Anthropogenic Influences on Dissolved Organic Matter in Three Coastal Bays, North China

Coastal bays bear anthropogenic influence strongly, and thus dissolved organic matter (DOM) in coastal bays, which is an important component of global carbon cycling, could be heavily affected by anthropogenic inputs. Utilizing absorbance, fluorescence spectroscopy, and stable carbon isotope (δ13C), this study analyzed the characteristics and distribution of DOM in three coastal bays (Jiaozhou Bay, Sishili Bay, and Taozi Bay), located in North China. The results showed that there was always a high concentration of DOM near the river mouth in all three bays and the DOM concentration decreased along the salinity gradient in Jiaozhou Bay, indicating the riverine inputs are the main factor that causes the variation of DOM in these coastal bays. The effects of inflowing rivers on DOM in coastal bays differed with their watershed characteristics (i.e., agricultural/urban). In addition, humic-like DOM components were found to be positively correlated with the apparent oxygen utilization, suggesting microbial activities could contribute to the DOM in this region. There was generally a higher averaged concentration of fluorescent DOM in surface water than that in bottom water in Jiaozhou Bay. In contrast, higher humic-like DOM was found in bottom water than that in surface water in Sishili Bay and Taozi Bay, which could be attributed to aquaculture activities and biological production. Moreover, photodegradation/photobleaching, dumping, and sewage discharge had their effects on DOM in coastal bays. This study demonstrates that DOM in coastal bays is regulated by multiple sources (rivers, aquaculture, dumping, and sewage) and processes (biological production and photodegradation), and anthropogenic activities have their influences on optical and isotopic characteristics of DOM in coastal bays.


INTRODUCTION
Estuarine and coastal areas link the land and the sea, thus playing important roles in carbon transport between these two largest carbon pools (Cai, 2011). Large amounts of terrestrial organic matter are transported to estuaries and coastal areas by rivers (Bianchi, 2011). Moreover, about 40% of the world's population lives within 100 km of the coast. Anthropogenic-derived nutrients and organic matter have been released to coastal areas, exerting a strong influence on aquatic ecosystems (eutrophication, algal bloom, hypoxia, etc.,) and carbon cycling . Despite the importance of coastal areas in global carbon cycling, coastal areas received less attention than estuarine regions (Li et al., 2015;He et al., 2020).
Dissolved organic matter (DOM) is the largest organic carbon pool in the ocean and has drawn increasing attention in the past decades (Hansell, 2015). Meanwhile, as the most reactive carbon pool, DOM could be significantly affected by anthropogenic activities and multiple biogeochemical processes (e.g., biodegradation and photodegradation) (Li et al., 2015;He et al., 2020;Lønborg et al., 2020). Previous studies have shown that in estuarine and coastal areas, riverine input would be one of the most important sources of DOM (Raymond and Bauer, 2001;Li et al., 2015;Guéguen et al., 2016). Besides, the input of riverine nutrients could stimulate primary production and thus influence the quantity and quality of DOM (Zhang and Yang, 2013;Asmala et al., 2018). Moreover, most marine aquaculture pools locate in coastal areas, releasing DOM into nearby waters (Zhao et al., 2001;Mahmood et al., 2017). Compared to large river estuaries, the water retention time in coastal areas, especially coastal bays, is relatively long due to hydrodynamic and geomorphological characteristics, allowing coastal areas to act as "reactors" as well as "pipes" of DOM (Casas-Ruiz et al., 2017;Liu et al., 2019).
Considering its complexity, multiple methods have been developed to study the quantity, quality, sources, and fates of DOM (Moran et al., 2016). Due to high sensitivity, fast response, easy use, and low cost, spectroscopic methods have been widely utilized to study DOM in the past several decades (Coble et al., 2014;Li and Hur, 2017). The ultraviolet-visible (UV-Vis) absorbance spectroscopy has been used to indicate the abundance and property of chromophoric DOM (CDOM) (Li and Hur, 2017). With high sensitivity, fluorescence spectroscopy has been extensively utilized to quantify and qualify DOM. Combining with parallel factor analysis (PARAFAC), fluorescence excitation-emission matrix (EEM) could be separated into independent fluorescent components. Besides, several indices were developed to indicate the source and property of fluorescent DOM (FDOM) (Zsolnay et al., 1999;McKnight et al., 2001;Cory and McKnight, 2005;Huguet et al., 2009). In addition to spectroscopic methods, stable carbon isotope analysis was also used to differentiate DOM from different sources (Cawley et al., 2012).
Due to rapid socio-economic development, most coastal areas in China are under increasing anthropogenic pressure. A significant fraction of treated domestic sewage has been discharged into coastal areas, increasing the concentration of organic matter and nutrients, and thus resulting in eutrophication and hypoxia (Strokal et al., 2014;Wang et al., 2018). Besides, heavy aquaculture activities occur in coastal areas, influencing the DOM in coastal bays , and also intensifying eutrophication and hypoxia Liang et al., 2018;Zhang et al., 2018). One recent study shows that the concentration of humic-like DOM could be elevated in the kelp mariculture zone in Ailian Bay, China . In addition, both eutrophication and hypoxia have a non-negligible effect on DOM dynamics in coastal areas (Zhang and Yang, 2013;Asmala et al., 2018;Zhang et al., 2018;Zhao et al., 2021). For example, Asmala et al. (2018) found that eutrophication could lead to the accumulation of recalcitrant DOM in coastal environments. Besides, hydrodynamics and geomorphology in coastal areas would affect the water exchange, thus exerting their effects on DOM dynamics (Jennerjahn and Ittekkot, 2002).
Coastal areas in Shandong Province are seriously affected by anthropogenic activities due to the fast development of the economy and aquaculture . Therefore, three coastal bays (Jiaozhou Bay, Sishili Bay, and Taozi Bay) in Shandong Province were selected to explore anthropogenic influences on DOM in coastal areas. Combining UV-Vis absorbance spectroscopy, fluorescence spectroscopy, and stable carbon isotope, this study aims to 1) investigate the spatial variations of DOM in coastal bays; 2) elucidate the sources of DOM in coastal bays; 3) evaluate the effects of anthropogenic activities on DOM in coastal bays.

