Hydrogen and carbon isotope fractionation factors of aerobic methane oxidation in deep-sea water

. Isotope fractionation factors associated with various biogeochemical processes are important in ensuring the reliable use of isotope tracers in biogeosciences at large. Methane is a key component of the subsurface biosphere and a notable greenhouse gas, making the accurate evaluation of methane cycles, including microbial methanotrophy,

Abstract. Isotope fractionation factors associated with various biogeochemical processes are important in ensuring the reliable use of isotope tracers in biogeosciences at large. Methane is a key component of the subsurface biosphere and a notable greenhouse gas, making the accurate evaluation of methane cycles, including microbial methanotrophy, imperative. Although the isotope fractionation factors associated with methanotrophy have been examined under various conditions, the dual-isotope fractionation factors of aerobic methanotrophy in oxic seawater remain unclear. Here, we investigated hydrogen and carbon isotope ratios of methane as well as the relevant biogeochemical parameters and microbial community compositions in hydrothermal plumes in the Okinawa Trough. Methanotrophs were found to be abundant in plumes above the Hatoma Knoll vent site, and we succeeded in simultaneously determining hydrogen and carbon isotope fractionation factors associated with the aerobic oxidation of methane (ε H = 49.4 ± 5.0 ‰, ε C = 5.2 ± 0.4 ‰) -the former being the first of its kind ever reported. This ε H value is comparable with values reported from terrestrial ecosystems but clearly lower than those from aerobic and anaerobic methanotroph enrichment cultures, as well as incubations of methanotrophic isolates. The covariation factor between δ 13 C CH 4 and δD CH 4 , (9.4 or 8.8 determined us-ing two different methods), was consistent with those from methanotrophic isolate incubations. These values are valuable for understanding dynamics of methane cycling in the marine realm, and future applications of the approach to other habitats with methanotrophic activity will help reveal whether the small ε H value observed is a ubiquitous feature across all marine systems.

