Characterization of Marine-Surface-Dissolved Organic Matter via Amino Acid Enantiomers and Its Implications Based on Diel and Seasonal Observations

: Due to the essential roles of dissolved organic matter (DOM) in both microbiol food loop and marine carbon cycling, changes in marine DOM composition have an important impact on the marine ecosystem and carbon cycling. In October 2014 and June 2015, two ﬁeld investigations for the DOM in the upper 200 m were conducted in the slope region of the northern South China Sea to characterize the DOM composition via amino acid enantiomers. In June, our sampling locations were under upwelling impact induced by an eddy-pair event, whereas in October there were no eddies. High-frequency sampling (a few hours interval) over 24 h reveals that the variability of the amino acid carbon yield (min. 0.2%) and the D/L alanine ratio (min. 0.02) is larger than its corresponding analytical and propagated errors, suggesting solid short-term changes for these two molecular-based indicators. Section samples from June showed a lower D/L alanine ratio (0.43 vs. 0.53) and a GABA mol% (1.0% vs. 1.6%) relative to the section samples from October, suggesting that DOM in June is more fresh (less degraded) compared to that in October. A higher serine mol% (19.5% vs. 13.2%) and lower D/L serine ratio (0.06 vs. 0.24) from the diel observation in June relative to October further indicates that phytoplankton, rather than bacteria, plays an more important role in DOM composition alternation. This is consistent with the higher phytoplankton biomass found in June, promoted by the eddy-pair. , samplings


