Zero to moderate methane emissions in a densely rooted , pristine Patagonian bog – biogeochemical controls as revealed from isotopic evidence

Peatlands are significant global methane (CH4) sources, but processes governing CH4 dynamics have been predominantly studied in the Northern Hemisphere. Southern hemispheric and tropical bogs can be dominated by cushionforming vascular plants (e.g. Astelia pumila, Donatia fascicularis). These cushion bogs are found in many (mostly southern) parts of the world but could also serve as extreme examples for densely rooted northern hemispheric bogs dominated by rushes and sedges. We report highly variable summer CH4 emissions from different microforms in a Patagonian cushion bog as determined by chamber measurements. Driving biogeochemical processes were identified from pore water profiles and carbon isotopic signatures. Intensive root activity throughout a rhizosphere stretching over 2 m in depth accompanied by molecular oxygen release created aerobic microsites in water-saturated peat, leading to a thorough CH4 oxidation (< 0.003 mmol L−1 pore water CH4, enriched in δC-CH4 by up to 10 ‰) and negligible emissions (0.09± 0.16 mmol CH4 m−2 d−1) from Astelia lawns. In sparsely or even non-rooted peat below adjacent pools pore water profile patterns similar to those obtained under Astelia lawns, which emitted very small amounts of CH4 (0.23± 0.25 mmol m−2 d−1), were found. Below the A. pumila rhizosphere pore water concentrations increased sharply to 0.40±0.25 mmol CH4 L−1 and CH4 was predominantly produced by hydrogenotrophic methanogenesis. A few Sphagnum lawns and – surprisingly – one lawn dominated by cushion-forming D. fascicularis were found to be local CH4 emission hotspots with up to 1.52±1.10 mmol CH4 m−2 d−1 presumably as root density and molecular oxygen release dropped below a certain threshold. The spatial distribution of root characteristics supposedly causing such a pronounced CH4 emission pattern was evaluated on a conceptual level aiming to exemplify scenarios in densely rooted bogs. We conclude that presence of cushion vegetation as a proxy for negligible CH4 emissions from cushion bogs needs to be interpreted with caution. Nevertheless, overall ecosystem CH4 emissions at our study site were probably minute compared to bog ecosystems worldwide and widely decoupled from environmental controls due to intensive root activity of A. pumila, for example. Published by Copernicus Publications on behalf of the European Geosciences Union. 542 W. Münchberger et al.: CH4 emissions in a densely rooted, pristine Patagonian bog


Introduction
Peatland ecosystems are significant natural methane (CH 4 ) sources on the global scale responsible for about 10 % of global annual CH4 emissions (Aselmann and Crutzen, 1989;Mikaloff Fletcher et al., 2004).CH 4 has a 28-fold higher global warming potential over a 100-year time horizon compared to carbon dioxide (CO 2 ) (IPCC, 2014) and, thus, processes governing CH 4 dynamics in peatlands receive much attention to estimate global greenhouse gas emissions.
The slowdown of decomposition and accumulation of C in rainwater-fed peatlands (bogs) results from a high water table, anoxic conditions, recalcitrant peat-forming litter, low temperatures, inactivation of oxidative enzymes and an accumulation of decomposition end products (Beer and Blodau, 2007;Limpens et al., 2008;Bonaiuti et al., 2017).Once methanogenic conditions are established, CH 4 production is mainly controlled by a substrate supply of labile organic matter (Hornibrook et al., 1997;Whalen, 2005).Then, slow diffusive transport along the concentration gradient from deep anoxic peat layers to the atmosphere leads to CH 4 release.In upper, unsaturated and oxic peat zones that typically extend only about a few decimetres (Whalen, 2005;Limpens et al., 2008), CH 4 might get consumed by methanotrophic microbes (Chasar et al., 2000;Whalen, 2005;Berger et al., 2018).Consequently, water table position and fluctuations strongly control the amount of emitted CH 4 (Blodau and Moore, 2003;Whalen, 2005).While CH 4 oxidation suppresses CH 4 emissions in diffusion-dominated systems, ebullition by fast release of gas bubbles or plant-mediated transport by aerenchymatic roots can substantially increase CH 4 emissions (Fechner-Levy and Hemond, 1996;Joabsson et al., 1999, and references therein;Chasar et al., 2000;Colmer, 2003;Whalen, 2005;Knoblauch et al., 2015;Burger et al., 2016;Berger et al., 2018).With a high share of ebullitive fluxes, vegetated or open water pools in peatlands are considered to be strong CH 4 emitters (Hamilton et al., 1994;Blodau 2002; Burger et al., 2016) that have, however, received less attention than the vegetated surfaces (Pelletier et al., 2014).Pools can even turn the peatlands' C balance into a source (Pelletier et al., 2014), but examples of low-emission pools have also been reported (Knoblauch et al., 2015).
The majority of studies elucidating CH 4 dynamics in peatlands have been conducted on the Northern Hemisphere (e.g.Blodau, 2002;Limpens et al., 2008;Yu, 2012, and references in these), whereas only little research deals with CH 4 dynamics in bogs on the Southern Hemisphere (hereafter termed southern bogs) (Broder et al., 2015).Contrary to northern bogs, southern bogs are not only vegetated by Sphagnum mosses, dwarf shrubs or a few graminoids such as the rushlike species Tetroncium magellanicum, but also by vascular plants (Kleinebecker et al., 2007) that can densely root upper peat layers (Fritz et al., 2011;Knorr et al., 2015).Another feature that can be found in both systems is pools with sparse or no vegetation (Fritz et al., 2011;González Garraza et al., 2018).The vascular plants, for instance of the genus Astelia (Asteliaceae), Donatia (Stylidiaceae) or Oreobolus (Cyperaceae), independently developed a cushion life form to protect from harsh climate in cold environments (Gibson and Kirkpatrick, 1985;Boucher et al., 2016).Bogs dominated by cushion-forming vascular plants (hereafter termed cushion bogs) can be found in many (mostly southern) parts of the world, for instance all along the high Andes (Coombes and Ramsay, 2001;Benavides et al., 2013;Fonkén, 2014) down to southernmost Patagonia (Ruthsatz and Villagran, 1991;Heusser, 1995;Kleinebecker et al., 2007;Grootjans et al., 2014), in the highlands of eastern Africa (Dullo et al., 2017) as well as in New Guinea (Hope, 2014), New Zealand, andTasmania (Gibson andKirkpatrick, 1985;Ruthsatz and Villagran, 1991) and some sub-Antarctic islands (Ruthsatz and Villagran, 1991).These cushion bogs could also be regarded as extreme examples for densely rooted northern bogs that are dominated by rushes and sedges.
To our knowledge, only two studies reported CH 4 emissions from southern vascular-plant-dominated bogs with quite inconsistent results.Over a New Zealand bog domi-nated by the evergreen "wire rush" Empodisma robustum, CH 4 emissions exceeded those commonly reported from northern bogs (Goodrich et al., 2015), while emissions from a Patagonian cushion bog dominated by Astelia pumila were negligible (Fritz et al., 2011).Both, high and negligible CH 4 emissions, were mainly explained by the extensive aerenchymous roots of the prevailing plant species: the aerenchyma might function as a conduit for CH 4 from deep peat layers to the atmosphere in the first case, whereas in the second case pronounced O 2 supply by roots might have resulted in a thorough oxidation of pore water CH 4 and enhanced organic matter decomposition.However, Fritz et al. (2011) were the first to investigate CH 4 dynamics below Astelia lawns and determined emissions only sporadically without considering all predominant microforms and plant communities.Thus, general conclusions about CH 4 emissions from this type of ecosystem are so far restricted.Furthermore, more knowledge on biogeochemical processes throughout the peat column of Astelia lawns as a model plant for cushion-forming aerenchymous plants and in other microforms is needed to extend our limited understanding of these first insights into CH 4 dynamics in cushion bogs.Despite the variety of bogs dominated by densely rooted aerenchymous plants in the Northern Hemisphere and Southern Hemisphere, there is some evidence from the latest research that very low CH 4 pore water concentrations could be a more common phenomenon in those ecosystems (Knorr et al., 2015;Dullo et al., 2017;Agethen et al., 2018).In general, promotion or suppression of CH 4 production should be determined by the ratio of root density and activity associated with O 2 release versus presence of labile root organic matter or exudates accompanied by O 2 consumption (Blodau, 2002;Agethen et al., 2018).
Here, we quantified austral summer CH 4 emissions from dominant microforms in a Patagonian cushion bog and examined possible environmental and biogeochemical controls.We hypothesized that (i) emissions from cushion plant vegetation dominated by Astelia pumila or Donatia fascicularis are negligible while pools and Sphagnum lawns emit CH 4 in considerable amounts; (ii) pore water CH 4 and dissolved inorganic carbon (DIC) concentration profiles and their carbon isotopic signatures below densely rooted Astelia lawns reflect a distinct CH 4 oxidation effect contrary to sparsely or even non-rooted peat below pools; and (iii) in addition to by methanotrophy, isotopic composition of pore water CH 4 could also be affected by a predominance of acetoclastic methanogenesis throughout the rhizosphere of Astelia lawns shifting towards CO 2 reduction with peat depth and in sparsely or even non-rooted peat below pools.Biogeochemical controls of CH 4 emissions from Sphagnum lawns were not the subject of this study, as they have been intensively studied in northern bogs (e.g.Hornibrook et al., 1997;Beer and Blodau, 2007;Steinmann et al., 2008;Corbett et al., 2013) and even southern bogs with mixed vegetation in the absence of notable cushion plant coverage (Broder et al., 2015).

