Sub-oxycline methane oxidation can fully uptake CH4 produced in sediments: case study of a lake in Siberia.

It is commonly assumed that methane (CH4) released by lakes into the atmosphere is mainly produced in anoxic sediment and transported by diffusion or ebullition through the water column to the surface of the lake. In contrast to that prevailing idea, it has been gradually established that the epilimnetic CH4 does not originate exclusively from sediments but is also locally produced or laterally transported from the littoral zone. Therefore, CH4 cycling in the epilimnion and the hypolimnion might not be as closely linked as previously thought. We utilized a high-resolution method used to determine dissolved CH4 concentration to analyze a Siberian lake in which epilimnetic and hypolimnetic CH4 cycles were fully segregated by a section of the water column where CH4 was not detected. This layer, with no detected CH4, was well below the oxycline and the photic zone and thus assumed to be anaerobic. However, on the basis of a diffusion-reaction model, molecular biology, and stable isotope analyses, we determined that this layer takes up all the CH4 produced in the sediments and the deepest section of the hypolimnion. We concluded that there was no CH4 exchange between the hypolimnion (dominated by methanotrophy and methanogenesis) and the epilimnion (dominated by methane lateral transport and/or oxic production), resulting in a vertically segregated lake internal CH4 cycle.

Introduction

Methane (CH4) released from lakes to the atmosphere is generally assumed to be produced in anoxic sediment and transported by diffusion or ebullition through the water column[1,2], where it is subject to oxidation, generally identified as a major sink[2-5]. In stratified lakes, the diffusion barrier formed by the thermocline promotes CH4 storage in the hypolimnion. The stored CH4 is released mainly to the epilimnion and then to the atmosphere during water column overturn[6,7]. Lateral transport of CH4 from the littoral to the pelagic zones may also substantially modify the CH4 balance at the epilimnion of the lakes[1,8,9]. Although this has been the prevailing theory, several other CH4 cycle processes have been discovered. For instance, CH4 production in the epilimnion under aerobic conditions has been observed in several aquatic ecosystems[8,10-13]. Another process previously reported is CH4 oxidation below the oxycline, which is typically associated with microaerophilic conditions[4], oxygenic photosynthesis[4,14], or, in some cases, assumed to be anaerobic[15]. These processes explain why deviations to standard CH4 concentration profiles (usually decreasing from the bottom to the surface) are often observed as concentrations increase or decrease in local areas.

In this study, we observed an atypical CH4 profile in a strongly stratified lake (Sila Lake) in north-central Siberia. From the bottom to the surface, the profile showed a sharp decrease of the CH4 concentration below the detection limit of our method (5 nmol L−1) in the benthic zone of the water column, well below the oxycline and the photic zone, and showed an increase in CH4 concentration in the aerobic epilimnion. This profile was created using a high resolution, high sensitivity method[16,17] which allowed the net methane production and oxidation rates to be determined using a diffusion-reaction model, showing a short distance transition between both processes in the hypolimnion. This CH4 concentration profile was correlated to methanotroph and methanogenic archaea abundances determined by pmoA and mcrA gene qPCR, respectively.

Results and Discussion

The depth measurements collected were used to develop a bathymetric map (Fig. S1) that shows a maximum depth of 12 m in the southeastern branch of the lake. Satellite imagery (Google Earth Engine) was used to determine that the area of the lake was 3.6 Ha. Combining the area of the lake and the bathymetric map, we constructed a hypsometric histogram to estimate the total water volume of the lake to be 168,000 m3. The lake was strongly stratified, with a mean surface temperature of 16.5 ± 0.7 °C (mean ± one standard deviation), a bottom temperature of 3.4 ± 0.3 °C, a thermocline 2–4 m deep, and a maximum gradient of 8.6 °C m−1 (Fig. S2). The oxycline matched the thermocline and anoxic conditions (i.e., dissolved oxygen (DO) concentration below 0.3 µmol L−1) were found at depths below 4 m. The Secchi depth was 1.8 ± 0.15 m, meaning that the euphotic depth (Z1%), i.e., depth at which the photosynthetically available radiation was 1% of its surface value, was estimated at 4.9 ± 0.4 to 5.6 ± 0.4 m.

