Variation of Isoprenoid GDGTs in the Stratified Marine Water Column: Implications for GDGT-Based TEX86 Paleothermometry

Isoprenoid glycerol dialkyl glycerol tetraethers (isoGDGTs) derived from archaea are lipid biomarkers that exhibit high sensitivity to changes in water temperature, leading to the widespread application of the isoGDGT-based tetraether index of 86 carbon atoms (TEX86) in surface seawater temperature (SST) reconstruction. However, there remain some uncertainties regarding the robustness of TEX86 under changing water conditions (e.g., variations in water depth, oxygen and pH). Here, we analyzed isoGDGTs in suspended particles at different depths of the East China Sea (ECS) during summer 2020, aiming to constrain the applicability of the TEX86 proxy in coastal waters. Our data showed that the isoGDGTs were mainly derived from planktonic Thaumarchaeota, as revealed by the low ratio of GDGT-0/crenarchaeol (<0.5). The vertical distribution of isoGDGT concentration depicted a downward increase from the surface to the bottom. This observation was likely shaped by Thaumarchaeota, which regulate the extent of ammonia oxidation based on the availability of ammonium. The occurrence of maximal isoGDGT concentrations in the bottom layer suggests that the isoGDGTs in sediments are mainly controlled by bottom archaeal production rather than surface archaeal production. By reanalyzing the published isoGDGT data of surface sediments in the ECS inner shelf, we found that the sedimentary TEX86 relates much better to the annual mean bottom seawater temperature (BST) than to the annual mean SST, indicating that sedimentary TEX86 is more inclined to be a proxy for the BST in the shallow ECS. In addition, the positive bias of TEX86 driven by low dissolved oxygen and low pH was observed under the pycnocline, indicating that the application of TEX86 to reconstruct seawater temperature should be carefully appraised in coastal environments with strong water column stratification.

Isoprenoid glycerol dialkyl glycerol tetraethers (isoGDGTs) derived from archaea are lipid biomarkers that exhibit high sensitivity to changes in water temperature, leading to the widespread application of the isoGDGT-based tetraether index of 86 carbon atoms (TEX 86 ) in surface seawater temperature (SST) reconstruction. However, there remain some uncertainties regarding the robustness of TEX 86 under changing water conditions (e.g., variations in water depth, oxygen and pH). Here, we analyzed isoGDGTs in suspended particles at different depths of the East China Sea (ECS) during summer 2020, aiming to constrain the applicability of the TEX 86 proxy in coastal waters. Our data showed that the isoGDGTs were mainly derived from planktonic Thaumarchaeota, as revealed by the low ratio of GDGT-0/crenarchaeol (<0.5). The vertical distribution of isoGDGT concentration depicted a downward increase from the surface to the bottom. This observation was likely shaped by Thaumarchaeota, which regulate the extent of ammonia oxidation based on the availability of ammonium. The occurrence of maximal isoGDGT concentrations in the bottom layer suggests that the isoGDGTs in sediments are mainly controlled by bottom archaeal production rather than surface archaeal production. By reanalyzing the published isoGDGT data of surface sediments in the ECS inner shelf, we found that the sedimentary TEX 86 relates much better to the annual mean bottom seawater temperature (BST) than to the annual mean SST, indicating that sedimentary TEX 86 is more inclined to be a proxy for the BST in the shallow ECS. In addition, the positive bias of TEX 86 driven by low dissolved oxygen and low pH was observed under the pycnocline, indicating that the application of TEX 86 to reconstruct seawater temperature should be carefully appraised in coastal environments with strong water column stratification.

