Impacts of temperature and soil characteristics on methane production and oxidation in Arctic tundra

Rapid warming of Arctic ecosystems accelerates microbial decomposition of soil organic matter and leads to increased production of carbon dioxide (CO2) and methane (CH4). CH4 oxidation potentially mitigates CH4 emissions from permafrost regions, but it is still highly uncertain whether soils in high-latitude ecosystems will function as a net source or sink for CH4 in response to rising temperature and associated hydrological changes. We investigated CH4 production and oxidation potential in permafrost-affected soils from degraded ice-wedge polygons on the Barrow Environmental Observatory, Utqiaġvik (Barrow), Alaska, USA. Frozen soil cores from flat and high-centered polygons were sectioned into organic, transitional, and permafrost layers, and incubated at −2, +4 and +8 C to determine potential CH4 production and oxidation rates. Significant CH4 production was only observed from the suboxic transition layer and permafrost of flat-centered polygon soil. These two soil sections also exhibited highest CH4 oxidation potentials. Organic soils from relatively dry surface layers had the lowest CH4 oxidation potential compared to saturated transition layer and permafrost, contradicting our original assumptions. Low methanogenesis rates are due to low overall microbial activities measured as total anaerobic respiration and the competing iron-reduction process. Our results suggest that CH4 oxidation could offset CH4 production and limit surface CH4 emissions, in response to elevated temperature, and thus must be considered in model predictions of net CH4 fluxes in Arctic polygonal tundra. Future changes in temperature and soil saturation conditions are likely to divert electron flow to alternative electron acceptors and significantly alter CH4 production, which should also be considered in CH4 models. Copyright statement. This manuscript has been authored by UTBattelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).


Introduction
Arctic ecosystems store vast amounts of organic carbon in active layer soils and permafrost (Hugelius et al., 2014;Shiklomanov et al., 2010).Rising temperatures, increased annual thaw depth, and a prolonged thaw season accelerate microbial degradation of this carbon reservoir (Shiklomanov et al., 2010;Schuur et al., 2015;Schuur et al., 2013).The potential carbon loss due to these direct effects is estimated to be Published by Copernicus Publications on behalf of the European Geosciences Union.
92 ± 17 Pg carbon over the coming century (Schuur et al., 2015).The extent of soil organic matter (SOM) decomposition and partitioning between CO 2 and CH 4 emissions highly depend upon soil saturation conditions.Thawing of ground ice and ice-wedge degradation cause ground subsidence and significant changes in soil water saturation (Liljedahl et al., 2016) generating heterogeneous surface CH 4 fluxes (Schädel et al., 2016), with a current estimate of net CH 4 exchange from tundra to the atmosphere ranging widely from 8 to 29 Tg C yr −1 (McGuire et al., 2012).Understanding the factors that control CH 4 fluxes is key to reducing model uncertainties and predicting future climate feedbacks.
The unique polygonal ground in Arctic coastal plain tundra creates natural gradients in hydrology, snow pack depth and density, and soil organic carbon storage (SOC) that control CH 4 fluxes (Liljedahl et al., 2016;Lara et al., 2015).Thermal contraction processes create cracks in the tundra soil, which can fill with water that freezes to produce massive ground ice (French, 2007).This ice forms wedges that create the borders of three dominant polygon types, defined by their surface relief and subsurface hydrology: low-centered polygons (LCPs), flat-centered polygons (FCPs), and highcentered polygons (HCPs) (MacKay, 2000).Poorly drained LCPs are characterized by wet centers bordered by raised, relatively dry rims and wet troughs.FCPs lack the rims of LCPs and are drier (Wainwright et al., 2015).When ice wedges erode and water drains from the polygons, troughs subside and rims disappear to form drier HCPs.Methane emissions from wet and inundated LCP sites were 1-2 orders of magnitude larger than emissions from drier FCP and HCP sites (Vaughn et al., 2016;Sachs et al., 2010).Although a number of factors -including vegetation height and plant composition (von Fischer et al., 2010), soil inundation (Sturtevant et al., 2012), thaw depth (Sturtevant and Oechel, 2013;Grant et al., 2017), and season (Chang et al., 2014) -were suggested as explanatory factors for CH 4 flux variations, the 100-fold difference in CH 4 flux between polygon types could not be fully explained by variations in moisture or temperature (Sachs et al., 2010;Vaughn et al., 2016).Measurements of dissolved CH 4 concentrations in soil pore waters suggested a huge disconnect between > 100 µM concentrations of dissolved CH 4 in the active layer below 20 cm and negligible dissolved CH 4 at 10 cm soil depth (∼ 1 µM).This CH 4 gradient is particularly steep in FCPs.δ 13 C-CH 4 data suggested CH 4 oxidation played an important role in mitigating CH 4 in soil pore water and limiting surface CH 4 emissions (Vaughn et al., 2016).However, CH 4 oxidation potential is rarely studied in permafrost-affected soils.Therefore, we investigated CH 4 cycling in FCP soil for this study.
Methane oxidation mitigates terrestrial CH 4 emission.Up to 90 % of CH 4 produced in the soil is consumed in the upper dry layers of soil by aerobic CH 4 oxidizing bacteria (methanotrophs) before reaching the atmosphere (Le Mer and Roger, 2001).Methane oxidation rates are usually greatest in oxic, surficial soils, although methane oxidation is known to occur under oxygen-limiting conditions as well (Roslev and King, 1996).In the classical model of CH 4 oxidation profiles, there is usually a vertical gradient of decreasing O 2 concentration in the top centimeters of the soil column that is inversely correlated with an increasing gradient of CH 4 through the suboxic/anoxic active layer.Methanogens are relatively abundant in the deeper layer where oxygen is limiting (Lee et al., 2015;Kim and Liesack, 2015).Methane oxidation potential in wetlands and peat bogs is highest near the water table, while most CH 4 is produced below the water table (Whalen and Reeburgh, 2000).In contrast, methanogenic and methanotrophic communities can overlap in the rhizosphere (Liebner et al., 2012;Knoblauch et al., 2015), where roots or Sphagnum create an oxic/anoxic interface providing a substantial amount of oxygen for methanotrophs (Laanbroek, 2010;Parmentier et al., 2011) and organic acid substrates for methanogenesis.
Soil CH 4 fluxes result from the net effect of microbial CH 4 production and oxidation, coupled with transport processes.The rates of CH 4 oxidation are mainly governed by the abundance and composition of methanotrophic microbial communities and environmental factors including CH 4 and O 2 availability, soil air-filled porosity, and soil-water content (Preuss et al., 2013).Previous studies of boreal lakes and wetlands showed that CH 4 production is more sensitive to temperature changes than CH 4 oxidation, as CH 4 oxidation rates respond more strongly to CH 4 availability than temperature increase (Liikanen et al., 2002;Segers, 1998).Soils at the Barrow Environmental Observatory in Utqiaġvik (Barrow), Alaska, experience a wide range of arctic temperatures, from −20 to +4 • C (Shiklomanov et al., 2010).Soil respiration and methanogenesis continue at low temperatures close to 0 • C, even after the soil surface freezes trapping gas under ice.Therefore, substantial annual CH 4 and CO 2 emissions from the Alaskan Arctic occur during the spring thaw (Commane et al., 2017;Raz-Yaseef et al., 2017;Zona et al., 2016).However, it is unclear how accelerated warming in Arctic soils affects the opposing processes of CH 4 production and oxidation due to their nonlinear response to temperature changes (Treat et al., 2015).
In this study, we investigated the rates and temperature sensitivities of CH 4 production and oxidation from permafrost-affected soils in Utqiaġvik (Barrow), Alaska.Although various studies have identified significant and frequently correlated factors affecting CH 4 and CO 2 production in permafrost ecosystems upon thawing, the oxidation of CH 4 is not considered in most incubation studies (Treat et al., 2015).Here we used the natural geomorphic gradient of FCP and HCP soils to represent degraded polygon tundra soils with relatively oxic active layers that will potentially act as CH 4 sinks.Methane production and oxidation assays were performed separately using anoxic or oxic incubations at three temperatures.We also measured additional mechanisms influencing CH 4 production, including accumulation of organic acids and competing anaerobic respiration through iron reduction (Lipson et al., 2010).Specifically, we tested the following hypotheses regarding CH 4 dynamics in FCP and HCP soils: (I) CH 4 production is localized in the more reduced subsurface while CH 4 oxidation occurs at the soil surface; (II) CH 4 production is more sensitive to temperature increase than CH 4 oxidation and will likely exceed CH 4 oxidation in wet areas in response to warming.