Study Area
All three bays (Jiaozhou Bay, Sishili Bay, and Taozi Bay) in this study were situated in Shandong Province of northern China, and on the western coast of the Yellow Sea (Figure 1; Schlitzer, 2018). As a semi-closed bay, Jiaozhou Bay has a narrow mouth of 2.5 km and a shallow average depth of 7 m. Jiaozhou Bay is surrounded by Qingdao City and more than 10 rivers flow into the bay. On the western coast, Dagu River, Yang River, and Lianwan River flow into the bay. On the eastern coast, more urban rivers (Moshui River, Baisha River, Loushan River, Licun River, Haipo River, etc.,) flow into the bay, and these rivers served as sewage discharge channels for many years (Hu et al., 2020). Besides, many marine aquaculture pools (e.g., fish and shellfish) locate at the western area of Jiaozhou Bay.
Compared with Jiaozhou Bay, Sishili Bay, and Taozi Bay are relatively open bays near Yantai City, and are divided by Zhifu Island. Guangdang River, Xinan River, and Xiaoyuniao River flow into Sishili Bay. Jia River, Huangjin River, and Baiyin River flow into Taozi Bay. Meanwhile, both Sishili Bay and Taozi Bay are important marine aquaculture (e.g., scallop) areas in Shandong Province, China.
The study area experiences subtropical dry and wet seasons, and the wet season starts from June to September. In most years, there is the highest precipitation in July and August, therefore most of the rivers in our study have their highest discharge during this period.
frozen at -20°C before analysis. Temperature and salinity were obtained through a CTD (Conductivity, Temperature, Depth; SBE 25 plus) recorder (Seabird, United States). Dissolved oxygen (DO) concentration was obtained from CTD equipment on board. The DO sensor was calibrated against water sample analyses conducted by the Winkler titration method. Concentrations of Chl-a were measured following the method described by Arar and Collins (1997).

Dissolved Organic Carbon and Physicochemical Parameters
Dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) were measured on a total organic carbon (TOC) analyzer equipped with a Total Nitrogen Module (TOC-L, Shimadzu, Japan). The stable carbon isotopic composition (δ 13 C) was measured based on solid-phase extracted DOM (SPE-DOM) on an isotope ratio mass spectrometer (Finnigan MAT 253, Thermo Scientific, United States). The solid-phase extraction of DOM was conducted using PPL cartridges (Agilent Bond Elut PPL, 200 mg, 3 ml) according to the method developed by Dittmar et al. (2008). The stable carbon isotopic values are reported in δ 13 C notation (‰) relative to the Vienna Pee Dee Belemnite with precision and accuracy within 0.2‰.