Introduction
Differences in the composition of dissolved organic matter (DOM) among samples provide an indication of sample origins and reflect the refractory or labile nature of DOM [1,2]. Therefore, understanding the composition of DOM is an essential requirement in studies of carbon and nitrogen cycling, e.g., [3,4]. Amino acids (AAs) account for a small fraction of the bulk DOM, but they offer a powerful tool for determining the composition of organic matter and are widely used in contemporary organic matter studies, e.g., [1,[5][6][7][8]. AAs are derived mainly from living organisms, and variations in both their total amounts and relative abundances are tightly coupled to the metabolism of organisms. Therefore, both the concentrations and relative abundances of AAs are dependent on the status or progress

Field Sampling
Two cruises were undertaken in R/V Nan Feng on the northern slope region of the SCS during 10-30 October 2014 and 10-30 June 2015. The temperature (T), salinity (S), fluorescence, and photosynthetically active radiation (PAR; in June 2015 only) were measured with a CTD (Seabird, USA). Discrete water samples were collected in depth profiles (top 200 m) along a section that covers the SCS northern slope area ( Figure 1). In addition, time-series stations (TS2014 in 2014 and TS2015 in 2015) were also observed. For the timeseries stations, sampling began at 06:00 and lasted for 1 day at TS2014. At TS2015, sampling began at 12:35 and also lasted for 1 day. In both cruises, the time-series observations were carried out 1 week later when compared to the corresponding section observations. Up to five profile samplings were conducted during the time-series stations in each of the two seasons. The sampling interval was approximately 6 h.
Immediately after collection, the seawater samples for the analysis of AAs and DOC were filtered with a clean nylon membrane (pore size 0.45 μm) and a syringe, and frozen (−20 °C) before analysis in the laboratory. Water samples for the analysis of dissolved nutrients were first filtered through acid-cleaned acetate cellulose filters (pore size 0.45 μM), treated with HgCl2, and stored at 4 °C in the dark before analysis. Nutrient samples were not collected during the last two samplings at TS2014.

Background
The benefits of multidisciplinary investigations and the physical, chemical, and biological background of this work are basically revealed [25]. In June 2015, an eddy-pair was found during our diel observations, whereas in October 2014, no mesoscale process was detected ( Figure 1) [26]. The observation sites were impacted by weak upwelling (Figure 1) [26]. Consequently, the surface phytoplankton biomass (0.12-0.48 µg L −1 ) was higher in June 2015 than in October 2014, when the biomass was as low as 0.02 µg L −1 [27]. In June 2015, the phytoplankton community was dominated by Chaetoceros and Thalassionema, both of which are diatoms [27]. Although their abundance remained unclear in previous work, when the bacterial diversity was studied in October 2014 the community was dominated by heterotrophs, and the abundance of autotrophic bacteria was rather low [28].

Field Sampling
Two cruises were undertaken in R/V Nan Feng on the northern slope region of the SCS during 10-30 October 2014 and 10-30 June 2015. The temperature (T), salinity (S), fluorescence, and photosynthetically active radiation (PAR; in June 2015 only) were measured with a CTD (Seabird, USA). Discrete water samples were collected in depth profiles (top 200 m) along a section that covers the SCS northern slope area ( Figure 1). In addition, timeseries stations (TS2014 in 2014 and TS2015 in 2015) were also observed. For the time-series stations, sampling began at 06:00 and lasted for 1 day at TS2014. At TS2015, sampling began at 12:35 and also lasted for 1 day. In both cruises, the time-series observations were carried out 1 week later when compared to the corresponding section observations. Up to five profile samplings were conducted during the time-series stations in each of the two seasons. The sampling interval was approximately 6 h.
Immediately after collection, the seawater samples for the analysis of AAs and DOC were filtered with a clean nylon membrane (pore size 0.45 µm) and a syringe, and frozen (−20 • C) before analysis in the laboratory. Water samples for the analysis of dissolved nutrients were first filtered through acid-cleaned acetate cellulose filters (pore size 0.45 µM), treated with HgCl 2 , and stored at 4 • C in the dark before analysis. Nutrient samples were not collected during the last two samplings at TS2014. At TS2014, water samples were further collected for the detection of bacterial abundance at 15:00 and 21:00 on 26 October (i.e., 3 h after the molecular indicators were sampled), whereas in late spring (June) at TS2015, the water samples for bacterial abundance were collected at the same time as the samples used to analyze the molecular indicators (one sampling point was missed: the 2nd sampling point). A volume of 47.5 mL of seawater was collected and fixed with 2.5 mL of formaldehyde for 15 min. The samples were stored at 2-8 • C and analyzed in the laboratory.

Laboratory Measurements
All the measured AAs (including both L-and D-enantiomers) and their abbreviations are shown in Table 1. The THDAAs were analyzed with a previously described method [29]. Briefly, the water samples were first hydrolyzed in sealed ampoules with 6 M HCl at 110 • C under N 2 . For samples with elevated nitrate levels, ascorbic acid was also added [30]. After hydrolysis, the samples were neutralized and the pH adjusted to 8.5. After precolumn derivatization with o-phthaldiadehyde (OPA) and N-isobutyryl-L/D cysteine, the AA enantiomers were separated and measured in the hydrolysates with high-performance liquid chromatography (HPLC; 1200 series Agilent, USA), with a fluorescence detector. During hydrolysis, asparagine (Asn) and glutamine (Gln) were deaminated to aspartic acid and glutamic acid, respectively. Therefore, "Asx" represents Asn plus aspartic acid (Asp), and "Glx" represents glutamine plus glutamic acid (Glu) ( Table 1). The DOC concentrations were measured with a total organic carbon (TOC) analyzer (TOC-L CPH , Shimazu, Japan). Ammonium was measured manually by the sodium hypobromite oxidation method, with an analytical precision of 0.04 µM. Concentrations of the other four nutrients were determined with an auto-analyzer (AA3: SEAL Analytical, Mequon, WI, USA), with precisions for nitrate, nitrite, dissolved inorganic phosphorus (DIP, PO 4 3− ), and silicate (SiO 3 2− ) of 0.01, 0.003, 0.005, and 0.02 µM, respectively. Milli-Q water was used in all laboratory processing, and the chemicals were purchased from Sigma-Aldrich (Missouri, USA) or Merck (USA) and were of HPLC grade or above. Table 1. Measured amino acids and abbreviations used in this study (L-and D-enantiomers are not listed). Total hydrolysable dissolved amino acids  THDAAs  Alanine  Ala  Arginine  Arg  Asparagine  Asx  Aspartic acid  Glutamine  Glx  Glutamic acid  Glycine  Gly  Isoleucine  Ile  Leucine  Leu  Lysine  Lys  Methionine  Met  Phenylalanine  Phe  Serine  Ser  Threonine  Thr  Tyrosine  Tyr  Valine Val γ-aminobutyric acid GABA 2-aminobutyric acid Aba

Name Abbreviation
To determine bacterial abundance, 1 mL 4 ,6-diamidino-2-phenylindole solution (DAPI; 10 µg mL −1 ) was added to 30 mL of each sample, and the mixtures were incubated in the dark for 5-10 min. Each sample was then filtered through 25 mm, and 0.2 µm black Whatman ® Nuclepore™ Polycarbonate Filters were mounted on pre-wetted Water 2021, 13, 685 5 of 16 0.8 Nuclepore Membra-fil ® cellulosic filters [31]. The filters were then set onto glass slides with a drop of Leica type A immersion oil. The bacterial abundances were determined with an epifluorescence microscope (Leica, DM5000B).

The AA-Based Molecular Indicators and Propagated Error Quantification
Among the AA-based molecular indicators, AA carbon yield is the AA carbon (in nM) divided by bulk DOC (in nM). AA% is the mol percentage of a given AA in THDAA. The degradation index (DI) is a measure that tells the overall degradation status of organic matter. The DI, and its calculation equation, was first proposed by Dauwe and Middelburg (1998) [32] for use in particulate and sediment studies. Later, it was applied to DOM with a revised constant derived from updated DOM samples [21]. In this work we calculated DI following the method of Dauwe and Middelburg (1998) [32], using the revised coefficients derived from DOM sample pools [21,33].
To check the errors in the measurements of the AA enantiomers, a water sample (10L) obtained at a depth of 5 m near TS2014 was sealed and stored for 3 years at room temperature in the dark to ensure it reached a stable and homogeneous state. Three subsamples were then taken from the 10 L sample. The concentrations of the THDAAs in the subsamples were determined, using the sampling and laboratory procedures described above. The standard deviation of the three subsamples was regarded as the error of the primary measurement of DCAAs.
To quantify the errors in measuring DOC, we used the deep-sea water standard (CRMbatch 9) from the Consensus Reference Materials project led by Prof. Dennis Hansell, at the University of Miami (FL, USA). Up to eight samples were analyzed. Similarly, the standard deviation of the eight DOC values was regarded as the error of the primary measurement of DOC.
When deriving the AA-based molecular indicators (D/L ratio, Gly%, percentage γaminobutyric acid (GABA%), AA carbon yield, and DI) from the primary measurements (i.e., individual AA concentrations and/or DOC concentration), the errors within the primary AA and/or DOC measurements were propagated to these molecular indicators. To obtain the equations for the propagation of uncertainties (as standard deviations), a statistical inference was conducted for each molecular indicator. The equations used to calculate the propagated errors are provided in the Supplementary Materials. In October 2014, a fluorescence peak was identified at 66 m depth at noon, decreasing both upwards and downwards. In June 2015, the fluorescence peak occurred at a depth of 40-54 m, and the depth at which the PAR value was 1% of that at the surface (i.e., the depth of the euphotic zone) was 56 m at that time.
The nutrient levels differed between June and October. In the surface waters, dissolved inorganic nitrogen (DIN) was slightly higher in June than in October (0.3 µM and 0.1 µM, respectively; Table 2). The DIP concentration was higher in October 2014 (mean 2.1 µM at the surface) than in June 2015, when it ranged from 0.03 µM to below the detection limit. A similar seasonality was observed at 200 m depth, with higher DIN and lower DIP in June but lower DIN and higher DIP in October (Table 2). Silicate was also clearly higher in June (2.0 µM) than in October (0.03 µM) in the surface waters and at a depth of 200 m (20 µM and 0.8 µM, respectively ( Table 2). The DIN/DIP ratio in the surface waters was very low (0.05) in October and higher (7.6) in June. At 200 m depth, the DIN/DIP ratio was 0.8 in  (Table 2). Over the entire upper 200 m water column (i.e., extending between 0 and 200 m), the DIN/DIP ratio was 0.48 and 16.66 over our diel observation period in October and June, respectively. In the surface waters, the bacterial abundance in October was 10 times that in June, with average values of 2.3 × 10 8 cell L −1 and 2.3 × 10 7 cell L −1 , respectively ( Table 3). The bacterial abundance decreased vertically in October; at 200 m depth, it was only half or one quarter that in the surface waters. In June, when the bacterial abundance in the surface waters remained low, it was similar from the surface to 200 m depth (or even increased).  [27]; the key algal groups in both months were diatoms.

Temporal Variability of DOM over 24 h
The levels of DOC in surface waters were similar in October and June (mean: 83 µM), but their variability over 24 h differed. In October, DOC ranged between 80 µM and 86 µM, whereas in June, it ranged between 75 µM and 87 µM ( Table 4). The ranges (6 µM in October and 12 µM in June) and the diel variability (standard deviations: 3 µM for October and 5 µM for June; Table 4) indicate that DOC varied more weakly in October than in June, and both values are clearly larger than our measurement precision (1 µM; Table S2). THDAAs in surface waters averaged 366 ± 37 nM in October, much lower than in June (584 ± 57 nM; Table 4). Similarly, the AA carbon yield was also lower in October (1.6% ± 0.2%) than in June (2.3% ± 0.3%). The diel variability in surface THDAAs and the AA carbon yield were smaller in October (37 nM and 0.2%, respectively) than in June (57 nM and 0.3%, respectively), as those of DOC (Table 4).
The diel THDAA pattern seems to have been decoupled from the variation in DOC. In October, THDAAs in the surface waters remained almost unchanged from 06:00 to midnight (24:00), when the concentration was slightly increased (Figure 2b). In June, the THDAA pattern did not mirror the corresponding DOC pattern either, but dropped steadily from 12:00 till the next dawn (06:00) and then increased again to 12:00 on the following day ( Figure 2b). The AA carbon yield pattern was very similar to that of THDAAs (Figure 2c).
Unlike the decoupling of THDAAs and DOC, the variations in the molecular indicators seemed to reflect some variations in DOC in October. At 18:00, when DOC in the surface water was highest, Gly%, D/L Ala, and DI in the surface waters were also at their maxima over the observation period (Figure 2d-f). Moreover, at night (24:00) and dawn (06:00), when DOC was low, Gly%, D/L Ala, and DI were also relatively low (Figure 2d-f). However, in June, there was no comparable coupling between DOC and the molecular indicators; instead, the indicators showed the reverse pattern relative to DOC (e.g., Gly% was elevated at 18:00 when DOC was low, and remained depleted at 06:00 when DOC was high; Figure 2d).
surface water was highest, Gly%, D/L Ala, and DI in the surface waters were also at their maxima over the observation period (Figure 2d-f). Moreover, at night (24:00) and dawn (06:00), when DOC was low, Gly%, D/L Ala, and DI were also relatively low (Figure 2df). However, in June, there was no comparable coupling between DOC and the molecular indicators; instead, the indicators showed the reverse pattern relative to DOC (e.g., Gly% was elevated at 18:00 when DOC was low, and remained depleted at 06:00 when DOC was high; Figure 2d).

Section Observations
In surface waters, vertical mixing was deeper in October than in June, as indicated by the density distributions [quantified as σ in kg/m 3

Section Observations
In surface waters, vertical mixing was deeper in October than in June, as indicated by the density distributions [quantified as σ in kg/m 3 , = 1000 × (seawater density-1.  (Figure 3). The DCM belt roughly followed the waters with density 23-23.5 kg/m 3 in both months (Figure 3), occurring at a depth between 50 to 100 m. The fluorescence signal was more intensive in June relative to that in October. in June 2015 (Figure 3). A clear deep chlorophyll maximum (DCM) belt can be found along the section in both October and June (Figure 3). The DCM belt roughly followed the waters with density 23−23.5 kg/m 3 in both months (Figure 3), occurring at a depth between 50 to 100 m. The fluorescence signal was more intensive in June relative to that in October.
(a) (b) Overall, THDAA was higher in June relative to October. In October, THDAA ranged from 181 nM to 317 nM, with an average value of 252 nM, whereas in June the mean value was 452 nM (ranging from 257 to 596 nM) ( Table 5). In June, the elevated THDAA values (e.g., >452 nM) basically appeared in the upper 100 m waters. The D/L Ala ratio in June was lower (0.43) when compared to that in October (0.53) ( Table 5). In October, elevated D/L Ala ratios roughly overlapped with the DCM zone, but in June such overlap pattern disappeared. Though the discrete samples cannot ideally cover the DCM belt, as revealed by the CTD probe, the D/L Ala ratio was low (e.g., <0.4) on the DCM belt zone in June (Figures 3 and 4). Gly% and GABA% in June were also lower than in October ( Table 5). As shown in Table 5, Gly% (GABA%) was 6% (0.6%) lower in June relative to October, judged by the mean values. Spatially, both Gly% and GABA% indicate a good DOM composition vertical mixing in June, whereas in October the DOM composition showed a horizontally laminated distribution feature ( Figure 4). As for DI values, a higher DI was observed in June (mean: 0.49) relative to October (mean: 0.20) (Table 5). Overall, THDAA was higher in June relative to October. In October, THDAA ranged from 181 nM to 317 nM, with an average value of 252 nM, whereas in June the mean value was 452 nM (ranging from 257 to 596 nM) ( Table 5). In June, the elevated THDAA values (e.g., >452 nM) basically appeared in the upper 100 m waters. The D/L Ala ratio in June was lower (0.43) when compared to that in October (0.53) ( Table 5). In October, elevated D/L Ala ratios roughly overlapped with the DCM zone, but in June such overlap pattern disappeared. Though the discrete samples cannot ideally cover the DCM belt, as revealed by the CTD probe, the D/L Ala ratio was low (e.g., <0.4) on the DCM belt zone in June (Figures 3 and 4). Gly% and GABA% in June were also lower than in October ( Table 5). As shown in Table 5, Gly% (GABA%) was 6% (0.6%) lower in June relative to October, judged by the mean values. Spatially, both Gly% and GABA% indicate a good DOM composition vertical mixing in June, whereas in October the DOM composition showed a horizontally laminated distribution feature ( Figure 4). As for DI values, a higher DI was observed in June (mean: 0.49) relative to October (mean: 0.20) (Table 5).

Solid Short-Term Changes of AA-Based Molecular Indicators
Overall, our results show that the propagated error played a minor role in most AAbased molecular indicators (Table S3). Though the propagated error for DI can be very high, this is not strange, as the DI value itself can be very close to zero [20,33,34]. Instead, we suggest a careful consideration of GABA% (Table S3). The absolute value of the propagated error for GABA% was of the order of 0.3-0.4%, which was similar to, or even larger than, the GABA% value itself, which was 0.3-2.0% in this study. Therefore, in studies using GABA% as a molecular indicator of surface DOM, this error should be considered. In deep waters, however, the impact of the propagated error on GABA% should be minor because GABA% increases with depth and can be >15% at water depths of >800 m [20].
The precision of the DOC and THDAA measurements was 1 μM and 20 nM, respectively, far lower than the corresponding short-term changes over the 24 h observation pe-

Solid Short-Term Changes of AA-Based Molecular Indicators
Overall, our results show that the propagated error played a minor role in most AAbased molecular indicators (Table S3). Though the propagated error for DI can be very high, this is not strange, as the DI value itself can be very close to zero [20,33,34]. Instead, we suggest a careful consideration of GABA% (Table S3). The absolute value of the propagated error for GABA% was of the order of 0.3-0.4%, which was similar to, or even larger than, the GABA% value itself, which was 0.3-2.0% in this study. Therefore, in studies using GABA% as a molecular indicator of surface DOM, this error should be considered. In deep waters, however, the impact of the propagated error on GABA% should be minor because GABA% increases with depth and can be >15% at water depths of >800 m [20].
The precision of the DOC and THDAA measurements was 1 µM and 20 nM, respectively, far lower than the corresponding short-term changes over the 24 h observation period (Figure 2a,b). For molecular indicators, a comparison of the short-term changes (expressed as standard deviations) and the corresponding seasonal propagated errors indicated that short-term changes were significant for AA carbon yield, Gly%, and D/L Ala ratio, whereas the other two molecular indicators were not (e.g., GABA%; Figure 5). The variability (as a standard deviation) in GABA% over 24 h was similar to, or even smaller than, the propagated error ( Figure 5).
The bulk DOC and THDAA changes over 24 h were clearly greater in June 2015 than in October 2014 (Figure 2a,b): for surface waters in June 2015, the standard deviations were 5 µM (DOC) and 57 nM (THDAAs), whereas in October 2014, they dropped to 3 µM (DOC) and 37 nM (THDAAs) ( Table 4). The short-term changes in THDAAs (maximum-minimum differences of 87 nM for October and 140 nM for June in surface waters) were similar to the diel changes in DFAAs (maximum-minimum difference 114 nM) in the Baltic Sea [14].
Water 2021, 13, x FOR PEER REVIEW 11 of 16 riod (Figure 2a,b). For molecular indicators, a comparison of the short-term changes (expressed as standard deviations) and the corresponding seasonal propagated errors indicated that short-term changes were significant for AA carbon yield, Gly%, and D/L Ala ratio, whereas the other two molecular indicators were not (e.g., GABA%; Figure 5). The variability (as a standard deviation) in GABA% over 24 h was similar to, or even smaller than, the propagated error ( Figure 5). The bulk DOC and THDAA changes over 24 h were clearly greater in June 2015 than in October 2014 (Figure 2a,b): for surface waters in June 2015, the standard deviations were 5 μM (DOC) and 57 nM (THDAAs), whereas in October 2014, they dropped to 3 μM (DOC) and 37 nM (THDAAs) ( Table 4). The short-term changes in THDAAs (maximumminimum differences of 87 nM for October and 140 nM for June in surface waters) were similar to the diel changes in DFAAs (maximum-minimum difference 114 nM) in the Baltic Sea [14]. In addition to their heterogeneous distribution, the short-term changes in DOM over 24h above the euphotic zone may be related to light-driven autotrophic contributions and related heterotrophic processes [13], whereas the short-term changes below the euphotic zone (i.e., 200 m in this study) may be attributable to either the heterogeneous distribution of microbes (and hence the DOM concerned) or the contribution of swimmers (vertically migrating zooplankton and fish). It has been suggested that the vertical migration of mesopelagic fish (and their excretions) plays a role in the diel variability in DOM (DOC and carbohydrates) during our sampling period [35]. The bulk DOC concentration and AAs showed a gradual decrease from the surface to deeper layers in June, but the statistical analysis of the respective diel observation data set-namely in June and October, respectively-indicates there are no seasonal relations between water masses (temperature, salinity, and density) and AA-based molecular indicators. Statistically, the bacterial abundance in June showed a negative relation with DOC, suggesting the potential influence of bacteria on DOM short-term changes. However, the relation is very weak (r = −0.55, respectively; figure not shown). Bacteria are periodically active in surface waters, even though their abundances are temporally similar [9]. Hence, it is reasonable that bacteria are playing an important role in constraining the short-term changes of DOM, but more sensitive biological proxies (e.g., ATP, RNA, and sequencing) are needed to quantify the biological role in DOM short-term changes. The mechanism of the short-term changes in DOM composition need further study.
The AA carbon yield, D/L ratio, DI, and nonprotein AA% have been intensively discussed previously, and have shown different patterns of change during DOM diagenesis [21,[36][37][38]. The AA carbon yield and the D/L ratio are most sensitive for tracking changes during the early stage of diagenesis (days to weeks) [20,36,38]. DI appears to be most effective during the intermediate stage of diagenesis (annual to decadal timescales), and nonprotein AA% (e.g., GABA% in this study) is the best indicator in the later stages of DOM diagenesis (e.g., >decades) [20]. Our results support these findings ( Figure 5) in that the daily changes in the AA carbon yield were apparent. The diel changes in D/L Ala were In addition to their heterogeneous distribution, the short-term changes in DOM over 24 h above the euphotic zone may be related to light-driven autotrophic contributions and related heterotrophic processes [13], whereas the short-term changes below the euphotic zone (i.e., 200 m in this study) may be attributable to either the heterogeneous distribution of microbes (and hence the DOM concerned) or the contribution of swimmers (vertically migrating zooplankton and fish). It has been suggested that the vertical migration of mesopelagic fish (and their excretions) plays a role in the diel variability in DOM (DOC and carbohydrates) during our sampling period [35]. The bulk DOC concentration and AAs showed a gradual decrease from the surface to deeper layers in June, but the statistical analysis of the respective diel observation data set-namely in June and October, respectively-indicates there are no seasonal relations between water masses (temperature, salinity, and density) and AA-based molecular indicators. Statistically, the bacterial abundance in June showed a negative relation with DOC, suggesting the potential influence of bacteria on DOM short-term changes. However, the relation is very weak (r = −0.55, respectively; figure not shown). Bacteria are periodically active in surface waters, even though their abundances are temporally similar [9]. Hence, it is reasonable that bacteria are playing an important role in constraining the short-term changes of DOM, but more sensitive biological proxies (e.g., ATP, RNA, and sequencing) are needed to quantify the biological role in DOM short-term changes. The mechanism of the short-term changes in DOM composition need further study.
The AA carbon yield, D/L ratio, DI, and nonprotein AA% have been intensively discussed previously, and have shown different patterns of change during DOM diagenesis [21,[36][37][38]. The AA carbon yield and the D/L ratio are most sensitive for tracking changes during the early stage of diagenesis (days to weeks) [20,37,38]. DI appears to be most effective during the intermediate stage of diagenesis (annual to decadal timescales), and nonprotein AA% (e.g., GABA% in this study) is the best indicator in the later stages of DOM diagenesis (e.g., >decades) [20]. Our results support these findings ( Figure 5) in that the daily changes in the AA carbon yield were apparent. The diel changes in D/L Ala were similar to (or seasonally larger than) the propagated error. Neither DI nor GABA% showed significant short-term changes compared with the propagated error.

Seasonal DOM Composition Difference and Its Implications
In the northern slope region of the SCS, Gly% in DOM is comparable to DOM samples from other regions in the northern SCS, like the offshore of the Pearl River Estuary [39], and from oligotrophic surface oceanic waters [37] (Figure 6). Diatom frustules and bacterial cell membrane macromolecules are the two DOM sources that are relatively enriched in Gly in the open sea [6]. In addition to Gly, bacterial cell membrane macromolecules contain a few other D-form AAs, whereas diatom frustules are enriched in L-form Ser and Thr [40,41]. During our diel observations, the elevated Ser% (mainly contributed by L-Ser) was as high as 10-23%, which was much higher than at BATS or HOT (3-12%) or in other oceanic DOM ( Figure 6; or see Table 4 in Jørgensen et al., 2014) [37]. However, our observed value (Ser%, 10-23%; Figure 6) was very similar to previously reported Ser% (9-21%) in the northern SCS (e.g., offshore waters from the Peal River Estuary) [39], suggesting that the elevated Ser% in the surface of the northern SCS DOM was a reproducible feature.
Water 2021, 13, x FOR PEER REVIEW 12 of 16 similar to (or seasonally larger than) the propagated error. Neither DI nor GABA% showed significant short-term changes compared with the propagated error.

Seasonal DOM Composition Difference and its Implications
In the northern slope region of the SCS, Gly% in DOM is comparable to DOM samples from other regions in the northern SCS, like the offshore of the Pearl River Estuary [39], and from oligotrophic surface oceanic waters [37] (Figure 6). Diatom frustules and bacterial cell membrane macromolecules are the two DOM sources that are relatively enriched in Gly in the open sea [6]. In addition to Gly, bacterial cell membrane macromolecules contain a few other D-form AAs, whereas diatom frustules are enriched in L-form Ser and Thr [40,41]. During our diel observations, the elevated Ser% (mainly contributed by L-Ser) was as high as 10-23%, which was much higher than at BATS or HOT (3-12%) or in other oceanic DOM (Figure 6; or see Table 4 in Jørgensen et al., 2014) [37]. However, our observed value (Ser%, 10-23%; Figure 6) was very similar to previously reported Ser% (9-21%) in the northern SCS (e.g., offshore waters from the Peal River Estuary) [39], suggesting that the elevated Ser% in the surface of the northern SCS DOM was a reproducible feature.  [39]; HOT and BATS data are from Kaiser and Benner (2009) [21]. October 2014 and June 2015: this study. Note that the time-series observations are shown as pooled dots with error bar (stdev). D forms of AAs have been widely used as an indicator of bacterial contribution to organic matter [38,[41][42][43]. The DOM samples had higher D/L Ala ratios in October than in June (Table 4, 5), which is consistent with the findings that bacteria were more abundant in October than in June (Table 3). As aforementioned, serine is contained in both diatom frustules (as the L-form) and bacteria membrane detritus (as D-Ser). For the TS2014 and TS2015, which showed contrasting Ser% values (Figure 6), the Ser D-and L-forms were 8.6 ± 3.3 nM (5.2 ± 1.5 nM) and 34 ± 9.3 nM (84 ± 20 nM) in TS2014 (TS2015), respectively. The D/L-Ser ratio in October (TS2014) was 0.25 ± 0.06, much higher than that in June (TS2015; 0.06 ± 0.009). Differences between the concentrations of Ser enantiomers, as well as the corresponding D/L ratios, suggest a stronger contribution from phytoplankton (e.g., diatom frustules) than from bacteria at TS2015 than at TS2014, which is consistent with the seasonal bacterial and phytoplankton results (Table 3). We should also bear in mind that the serine enrichment in phytoplankton detritus [40,41] is based on a particulate approach D forms of AAs have been widely used as an indicator of bacterial contribution to organic matter [38,[41][42][43]. The DOM samples had higher D/L Ala ratios in October than in June (Tables 4 and 5), which is consistent with the findings that bacteria were more abundant in October than in June ( Table 3). As aforementioned, serine is contained in both diatom frustules (as the L-form) and bacteria membrane detritus (as D-Ser). For the TS2014 and TS2015, which showed contrasting Ser% values (Figure 6), the Ser D-and L-forms were 8.6 ± 3.3 nM (5.2 ± 1.5 nM) and 34 ± 9.3 nM (84 ± 20 nM) in TS2014 (TS2015), respectively. The D/L-Ser ratio in October (TS2014) was 0.25 ± 0.06, much higher than that in June (TS2015; 0.06 ± 0.009). Differences between the concentrations of Ser enantiomers, as well as the corresponding D/L ratios, suggest a stronger contribution from phytoplankton (e.g., diatom frustules) than from bacteria at TS2015 than at TS2014, which is consistent with the seasonal bacterial and phytoplankton results (Table 3). We should also bear in mind that the serine enrichment in phytoplankton detritus [40,41] is based on a particulate approach study. Organic matter changes during particulate-dissolved phase transformation and its impact on the serine enrichment pattern are unclear and need further study.
In the oceans, the presence and enrichment of non-protein AA GABA indicates the aging and advanced diagenesis status of the DOM [44]. Moreover, due to the contribution from bacteria in DOM diagenesis, the D/L ratio of DOM is also expected to increase along with the early diagenesis progress. Though with an outlier, of all the time-series and section samples, the relation between GABA% and D/L Ala is statistically negative (r = 0.32, p = 0.005, n = 75). When GABA% was plotted against the D/L Ala ratio, samples from TS2015 showed the lowest GABA% and D/L Ala ratio relative to all other samples, being located in the left bottom corner (Figure 7). The lower GABA% and higher D/L Ala ratio of TS2015 is consistent with its corresponding higher Ser% content ( Figure 6). Both figures (Figures 6 and 7) support the idea that the DOM of samples from TS2015 was likely affected by phytoplankton instead of by bacterial rework.
Water 2021, 13, x FOR PEER REVIEW 13 of 16 study. Organic matter changes during particulate-dissolved phase transformation and its impact on the serine enrichment pattern are unclear and need further study. In the oceans, the presence and enrichment of non-protein AA GABA indicates the aging and advanced diagenesis status of the DOM [44]. Moreover, due to the contribution from bacteria in DOM diagenesis, the D/L ratio of DOM is also expected to increase along with the early diagenesis progress. Though with an outlier, of all the time-series and section samples, the relation between GABA% and D/L Ala is statistically negative (r = 0.32, p = 0.005, n = 75). When GABA% was plotted against the D/L Ala ratio, samples from TS2015 showed the lowest GABA% and D/L Ala ratio relative to all other samples, being located in the left bottom corner (Figure 7). The lower GABA% and higher D/L Ala ratio of TS2015 is consistent with its corresponding higher Ser% content ( Figure 6). Both figures (Figures 6,7) support the idea that the DOM of samples from TS2015 was likely affected by phytoplankton instead of by bacterial rework. When comparing both section and time-series observations, DOM in surface waters in June 2015 showed less GABA% and D/L Ala ratio when compared to that in October, 2014 (Figure 7), suggesting that DOM in June 2015 was overall more fresh (less degraded) when compared to that in October 2014. We note that such seasonal DOM degradation pattern is consistent with the seasonal bacterial and phytoplankton biomass pattern ( Table  3). The contribution of the eddy-pair in June during our observation is highly suspicious. In June, our section and time-series observations were under a mild upwelling effect, which can be judged from both the sea level anomaly (Figure 1) and the section density profiles ( Figure 4). Accordingly, higher chlorophyll was found in June relative to October (Table 3). Though the upwelling of deeper waters also brings up highly degraded DOM from the deep, the upwelling-promoted algae growth and the addition of newly assimilated organic matter into the surface DOM pool seem to play a more important role in determining the overall apparent DOM degradation status in the upper ocean.

Conclusion and Perspective
The THDAAs and their related molecular indicators reveal a bulk DOM diagenesis status that can correspond to thousands of years, e.g., [21]. Against this background, shortterm changes in THDAAs and AA-based molecular indicators, if they exist, would interfere with our understanding of long-term carbon cycling. Hence, the study of THDAAs and their molecular indicators' short-term changes has a potential significance on longtime scale materials cycling. When comparing both section and time-series observations, DOM in surface waters in June 2015 showed less GABA% and D/L Ala ratio when compared to that in October, 2014 (Figure 7), suggesting that DOM in June 2015 was overall more fresh (less degraded) when compared to that in October 2014. We note that such seasonal DOM degradation pattern is consistent with the seasonal bacterial and phytoplankton biomass pattern ( Table 3). The contribution of the eddy-pair in June during our observation is highly suspicious. In June, our section and time-series observations were under a mild upwelling effect, which can be judged from both the sea level anomaly (Figure 1) and the section density profiles ( Figure 4). Accordingly, higher chlorophyll was found in June relative to October (Table 3). Though the upwelling of deeper waters also brings up highly degraded DOM from the deep, the upwelling-promoted algae growth and the addition of newly assimilated organic matter into the surface DOM pool seem to play a more important role in determining the overall apparent DOM degradation status in the upper ocean.

Conclusion and Perspective
The THDAAs and their related molecular indicators reveal a bulk DOM diagenesis status that can correspond to thousands of years, e.g., [21]. Against this background, shortterm changes in THDAAs and AA-based molecular indicators, if they exist, would interfere with our understanding of long-term carbon cycling. Hence, the study of THDAAs and their molecular indicators' short-term changes has a potential significance on long-time scale materials cycling.
Overall, our work indicates that care should be taken when using THDAAs and their related molecular indicators to determine the DOM's degradation status. This is because some THDAA molecular indicators showed solid short-term changes within 1 day. The scale of variability within 1 day can be viewed as a resolution of the corresponding molecular indicators when elucidating the DOM composition changes over longer time scales. Such molecular indicators include the AA carbon yield and D/L Ala, which show 0.3-0.4% and 0.02-0.13 short-term changes, respectively.
Given these revealed short-term changes, we applied our findings to DOM seasonal composition comparison. We conclude that DOM in June 2015 is more fresh when compared to that in October 2014. This is quantified by the higher Ser% and lower GABA% and D/L Ala in DOM's composition in June relative to October. The seasonal differences of these molecular indicators are larger than their corresponding short-term changes. The reason for such DOM composition seasonality is the fact that there was an eddy-pair event going on during our observation in June 2015. Accordingly, more upwelled deeper waters, better nutrient concentration, and higher phytoplankton standing stock are reported in June relative to October. Our finding of such seasonal DOM composition differences, confirms the seasonal DOM's lability difference, which sheds light on the microbiol food loop and marine carbon cycling studies in the SCS.