Site description and soil sampling
The study site is located at the Barrow Environmental Observatory (BEO), Utqiaġvik (Barrow), Alaska, as part of the intensive study site areas B (high-centered polygon, HCP) and C (flat-centered polygon, FCP) of the Next Generation Ecosystem Experiments in the Arctic project.The centers of HCPs in area B were covered by lichens, moss, and dry tundra graminoids, while centers of FCPs in area C hosted wet tundra graminoids, mosses, and bare ground (Langford et al., 2016).Intact frozen soil cores from the centers of a water-saturated FCP (71 • 16.791 N, 156 • 35.990W) and a well-drained HCP (71 • 16.757 N, 156 • 36.288W) were collected with a modified SIPRE auger containing a sterilized liner (3 in.inner diameter), driven by a hydraulic drill during a field campaign in April 2012 (Herndon et al., 2015a, b;Roy Chowdhury et al., 2015).All samples were kept frozen during core retrieval, storage, and shipment to Oak Ridge National Laboratory (Oak Ridge, TN, USA).The frozen cores were stored at −20 • C until processing.The thaw depth measured in September 2012 in the HCP center was 40 cm, and thaw depths in the FCP varied from 41 to 47 cm.
The frozen soil cores were inspected and processed inside an anaerobic chamber (Coy Laboratory Products, MI, USA; H 2 ≤ 2 % and O 2 < 1 ppm).Both cores were sectioned into 10 cm segments, and each segment was inspected for evidences of roots and undecomposed organic matter.Soil Munsell color was recorded to qualitatively infer redoximorphic conditions.The soil core collected from the FCP center showed evidence of buried and discontinuous organic matter at approximately 50 cm depth, below the active, organic layer and above the permafrost.This transition layer was attributed to episodic thawing and cryoturbation (Schuur et al., 2008;Bockheim, 2007).Soil geochemical properties were measured for each 10 cm segment to assess soil geochemical depth profiles.Specifically, we measured soil gravimetric water content and soil pH, using the 1 : 2 soil slurry method within 1 M KCl.Fe(II) concentrations were measured to characterize redox conditions in the soil segments.Briefly, a soil subsample (∼ 2 g) was extracted with 0.1 M KCl solution for 30 min in an anaerobic chamber.The extracts were filtered with 2 µm syringe filters and analyzed immediately using the colorimetric 1,10-phenanthroline method (Hach method 8146).Absorbance was determined at 510 nm using a DU 800 spectrophotometer (Beckman Coulter, CA, USA).