Introduction
Stable isotope ratios have been widely used for tracing biogeochemical cycles and microbial activities (e.g., Ohkouchi et al., 2010;Sharp, 2017). To ensure the reliable use of isotope tracers in biogeoscience, it is important to understand the isotope fractionation factors associated with expected processes in the environment of interest. The isotope fractionation factor associated with microbial metabolism is known to be variable due to physiological responses to environmental conditions such as temperature and substrate availability (e.g., Valentine et al., 2004). As the complex biogeochemistry of a natural system cannot be practically rebuilt in full in either laboratory or models, determination of the fractionation factor in the natural environment is most appropriately done via observation.
Measuring multiple-isotope ratios of a compound is a powerful tool to trace a specific process. Currently available multi-isotope-ratio approaches include mass-independent fractionations (e.g., Farquhar et al., 2000;Michalski et al., 2003), position-specific differences (Blair et al., 1987;Yoshida and Toyoda, 2000), multiple substitution "clumped" isotopes Ghosh et al., 2006), and multi-element isotope ratios of a compound (e.g., Vogt et al., 2016). Among these, the two-dimensional analysis of dual-element isotope ratios in compounds, such as δ 15 N and δ 18 O of nitrite (Casciotti, 2016) and δ 13 C and δ 37 Cl of methyl chloride (Konno et al., 2013), is the most conventional method. As isotope effects associated with molecular dynamics, such as diffusion, occur to the same extent for each element in a compound, the two-dimensional analysis is useful in distinguishing the processes producing and consuming the compound.
Methane (CH 4 ) is a representative compound of the subsurface environment (e.g., Ijiri et al., 2018) and is also a notable atmospheric compound due to its strong greenhouse effect (e.g., Saunois et al., 2020). As such, the accurate evaluation of the marine CH 4 cycle is a matter of urgency (Reeburgh, 2007;Dean et al., 2018). Two-dimensional analysis with carbon and hydrogen isotope ratios has been applied to trace the origins and behaviors of CH 4 (Whiticar et al., 1986;Sugimoto and Wada, 1995;Nadalig et al., 2013). The kinetic isotope effect on the cleavage of the C-H bond in CH 4 causes isotope fractionation in the remnant CH 4 . Carbon and hydrogen isotope fractionation factors associated with the microbial consumption of CH 4 , methanotrophy, mediated by enzymes such as methane monooxygenase, have been determined through incubation of isolates (Feisthauer et al., 2011), incubation of enrichment cultures (Coleman et al., 1981;Kinnaman et al., 2007;Holler et al., 2009;Rasigraf et al., 2012), and observations in the natural environment (Snover and Quay, 2000;Kessler et al., 2006). While wide variations in isotope fractionation have been observed in each, with a factor between 2 ‰-36 ‰ for carbon and 36 ‰-320 ‰ for hydrogen, the ratio between carbon and hydrogen fractionation factors falls into a narrow range of 6-15 regardless of magnitudes of the fractionations (Feisthauer et al., 2011;Rasigraf et al., 2012, and references therein). The dual-isotope fractionation factors associated with methanotrophy have been examined under a variety of environmental and physiological properties such as aerobic-anaerobic, terrestrial and marine, temperatures, and forms of key enzyme methane monooxygenase (reviewed by Feisthauer et al., 2011;Rasigraf et al., 2012). Despite these efforts, the dual-isotope fractionation factors of methanotrophy in oxic seawater have never been determined. Only the carbon isotope fractionation factors of aerobic methanotrophy have been determined in oxic seawater, through observations of CH 4 -rich hydrothermal plumes (Tsunogai et al., 2000;Gharib, 2005;Gamo et al., 2010;Kawagucci et al., 2010).
Here, we studied both hydrogen and carbon isotope ratios of CH 4 , as well as the relevant biogeochemical processes and microbial community compositions, in plumes above two hydrothermal vents in the Okinawa Trough, at the Hatoma Knoll site and the Daisan-Kume Knoll ANA site, in order to simultaneously determine hydrogen and carbon isotope fractionation factors associated with aerobic oxidation of CH 4 in the marine environment. High-resolution vertical samplings inside and above seamount craters, at the base of which the vents are located, allowed us to capture gradual biogeochemical alteration of the plumes. Since vigorous venting of fluids with high CH 4 concentration (> 10 mM) with CH 4 -bearing boiling-derived bubbles at Hatoma Knoll (Toki et al., 2016) contrasts with relatively weak discharge of fluids with lower CH 4 concentration (< 1 mM) at the ANA site (Makabe et al., 2016), our initial intention was to reveal critical mechanisms controlling the fractionation factors by comparing them. We successfully determined the first hydrogen isotope fractionation factor of aerobic methanotrophy in seawater column at the Hatoma site, along with carbon isotope fractionation factors consistent with those reported to date (Tsunogai et al., 2000;Gharib, 2005;Gamo et al., 2010;Kawagucci et al., 2010).  (Toki et al., 2016) and the Daisan-Kume Knoll ANA site (depth: 1079 m) (Makabe et al., 2016;Uyeno et al., 2020). Hatoma Knoll is characterized by vigorous venting of fluids at 323 • C with boilingderived bubbles. The galatheoid squat lobster Shinkaia crosnieri, exhibiting ectosymbiosis with a microbial community of methanotrophs and thiotrophs (Watsuji et al., 2014), is dominant and densely distributed across the hydrothermally active area of Hatoma Knoll. Discovered in 2015 near Kumejima, the ANA site takes its name from an acronym of "Acoustic anomaly from Narrow pit Areas on caldera floor" (Nakamura et al., 2015), with a double meaning where ana means "hole" in Japanese. The maximum measured vent fluid temperature at the ANA site was 229 • C. The total vent fluid flux of the ANA site appeared to be more than an order of magnitude lower than that of Hatoma Knoll (unpublished data). Active venting at both sites is located on the bottom of the caldera at the top of the knolls, with the surrounding caldera walls being approximately 200 m high ( Fig. 1).
High-resolution vertical sampling with 23 seawater samples collected within 400 m altitude from the seafloor was .158 E -1079 m) using a conductivity-temperature-depth (CTD) profiler with a Carousel Multiple Sampling (CMS) system. The CTD-CMS system deployed during the MR17-03C cruise consisted of a CTD sensor (SBE911plus, Sea-bird Scientific), a SBE32 Carousel water sampler (Sea-bird Scientific) for 36 Niskin-X bottles, a dissolved oxygen (DO) sensor (RINKO), and a turbidity meter (Seapoint turbidity meter, Sea-Bird Scientific).

