Substantial Biogeochemical and Biomolecular Processing of Dissolved Organic Matter in an Anticyclonic Eddy in the Northern South China Sea Down to Bathypelagic Depths

Solid-phase extracted dissolved organic matter (SPE-DOM) was isolated from two depth profiles at the core and at the edge of an anticyclonic eddy (ACE) in the northern South China Sea. Non-target nuclear magnetic resonance (NMR) spectroscopy and Fourier transform ion cyclotron mass spectrometry (FTICR/MS) of SPE-DOM revealed a higher uniformity of DOM molecules within the ACE than at the edge of the ACE. Small-scale upwelling of external nutrients may have contributed to higher productivity and production of fresher DOM, with higher proportions of CHNO and CHNOS compounds and low molecular weight species at the edge of the eddy. Common SPE-DOM molecules of supposedly biological origin such as carbohydrates and olefins were most abundant in the chlorophyll maximum layer in both stations. An unusual suite of ~10 abundant and ~35 less abundant tert-butyl benzene derivatives with potential to act as endocrine disruptors within a marine food chain and ~two dozen ketones of putative bacterial origin was recognized at meso- and bathypelagic depths in single-digit micromolar concentrations, with a distinct maximum at 1000 m depth at the edge of ACE. Downwelling might bring temporary large volumes of productive marine waters into deep waters, with micromolar concentration of abundant, microbial food web-specific metabolites (e.g. 2,4-di-tert-butylphenol et al.). In our study, these eventually added up to one quarter of common background biogeochemical marine organic matter even at bathypelagic depths and beneath and are significant food and energy sources for marine biota. Mesoscale chemical heterogeneity of marine water columns might extend to larger depths than currently anticipated and may create activity hotspots influencing biota, processing of DOM, and cycling of nutrients and trace elements.

Solid-phase extracted dissolved organic matter (SPE-DOM) was isolated from two depth profiles at the core and at the edge of an anticyclonic eddy (ACE) in the northern South China Sea. Non-target nuclear magnetic resonance (NMR) spectroscopy and Fourier transform ion cyclotron mass spectrometry (FTICR/MS) of SPE-DOM revealed a higher uniformity of DOM molecules within the ACE than at the edge of the ACE. Small-scale upwelling of external nutrients may have contributed to higher productivity and production of fresher DOM, with higher proportions of CHNO and CHNOS compounds and low molecular weight species at the edge of the eddy. Common SPE-DOM molecules of supposedly biological origin such as carbohydrates and olefins were most abundant in the chlorophyll maximum layer in both stations. An unusual suite of~10 abundant and~35 less abundant tert-butyl benzene derivatives with potential to act as endocrine disruptors within a marine food chain and~two dozen ketones of putative bacterial origin was recognized at meso-and bathypelagic depths in single-digit micromolar concentrations, with a distinct maximum at 1000 m depth at the edge of ACE. Downwelling might bring temporary large volumes of productive marine waters into deep waters, with micromolar concentration of abundant, microbial food web-specific metabolites (e.g. 2,4-di-tertbutylphenol et al.). In our study, these eventually added up to one quarter of common background biogeochemical marine organic matter even at bathypelagic depths and beneath and are significant food and energy sources for marine biota. Mesoscale chemical heterogeneity of marine water columns might extend to larger depths than currently anticipated and may create activity hotspots influencing biota, processing of DOM, and cycling of nutrients and trace elements.