Soil pore water gas measurements
Dissolved gas (CO 2 and CH 4 ) concentrations in soil pore water were determined at each 10 cm depth interval from the FCP and HCP cores.A 1 : 1 (w : v) soil slurry was prepared by mixing 10 g wet soil in 10 mL deionized and degassed water under anoxic conditions inside an anaerobic chamber.The samples were then placed in 15 mL crimp-sealed serum vials (Wheaton, NJ, USA).Vials were inverted and shaken at 4 • C for 12 h to allow for exchange between the dissolved and soil gas phases.Then, ∼ 5 mL of the aqueous phase was exchanged with Ar (99.9 % purity) using a gastight syringe (Hamilton, NV, USA), and samples were manually shaken vigorously for 5 min to allow for equilibration between the aqueous and gas phases.Subsequently, a 500 µL headspace sample was drawn and immediately analyzed with an SRI 8610C gas chromatograph using the method previously described (Roy Chowdhury et al., 2015).The detection limit for CH 4 was 1 ppm v .Concentrations of CH 4 and CO 2 were corrected for dissolved gases based on temperature and pHdependent Henry's law constants (Sander, 2015).

Low temperature respiration experiments
To investigate the temperature response of CH 4 production and overall organic carbon mineralization rates, samples of FCP and HCP soil cores were incubated at −2, +4, or +8 • C for approximately 90 days.Based on measured geochemical similarities, core segments were combined to represent the organic layer that comprises most of the active layer, transition layer (only present in FCP), and permafrost.Soils from each layer were homogenized inside an anaerobic chamber using sterile tools and equipment to establish microcosms (Roy Chowdhury et al., 2015).Our homogenization method using an oscillating power tool inside a glove box disrupts large, frozen soil clumps and removes gravel or litter, allowing us to place representative soil samples inside the serum vials.This method does not affect the soil microaggregate structure or expose the samples to oxidation, drying, or significant warming that could disrupt anaerobic microbial processes.Gases for the microcosm headspaces were selected based on gravimetric water contents and concentrations of reduced Fe(II) (Howeler and Bouldin, 1971) to best represent field conditions.Thus, organic layer soils from both FCP and HCP cores were incubated under air, while the transition layer and permafrost of FCP and HCP cores were incubated under anoxic conditions with N 2 headspaces.All vials were sealed with butyl rubber septa and crimp sealed to prevent evaporation and gas exchange.Anoxic microcosms were flushed with N 2 three times after sealing to remove residual H 2 and O 2 from the headspace.It is important to note that at −2 • C soil water remained unfrozen in these sam-ples due to freezing point depression (Romanovsky and Osterkamp, 2000).
A subsample of the combined, homogenized soil (∼ 10 g) was placed in a sterile 70 mL serum bottle sealed with a blue butyl stopper and aluminum crimp seal to form a soil microcosm.For each homogenized soil layer, nine replicate microcosms were constructed at each incubation temperature (Fig. 1).Headspace CO 2 and CH 4 concentrations were measured at 2-to 15-day intervals using gas chromatography (see Sect. 2.2).After 5, 10, or 20 days of incubation, three replicated microcosms were destructively sampled for methane oxidation assays (see Sect. 2.4 and Fig. 1) and additional soil geochemical analysis (Zheng et al., 2017).
Subsamples from microcosms opened after 10, 20, and 90 days incubation time were further processed for analysis of pH, Fe(II), and total organic acids.The pH and Fe(II) concentrations were measured as described in Sect.2.1.Iron reduction rates were estimated by the changes in measured Fe(II) concentration.Organic acids were analyzed using NH 4 HCO 3 extracts.Briefly, NH 4 HCO 3 extractions (12 h) were centrifuged for 15 min at 6500 g, then the supernatants were filtered through 0.2 µm membrane filters before analysis.Filtered samples were analyzed for low-molecularweight organic acids using a Dionex ICS-5000 + system (Thermo Fisher Scientific, MA, USA) equipped with an Ion-Pac AS11-HC column with a KOH mobile phase.Total organic acid concentrations (T OA ) were calculated by summing molar C equivalents for each measured species, normalized per gram of dry soil.

Methane oxidation potential assay
Soil samples were incubated in oxic conditions supplemented with ample CH 4 substrate to measure CH 4 oxidation potential (Roy Chowdhury et al., 2014).Methanogenesis is expected to be negligible under these oxic conditions.The methane oxidation assays (MOAs) were constructed using both freshly thawed (labeled as 0 day) and pre-incubated (labeled as various days of pre-incubation) samples to account for potential delays in the overall microbial activities (Fig. 1).Replicated samples (about 2 g) were slurried in a 1 : 1 (w : v) ratio with autoclaved deionized water in 26 mL serum bottles under ambient conditions.A 1 % CH 4 headspace was introduced into each crimp-sealed bottle by replacing 0.23 mL of headspace with 99.99 % CH 4 (Scott, Air Liquide).For freshly thawed samples, nine replicates were constructed for each incubation temperature and each soil layer.For pre-incubated samples at the designated sampling day, three replicated incubations (∼ 10 g per sample) were destructively sampled to construct nine replicates of MOAs to incubate at the same incubation temperature they were pre-incubated (Fig. 1).Due to the limited number of shaking incubators, only the 4 • C MOA from FCP and 8 • C MOA from HCP were shaken to minimize potential gasliquid phase transfer limitations.Headspace CO 2 and CH 4 concentrations in MOAs were measured at 2-to 15-day time intervals using gas chromatography (Sect.2.2).