Analytical procedure
Concentrations and dual-isotope ratios of CH 4 and nitrous oxide (N 2 O) in the seawater samples were measured by a custom-made purge-and-trap system with a MAT253 isotope-ratio mass spectrometer (IRMS), following previously published methods (Hirota et al., 2010) and modifications (Okumura et al., 2016;Kawagucci et al., 2018). Immediately after the recovery of the CTD-CMS system on deck, each seawater sample in the Niskin bottle was subsampled into a 120 mL glass vial, sealed with a butyl rubber septum and aluminum cap after the addition of 0.2 mL mercurychloride-saturated solution, and stored at 4 • C until shorebased measurement in laboratories at the Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Yokosuka.
Each seawater sample was transferred by helium stream (rate: 100 mL min −1 ) into a 250 mL purge bottle to extract dissolved gases by helium bubbling and magnetic stirring. CH 4 and N 2 O in the stripped gas were purified by passing through a stainless-steel tubing coil trap held at −110 • C (ethanol / liquid-N 2 bath) as well as a chemical trap filled with magnesium perchlorate (Mg(ClO 4 ) 2 ; Merck KGaA) and Ascarite II (sodium-hydroxide-coated silica; Thomas Scientific), followed by a stainless-steel tubing trap filled with HayeSep-D porous polymer (60/80 mesh, Hayes Separations Inc.) held at −130 • C (ethanol / liquid-N 2 bath). Then, the CH 4 and N 2 O on the HayeSep-D trap were released into another helium stream (rate: 1.0 mL min −1 ) and again condensed on a capillary trap made of PoraPLOT Q (20 cm long, 0.32 mm i.d.) held at −196 • C (liquid-N 2 bath) for cryofocus, and finally released at room temperature.
After the complete separation of CH 4 and N 2 O from the other molecules using a GS Carbon-PLOT capillary column (30 m long, 0.32 mm i.d.) at 40 • C, the effluent CH 4 was put through a 960 • C combustion unit (Thermo Fisher Scientific) to convert into CO 2 prior to introduction into MAT253 via an open-split interface (GCC III, Thermo Fisher Scientific). N 2 O was introduced into GCC III and IRMS without conversion. Carbon, nitrogen, and oxygen isotope ratios were obtained through simultaneous monitoring of CO + 2 and N 2 O + isotopologues at m/z = 44, 45, and 46 by Faraday cups with 10-times-higher resistors, for increased sensitivity compared to commercially available general settings . For the determination of hydrogen isotope ratio of CH 4 , another vial of the same sample was used, with the analytical procedure being almost identical to the above except a 1440 • C pyrolysis unit (Thermo Fisher Scientific, Massachusetts, USA) was applied for conversion into H 2 , m/z = 2 and 3 (Okumura et al., 2016). Isotope ratios are represented by conventional δ notation and presented at a per mill scale. Errors for the analyses of deep-sea water conducted during the present study were estimated from repeated analyses of a sample and were within 10 % for N 2 O concentration, 0.2 ‰ for δ 15 N N 2 O , 0.5 ‰ for δ 18 O N 2 O , 20 % for CH 4 concentration, 0.3 ‰ for δ 13 C CH 4 , and 5 ‰ for δD CH 4 . Table S1 in the Supplement.

