The inﬂuence of near-surface sediment hydrothermalism on the TEX 86 tetraether-lipid-based proxy and a new correction for ocean bottom lipid overprinting

. The diversity and relative abundances of tetraether lipids produced by archaea and bacteria in soils and sediments are increasingly used to assess environmental change. For instance, the TetraEther indeX of 86 carbon atoms (TEX 86 ), based on archaeal isoprenoidal glycerol dialkyl glycerol tetraether (iGDGT) lipids, is frequently applied to reconstruct past sea-surface temperatures (SSTs). Yet, it is unknown how the ratio fully responds to environmental and/or geochemical variations and if the produced signals are largely the adaptive response by Thaumarchaeota to oceanographic effects associated with climate or seasonal temperature changes in the upper water column. We present the results of a four push-core transect study of surface sediments collected along an environmental gradient at the Cathedral Hill hydrothermal-vent system of The transect crosses a region where advecting hydrothermal ﬂuids reach 155 ◦ C within the upper 21 cm below the seaﬂoor (cm b.s.f.) close to the vent center to near-ambient conditions at the vent periphery. The recovered iGDGTs closest to the vent center experienced high rates of turnover with up to 94 % of the lipid pool being lost within the upper 21 cm b.s.f. Here, we show that the turnover is non-selective across TEX 86 GDGT lipids and does not affect the ratio independently. However, by TEX 86 highly Hill

Abstract. The diversity and relative abundances of tetraether lipids produced by archaea and bacteria in soils and sediments are increasingly used to assess environmental change. For instance, the TetraEther indeX of 86 carbon atoms (TEX 86 ), based on archaeal isoprenoidal glycerol dialkyl glycerol tetraether (iGDGT) lipids, is frequently applied to reconstruct past sea-surface temperatures (SSTs). Yet, it is unknown how the ratio fully responds to environmental and/or geochemical variations and if the produced signals are largely the adaptive response by Thaumarchaeota to oceanographic effects associated with climate or seasonal temperature changes in the upper water column. We present the results of a four push-core transect study of surface sediments collected along an environmental gradient at the Cathedral Hill hydrothermal-vent system in Guaymas Basin, Gulf of California. The transect crosses a region where advecting hydrothermal fluids reach 155 • C within the upper 21 cm below the seafloor (cm b.s.f.) close to the vent center to near-ambient conditions at the vent periphery. The recovered iGDGTs closest to the vent center experienced high rates of turnover with up to 94 % of the lipid pool being lost within the upper 21 cm b.s.f. Here, we show that the turnover is non-selective across TEX 86 GDGT lipids and does not affect the ratio independently. However, as evident by TEX 86 ratios being highly correlated to the Cathedral Hill vent sediment porewater temperatures (R 2 = 0.84), the ratio can be strongly impacted by the combination of severe lipid loss coupled with the addition of in situ iGDGT production from archaeal communities living in the vent sediments. The resulting overprint produces absolute temperature offsets of up to 4 • C based on the TEX H 86 calibration relative to modern climate records of the region. The overprint is also striking given the flux of iGDGTs from the upper water column is estimated to be ∼ 93 % of the combined intact polar lipid (IPL) and core GDGT lipid pool initially deposited on the seafloor. A model to correct the overprint signal using IPLs is therefore presented that can similarly be applied to all near-surface marine sediment systems where calibration models or climate reconstructions are made based on the TEX 86 measure.