Description of the study site
The study was conducted at a cushion bog on the Península Mitre in southernmost Patagonia, Tierra del Fuego (Moat, 54 • 58 S, 66 • 44 W; Fig. 1).Peatlands of the Península Mitre provide the today rare opportunity to study carbon dynamics under pristine conditions as they can be considered largely undisturbed by human activities such as drainage, agriculture or elevated atmospheric nitrogen deposition (Kleinebecker et al., 2008;Fritz et al., 2011;Grootjans et al., 2014;Paredes et al., 2014).About 45 % of this area is covered by peatlands (Iturraspe, 2012) while little is known about these ecosystems because of their poor accessibility (Iturraspe, 2012).
The study site at Moat belongs to a complex system of sloping mires, blanket bogs, fens and -in coastal areas -cushion bogs (Iturraspe, 2012;Borromei et al., 2014).The studied cushion bog is located in exposed proximity to the Beagle Channel with harsh winds (Grootjans et al., 2014), resulting in an oceanic climate with average daily temperatures of 6 • C and annual precipitation of 500 mm (Fritz, 2012).
Large areas of the cushion bogs in Moat are characterized by a comparatively flat microrelief and covered by the cushion-forming plants Astelia pumila and Donatia fascicularis.The vegetation of the studied bog is composed of a mosaic of different plant communities characterized by specific plant species with lawns of A. pumila being the prevailing microform.Sphagnum magellanicum or D. fascicularis grow in small lawns (patches of a few square metres) at pool margins where S. magellanicum supposedly benefits from protection against desiccation.Tetroncium magellanicum, a rush-like herb that does not form cushions but presumably establishes aerenchymous roots as it tolerates inundation (von Mehring, 2013), is associated predominantly with Sphagnum lawns but also with Astelia and Donatia lawns, or it sporadically inhabits pools.These sparsely vegetated pools of different size (< 0.5 m 2 -∼ 10 m 2 ) and ∼ 0.5 m in depth are embedded in Astelia lawns.Pool sediments are inhabited by cyanobacteria (Arsenault et al., 2018;González Garraza et al., 2018) and submerged Sphagnum mosses such as Sphagnum cuspidatum (Kleinebecker et al., 2007;Kip et al., 2012).Peat formation started ∼ 11 000 years ago (Borromei et al., 2014), and previously the bog was dominated by S. magellanicum while A. pumila invaded the area around 2600 years before present as determined by pollen analyses (Heusser, 1995).The peat below Astelia lawns is densely rooted with an average of 2.15 g DW L of peat −1 down to 1.7 m (Fritz et al., 2011).Mean root lifetimes of A. pumila have been estimated to be ∼ 3-4 years, indicating a high root turnover (Knorr et al., 2015).

Sampling and analysis of solid peat and root biomass
Peat coring in Astelia and Sphagnum lawns (data not shown) was carried out using a Russian peat corer (Eijkelkamp Agrisearch Equipment, Giesbeek, the Netherlands) down to a maximum depth of 7.5 m in Astelia lawns.Roots of A. pumila in Astelia lawn cores were sorted out and treated separately.
The density of D. fascicularis root biomass was determined by cutting three 0.1 × 0.1 m sods down to a depth of 0.4 m.Peat and root biomass samples were oven dried at 70 • C until a constant weight to calculate peat bulk density and porosity as well as D. fascicularis root density.Total C and N contents together with abundance of 15 N and 13 C in the solid peat were determined using an elemental analyser (EA 3000, EuroVector, Redavalle, Italy) connected to an isotope ratio mass spectrometer (IRMS, Nu Instruments, Wrexham, UK).

Environmental variables
Environmental A cylindrical, transparent chamber with a basal area of 0.13 m 2 and a height of 0.4 m was used for soil-atmosphere flux measurements.The chamber was equipped with a fan to ensure mixing of the headspace air, a PAR sensor (HOBO S-LIA-M003, Onset, USA) and a temperature sensor (HOBO S-TMB-M0x, Onset, USA).A second temperature sensor recorded ambient temperature (HOBO S-TMB-M0x, Onset, USA) and all sensor data were logged in 1 s intervals to a data logger (HOBO U30 NRC, Onset, USA).This set-up allowed us to control air temperature inside the chamber within an approximately 3 • C deviation of the ambient air temperature.If necessary, ice packs were installed inside the chamber to avoid temperature increase.Due to strong wind conditions we decided against the installation of a vent tube.An opening in the top of the chamber avoided overpressure during chamber placement.This opening was closed immediately after the chamber was gently placed on a collar for at least 3 min to conduct measurements.Collars were equipped with a water-filled rim to ensure a gas-tight seal between chamber and collar during measurements.Pool fluxes were performed with a floating chamber of identical dimension and design on four pools located at each platform.The chamber wall extended approximately 0.04 m into the water during measurements.
All measurements were performed under a broad range of irradiance between 07:00 and 22:00 (local time) to determine diurnal variations in CH 4 fluxes.At least once during each sampling period of 2 to 3 consecutive days, all selected microforms and pools were measured under dark conditions.Opaque, reflective material was used to cover the chamber for dark conditions.Measurements were taken every 2 weeks during the first campaign in 2015 (N = 405 measurements) and on two occasions in 2016 (N = 132 measurements).A more regular measurement routine throughout the whole season was constrained by harsh weather conditions and the remoteness of our study site.
The chamber was connected by a 2 mm inner diameter polyethylene tubing to a gas analyser (Los Gatos ultraportable greenhouse gas analyser 915-001, Los Gatos Research) to record the increase in CH 4 concentrations over time at a rate of 1 Hz.Between measurements, atmospheric background concentrations were achieved inside the chamber.The gas analyser was equipped with an external pump providing a flow rate of 2 L min −1 and pumping the analysed gas back to the chamber, creating a closed system.Prior to each campaign, the instrument was recalibrated.
We estimated the cover of plant species within each collar to the nearest 5 % and counted the number of T. magellanicum shoots to further characterize the different microforms (Table 1).In 2015, vegetation surveys were conducted at the beginning and at the end of the measurement campaign.Since no change in vegetation was observed, surveys were conducted only once during the campaign in 2016.For each collar, the position above the water table was determined once during each sampling period of consecutive days to account for different positions in the microrelief and to determine a specific water table position.Additionally, we determined the area of each individual microform patch in which a collar had been installed.