Results of the N 2 O analysis are provided in
Isotope fractionation factor caused by kinetics, ε (‰), can be expressed by the following equation: where k is reaction constant while ( h k/ l k) is the kinetic isotope effect and superscripts "h" and "l" represent the heavier isotopologue (e.g., 13 C) and the lighter one (e.g., 12 C), respectively. Observations of fractional yield and isotope composition of the reactant in a closed system can be used to determine the fractionation factor using the following equation (Mariotti et al., 1981): where the subscripts designate initial reactant (r 0 ) and reactant remaining in a sample (r x ), and f represents the fraction of the reactant consumed after the reaction. When δr x is plotted as a function of ln(1 − f ), the ε value is given by the slope of the line. A factor expressing covariation between δ 13 C CH 4 and δD CH 4 , , is defined (Elsner, 2010;Feisthauer et al., 2011) by the following equation: Assuming variations in the isotope ratios of CH 4 in the environment are attributable solely to CH 4 oxidation, can also be calculated using another equation without the parameter f from discriminations of isotope ratios observed (Tsunogai et al., 2020), as follows: The manganese concentration was measured in a laboratory of Kaiyo-Keisoku Co. Ltd. (Kei Okamura) using 100 mL of unfiltered seawater via the luminol-H 2 O 2 chemiluminescence detection method (Ishibashi et al., 1997; calibrated with the international standard NASS-5. Ammonium concentrations were determined on board the research vessel with the conventional autoanalyzer method. The concentration of particulate ATP (pATP) was determined by a luciferin-luciferase assay on board the ship. Seawater samples collected using the CTD-CMS were transferred aseptically to clean plastic tubes, and an aliquot of this subsample was immediately filtrated through a sterilized membrane filter unit (0.2 µm pore size). Concentrations of ATP in both unfiltered and filtered aliquots were measured with a simplified quantification assay using the ATP assay kit (CheckLight HS, Kikkoman Biochemifa) without pre-concentration and extraction process. To start, 100 µL of cell lysis reagent was directly added to 100 µL of seawater sample for ATP extraction in a disposable test tube. After leaving the lysate for at least 30 min at room temperature to minimize the effect of temperature differences between the samples, 100 µL of the luciferin-luciferase mixture reagent was added to the lysate, and the luminescence intensity was immediately measured using a desktop photodetector (Luminescencer Octa AB-2270, ATTO). The ATP concentrations from unfiltered and filtered subsamples were regarded as total ATP and dissolved ATP, respectively. The pATP concentration was calculated by subtracting the dissolved ATP concentration from the total ATP concentration.
Samples for measuring the microbial cell density were fixed with 0.5 % (wt / v) glutaraldehyde (final concentration) in 2 mL cryo-vials on board and stored at −80 • C until further analysis. For cell density measurements, 200 µL of each sample was stained by SYBR Green I nucleic acid gel stain (Thermo Fisher Scientific, ×5 of manufacture's stock) at room temperature for > 10 min. The total microbial cell abundance in 100 µL of sample was determined using an Attune NxT acoustic focusing flow cytometer (Thermo Fisher Scientific) by their signature in a plot of green fluorescence versus side scatter (Brussaard, 2004;Giorgio et al., 1996).
To obtain microbial community structure based on amplicon sequencing of the 16 rRNA gene, 2-4 L seawater samples were filtered with 0.2 µm pore-size cellulose nitrate or acetate membrane filters. The filters were stored at −80 • C until environmental DNA extraction, following previously published methods (Hirai et al., 2017). The 16S SSU rRNA gene was amplified from the extracted DNA with the primer mixture of 530F and 907R, using the LA Taq polymerase with GC buffer (Takara Bio) as described previously (Nunoura et al., 2012;Hiraoka et al., 2020). The amplicon sequencing libraries were sequenced using an Illumina MiSeq high-throughput sequencer (2 × 300 paired-end platform) at JAMSTEC. The sequence data generated herein are publicly available in the DDBJ sequence read archive (DRA) under the BioProject PRJDB11835.

Hydrothermal plume signature
Vertical profiles of biogeochemical parameters drawn by altitude from the seafloor are presented in Fig. 2, and with all data provided as Table S1. Temperatures of seawater collected ranged between 3.5 and 7.0 • C (Fig. 2a). Water columns showed simultaneous increases in turbidity, manganese, and CH 4 and demonstrated the presence of hydrothermal plumes above both the Hatoma Knoll and ANA sites. Temperature changes (from the bottom water temperature) fluctuated at both sites along the vertical profile, with no clearly identifiable peaks. Dissolved oxygen profiles at both sites demonstrated sufficiently oxic seawater for aerobic metabolisms through the water column observed. The turbidity profile above Hatoma Knoll exhibited two vertically broad turbid water masses, one near the seafloor and another at 100-200 m altitude. At the ANA station the turbidity profile displayed a single peak centered around 130 m altitude. The maximum turbidity seen above Hatoma Knoll was an order of magnitude lower (0.6 FTU) than at the ANA site. Drastic change in temperature and decline of turbidity above 200 m altitude in the Hatoma Knoll plume suggests the height of the caldera rim. Vertical patterns of manganese concentrations at each station were generally similar to turbidity. The maximum manganese concentrations observed were lower at Hatoma Knoll (21 nM) than the ANA site (119 nM).
Vertical distribution of CH 4 concentrations also showed similar patterns to those of turbidity and manganese, but they differed from them in that the maximum CH 4 concentrations observed were clearly higher at Hatoma Knoll (935 nM) compared to ANA (254 nM). The CH 4 /Mn ratios varied between 1 and 45 above Hatoma Knoll but were roughly constant around 1.5 in the water column above ANA. Ammonium was more enriched in the water column above Hatoma Knoll compared to the ANA site. The vertical pattern of ammonium concentrations above Hatoma was similar to that of CH 4 concentrations.

Methane isotope composition
Carbon and hydrogen isotope ratios of CH 4 were distinct between the two vent sites (Figs. 2 and 3). Above Hatoma Knoll, vertical patterns of CH 4 concentrations and the δ 13 C CH 4 values were mirror images of each other, with the δ 13 C CH 4 values being the lowest (approximately −47 ‰) at 107 m altitude where the CH 4 concentrations were high, while the highest δ 13 C CH 4 value (−35.6 ‰) appeared at 77 m altitude where there was a minimum in CH 4 concentrations. The vertical pattern of δD CH 4 values above Hatoma was similar to that of the δ 13 C CH 4 value, peaking at 77 m altitude with +17 ‰. Above the ANA site, however, the δ 13 C CH 4 and δD CH 4 values were almost constant across the depth gradient at −28 ‰ and −110 ‰, respectively. The 13 C-D diagram for CH 4 observed above Hatoma Knoll demonstrated a linear trend (Fig. 3). Isotopic composition of the low-δ 13 C CH 4 root of the linear distribution is consistent with the δ 13 C CH 4 values of Hatoma vent fluids determined by direct sampling (−54 ‰ to −49 ‰) (Toki et al., 2016), and the δD CH 4 values determined from all known deep-sea vent fluids (approximately −120 ‰) (Proskurowski et al., 2006). A least-square linear fitting for the dataset obtained above Hatoma Knoll yielded a slope of 8.8 (Fig. 3), which corresponds to according to Eq. (4).

Microbiological characteristics
Total cell density (cells mL −1 ) above Hatoma Knoll generally fell within a narrow range between 1.9-7.6 × 10 4 across the depth gradient (Fig. 2j). The cytogram obtained from cell counting using flow cytometry exhibited two clusters, each of which represented the typical microbial cells and cells with typical fluorescence signals with significantly higher side-scatter signals, named the high-side-scatter population (HSS). The HSS is attributable to large cells and/or cell aggregates. Only two samples between 200-230 m showed sig- nificant HSS cell abundances at Hatoma Knoll. At the ANA site, the total cell density below 250 m altitude was an order of magnitude higher (≤ 4.9 × 10 5 cells mL −1 ) than that above 250 m (4.0 × 10 4 cells mL −1 ). Between 40-160 m altitude, HSS cells occupied 20 %-58 % of the total cells. The cell densities in these samples are possibly underestimated because a HSS signal possibly consists of aggregated cells.
pATP concentrations above Hatoma Knoll increased with depth and were nearly constant below 200 m (the height of the caldera rim) at an order of 10 1 pmol L −1 (Fig. 2k). Above the ANA site, pATP concentrations were the highest between 33-140 m altitude at 3.0 × 10 2 pmol L −1 . pATP concentrations above 250 m altitude were comparable between the two vent sites on the order of 10 0 . Cellular ATP contents (ng of ATP per cell) below 10 −7 were observed above 300 m altitude at both the Hatoma Knoll and ANA sites (Fig. 4). The level observed is consistent with values previously reported from open ocean water (e.g., Winn and Karl, 1986), suggesting that they represent background levels without any effect from hydrothermal input. Gradual increase in the cellular ATP content along with increased depth was found from the caldera rim depth (200 m altitude) to the seafloor at Hatoma as well as across the depth gradient above the ANA site. Cellular ATP contents above 10 −7 ng of ATP per cell, despite the low cell density (< 10 5 cells mL −1 ), were observed below 200 m altitude at Hatoma Knoll.
Microbial community analysis using amplicon sequences of the 16S rRNA gene revealed distinct communities between the Hatoma Knoll and ANA sites (Fig. 5). In the water column from the seafloor to 200 m altitude above Hatoma, more than half of the community consisted of members of SUP05 (composed of three ASVs: SUP05_1, _2, and _3), a well-known group of chemolithotrophic sulfuroxidizing Gammaproteobacteria frequently detected around hydrothermal vents and within hydrothermal plumes (Dick et al., 2013). Approximately 10 % of the community was occupied by the aerobic methanotrophic Gammaproteobac- teria family Methylococcaceae (composed of two ASVs: Methylococcaceae_1 and _2) and > 5 % by the ammoniumoxidizing archaeal family Nitrosopumilaceae (composed of two ASVs: Nitrosopumilaceae_1 and _2). Occupancies of SUP05 decreased with increasing altitude from the 200 m mark. The microbial community in the water column above ANA consisted of over 50 % SUP05 (composed of three ASVs: SUP05_1 and _2) and 19 % Sulfurovaceae (sulfuroxidizing Epsilonproteobacteria family) from the seafloor to 250 m altitude. In contrast to Hatoma Knoll, methanotrophic lineages were not detected at the ANA site.

Active aerobic methanotrophy in Hatoma Knoll plume
Multiple lines of evidence point to the occurrence of microbial oxidation of CH 4 in the hydrothermal plume above Hatoma Knoll, between the water column at the seafloor to 200 m altitude, but not at the ANA site. Manganese concentration has been utilized as an indicator for the dilution of Mn-rich hydrothermal vent fluid with Mn-depleted ambient seawater (e.g., Kawagucci et al., 2008), owing to the long life of manganese compared to particles and aerobically energetic molecules like CH 4 (Kadko et al., 1990). In principle, a two-component mixing is represented by a straight line on the Mn plot, and a downward deviation from the ideal mixing line suggests significant removal of the counterpart compound (German and Seyfried, 2014). Indeed, the Mn plot of waters above Hatoma Knoll exhibited downward convex curves and not a straight line for CH 4 concentrations, ammonium concentrations, and turbidity (Fig. 6). The NH 4 -Mn dataset of the Hatoma plume observed is shifted downwards from the ideal mixing line, assuming NH 4 /Mn ratios of 11-18 observed in the Hatoma vent fluid (Toki et al., 2016), confirming the removal of ammonium from the plume. Despite the lack of available CH 4 /Mn data for the estimated endmember fluid at the Hatoma Knoll vent site ASVs, in terms of total read abundance in this study, were defined as the major group. Read abundance of the major group corresponds to 83 % of the total read abundance of the samples. The names and colors in the legend are identical in these two panels. (Toki et al., 2016), the downward convex curve of CH 4 concentrations on the Mn plot strongly suggests CH 4 removal. From the viewpoint of microbial ecology, the ATP-rich microbial cells (Fig. 4) and the significant appearance of aerobic methanotrophic lineage (Fig. 5) observed in the plume above Hatoma Knoll are both strongly indicative of active microbial methanotrophy within the plume. In contrast, the rather uniform CH 4 /Mn ratios above the ANA site regardless of concentrations suggest negligible CH 4 consumption in the plume. No dominant methanotrophic microbial groups were detected by the 16S rRNA gene community analysis above ANA, supporting a lack of methanotrophic activity. The contrasting microbial community composition and methanotrophic activity seen between plumes above the Hatoma Knoll and ANA sites are likely attributable to differences in the vent fluid chemistry. High-temperature hydrothermal fluids directly collected from vent orifices at Hatoma Knoll showed relatively high CH 4 concentrations sometimes > 10 mM (Toki et al., 2016), in line with significant methanotrophic activity being detected in the Hatoma plume. In the plume originating from the CH 4 -rich Guaymas Basin site (McDermott et al., 2015), the presence of Methylococcaceae and the high expression of the methane oxidation gene (pmo) were revealed (Lesniewski et al., 2012). In contrast, high-temperature vent fluid collected from the ANA site only contained CH 4 at concentrations below 1 mM (Makabe et al., 2016), which may explain the lack of significant methanotrophic activity in the plume above this site. The increase in total cell abundance (Fig. 2) and the dominance of the plume-associated sulfur-oxidizing bacteria group SUP05 (Fig. 5) (Sunamura et al., 2004) in the ANA plume are evidence for significant shifts of the entire microbial community supported by hydrothermal fluids. Although available CH 4 likely fuels the methanotrophs, the low concentration only allows for slow methanotrophic activity compared to the residence time of the plume. This could make it difficult to de-  tect signatures of methane consumption using concentrations and isotope ratios.

Determination of isotope fractionation factors at Hatoma Knoll
Hereafter, we assume the CH 4 /Mn ratios of the plume above Hatoma vary only by aerobic CH 4 oxidation, for the evaluation of isotope fractionation factors. When the δ 13 C CH 4 and δD CH 4 values of the Hatoma Knoll plume are plotted as a function of the CH 4 /Mn ratio, there are observable increments of δ values along decrements of the CH 4 /Mn ratio (Fig. 7). This phenomenon is well explained by the kinetic isotope effect on aerobic CH 4 oxidation which causes the enrichment of heavier isotopologues, 13 CH 4 and CH 3 D, in the remnant reactant. Some large scatters of the δ 13 C CH 4 and δD CH 4 values are seen, particularly at low CH 4 /Mn ranges (e.g., < 10) -these data include two distinct seawater masses below 200 m altitude (the height of the caldera wall), corresponding to the high pATP water inside the caldera (Fig. 2) and waters coming in from above 200 m (Fig. 1). As our aim is to determine the isotope fractionation factors from the observed values, the gradual changes in isotope composition within the water mass below 200 m altitude were used for further analysis.
According to previous studies (e.g., Gamo et al., 2010), we applied the CH 4 /Mn ratio instead of f for Eq. (2), as follows: where a represents carbon (C) or hydrogen (H) isotopes. The δ − ln[CH 4 /Mn] plot analysis for the plume above Hatoma Knoll yielded slopes representing ε C and ε H values of 5.2 ± 0.4 ‰ and 49.4±5.0 ‰, respectively (Fig. 7). The ε C and ε H values in turn yielded a value of 9.4 according to Eq. (3). This ε-based value calculated from the plume sample is similar to the δ-based value calculated using the entire dataset from Hatoma Knoll (8.8) (Fig. 3).
The ε C values obtained through observations of deep-sea hydrothermal plumes reported to date are comparable among the sites and regions, including the 27.5-32.5 • S area on the East Pacific Rise (4 ‰-6 ‰) (Gharib, 2005), the Myojin Knoll field on the Izu-Ogasawara Arc (5 ± 1) (Tsunogai et al., 2000), the Daiyon-Yonaguni site in the Okinawa Trough (5 ‰ and 12 ‰) , Izena Hole also in the Okinawa Trough (< 7 ‰) (Kawagucci et al., 2010), and Hatoma Knoll (5.2 ± 0.4 ‰) (this study). These ε C values from the deep-sea plumes are lower than those estimated from methanotrophic isolate incubations (18.8 ‰-27.9 ‰) and methanotrophic communities (7.9 ‰-26.6 ‰) when not considering some exceptional values (e.g., Feisthauer et al., 2011). Regardless of the ε H and ε C values, the value of the hydrothermal plume (9.4 and 8.8) reported herein is consistent with those from methanotrophic isolate incubations (7.3-10.5), including the cultivation of M. capsulatus (Feisthauer et al., 2011). The consistencies in both ε C values among deep-sea plumes and values among isolates and plumes suggest that the ε H , ε C , and values determined in this study are appropriate. However, the reasons why the ε H and ε C values are smaller in the natural environment including deep-sea water columns compared to others remain unclear at this point. Although temperature dependency of isotope fractionation factors on the methanotrophy was discussed (Chanton et al., 2008), the small ε values at 4 • C obtained in this study oppose the trend reported so far. A possible explanation, originally proposed for the case of benzene biodegradation (Fischer et al., 2009), is that the cells take up only a limited amount of methane, which is then virtually all consumed, leading to little change in the isotopic ratios of water column methane. Future applications of the approach to the other marine habitats with methanotrophic activity will reveal whether or not the small ε H value observed here is ubiquitous in the marine realm.

Concluding remarks
The ε H , ε C , and values associated with aerobic CH 4 oxidation in seawater above hydrothermally active areas are useful for understanding marine CH 4 dynamics. The δD CH 4 values of seafloor hydrothermal vent fluids and hydrocarbon seep fluids are expected to be −130 ‰ and −180 ‰, respectively (Whiticar, 1986;Okumura et al., 2016). As such, observations of the δD CH 4 values in the water column above the geofluid sites, combined with the ε H value of 49.4 ± 5.0, enable us to estimate how CH 4 oxidation has progressed through Eq. (2). The estimation of the fraction as well as the δ 13 C CH 4 values of vent plumes further allow us to estimate the δ 13 C CH 4 value of the endmember effluent CH 4 . The estimated δ 13 C CH 4 value of the endmember geofluid denotes the origin of CH 4 there, allowing "sneak peeks" of the ongoing subseafloor process.
Our approach to determine isotope fractionation factors in the seawater column environment by high-resolution hydrothermal plume sampling can be applied not only to conventional carbon and hydrogen isotope ratios of CH 4 but also to "clumped" isotope composition of CH 4 (Wang et al., 2016;Ono et al., 2021). The same approach at Hatoma Knoll is also applicable to determining isotope fractionation factors for ammonium oxidation in marine environment because of the decreases in NH 4 /Mn ratios (Fig. 6) and significant appearance of the ammonium-oxidizing archaeal family Nitrosopumilaceae (Fig. 5). If the focus was on sulfur isotope fractionation of aerobic sulfide oxidation, then the same approach at the ANA site would be appropriate because sulfur-oxidizing microbes dominated the microbial community (Fig. 5). Drastic changes in N 2 O and H 2 concentrations in plumes at the Izena Cauldron (Kawagucci et al., 2010) also allow us to determine the isotope fractionation factors associated with their metabolisms in the water column by the same approach. These are foci in future studies.
Data availability. The dataset reported is available in Table S1 in the Supplement.
Author contributions. SK and TY designed the study. YM, AM, TF, YO, TN, and TY conducted chemical and microbiological analyses. SK made a draft. All authors contributed to sampling and gave final approval for submission and publication.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.