In addition to each sampling and monitoring station, the profiles of dissolved CH4 and carbon dioxide (CO2) concentration (CCH4 and CCO2, respectively) were measured in triplicate at the center of the lake (P1, see materials and methods, Fig. 1A). Epilimnetic CCH4 and CCO2, at depths between 0 and 2.5 m, were 0.56 ± 0.05 μmol L−1 and 44 ± 7 μmol L−1, respectively. Below 2.5 m, CCH4 rapidly decreased and was no longer detected below depths of 2.96 ± 0.13 m (until depths greater than 7 m). On the contrary, an increase of CCO2 was observed to 209 ± 16 μmol L−1. At a depth of 7.12 ± 0.14 m, well below the oxycline and Z1%, CCH4 rapidly increased and reached 196 ± 18 μmol L−1 at 10 m. Higher values were observed below 10 m but these measurements were discarded as potentially being the result of sediment disturbances, at a depth of 10.4–10.7 m. Inversely, CCO2 decreased in the deepest section of the hypolimnion to 135 ± 11 μmol L−1 at 10 m. The section of the water column where CH4 was not detected (3–7 m), called the “methane minimum zone” (MMZ), expanded along the transversal and longitudinal transects (Fig. 2), and acted as a segregation zone between the hypolimnetic and the epilimnetic CH4.

Figure 1

(A) Depth profiles of dissolved CH4 (CCH4) and CO2 (CCO2) in logarithmic scale; (B) net methane production rates of one profile (for clarity purposes, replicate profiles are shown in Fig. S3); and relative abundance, in reference to total prokaryotes, of pmoA genes (green bars), mcrA genes (blue bars) and total archaea genes (pink bars). Epilimnion NMPR are multiplied by 1000; (C) δ13C-DIC (brown square), δ13C-CH4 (green circles), and δ2H-CH4 (blue triangles).

Figure 2

Longitudinal (east-west; A) and transversal (north-south; B) transectional maps of dissolved CH4 concentration showing the expansion of the minimum methane zone (MMZ). ND stands for not detected.

To better characterize the CH4 profiles observed, a diffusion-reaction model was applied to estimate net methane production rates (NMPR; Fig. 1B, S3D–F) and vertical fluxes through the water column (Fig. S3A–C). As shown in the diagram, a peak of negative NMPR, i.e., CH4 sink, was observed below the MMZ, immediately above a peak of positive NMPR; i.e., CH4 source. The same pattern was observed in the three replicates (Fig. S3), although with differences in depths (maximum difference of 0.54 m between the three negative NMPR peaks), which were likely caused by boat motion and inaccurate positioning, evidenced by differences in water column depth varying between 10.4 and 10.7 m. It is noteworthy that the depth intervals between CH4 sink and source peaks were consistently small, 0.51 ± 0.11 m, indicating a sudden shift in the dominant CH4 process. The average minimum and maximum NMPR corresponding to the first methanotrophic and the methanogenic peaks, in a downward direction, were –1.82 ± 0.34 and 1.69 ± 0.68 μmol L−1 h−1, respectively, which suggests that the magnitude of both processes was comparable.

The isotopic signature of CH4 and the CH4 production fractionation factor (α; see materials and methods), show two homogeneous zones (Fig. 1C). In the hypolimnetic layer, depleted δ13C-CH4 (–78.6 to –82.2‰) and δ2H-CH4 (–360 to –383‰), coupled to relatively high α values (1.072 to 1.077), strongly support hydrogenotrophic production of CH4 from CO2[12,18]. On the contrary, at the epilimnion, enriched values of δ13C-CH4 (–49.8 to –50.8‰) and δ2H-CH4 (–274.2 to –287.8‰), coupled with relatively low α values (1.046 to 1.048), are consistent with a predominance of acetoclastic production of CH4[12,18]. Indeed, according to Whiticar and Faber[19], α from 1.055 to 1.09 most likely corresponds to hydrogenotrophic methanogenesis, while α from 1.04 to 1.055 corresponds to acetoclastic methanogenesis.