CH 4 flux calculation and statistical analyses
Soil-atmosphere CH 4 fluxes were calculated from the gas concentration increase over time within the chamber using the software package (MATLAB Release R2015a) routine described in Eckhardt and Kutzbach (2016) and Kutzbach et al. (2007).Gas concentration was modelled as either a linear or exponential function of time.Model performance was compared using Akaike's information criterion (AIC) as a measure of goodness of fit.This routine resulted in 90 % of cases in a linear function of time (N = 485) and in 10 % of cases in an exponential function of time (N = 52).
We visually inspected all concentration increases over time and none of the measurements showed a stepwise concentration increase indicative of ebullient events.Of the 3 min measurement time, 50-180 s (mostly around 90 s) was selected for the CH 4 flux calculation to exclude unstable conditions due to chamber placement, particularly at the beginning of the measurement.Fluxes with a slope not significantly different from 0 (tested with an F test, p > 0.05) were set to zero (20 % of cases, N = 105).
Normal distributions of CH 4 flux datasets were checked by a Kolmogorov-Smirnov test.Due to non-normality, the Kruskal-Wallis analysis of variance followed by multiple comparison U tests adjusted by the Bonferroni method were applied to determine significant differences among CH 4 fluxes of different microforms and collars.The relationships between CH 4 fluxes and environmental variables (soil temperature, water table depth fluctuations) as well as additional features of individual collars were evaluated by Spearman's rank correlation (Table 2).A response of CH 4 emissions to the investigated environmental variables was identified only for a few individual collars (see the Supplement).Therefore, no seasonal or annual ecosystem emissions were calculated.Statistical analyses were performed using MATLAB (MAT-LAB Statistics Toolbox Release R2017a).All flux data were reported as mean ± standard deviation.

Sampling and field measurements
Multilevel piezometers (MLPs), as described in Beer and Blodau (2007), were installed during both measurement campaigns in austral summer 2015 and 2016 to determine pore water concentrations in depth profiles.The MLPs provided a spatial sampling resolution of either 0.1 or 0.2 m and were equipped with diffusive equilibration samplers made of permeable silicone tubes providing a 4 mL gas sample volume together with a 5 mL crimp vial filled with deionized water and covered with a permeable membrane (Supor 200 0.2 µ, Pall Corporation, Pall Life Sciences).To collect gas samples, MLPs were stepwise retrieved and the gas volume extracted from the equilibration sampler with a 3 mL syringe was transferred into nitrogen-flushed 5 mL crimp vials.The www.biogeosciences.net/16/541/2019/Biogeosciences, 16, 541-559, 2019 crimp vials were capped with a butyl stopper with an aluminium crimp seal.MLPs were installed in Astelia lawns and pools in three replicates at each platform where closed chamber measurements were performed.
In January 2015, MLPs were installed in Astelia lawns (0.2 m resolution, 2 m depth) and pools (0.1 m resolution from 0 to 1 m in depth, 0.2 m resolution below down to a maximum depth of 4 m) for 5 weeks of equilibration time.From gas samples pore water concentrations of CH 4 , DIC and hydrogen (H 2 ) were determined.To measure sulfate concentrations, which are potential electron acceptors originating from sea spray and affecting competitiveness of methanogenesis (Broder et al., 2015), a pore water subsample obtained from each crimp vial was transferred to a storage vial and stored frozen until analysis.An aliquot of each remaining pore water sample was measured in situ with a pH / EC (electrical conductivity) meter (Combo HI 98129/130, Hannah Instruments, Germany).
In January 2016, MLPs were installed to a depth of 3 m in pools (resolution as in 2015) for 3 weeks of equilibration time and afterwards again in Astelia lawns (0.1 m resolution in transition zone from rooted to non-rooted peat at 1-2 m in depth) for 9 weeks of equilibration time.In 2016, equilibration samplers were used to collect CH 4 and CO 2 samples and determine their carbon stable isotopic signatures.Pore water samples obtained from crimp vials were used to measure O 2 concentrations in situ with a planar trace oxygen minisensor (Fibox 3, PreSens Precision Sensing GmbH, Germany).

Analytical procedures
Gas samples were transported to Germany and analysed within 4 weeks.Gaseous CH 4 and CO 2 concentrations were measured with a gas chromatograph (8610C, SRI Instruments, USA) equipped with a methanizer and flame ionization detector (FID).Hydrogen concentrations were analysed on a H 2 analyser (Ametek ta3000 H 2 analyser, Trace Analytical TA 3000r).Sulfate concentrations were obtained by ion chromatography (883 Basic IC plus, Metrohm, Herisau, Switzerland).
Stable carbon isotopic signatures of CH 4 and CO 2 were simultaneously determined by cavity ring-down spectroscopy (CRDS; Picarro G2201-i connected to a small sample isotope module (SSIM, model A0314), Picarro Inc., USA).As this set-up required a minimum sample volume of 20 mL, samples were diluted with nitrogen prior to measurements.The instrument was calibrated in the beginning of every measurement day using two working standards of CH 4 (1000 ppm, −42.48 ‰) and CO 2 (1000 ppm, −31.07 ‰).Isotopic signatures are given in δ notation relative to Vienna Peedee Belemnite (VPDB).As samples were stored for several weeks and δ 13 C-CO 2 values were biased in case of high CH 4 concentrations in the sample, a correction procedure was applied following Berger et al. (2018).

Calculations and statistical analyses
Pore water concentrations of CH 4 , DIC and H 2 were recalculated from gaseous concentrations obtained from equilibration samplers by applying Henry's law following Eq.( 1): where c is the concentration in moles per litre, p the pressure in atmospheres and K H the Henry constant that was corrected according to Sander (1999) for mean soil temperatures of 10 • C in February 2015 and 2016 as well as 7 • C in April at a depth of 0.5 m.DIC concentrations at prevalent pH conditions were calculated with respect to carbonate speciation considering equilibrium constants following Stumm and Morgan (1996).All pore water data were reported as mean ± standard deviation (N = 3).Zones of CH 4 production and consumption in the peat column were visually identified based on observed pore water concentration gradients.For a quantitative evaluation of pore water concentration gradients, steady-state conditions as well as a dominance of diffusive gas transport would have to be assumed; this may only partly apply to the system under  study here, given the large root biomass.Nevertheless, we additionally applied a software routine using inverse modelling (PROFILE, Berg et al., 1998) for zone identification.This modelling approach supported the results obtained visually, but provided further rough estimates for production and consumption zones due to the complex diffusivity in the rhizosphere of highly rooted peat and, thus, results are presented in the Supplement only.