Ultraviolet-visible Absorbance Spectra
The UV-Vis absorbance and fluorescence spectra were determined on an Aqualog (Horiba) fluorescence spectrometer equipped with a 1 × 1 cm quartz cuvette. The absorbance spectra were measured from 240 to 600 nm with an interval of 3 nm. Absorbance spectra of milli-Q water were used as a reference and average absorbance between 580 to 600 nm was used for baseline correction. The absorption coefficient (a 351 ) was calculated to quantify the CDOM abundance. The spectral slope (S 275-295 ) was calculated through the linear fitting of log-transformed absorption coefficient according to Helms et al. (2008) to indicate the source and quality of CDOM. The low S  suggests that the DOM contains a high proportion of terrestrial-derived organic matter and S 275-295 would increase due to photobleaching (Fichot and Benner, 2012).

Fluorescence Spectra
The three-dimensional fluorescence spectra (excitation-emission matrices, EEMs) were scanned across 240-600 nm at 3 nm intervals for excitation wavelengths (Ex) and 248.1-825.1 nm at 3.7 nm intervals for emission wavelengths (Em). The integration time was 0.5 s. Milli-Q water was scanned every day for blank subtraction and Raman normalization. All spectra were blank subtracted and corrected for inner-filter effect with the Horiba Aqualog software. Also, fluorescence index (FI, the ratio of the fluorescence intensity at 470 and 520 nm when Ex 370 nm), biological index (BIX, the ratio of the fluorescence intensity at 380 and 430 nm when Ex 310 nm), and humification index (HIX, the ratio of the sum fluorescence intensity in the 300-345 nm and 435-480 nm when Ex 254 nm, and Ex 255 nm is used herein) were calculated based on fluorescence spectra to indicate the source and quality of DOM (Zsolnay et al., 1999;McKnight et al., 2001;Cory and McKnight, 2005;Huguet et al., 2009). The FI has a value of ∼1.55 for microbially derived fulvic acid and ∼1.21 for terrestrially derived fulvic acids (Cory et al., 2010). The BIX could indicate the contribution of autochthonous organic matter and the high BIX value suggests a high contribution of autochthonous organic matter Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 697758 (Huguet et al., 2009). The HIX increases with increasing humification degree, and could indicate the contribution of terrestrial organic matter (Zsolnay et al., 1999;Zhang et al., 2010).

Parallel Factor Analysis
A fluorescence dataset of 173 EEMs was analyzed by PARAFAC. The PARAFAC modeling was conducted using the DOMFluor toolbox on MATLAB R2020b (Mathworks, United States) according to the protocol described by Stedmon and Bro (2008). Before PARAFAC modeling, Raman normalization was conducted to enable comparison among different studies (Murphy et al., 2010), therefore the fluorescence intensity in this study was reported as Raman Unit (RU). In addition, to facilitate the modeling speed and to avoid the influence of scattering, the EEMs for PARAFAC modeling only covered the excitation range of 252-399 nm and the emission range of

Physicochemical Parameters
The physicochemical parameters in the three bays are summarized in Table 1. The salinity ranged from 26.01 to 31.59 in Jiaozhou bay, while the salinity range in Sishili Bay and Taozi Bay was much narrower ( Table 1). In Jiaozhou Bay, the lowest salinity was found in the northeastern area where the bay received discharges from urban rivers ( Figure 2). In Sishili Bay (31.24-32.01) and Taozi Bay (31.49-31.95), the lowest salinity was also found near the river mouth along the coast (Figure 3). In contrast, the temperature in Jiaozhou Bay was relatively stable from surface to bottom, while the temperature showed a depth gradient in Sishili Bay and Taozi

Dissolved Organic Matter and Chromophoric Dissolved Organic Matter
The DOC concentration had a range of 147.7-267.5 μmol L −1 in our study (Table 2), similar to the results reported in the same study area (Zhao et al., 2001;Zhang et al., 2009;Zhang et al., 2018). The SPE-DOM δ 13 C had a narrow range from −24.8‰ to −22.9‰. There was higher DOC concentration and lower SPE-DOM δ 13 C in bottom water than in surface water in Taozi Bay. There was higher CDOM abundance (a 351 ) in surface water than bottom water in Jiaozhou Bay, while there was no significant difference (p > 0.05) in CDOM abundance among different layers in Sishili Bay and Taozi Bay. There was lower S 275-295 in surface water than bottom water in Jiaozhou Bay and Taozi Bay, while there was higher S 275-295 in surface water than bottom water in Sishili Bay.