INTRODUCTION
As an important life form in the ocean, microorganisms play a significant role in the ocean matter cycle (Copley, 2002). Archaea are ubiquitous in the marine environment (Lipp et al., 2008;Hoshino and Inagaki, 2019). In view of their unique living habits and diverse metabolic types, archaea regulate the carbon and nitrogen cycles and energy exchange in the ocean (Copley, 2002;del Giorgio and Duarte, 2002;Francis et al., 2005;Martens-Habbena et al., 2009;Jiao et al., 2010;Zhang et al., 2020). Especially under global change, the role of archaea in biogeochemical processes has become increasingly important (Zhang et al., 2018).
The evolution of environmental conditions can be recorded by archaea through the production of specific organic molecules (Bianchi and Canuel, 2011), including glycerol dialkyl glycerol tetraethers (GDGTs). GDGTs are an important component of cell membrane lipids that can be preserved in the environment for a long time, thereby becoming an excellent paleoarchive Hopmans et al., 2004;Pearson and Ingalls, 2013). GDGTs derived from archaea are commonly referred to as isoprenoid GDGTs (isoGDGTs). The isoGDGT carbon chain skeleton has a unique isoprenoid structure and is connected to the glycerol molecule by ether bonds. GDGTs are mainly composed of compounds containing 0-4 cyclopentane rings (GDGT-0 to GDGT-4). However, crenarchaeol has 4 cyclopentane rings and 1 cyclohexane ring and is generally considered to be a unique biomarker for Thaumarchaeota (Figure 1; Sinninghe Weijers et al., 2006;Pester et al., 2011;Schouten et al., 2013). The study of isoGDGTs has spawned a series of environmental proxies. One of these proxies is the tetraether index of 86 carbon atoms (TEX 86 ) established based on the relative abundance of the isoGDGT component, which has been widely used in the reconstruction of seawater surface temperature, and the time scale spans more than 100 million years of geological history (Schouten et al., , 2013Kim et al., 2010;Jenkyns et al., 2012;Lü et al., 2019b). However, the use of TEX 86 in sedimentary cores to retrieve surface temperature information is based on the assumption that the isoGDGTs are mainly derived from surface archaeal production, but this assumption still prompts controversy. Previous research has found that Thaumarchaeota prevail in the deeper part of the Mediterranean Sea water column, and the temperature retrieved by the sedimentary TEX 86 cannot reflect the annual mean surface temperature (Besseling et al., 2019). There remain several uncertainties regarding the utility and reliability of TEX 86 under changing water conditions (e.g., changes in water depth, oxygen and pH) (Huguet et al., 2009;Blaga et al., 2011;Boyd et al., 2011;Zhu et al., 2011;Jia et al., 2012;Lü et al., 2014;Qin et al., 2015;Ho and Laepple, 2016;Wang et al., 2019). Especially in summer, the stratification of marginal seas under the influence of large rivers develops extensively, and the bottom layer in the water column usually has a lower dissolved oxygen (DO) content and pH. However, the effect of seawater stratification on the efficacy of TEX 86 is still unclear. Recent studies have shown that in stratified lakes, TEX 86 is primarily controlled by the DO content, not the temperature (Zhang et al., 2016). In addition, the study of isoGDGTs in the water column of the marginal sea did not find correlation between bottom TEX 86 values and the in situ temperature during summer (Wang et al., 2019).
The vertical variation of isoGDGT producers in the water column directly affects the response of temperature signals incorporated in the sedimentary TEX 86 to water depth. Thaumarchaeota, a phylum in the ammonia-oxidizing archaea (AOA) class, are the major source of isoGDGTs (Brochier-Armanet et al., 2008;Pitcher et al., 2011;Pearson and Ingalls, 2013). As one of the predominant ammonia oxidizers, Thaumarchaeota play an important role in marine nitrification (Könneke et al., 2005;Pitcher et al., 2011;Schouten et al., 2013). The prevalence of autotrophic nitrifying archaea in the water column is controlled by the availability of ammonium, which further drives the change in isoGDGT abundance (Pitcher et al., 2011). The coupling of AOA and crenarchaeol has been confirmed in many experiments, especially the excellent correspondence between the gene copy number of amoA (the functional gene for ammonia oxidation) and crenarchaeol abundance, indicating that crenarchaeol could be used as a tracer of ammonia oxidation (Leininger et al., 2006;Schouten et al., 2008;Pitcher et al., 2009Pitcher et al., , 2011.
Clarifying the vertical distribution pattern of isoGDGTs and their driving factors can better constrain the application of TEX 86 . However, in marginal seas, the variability of isoGDGT concentration is still unclear due to the dynamic interaction of physical and biochemical processes, as well as the source and transformation of nitrogen are extremely complicated.
Here, we collected samples of suspended particulate matter from the East China Sea (ECS) and measured a series of biochemical parameters to advance the understanding of the vertical distribution pattern of isoGDGTs in the stratified water column and reassess the records and robustness of TEX 86 .

Study Area and Sampling
The ECS is a typical marginal sea influenced by large rivers and characterized by complex hydrological conditions. The world's third-largest river, the Changjiang River, carries a large amount of fresh water and nutrients into the ECS annually. Driven by the East Asian monsoon, multiple currents converge in the ECS and interact in a dynamic manner (Guo et al., 2018). Affected by the plume of the Changjiang River, the water column of the ECS becomes stratified during summer, and the surface layer possesses high primary production due to favorable light and temperature and sufficient nutrients (Gong et al., 2003).
Sampling was conducted aboard the R/V Kexue 3 from August 8 to August 20, 2020. Water samples were collected from discrete depths (2, 10, 20, 30, 50 m and bottom layer) at five stations (Figure 2). A Seabird conductivity-temperaturedepth (CTD) system incorporated with 12-L Niskin bottles was used for seawater collection and to obtain the in situ salinity and temperature. Water samples for particulate isoGDGTs and particulate organic nitrogen (PON) were filtered through Whatman GF/F membranes (0.7 µm pore size; precombusted at 450 • C for 5 h) immediately after collection. The filtrate was collected in duplicate in acid-cleaned, 60-mL high-density polyethylene (HDPE) bottles and stored frozen (-20 • C) for further analysis of nutrients. Samples for chlorophyll-a (Chla) analysis were collected using cellulose acetate fiber filters (0.45 µm pore size; Whatman) and stored in the dark at -20 • C. Water samples for heterotrophic bacterial abundance (HBA) analysis were collected in cryopreservation tubes and immediately fixed with 1 mL of 1% paraformaldehyde. Then, the samples were stored in liquid nitrogen until analysis in the laboratory. The surface sediments (top 0-2 cm) were obtained using a box sampler, and the samples were stored at -20 • C.

DO, pH, and Chl-a
Dissolved oxygen and pH were determined by electrochemical probes immediately following sample collection (SevenExcellence, Mettler Toledo) (Guo et al., 2020). The relative accuracy of the pH measurement was ±0.002, and the pH was corrected according to the in situ temperature. Traditional Winkler titration (Bryan et al., 1976) was performed to calibrate the DO values obtained by the probe, and the two DO results agreed well (R 2 = 0.98, p < 0.01).
The samples (cellulose acetate fiber filters) for the determination of Chl-a were extracted in 10 mL of N,N-dimethyl formamide for 12 h and then centrifuged. The supernatant was subjected to fluorescence measurement with a fluorescence spectrophotometer (F-4500, Hitachi Co, Japan). The detection limit was 1.5 µg L −1 , and the standard deviation (SD) was <1% (n = 11; Guo et al., 2020).

HBA Analysis
The HBA was determined by fluorescence staining (Zhao et al., 2011). Briefly, the water samples were diluted six times with TE buffer (Tris-EDTA, 100 mmol L −1 Tris-Cl, 10 mmol L −1 EDTA, pH = 8.0, Sigma) and stained with the nucleic acid dye SYBR Green I (Molecular Probes). The reaction continued for 20 min at room temperature in the dark and was then subjected to flow cytometry (BD FACSVantage SE, Becton, Dickinson, United States).

Inorganic and Organic Nitrogen
Nitrite (NO 2 -N) and ammonium (NH 4 -N) were determined spectrophotometrically on the ship immediately after pretreatment, whereas nitrate (NO 3 -N) was analyzed using a Bran-Lubbe Quaatro-SFA autoanalyzer in the laboratory (Duan et al., 2016). Total dissolved nitrogen (TDN) was oxidized by potassium sulfate to NO 3 -N and then measured with an autoanalyzer (Yuan et al., 2018). The SD of three injections was typically <3%. The dissolved organic nitrogen (DON) was calculated as the difference between TDN and the sum of NO 3 -N, NO 2 -N, and NH 4 -N. Concentrations of PON were determined using an autoanalyzer after alkaline potassium persulfate oxidation and with an SD less than 5% for triplicate measurements (Yuan et al., 2018).

Extraction and Analysis of isoGDGTs
The methods of lipid extraction were described by Duan et al. (2019). Briefly, freeze-dried GF/F filters or powdered sediments with an added internal standard (i.e., C 46 glycerol trialkyl glycerol tetraether, C 46 GTGT) (Huguet et al., 2006) were ultrasonically extracted four times using a mixture of dichloromethane (DCM) and methanol (MeOH) (v/v = 3:1). The extracted samples were dried with nitrogen (N 2 ) and then saponified with 6% KOH/MeOH (w/w) overnight. Then, the polar and nonpolar fractions were separated in a silica gel column by elution with n-hexane and DCM/MeOH (95:5, v/v), respectively. The polar eluate was concentrated with N 2 and redissolved in n-hexane/isopropanol (v/v = 99:1) for instrumental analysis.
The TEX 86 index was calculated using the following equations proposed by Schouten et al. (2002):
The temperature and salinity values in the upper 20 m of the nearshore stations (DH4-0 and DH5-1a) changed dramatically, whereas the salinity was vertically uniform at the offshore stations (DH5-2, DH5-2a, and DH5-3). The pH showed a narrow range (7.98-8.24) and was distributed relatively uniformly throughout the water column ( Figure 3C). Concentrations of DO ranged from 125.6 to 238.2 µmol L −1 and shared a similar distribution pattern with the temperature (Figure 3D). Compared with those of offshore stations (DH5-2, DH5-2a, and DH5-3), the values of temperature, salinity, pH and DO at nearshore stations (DH4-0 and DH5-1a) were relatively low. In contrast, the Chl-a concentrations (0.08-3.63 µg L −1 ) were higher at nearshore stations (DH4-0 and DH5-1a; mean ± SD: 0.98 ± 1.09) than at offshore stations (DH5-2, DH5-2a, and DH5-3; mean ± SD: 0.37 ± 0.32) (Figure 3E). A remarkable vertical gradient of Chl-a concentrations was observed at nearshore stations DH4-0 and DH5-1a. The depth profile showed that the HBA values were elevated near the surface (upper 20 m) and decreased with increasing depth (Figure 4F). The concentrations of NH 4 -N, NO 2 -N, and NO 3 -N ranged from 0.47 to 2.92, 0.01-0.91, and 0.01-12.87 µmol L −1 , respectively (Figures 3F-H). The vertical distributions of NH 4 -N and NO 2 -N were highly variable. In comparison, the NO 3 -N concentration generally exhibited an increasing trend from the surface to the bottom except for the obviously high value at the surface of station DH4-0. Concentrations of DON ranged from 23.69 to 51.82 µmol L −1 (average ± SD: 33.91 ± 6.49 µmol L −1 ) but displayed no apparent pattern with increasing depth (Figure 5A). The vertical distribution of PON was similar to that of Chl-a and varied between 0.39 and 2.51 µmol L −1 (Figure 5B).

Concentrations and Distributions of isoGDGTs in the Water Column
The concentrations of isoGDGTs were between 3.66 and 52.06 ng L −1 (average ± SD: 19.63 ± 14.75 ng L −1 ) and gradually increased with increasing depth (Figure 4A). The isoGDGT concentrations in this study were comparable to those (8.30-237.10 ng L −1 ) reported by Wang et al. (2019) in the same season but lower than those (27.60-6,230 ng L −1 ) reported by Lü et al. (2019a). Similar to NO 3 -N, relatively high concentrations of isoGDGTs were observed at nearshore stations DH4-0 and DH5-1a. Cren and GDGT-0 accounted for 66.2 ± 6.1 and 15.5 ± 5.2% of the total isoGDGT content, respectively (Figures 4B,C). Cren and GDGT-0 showed similar distribution patterns, and the ratio of GDGT-0/Cren ranged from 0.08 to 0.43 (mean ± SD: 0.24 ± 0.10) among all samples ( Figure 4D). The TEX 86 ranged from 0.50 to 0.86 and displayed considerable variability in the water column ( Figure 4E).

Source and Variation of isoGDGTs in the Stratified Water Column
The origins of isoGDGTs can be evaluated by the ratio of GDGT-0/Cren. A low ratio value (<2.0) of GDGT-0/Cren indicates a main source from Thaumarchaeota, while a high ratio value (>2.0) suggests the existence of methanogenic archaea (Blaga et al., 2009). In the present study, GDGT-0/Cren values between 0.08 and 0.43 suggest the dominance of Thaumarchaeota-derived isoGDGTs in the ECS. This is also supported by the archaeal 16S rRNA gene data (Zeng et al., 2007).
The vertical distribution of isoGDGT concentration depicted an increasing trend from the surface to the bottom ( Figure 4C). This vertical trend was coupled with the activity of Thaumarchaeota in the water column (Pitcher et al., 2011). Previous studies on the ECS ammonia monooxygenase gene showed that the amoA in bottom water was significantly higher than that in the surface layer . This result is consistent with our findings, indicating that Thaumarchaeota prefer to thrive in the deeper layers of the ECS. To further analyze the influence of various physical, chemical and biological parameters on the distribution of isoGDGTs, redundancy analysis (RDA) was performed in this study. The RDA results showed that the first principal component explained 99.09% of the variation in the isoGDGT data (Figure 6). Temperature, DO and pH, as well as nitrogen, presented a strong influence on the distribution of isoGDGTs. The relatively low oxygen and low pH microenvironment of the bottom particles are conducive to the ammonia oxidation process of Thaumarchaeota, which has been confirmed by previous studies (Hanaki et al., 1990;Hsiao et al., 2014). As a reactant in ammonia oxidation, the availability of NH 4 -N plays a leading role in the prevalence of Thaumarchaeota (Könneke et al., 2005;Yan et al., 2012;Evans et al., 2018). In summer, stratification is widely developed in the ECS, especially at the nearshore DH4-0 and DH5-1a stations. The pycnocline separates the vertical exchange of nutrients, resulting in two different patterns of nitrogen cycling in the water column. Above the thermocline, the high phytoplankton production, as revealed by the high Chl-a content (Figure 3E), consumed a large amount of NH 4 -N. At the same time, the fresh organic matter produced was quickly decomposed by heterotrophic bacteria, which also stimulated heterotrophic bacterial blooms at the surface ( Figure 4F). The NH 4 -N released by mineralization was utilized by phytoplankton. This rapid NH 4 -N cycling made it difficult for nitrogen to accumulate and nitrify at the surface, showing some characteristics oligonitrates (excluding the DH4-0 surface NO 3 -N concentration, which was severely affected by terrestrial input). In this pattern, Thaumarchaeota can hardly compete with phytoplankton for NH 4 -N (Murray et al., 1999;Galand et al., 2010). Furthermore, the intense radiation in the surface water in summer also exerted a negative effect on the growth of Thaumarchaeota (Merbt et al., 2012). The NO 3 -N concentration at the surface of station DH4-0 was remarkably high (Figure 3H), which was attributed to terrestrial input, as revealed by the low salinity ( Figure 3B). However, below the pycnocline, the stress from phytoplankton competition gradually decreased due to light limitation. The NH 4 -N released by the mineralization of sedimentation particles can be used by Thaumarchaeota to mediate the nitrification process and promote the accumulation of NO 3 -N . Even when the concentrations of NH 4 -N in the ambient environment are low, archaea can still maintain thriving growth conditions (Martens-Habbena et al., 2009). Although vertical variations in isoGDGTs were observed, the expected increase in NO 2 -N with depth was not found, which may be due to its active chemical properties (such as rapid oxidation). However, by excluding the DH4-0 surface data, the strong correlations between isoGDGTs and NO 3 -N (R 2 = 0.75, p < 0.01) and between Cren and NO 3 -N (R 2 = 0.74, p < 0.01) indicate that archaeal ammonia oxidation plays an important role in the ECS nitrogen transformation. Hsiao et al. (2014) found that the nitrification rate was significantly correlated with PON but not with DON in the ECS, implying that PON had a greater potential to provide NH 4 -N than DON. The coupling of relatively high PON and isoGDGTs at DH4-0 and DH5-1a (Figures 4A, 5B) and a significant correlation between DO and PON (R 2 = 0.52, p < 0.01) further indicated that PON as  a remineralization substrate was more effective in providing NH 4 -N. However, for DON, its vertical distributions were highly variable (Figure 5A), and no significant relationships between DON and DO and isoGDGTs were found (p > 0.05), which may be related to its low bioavailability and diverse sources Carlson and Hansell, 2015).