Characteristics of solid peat and root biomass
Roots of A. pumila were present in the upper profile down to a depth of 1.8 m (rhizosphere) in Astelia lawn cores.Throughout the rhizosphere, we observed a highly decomposed and amorphous peat supposedly originating from A. pumila that developed above Sphagnum peat.Peat below Sphagnum lawns was continuously formed predominantly by Sphagnum.Total C and N contents in the peat together with natural abundance of 15 N and 13 C are given in Table 3. Root density of D. fascicularis was 0.014 ± 0.010 g DW L of peat −1 at 0.35 m in depth while root density of A. pumila was not determined once again as it has been already quantified with > 4 g DW L of peat −1 at around 0.5 m in depth (Fritz et al., 2011).

Environmental conditions and potential controls on CH 4 fluxes
The study site was characterized by daily mean PAR values up to ∼ 700 µmol m −2 s −1 with maximum values exceeding 2000 µmol m −2 s −1 and daily mean air temperatures ranging from −0.5 • C to 17 • C during the measurement campaigns in austral summer from January to April (Fig. 2).Soil temperature in summer reached maximum daily average values of 15 and 12  temperature fluctuations could be identified (Table 2).Yet, significant relationships among water table fluctuations or variations in soil temperatures were established for those few lawns with considerable emissions (Table 2).

Soil-atmosphere CH 4 fluxes and features of microforms possibly affecting CH 4 emissions
Summer CH 4 emissions in a south Patagonian cushion bog were highly variable among the dominant microforms Astelia lawns, Sphagnum lawns and Donatia lawns characterized by specific plant species as well as pools and showed a pronounced spatial pattern even within microforms.Compared to all other microforms, only Sphagnum lawns showed considerable emissions with 1.52 ± 1.10 mmol CH 4 m −2 d −1 while Astelia lawns emitted nearly no CH 4 with 0.09 ± 0.16 mmol m −2 d −1 (Fig. 3).Methane emissions from pools were also low (0.23 ± 0.25 mmol m −2 d −1 ) and not significantly different from intermediate Donatia lawn fluxes (0.66 ± 0.96 mmol CH 4 m −2 d −1 ).During some occasions, we determined even significant, negative CH 4 fluxes, mainly from Astelia lawns (Fig. 3a, Table 2).Most of these negative fluxes as well as the majority of fluxes (41 %) that were set to zero were obtained from one individual collar throughout  the whole measurement campaign precluding the possibility of a measurement artefact.
As indicated by high standard deviations, CH 4 emissions among replicates of all microforms were not consistent.Field observations suggested that CH 4 emissions of a microform were not only controlled by the predominant plant species, but might also have been associated with additional features of the respective microform.Therefore additional microform features (i.e.number of T. magellanicum shoots, cover of D. fascicularis, extent of the microform in which a specific collar was installed) were assessed.However, as we did not expect such small-scale spatial variability among emissions of individual collars, the results of this survey given in Table 2 are of rather explorative, preliminary character as the low number of replicates (three to five collars) in combination with several outliers barely allows any sound statistical analysis.Nevertheless, collars with elevated or comparatively high emissions had some fea-

Pore water CH 4 and DIC concentration profiles
Pore water profiles in peat columns below Astelia lawns and pools showed similar trends in pore water concentrations and during both sampling years.CH 4 concentrations were almost zero in the upper pore water profile below Astelia lawns (< 0.003 mmol L −1 ) down to a depth of around 1.5 m, which corresponds to the zone where the rhizosphere was most pronounced (Fig. 4a).Below this depth, CH 4 concentrations sharply increased up to 0.25 ± 0.08 mmol L −1 at 3 m in depth.CH 4 concentrations in pore water profiles below pools resembled profiles obtained under Astelia lawns on elevated levels.Throughout the upper profile down to 1.5 m, CH 4 concentrations were ∼ 0.03 mmol L −1 with a peak around 0.3 m in depth (0.06 ± 0.06 mmol L −1 ) and increased steeply to 0.40 ± 0.25 mmol L −1 at 3 m in depth.Maximum concentrations in comparable depths reached similar levels in both sampling years below Astelia lawns, but maximum concentrations below pools were about 3 times higher in 2015 compared to 2016.DIC predominantly occurred as dissolved CO 2 because of the low pH ranging from 3.36 to 4.77.Contrary to CH 4 , DIC concentrations increased constantly with depth from around 1 mmol L −1 near the surface to 2.60±1.0mmol L −1 at 2 m in depth below Astelia lawns and 2.94±1.1 mmol L −1 at 3 m in depth below pools (Fig. 4b).With depth, DIC converged with CH 4 concentrations, which was reflected by DIC : CH 4 ratios that were extremely high in the rhizosphere below Astelia lawns exceeding 100 (Fig. 4c).Beneath the rhizosphere, ratios steeply approached values below 40.Under pools, ra-tios slightly decreased with depth but were mostly around 40 down to 1.5 m with two distinct peaks of very little CH 4 at the surface and around 0.8 m in depth.In deep peat layers below the rhizosphere, DIC : CH 4 ratios detected under Astelia lawns and pools converged and reached the lowest ratios of ∼ 10 and ∼ 5, respectively.
Visual inspection of concentration profiles but also modelling of CH 4 production and consumption rates indicated a predominance of CH 4 consumption throughout the whole rhizosphere below Astelia lawns and in roughly corresponding depths under pools (Figs.4a, 6).Maximum CH 4 consumption was identified in the lower rhizosphere around 1.5-2 m where the increase in CH 4 concentration was most pronounced.Deep peat layers below the rhizosphere of Astelia lawns and similar depths below pools were considered to be CH 4 sources.In these depths, the high concentration levels of CH 4 (up to 0.40±0.25 mmol L −1 ) sustain substantial upward diffusion of CH 4 following the concentration gradient into the consumption zone throughout the rhizosphere of Astelia lawns and in corresponding depths below pools.

Carbon isotopic values in pore water and apparent fractionation
Values of δ 13 C in CH 4 did not show a clear depth trend below Astelia lawns and ranged in the upper profile down to 2 m in depth between −87.2 ± 10.1 and −72.1 ± 10.3 ‰ (Fig. 4d).Isotopic δ 13 C-CH 4 values below pools became less negative with depth from −93.4 ± 9.2 to −73.7 ± 1.2 ‰, with a minimum peak in the upper peat layers around 0.3 m of −101.0 ± 5.6 ‰.In the upper metre of the profile, values of δ 13 C-CH 4 were thus up to 25 ‰ more negative below pools compared to Astelia lawns while the isotopic values of both profiles converged in deep peat layers below 2 m at −80 ‰.Values of δ 13 C in CO 2 below Astelia lawns became more negative with depth in the upper profile down to 1 m, reaching values as low as −32.4 ± 0.8 ‰, and increased to −19.5 ± 1.0 ‰ at 3 m in depth.Below pools, values became less negative with depth throughout the whole profile from −21.2 ± 0.4 to −6.0 ± 1.5 ‰ (Fig. 4e).Apparent fractionation factors varied only slightly and showed opposite trends below the two microforms.Values increased slightly from 1.060 to 1.066 below Astelia lawns, peaking at 1.045 at 0.5 m in depth.Below pools, apparent fractionation factors decreased with depth from 1.079 to 1.073 with a maximum of 1.091 around 0.3 m in depth (Fig. 4f).