INTRODUCTION
Marine dissolved organic matter (DOM) represents one of the largest pools of organic carbon on earth, comprises more than 90% of organic carbon in the ocean, and plays an essential role in the global carbon cycle. Molecular characterization of DOM will provide relevant information about its precursors (Hertkorn et al., 2013;Igarza et al., 2019), the biological and biochemical processes relating to DOM formation and degradation, including its reactive components (Gonsior et al., 2014;Padmavathi et al., 2014;Ali and Tremblay, 2019;Powers et al., 2019;Seidel et al., 2022), and might enable discrimination of biotic and abiotic processing of DOM in various ocean layers.
Ocean mesoscale eddies are ubiquitous in near coastal and open oceanic waters (McGillicuddy Jr, 2016), and strongly affect the transport of waters, abundance and biodiversity of marine organisms and marine DOM (Doblin et al., 2016;Sarma et al., 2019). Mesoscale eddies not only influence advective transport and vertical water mixing, but also affect biogeochemical cycles and the ocean biological pump (Garcon et al., 2001;McGillicuddy Jr, 2016, Siegel et al., 2011. In the northern hemisphere, most anticyclonic systems are oligotrophic, with nutrient-poor surface waters and restricted primary production. However, anticyclonic eddies (ACEs) are usually associated with relatively abundant microorganisms, an appreciable growth rate and DOC (dissolved organic carbon) accumulation (Bidigare et al., 2003;Ewart et al., 2008;Baltar et al., 2010), possibly triggered by downwelling in early-stage anticyclonic systems (Lasternas et al., 2012). DOC-rich surface waters and particulate organic carbon (POC) in ACEs were downwelled into deeper layers where extensive biogeochemical processing greatly affected its composition and the contributions of DOC to the oceanic carbon pump (Kolasinski et al., 2012;Lasternas et al., 2012). Intensification of ageostrophic secondary circulation produces small-scale hotspots of upwelling on the periphery of ACEs, which serve as a frontal zone between the eddy and the surrounding waters (Klein andLapeyre, 2009, Siegel et al., 2011). As a result, both episodic nutrients and deep phytoplankton are driven upwards into euphotic zone, leading to the stimulation of primary production (PP) at the rim of anticyclonic eddies (Samuelsen et al., 2012;Zhou et al., 2013;Wang et al., 2018).
The combination of optical spectroscopy with non-target nuclear magnetic resonance (NMR) spectroscopy and Fourier transform ion cyclotron resonance mass spectrometry (FTICR/ MS) allows for an in-depth assessment of molecular diversity in marine and other DOM, which range among the most complex mixtures of organic molecules on earth (Hertkorn et al., 2007). This includes FTICR/MS-based assessment of heteroatoms CHNOS, distinction of chemical environments of hydrogen and carbon by NMR spectroscopy, and classification of chromophoric dissolved organic matter (CDOM) (Hertkorn et al., 2013;Hertkorn et al., 2016;Gonsior et al., 2017). Spectroscopic characterization of marine CDOM demonstrated a marked diversity between DOM in different water masses and allowed to distinguish putative sources of DOM (Kellerman et al., 2015;Wagner et al., 2015;Gonsior et al., 2017;Martıńez-Peŕez et al., 2017;Brogi et al., 2018). So far, only limited studies have been carried out to decipher the chemical diversity of DOM in dynamic, upwelling and eddy marine systems; one example describes biotic and abiotic processing of DOM in the Peruvian upwelling system (Igarza et al., 2019).
The South China Sea (SCS) is the largest semi closed marginal sea between the Western Pacific and Eastern Indian Ocean, and mesoscale phenomena in SCS are ubiquitous (Xiang et al., 2016;Zhang et al., 2016;He et al., 2018), with a frequency of more than 30 per year including~48% ACEs (Xiu et al., 2010). Wang et al. (2003) reported that most of the mesoscale eddies were distributed in the northern and eastern SCS, which modulated biogeochemical processes of SCS to remarkable extents.
To assess the impact of ACE on DOM composition and structure in SCS, two depth profiles of SPE-DOM from the edge and within an ACE from the northern SCS were isolated ( Figure 1; Table S1), comprising samples obtained from surface (5 m), chlorophyll maximum zone (80-90 m), upper mesopelagic zone (200 m), mesopelagic zone (400-450 m), oxygen minimum zone (OMZ, 700-1000 m) with oxygen concentrations varying from 1.5-1.8 mg/L, deep layer (1400-1900m), with about 40-70% carbon yield. While effects of temporal isolation of water masses within ACE on microbial communities have been reported (Baltar et al., 2010, Fernańdez et al., 2008, we have determined bulk and spectroscopic parameters of these SPE-DOM samples with high spatial resolution to reveal whether temporal isolation of water masses within ACE will affect DOM composition and structure, in particular with regard to DOM at deeper waters. These data will allow a more comprehensive integration of mesoscale marine currents into the global carbon and other element cycles (e. g. N and S).

MATERIALS AND METHODS
Sampling DOM samples were collected during March 2017 in the northern slope of SCS ( Figure 1). An anticyclonic eddy born in January 2017, moving southwest was observed in the study area. Sampling points L05, L11 and L15 (from southwest to northeast) were across the eddy, with L05 and L15 on the edge and L11 in the center of the ACE (Figure 1). The downward contours of temperature, salinity and density (s t ) in L11 reflected the impact of downwelling of ACE in the upper layers ( Figure S1).
Seawater for each sample was retrieved with acid-cleaned 12 L Niskin bottles attached to the CTD rosette, which recorded temperature, salinity, s t , chlorophyll fluorescence. Dissolved oxygen (DO) was measured immediately after sample collection with a DO analyzer (Jenco, San Diego, California, USA, Model 9173), which was calibrated with traditional iodometric method (r 2 = 0.9, slope = 0.97). Apparent oxygen utilization (AOU) was computed as the difference between oxygen saturation and the observed oxygen concentration (Garcia et al., 1992).
Seawater for DOM extraction was acidified to pH = 2 with HCl (Emsure, Merck) after filtration with GF/F glass filters (Whatman, 0.7 mm pore size; precombusted at 450°C for 5 h). A solid phase extraction (SPE) procedure was used to isolate DOM from the water with PPL cartridges (1 g, Agilent Bond Elut) as previously described (Dittmar et al., 2008). After conditioning of the cartridges with 10 mL methanol (GC-grade, Sigma-Aldrich) and Milli-Q water (pH = 2), 10 L acidified seawater for each sample were extracted with a cartridge. For removing the salt effect, the cartridges were rinsed with more than 20 mL Milli-Q water (pH = 2) and stored at -20°C in the dark until analyses. Before measurement, the cartridges were dried by flushing with high purity N 2 and eluted with methanol (GC-grade, Sigma-Aldrich). The recovery of SPE-DOM ranged from 40-70% (Table S1).

NMR Spectroscopy
For NMR spectroscopy, aliquots of original methanol solutions of marine SPE-DOM samples were evaporated in vacuo to dryness and deuterated methanol (CD 3 OD) was added; this cycle was repeated three times to largely exchange methanol-h 4 with methanol-d 4 . Here,~1-3.5 mg yellow solid marine SPE-DOM was obtained. The 1 H NMR spectra of all the South China Sea SPE-DOM samples were acquired with a Bruker Advance III NMR spectrometer operating at 800.35 MHz (B 0 = 18.8 Tesla) and Topspin 3.5.7 software with samples from redissolved solids in sealed 2.0 mm Bruker Match tubes. Proton detected NMR spectra were acquired at 283 K with an inverse geometry 5 mm zgradient 1 H/ 13 C/ 15 N/ 31 P QCI cryogenic probe (quaternary cryogenic inverse; 90°excitation pulses: 13 C~1 H~10 ms). NMR chemical reference: 1 H/ 13 C NMR, CD 3 OD: 3.30/49 ppm. 1-D 1 H NMR spectra were recorded with a spin-echo sequence (10 ms delay) to allow for high-Q probe ring-down (Q: high quality factor), and classical presaturation to attenuate residual water present: noesypr1d, typically 512-2048 scans (5 s acquisition time, 5 s relaxation delay (d1), 1 ms mixing time; 1 Hz exponential line broadening). Absolute value J-resolved (JRES), COSY and phase-sensitive TOCSY/NOESY NMR spectra (correlation spectroscopy; total correlated spectroscopy, nuclear Overhauser spectroscopy) used acquisition times of 1 s at a spectral width of 9615 Hz, with 50 Hz sweep with in F1 (JRES), 70 ms dipsi-2 mixing sequence (TOCSY with solvent suppression: dipsi2etgpsi19); NOESY mixing time 350 ms (noesygpphpp). Carbon decoupled 1 H, 13 C HSQC spectra (heteronuclear single quantum coherence; hsqcedetgpsisp2.2) were acquired under the following conditions: acquisition time: 250 ms; 13 C 90°decoupling pulse, GARP (70 µs); 50 kHz WURST 180°1 3 C inversion pulse (wideband, uniform, rate, and smooth truncation; 1.2 ms); F2 ( 1 H): spectral width (SW) of 9615 Hz, 1 J CH = 150 Hz, 1.25 s relaxation delay; F1 ( 13 C): SW = 33784 Hz (167.9 ppm); number of scans (F2)/F1-increments ( 13 C frequency) for SPE-DOM (400-1800/137-142). HMBC s p e c t r a ( h e t e r o n u c l e a r m u l t i p l e b o n d c o h e r e n c e ; hmbcacbigplnd) were acquired with aq = 1 s, d1 = 0.5 s, F1 ( 13 C): SW = 52083 Hz (258.8 ppm). 13 C NMR spectra were acquired with a Bruker Advance III NMR spectrometer operating at 500.13 MHz (Bo = 11.7 T) and used inverse gated WALTZ-16 decoupling (wideband alternating-phase low-power technique for zero residual splitting; 19 s relaxation delay; with an acquisition time (aq) of 1 s; 13 C DEPT NMR spectra used aq = 1 s and d1 = 2 s; the one bond coupling constant 1J(CH) was 150 Hz; exponential line broadening was lb = 12.5 Hz. Other acquisition conditions are provided in Table S2. Compound M2* (3-(3,5-Di-tert-butylphenyl) propionic acid was purchased from Abcr GmbH, Karlsruhe, Germany, and dissolved as provided in CD 3 OD for acquisition of NMR spectra.

FTICR/MS Analysis
Negative electrospray ionization [(-)ESI] Fourier transform ion cyclotron resonance mass spectra (FTICR/MS) were acquired using a 12 T Bruker Solarix mass spectrometer (Bruker Daltonics, Bremen, Germany) and an Apollo II electrospray ionization source in negative mode. Nebulizer gas pressure, drying gas pressure and the source heater temperature were 138 kPa, 103 kPa and 200°C, respectively. The spectra were acquired with a time domain of 4 MW and 300 scans were accumulated for each mass spectrum. Computation of elemental formulas for each peak was done in a batch mode by an in-house written software tool. The generated formulae were validated by setting sensible chemical constraints (N rule, O/C ratio ≤ 1, H/C ratio ≤ 2n + 2 (C n H 2n+2 ), element counts: C ≤ 100, H ≤ 200, O ≤ 80, N ≤ 3, S ≤ 2, P < 1 and mass accuracy window, set at ± 0.5 ppm). Final formulae were generated and categorized into groups containing CHO, CHNO, CHOS or CHNOS molecular compositions, which were used to reconstruct the group-selective mass spectra (Schmitt-Kopplin et al., 2010a;Schmitt-Kopplin et al., 2010b). The computed average values for H, C, N, O and S (atom %), elemental ratios, DBE-related parameters and average mass (m/z) given in Table 1 were based upon intensity-weighted averages of area-normalized mass spectra using assigned molecular formulae.

Statistical Analysis
Hierarchical Cluster Analysis (HCA) was performed based on 0.01 ppm buckets of 1 H NMR spectra (from 0.5-10 ppm, with exclusion of HDO and CD 2 HOD) and assigned, intensityweighted mass peaks (n = 7960) in FTICR mass spectra, using the Hierarchical Clustering Explorer 3.0 (HCE; http://www.cs. umd.edu/hcil/multi-cluster/), with clustering of the dataset using the Average Linkage (UPGMA) method and the Pearson correlation coefficient (Pearson's R) as the similarity/distance measure. Based on the HCA, we used the "profile search" tool from HCE 3.0, choosing a search method (model-based), a distance measure (Pearson's R) and a threshold (0.9). Principal Component Analysis (PCA) was done with the software SIMCA-P 9.0 (Umetrics AB, Sweden). FTICR mass spectra were arranged with the samples as observations and the peak areas of the assigned FTICR mass peaks as the response variables. 1 H NMR spectra were arranged with the samples as observations and the NMR resonances as the response variables (800 MHz 1 H NMR, CD 3 OD, area-normalized from 0.5-10.0 ppm; 0.01 ppm bucket resolution; with the exclusion of residual water and methanol NMR resonances). Before multivariate statistics were performed, the response variables were centered and scaled to unit variance. The based weight was computed as 1/sqrt (standard deviation of the response variables). SPSS 23.0 was used to compute the correlation between assigned molecular series and optical parameters (Spearman's R) within 0.95 confidence limit.

Chemical Diversity of SPE-DOM in SCS
The proton NMR spectra of the SPE-DOM samples in SCS showed smooth, broad bulk envelopes familiar from highly processed biogeochemical marine DOM that represent superpositions of millions of atomic environments grouped in key substructures (Figures 2, S2-S4) (Hertkorn et al., 2013). In general, the abundance of sharp resonances in 1 H NMR spectra declined in the order of   (Figures 2, S2, S3). Massive projection of millions of individual atomic environments of 1 H atoms produces sweeping intrinsic averaging in NMR spectra of SCS organic matter (Hertkorn et al., 2007), and extensive molecular changes are required to produce visible differences of bulk NMR spectral features, which reflect multiple layers of biogeochemical processing in marine foodwebs. In comparison, the chemical diversity of biomolecules from known metabolic networks is more limited than that of biogeochemical molecules like DOM.
Recognizable biomolecules in marine SPE-DOM are typically much more abundant than (non-assigned and largely unknown) individual biogeochemical molecules (~10 -4 compared with ≤ 10 -7 of 1 H NMR integral); they share common structural features (Hoffmann et al., 2020) and project on narrow ranges of 1 H NMR chemical shift. Hence, even minor operating biochemical processes will produce readily recognizable features in 1 H NMR spectra of marine SPE-DOM. The differences between 1 H NMR spectra of each sampling depth indicated various physical and biochemical processes in different layers (Figures 2, S2-S4). While many methyl resonances terminating alkyl groups resonated near d H~0 .9 ppm, the resonances of upper layer (< 200 m) for L05 and L11 were more smooth than those of deep layers (> 400 m) (Figures 2-a2, b2). Similarly, the intensity maxima at d H~1 .3 ppm at deep water layers, representing methylene, alkyls and other purely aliphatics with HC ≥ 4 bonds away from the next heteroatoms, were higher than those in upper layers (Figures 2-a1, b1). Combined with the progressively increasing downfield shoulders ranging d H~1 .4-1.7 ppm, which represented certain remotely oxygenated aliphatics and alicyclic rings (Figures 2-a2, b2), higher proportions of aliphatic protons (d H~0 .5-1.9 ppm) were observed in the deep layers rather than in the upper layers ( Figure S4). The resonances coincident with acetyl and acetate derivatives (d H~2 .1 ppm) decreased from surface to OMZ, and then increased slightly (Figures 2-a5, b5). Functionalized aliphatics including those with remote carboxylic substitution typical of CRAM (heteroatoms ≥ 2 bonds away from protons; d H~2 .3-3.1 ppm) increased from surface to bottom, except SPE-DOM of L05_1000 m which showed comparatively low abundance of these units ( Figure S4). The superimposed small NMR resonances indicating the diversity of CRAM protons and related biomolecules were more significant in deep layer (> 400 m) (Figures 2-a5, b5). Isolated double bonds (d H 5.2-6.0 ppm) were characteristic of natural products (Pisani et al., 2018), and commonly declined from surface to depth, in the order of 80 m > 5 m > 200 m > 400 m ∼ 1898 m > 1000 m in L05 station and 5 m > 90 m > 200 m > 450 m > 700 m ∼1490m in L11 station ( Figures 2-a6, b6). Oxygenated aliphatic (OCH) units also decreased from surface to deep waters ( Figure S4) (Hertkorn et al., 2013;Benner and Amon, 2015).
The relative abundance of pure aliphatic protons (HC-C-C-; d H~0 .5-1.9 ppm) in marine DOM always exceeded 40%, with higher abundance in deep (~48%) than in surface (~42%) waters ( Figure S4). In both L05 and L11, the relative abundance of the functionalized protons (HC-N, HC-C-O; d H~1 .9-3.1 ppm) increased from surface to bottom layers. In OMZ (1000 m of L05, 700 m of L11), however, the relative abundance of functionalized protons decreased slightly. The abundance of carbohydrates, alcohols, esters, ethers and other singly oxygenated units (HCO; d H~3 .1-5.0 ppm) at L05 and L11 stations declined from surface to deep ( Figure S4), and more strongly in the upper 200 m, demonstrating that oxygenated aliphatic molecules were produced primarily within the sunlit layer of the surface ocean and consumed primarily within the upper mesopelagic zone (Benner and Amon, 2015).
The total contributions of aliphatic HC-C-C-, HC-N, HC-C-O and HCO units were higher than 90% in marine SPE-DOM. Olefins (HC=C-; d H~5 .0-7.0 ppm) and aromatics (H ar ; d H7 .0-10.0 ppm) contributed to the unsaturated protons, with <10% of total 1 H NMR integral and slight variance in the water profile ( Figure S4). The relative abundance of olefinic protons in upper layers was higher than that in the deep layers, in line with the supposed biological origin of olefins (HC = C, some anomeric positions HCO 2 of supposedly biotic origin also resonated here) (Hertkorn et al., 2013;Pisani et al., 2018), ( Figure S4). In the OMZ of L05, however, the olefins were also abundant ( Figure 2).
The content of aromatic proton (C ar H) was highest in the OMZ, which could result from the decreased activity of microbes when lacking oxygen. The ratio between olefinic and aromatic protons (H olefinic /H aromatic ) varied between 1.7 and 2.8, indicating that the relative abundance of olefinic protons exceeded that of aromatic protons in marine SPE-DOM, in accordance with previous results (Hertkorn et al., 2013). The ratio H olefinic /H aromatic was highest in the chlorophyll maximum layer and lowest in the OMZ ( Figure S4), indicative of a production of olefins and aromatic substances by different biochemical processes in marine water profiles. Ratios of aliphatic to aromatic protons (H aliphatic /H aromatic ) were lowest in OMZ of both stations, suggesting the decline of microbialderived aliphatic rich DOC (Hertkorn et al., 2016).
FTICR mass spectrometry with negative electrospray ionization of SCS SPE-DOM provided several thousands of mass peaks ( Figure S5), of which many hundreds could be assigned to extended CHO, CHOS, CHNO, and CHNOS molecular series (Table 1; Figure S5) based on excellent mass accuracy and mass resolution Schmitt-Kopplin et al., 2010b). In the upper 200 m, higher numbers of molecular formulas were assigned in L11 than in L05, while more mass peaks were assigned in the deep layer for L05 than for L11 samples ( Table 1).
Visual inspection of the spectra showed smooth, slightly skewed near Gaussian distributions of mass peaks, with prominent mass spacings caused by ± CH 2 (Dm = 14.0565 Da) and ± H 2 (Dm = 2.00121 Da), accompanied by sets of abundant low mass peaks with declining relative abundance in the order L11_700 m ≈ L05_1000 m > L11_1490 m ≈ L05_200 m > L11_90 m ≈ L11_5 m ( Figure S5); these irregular series often represent relatively abundant biomolecules, but also impurities like pollutants or surfactants might have contributed. No clear mass shift of the most abundant mass peaks was observed, but a relative depletion of low mass peaks from top to bottom layers in L05 and less so in L11 took place (m/z < 350; Figure S5), reflecting progressive degradation of low m/z molecules along the depth profile (Flerus et al., 2012;Hertkorn et al., 2013). The OMZ at L05_1000 m and L11_700 m showed the lowest average mass (Table 1; Figure S5).
The FTICR/MS-derived relative weight percentage of oxygen in SPE-DOM decreased slightly from surface to deep layers in station L05 and showed no obvious trend in station L11 (Table 1). Meanwhile, the weight percentage of carbon increased from surface to deep, and hence, the computed O/C ratio decreased steadily from upper to deep layers, probably caused by a combination of abiotic (e.g. redox chemistry, photochemistry, mineral association) and biotic processes (e.g. metabolism of organisms in foodwebs, lysis and exudation) subject to upwelling/downwelling propargation, resulting in less oxygenated (i.e. defunctionalized) organic matter following loss of H 2 O and CO 2 upon processing, which then might become less degradable and acts as carbon sink in the ocean (Table 1) (Nebbioso and Piccolo, 2013). The average H/C ratio very slightly decreased from surface (~1.23) to bottom waters (~1.20) in L05 and L11. The total weight percentage of assigned nitrogen and sulfur in organic compounds was < 1% in the water profile. CHO-molecular series comprised > 40% in all SCS SPE-DOM, with a relative maximum at the OMZ for both stations (44.4% of L05 and 43.7% of L11, respectively), and the count of CHO compounds increased from surface to bottom in L05 but decreased in L11. CHNO-molecular series were the second most abundant compounds in SCS SPE-DOM, with contributions of 35-41% to total assigned molecules. The relative proportions of both CHO and CHNO molecular series increased with depth in L05 and decreased in L11, and the highest percentage of CHNO compounds was observed in the chlorophyll maximum zone of both stations. The contributions of assigned CHOS-and CHNOS-molecular series were about 20% with less variation in the water profiles. Double Bond Equivalent (DBE), an empirical criterion to provide counts of C=C and C=O double bonds as well as numbers of alicyclic rings present in molecular compositions and DBE per carbon (DBE/C) were proposed for identification of (condensed) aromatic structures from FTICR mass spectra (Koch and Dittmar, 2006). DBE/C in SPE-DOM increased from surface to deep (Table 1), in line with previous studies (Benner et al., 1992;Flerus et al., 2012;Repeta, 2015), reflecting a progressive increase of relative unsaturation in deep water DOM, which may impart some resistance to microbial degradation of "residual" DOM.
The highly unsaturated compounds (HU; for definition, see methods) contributed most to the SCS SPE-DOM (86.2%) ( Figure S6; Table S3), and was higher than in the Pacific ocean (Igarza et al., 2019). HU showed highest average molecular weight and largely comprised CHO compounds. The relative abundance of HU increased from surface to deep waters, but was lowest in the OMZ ( Table S3). The unsaturated compounds (UC) represented 9.5% of SPE-DOM, including 2.4% of peptide-like molecular series. Although no sulfur or nitrogen atoms were present in more than 74% of UC, the percentage of S-containing molecular series was higher than that in HU ( Figure S6). The UC molecular series decreased with depth except at the OMZ, in which more UC were observed. UC in the chlorophyll a maximum zone was slightly more abundant in L11 (11.7%) than in L05 (10.5%) ( Table S3). The average molecular weight of CHNO compounds (418 Da) was higher than other UC, which decreased from surface to deep. Meanwhile, more CHNO compounds were present in L11 chlorophyll a maximum zone than at L05 ( Table S3). The MSbased percentages of condensed compounds (CC), saturated compounds (SC) and sugars-like molecular series were less than 1% respectively ( Figure S6). More than half of condensed molecular series had no sulfur atom, while most contained one or more nitrogen atom ( Figure S6), also shown in van Krevelen diagrams ( Figure S7). The percentage of CC in chlorophyll a maximum zone of L05 was much lower than others, indicating that fresher DOM was produced ( Table S3). The average molecular weight of SC was lowest (343 Da), most of which contained one sulfur atom ( Figure S6). The distribution of SC in L05 and L11 was reversed (Table S3). About 80% of carbohydrate-like compounds had one or more nitrogen, only 20%, however, were sulfur-bearing ( Figure S6). In chlorophyll a maximum zone, sugars-like compounds were more abundant in L05 than L11, congruent with NMR results.
In L05 and L11, CHO and CHNO compounds formed contiguous lobes in van Krevelen and mass-edited H/C ratio diagrams, whereas CHOS compounds formed single contiguous lobes in mass-edited H/C diagrams and two distinct lobes in van Krevelen diagrams, indicative of two different classes of CHOS compounds being present in L05 and L11 SPE-DOM ( Figure  S7). The CHOS and CHNOS compounds with high H/C and O/ C ratios progressively decreased with water depth ( Figure S7). CHNOS compounds were more scattered in all these diagrams, with a minor contiguous lobe (H/C~1.0 -1.6; O/C~0.45-0.7) in van Krevelen diagrams of L11 ( Figure S7). Contiguous lobes are indicative of processed biogeochemical materials, with a steep decline of relative abundance at the edges that originate from intrinsic averaging in extremely complex mixtures of organic molecules (Hertkorn et al., 2007). In comparison, scatter typically results from a less deep coverage of the compositional space; this applies here as well because CHNOS compounds were much fewer than the other molecular series (Table 1). Considering, that entropy-driven biogeochemical processing of organic matter frequently results in continual distributions of bulk properties, with strong gradients toward the edge, the conspicuous scatter of CHNOS compounds in both van Krevelen and mass-edited H/C diagrams at first is in line with the lower relative counts and abundance (Table 1) but also testifies for significant contributions of biological processes in their formation (Hertkorn et al., 2016). The CHO compounds in FTICR/MS-derived van Krevelen diagrams of L05 and L11 were positioned in accordance with CRAM (O/C = 0.3-0.6, H/C = 0.7-1.6; (Hertkorn et al., 2006)) ( Figures S7-e1, e5). CHNO compounds occupied larger areas than CHO compounds in van Krevelen diagrams and mass-edited H/C ratios ( Figure S7), indicative of a larger chemical diversity (Andersson et al., 2019) for CHNO compounds of yet unknown provenance, with more saturated (H/C > 1.6) and more oxygenated (O/C > 0.8) molecules than observed in CHO compounds. Highly oxygenated CHNO compounds (O/C > 0.7), which were more abundant in the upper than in the deep layers, showed higher relative saturation (H/C > 1.2) than the other CHNO compounds and possibly comprised polyoxygenated aliphatic compounds with N-containing functional groups of yet unknown provenance. The assigned CHOS molecular series showed appreciable unsaturation, with H/C ratios ranging from 1.0-1.4 ( Figure S7), and less (O/C < 0.4) and more oxygenated (O/C > 0.4) compounds ( Figure S7) resulted from different biogeochemical processes. In the upper layers of L05 (5 m and 80 m), S-bearing polyoxygenated compounds were present ( Figures S7-b, S7-a2, b2, section b), which were not obvious in L11; these were probably related to primary production at the edge of ACE. Higher proportions of less oxygenated CHOS molecular series were observed in the deep layer of L05 ( Figure S7-f2, section c).

Spatial Discrimination of SPE-DOM in the ACE
Nutrients when supplied to the upper layers with the small-scale upwelling at the edge of ACE might stimulate primary production, with synthesis of rather fresh organic matter (Samuelsen et al., 2012;Karstensen et al., 2017). NMR and FTICR mass spectra were supportive of biotic process participating in the formation of DOM. When the 1 H NMR spectral intensities were scaled to 100% total integral within the entire region of chemical shift (d H : 10-0.5 ppm), the classical carbohydrates (d H~3 .4-4.1 ppm) including methoxy peaks (d H~3 .65 ppm) declined from upper to deep layers (Figures 2-a3, 2-b3). The bioavailable carbohydrates in the chlorophyll maximum zone at the edge of the ACE (L05_80 m) were more abundant than in the core (L11_90 m) ( Figure  S4). In line with the supposed biological origin of olefins (HC=C-) (Hertkorn et al., 2013;Pisani et al., 2018), the relative abundance of olefinic protons in L05_80 m was higher than in L11_90 m ( Figure S4). Here, lower nutrient availability at the core of the ACE could have attenuated photosynthetic and bacterial phytoplankton (Lasternas et al., 2013). In L11 station, the weight percentage of nitrogen (N%) decreased in the deep euphotic layer (~200 m) ( Table 1), probably related to the fast recycling and remineralization of N-rich DOM by heterotrophs (Kujawinski, 2011) in SCS, where phytoplankton growth was limited by nitrogen availability (Chen et al., 2007). An analogous deficit of nitrogen content was not observed in L05 station, which could be explained by the addition of nutrients from the smallscale hotspot upwelling parts at the edge of the ACE (Samuelsen et al., 2012). The increase of C/S ratio at the chlorophyll maximum zone was more apparent in L05 station (211) ( Table 1), indicating that depletion of sulfur was much stronger at the edge than at the core of the ACE (Zhou et al., 2013).
Principal component analysis (PCA) derived from 0.01 ppm buckets of 1 H NMR spectra indicated that SPE-DOM of L11 formed a more homogeneous group than SPE-DOM of L05 ( Figure S9A), suggesting that organic matter inside the eddy was chemically less diverse, likely resulting from a limited exchange of DOM in ACE with outside waters. PC1 represented more than half of influencing factors, with opposing horizontal ends indicating DOM in upper and deep layers, respectively. SPE-DOM in L11_450 m was more similar to DOM in upper layers than L05_400 m; possibly, downwelling in the ACE center had impact on DOM down to 450 m depth. The distribution of temperature, salinity and s t on section X also demonstrated a depth of eddy in excess of 450 m ( Figure S1). In addition, the distribution of EEM-derived fluorescence components (discussed in supplementary material; Figure S10) as well as other bulk parameters of DOM (e.g. DOC, SP/SH, A/C, M/C; discussed in supplementary material; Figures S11, S12, S13, S14; Tables S4, S5) on the section X provided properties of water masses within and at the edge of the ACE ( Figure S11). The ACEs in SCS influenced the waters down to more than 500 m (Xiang et al., 2016;He et al., 2018), and even more than 2000 m (Zhang et al., 2016). For the absence of sampling between 450 m and 700 m at L11, the exact depth that was visibly influenced by downwelling in the ACE center remained undefined.
NMR-derived hierarchical cluster analysis (HCA) produced two groups for both L05 and L11, comprising an upper layers group (≤ 200 m) and a deep layers group (≥ 400 m) (Figures S9b, S9c). Within the former, the SPE-DOM at 200 m depth was distinct from the surface samples in agreement with previous work (Hertkorn et al., 2013) that had shown maximum primary production at fluorescence maximum (L05_80 m/L11_90 m) whereas the surface sample (5 m) was affected by photodegradation.
FTICR/MS-derived HCA produced four groups in L05, discriminating upper layers (5 m and 80 m), mesopelagic zones (400 m and 1000 m), top mesopelagic layer (200 m) and deep layer (1898 m), whereas only two groups (upper layers ≤ 200 m and deep layers ≥ 450 m) were distinguished within L11 (Figures 3-a1, 3-a2). This indicated more distinctive differences of organic matter composition in the water profile at the edge of the eddy than at the center. FTICR/MS-derived HCA demonstrated that while in L05, small molecules (< 350 Da) were more abundant in the upper and in the mesopelagic layers (Figures 3-b2, 3-c2), large molecules (> 450 Da) were more abundant in the deep layers (Figure 3-e2). However, both small and large molecules were observed at 200 m depth (Figure 3-d2). The synchronicity of molecular distributions in analogous depth zones could be explained by analogous relationships of DOM production and degradation. In L11, comparative HCA-based analysis of mass spectra for surface and deep waters led to a remarkable discrimination: a tight cluster of rather unsaturated (H/C: 0.9-1.4) compounds from 300-470 Da was abundant at depth (Figure 3-g2) whereas more saturated compounds (H/C: 1.0-1.7) of considerable mass range from 230-830 Da were abundant in surface waters (Figure 3-f2). Both sets were mutually exclusive, producing a characteristic "dent" in the mass-edited H/C diagram of the surface SPE-DOM (Figure 3-f2). CRAM was commonly distributed from surface to deep (Figure 3), somewhat enriched in the deep layers (Figures 3-e3,  3-g3), indicating refractory DOM contributed more to total DOM in the deep layers, also in accordance with 1 H NMR spectra (Figure 2-a5, b5). Carbohydrates and tannin-like compounds were abundant in upper 200 m ( Figures. 3-b3, d3,  f3), while the less oxygenated CHOS molecular series were more abundant mainly in the mesopelagic and deep layers (Figure 3-e3).

Metabolites and Biogeochemical Organic Matter in ACE of SCS
All 1 H NMR spectra of L05 and L11 SCS SPE-DOM, including those in deep layers of water, showed a multitude of intense sharp NMR resonances, especially for C sp2 H units and within a narrow range of aliphatic protons from d H 0.9-1.5 ppm (Figures 2, 4, 5, S2, S3; cf. Chemical Diversity of SPE-DOM in SCS). These represented relatively abundant small molecules of putative biological (metabolites) or anthropogenic (pollutants) origin. Admitting a limited availability of meaningful NMR data (Hertkorn et al., 2013;Arakawa et al., 2017;Gonsior et al., 2017;Powers et al., 2019), 1 H NMR spectra of typical PPL-based SPE-DOM isolated from mesopelagic and lower oceanic waters commonly display common broad bulk envelopes with rather marginal fine structure. Attenuation of sharp NMR resonances in 1 H NMR spectra of methanolic solution of open ocean SPE-DOM arises from the presence of millions of atomic environments  in a highly processed biogeochemical material (cumulative projection of Lorenztian lines with finite line width onto NMR spectra with a limited overall bandwidth) which distribute in highly functionalized, oxygen-and often carboxyl-rich molecules (Hertkorn et al., 2006). This considerable functionalization facilitates mutual interactions between DOM molecules in solution, causing acceleration of transverse NMR relaxation and line broadening (Hertkorn, 2014), further lowering disposable S/N ratio and observability of concerned 1 H NMR resonances. In contrast, photochemical degradation of isolated Pacific and Atlantic Ocean surface water in the lab had produced a suite of~2% polyols within hours which showed sharp NMR resonances and distinctive cross peaks in 2D NMR spectra (Gonsior et al., 2014). In SCS SPE-DOM, the relative abundance of analogous sharp NMR resonances was particularly pronounced in SPE-DOM L05_1000 m, sufficient to produce a remarkable distinction in the visual appearance of its 1 H NMR spectrum (Figures 2, 5, S2, S3), in NMR-derived HCA, and, particularly strong, in PCA ( Figure S9). The large S/N ratio of the intense sharp NMR resonances in the 1D 1 H NMR spectrum of L05_1000 m relatively attenuated the (otherwise dominant) broad signature of common organic matter ( Figure S2). The 1 H NMR spectrum of SPE-DOM L05_1000 m showed an array of~15 very and~20 less intense singlet resonances at d H~0 .9-1.5 ppm (Figures 4, 5, S15A1, S16A1) and strong NMR resonances in the aromatic region from d H~6 .6-7.6 ppm (Figures 2, 5). The most abundant aromatic compound showed J-coupling patterns characteristic of a 1,2,4-substituted benzene of which the upfield proton at d H~6 .6 ppm showed a single vicinal coupling 3 J HH~7 .5 Hz (H6 in Figures 5, S15A1, S16A1), suggesting an adjacent oxygen atom in ortho-position. The FIGURE 3 | Comparative analysis of intensity weighted negative ESI FTICR mass spectra of marine SPE-DOM in stations L05 and L11 produced hierarchical cluster analysis (HCA) diagrams (a1, a2). From top to bottom rows: FTICR mass spectra from m/z 150-950 with the assigned molecular compositions and relative proportions of color coded CHO (blue), CHOS (green), CHNO (orange) and CHNOS (red) molecular series, shown as van Krevelen diagrams and mass edited H/C ratios of the relatively more abundant compositions found in near surface layers (L5_5m and L5_80m) (b1-b3), mesopelagic layers (L5_400m and L5_1000m) (c1-c3), top mesopelagic layer (L5_200m) (d1-d3), L05 deep (e1-e3), upper layers (L11_5m, L11_90m and L11_200m) (f1-f3), and deep layers (L11_450m, L11_700m and L11_1490m) SPE-DOM. Section a: CRAM; section b: Carbohydrates and tannin-like compounds; section c: Less oxygenated CHOS molecular series.

NMR integral of the large aliphatic singlet NMR resonances was
10-fold compared with that of the (single) aromatic protons H3, H5 and H6 (Figures 5, S15), which therefore represented tert-butyl (C 4 H 9 ) rather than methyl (CH 3 ) groups. This identification also conformed to the observed 1 H NMR chemical shift of the singlet NMR resonances at d H~1 .2 ± 0.3 ppm, because C ar -CH 3 groups would resonate at d H > 2.1 ppm. Furthermore, 2-D NMR spectra of L05_1000 m provided remarkable richness in detail and aided in improving preliminary assignment proposals from the onedimensional 1 H spectra. 1 H, 1 H -resolved (JRES) NMR spectra confirmed 1,2,4 substitution of benzene ( 3 J HH , 4 J HH ) and the singlet character of the main NMR resonances of aliphatic protons from d H 0.9-1.5 ppm ( Figure S15A). Strong cross peaks in the section from d H ∼0.9-1.5 ppm were observed in 1 H, 13 C HSQC NMR spectra at d C ∼ 26 ppm ( Figure 4A, red) and in 1 H, 13 C HMBC NMR spectra at d C ∼35 ppm (Figure 4, blue), again corroborating presence of tert-butyl groups in abundance. The 1 H, 13 C HMBC cross peak at d H/C 2.79 and 2.55 ppm/~175 ppm was assigned to a C ar -CH-CH-COOH unit (Figure 4). The combination of 1-D and 2-D NMR spectra showed the presence of two major molecules: 2,4-di-tertbutylphenol (M1; C 14 H 22 O; NMR data in accordance with (Yoon et al., 2006)) and 3,3'-(3-tert-butylbenzene-1,5-diyl) dipropanoic acid (M2; C 16 H 22 O 4 ), a not yet described compound; SciFinder July 2020) and their analogues, which comprised 15% of the total FIGURE 4 | 1 H NMR detected NMR spectra of CCCH 3 groups of L05_1000m SPE-DOM; (A): most intense NMR cross peaks indicative of tert-butyl containing abundant molecules, with cross peaks from 1 H, 13 C DEPT-HSQC NMR spectra in red and those of 1 H, 13 C HMBC in blue color; (B): 1 H, 13 C DEPT-HSQC (red: CH 3 ; green CH 2 ) and 1 H, 13 C HMBC (blue) NMR spectra of L05_1000m SPE-DOM, shown down to lower cross peak amplitude; the common, contiguous cross peaks of biogeochemical organic matter show up in this display labelled by ellipsoids (cf. Figures S15, S16);. 1 H, 13 C HMBC cross peaks (blue color) represent common aliphatic branching mofis of DOM and relatively abundant molecules with considerable extent of aliphatic branching as exemplified by d C ranging down to d C~5 5 ppm (CCCH units). (C): 1 H, 13 C HMBC cross peaks (blue color) demonstrating abundance of aromatic compounds (C ar H and C ar O units; d C~1 20-160 ppm), a few carboxylic acids (d C~1 65-185 ppm), and several ketones (d C~1 95-220 ppm). Some of these units could represent to benzylic units of general character C 6 H x Z-Ca(CH 3 ) n -Y 3-n in which Z would represent aromatic substitution and Y oxygenated aliphatic and aromatic units, including carboxylic acids and ketones. The benzylic position Ca would not carry directly bound protons, and tert. butyl C 4 H 9 then represents the most basic unit of such benzylic position. 13 C NMR spectra are 13 C single pulse (black) and CH 2selective (blue) 13 C DEPT NMR spectra.
1 H NMR integral in L05_1000 m ( Figure S17). The relative abundance of compound M1 declined in the order L05_1000 m > L05_400 m > L11_200 m > L5_5 m ≈ L11_450 m > L5_1898 m > L11_1898 m > others, whereas that of M2 declined in the order L05_1000 m > L11_450 m > L5_400 m > L5_1898 m ≈ L11_700 m > others ( Figure 5). Molecules M1 and M2 were accompanied by about a dozen of compounds with similar 1 H NMR properties (Figure 4), several of those showing tert-butyl groups as well, which frequently were connected to aromatic rings. The 1 H, 13 C HMBC NMR spectrum of L05_1000m presented remarkably strong cross peaks derived from ketones which very probably are marine derived natural products (Figures 4, S15). Tert-butylphenols readily produce phenoxyl radicals and are used for stabilizing plastics; they are feedstock for the industrial synthesis of more complex antioxidants (Gerritsen et al., 1998;Remberger et al., 2004;Dhawan and Cholli, 2006;Yoon et al., 2006;Tollefsen et al., 2008, Zhang et al., 2013Papa et al., 2014, Zhao et al., 2020. Additionally, they are found in inflows and effluents of municipal sewage treatment plants (Remberger et al., 2004). They also contribute to the endocrine disrupting potential in natural waters (Tollefsen et al., 2008), and are subject to bioaccumulation in fish (Papa et al., 2014) and other organisms (Gerritsen et al., 1998). The uptake, transformation and elimination within a food web depend on individual characteristics of organisms. In addition, 2,4-di-tertbutylphenol is a remarkably widespread metabolite of plants and bacteria also in the ocean (Chawawisit et al., 2015;Padmavathi et al., 2015), overall being produced by almost 200 species (Zhao et al., 2020). It shows a distinct potent broadband toxicity against almost all testing organisms, including the producing species itself (Zhao et al., 2020), and inhibits microbial quorum sensing in vitro (Chawawisit et al., 2015;Padmavathi et al., 2015). Photochemical a-cleavage of di-tbutylketone in the presence of phenols leads to ortho, para-tbutylphenols, and t-butyl phenyl ethers (Galindo et al., 1998). Tert-butylbenzene was a key oil contamination, which could be degraded by bacteria in deep water (Hazen et al., 2010;Lu et al., 2012). Bergauer et al. (2018) indicated that carboxylic acids could be one of essential substrates for the microbial community in the bathypelagic, and aromatic compounds were more likely to become utilized by heterotrophic prokaryotic community. The observed depth distribution of compounds M1, M2 and associated molecules (Figures 2, 5, S2, S3, S17) was remarkable because relative proportions were maximal at meso-and bathypelagic depths and significantly exceeded near surface concentrations; specifically, M1 and M2 were virtually absent in the chlorophyll fluorescence maximum of both sampling stations L05_80 m and L11_90 m. Such specific low molecular weight phenols and carboxyl aromatics have not been previously reported in deep marine waters (1000 m), suggesting possible intermediate products in microbial food webs. Further studies are essential to demonstrate whether the abundance of M1 and M2 at the meso-and bathypelagic depths is related to the mesoscale eddies. 1 H, 1 H COSY NMR spectra of marine SPE-DOM isolated from waters with high biological activity are very rich in detail and depict atomic environments of hundreds of organic molecules (Hertkorn et al., 2013). Common COSY cross peaks in marine DOM originate from relative abundant molecules, and direct comparison of those information-rich NMR spectra allows assessment of relative similarities and deviations of putative biological production in different ocean waters. It was therefore instructive to compare 1 H, 1 H COSY NMR spectra of SCS SPE-DOM L05_1000m with that of another previously recorded marine system, Atlantic Ocean fluorescence maximum SPE-DOM obtained at 48m depth (Figures 6 and Figure S18). 1 H, 1 H TOCSY NMR spectra with short mixing time (e.g. 15 ms) reveal short range J-coupled spin systems ( 2/ 3 J HH ) alike 1 H, 1 H COSY NMR spectra, but with substantially higher sensitivity (Hertkorn, 2014) making these an obvious choice because only limited quantities of SPE-DOM L05_1000m were available. (Bruch, 1996;Keeler, 2011). Gonsior et al. (2014) observed very similar 1 H NMR spectra of Atlantic and Pacific surface ocean SPE-DOM, indicative of common features governing the formation and degradation of highly processed biogeochemical DOM in these open ocean waters. In contrast, comparison of SPE-DOM isolated from Atlantic Ocean fluorescence maximum and South China Sea L05_1000m waters revealed presence of fundamentally different but rather abundant organic molecules, in particular referring to aliphatic OCH and CCH units (Figures 6, S18). While both samples showed a higher relative concordance of their aromatic C ar H units, distinction of cross peaks was visible as well in the major sections with electron withdrawing substitution, polycyclic aromatic hydrocarbons and six-membered N-heterocycles (d H > 7.3 ppm), in the sections of C ar H with neutral (d H~7 -7.3 ppm) and electron donating substitution, as well as olefins and five membered heterocycles (d H < 7 ppm) ( Figures 6A, C). The two major molecules M1 (2,4-di-tert-butylphenol) and M2 (3,3'-(3tert-butylbenzene-1,5-diyl)dipropanoic acid) showed a few, not too influential strong cross peaks for C ar H units while tert-butyl groups did not show appreciable TOCSY cross peaks from long range couplings ( Figure S15B). The main distinction of both samples resided in the oxygenated aliphatic groups (OCH and OCCH units) and in the functionalization of aliphatic carboxylic acids; the latter were much more abundant and chemically diverse in SCS L05_1000m than in Atlantic Ocean FMAX SPE-DOM. Carbohydrate-type OCH cross peaks (-OCH-OCH units, and alkylated carbohydrates: -OCH-CH units) were much more abundant in Atlantic Ocean FMAX SPE-DOM, suggesting earlier stages of organic matter production than observed in SCS L05_1000m SPE-DOM. Here, tert-butyl derivatives of carboxylic acids were likely contributors to some but not to the majority of these TOCSY cross peaks (HOOC-CaH-CbH-Z units; with Z including remote tert-butyl groups C 4 H 9 ).
Taking the bulk parameter DOC of SCS samples into consideration, the total concentration of these specific compounds M1 and M2 as computed from 1 H NMR integrals (ratio of compound 1 H NMR integral to total 1 H NMR integral of SPE-DOM multiplied by DOC concentration) was maximal (~6 mmol L -1 ) in L05_1000 m, compared with the other depths. Such high abundance of putatively bioactive compounds at meso-and bathypelagic depths in SCS waters is truly remarkable. Since M1 and M2 were virtually absent in the chlorophyll fluorescence maximum of both the SCS water SPE-DOM (L05_80 m and L11_90 m), and the Atlantic Ocean water SPE-DOM (at 48 m), the enrichment of M1 and M2 in deep waters revealed the prominent effect of eddy on the production and depth distribution of these marine metabolites. Such production of such "unusual" molecules at depth has scarcely been reported, in part because only NMR spectroscopy provides straightforward access to these "unknown unknowns". Most importantly, tert-butylphenols are endocrine disruptive compounds, with distinct potent broadband toxicity. These molecules may accumulate in organisms along the food chain and have long-term effects to marine system, not only in surface but even in deep ocean waters (> 1000 m). Mesoscale marine chemical hotspots might induce chemotaxis of susceptible organisms (Raina et al., 2022) not only contributing to further processing of marine DOM but also influencing the productivity of marine ecosystems. In addition, physicochemical mechanisms of exporting unique molecules from and to meso-and bathypelagic depths are complex and ill constrained. Downwelling induced by ACE could penetrate down to 1000 m depth and residuals of unusual (bio)molecules could benefit from the "dilution effect" (Arrieta et al., 2015;Seidel et al., 2022). Zhang et al. (2014) observed deep-sea sediment transport processes driven by mesoscale eddies in the northern South China Sea at 2100 m water depth, close to our study zone, and deep-reaching eddies were shown to exert significant influence on the settling velocity of marine snow and other particulates (Purkiani et al., 2020). Better understanding of TEP (transparent exopolymer particles) -mediated biogeochemical and ecological processes in deep oceans (Nagata et al., 2021) is needed. An enhanced understanding of lateral advection (from shelves, continents, or islands) and the dissolution from sinking particulate organic matter are other conceivable contributors of relevance (Seidel et al., 2022).

CONCLUSION
Oceanic mesoscale eddies with horizontal scales of 50-300 km modulate the ocean biogeochemistry and are very influential in FIGURE 6 | Homonuclear two dimensional 1 H, 1 H NMR spectra and 1 H NMR projection NMR spectra above of SCS L05_1000 m SPE-DOM (A, B, TOCSY NMR with 15 ms mixing time) and of open Atlantic Ocean fluorescence maximum zone SPE-DOM at 48 m depth (C, D, COSY NMR spectrum) (Hertkorn et al., 2013) demonstrate considerable biomolecular diversity of these two marine DOM. NMR spectra indicate an earlier, carbohydrate-rich stage of organic matter processing following primary production in open Atlantic Ocean surface water, compared with that of South China Sea deep water, which is further oxidized showing a high proportion of carboxyl-rich alicyclic molecules CRAM (Hertkorn et al., 2006) (cf. text; for assignment of key substructures see the supporting online material). the transport of heat, nutrients, DOC and other organic and inorganic compounds, and they play a decisive role in the confined South China Sea. While cyclonic eddies often inject nutrients from deep waters and thereby increase primary production in near surface waters, anticyclonic eddies (ACE) are characterized by downwelling at their core regime, with supposedly minor associated biogeochemical effects. While ACE-induced marine currents were not expected to fundamentally change the nature of its marine DOM as an extremely complex mixture of organic molecules, with all associated challenges of its molecular and structural analysis, a certain confinement of water masses on week-long time and kilometer size scales apparently decreased the overall molecular diversity of DOM of a depth profile from surface to bottom waters within an ACE, compared with the edge of an ACE. While small molecules of putative biological origin, such as 2,4-di-tertbutylphenol (C 14 H 22 O) and 3,3'-(5-(tert-butyl)-1,3-phenylene) dipropionic acid (C 16 H 22 O 4 ), were present in all waters, NMR spectroscopy enabled identification of specific tert-butyl derivatives at meso-and bathypelagic depths in single digit micromolar concentrations, and offered lead structures for other, yet unidentified olefinic and ketone-derived marine natural products. 2,4-di-tert-butylphenol could be a widespread product of marine algae and bacteria with potential for bioaccumulation within marine food webs, and associated toxicity for plentiful organisms. 3,3'-(3-tert-butylbenzene-1,5diyl) dipropanoic acid derived radicals could be the degraded by heterotrophic prokaryotic communities in the bathypelagic ocean, implying a potential contribution of microbial carbon pump in the marine system. Deviations from smooth mass peak distributions in FTICR mass spectra provided rather indirect evidence of biological processes associated with ACE, and mathematical analysis revealed alterations of biogeochemical processing at different positions and water depths. Active microbial metabolism across the entire water depth column in regions of up-and downwelling might create more abundant organic molecules than currently anticipated. Our study highlights the necessity to deciphering the structure and reactivity of marine organic molecules in a kinetic context, including the microbial and physicochemical constraints especially in the deep ocean (Gonsior et al., 2022).
These findings again emphasize the need for in-depth NMR spectral analyses for identification of "unknown unknowns" in marine waters, because FTICR mass spectrometry alone could not have produced significant lead information. 1 H detected NMR spectroscopy could identify the structures of key molecules including quantification. Less superimposed 2D NMR spectra provided comprehensive information of key atomic chemical environments, demonstrating an active oceanic biogeochemistry of these unusual molecules. Future studies which combine investigation of spatially and temporally resolved microbial diversity, external (bulk) conditions and molecular features of DOM might select oceanic eddies as an attractive example for investigating the effect of intermediate size gradients on element, nutrient and molecule cycles. This will better constrain the relative proportions of biotic and abiotic interactions governing the global carbon and other element cycles. These studies will propose relationships between groups of organisms and arrays of molecules, which allow testable predictions about biogeochemical processing in ocean waters.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
MZ and YW contributed to the research design and sample collection. MZ, SL, and NH performed NMR acquisition, data analyses and interpretation. MZ, SL, and MH performed FTICR/MS acquisition and data analyses. PS-K provided support for the FTICR/MS analyses. DX provided the physical parameters. MZ, YW, SL, NH, and MH participated paper writing. All authors provided significant input on the final manuscript.