Fluorescent Dissolved Organic Matter
Four fluorescent components were identified by PARAFAC, including one protein-like (C1) and three humic-like (C2, C3, and C4) components ( Figure 4). All four components were matched with the components identified in previous studies based on the OpenFluor database. C1, resembling traditionally defined peak B, has been reported as one tyrosine-like fluorophore (Li et al., 2015). C2 has two excitation maxima at <252 and 309 nm with its emission maxima at 393 nm. C2 Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 697758 resembles traditionally defined peak M, and has been regarded as microbially derived humic-like substances (Li et al., 2016). C3 resembles a combination of traditionally defined peaks A and C, and has been reported as terrestrially derived humic-like substances (Li et al., 2015). C4 has its excitation maxima <252 nm with its emission maxima at 453 nm, resembling traditionally defined humic-like peak A (Coble, 2007). C4 has also been regarded as terrestrially derived humic-like substances.
The statistics of four identified PARAFAC components are shown in Table 2. There was a higher FDOM concentration (C1-C4) in surface water than bottom water in Jiaozhou Bay. In contrast, there were higher concentrations of humic-like FDOM components in bottom water than surface water in Sishili Bay and Taozi Bay.
FI had similar ranges in three bays, with a range of 1.46-1.69 in Jiaozhou Bay, 1.34-1.61 in Sishili Bay and 1.44-1.57 in Taozi Bay. BIX had a range of 0.89-1.48 in all samples from Jiaozhou Bay, Sishili Bay and Taozi Bay, suggesting a dominant autochthonous origin in the study area, which is supported by the low HIX (0.89-2.24) in these three bays. There was higher HIX and lower BIX in bottom water than in surface water in Sishili Bay.

Riverine Inputs of Dissolved Organic Matter
In Jiaozhou Bay, high DOM abundance (DOC, a 351 , and four fluorescent components) was found in the northeastern and western area where rivers (Dagu River, Yang River, Moshui River, Loushan River, Licun River, and Haipo River) flow into the bay, suggesting riverine inputs of DOM (Zhang et al., 2009;Zhang and Yang, 2013;Hu et al., 2020) (Figure 5). Moreover, negative correlations between DOM abundance and salinity were identified ( Table 3), suggesting that conservative mixing between freshwater and seawater controlled the distribution of DOM in Jiaozhou Bay (Zhang and Yang, 2013). In addition, the correlation coefficients varied among different fluorescent components (C1 to C4; Table 3), suggesting that the degree of the influence of mixing on DOM varied among different components.
In Sishili Bay and Taozi Bay, high DOM abundance was found in the nearshore area, especially in the area near the mouth of Jia River and Xin'an River, indicating a contribution of riverine input ( Figure 6). However, the riverine input had a limited influence on DOM abundance in the whole bay area, which could be corroborated by the narrow salinity ranges observed in these two bays (Table 1), consistent with a previous study . The limited influence of riverine discharge on DOM abundance could be possibly explained by the characteristics of the inflowing rivers and the bays. First, the inflowing rivers are relatively small and have low flow rates. Moreover, these two bays are relatively open, allowing an exchange of water between the inner and outer bay (Dong et al., 2019).
Riverine inputs had their signatures on the optical properties and stable isotopic properties of DOM in Jiaozhou Bay. Lower S 275-295 values were found in the northeastern and northwestern part of Jiaozhou Bay (Figure 7), which could be attributed to the input of terrestrial-derived DOM from inflowing rivers (Helms et al., 2008). Higher FI and BIX values were found in the eastern part of Jiaozhou Bay (Figure 7), indicating a stronger input of microbial-derived source or an autochthonous source (Huguet et al., 2009;Cory et al., 2010). Lower HIX was found in the eastern part of Jiaozhou Bay (Figure 7), suggesting a lower humification degree of DOM, corroborating with higher FI and BIX. In contrast, relatively higher HIX was found in the western part of Jiaozhou bay, where Dagu River and Yang River flow into the bay (Figure 7). On the western coast of Jiaozhou Bay, the land cover is dominated by agricultural land, while on the eastern coast, the land cover is dominated by urban areas (Ke et al., 2020). Previous studies have shown that rivers draining agricultural land had a higher concentration of humic-like components and higher HIX, while rivers flowing through the urban area had higher concentrations of protein-like and microbial-derived humic-like components, and higher FI and BIX values due to domestic wastewater input (Shang et al., 2018;Chen et al., 2019). In addition to optical properties, the values of SPE-DOM δ 13 C indicate a mixture of terrestrial-and marine phytoplanktonderived DOM (Cawley et al., 2012), and they differed in these two areas of Jiaozhou Bay. Lower SPE-DOM δ 13 C values were found in the northeastern part of the bay, which could be attributed to the input of 13 C-depleted terrestrial organic matter (e.g., C3 plants or freshwater phytoplankton) through riverine discharge. This pattern is similar to the distribution of particulate organic matter (POM) δ 13 C observed during summer in the same area (Ke et al., 2020).
In addition to land cover, sewage discharge, aquaculture and dumping also had their effect on DOM. DOC and three humic components (C2, C3, and C4) exerted high concentrations in the southeast side of Zhifu Island (Figure 6), which could be attributed to the discharge of domestic and industrial sewage (Wang et al., 2012). Meanwhile, in southeast stations of Sishili Bay, high concentrations of DOC, C1 and C2 could be possibly attributed to aquaculture activities, consistent with previous studies (Zhao et al., 2001;Li et al., 2018). Moreover, high CDOM (a 351 ) abundance was found in the northern part of Sishili Bay, which could be attributed to marine dumping activities (Wang et al., 2012).