The abundance of methane oxidizing bacteria (MOB) quantified through the relative abundance of pmoA gene (Fig. 1B) was minimal in the MMZ, representing <0.03% of the prokaryotic community, as expected from the absence of detectable levels of CH4. Higher abundance of MOB was detected within the anoxic water column (8–10 m), where MOB represented from 1.1% to 2.5% of the total community. Notably, the peak of MOB abundance was found at 8 m depth (1.6 ± 0.05 × 108 pmoA copies L−1), which coincided with the depth of maximum methanotrophic activity (Fig. 1B, S3), suggesting that MOB were potentially major contributors to the CH4 oxidation in this lake. This is not surprising since the possible involvement of aerobic methanotrophs in anaerobic methane oxidation has been previously suggested[4,20,21]. Although MOB activity was not determined, this finding indicates that MOB might be active in anoxic waters. Quantitative PCR also revealed that the amount of mcrA gene was minimal in the MMZ and remained low where the methanotrophic peak was observed, being two orders of magnitude less abundant than the pmoA gene (Fig. 1B). This result suggests that mcrA-carrying anaerobic methane oxidizing archaea (ANME) were not major contributors to the strong methane oxidation observed in this anoxic layer. The mcrA gene abundance only increased significantly at the bottom of the water column (1.1 × 107 ± 0.3 × 107 pmoA copies L−1 at 10 m depth); similar results were observed with total archaea (quantified through their 16 S rRNA gene). In addition, a local mcrA maximum (2.3 × 105 ± 0.4 × 105 copies L−1) was observed in the oxic epilimnion at a depth of 2 m, which indicates the presence of methanogens despite the oxic environment.

The major CH4 oxidation observed in the anoxic epilimnion raises the important question of the oxic or anoxic nature of the process. Several arguments suggest the absence of oxygen, either diffused from the top or locally produced. First, the euphotic depth of 4.9–5.6 m was 2.9–3.6 m above the peak of methanotrophic activity, where the photosynthetic available radiation was 0.03–0.08% of its surface value. Second, the peak of methanotrophic activity of 1.82 ± 0.34 μmol CH4 L−1 h−1 would require 3.64 ± 0.68 μmol O2 L−1 h−1, i.e., two moles of O2 required to oxidize one mole of CH4. That oxygen requirement is within the higher range of primary production reported for the euphotic or epilimnion zones of 118 lakes[22] worldwide, thus unlikely to occur 2.9–3.6 m below the euphotic depth. Third, the average methanotrophic activity over the entire water column, determined by integration of negative NMPR, corresponded to 611 ± 80 μmol CH4 m−2 h−1. That figure would require a counter flux of 1,222 ± 160 μmol O2 m−2 h−1. According to a maximum hypolimnetic diffusivity of 2.819 m2 h−1[23], the DO gradient along the water column, required to sustain aerobic methanotrophy, corresponds to 433 μmol L−1 per meter of water column depth, which is incompatible with the measured DO concentration in the hypolimnion; i.e., below 0.3 μmol L−1. Fourth, the proximity of methanogenic peaks, observed 0.51 ± 0.11 m below the methanotrophic peak, does not support a sudden shift from an active oxic to a strict anaerobic process. The unlikeliness of aerobic CH4 oxidation raises the question of the electron acceptor as an alternative to O2. No evidence arises from nitrate or nitrite (Table S2), as none of these electron acceptors for anaerobic methane oxidation[20,24] showed a significant concentration change in the region where methanotrophic activity was found, and both nitrate or nitrite in the lake water column were below or at the lower range of concentration to support AOM coupled to denitrification[25]. Thus, the present work does not reveal the possible electron acceptors for AOM in the epilimnion of the lake, since sulfate, oxidized metals, and organic matter[26] were not tested.