Hydrogen, oxygen and sulfate concentrations in pore water profiles
Hydrogen concentrations were mostly below 3 nmol L −1 below Astelia lawns while they were elevated below pools, ranging between 3 and 10 nmol L −1 (Fig. 5a).Maximum concentrations were reached in the upper profile of Astelia lawns with up to 11.83 ± 17.88 nmol L −1 and, in addition to a peak at the surface, below pools from 1 to 2 m in depth with up to 41.61 ± 64.63 nmol L −1 .Molecular oxygen concentrations were mostly below 5 % saturation and tended to be higher in the upper profile down to 2 m below Astelia lawns compared to pools (Fig. 5b).Sulfate concentrations were similar below both, Astelia lawns and pools, and reached the highest values near the surface with up to 16.43 ± 11.85 µmol L −1 .With depth, concentrations decreased and approached zero below 1 m (Fig. 5c).

Discussion
This study is among the first dealing with CH 4 dynamics in an austral, vascular-plant-dominated cushion bog.We aimed to reveal patterns of CH 4 emissions and their environmental as well as potential below-ground biogeochemical controls.We furthermore attempted to elucidate the CH 4 dynamics in the peat below the four predominant microforms on a conceptual level based on the results presented here and in previous studies.

Environmental controls on CH 4 emissions
Summer CH 4 emissions of a south Patagonian cushion bog dominated by A. pumila were low, albeit spatial heterogeneity was pronounced with small patches of Sphagnum or Donatia lawns being local emission hotspots.The variation in CH 4 fluxes was mainly controlled by the dominant plant species of each lawn microform, which emitted significantly different amounts of CH 4 .We found only a few significant responses of CH 4 fluxes to water table fluctuations or variations in soil temperatures (Table 2).Such weak coupling between environmental variables and CH 4 emissions is not surprising as fluxes from most individual collars were relatively low.Furthermore, the water table fluctuated only slightly near the surface (Fig. 2) and may thus not serve as a primary control on CH 4 emissions.Stronger effects may only be expected when amplitude and duration of fluctuations affect peat redox conditions substantially (Blodau and Moore, 2003;Knorr et al., 2009).Under such stable moisture conditions and as long as methanogenesis is not limited by substrate supply, the sensitivity of methanogenic microbial consortia to temperature should become apparent (Whalen, 2005).Variations in landscape CH 4 fluxes have indeed been related to soil temperature previously (Rinne et al., 2007;Jackowicz-Korczyński et al., 2010;Goodrich et al., 2015), which is in accordance with the present study.The temporal variability in CH 4 emissions over the two investigated seasons was rather low (data not shown) due to low seasonality in both, temperature and precipitation, in the study region.Measurement campaigns extending over shoulder seasons and comparing years of contrasting weather conditions would be necessary to reveal the impact of more pro-nounced water table fluctuations or temperature regimes on CH 4 fluxes in south Patagonian cushion bogs.

Astelia lawns -zero emission scenario
Emissions from Astelia lawns were minute (0.09 ± 0.16 mmol m −2 d −1 ) and significantly lower than from all other microforms (Fig. 3), verifying our first hypothesis.The near-zero emissions were well explained by near-zero CH 4 concentrations (< 0.003 mmol L −1 ) in the rhizosphere down to around 1.8 m.Only below, concentrations increased sharply to 0.2 mmol L −1 at 2 m in depth (Fig. 4a).The low CH 4 concentrations even turned some Astelia lawns into a small sink for atmospheric CH 4 as shown by negative fluxes that were obtained sporadically throughout the whole measurement campaign, confirming earlier results (Kip et al., 2012).These findings agree with previous research by Fritz et al. (2011), who observed CH 4 emissions from Astelia lawns at a low magnitude compared to that presented here and presented a similar CH 4 concentration depth profile.Extremely low CH 4 concentrations have also been obtained by Dullo et al. (2017) in upper peat layers of an Ethiopian cushion bog.Nevertheless, in deep peat layers below the rhizosphere, CH 4 concentrations of our study approached magnitudes reported for Chilean bogs with mixed vegetation consisting of Sphagnum mosses and A. pumila (Broder et al., 2015).This is in the lower range of concentrations described from northern bogs mostly around 0.5 mmol L −1 (e.g.Blodau and Moore, 2003;Beer and Blodau, 2007;Beer et al., 2008;Corbett et al., 2013) but even reaching up to 5 mmol L −1 at comparable depths (Steinmann et al., 2008).
A CH 4 concentration profile as observed at our site suggests that in deep peat layers below the rhizosphere (> 2 m) CH 4 was produced at substantial rates and comparably low oxidation, as indicated by DIC : CH 4 ratios notably smaller than 40 (Fig. 4c).Subsequently, CH 4 became almost completely oxidized on its diffusive way to the atmosphere as visualized on a conceptual level (Fig. 6).Ratios of DIC : CH 4 under methanogenic conditions in northern bogs are typically around 5 (e.g.Hornibrook et al., 1997;Blodau and Moore, 2003;Corbett et al., 2013), but can be as high as 30 (Steinmann et al., 2008).Notably higher ratios throughout the rhizosphere suggested methanotrophic activity and availability of either alternative electron acceptors or O 2 from aerenchymous roots suppressing methanogenesis (Colmer, 2003;Mainiero and Kazda, 2005;Knorr et al., 2008a, Fritz et al., 2011;Dullo et al., 2017).Diffusion and plant-mediated transport as pathways for upward CH 4 transport would have caused higher emissions and ebullition requires supersaturation of at least ∼ 350 mmol CH 4 L −1 in the pore water (Fechner-Levy and Hemond, 1996;Beer et al., 2008).The sharp drop in CH 4 concentrations from −2 to −1.5 m in depth characterized this peat layer as a zone of the strongest methanotrophic activity, as supported by modelled consumption zones.
Carbon isotopic values obtained from pore water profiles reflected rhizosphere processes, although the δ 13 C-CH 4 signature alone did not provide a clear indication for oxidative effects in the rhizosphere as proposed in our second hypothesis.A strong effect of only methanotrophy should result in a shift of the mean δ 13 C-CH 4 signature to more enriched values and distinct changes in the fractionation factor from lower to upper horizons.Methane isotopic signatures were close to −80 ‰ in deep peat layers while above 2 m in depth values ranged from −90 to −70 ‰ throughout the rhizosphere (Fig. 4d) while the fractionation factor increased slightly (Fig. 4f).These surprisingly negative δ 13 C-CH 4 values throughout the rhizosphere along with a comparatively small shift and a wider standard variation in replicate samples are an unexpected result that needs further explanation.
We can only explain the pronounced variation in replicate δ 13 C-CH 4 values to be associated with the presence of aerenchymous roots and the mean δ 13 C-CH 4 depth pattern by a coexistence of aerobic and anaerobic microsites that would both support such strong variation within short spatial scales.The occurrence of small-scale aerobic microsites attached to roots releasing O 2 into water-saturated (peat) layers has been previously suggested (Colmer, 2003), also for peatlands in the Northern Hemisphere (Popp et al., 1999;Knorr et al., 2008a;Corbett et al., 2013).Life time and activity of A. pumila roots are highly dynamic (Knorr et al., 2015) and supposedly governed temporal and spatial expansion of microsites.Within these microsites, labile organic matter and CH 4 were rapidly consumed by aerobic respiration and methanotrophy, resulting in CH 4 enriched in 13 C by up to 10 ‰ (Fig. 4d).This δ 13 C-CH 4 signature and the apparently small 12 C/ 13 C fractionation observed throughout the rhizosphere (Fig. 4f) revealed the influence of a strongly fractionating process such as CH 4 oxidation.Although based on 13 C isotope fractionation alone CH 4 oxidation cannot be clearly separated from acetoclastic methanogenesis (Whiticar et al., 1986;Conrad, 2005), oxidative effects in the rhizosphere were very likely since diffusion or plant-mediated transport were excluded as sinks for CH 4 .Methanotrophy leaves CH 4 in the pore water enriched by ∼ 10 ‰ while produced CO 2 becomes depleted (Chasar et al., 2000;Popp et al., 1999) as observed here throughout the rhizosphere with fewer negative signatures occurring only at greater depths (Fig. 4e).Also, values of δ 13 C-CO 2 of less than −30 ‰ that were substantially lower than the source organic matter material around −26 ‰ (Table 3) can only be explained by the influence of 13 C-depleted carbon from oxidation of 13 C-depleted CH 4 .Substantial rates of root respiration were a further source for depleted CO 2 in the pore water (Corbett et al., 2013) as suggested by high DIC : CH 4 ratios mostly exceeding 100 throughout the prevailingly watersaturated rhizosphere (Fig. 4c).Nevertheless, O 2 release by roots still exceeded consumption, as demonstrated by O 2 saturation reaching up to 10 % (Fig. 5b).This O 2 supply may thus not only fuel CH 4 oxidation but also heterotrophic respiration using molecular oxygen or other electron acceptors regenerated by this O 2 input (Colmer, 2003;Mainiero and Kazda;2005).Likely further sources of CO 2 could have been dissimilatory sulfate reduction in anaerobic microsites of surface peat layers down to 0.4 m (Fig. 5c) with sulfate originating from sea spray (Kleinebecker et al., 2008;Broder et al., 2015).Taking these results obtained from CH 4 and DIC concentration profiles and their carbon isotopic signatures together, the second hypothesis stating that pore water concentrations below Astelia lawns reflect a distinct CH 4 oxidation effect was confirmed.
Enriched CH 4 being only a leftover from CH 4 oxidation should result in a δ 13 C-CH 4 depth pattern following that of the δ 13 C-CO 2 profile (Hornibrook et al., 2000), but such a pattern was not observed in our study (Fig. 4d and e).The surprisingly negative δ 13 C-CH 4 signature throughout the rhizosphere is likely explained by production of small amounts of CH 4 at anaerobic microsites in the absence of living roots.Either hydrogenotrophy or acetogenesis and subsequent acetoclastic methanogenesis cause large fractionation (Steinmann et al., 2008) and could add depleted CH 4 to the rhizosphere.As we did not quantify other parameters such as labile organic root matter or root exudates and acetate concentrations or their carbon isotopic signatures, we cannot clearly determine the relative importance of both pathways for our study.Elevated H 2 concentrations in surface peat layers (Fig. 5a) indicated high H 2 production by fermentation compared to consumption by sulfate reduction or hydrogenotrophic methanogenesis, presumably due to high root litter input and root exudates (Knorr et al., 2009;Estop-Aragonés et al., 2013).But as acetogens can outcompete hy- drogenotrophic methanogens at low temperatures, especially when H 2 is available, (Kotsyurbenko et al., 2001), we propose that acetoclastic methanogenesis was the more important pathway here.Together with the probably high amounts of labile organic matter available to serve as substrate for methanogenesis throughout 2 m of rhizosphere, this explanation would support our third hypothesis.With depth, root density presumably decreased (Fritz et al., 2011) and a slight increase in 12 C/ 13 C fractionation indicated a to a predominance of hydrogenotrophic methanogenesis and, below 2 m in depth, absence of microsites (Fig. 4f) as α c values of replicate samples varied less and increased towards 1.065 (Whiticar et al., 1986;Conrad, 2005).