Multiple Mechanisms Controlling the Vertical Distribution of Dissolved Organic Matter
In estuarine and coastal regions, DOM concentration is generally dominated by the mixing between freshwater and seawater. Therefore, there is often higher DOM in low salinity surface water than that in bottom water due to freshwater input (Li et al., 2015). As expected, high DOM abundance (DOC, C1-C4) was observed in surface water in most stations of Jiaozhou Bay (Table 2; Figure 8). Considering similar salinity observed in surface and bottom water in Taozi Bay and Sishili Bay, DOM concentration in surface and that in bottom water should also be similar. However, there was a higher average DOM abundance in bottom water than in surface water in Sishili Bay and Taozi Bay (Table 2; Figure 8). Higher DOC and FDOM concentrations were found in bottom water than that in surface water in Taozi Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 697758 Bay. In Sishili Bay, three humic-like components were also found to be higher in bottom water than that in surface water. Theoretically, two possibilities could cause higher DOM abundance in bottom water than in surface water. First, DOM in surface water decreased due to a series of removal processes (e.g., photodegradation and microbial degradation). Second, DOM in bottom water increased due to additional inputs (e.g., sediment resuspension and subsequent release) (Li et al., 2015).
To check which is the main reason behind the vertical  Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 697758 8 distribution of DOM in these two bays, a cross-bay comparison has been made for all three bays of this study. A similar concentration of humic-like components (C2-C4) was found in surface water samples with similar salinity in all three bays, while higher humiclike components were found in bottom water samples of Sishili Bay and Taozi Bay than that in bottom water samples of Jiaozhou Bay (Figure 9). These results suggest that the second one (additional inputs at bottom layer) should be the main underlying reason.
Previous studies suggested that sediment resuspension could be an important source for DOM in estuarine and coastal regions (Li et al., 2015). Due to shallow depth, sediment resuspension triggered by wind or tidal forces frequently occurs in estuarine and coastal regions, driving the release of DOM from sediment pore water into the overlying water, and thus resulting in high DOM concentration in bottom water. The DOM abundance, especially protein-like components, was found to be elevated in the Yangtze estuary due to sediment resuspension (Li et al., 2015). However, a previous study showed that the tidal current inside Sishili Bay is weak (Li et al., 2013). Besides, the wind is not strong in summer, which could be supported by the depth gradient of temperature in Sishili Bay and Taozi Bay. Moreover, scallops are cultivated in suspended lantern nets, not disturbing sediments. Therefore, wind or tidal forces or aquaculture could not be able to drive sediment resuspension, and thus sediment resuspension and subsequent release of DOM should not be the main mechanism for high DOM concentration in bottom water in this study.
Aquaculture activities could be a possible reason behind high DOM abundance in bottom water in Sishili Bay and Taozi Bay. High DOC concentration has been reported in the aquaculture area of Sishili Bay (Zhao et al., 2001). Previous studies have found high DOC concentration at bottom water of aquaculture area and attributed it to accumulation due to long-term aquaculture activities (Mahmood et al., 2017;Li et al., 2018). Besides, elevated humic-like components have also been found in bottom water in areas of kelp culture , similar to the distribution observed in Sishili Bay and Taozi Bay. However, it should be noted that the effects of aquaculture on DOM depend on the species of culture and the growth phase (Mahmood et al., 2017). In this study, there is no significant difference between surface DOC and bottom DOC in Sishili Bay. In Sishili Bay and Taozi Bay, scallops are the main cultural species. A previous study has shown that scallops could either increase or decrease DOC concentration (Mahmood et al., 2017). The increase of DOC was mostly attributed to DOC excretion, while the decrease was mostly caused by DOC uptake through filter-feeding (Mahmood et al., 2017). Besides sediment resuspension and aquaculture activities, biological production could also be another important source of DOM at the bottom layer. Humic-like components could be produced in deep water through biological production based on a significant correlation (p < 0.001) between humic-like components and AOU (Yamashita and Tanoue, 2008). In this Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 697758 10 study, significant correlations (p < 0.05) between AOU and humic-like components (C3 and C4 in Jiaozhou Bay, C2-C4 in Sishili Bay) are also found in the bottom water of Jiaozhou Bay and Sishili Bay ( Table 4), indicating that biological production is an important source of DOM and could be the main reason behind the vertical distribution of DOM in our study area. It should be noted that there are only three bottom samples available in Taozi Bay, which could be the reason that there is no correlation found between humic-like components and AOU in Taozi Bay.
Although aquaculture and biological production could be the main reasons behind high DOM abundance at bottom water in Taozi Bay and Sishili Bay, photodegradation/ photobleaching also has its effect on the distribution of DOM. There are some stations with higher humic-like components (C2-C4) in bottom water than in surface water in Jiaozhou Bay, although there is higher DOM abundance in surface water than in bottom water in most stations of Jiaozhou Bay. These results are consistent with that humic-like components are more susceptible to photodegradation (Lønborg et al., 2015;Yang et al., 2020). Moreover, there are higher S 275-295 in surface water than that in bottom water in most stations of Sishili Bay, indicating photobleaching is an important factor influencing surface DOM.