Rate estimation, temperature sensitivity, and statistical analyses
Concentrations of CH 4 and CO 2 from soil microcosms were fitted with hyperbolic, sigmoidal, or linear functions (Roy Chowdhury et al., 2015).Rates of CO 2 production were calculated using derivatives of the best curve-fitting equations with parameters listed in Supplement Table S3.Methane oxidation rates were calculated from the loss of headspace CH 4 , which were best fitted with simple linear regression.All rate calculations are reported on per gram soil dry weight basis (denoted dwt in figures).The temperature dependence was calculated using the conventional Q 10 relationship by taking the ratio of maximum production or oxidation rates at 8 and −2 • C based on triplicate measurements.Changes in soil physicochemical properties were evaluated with one-way ANOVA, Tukey's honest significant difference (HSD) test.The effect of soil layers (organic, transition layer, and permafrost) and incubation temperature (−2, 4, and 8 • C) were examined with Tukey's HSD test.All curve fittings and statistical analyses are performed with R 3.4.0(The R Foundation for Statistical Computing) and validated with Prism (GraphPad Software, ver.7.0a).

Calculation of net CH 4 emission
To evaluate the net result of CH 4 production and oxidation and how this result changes in response to temperature increase, we applied a simplified model simulation.Representation of the CH 4 oxidation rate (R oxi ) is based on Michaelis-Menten kinetics with linear dependence on the biomass of methanotrophs (Xu et al., 2015), while the CH 4 production rate (R pro ) is calculated from measurements directly where  3 Results