Pools -low-emission scenario
Pools exhibited, with 0.23 ± 0.25 mmol CH 4 m −2 d −1 , the second lowest emissions of the microforms under study (Fig. 3).Such low emissions were contrary to our first hypothesis and somewhat surprising since they have been previously reported to be similar to Sphagnum lawns (Fritz et al., 2011) and were expected to be relevant CH 4 sources because of prevailing anoxic conditions (Blodau, 2002).Moreover, other studies identified pools as hotspots of CH 4 emissions from peatlands (Hamilton et al., 1994;Pelletier et al., 2014;Burger et al., 2016).This was clearly not the case here, as emissions were of the same order of magnitude as pool fluxes in south Patagonian Sphagnum bogs (Lehmann et al., 2016).
Low emissions were associated with CH 4 pore water concentrations surprisingly resembling those below Astelia lawns on slightly elevated levels (Fig. 4a).These low pore water concentrations in sparsely or even non-rooted peat below pools were not supported by our second hypothesis and require another explanation.The δ 13 C-CH 4 signature was depleted compared to Astelia lawns down to 2 m in depth and largely followed that of δ 13 C-CO 2 (Fig. 4d, e) and, thus, did not indicate a methanotrophic effect.As upward diffusion of CH 4 against the concentration gradient is impossible, we can only explain low pore water concentrations by lateral exchange of CH 4 .We suggest that root O 2 release by A. pumila controlled even CH 4 dynamics below adjacent microforms as visualized on a conceptual level in Fig. 6.The thorough CH 4 oxidation below Astelia lawns may establish a lateral concentration gradient of CH 4 .Thereby, CH 4 from adjacent microforms may diffuse to the rhizosphere of Astelia lawns where it was consumed.As the microrelief is not pronounced at our study site, we propose that there is a large depth zone throughout the 2 m rhizosphere of Astelia lawns and in corresponding depths below pools with negligible water flow and a predominance of diffusive transport.In this depth zone, we observed a highly decomposed and amorphous peat as previously described for cushion bogs (Ruthsatz and Villagran, 1991;Fritz et al., 2011), indicating a very high water permeability that would together with a high water table position facilitate such lateral exchange (Baird et al., 2016).The assumption of lateral exchange was supported by a peak in DIC : CH 4 ratios accompanied by very low CH 4 concentrations around 0.8 m below pools (Fig. 4c), inhibiting upward diffusion.Lateral root ingrowth and O 2 supply from adjacent Astelia lawns could furthermore suppress methanogenesis and keep CH 4 concentrations below pools low, in particular if pools are small.
Patterns of carbon isotopic signatures, levels of H 2 concentrations of ∼ 10 nmol L −1 sufficient to maintain hydrogenotrophic methanogenesis (Heimann et al., 2010) and low DIC : CH 4 ratios suggested methanogenic conditions over the whole profile (Figs.4c, 5c).Methanogenesis was especially indicated by a slight peak of CH 4 concentrations associated with lower DIC : CH 4 ratios and elevated H 2 ratios around 0.3 m in depth.However, even high activity of fermentation processes as indicated by peaking H 2 levels from 1 to 2 m in depth did not translate to changes in fractionation factor (Figs. 4f, 5a) but suggest that CH 4 was produced only very locally if at all.Low substrate supply from highly decomposed organic matter and thus unfavourable thermodynamic conditions seemed to limit CH 4 production as suggested by very negative δ 13 C-CH 4 signatures (Fig. 4d) accompanied by peaking α c values in the upper profile (Fig. 4f) (Hornibrook et al., 1997;Knorr et al., 2008b).Instead, methanogenesis might have been outcompeted by other electron-accepting processes, such as sulfate reduction (Fig. 5c).If produced at all, CH 4 was rapidly consumed around 0.3 m in depth by methanotrophic bacteria inhabiting submerged Sphagnum mosses that provided recalcitrant organic matter on pool bottoms (Kleinebecker et al., 2007;Kip et al., 2012;Knoblauch et al., 2015) or were associated with cyanobacteria (Arsenault et al., 2018;González Garraza et al., 2018).Within the open water of pools, which was well mixed by strong winds, CH 4 production and excess of H 2 were diminished due to electron acceptor availability for consumption (e.g.O 2 , sulfate).Thus, the pools at our study site belonged to examples of low-emission pools as described in Knoblauch et al. (2015).