Introduction
Archaeal and bacterial tetraether cellular membrane lipids represent a group of common and structurally diverse compounds frequently used to track the presence of living and dead microorganisms as well as geochemical and physical conditions within present-day and paleoenvironments (e.g., Schouten et al., 2002Schouten et al., , 2004Schouten et al., , 2013Sturt et al., 2004;Hopmans et al., 2004;Weijers et al., 2007Weijers et al., , 2014Hollis et al., 2012;O'Brien et al., 2017). In this regard, the proportional abundances of these lipids form various prominent proxies for assessing environmental change through time. For example, TEX 86 (TetraEther indeX of 86 carbon atoms; Schouten et al., 2002) is a widely used archaeal-lipid-based paleotemperature proxy for marine environments. The ratio measures variations in the number of cyclopentyl rings for a select group of archaeal core lipids (CLs) (Fig. S1 in the Supplement) following the assumption that biphytanyl cyclization is an organismal response to changing sea-surface temperatures (SSTs). The proxy is therefore used in many regions around the world with TEX 86 values typically ranging from 0.2-0.9 in marine settings (e.g., Huguet et al., 2006;Kim et al., 2008;McClymont et al., 2012;Tierney, 2014). The utility of TEX 86 rests on the premise that iGDGTs found in ocean bottom sediments are almost exclusively produced by marine planktonic archaea that inhabit the epipelagic zone Tierney, 2014;Besseling et al., 2019Besseling et al., , 2020. Lipids are therefore required to be efficiently and continually transported from the upper water column to the underlying ocean floor to produce a fossil chemostratigraphic record of microbial response to changing SST conditions with time (Wuchter et al., 2005).
Since its introduction, the reliability of TEX 86 to accurately track paleoclimate variations has been questioned. TEX 86 -based SST estimates have been observed to substantially deviate from other temperature proxies (e.g., Huguet et al., 2006;Rommerskirchen et al., 2011;Seki et al., 2012Seki et al., , 2014. For example, over the past decade, considerable effort has been made to reconstruct the early Paleogene greenhouse climate system. However, TEX 86 appears to significantly overestimate reconstructed SSTs (Hollis et al., 2012) relative to other proxies such as Mg/Ca or clumped isotopic compositions of foraminiferal calcite, as well as various climate models based on partial pressure of carbon dioxide (pCO 2 ) predictions (Lunt et al., 2012;Naafs et al., 2018). For late Neogene climate reconstructions, TEX 86 has been shown to underestimate warming trends relative to the U k 37 alkenone-index-derived (Brassell et al., 1986) temperatures (Lawrence et al., 2020). The apparent SST offsets have been attributed to how the proxy's associated lipids change in relation to their environment and if these changes are regulated by internal adaptations within the archaeon or by an overarching community succession. In this regard, the debate surrounding these discrepancies largely centers on establishing responses to seasonal biases (e.g., Herford et al., 2006;Wuchter et al., 2006;Huguet et al., 2011); the development of adequate calibration methods (e.g., Kim et al., 2010;Pearson et al., 2013;Tierney et al., 2014); identifying lipid sourcing effects, including subsurface sediments origins for those used with the calculation of TEX 86 (e.g., Lipp and Hinrichs, 2009); and physical, chemical, and ecological controls for archaeon iGDGTs cyclization (e.g., Elling et al., 2015;Qin et al., 2015;Hurley et al., 2016).
For non-thermal influences, the primary concern is what archaeal taxa produce iGDGTs and where they are sourced. To this end, most TEX 86 lipids are thought to be produced by Marine Group I (MGI) planktonic Thaumarchaeota (Brochier-Armanet et al., 2008), which are most abundant below the photic and epipelagic zone (e.g., Karner et al., 2001). Within this context, many regions of the ocean floor may become highly impacted by colder, deeper water column inputs (e.g., Karner et al., 2001;Huguet et al., 2007;Lopes dos Santos et al., 2010;Kim et al., 2012aKim et al., , b, 2015Pearson et al., 2013;Ho and Laepple, 2016;Hurley et al., 2016;Lui et al., 2018;Sinninghe Damsté et al., 2018). Other non-thermogenic driving forces impacting the production, cyclization, and relative abundance of TEX 86 -based lipids include organismal selectivity to specific growth phases and growth rates (Elling et al., 2014;Hurley et al., 2016), redox conditions (Qin et al., 2015), and the incorporation of iGDGT (isoprenoidal glycerol dialkyl glycerol tetraether) from archaeal communities living in the ocean floor sediments. With respect to the latter, Lipp and Hinrichs (2009) demonstrated that the production of intact polar lipid GDGTs ( IPL GDGTs) by ocean floor sediment microbial communities collected in the Peru Margin were distinctly different from upper-water-column-sourced CLs and that the conversion of this living pool to fossil lipids would shift TEX 86 ratios to higher values. However, the overall impact may not be substantial, as Umoh et al. (2020) found little effect to the TEX 86 paleoclimate ratio when examining surface sediments near hydrothermal-vent sites on the Southeast Indian Ridge in the southern Indian Ocean. Lengger et al. (2012Lengger et al. ( , 2014 also reported no significant deviation between the TEX 86 values in sediment cores collected near the oxygen minimum zone from that of the overlying water column in the Arabian Sea with near-linear degradation rates of both IPLs and CLs. All together, the iGDGT abundances recorded in a TEX 86 sediment value may ultimately constitute a multi-variable data point -mixing lipid components that are themselves responses to temperature, organismal substrate, and metabolism dynamics, biozone niche partitions spanning from the ocean surface to the shallow sediment archaeal community, which ultimately become further attenuated by depositional and diagenetic processes. While not an ideal location to create SST reconstructions, hydrothermal vents of sedimented ocean basins do represent an anomalous endmember to the vast expanse of ambient ocean floor sediment where paleoclimate reconstructions are commonly produced. The sedimented vent systems of Guaymas Basin, Gulf of California (Fig. 1), is one such site. The basin experiences high sedimentation rates ranging from 0.4-0.2 cm yr −1 (Gieskes et al., 1988) due in part to the high productivity of the upper water column. The ocean floor hydrothermally impacted surface sediments are also a location of active and diverse microbial communities with vents that are often covered by Beggiatoa-dominated microbial mats (e.g., McKay et al., 2012;Meyer et al., 2013;Teske et al., 2016). These sites should, in principle, enable a high-resolution archaeal-lipid stratigraphic record that provides optimal conditions for studying potential shallow diagenetic and subseafloor interferences to common archaeallipid-based environmental proxies. The region further offers an ideal setting to compare TEX 86 proxy responses to in situ lipid production from thermophilic sedimentary ar-chaea that differ from the pelagic background communities (e.g., Schouten et al., 2003). Recently, Bentley et al. (2022) produced a survey of the source and diagenetic and catagenetic alteration of archaeal lipids from the Cathedral Hill hydrothermal-vent complex ( Fig. 1) in Guaymas Basin, Gulf of California. Within the investigation, it was observed that most iGDGTs are sourced from the overlying water column. Building on the results of Schouten et al. (2003), it was observed that these lipids can become heavily turned over in the hotter portions of the vent site where they rarely survive long enough to become cracked into hydrocarbon biomarkers such as biphytanes and derivatives of biphytanes. For this study, we further examine the iGDGT lipid distributions in these near-surface ocean floor sediments to determine if paleoclimate proxy signals can be impacted by the presence of subsurface archaeal populations. The distribution of iGDGTs and their corresponding environmental proxy signals were measured within the sediments along a transect at the complex. In this regard, this site offers the unique opportunity to evaluate the response of TEX 86 and other tetraether-lipid proxies within a microbially diverse sedimentary environment that is exposed to high-temperature vent fluids.

Study location and sampling
Four sediment push cores were collected using HOV Alvin (dive 4462; 22 October 2008; human-occupied vehicle) at the Cathedral Hill hydrothermal-vent site, located at a water depth of 1996 m in the southern trough of Guaymas Basin, Gulf of California (27 • 0.629 N, 111 • 24.265 W) (Fig. 1). The push cores, labeled 1 to 4, were taken along a transect with ∼ 2 m spacing extending outwards from microbial-matcovered sediments near the sulfide chimney complex to just outside of the microbial-mat area in ambient seafloor sediment. Thermal-probe measurements were sequentially taken beside each core (Table 1). Once the push cores were brought to the surface, the sediments were subsampled into 2-3 cm thick depth intervals, transferred to combusted glass vials, and immediately stored at −40 • C (on board the ship) before being shipped under dry ice to the laboratory and later freeze-dried and stored at −80 • C.

Lipid extraction
Lipid extractions followed a modified Bligh and Dyer protocol laid out in Bentley et al. (2022) and following Sturt et al. (2004). A subsample of freeze-dried sediment was added to a Teflon© centrifuge tube followed by the addition of 6 ml of mix A solvent solution comprising of 2 : 1 : 0.8 v/v/v methanol (MeOH), dichloromethane (DCM), and phosphate buffer (5.5 g L −1 Na 2 HPO 4 ; Avantor Performance Materials; adjusted to pH of 7.4 with HCl; Anachemia Canada Inc.). The solvent sediment mixture was further spiked with 1-alkyl-2-acetoyl-sn-glycero-3-phosphocholine (PAF) recovery standard purchased from Avanti Polar Lipids, Inc. The slurry was sonicated for 5 min then centrifuged for 5 min at 1250 rpm. The resulting supernatant was added to a separatory funnel. This procedure was performed twice before being joined by two replicate extractions using mix B, a 2 : 1 : 0.8; v/v/v solution of MeOH, DCM, and trichloroacetic acid buffer (50 g L −1 C 2 HCl 3 O 2; Avantor Performance Materials, LLC; pH of 2), and finally two replicate extractions using mix C, a 5 : 1 v/v solution of MeOH and DCM. Once complete, this combined A, B, and C. For each step, the organic fraction was collected in a beaker, and the combination of mix A, B, and C was subjected to 10 ml of DCM and H 2 O (Milli-Q) to achieve separation. The organic phase was drawn off, and the water was extracted using three DCM washes, drawing off the organic phase after each wash. The organic phase was back-extracted with H 2 O to ensure purity. The resulting organic phase was then evaporated to dryness at 60 • C under dry nitrogen. The resulting total lipid extract (TLE) was transferred to pre-weighed autosampler vials using DCM : MeOH 1 : 1 v/v, spiked with 1, 2-diheneicosanoyl-sn-glycero-3-phosphocholine (C 21 -PC; Avanti Polar Lipids, Inc.) and stored at −20 • C.

High-performance liquid chromatography-mass spectrometry (HPLC-MS)
Mass spectrometric analyses were performed on an Agilent Technologies 1260 Infinity II HPLC coupled to an Agilent Technologies 6530 quadruple time-of-flight mass spectrometer (qToF-MS) operated in positive mode. Chromatographic separation used a reverse-phase method outlined by Zhu et al. (2013). The HPLC was fitted with an Agilent Technologies ZORBAX RRHD Eclipse Plus C18 (2.1 mm × 150 mm × 1.8 µm) reverse-phase column and guard column maintained at 45 • C. The sample injection solvent was methanol. An aliquot of each sample representing 1 % of the TLE was analyzed. A 0.25 mL min −1 flow rate was established with mobile phase A consisting of methanol / formic acid / ammonium hydroxide (100 : 0.04 : 0.10 v/v/v) held at 100 % for 10 min, thereafter mixed following a linear gradient with mobile phase B of propan-2-ol / formic acid / ammonium hydroxide (100 : 0.04 : 0.10 v/v/v) to 24 %, 65 %, and 70 % over 5, 75, and 15 min intervals, respectively. Each sample run was finished by re-equilibrating the system with 100 % mobile phase A for 15 min. The effluent was ionized by an electrospray ionization source with a gas temperature of 300 • C, a 3 L min −1 drying gas flow, and a 5.33 µA source current. The mass spectrometer was set to a 100-3000 m/z scan range in positive mode in an untargeted method with 10 ppb resolution to simultaneously resolve both archaeal IPLs and CLs. Analyte identification was achieved by accurate mass resolution, mass spectral analysis using Agilent Technologies's MassHunter software, and comparison of fragmenta- tion patterns with the literature (e.g., Knappy et al., 2009;Liu et al., 2010;Yoshinaga et al., 2011;see Bentley et al., 2022, for further details). Mass fragments consistent with the loss of a biphytane (m/z 743.7) were screened for all archaeal lipids. Quantification was achieved by summing the integration peak areas of [M+H] + , [M+NH 4 ] + , and [M+Na] + adducts for the respective IPLs and CLs of interest. Concentration values were obtained relative to the internal C 21 -PC standard and reported in microgram per gram dry sediment weight. Response factors were determined by a series of injections of a standard solution containing PAF; C 21 -PC; 1,2-diacyl-3-O-(α-D-galactosyl1-6)-β-D-galactosyl-sn-glycerol (DGDG); 1,2-diacyl-3-O-β-D-galactosyl-sn-glycerol (MGDG); 1,2-di-O-phytanyl-snglycerol (archaeol); 1 ,3 -bis[1,2-dimyristoyl-sn-glycero-3phospho]-glycerol (14 : 0 cardiolipin) from Avanti Polar Lipids, Inc., USA; and 2,2 -di-O-decyl-3,3 -di-O-(1 ,ωeicosanyl)-1,1 -di-(rac-glycerol) (C 46 -GTGT) from Pandion Laboratories, LLC in amounts ranging from 100 pg to 30 ng. Response factors were calculated relative to the C 21 -PC, and the appropriate correction factor was then applied to the lipid class of interest.
A series of samples were re-run to identify or confirm deviations in the data set. The variations between the concentrations of GDGTs in the re-run and the initial runs yielded a maximum difference of ∼ ± 4 µg g −1 per GDGT compound, providing confidence in the initial results and confirming the presence of two outliers in the data set (Bentley et al., 2022). These outliers are core 4 at 8-10 cm, with abnormally low concentrations of all compounds that are likely ion suppres-sion from a sample heavily impregnated with oil, and core 3 at 15-18 cm, which contains relatively high lipid concentrations that are yet to be explained.

Archaeal-lipid diversity and turnover
The Cathedral Hill transect sediments have iGDGTs containing 0-4 cyclopentyl (GDGT-0-4) as well as crenarchaeol (Cren) and the isomer of crenarchaeol (Cren') that contains five rings (four cyclopentyl and one cyclohexyl moiety) (Table S1). Branched GDGTs (brGDGTs) including Ia-c, IIa-c, and IIIa were found to have discontinuous and/or low absolute abundances, with some compound classes not being detected (i.e., brGDGT-IIIb; Table S2). The brGDGTs are therefore not further examined in this study. For cores 1 to 3 the concentrations of nearly all iGDGT compounds systematically decrease with depth ( Fig. 2). Bentley et al. (2022) established the sedimentation of archaeal lipids from the upper water column as being uniform both in terms of spatial loading across the length of the transect as well as over an inferred 52.5-105 years of sedimentation as penetrated by the length of the push core (based on sedimentation rates). From this, it is estimated that ∼ 70.6 ± 23.5 µg iGDGT g −1 sed yr −1 is being deposited on the seafloor from the overlying water column. However, for cores closest to the vent site, lipid abundances exhibited a much sharper decrease with depth, which Bentley et al. (2022) attribute to the turnover of archaeal lipids cou- pled to, but not directly caused by, hydrothermalism. For cores 1 and 2, losses reach as high as 94 % within the upper 21 cm b.s.f. (centimeters below seafloor). The lipid loss is less severe for core 3 at ∼ 60 %. For the ambient core 4, iGDGTs have similar downcore stratigraphic trends with a near-consistent average of 400 µg g −1 sed concentration and no systematic loss of lipids. Due to the high-temperature conditions of the vent fluids at Cathedral Hill, the identified archaeal iGDGT-based IPLs within the sediments most likely represent the composition of cellular membrane material from archaeal communities living in the sediments. These lipids have exclusively monoglycosyl (1G) or diglycosyl (2G) head groups linked to a 2,3sn-glycerol. Within the pyrolytic environment, the transformation of IPL iGDGTs could hypothetically add to the core iGDGT lipid pool. Similar to CLs, the 1G-GDGTs contain 0-4 cyclopentyl moieties and include Cren and Cren'. Surface concentrations of these lipids are ∼ 15 µg g −1 sed in cores 1 to 3 (residing within the microbial mat) and 11 µg g −1 sed for core 4 (Table S2). Also similar to the CLs, the archaeal IPL concentrations decrease downcore and are closely coupled to increasing porewater temperatures (Table S2). For cores 1 and 2, the maximum depths for detectable 1G-GDGTs are 15-18 and 12-15 cm b.s.f., corresponding to vent porewater temperatures of 145 and 87 • C, respectively. In core 3, 1G-GDGTs persist downcore with a consistent lipid depletion that reaches its lowest concentration of 5.2 µg g −1 sed in the bottom of the core at 18-21 cm b.s.f. sediment depth where porewater temperatures rise to 80 • C. In core 4, which is most similar to the ambient ocean bottom conditions and falls outside of the area covered by the microbial mat, the lipid concentrations average is ∼ 8 µg g −1 sed across the depth of the core. The 2G-GDGTs have zero to two cyclopentyl rings that for cores 1 and 2 are restricted to the upper 4 to 6 cm b.s.f. These lipids are not further investigated in this study, as 2G-GDGTs are of limited abundance (max summed concentrations <2 µg g −1 sed) and their structural diversities negligibly affect isoprenoid-based proxies.  Bentley et al. (2022). 2 Sum of GDGT-1, -2, -3. -4, -5, and -5 (Table S1). 3 Sum of all detected 1G-and 2G-GDGTs (Table S3). a Collected core numbers are relabeled in the sample name to reflect a relative transect position (1-4). b Sediment lithology based on freeze-dried sediments. c TEX 86 = (GDGT-2 + GDGT-3 + GDGT-5 )/(GDGT-1 + GDGT-2 + GDGT-3 + GDGT-5 ) (Schouten et al., 2002) applied to both core GDGTs and 1-glycosyl-GDGTs (also referred to as M TEX 86 in Sect. 3.4). d TEX H 86 = log ((GDGT-2 + GDGT-3 + GDGT-5')/(GDGT-1 + GDGT-2 + GDGT-3 + GDGT-5')) for sediments outside low latitudes (Kim et al., 2010).   (Bentley et al., 2022). The gray background indicates transect regions where porewater temperatures exceeded 123 • C, marking the known upper thermal limit of life (Kashefi and Lovley, 2003). Yellow fields are zones where oil generation and hydrocarbon degradation occur (Dalzell et al., 2021).
Lipid-based proxies for the calibration or reconstruction of paleoclimate records such as TEX 86 are based on environmentally scaled contributions of select GDGT compounds. These proxies could be negatively impacted should other ocean floor sediment systems experience high rates of lipid turnover (Lengger et al., 2014). To evaluate whether downcore depletions of lipid concentrations impacted tetraetherbased proxies, the concentrations of the highly abundant GDGT-0 was plotted relative to the TEX 86 ratio lipids (iGDGT-1, iGDGT-2, iGDGT-3, and Cren') (Fig. 3a). For Fig. 3a, straight lines in the logarithmic plot indicate nearequal depletion rates between the paired xand y-axis lipid classes. Similarly, parallel slopes for the various lipid pairs also indicate near-equal depletion rates, with vertical offsets between pairs marking different initial starting abundances of the compared lipid. In this regard, iGDGT-0, iGDGT-1, iGDGT-2, and Cren' have undergone the same rate of turnover. However, the depletion rate of iGDGT-3's is lower than that of other lipid classes for cores 1 and 2. Although, this may represent a distinct resilience to turnover, we suggest it instead results from overprinting by the subsurface hyperthermophilic-archaeal community (see below).
To better track changes across each core, the degradation rate constants (k ) of TEX 86 lipid classes were calculated for each push core ( Fig. S2; Table S3) using a first-order kinetic model: in which C t and C i are concentration at time (t) and the initial concentration, respectively (e.g., Schouten et al., 2010). Rearranging Eq. (1), the k values were calculated as from which data, it is evident that the downcore concentrations of each lipid decrease at equivalent rates (i.e., they have the same slopes for their rates of decay; s 2 = 0.2). The exception to this is core 2, which independent of two outliers has different decay paths for GDGT-3 and GDGT-5. This is consistent with the TEX 86 iGDGT lipid classes largely being removed from the sediment lipid pool in a non-selective manner.
Based on these results, the TEX 86 , ring index (RI), and methane index (MI) values were plotted against their respective summed iGDGTs lipid concentrations (Fig. 3b-d). For samples located within the habitable zone (having porewaters ranging from 0-123 • C; Kashefi and Lovley, 2003), no correlation is observed between the lipid abundances and proxy ratios of TEX 86 , RI, or MI ( Fig. 3b-d). This further suggests these proxies are not affected by turnover in the habitable zone. However, once sediment burial reaches beyond the habitable zone, TEX 86 ratios trend to higher values (similarly also reflected in GDGT-3 concentration trends of Fig. 3a). Collectively, these data strongly indicate that archaeal-lipid turnover is largely nonselective of the TEX 86 lipid classes and will therefore theoretically not in and of themselves significantly impact archaeal-lipid paleoclimate proxy reconstructions.
Apart from paleoclimate reconstructions, archaeal-lipid CLs are sometimes used to resolve aspects of localized biogeochemical cycles within sediments. To this end, the location and degree of anaerobic oxidation of methane (AOM) is determined by methane and archaeal-lipid carbon isotope measures (e.g., Boetius et al., 2000;Schouten et al., 2003;Stadnitskaia et al., 2008;Biddle et al., 2012) as well as by the proportional abundances of core GDGTs (cGDGTs) in the form of the MI (Zhang et al., 2011;Carr et al., 2018;Petrick et al., 2019). With respect to the latter, the MI proxy is used to differentiate regions of normal marine (with values between 0-0.3) and active AOM conditions in and around cold seeps (where values > 0.5-1 are reported for gas hydrate impacted sediments and subsurface environments with high AOM levels). To our knowledge, the use of this proxy for hydrothermal-vent systems has not been thoroughly investigated even though this microbial process has been well documented at Guaymas Basin. For example, highly 13 C-depleted CLs reaching up to −70 ‰ in hydrothermal-vent sediments with porewater temperatures as high as 95 • C indicate thermophilic archaea actively engaging in AOM . Biddle et al. (2012) through the detection of relevant archaeal communities by 16S RNA in conjunction with highly depleted methane carbon isotope values determined active AOM spanning 35 to 90 • C porewater conditions. AOM is not likely to be the dominant form of carbon and sulfur metabolism, as it generally accounts for less than 5 % of sulfate reduction (Kallmeyer and Boetius, 2004). When applying the MI to the Cathedral Hill push-core transect survey, low values (ranging from 0.2-0.38; Table 1) are recorded with no correspondence to thermal controls across the vent transect (Fig. 4). Although, it could be considered that the low values arise from a lack of AOM within these sediments, the low MI values are consistent with a high upper-watercolumn iGDGT loading as estimated by Bentley et al. (2022).

TEX 86 and reconstructed SSTs
McClymont et al. (2012) reported a GDGT-based reconstructed annual SSTs of 16-18 • C from particulate organic matter collected in ambient sediment traps in Guaymas Basin during an annual cycle from 1996-1997. The reconstructed temperatures followed the calibration model for sediments outside of polar regions proposed by Kim et al. (2010). These authors demonstrated the temperatures derived from the TEX 86 reconstruction were significantly lower than those produced by the closely co-varying U k 37 paleoclimate proxy and satellite-measured estimates that jointly estimated a mean annual sea-surface temperature (MASST) of 23 • C. The longer 21-year (1982-2004) satellite-derived MASST is also reported to be higher at 24 • C (Herrera- Cervantes et al., 2007). Although, the sites and timeframes of these surveys do not match that of the Cathedral Hill survey, they do provide context to what our reconstructed TEX 86 values should record.
The high sedimentation rate at Cathedral Hill has resulted in near-homogenous inputs of organic matter from the upper water column across the transect area (Dalzell et al., 2021;Bentley et al., 2022). Therefore, TEX 86 reconstructions should produce equivalent cross-transect trends with sediment depth. Nonetheless, as with changes in the archaeal-lipid concentrations, the profiles of iGDGT proxies TEX 86 and RI of the transect similarly have downcore trends ( Fig. 2; Bentley et al., 2022). For core 4, TEX 86 spans a narrow range of values (n = 7; 0.52-0.54, average: 0.53 ± 0.01; Fig. 4a) across a period of ∼ 37.5 to 75 years, corresponding to the depth of the cores. To a slightly lesser degree, the core top (0-2 cm b.s.f.) across the transect also display near-equal values to core 4 (n = 4; 0.56-0.54; average: 0.55 ± 0.01). These values mark a TEX H 86 reconstructed mean annual SST of 19.3-20.4 • C following the Kim et al. (2010) calibration model (Table 1) Fig. 4a), corresponding to a TEX H 86 reconstructed SST change of 3.1 • C, marking a range from 21.8 to 24.9 • C ( Table 1). The fundamental driver for the proxy is likely influenced by the archaeal community composition  Kim et al. (2010). Blue field indicates MI values for normal marine conditions (Zhang et al., 2011. that is responding to their exposure to in situ vent fluid temperatures (Fig. 4).
Two mechanisms are considered for the observed proxy variations. The first is that progressive ring loss due to carbon-carbon bond cleavage of pentacyclic rings moieties by exposure to the sharp geothermal gradient acts to systematically attenuate the iGDGT lipid pool. Hydrous pyrolysis experiments conducted by Schouten et al. (2004) demonstrated that at extreme temperatures (ca. > 160 • C), TEX 86 values become negatively impacted by the preferential destruction of polycyclic GDGTs. Such losses produce progressively lower ratio values. Although the transect sediment porewaters do not reach the pyrolytic temperatures of the Schouten et al. (2004) experiment, they are high enough to generate hydrocarbons (Dalzell et al., 2021) and thermochemically degrade iGDGTs in the hottest regions of the transect; they are also more long-lived than what is produced from a laboratory experiment. However, the observed strati-graphic TEX 86 trends do not match those of predicted ring loss, as the values increase rather than decrease in relation to elevated porewater condition. Nonetheless, the thermochemical oxidative loss of GDGTs and its effect on the TEX 86 ratio is further explored below (Sect. 3.4).
The second mechanism is that subsurface microbial communities donate enough core GDGTs to overprint the detrital signal source. The RI (Fig. 4b) values were similarly compared to recorded porewater temperatures to better interpret the TEX 86 trends and to ensure that the Cathedral Hill reconstructed temperatures are influenced by the subsurface microbial community. In this regard, RI is used to monitor the adaptive response of an archaeal community at the hydrothermal-vent site. Lipid cyclization is an adaptive response to changing environmental temperature or acidity in which an archaeon increases its rigidity by decreasing the fluidity and permeability of its cellular membrane that, therefore, also further regulates the flow of solutes and nutrients in and out of the cell (Gliozzi et al., 1983;De Rosa and Gambacorta, 1988;Uda et al., 2001;Schouten et al., 2002;Macalady et al., 2004;Boyd et al., 2013). Both cores 1 and 2 have RI values highly correlated to temperature (R 2 = 0.87 and 0.75, respectively), consistent with heat stress adaption. This same was also observed in Guaymas Basin by , who reported an increase in the RI of core lipid GDGTs with in situ temperature. As such, a significant proportion of the measured iGDGTs likely emanates from archaeal communities living in the shallow sediments of Cathedral Hill. As such, the lipid cyclization pattern may reflect stratigraphically discrete thermophilic to hyperthermophilic communities that are selectively adapted to more extreme temperature conditions (see Bentley et al., 2022, for further discussion on the lipid-based taxonomic makeup of the vent site).

Lipid signal sourcing
To evaluate the sources of measured archaeal lipids, CL and IPL TEX 86 (the ratio applied to IPLs that contain equivalent core lipids) indices were compared as signal responses from their respective pools of living and dead cellular debris (Fig. 5). For cores 1, 2, and 3 the 1G-iGDGT IPL TEX 86 measures are positively correlated with temperature (R 2 = 0.46, 0.74, and 0.66, respectively; Fig. 5a). In this regard, 1G-iGDGT IPL TEX 86 ratio appears to be largely influenced by in situ porewater temperatures as well as potentially the archaeal community ecology of the vent system. Factors such as community composition and adaptation may further impact the IPL TEX 86 ratio as the rates of changes between cores 1-3 are not the same. Similar to the CL TEX 86 values, the IPL TEX 86 is not correlated to their summed TEX 86 lipid abundances (Fig. 5b). Such a condition is largely consistent with the living lipid pool being modified by the archaeal community's response to thermal stress and not by subse-quent thermal-oxidative transformation occurring shortly after cell death.
The IPL and CL lipids of transect samples can be further grouped into three clusters (A, B, C), suggesting a mixed signal for the sourcing of archaeal GDGTs from both the living and dead pools of archaea (Fig. 5c), closely tracking temperature. In this plot, we assume that clusters falling on the 1 : 1 line indicate the living biota can equally contribute to the dead pool of total recovered GDGTs. Those off axis contribute either less or more to one or the other lipid pool. The three clusters mark unique thermal zones within the transect area, with cluster A being composed of the ambient core 2 to 4 seafloor surface samples, cluster B marking a mix of intermediate temperature samples from all cores, and cluster C being composed of high-temperature samples. The lipid groups likely mark distinct archaeal communities. As cluster B resides on the 1 : 1 line, the TEX 86 core lipids likely have a mix of detrital and in situ inputs. Cluster C, however, appears likely dominated by in situ lipid production. The thermal zonation and equivalent directionality of the resulting ratios (i.e., both CL and IPL TEX 86 ratios increase with porewater temperature) further support overprinting of the original CL TEX 86 sea-surface signal by the ocean bottom sediment archaeal community as a mechanism for the observed CL TEX 86 trends.
Collectively, these results suggest the source of the archaeal CLs measured in the TEX 86 and RI indices progressively become more dominated by subsurface microbial communities adapted to the hotter hydrothermal-vent fluids. Our results also indicate that in select natural environments, such as hydrothermal-vent complexes, the TEX 86 SST proxy may entirely record ocean bottom sediment porewater temperatures. To our knowledge, a clear case of overprinting to this level has not yet been demonstrated.

TEX 86 overprint corrections
The measured TEX 86 ( M TEX 86 ) value of the Cathedral Hill sediments is herein considered to be a weighted sum of a sea-surface TEX 86 ( SS TEX 86 ) value acquired from lipids sourced in the upper water column that is further modified by a component of the deeper-water-column-sourced core lipids ( WC TEX 86 ) as well as by additions of archaeal lipids from the benthic and subsurface microbial communities ( Sed TEX 86 ). These ratio loadings are collectively also potentially further modified by diagenetic influences in the ocean bottom sediments. Over the cumulative sediment burial period and in consideration of the measured porewater temperatures of the Cathedral Hill push-core sediments, these influences include the selective loss of lipids by their binding into protokerogen (K) and by potential changes due to the loss of lipid by turnover (ϕ; Sect. 3.1 ϕ). Additional catagenetic effects from thermochemical alteration of lipids (θ ) may also attenuate the sum of sedimentary core lipids by their exposure to high-temperature vent fluids. Collectively, these ef-fects are considered to form the following relationship: where a, b, and c, are measured scaling parameters for lipid loading and ϕ, K, and θ are diagenetic and catagenetic alteration parameters. Solving for SS TEX 86 , In this regard, a portion of the archaeal community from the upper water column, presumably initially sourced of IPLs, and an additional community inhabiting the ocean floor sediments were assumed to eventually die, with their respective IPLs gradually hydrolyzing, joining the CL pool where they further contribute to the observed M TEX 86 value. For this study, no data were collected to calculate b WC TEX 86 , and its potential impact on M TEX 86 cannot be further considered in this study. However, it is highly likely, given the longer residence times for glycosidic-based head groups of the identified archaeal IPLs and their relatively short settling time through the water column (Lengger et al., 2012;Xie et al., 2013), that a component of this lipid source was already mixed with the Sed TEX 86 contribution. For this study, Sed TEX 86 is an IPL TEX 86 ratio based on detected 1G-GDGT-1, 1G-GDGT-2, 1G-GDGT-3, Cren', 2G-GDGT-1, and 2G-GDGT-2, as present in the original paleoclimate proxy (Table 1; Fig. 6). Testing the removal of 2G-GDGT lipids, which have a low absolute concentration (< 2 µg g −1 sed) and shallow stratigraphic zones of occurrence (Sect. 3.1; Table S2), yielded a negligible <1 • C change in the summed average reconstructed SST.
using the accumulated summed concentrations of 1G-GDGTs and 2G-GDGTs that have the potential to become converted to cGDGTs by progressive burial diagenesis and d 0−n marking the range of sampled sediment depths, with 0 being the 0-2 cm b.s.f. core top and n being the deepest point of sediment burial. These intervals are divided by the water column input of TEX 86 lipids ([GDGTs CL− TEX86 lipids ] 0−2 cm ) estimated to be 120 µg g −1 sed based on their average measured concentration across the four-core transect. The function assumes the surface sediment does not hydrolyze its IPL GDGTs to CLs (Table 2). When applied to Eq. (4) and further excluding ϕ, K, and θ , the SS + WC TEX H 86 reconstructed SSTs average 19.68 ± 0.79 • C (Table 2; Fig. 6a), with the total samples having an unchanging depth profile that mirrors the range of values measured in the ambient sediments of core 4 (Fig. 2).   (c) TEX 86 proxy for core GDGTs versus 1G-GDGTs. Clusters A-C may represent different archeal communities that are providing varying inputs of iGDGT to the core GDGT lipid pool. The dotted trend line is the partial least-square regression of the complete core lipid TEX 86 data set. The solid line marks the 1 : 1 CL-to-IPL proxy correspondence indicating both allochthonous and autochthonous sources contribute equally to the core GDGT lipid pool.
The selective lipid removal by diagenetic and catagenetic processes theoretically may also affect the TEX 86 value; however, their perspective impact on the directionality and magnitude of the ratio are difficult to predict and equally hard to discretely measure. Although the loss of GDGTs to protokerogen formation could potentially impact the ratio, it was shown to be a negligible sink for the lipids (Bentley et al., 2022). As such, the K parameter in Eqs. (3) and (4) was therefore assigned a 0 value. Due to the high geothermal gradient at Cathedral Hill, some of the transect push-core sediments resided within zones of active catagenesis ( Fig. 2; Dalzell et al., 2021). The degradation rates of each TEX 86 lipid were independently measured for the four push cores (Eq. 2; Fig. S2). As the abundance of both CLs and IPLs differentially decreases through the various core sediment profiles with turnover rates that appear to be constrained by porewater temperature changes (Sect. 3.1), the degradation rates must also record the effects of thermochemical oxida-tive weathering (Fig. 3b). In this case, ϕ and θ are treated as grouped parameters. To determine if individual lipid classes were selectively removed during degradation, the variance (s 2 ) of the rate change as measured from its respective regression slope (i.e., m logk ) from the TEX 86 lipids ( Fig. S2; Table S4 from Eq. 2) was calculated. For the Cathedral Hill transect, the calculated m logk s 2 is 0.20, which is due to accelerated degradation rates for higher-ring lipids, GDGT-3 and Cren', in samples from cores 1 and 2, where high vent temperatures resulted in hydrocarbon generation of the sediments (Dalzell et al., 2021). A weighing function for the degree of lipid class selectivity during turnover is proposed: When applied to Eq. (4), the corrected data series produces an average transect SS + WC TEX H 86 reconstructed SST of 23.66 ± 0.59 • C with a near-zero partial least-square re- gression slope (Table 2; Fig. 6b). These modeled values are within the 23-24 • C obtained for the 21-year (1982-2004) satellite-derived MASST data for the Guaymas Basin region (Herrera-Cervantes et al., 2007). Based on these calculations, nearly all M TEX 86 attenuation can be attributed to sediment microbial overprinting coupled to diagenetic and catagenetic loss of lipids consistent with prior observations at Guaymas Basin Zhang et al., 2011). The high degree of influence this has on the TEX 86 proxy is striking given that the upper water flux of GDGTs at Cathedral Hill is estimated to represent up to 93 % of the total intact polar and core GDGT lipid pool within these sediments. Although, this study demonstrates the benthic microbial community can influence TEX 86 values in anomalous, endmember environments, the above model has not yet been tested across conventional ocean shelf environments.

Conclusions
In this study, we demonstrate a pronounced overprint of cGDGTs sourced from the ocean floor sedimentary archaeal community at the Cathedral Hill vent site in Guaymas Basin. The overprint is marked by lipids with more cyclized ring moieties, marking an adaptive response by archaea to rigidify the cellular membranes against localized heat stress. This in turn has resulted in the commonly used TEX 86 paleoclimate proxy partially recording advecting porewaters temperatures. As the vast majority of cGDGTs in these sediments is sourced from the overlying water column, the impact on the TEX 86 ratio is further the product of rapid lipid turnover rates and diagenetic and catagenetic alteration processes potentially unique to the hydrothermal system. Together, these factors resulted in absolute TEX H 86 temperature offsets of up to 4 • C based on calibrations closely suited to the latitudinal position of Guaymas Basin. To untangle the impact of these coupled drivers on the TEX 86 proxy, we further present a method to correct the overprints by both the water column and subsurface archaeal community using IPLs extracted from both of these sources. Although, we have not been able to test this model with lipid inputs from the overlying water column, we have demonstrated its effectiveness at removing sediment-sourced overprints, which may not be unique to hydrothermal systems. This approach should be capable of being extended to all near-surface marine sediment systems and may improve the quality of calibration models or climate reconstructions that are based on modern TEX 86 measures.
Data availability. Data files are available by request.
Author contributions. JNB and GTV collected and provided analysis of the acquired data and carried out the conceptualization of the study. JNB was responsible for all chemical preparation of the samples. JNB, GTV, CCW, SMS, and JSS wrote the paper. All authors were involved with the review of various drafts as the paper developed across the peer-review process. GTV, SMS, and JSS provided funding for this project.

Acknowledgements.
A special thank you is extended to associate editor Jack Middelburg and the Biogeosciences reviewers who provided highly constructive feedback. We are grateful to Carl Peters, formally at Saint Mary's University, who provided considerable feedback and advice during the course of this study. We further thank the officers, crew, and pilots of R/V Atlantis and HOV Alvin for their expert help at sea and their outstanding efforts acquiring the samples for this study. Julius Lipp, Florence Schubotz, and Kai-Uwe Hinrichs of MARUM (Center for Marine Environmental Sciences) assisted our lab in the development of lipidomic analytical techniques. Special thanks is extended to Clarissa Sit for the use of her HPLC-qToF-MS. Sean Sylvia assisted with the preparation of push cores used in sampling. Review statement. This paper was edited by Jack Middelburg and reviewed by three anonymous referees.