Soil attributes and pore water characteristics
Soil cores from FCP and HCP center positions showed distinct vertical profiles of soil moisture expressed as gravimetric water content (g g −1 dry soil).The soil core from FCP was characterized by a wet surface within the top 10 cm below ground, a much drier organic layer between 10 and 40 cm, and a bottom layer below 40 cm with significantly higher water content.In the HCP core, soil moisture gradually increased from the top to the bottom (Fig. 2).Soil bulk density (g cm −3 ) is negatively correlated with gravimetric water content in both FCP and HCP along the depth profile (R 2 = −0.93 and R 2 = −0.86,respectively).A similar water distribution has been recorded by continuous field measurements of volumetric water content at the nearby NGEE_BRW_C soil pit monitoring site (http://permafrostwatch.org, last access: 3 November 2018.).Fe(II) concentration showed a strong positive correlation with gravimetric water content in both FCP and HCP cores (R 2 = 0.81 and R 2 = 0.91, respectively).The soil pH in FCP increased steadily with soil depth (R 2 = 0.95), with an average of 4.7.In the HCP soil core, soil pH varied by 1.5 pH units, with an average of 5.4.Overall, soil moisture, Fe(II) concentration, and pH increased with depth in the centers of both polygon types.
Dissolved CO 2 in soil pore water showed a similar general trend in both FCP and HCP cores.The concentration of dissolved CO 2 increased from 100 µM in the surface soil (0-10 cm) to approximately 950 µM at 30 to 40 cm depth in FCP, and it decreased below 40 cm.In HCP, approximately 400 µM dissolved CO 2 was measured in the surface soil.The concentration was 400-500 µM in the top 30 cm and also decreased significantly below 40 cm.A strong correlation between dissolved CO 2 concentration and soil bulk density was observed for both FCP and HCP (R 2 = 0.91 and R 2 = 0.85, respectively).
The highest dissolved CH 4 concentration (about 85 µM) was found between 30 and 40 cm in soil pore water of FCP, approximately 10 times the CH 4 concentration measured from the top 10 cm and 2-4 times higher than the CH 4 concentration measured below 40 cm.In HCP, significant CH 4 accumulation in soil pore water was found below 50 cm of the HCP core, while no dissolved CH 4 was detected above 50 cm.
Soil cores of FCP and HCP were divided into organic, transitional, and permafrost layers to facilitate incubation setup.The top 10 cm of both cores contained mostly plant material, litter, and ice or snow: these sections contained little soil and were not studied further.The organic layers of FCP and HCP cores were both oxic, with low Fe(II) concentrations and minimal dissolved CH 4 in soil pore water, while deeper layers were more reduced with 5-to 7-fold higher Fe(II) concentrations and more dissolved CH 4 (Figs. 1  and S1).Measured total carbon and nitrogen content showed distinct patterns in FCP and HCP cores (Table S1).The total carbon content of the FCP permafrost (31 %) was nearly twice as large as the organic layer.The FCP transition layer contained much less carbon than the adjacent layers (20 % of FCP permafrost), leading to a low C/N ratio of 16.For the HCP, the total carbon contents of organic and permafrost layers were 21 % and 17 %, respectively, significantly lower than that of the FCP permafrost.Inorganic carbon quantified as CO 2 released upon acid treatment was less than 0.001 % for each layer of the FCP and HCP cores.

Temperature responses of CH 4 production and oxidation
Thawed FCP and HCP soil samples were incubated in microcosms at fixed temperatures to assess methanogenesis rates.CH 4 production was only observed in microcosms from the transition layer and permafrost of FCP, which were incubated under anoxic conditions.CH 4 production started within 5 days after the anoxic incubations were set up (Fig. 3).Cumulative CH 4 concentrations at all temperatures were best fitted with a linear model (Table S2).Soils from the transition layer yielded about 10 times more CH 4 than permafrost at same incubation temperatures.Transitionsoil showed a stronger temperature effect than permafrost.CH 4 production rate increased by 1.6 and 3.1 times as the temperature increased from −2 to 4 and 8 • C, respectively.Measurements of CH 4 concentrations in the headspace of HCP permafrost were mostly below the quantification limit of the gas chromatograph flame ionization detector (1 ppm v , Fig. S2 in the Supplement).Organic soils from both FCP and HCP were incubated with air to best represent the field condition, and methanogenesis was unlikely to occur under oxic conditions.Headspace O 2 was not completely consumed after 90 days incubation (calculations not shown).Potential rates of aerobic CH 4 oxidation were measured in freshly thawed and pre-incubated soils to minimize total mi- crobial growth limitations.Pre-incubated soils showed much higher CH 4 oxidation potential compared to corresponding freshly thawed soil in both FCP and HCP.In FCP soil, CH 4 oxidation potentials measured in soils from the transition layer and permafrost were significantly higher than those measured in organic soils (Fig. 4a, b).Similarly, permafrost from HCP showed higher CH 4 oxidation potentials than organic soils (Fig. 4c, d).Overall, CH 4 oxidation rates in HCP soils were 80 %-90 % lower than rates from the equivalent FCP soil layers.
Rates of both CH 4 production and oxidation responded positively to temperature increase (Fig. 5).In the transition layer of FCP, CH 4 production showed much higher temperature sensitivity than CH 4 oxidation from both freshly thawed and pre-incubated soils, with an estimated Q 10 value of 4.1.The Q 10 value of CH 4 oxidation was 2.0 from freshly thawed soils and only 1.1 in pre-incubated soils.Similarly, the Q 10 value for CH 4 oxidation in permafrost also dropped from 1.7 in freshly thawed soils to 1.0 in pre-incubated soils.However, CH 4 production in the permafrost responded slowly to temperature increase, with an estimated Q 10 value of 1.7.Overall, permafrost showed significantly lower CH 4 production rates than soils from the upper transition layer, and also had a much lower temperature sensitivity for CH 4 production.However, the measured CH 4 oxidation potentials were similar in both soils, with similar temperature responses.

Soil respiration in response to rising temperature
Total soil respiration was evaluated using CO 2 production to characterize the observed variation in CH 4 production.Soils in all of the microcosm incubations produced CO 2 , by aerobic respiration under oxic conditions or by anaerobic respiration and fermentation under anoxic conditions.CO 2 production started immediately in microcosms of FCP samples, including organic soils that were incubated under oxic conditions and soils from transition layer and permafrost samples that were incubated under anoxic conditions (Fig. 6).CO 2 accumulation was best modeled by a hyperbolic function, except the organic layer soil incubated at −2 or +4 • C, where the best fit was a linear function (Table S3, Fig. 6).This exception is likely due to continuous aerobic respiration at lower temperatures, indicating that substrate limitation was not reached within 90 days.The transitional layer had the slowest CO 2 production rates and least carbon loss via CO 2 , in contrast to its relatively high rate of methanogenesis.
CO 2 production in microcosms of HCP samples was significantly delayed with much lower production compared to FCP soils.The HCP cumulative CO 2 production profiles were best fitted with a sigmoidal model, compared to the hyperbolic model that best fit FCP data (Table S3).A prolonged delay in CO 2 accumulation was observed in both HCP organic and permafrost samples.CO 2 production started about 10 days after the microcosm setup in the organic layer and reached a maximum rate at 30 days (+4 and +8 • C) or 75 days (−2 • C).The delay was longer in permafrost incubations, with a rapid increase after 40 to 50 days followed by a plateau.Therefore, microorganisms mineralized more carbon from FCP soils than HCP soils, and CO 2 production began sooner in FCP than HCP soil incubations.
Temperature showed significant effects on CO 2 production from the organic and transitional layers of FCP during 90-day incubations (p < 0.01 for each layer, Fig. 6).FCP soils incubated at +8 • C produced substantially more CO 2 than those incubated at lower temperatures.In the permafrost layer of FCP, CO 2 production was significantly higher at +8 • C compared to 2 • C (p < 0.05, ANOVA with Tukey's multiple comparisons test).Anaerobic CO 2 production in FCP permafrost showed much higher temperature sensitivity than that from FCP transition layer (Q 10 = 2.3-3.3 and 1.2-1.3,respectively, Table S4).

Organic acids production and iron reduction in FCP soils
Organic acids were produced as intermediate metabolites during microbial degradation of organic matter, and they probably fueled methanogenesis and iron reduction.Specific organic acids were analyzed from soil extracts from FCP transitional and permafrost layers (Table S5).and oxalate, consistent with previous analyses from LCP soils (Herndon et al., 2015a).Trace amount of lactate, pyruvate, and succinate were detected, with concentrations less than 0.05 µmol g −1 soil.To compare the changes in organic acids over time, total organic carbon contained in formate, acetate, propionate, butyrate, and oxalate products was calculated (T OA , µmol C g −1 ; Fig. 7a, b).Concentrations of dominant organic acids measured in the permafrost were approximately 10 times higher than those measured in the transition layer.The difference could be partly explained by much lower SOC content (5.8 %) in the transitional layer than that in the permafrost (30.8 %).T OA increased by 23 %, 70 %, and 65 % at −2, 4, and 8 • C, respectively, in transitional soils during 90-day anoxic incubations.Permafrost initially contained a much higher concentration of organic acids, and T OA increased by a lower percentage, by 1 %, 15 %, and 25 % at −2, 4, and 8 • C, respectively.Significant increases in Fe(II) concentrations were observed from anoxic incubations of FCP samples over the incubation period (Fig. 7c, d).In soils from the transition layer, Fe(II) concentrations stayed at a similar level during the first 20 days of anoxic incubation, and then increased significantly from ∼ 50 to ∼ 100 µmol g −1 between 20 and 90 days of incubation at −2, 4, and 8 • C. In the permafrost, the highest level of iron reduction within the first 20 days was observed in samples incubated at −2 • C. Between 20 and 90 days, the highest iron reduction rates were observed in samples incubated at 8 • C. The estimated Q 10 values were 1.2 and 1.3 for transition layer and permafrost, respectively.If we assume that iron reduction is coupled to acetate oxidation to produce CO 2 , stoichiometric calculations suggest that iron reduction could account for 96 % and 70 % of CO 2 produced in the transitional and permafrost layers at −2 • C. At 8 • C iron reduction could account for 74 % and 61 % of acetate oxidation in the transitional and permafrost layers.

Discussion
The widespread ice-wedge degradation in the Arctic causes morphological succession and hydrological changes in tundra ecosystems (Liljedahl et al., 2016).FCP and HCP features represent successively more degraded polygons.Despite clear geomorphological differences between polygon types that affect drainage, vegetation, and snow cover, the frozen FCP and HCP organic layers share similar gravimetric water contents, pH, SOC, and Fe(II) concentrations.Permafrost from both polygons contains more water, dissolved CH 4 , and Fe(II) than organic layer soils indicating more reducing environments.Concentrations of CH 4 measured in soil pore water from FCP and HCP cores increased with depth.These results are consistent with field measurements of CH 4 dissolved in soil water sampled from the sites during the thaw season (Herndon et al., 2015b).This CH 4 gradient suggests CH 4 oxidation in the upper organic layer.Therefore, we developed the following hypotheses for CH 4 production and oxidation, which were informed by previous studies of methane cycling in temperate ecosystems but untested in the Arctic.(1) CH 4 is produced in the more reduced subsurface, and consumed by methane oxidizers at the upper section of the soil column where O 2 is available.(2) Methane production has a higher temperature sensitivity than CH 4 oxidation, and is likely to exceed the CH 4 consumption rate in wet areas in response to warmer temperature.
Aerobic CH 4 oxidation is usually assumed to be limited by O 2 and CH 4 diffusion.Therefore, upper layers of soil, the rhizosphere, and soil at the water table would be expected to have the highest CH 4 oxidation activities (Shukla et al., 2013;Gulledge et al., 1997).The abundance of the pmoA marker gene for CH 4 oxidation decreased with soil depth in HCP trough soils (Yang et al., 2017) and in permafrostaffected soils from the Canadian Arctic (Frank-Fahle et al., 2014), consistent with that conceptual model.Thus, the organic layers of FCP and HCP in the BEO tundra would have higher potential for CH 4 oxidation.However, the highest CH 4 oxidation potentials were observed in the transitional and permafrost layers of FCP, the only layers with active methanogenesis.This result suggests that most active methanotrophs are found in these deeper soil layers, where CH 4 is available as a carbon and energy source.
Our results demonstrated that CH 4 oxidation might not be primarily O 2 diffusion limited, but rather limited by the availability of CH 4 in this system.The highest CH 4 oxi- were calculated from the concentrations of individual organic acids (Table S5).Error bars for Fe(II) concentrations are ±1 standard deviation from three replicate incubations.dation potentials were measured below the rhizosphere, in suboxic layers where CH 4 has accumulated.Given that halfsaturation constants for CH 4 and O 2 used in methanotrophy models vary over 1-2 orders of magnitude (Segers, 1998;Riley et al., 2011), aerobic CH 4 oxidation could occur throughout much of the soil column, as advective, diffusive, or plantmediated transport processes introduce O 2 into the soil.Others have observed deep soil CH 4 oxidation activity in peatlands (Hornibrook et al., 2009), fens (Cheema et al., 2015), and wet tundra (Barbier et al., 2012), often correlated with water table depth (Sundh et al., 1994).
The water table in the center of Barrow LCPs and FCPs varies somewhat during the thaw season but remains close to the surface (< 10 cm below surface for most of the thaw season; Liljedahl et al., 2015Liljedahl et al., , 2016)).Precipitation balances evapotranspiration during the thaw season, with little lateral runoff (Dingman et al., 1980), and volumetric water contents remain constant for these features (http://permafrostwatch. org, last access: 3 November 2018).In HCPs, the water table drops up to 20 cm below the surface following snowmelt (Liljedahl et al., 2015), and the soils have a lower volumetric water content.Due to limited drainage in the flat coastal plain, the frozen cores analyzed here are representative of field conditions for much of the thaw season.Water isotope analysis demonstrated that most water in the deep active layer comes from summer precipitation rather than seasonal ice melt (Throckmorton et al., 2016  age for Utqiaġvik (http://climate.gi.alaska.edu/,last access: 3 November 2018), suggesting a high water table and limited gas diffusion during the winter freeze-up before we collected soil cores in early 2012.As a result of these annual and transient changes in saturation, methanotrophs could colonize a broad range of the soil column below the rhizosphere.A comparison of methanogenesis and methane oxidation potential in peat bogs demonstrated that methanotrophs survived temporary exposure to anoxic conditions, suggesting these organisms can tolerate rapid changes in the water table and redox potential (Whalen and Reeburgh, 2000).These observations argue against hypotheses that CH 4 oxidation occurs primarily at the surface layer or at the water table interface.
The net effect of CH 4 production and oxidation in response to temperature change determines the sign of surfaceatmosphere CH 4 flux.The 4.1-fold increase in CH 4 production in the transition layer due to a 10 • C rise in temperature (Q 10 ) was similar to the average value of 4.26 reported in a recent meta-analysis of permafrost-affected soils (Schädel et al., 2016).These values are substantially higher than the temperature sensitivity of CH 4 oxidation from both freshly thawed and pre-incubated samples (Q 10 = 1 to 2).While both CH 4 production and oxidation respond positively to increased temperature, CH 4 production rates are predicted to increase more rapidly with higher temperature at this critical interface between the organic layer and permafrost.This difference in temperature sensitivity of CH 4 production and oxidation was also found in Arctic lakes (Lofton et al., 2014).CH 4 oxidation potential from freshly thawed FCP soils showed a temperature sensitivity coefficient (Q 10 ) between 1.7 and 2.0, which is consistent with reported values from peat (Segers, 1998).Higher Q 10 values for CH 4 oxidation were reported in drier mineral Arctic cryosols with low organic carbon content (Jørgensen et al., 2015;Christiansen et al., 2015).Q 10 relationships are strictly empirical representations of a processes' temperature sensitivity.Microbial community changes, differences in gene expression and enzyme production, and enzyme kinetics are all grouped together in this term, along with any changes in sorption or transport rates.Current process-enabled models do not distinguish between these factors, but future field experiments and model structures could provide a framework to differentiate mechanisms.
We simulated the net effect of CH 4 production and oxidation based on rate measurements from the transitional layer of FCP that exhibited the highest CH 4 oxidation potential.Methane oxidation rates are 14, 9, and 7 times the methanogenesis rate at −2, 4, and 8 • C, respectively.It is quite likely that methanogenesis will outpace CH 4 oxidation under much warmer temperature if there is no change in soil water content.We introduced model representations to explore this possibility, as previous studies suggest CH 4 oxidation rate is strongly regulated by the CH 4 supply (Liikanen et al., 2002;Lofton et al., 2014).Excluding diffusion in this thin layer, we modeled the distribution of the active biomass ratio between methanogens and methanotrophs (B methanotrophs /B methanogens ), which would cause Arctic soils to act as CH 4 sink or source in response to rising temperature (Fig. 8).Given the large difference in the rates of CH 4 production and oxidation, CH 4 oxidation would easily exceed methane production even with an active biomass ratio, B methanotrophs /B methanogens , lower than 1.A much lower biomass ratio, B methanotrophs /B methanogens (0.07-0.26), would make the soil a net source of CH 4 .Studies of functional genes involved in methane production (mcrA) and oxidation (pmoA) in the active layers in the western Canadian Arctic region suggested substantial variation in the pmoA/mcrA abundance ratio: the range is between 8.5 × 10 −5 and 7.6 × 10 2 (Frank-Fahle et al., 2014).Thus, accurate prediction of soil CH 4 production requires quantification of methanogen and methanotroph populations as model constraints.These results suggest the importance of parameterizing the temperature response function and biomass growth function specifically for methanogenesis and methane oxidation in model simulations to determine if the rates of methanogenesis and methane oxidation offset each other.These results support the hypothesis that methanogenesis can be offset by high methane oxidation rates in degraded tundra, although comprehensive simulations of surface gas fluxes will need to account for gas transport through the full soil column.
The incubations of thawed FCP and HCP soils from all layers revealed substantial differences in the temporal dynamics of CO 2 production.An initial lag was observed from HCP samples at all temperatures, suggesting low initial microbial activity, which might be due to low initial microbial biomass or substrate limitation.In contrast, rapid CO 2 production in FCP soils was observed from both the organic layer incubated under oxic conditions and the transitional and permafrost layers incubated under anoxic conditions.This temporal pattern of rapid accumulation after thawing was also observed for anaerobic respiration from LCP soils and HCP trough soils (Roy Chowdhury et al., 2015;Yang et al., 2016).
Coupled iron reduction and organic carbon oxidation processes made substantial contributions to total anaerobic CO 2 production relative to fermentation and methanogenesis.Acetate was the most abundant intermediate and exhibited the most dynamic concentration changes among individual organic acids measured from the soils.Using reaction stoichiometry for acetoclastic methanogenesis and anaerobic respiration through iron reduction (Istok et al., 2010), we estimated the amount of acetate being consumed by these parallel processes in soils from the transition layer and permafrost (Fig. 9).In transition-layer soils, about half of the available acetate was consumed by methanogenesis during the first 20 days of incubation (1.34 µmol g −1 out of 2.54 µmol g −1 ), while Fe(III) reduction consumed approximately 23 % of available acetate (0.60 out of 2.54 µmol g −1 ).In contrast, the partitioning of methanogenesis and Fe(III) reduction inverted during 20 to 90 days of anoxic incubation.Methanogenesis and Fe(III) reduction consumed 31 % and 59 % of available acetate, respectively.Permafrost con-tained a higher level of available acetate at the beginning of the incubation, and over 60 % of the available acetate was rapidly consumed by Fe(III) reduction during the first 20 days of incubation.From 20 to 90 days of incubation, iron reduction still consumed over half of the depleted acetate in the permafrost.This estimate is consistent with our initial characterization of the FCP core, where the soils from the transition layer contained the highest dissolved CH 4 concentrations, and permafrost was associated with significantly higher Fe(II) concentrations (Fig. S1).If hydrogen or other organic anions such as formate or propionate were oxidized by methanogens or iron reducers, then estimated acetate production levels would decrease slightly.These simulations indicate that Fe(III) reduction is responsible for most of the acetate mineralized and CO 2 produced in these soil incubations.As a significant anaerobic respiratory process in Arctic soils, Fe(III) reduction should also be included in CH 4 models to better predict greenhouse gas production in response to a changing climate.
Based on these findings, we propose the following scheme of soil carbon biogeochemistry in the FCP (Fig. 10).(1) Increasing temperature facilitates aerobic decomposition of organic carbon in the organic layer, and accelerates anaerobic carbon decomposition in the lower active layer and transition layer through fermentation, iron reduction, and methanogenesis to form CH 4 and CO 2 .(2) CH 4 produced in the transitional and permafrost layers is oxidized close to the site of production or transported to the atmosphere.(3) Fe(III) reduction is the primary anaerobic process responsible for the depletion of acetate, the major SOC decomposition intermediate.In this scheme, the oxic-anoxic interface could dynamically move in the soil column with changes in the water table and pore water distribution.Although the transitional layer contained much less carbon than the permafrost layer, the total carbon loss as the sum of CO 2 and CH 4 was comparable to that from permafrost.Laboratory measurements suggested that acetate accumulated in the organic layer could be transported into deeper layers to support iron reduction and methanogenesis (Yang et al., 2016).This transport might occur through vertical movement of dissolved organic compounds or mixing through cryoturbation (Drake et al., 2015).We will use results from these incubation experiments to structure and parameterize a thermodynamically based microbial growth model that improves simulations of anaerobic organic matter decomposition with CO 2 and CH 4 production.

Conclusions
Increased warming is predicted to accelerate the transition from wet LCPs to drier FCPs and HCPs in Arctic tundra, which probably function as potential CH 4 sinks.This study demonstrated that CH 4 oxidation capacity was tightly linked to CH 4 availability.Thus, the zone of highest CH 4 oxida-tion potential is at the suboxic area near the FCP transition layer and the upper permafrost.The measured CH 4 oxidation potential is an order of magnitude higher than the methanogenesis rate.With higher CH 4 residence time in the soil column due to limited gas diffusion in the field, CH 4 oxidation could easily consume CH 4 produced in deep permafrost soil at warming temperature.Given that iron reduction-coupled respiration predominates anaerobic organic carbon decomposition, CO 2 is likely to remain the major form of carbon emission from degraded polygons.This finding provides critical information about the dynamics of CH 4 production and oxidation with increasing temperature that needs to be incorporated into Arctic terrestrial ecosystem models for better predictions.

Figure 1 .
Figure1.Soil incubation and methane oxidation assay (MOA) experimental design.Soils from the transitional layer of a FCP core were used here to illustrate how incubations and MOAs were constructed at a given temperature.Incubation replicates were destructively sampled after 10 or 20 days to set up MOA with pre-incubated soils.For example, nine replicate microcosms were created for the FCP transition-layer soil, incubated under anoxic conditions.CO 2 and CH 4 in the headspace were cumulatively measured.After 10 days, three were opened and subsampled for triplicate MOAs.After 20 days, another three were opened and subsampled for triplicate MOAs.

Figure 2 .
Figure 2. Depth profile of gravimetric water content, bulk density, Fe(II) concentration, pH, and soil pore water dissolved CO 2 and CH 4 concentrations from FCP (blue) and HCP (red) cores.Replicate measurements were plotted as high and low values in each soil section as blue filled and empty circles (FCP) and red triangles (HCP).Trend lines were plotted based on the average values.

Figure 3 .
Figure 3. CH 4 production in soil microcosms from (a) transition layer, and (b) permafrost layers of FCP at indicated temperatures.The temporal profiles of CH 4 production were best fitted with the linear regression model, and the shaded area represents the 95 % confidence interval.

Figure 4 .
Figure 4. CH 4 oxidation potential measured from soils incubated at the indicated temperatures from (a) FCP after 0 days, (b) FCP after 20 days, (c) HCP after 5 days, and (d) HCP after 10 days.Error bars indicate ±1 standard deviation from three replicate incubations.

Figure 5 .
Figure 5. Temperature sensitivity of CH 4 production and oxidation measured from the transition and permafrost layers of FCPs.Temperature responses of CH 4 oxidation rates from freshly thawed (O transition and O permafrost ) and pre-incubated (O20 transition and O20 permafrost ) soils were both estimated.

Figure 6 .
Figure 6.Cumulative CO 2 production in soil microcosms from FCP and HCP samples at indicated temperatures.

Figure 7 .
Figure 7. Changes in total organic acids carbon (a, b) and Fe(II) concentrations (c, d) in soils from transition layer and permafrost of FCP during anoxic incubations.Total organic acids (µmol C g −1 )were calculated from the concentrations of individual organic acids (TableS5).Error bars for Fe(II) concentrations are ±1 standard deviation from three replicate incubations.
Figure 8. Simulations of active biomass ratio B methanotrophs /B methanogens distribution for Arctic soils to act as CH 4 sink or source in response to rising temperature.Dissolved CH 4 and O 2 concentrations of 0.1 mM in soil pore water are assumed.Half saturation constants (K m,CH 4 and K m,O 2 ) represent the baseline value (grey), high (blue), and low (red) range for sensitivity analysis.Biomass ratios below the curves indicate the system could be a CH 4 source, while biomass ratios above the curves suggest a CH 4 sink.

Figure 9 .
Figure 9. Changes in acetate concentrations associated with production (white bars) and consumption by iron-reducing bacteria (red bars) or methanogens (blue bars) were estimated using stoichiometric calculations based on measurements of methane and Fe(II) during incubations from 0 to 20 days and from 20 to 90 days at the indicated incubation temperatures.

Figure 10 .
Figure 10.Conceptual model of aerobic and anaerobic soil organic carbon decomposition pathways and the release of CO 2 and CH 4 from a flat-centered polygon (FCP).Lines on the left marked as Fe(II), CO 2 , and CH 4 represent measurements across the soil column.
methanotrophs and B methanogens represent the estimated biomass of methanotrophs and methanogens, respectively.K m,CH 4 and K m,O 2 are the half-saturation coefficients (mM) with respect to CH 4 and O 2 concentrations, respectively.
2 concentrations were calculated from soil pore water dissolved gas measurement and soil air-filled porosity estimations.With the above parameters, we estimated the biomass ratio of methanogens to methanotrophs (B methanogens /B methanotrophs ) under both net CH 4 production ). Precipitation during September and October 2011 was above aver-