Sphagnum and Donatia lawns -low-to moderate-emission scenario
CH 4 emissions from Sphagnum and Donatia lawns were most variable among microforms ranging from near zero to local emission hotspots.Sphagnum lawns showed significantly higher CH 4 emissions of 1.52 ± 1.10 mmol m −2 d −1 than all other microforms (Fig. 3).Substantially higher fluxes from Sphagnum lawns at our study site have been previously described by Fritz et al. (2011) and emissions reported here were on the same order of magnitude as determined for Sphagnum lawns in south Patagonian bogs dominated by Sphagnum mosses (Broder et al., 2015;Lehmann et al., 2016).The highest emissions obtained in the present study from Sphagnum and Donatia lawns were at intermediate levels compared to northern bogs (Blodau, 2002;Laine et al., 2007;Limpens et al., 2008).
We present the first CH 4 emission data obtained from Donatia lawns.Contrary to our expectations, Donatia lawns emitted 0.66±0.96mmol CH 4 m −2 d −1 and thus had the second highest CH 4 emissions of our study.One lawn even exceeded emissions of the strongest Sphagnum lawn source (Fig. 3c, Table 2) on all measurement occasions.As these comparatively high emissions could be reproduced well, we conclude that an artefact can be excluded.Therefore, our first hypothesis, that cushion plants emit negligible amounts of CH 4 , could retrospectively only be partly confirmed, as it depends on the specific cushion-forming plant species.We conclude that the presence of cushion plants as a proxy for negligible CH 4 emissions from cushion bogs needs to be carefully interpreted.
The spatial distribution of vascular plant root characteristics and its oxidative effects apparently controlled CH 4 dynamics at the study site as suggested by the pronounced spatial variability in emissions associated with our preliminary observations of microform features (Table 2).D. fascicularis established a deep and aerenchymous, but with < 0.05 g DW L of peat −1 , a substantially less dense rooting system compared to A. pumila with > 4 g DW L of peat −1 at around 0.5 m in depth (Fritz et al., 2011).The lower root density of D. fascicularis was probably below a specific threshold at which O 2 release did not suppress emissions any more.Instead, the roots might have accelerated CH 4 production through the presence of labile organic matter while facilitating emissions (e.g.Joabsson et al., 1999;Blodau, 2002;Colmer, 2003;Whalen, 2005;Agethen et al., 2018).As root O 2 release promotes nutrient mobilization (Colmer, 2003), this assumption of little if any O 2 release by D. fascicularis is in line with Schmidt et al. (2010), who found nutrient concentrations in D. fascicularis to be almost as low as in S. magellanicum while A. pumila contained significantly higher amounts of nutrients.
Sphagnum and Donatia lawns that showed substantial CH 4 emissions were beyond the sphere of influence by A. pumila roots.Although not measured here, pore water CH 4 concentrations below Sphagnum and Donatia lawns with remarkably higher emissions probably increased steeply below the water table (Fig. 6), like common depth profiles in northern (e.g.Beer and Blodau, 2007;Limpens et al., 2008;Corbett et al., 2013) or southern bogs (Broder et al., 2015).This explanation implies that the pore water below those lawns was not well connected to pore water influenced by A. pumila roots.A spatial variability in peat physical properties on the ecosystem scale was suggested by water table records (Fig. 2), demonstrating that one site responded faster than the other to, for instance, precipitation events.
We speculate that the biogeochemistry below Sphagnum and Donatia lawns with near-zero emissions was, similar to pools, influenced by root activity of nearby A. pumila (as visualized in Fig. 6).This assumption was supported by previous research at our study site.Fritz et al. (2011)  The highest summer emissions in the present study were considerably lower compared to maximum austral summer fluxes of ∼ 150 mg CH 4 m −2 d −1 (9.4 mmol m −2 d −1 ) determined for a New Zealand bog dominated by the vascular evergreen wire rush, Empodisma robustum (Goodrich et al., 2015).The authors explained such high CH 4 emissions by comparatively wet conditions and a high density of aerenchymous vegetation providing a gas conduit for CH 4 transport in their study site.Compared to CH 4 emissions of northern bogs ranging from 3 to 53 mg m −2 d −1 at intermediate levels in wet bogs to up to 80 mg m −2 d −1 (Blodau, 2002;Laine et al., 2007;Limpens et al., 2008; 0.2 to 2.2 and up to 5 mmol m −2 d −1 ), the highest summer emissions in the present study represent a pronounced CH 4 source.Nevertheless, overall ecosystem CH 4 emissions would probably be among the lowest described for pristine bog ecosystems worldwide when taking into account the high proportion of the surface area where Astelia lawns prevailed.While Donatia lawns covered only small parts of our study site, we observed other parts of the complex bog system in the whole study area to be widely dominated by D. fascicularis instead of A. pumila, which is in line with previous studies (Heusser, 1995;Grootjans et al., 2014) that described the cushion bogs in the Moat landscape to be covered by both, A. pumila and D. fascicularis.Thus, amounts of emitted CH 4 could also be significant on the ecosystem scale once D. fascicularis dominates.For instance, a predominance of D. fascicularis is typical for later vegetation succession stages in cushion bogs of Tierra del Fuego (Heusser, 1995) or Chile (Ruthsatz and Villagran, 1991;Kleinebecker et al., 2007).

Conclusion
We conclude from our study that the spatial distribution of cushion-forming A. pumila root density and activity and associated O 2 supply strongly controlled CH 4 production and consumption in a pristine Patagonian cushion bog.Thereby, CH 4 emissions were reduced to near-zero levels and largely decoupled from environmental controls.The high root density even regulated CH 4 dynamics below adjacent microforms with less or non-rooted peat such as pools, Sphagnum lawns or lawns dominated by cushion-forming D. fascicularis by maintaining lateral concentration gradients in upper peat layers.Nevertheless, when root density dropped below a certain threshold, CH 4 production might have been acceler-ated by the presence of labile root organic matter or exudates accompanied by O 2 consumption along with facilitated CH 4 transport by aerenchymatic roots of D. fascicularis.Under such circumstances, CH 4 emissions increased to intermediate levels compared to northern bogs.Therefore, the presence of cushion plants as a proxy for negligible CH 4 emissions from cushion bogs should be interpreted carefully.As cushion bogs can be found in many (mostly southern) parts of the world and only very limited knowledge about CH 4 dynamics from these systems exists, future research should take into account a possibly high spatial variability in CH 4 emission from bogs dominated by cushion plants.We demonstrated an extreme scenario for how a spatial distribution of root density in the peat can lead to a pronounced pattern of CH 4 emissions.Yet, the underlying ratio between root characteristics and O 2 supply determining this emission pattern should be applicable to other densely rooted peatlands in general.Further research should be undertaken to prove relationships that have been developed here on a conceptual level in order to extend our knowledge on CH 4 dynamics in under-researched cushion bogs.
Sabrina Knaack, Sina Berger, Simona Bonaiuti and Leandra Praetzel is greatly acknowledged.We thank Peter Sulmann for technical assistance.We are grateful for the valuable comments of the two anonymous referees.The support by Open Access Publication Fund of the University of Münster is acknowledged.
Edited by: Luo Yu Reviewed by: two anonymous referees