Our results show that the MMZ segregates two different CH4 zones of the water column with no diffusive exchange between them. The existence of an MMZ segregation zone is supported not only by the CH4 profiles, discarding diffusive transfer between the epilimnion and the hypolimnion, but also by the CH4 isotopic signature, suggesting different CH4 origins in each zone. However, despite segregation between the epilimnion and the hypolimnion, bubbles formed in the sediments might transfer CH4 to the epilimnion during their migration to the surface of the lake. Nevertheless, at each sampling and monitoring station, CH4 and CO2 fluxes were determined in triplicate with a dynamic closed chamber coupled to a greenhouse gas analyzer (UGGA 30 P, Los Gatos Research, CA, USA; data acquisition frequency of 1 s−1 and CH4 sensitivity of 30 ppb)[17]. During a total of 54 flux measurements over a total time of 4.5 h, no evidence of bubbling was found; i.e., peak increase of CH4 concentration within the chamber[17]. Bubbles were only occasionally visually observed in the littoral region at the western section of the lake. Thus, ebullitive transfer of CH4 from the hypolimnion to the epilimnion, is unlikely although not discountable. The CH4 found at the epilimnion might be considered as the product of local oxic production and/or lateral transport from the littoral zone[1,8,10-13].

We estimated the flux, downward to the MMZ, as 0.62 ± 0.15 μmol CH4 m−2 h−1 (Fig. S3 A–C). A triplicate measurement of the CH4 flux to the atmosphere, at each of the three locations where the CCH4 profiles were made, gave an average of 54 ± 14 μmol CH4 m−2 h−1. Thus, the transfer of CH4 to the MMZ was negligible and, assuming steady-state concentrations (i.e., those that did not change over time), we estimated the oxic production and/or lateral transport to equalize flux to the atmosphere, i.e., 54 ± 14 μmol CH4 m−2 h−1. Figure 3 shows the CH4 mass balance of Sila Lake at the time of characterization and indicates that about 92% of the total CH4 produced in or transported to the lake is oxidized. The steady-state assumption used to establish the CH4 mass balance might be a simplistic consideration regarding the epilimnion. Sila Lake is a northern lake with an ice-free period ranging from three to four months (personal communication with local inhabitants); this short period of summer stratification suggests relatively rapid changes of the water column structure, which is potentially contradictory to steady-state conditions. However, the total amount of CH4 present in the epilimnion, considering a depth layer of 2.3 m, was estimated to be 1.34 ± 0.12 mmol CH4 m−2. Dividing the amount of CH4 present in the epilimnion by the flux to the atmosphere mentioned above, we estimated the CH4 turnover time to range from 22–27 h. The latter indicates that the CH4 cycle in the epilimnion is dynamic and that the CH4 emitted to the atmosphere is rapidly replaced by lateral transport and/or oxic production.

Figure 3

CH4 mass balance in the water column of Sila Lake. Arrows indicate CH4 transport and production while numbers indicate the magnitude of these processes, all of which are expressed per unit of lake area (μmol CH4 m−2 h−1). The minimum methane zone is called MMZ.

This study shows that, in at least some cases, CH4 cycling in stratified lakes should be considered as two segregated systems with no exchange between them. While the hypolimnion is dominated by methanotrophy and hydrogenotrophic methanogenesis, the epilimnion is dominated by CH4 originating from acetoclastic production, in situ or transferred from surrounding terrestrial ecosystems. The segregation between the CH4 cycle in the epilimnion and the hypolimnion has been suggested by several reports[5,8,10,12,13] and used as a supporting evidence that epilimnetic CH4 is locally produced[8,10-13,27] or transported from the littoral zone[1,8,28]. In addition to similar conclusions, the present work was based on an innovative high resolution and sensitive method for determining CCH4, which allowed the NMPR profile to be created and demonstrated that methanotrophy can occur well below the oxycline and euphotic zone of a lake. This method also revealed the drastic shift, relative to depth, between the two opposing dominant CH4 processes of anaerobic methanogenesis and methanotrophy, with the latter being potentially anoxic or at least not associated with locally produced oxygen. Notably, pmoA-carrying bacteria were the major microbial contributors of methane oxidation in this lake and largely predominated (over potential ANMEs) at the depth of highest methane oxidation activity, leading to a complete exhaustion of produced methane at the hypolimnetic water column. The extent to which the behavior observed in Sila Lake can be extrapolated to other lakes is unknown but is certainly an important milestone and should be investigated further, along with the identity of the methane oxidizing bacteria potentially involved in AOM in this lake and others.