Re-evaluating the Records and Robustness of TEX 86
As a bioarchive of seawater temperature, TEX 86 has been widely used and reported (Pearson and Ingalls, 2013;Schouten et al., 2013;Lü et al., 2015Lü et al., , 2019a. Since TEX 86 was proposed, its statistical relationship with surface seawater temperature (SST) has undergone many modifications, but it is not feasible to obtain accurate results in different regions with one formula (Kim et al., 2010;Tierney and Tingley, 2014). Especially in the marginal seas, due to the complex hydrodynamic process and the disturbance of human activities, the records and robustness of TEX 86 still need to be further evaluated.
The reconstruction of paleo-SST by TEX 86 is based on the assumption that sedimentary isoGDGTs originated from surface archaea. However, our study found that the isoGDGT concentration increased with increasing depth in the water column. This indicated that, at least in summer, the isoGDGTs in ECS surface sediments are not derived from the surface and are mainly controlled by bottom archaeal production. The depth profile including the surface sediments showed that the surface sediments and the bottom layer exhibit similar isoGDGT concentrations (Figure 7), further supporting this hypothesis. From this perspective, the sedimentary TEX 86 cannot accurately retrieve the surface water temperature, and it is more likely to record the bottom temperature signals. This also explains why TEX 86 in surface sediments usually does not accurately represent surface temperature in summer (Duan et al., 2019). To further verify the suitability of sedimentary TEX 86 as a bottom temperature proxy, we collected and reanalyzed the previously published isoGDGT data of surface sediments in the ECS (Lü et al., 2014). We found that the sedimentary TEX 86 had a stronger relationship with the annual mean bottom seawater temperature (BST) (R 2 = 0.82, p < 0.01; Figure 8B) than with the annual mean SST (R 2 = 0.22, p < 0.01; Figure 8A), indicating that the sedimentary TEX 86 in the ECS is more suitable as a proxy for BST. The optimal relationship between sedimentary TEX 86 and annual average BST was also found by Wang et al. (2019) in the eastern China marginal sea, further confirming that sedimentary TEX 86 is a proxy for ECS bottom temperature. In this scenario, in coastal settings, sedimentary isoGDGTs could be mainly derived from bottom archaeal production, and sedimentary TEX 86 best reflects BST. However, in the open sea, the sedimentary  Lü et al. (2014). Annual mean seawater temperature data were obtained from NOAA website (WOA13; https://www.nodc.noaa.gov/). isoGDGTs mainly originate from the subsurface water, resulting in sedimentary TEX 86 reflecting the subsurface temperature (Huguet et al., 2007;Hernández-Sánchez et al., 2014).
In addition to the temperature records by TEX 86 , whether TEX 86 only depends on temperature is also debatable. In this study, the relationship between TEX 86 and temperature presented two scenarios (Figure 9). When the temperature was greater than 25.0 • C, TEX 86 was significantly related to temperature (R 2 = 0.58, p < 0.01), while this relationship disappeared when the temperature was below 25.0 • C. Previous studies have found that DO and pH could also affect TEX 86 , which was confirmed in pure archaeal culture experiments (Oger and Cario, 2013;Elling et al., 2015;Qin et al., 2015). Significant correlations were found between TEX 86 and DO (R 2 = 0.14, p < 0.01) and pH (R 2 = 0.36, p < 0.01), indicating that TEX 86 does not only serve temperature in the ECS during summer. Samples with temperatures below 25.0 • C have relatively low DO content and pH values driven by stratification, and TEX 86 increased compared to the relational equation for samples with a temperature greater than 25.0 • C. This phenomenon was consistent with previous findings, indicating that a decrease in DO or pH could lead to an increase in TEX 86 Zhang et al., 2016;Lü et al., 2019a;Cao et al., 2020). For samples with a temperature below 25.0 • C, a significant correlation was found between pH and TEX 86 (R 2 = 0.33, p < 0.05), but no significant correlation was found between DO and TEX 86 , indicating that the effect of pH on TEX 86 was greater than that of DO. However, more data are needed to confirm this hypothesis. Another argument for TEX 86 bias is the supply of NH 4 -N. Recent studies have shown that the low ammonia oxidation rate and energy limitation caused by a low supply of NH 4 -N resulted in a positive deviation of TEX 86 (Hurley et al., 2016;Evans et al., 2018). Lü et al. (2019a) investigated isoGDGT concentrations FIGURE 9 | Relationship between TEX 86 and temperature in the East China Sea (ECS) during summer. Inside the red rectangles are samples characterized by the relatively low dissolved oxygen (DO) and pH. The temperature is derived from in situ CTD measurement.
in the water column of the ECS and found that the deviation between the reconstructed temperature based on isoGDGTs and the in situ temperature gradually increased with increasing water depth in summer. They attributed this phenomenon to the reduction in the ammonia oxidation rate at the bottom layer. However, in this study, no significant relationship was found between NH 4 -N and TEX 86 , indicating that NH 4 -N had a limited effect on TEX 86 . Alternatively, the rapid turnover of NH 4 -N in the water column could have caused the statistical relationship between NH 4 -N and TEX 86 to be lost. Furthermore, the growth phase of archaea also has a profound effect on the composition of membrane lipids, and the bias of TEX 86 -derived temperature between the growth and stationary phases can be up to 9 • C (Elling et al., 2014). However, it is difficult to distinguish the contribution of various influencing factors in this study. Overall, the coupling of multiple factors leads to the divergence of TEX 86 to the summer temperature in the ECS. More seasonal data are needed to further understand this issue.

CONCLUSION
Our data presented in this study reveal that Thaumarchaeota were the predominant producers of isoGDGTs in the water column of the ECS. The isoGDGT concentration increased with increasing depth, which may be related to Thaumarchaeotal ammonia oxidation driven by the availability of NH 4 -N. Archaea prefer to produce in relatively deep water, indicating that sedimentary TEX 86 could not record summer surface temperatures in the ECS. By compiling previously published isoGDGT data for the ECS surface sediments, we found that sedimentary TEX 86 exhibited a stronger relationship with the annual mean BST (R 2 = 0.82, p < 0.01) than with the annual mean SST (R 2 = 0.22, p < 0.01). One major implication is that sedimentary TEX 86 is more reliable in reflecting BST than reflecting SST in the shallow ECS (<100 m). In addition, in environments where the DO content and pH were low, the relationship between TEX 86 and temperature varied, resulting in a bias toward warmer temperatures in the TEX 86 -derived temperature, indicating that the robustness of TEX 86 could be weakened in coastal settings when the water column is stratified during summer. Further work in combination with archaeal community and intact polar ether lipids analysis are needed to better decipher the production and variation mechanisms of isoGDGTs.

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
JG: conceptualization, formal analysis, investigation, and writing original draft. HY: methodology, resources, review and editing, and Supervision. JS: project administration, supervision, and funding acquisition. XL: data curation. LD: validation. BQ, JX, QW, YW, and NL: formal analysis, methodology, and investigation. All authors contributed to the manuscript revision, read, and approved the submitted version.