Primary and Microbial Production
In Sishili Bay, positive correlations have been identified between Chl-a and three humic-like components in surface water ( Table 4), suggesting that primary production could be an important source of DOM. Correlation between Chl-a and DOM abundance (DOC) has also been reported in the surface layer of Jiaozhou Bay (Zhang and Yang, 2013), partly corroborating our results. Also, Chl-a was found to be positively related to TDN, indicating TDN was a major controlling factor of primary production in this region (Wetz et al., 2017). Moreover, correlations have been found between AOU and three humic-like components in surface water. These correlations have been found in the deep ocean, and have been attributed to biological oxidation of organic matter and subsequent production of humic-like organic matter (Yamashita and Tanoue, 2008). In most cases, there was no correlation between AOU and humiclike components observed in the surface water due to the air-sea gas exchange. The correlation between AOU and humic-like components observed in the surface water of Sishili Bay corroborated that the tidal current within the bay is weak so that the oxygen consumed by microbial activities cannot be immediately replenished through the air-sea exchange. Therefore, we could conclude that the TDN Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 697758 concentration could affect, to some extent, the concentration of humic-like components in surface water in Sishili Bay.

CONCLUSION
Multiple sources and processes regulate DOM in coastal bays in this study. Riverine inputs are one of the most important DOM sources in coastal bays, and the effects vary with the discharges and watershed characteristics (agricultural and urban activities) of inflowing rivers, and also the geomorphology of the bays. There is a higher DOM abundance, especially humic-like FDOM components, in the bottom water of Sishili Bay and Taozi Bay, which could be attributed to aquaculture activities and biological production at the bottom layer. Besides, significant correlations have been found between physicochemical parameters (TDN, Chl-a, and AOU) and DOM abundance (humic-like components), indicating eutrophication and subsequently enhanced primary and microbial production are important factors controlling DOM abundance and distribution in coastal bays. Moreover, photodegradation/photobleaching, dumping and sewage discharge had their effects on DOM in this study. This study demonstrates that DOM in coastal bays is strongly influenced by anthropogenic activities and that spectroscopic and isotopic methods are useful for elucidating the source