Materials and Methods

A small, Y-shaped unnamed lake, located in a discontinuous permafrost area of northern Siberia[29] 5 km north of the town of Igarka, Krasnoyarsk Kray, Russia (Fig. S1), was selected and called Sila for the purpose of this study (Lat. 67.5138, Long. 86.5915). Sila Lake is of glacial origin and heavily influenced by thermokarst. It is surrounded by northern forest and peatlands, which is the dominant landscape of the region. About one-third of the perimeter (600 m), mostly on the west, is bordered by shallow peatland, while mixed forest, i.e., larch (Larix sibirica), birch (Betula pendula), and Siberian pine (Pinus sibirica) borders about two-thirds of the perimeter (1100 m). In August 2016, sampling and monitoring stations were established at 13 locations along a west-east longitudinal transect and 5 locations along a transversal north-south transect (Fig. S1). CCH4 and CCO2 concentrations along the water column were determined using a membrane-integrated cavity output spectrometry method (M-ICOS)[16]. This method, described in more detail in the supporting information, allowed for the continuous measurement of dissolved gas at a frequency of 1 s−1, which corresponded to approximately 50 dissolved gas concentration data points per meter of the water column. The lower detection limit of the method under the present configuration was 5 nmol L−1 for CCH4 and 4 μmol L−1 for CCO2. Vertical CH4 fluxes and NMPR within the water column were derived from the estimation of turbulent diffusion of CH4 across the concentration gradient according to the method established in Kankaala et al.;[30] details of the method are provided in the supporting information. At each location, water column parameters including pH, temperature, DO, redox potential, and conductivity were also determined at 1 m depth intervals using multi-parametric probes (HI 9828, Hanna Instrument, Mexico). Water transparency was measured with a 30 cm Secchi disk. The euphotic depth (Z1%) was estimated according to methods established by French et al.[31] and LaPierre and Edmundson[32].

At the deepest location, water samples were taken at 1, 2, 4, 5, 6, 8, 9, and 10 m with a Van Dorn Bottle. These samples were used for molecular biology, stable isotope analysis of dissolved inorganic carbon (δ13C-DIC) and CH4 (δ2H-CH4 and δ13C-CH4). All stable isotope samples were analyzed in replicates of three and standard deviation was typically 0.2‰. The δ13C-DIC was analyzed by an Isoprime 100 unit (MultiFlow-Geo, Elementar, UK). Stable isotopic analysis of CH4 (δ2H and δ13C) was completed at UC Davis’ Stable Isotope Facility with a Thermo Scientific PreCon unit interfaced to a Thermo Scientific Delta V Plus Isotope Ratio Mass Spectrometer (Bremen, Germany). The fractionation factor α indicates the magnitude of isotopic separation between the δ13C values of ΣCO2 (total dissolved inorganic carbon) and CH4 in anaerobic environments. Its value reflects the dominant CH4 production pathways. It also shows systematic shifts in CH4 oxidation processes[33]. It was calculated as: α = (δ13C-DIC + 1000)/(δ13C-CH4 + 1000)[18]. For molecular biology characterization, water samples were filtered as soon as possible after sampling (i.e. <24 h). DNA extraction was completed and quantitative PCR (qPCR) was performed to assess the abundance of the following genes: bacterial 16 S rRNA gene (total bacteria), archaeal 16 S rRNA gene (total archaea), pmoA gene (particulate methane monooxygenase) and mcrA gene (methyl coenzyme M reductase, indicative of methanogens and ANMEs). More details of these methods are provided in the supporting information.

Supplementary Information.