Figure 1 .
Figure 1.Location of the study area in southernmost Patagonia.The investigated cushion bog in Moat is located 130 km east of Ushuaia (Argentina).The map was created at http://www.simplemappr.net(last access: 11 December 2017) and in MATLAB (MATLAB Mapping Toolbox Release R2017a).
CH 4 fluxes determined from individual collars installed in dominant microforms of a Patagonian cushion bog.Number of flux measurements, their mean and standard deviation are presented.Significantly different fluxes from collars within one microform were identified by Kruskal-Wallis analysis of variance and multiple comparison U tests with Bonferroni adjustment.Results of computing Spearman's rank correlation coefficient ρ for the relationships between CH 4 fluxes and possible environmental controls, namely water table depth determined separately for each collar and soil temperature at 0.05 m in depth, are given.Furthermore, characteristics of collars and results of computing Spearman's rank correlation for the relationships among mean CH 4 fluxes during the second campaign and each feature are shown.Microform size refers to the extent of a microform in which a specific collar was installed.P values < 0.05 indicate a correlation significantly different from zero..: not significant; n.d.: not determined.a Despite a (relative) high absolute value of Spearman's rho, p values did not indicate a significant correlation, probably because of outliers.www.biogeosciences.net/16/541/2019/Biogeosciences, 16, 541-559, 2019

Figure 2 .
Figure 2. Time series of environmental variables during the study period.Photosynthetic active radiation (PAR, a, b) and air temperature (c, d) are presented as half-hourly and daily means during two austral summer periods.Soil temperature (e) and water table position below the surface (f) were continuously recorded and are shown as half-hourly means.Note the difference in response time between platform a and platform b where water table fluctuations were recorded.

Figure 3 .
Figure 3. Mean CH 4 emissions and their standard deviation determined from dominant microforms (a) and individual collars during the two measurement campaigns (b, c) in a Patagonian cushion bog.Different letter superscripts denote significantly different amounts of emitted CH 4 among microforms despite variability within microforms being substantial (a) or from collars within one microform (b, c) that were identified by Kruskal-Wallis analysis of variance and multiple comparison U tests with Bonferroni adjustment.

Figure 4 .
Figure 4. Pore water composition in depth profiles obtained from MLPs (N = 3) installed in Astelia lawns and pools during two sampling campaigns in austral summer in 2015 and 2016.Mean values and their standard deviation are presented.CH 4 (a) and DIC (b) concentrations as well as related DIC : CH 4 (c) ratios were determined during both sampling events while carbon isotope values of CH 4 (d) and DIC (e) and corresponding fractionation factors (f) were investigated only during the second campaign.DIC : CH 4 ratios are shown on a logarithmic scale without negative standard deviation for Astelia lawns and pools.Displayed depths represent the centre of a MLP segment and had to be standardized because of different sampling resolution in both years.The water table position and approximate maximum rooting depth (rhizosphere) refer to Astelia lawns only.Note the different scale of the y axis in the upper and lower three panels.

Figure 5 .
Figure 5. Profiles of H 2 (a), O 2 (b) and sulfate (c) pore water concentrations obtained from MLPs (N = 3) installed in Astelia lawns and pools during two sampling campaigns in austral summer in 2015 and 2016.Mean values and their standard deviation are presented.Displayed depths represent the centre of a MLP segment and had to be standardized because of different sampling resolution in both years.The water table position and approximate maximum rooting depth (rhizosphere) refer to Astelia lawns only.

Figure 6 .
Figure 6.Schematic CH 4 concentration profiles below the four dominant microforms in a Patagonian cushion bog.The shape of profiles was derived from data obtained in the present study together with those by Limpens et al. (2008), Fritz et al. (2011) and Broder et al. (2015).The size of arrows represents the relative magnitude of either a CH 4 or CO 2 flux following the diffusion gradient.The water table position is displayed by a dashed line near the surface.CH 4 oxidation by root O 2 release resulted in near-zero CH 4 pore water concentrations in densely rooted peat below Astelia lawns.With increasing distance from densely rooted peat, CH 4 oxidation decreased while CH 4 concentration increased until reaching a depth profile typical for those of northern peatlands with fewer if any roots in water-saturated peat layers.Exemplary sequence of microforms, derived from observed amounts of emitted CH 4 .Other sequences of microforms are reasonable as microforms occur in various combinations in the field.

Table 1 .
Characteristic plant species of dominant microforms in a Patagonian cushion bog.Species composition was determined during the second measurement campaign within collars where closed chamber measurements were conducted.Mean cover (%) and its standard deviation of characteristic were given and plant species with a mean cover < 5 % were listed with +.

Table 2 .
• C at 0.05 m in depth and 0.5 m in depth, respectively.Concurrently, the water table fluctuated close to the surface between −0.03 and −0.23 m and differed consistently by about a few centimetres between the two water table measurement sites in Astelia lawns.CH 4 fluxes approached zero in many cases, but differed significantly from zero down to absolute values of 0.01 mmol CH 4 m −2 d −1 .Fluxes below absolute values of 0.01 mmol CH 4 m −2 d −1 did not differ significantly from zero and were therefore set to zero.For collars exhibiting low-magnitude fluxes neither seasonal trends (data not shown) nor a relationship to water tableand soil

Table 3 .
Summary of peat characteristics obtained from two peat cores taken in Astelia lawns in a Patagonian cushion bog.Data were averaged over the respective depths sampled at 0.1 m resolution.
Münchberger et al.: CH 4emissions in a densely rooted, pristine Patagonian bog tures in common: emissions from one Astelia lawn were with ∼ 0.3 mmol CH 4 m −2 d −1 , (significantly) higher in both sampling years compared to fluxes measured from the other Astelia lawns.This Astelia lawn collar with elevated emissions was characterized by (i) the presence of many T. magellanicum shoots together with a high share of D. fascicularis and (ii) placed in a lawn with small extent surrounded by a mosaic of small pools and D. fascicularis patches (both smaller than 1 m 2 ).Surprisingly, one Donatia lawn was also a substantial CH 4 source of 2.10 ± 0.14 mmol m −2 d −1 , even exceeding the highest emissions observed from Sphagnum lawns on all measurement occasions in 2016.Similar to the Astelia lawn with elevated emissions, the collars of Donatia and Sphagnum lawns with high emissions were characterized by (i) a high number of T. magellanicum shoots and (ii) placed in lawns surrounded by small pools and other patches of D. fascicularis with no A. pumila nearby.Emissions obtained from one Sphagnum lawn were, with less than 0.3 mmol CH 4 m −2 d −1 in both sampling years, significantly lower compared to emissions from other Sphagnum lawns.
www.biogeosciences.net/16/541/2019/Biogeosciences, 16, 541-559, 2019 W. Vice versa, this collar with low Sphagnum emissions was characterized by (i) the lowest number of T. magellanicum shoots among all Sphagnum lawn collars and (ii) installed in a small Sphagnum lawn surrounded by A. pumila.
showed that below Sphagnum lawns pore water CH 4 concentrations can www.biogeosciences.net/16/541/2019/Biogeosciences, 16, 541-559, 2019 W. Münchberger et al.: CH 4 emissions in a densely rooted, pristine Patagonian bog be substantially lowered throughout rhizosphere-influenced peat layers.Further research should be undertaken to investigate why CH 4 emissions from Sphagnum and Donatia lawns appear to be highly variable in order to understand CH 4 dynamics in under-researched cushion bogs.