Edinburgh Explorer Methane Emissions from a Grassland-Wetland Complex in the Southern Peruvian Andes

: Wet organic-rich mineral and peat soils in the tropical Andes represent a potentially signiﬁcant, but little studied, source of methane to the atmosphere. Here we report the results of ﬁeld and laboratory measurements of soil–atmosphere methane exchange and associated environmental variables from freely draining upland and inundation prone wetland soils in a humid puna ecosystem in the Southeastern Andes of Peru. Between seasons and across the landscape soil–atmosphere exchange varied between uptake and emission. Notable hotspots of methane emission, peaking during the wet season, were observed from both upland and wetland soils with particularly strong emissions from moss-accumulating topographic lows. This variability was best explained by the inﬂuence of oxygen concentration on methane production in superﬁcial soil horizons.


Introduction
The availability of satellite retrievals of atmospheric methane (CH 4 ) concentration at the beginning of the 21st century fuelled renewed interest in understanding variability in the global atmospheric budget of this radiatively important greenhouse gas [1][2][3]. This work re-confirmed the findings of previous inverse modelling simulations using airborne or ground sampling networks that the tropics are a larger source of CH 4 to the atmosphere than previously thought and that current bottom-up source-sink inventories based on scaling of field observation or process-models poorly characterize landscape-atmosphere exchange in these regions [4][5][6]. In this respect, there has been a drive to better constrain not only the dynamics of traditional wetland sources, such as the soils and waters of tropical swamp and seasonally inundated floodplain forests [7][8][9], but also less well understood elements of the tropical CH 4 cycle. Such elements include transport by wetland trees [10,11], emissions from wet upland soils [12][13][14], abiotic degradation of foliar pectin [15,16], and the function of cryptic environments like the leaf axes of canopy epiphytes [17]. In this context, the highlands of the tropical Andes are of particular interest as the presence of organic-rich mineral soils and peatlands in high altitude montane ecosystems potentially represent a poorly documented, but significant component, of the tropical South American CH 4 budget [18][19][20].
The tropical Andes, extending from latitudes of 11 • N to 23 • S and elevations of 600 to 6962 m above sea level (asl), represents a hotspot for biodiversity and endemism [21] covering between high CH 4 and O 2 availability close to the limit of the water-table or in the rhizosphere of aerenchymatous plants that transport O 2 below the water-table to maintain aerobic root function [44,45]. The remainder is either emitted to the atmosphere following diffusion or ebullition through the soil to the surface or passes into plant roots and is emitted above-ground through stems and leaves [46][47][48]. Reflecting the balance among these processes, wetland-atmosphere CH 4 exchanges are predominately controlled by water-table depth, soil temperature and vegetation cover [30]. In contrast, the uptake of atmospheric CH 4 by well-aerated upland soils is driven by communities of low-capacity-high-affinity methanotrophs that exist at close to ambient CH 4 and O 2 concentrations [39]. Under such conditions, oxidation is unsaturated with respect to CH 4 but not O 2 [49,50] and variations in soil-atmosphere CH 4 exchange primarily arise from the influence of soil texture, structure, and water content, integrated by terms such as soil water-filled pore space (WFPS), on the rate at which CH 4 can diffuse from the atmosphere to niches occupied by methanotrophs in the soil [51,52]. Additional to this physical constraint, the rate of uptake is also sensitive to the size and composition of the methanotrophic communities involved and factors, such as nitrogen availability, which can either inhibit or promote methanothrophy [53,54].
Whilst these models are generally successful in explaining soil-atmosphere CH 4 exchange, it is also apparent that methanogenic archaea are present not only in anoxic wetland soils, but are also ubiquitous [55,56] and active [12,57] in oxic environments. Indeed, field observations of CH 4 emissions from the oxic soils of wet uplands [52,58] and wetlands [59,60] have required refinement of the conceptual models which describe controls on soil-atmosphere exchange [61] and consideration from a global perspective [13]. Methanogenesis within, and emission of CH 4 from, oxic soils has been explained by the presence of anoxic microsites that form through limitations on diffusion of O 2 into soil pores occluded by aggregates and saturation combined with the demand of aerobic respiration [52,62]. As in wetlands, methanotrophic communities utilising elevated CH 4 concentrations close to the interface between anoxic and oxic niches can consume the majority of that produced [12]. In situations where this activity is sufficient to consume all of the CH 4 produced, soils may act as sinks for atmospheric CH 4 despite considerable flows of C passing through the below-ground CH 4 cycle [61,63,64]. The cryptic nature of CH 4 cycling in heterogeneous soils has made it difficult to identify the proximal controls on soil-atmosphere exchange; however, the emerging picture from tracer studies is that CH 4 production and emissions to the atmosphere are best explained by the methanogenic fraction of total C mineralisation as a proxy for the availability of C to methanogenic communities [59,61,64].
Here we report net CH 4 fluxes and associated environmental conditions from a grassland-wetland complex in the humid puna of Southeastern Peru. The first year of the data shown here has been previously reported and indicated that the study site is a significant source of atmospheric CH 4 when compared to surrounding montane forests, which function as net sinks for atmospheric CH 4 [19]. These preliminary data highlighted the presence of hotspots of CH 4 emission, most notably during the wet season, associated with anoxic conditions in localised topographic low points. Here we report data from longer-term measurements from January 2011 through June 2013 and seasonal intensive measurement campaigns from the 12th through 22nd of November 2011 and 12th through 21st of August 2012. These longer-term observations are used to (1) improve our assessment of the magnitude and seasonality of CH 4 emissions previously indicated [19] through analysis of a further 18 months of observations. To better understand the occurrence of emission hotspots, the intensive seasonal campaigns observations are used to (2) characterise the spatial controls on soil-atmosphere CH 4 exchange across and within landscape features in terms of edaphic and environmental conditions and (3) investigate the relationship between gross methanotrophic and methanogenic process rates in these soils. We hypothesise that: (H1) variations in net CH 4 flux will be best explained by soil O 2 concentration, rather than water-table depth or WFPS; and (H2) variations in methanogenic activity will be driven by C availability.

Study Sites
The study was carried out in the vicinity of Tres Cruces (13 • 07 19" S, 71 • 36 54" W), on the western border of Manu National Park in the Southeastern Peruvian department of Cusco (Figure 1a). This area falls within the watersheds dividing western xeric puna highlands and eastern steeply sloped evergreen tropical montane cloud forests that extend towards the Amazonian basin. The study area lies between the summits of ridges at approximately 4000 m asl and the treeline at approximately 3400 m asl, and has previously been described as humid puna grassland [29,65]. Such grasslands cover approximately 210 km 2 along the western border of Manu National Park or on the order of 20% of the land-cover above 600 m asl [65]. Total precipitation is 1900 to 2500 mm yr −1 and mean annual air temperature, at 3600 m asl, is 11 • C [65]. Precipitation is intensely episodic and there is a pronounced wet season between November and April ( Figure S1). In contrast, diurnal differences in air temperature are greater than seasonal variations ( Figure S1). This is particularly notable during the dry season, when cloud free nights commonly result in morning frosts. in net CH4 flux will be best explained by soil O2 concentration, rather than water-table depth or WFPS; and (H2) variations in methanogenic activity will be driven by C availability.

Study Sites
The study was carried out in the vicinity of Tres Cruces (13°07′19′′ S, 71°36′54′′ W), on the western border of Manu National Park in the Southeastern Peruvian department of Cusco ( Figure  1a). This area falls within the watersheds dividing western xeric puna highlands and eastern steeply sloped evergreen tropical montane cloud forests that extend towards the Amazonian basin. The study area lies between the summits of ridges at approximately 4000 m asl and the treeline at approximately 3400 m asl, and has previously been described as humid puna grassland [29,65]. Such grasslands cover approximately 210 km 2 along the western border of Manu National Park or on the order of 20% of the land-cover above 600 m asl [65]. Total precipitation is 1900 to 2500 mm yr −1 and mean annual air temperature, at 3600 m asl, is 11 °C [65]. Precipitation is intensely episodic and there is a pronounced wet season between November and April ( Figure S1). In contrast, diurnal differences in air temperature are greater than seasonal variations ( Figure S1). This is particularly notable during the dry season, when cloud free nights commonly result in morning frosts. The ecosystem at this site is characterised by upland summits and backslopes that are dominated by tussock grasses and transition into toeslopes containing wetland environments of up to 1 ha in size ( Figure S2). These wetlands consist of peat-forming depressions, moss-filled hollows and shallow lakes of varying degrees of permanency ( Figure S3). The study site has a history of cattle-grazing by local communities, and the landscape is susceptible to burning during dry periods [65,66]. Typical of paramo and wetter puna ecosystems, above-ground biomass is dominated by tussocks of Calamagrostis sp. and lower abundances of other grasses such as Scirpus sp., Festuca sp., and Juncus sp., in addition to mosses in moist locations and diverse herbs, shrubs and ferns [26,65,67]. The upland soils are 20 to 40 cm in depth, and consist of a thick organic rich mineral horizon overlying a thinner, stony layer. The superficial soils are acidic and typically have bulk densities on the order of 0.40 g cm −2 , C contents of approximately 15% C and C to nitrogen (C:N) The ecosystem at this site is characterised by upland summits and backslopes that are dominated by tussock grasses and transition into toeslopes containing wetland environments of up to 1 ha in size ( Figure S2). These wetlands consist of peat-forming depressions, moss-filled hollows and shallow lakes of varying degrees of permanency ( Figure S3). The study site has a history of cattle-grazing by local communities, and the landscape is susceptible to burning during dry periods [65,66]. Typical of paramo and wetter puna ecosystems, above-ground biomass is dominated by tussocks of Calamagrostis sp. and lower abundances of other grasses such as Scirpus sp., Festuca sp., and Juncus sp., in addition to mosses in moist locations and diverse herbs, shrubs and ferns [26,65,67]. The upland soils are 20 to 40 cm in depth, and consist of a thick organic rich mineral horizon overlying a thinner, stony layer. The superficial soils are acidic and typically have bulk densities on the order of 0.40 g cm −2 , C contents of approximately 15% C and C to nitrogen (C:N) ratios of approximately 14 [29]. Wetland peat soils range from 40 to over 100 cm in depth and have low bulk densities and C contents in excess of Soil Syst. 2019, 3, 2 5 of 21 40% C [29]. These peats are well humidified in toeslope settings where a mixture of moss and grass species occur, but may develop extensive accumulations of moss litter in persistently wet hollows. The composition of the stony sub-soil in this environment reflects the Palaeozoic shale-slate geology of the region [68].
Soil-atmosphere CH 4 exchange for 2011 has previously been reported for this site as part of a study investigating non-CO 2 trace gas fluxes along an Andean altitude transect [19]. These measurements indicated that the site acts as a considerable source of CH 4 to the atmosphere when compared to uptake by the surrounding montane forests [19,69]. Emissions from the site were highly seasonal and largely driven by hotspots at the transition between the foot-and toeslope, whilst the backslope and summit varied between source and sink activity. Greater fluxes from these soils than those of the surrounding forests was attributed to near-saturated soils and subsequent lower soil O 2 concentrations.

Sampling Approach
As previously described, this study focused on a northeast-southwest trending summit to toeslope transition ( Figure 1b) covering approximately four hectares at 3650 m asl [19,70]. The transition from a narrow ridge summit, through a backslope of approximately 15 • , to the footslope occurs over 200 m and a 50 m decrease in elevation. A relatively abrupt transition leads to a broad toeslope that extends southwest for a further 100 m with little change in elevation before terminating in an escarpment where slope angles increase once again and eventually lead into the forests occupying the steeply inclined valleys that drain this landscape. To investigate the controls on soil-atmosphere CH 4 exchange from this environment, measurements were taken at two temporal and spatial scales. In both cases sampling equipment was installed a minimum of four weeks prior to measurements to minimise the influence of disturbances. The large and dense nature of the grass tussocks, reaching heights and diameters of up to 60 and 15 cm respectively, which dominate the landscape precluded their direct inclusion in the sampling. For this reason the 20 cm diameter soil collars used in gas exchange measurements, which were inserted to 5 cm depth in order to minimise long-term root damage, typically encompassed bare soil, small herbs, individual grass stalks and mosses, but not mature examples of the principal tussock-forming grasses. As such the data collected is biased towards the soil surface and we are unable to assess the role of the dominant plants in transporting gas between the soil and atmosphere. Sampling order was rotated and sampling was conducted between the hours of 08:00 and 16:00 to minimise effects of temporal variability associated with sunrise and sunset at approximately 06:00 and 18:00 [71].
To investigate seasonal variability in soil-atmosphere CH 4 exchange, four plots were established for longer-term measurements on dominant landscape features encompassing the upper and lower reaches of the study area; specifically, a summit plot on the ridge, a backslope plot near the ridge, a wet footslope plot and a dry footslope plot both at the transition between the foot-and toeslope. Each plot was instrumented at five sampling stations approximately 5 to 10 m apart. In addition to soil collars inserted at each sampling station, soil gas equilibration chambers were buried at three stations per plot. At each sampling station, soil-atmosphere gas exchange, soil moisture, soil temperature, and, where soil-gas equilibration chambers were present, soil O 2 concentration were measured monthly. Measurements at the summit, backslope, and dry footslope plots were carried out, inclusively, from January 2011 to June 2013. The wet footslope plot was established at a later date than these plots with measurements beginning in August 2011. Measurements were not possible in July and December 2012 and February 2013 due to access restrictions.
To investigate spatial variability in soil-atmosphere CH 4 exchange and soil CH 4 cycling, intensive measurement campaigns were carried out in both the wet and dry season. A stratified sampling grid consisting of 24 sampling stations encompassing the topographic transition was established along six 75 m transects running perpendicular to the principle northeast-southwest trend in slope. Each transect consisted of four sampling stations with a footprint of 0.5 m 2 positioned 25 m apart. The transects had a down-slope separation of 50 m. Prior to each campaign, running from the 12th to 22nd of November 2011 in the wet and 12th to 21st of August 2012 in the dry season, each sampling station was equipped with a soil collar, soil gas equilibration chamber, and a piezometer. Sampling stations were visited daily to measure soil-atmosphere gas exchange, soil moisture, soil temperature, water-table depth, and soil O 2 concentration. On the final day of each campaign soil CH 4 concentration was also determined. The equipment was removed following the wet season campaign and replaced directly adjacent to its previous position on undisturbed ground prior to the dry season campaign. Following the wet season campaign, the footprints of each soil collar were sampled to characterise soil properties and provide material to investigate gross process rates. Following the dry season campaign a vegetation survey was conducted at each sampling station.

Soil-Atmosphere Gas Exchange and Environmental Conditions
Soil-atmosphere exchange of CH 4 and CO 2 and environmental conditions were determined following standard procedures as previously described in detail by Jones et al. [69]. Briefly, diffusive gas fluxes were determined based on the temporal change in concentration of CH 4 and CO 2 across four discrete 20 mL gas samples (12 mL Exetainer, Labco Ltd., Lampeter, UK) taken from the headspace of a static chamber over a period of approximately 35 min. Each chamber consisted of a cylindrical cap equipped with a gas sampling port, pressure equilibration port, thermocouple port and a small mixing fan. To initiate a measurement the cap would be carefully mounted, using a rubber sleeve, on the pre-installed soil collar at each sampling station to enclose a volume of 0.008 m 3 over a soil surface area of 0.031 m 2 . The cap was demounted following each measurement. Care was taken when mounting and demounting the caps to minimise disturbance of the soil in the vicinity of the sampling equipment and sampling ports were extended with 2 m lengths of tubing to avoid disturbance during gas sampling over the course of each measurement. Temperature (type k thermocouple, Omega Engineering Ltd., Manchester, UK) within the chamber and ambient air temperature at 5 cm above the soil surface and atmospheric pressure (Garmin GPSmap 60CSx, Garmin Ltd., Olathe, Kansas, USA) were recorded at each gas sampling interval. Concentrations of CH 4 and CO 2 were determined by gas chromatography (Thermo TRACE GC Ultra, Thermo Fisher Scientific Inc., Waltham, Massachusetts, USA). Diffusive fluxes were calculated from concentration time-series, after assessment for evidence of sampling and analysis artefacts such as ebullition events or sample leakage during storage and transport, based on linear and non-linear relationships with time. Fluxes were converted from a concentration to amount basis following the Ideal Gas Law and are reported in mg CH 4 -C m −2 d −1 and g CO 2 -C m −2 d −1 . Positive and negative fluxes respectively indicate emission and uptake by the soil surface. Soil environmental condition measurements coincided with each gas flux measurement. As described in Jones et al. [69], soil gas equilibration chambers, constructed from gas-permeable silicone rubber tubing (AP202/60, Advanced Polymers Ltd., Worthing, UK) and sampling tubes fitted with stop-cocks, were buried at 10 cm below the soil surface. Each chamber had an internal volume of 50 cm 3 and a surface area of 57 cm 2 . Soil O 2 concentration in each chamber was determined by using a pair of syringes and stop-cocks to remove 40 mL gas from a chamber, pass it through the flow-through head of an oxygen sensor (MO-200, Apogee Instruments Inc., Logan, Utah, USA) and then return it to the chamber. On the final day of each campaign, the gas was stored (12 mL Exetainer, Labco Ltd., Lampeter, UK) for determination of soil CH 4 concentration by gas chromatography rather than being returned to the chamber. Soil water content and soil temperature (type k penetration probe, Omega Engineering Ltd., UK) were measured in triplicate at each sampling station. During the longer-term measurements soil water content was determined in the upper 6 cm (ML2x ThetaProbe, Delta-T Ltd., Burwell, UK) and, as reliable plot estimates of porosity are unavailable, reported as volumetric water content (VWC). Soil temperature was measured at 5 cm depth. During the intensive seasonal campaigns soil water content was determined in the upper 20 cm (CS620 Hydrosense, Campbell Scientific Inc., Logan, Utah, USA) and, calculated from VWC and sampling station estimates of porosity, reported as WFPS. Soil temperature was measured at 5 cm and 10 cm depth. During the intensive campaigns, water-table depth from the surface was measured at each sampling station in Soil Syst. 2019, 3, 2 7 of 21 piezometers constructed from 35-mm diameter plastic tubing and installed to a maximum depth of 60 cm. Water-table depth is reported to a maximum depth of 20 cm below surface to account for variations in total soil depth across the study area.

Site Characterisation
The soils and vegetation of the intensive campaign sampling stations were characterised following the wet and dry season campaigns. Following the intensive wet season campaign, soil depth was measured and the soil sampled at each sampling station. Paired soil samples were taken from 0-5 and 5-15 cm within the soil collar footprint. To sample the surface soils between 0-5 cm, vegetation was clipped to the soil surface and samples consisting of a block 10 × 5 × 5 cm, with a volume of 250 cm 3 , removed using scissors. Soils from 5-15 cm depth were sampled using a 10 cm long, 5.0 cm diameter corer to yield a sample with a volume of 173 cm 3 . Soil depth was determined by inserting a metal rod until physical resistance to further insertion was encountered [29]. One sample from each depth pair was processed by gently homogenizing and removing root fragments for use in the incubation experiments and determination of soil chemical properties. Soil pH was determined (HANNA pHep 4, HANNA Instruments, Woonsocket, Rhode Island, USA) in a gravimetric 1:2 slurry of air dried soil and deionised water [72]. Soil C and N content were determined by elemental analysis after grinding to a fine powder (Costech ECS4010, Costech Analytical Technology Inc., Valencia, California, USA and Finnigan Deltaplus XP GC-IRMS, Thermo Fisher Scientific Inc., Waltham, Massachusetts, USA). The second sample from each depth pair was kept intact to determine bulk density after drying for 24 h at 105 • C. Bulk density was calculated as the oven dry mass per soil field volume of this sample. Subsequently, particle density was determined using 10 mL pyncometers [73]. Porosity was estimated from these data based on bulk and particle density estimates. A vegetation survey was carried out following the intensive dry season campaign. At each of the sampling locations two 0.25 m 2 quadrats were placed, perpendicular to the slope, at either side of each sampling station. Additionally, slope angle was measured in the principle downslope, northeast-southwest, orientation, and perpendicularly across slope using a clinometer and ranging poles. Within each quadrat every non-bryophyte was identified to the genus or species level. The biomass density, reported in g m −2 , associated with the dominant grass genera, Calamagrostis sp., Scirpus sp. and Juncus sp., was determined from measurements of tussock height, and basal and crown diameter, following the optimal allometric model previously developed for these grasslands [67].

Laboratory Incubations
Gross rates of CH 4 production and consumption were determined by incubation of the soils sampled from 0-5 and 5-15 cm depth at each of the sample stations at the end of the wet season intensive campaign using the trace gas pool dilution approach developed by Von Fischer and Hedin [57] and applied as described by Yang et al. [74], Yang and Silver [64] and Yang et al. [59]. Briefly, approximately 100 g of homogenised soil, at field water content, from each sampling station and depth were placed in 1 L Kilner jars and loosely sealed with screw-cap lids fitted with septa. After 48 h of pre-incubation, a jar was vented and fully sealed before being spiked with 1 mL of nitrogen carrier containing 0.3 ppm sulphur hexafluoride and 70 ppm 13 CH 4 . In doing so the initial total CH 4 concentration in the jar headspace, of approximately 2 ppm, was slightly elevated but close to that of the atmosphere. The jar headspace was mixed with a 60 mL syringe and pre-incubated for 30 to 60 min to allow these tracers to fully equilibrate within the headspace. Following the pre-incubation period, 100 mL of nitrogen was injected and the headspace mixed at four discrete times over the course of 10 to 12 h. The resulting 100 mL overpressure at each time-step was sampled and stored in evacuated 60 mL Wheaton bottles sealed with 20 mm butyl septa (Geo-microbial Technologies Inc., Ochelata, Oklahoma, USA). Sub-samples of 5 mL were taken from each bottle to determine the concentration of CH 4 , CO 2 , and sulfur hexafluoride by gas chromatography (Thermo TRACE GC Ultra, Thermo Fisher Scientific Inc., Waltham, Massachusetts, USA) and the remainder was used to determine the carbon isotope composition of the CH 4 by isotope ratio mass spectrometry (PreCon and Finnigan Deltaplus XP GC-IRMS, Thermo Fisher Scientific Inc., Waltham, Massachusetts, USA). Concentration measurements were corrected for the effect of sampling dilution based on the change in sulphur hexafluoride concentration and net fluxes of CH 4 and CO 2 calculated based on the temporal change in their concentrations. The change in the concentration of 13 CH 4 over time was used to estimate the rate of oxidation and the gross rate of CH 4 production was determined iteratively [59,64] by assuming the fractionation factor associated with oxidation to be 0.98 [57] and the C isotope composition of CH 4 produced, based on pilot anoxic incubations of these soils, to be −75‰ VPDB. The gross rate of consumption was calculated as the difference between the net CH 4 flux and the gross rate of production [59,64] and the total and methanogenic fraction of carbon mineralisation was calculated following Von Fischer and Hedin [61]. These incubations took place in the dark at 24 • C.

Data Processing and Statistics
Data processing and statistical analyses were conducted in R [75]. A generalised least squares approach was used to reduce the influence of heteroscedasticity in investigating the aggregate effect of season on monthly mean net CH 4 flux, net CO 2 flux, VWC, O 2 concentration, and soil temperature within each longer-term measurement plot [76,77]. The drivers of temporal variations in net CH 4 flux within the longer-term measurement plots were investigated after averaging monthly observations by plot. The drivers of spatial variations in net CH 4 flux among intensive campaign sampling stations were investigate after averaging each campaign by sampling station. Temporal and spatial relationships between field measurements of gas fluxes and environmental variables produced residuals with a high degree of non-normality when assessed with parametric approaches. For this reason, temporal relationships between monthly plot means for the longer-term measurements and spatial relationships between sampling location campaign means for the intensive campaign measurements were investigated using Spearman's rank correlation [78]. Relationships between gross process rates of CH 4 cycling and C mineralisation determined in the laboratory incubations were investigated by linear regression [78]. Significance is reported as p < 0.05 and summarised data are reported as means and, in parentheses, standard deviations, unless stated otherwise.

Longer-Term Measurements
Net CH 4 fluxes measured monthly at the summit, backslope and wet footslope longer-term plots were significantly greater during wet than dry season months, with all plots acting as net wet season sources to the atmosphere (Table 1). This behaviour was under-pinned by considerable within-plot variability with individual monthly sampling locations varying between sink and source activity throughout the measurement period ( Figure 2). This was particularly true of the wet footslope plot where CH 4 emissions had the potential to be one to two orders of magnitude greater than the fluxes measured from the summit, backslope or dry footslope plots. Net CO 2 fluxes showed little apparent difference between wet and dry season months and tended to be slightly greater on the footslope than summit or backslope plots. Soil O 2 concentrations were somewhat, but not significantly, lower during wet season months and indicate, mostly notably at the wet footslope, the presence of transient anoxic conditions. Soil temperatures were similar among plots with significantly warmer temperatures during wet season months, whilst, little variability in VWC was observed. The presence of different superscript letters a, b within plots indicates significant differences between wet and dry season for a given variable (generalised least squares, p < 005) Soil Syst. 2019, 3, x FOR PEER REVIEW 9 of 22 The presence of different superscript letters a, b within plots indicates significant differences between wet and dry season for a given variable (generalised least squares, p < 005) Temporal variations in mean monthly net CH4 flux were significantly positively correlated with soil temperature and consistently, but not significantly, negatively correlated with O2 concentration at all the longer-term measurement plots (Table 2). Clear patterns were not found with the other measured variables, for example, net CH4 flux was also significantly correlated with Temporal variations in mean monthly net CH 4 flux were significantly positively correlated with soil temperature and consistently, but not significantly, negatively correlated with O 2 concentration at all the longer-term measurement plots (Table 2). Clear patterns were not found with the other measured variables, for example, net CH 4 flux was also significantly correlated with net CO 2 flux at the summit and dry footslope plots but the direction of these relationships were opposing.

Site Characteristics
The transition from the top of the backslope to the toeslope dominated the environmental gradients captured by the intensive sampling campaigns (Table 3). Based on the characteristics of the topography, vegetation and soil, the 24 sampling stations were assigned to one of four morphological groups (Figure 1b). Higher elevation sampling stations on the steep ground of the backslope were defined as upper slope (n = 8) or lower slope (n = 9) sampling stations. Similarly sampling stations on the gently sloping or flat ground of foot-and toeslope were defined as depression (n = 3) or hollow (n = 4) sampling stations. Of twelve non-bryophyte genera identified, the grasses Calamagrostis sp., Scirpus sp., and Juncus sp. were, by far, the most abundant. The total abundance of these grasses generally decreased downhill from the upper slope to hollow sampling stations. The upper slope and hollow sampling stations were typified by inverse patterns in the proportion of total above-ground grass biomass contributed by Calamagrostis sp. and Juncus sp. In contrast, the lower slope and depression sampling stations were not clearly distinguishable with respect to grass biomass. Whilst not quantified, mosses were notably present in all the sampling stations on the foot-and toeslope but were rare on the slopes. These differences were also reflected by soil depth, with shallower soils on the backslope and deeper soils on the foot-and toeslope. The shallower backslope soils were typically organomineral in origin, whilst peats formed in the foot-and toeslope (Table 4). These organomineral soils had higher bulk and particle densities, and lower C contents and C to nitrogen ratios (C:N) than the peats. In the foot-and toeslope, peat soils were differentiated from each other by the presence of well-humidified peats in the depressions, whilst the hollows principally consisted of mosses and their poorly decomposed litter. This distinction is reflected by higher C:N in the hollows than depressions. Soils across this landscape were acidic.

Soil-Atmosphere Gas Exchange and Environmental Conditions
Uptake and emission of CH 4 by the soils of the intensive campaign sampling stations was observed in both the wet and dry season. Reflecting the temporal patterns seen in the longer-term measurements (Table 1), weaker sink and stronger source activity were found in the wet season than dry season with daily net CH 4 fluxes ranging from −0.8 to 291.9 mg CH 4 -C m −2 d −1 and −1.1 to 39.7 mg CH 4 -C m −2 d −1 , respectively. On average in both the wet (Figure 3a) and dry season (Figure 4a) campaigns hollow sampling stations acted as sources of CH 4 to the atmosphere, whilst the backslope and depression sampling stations acted as both sources and sinks. Soil-atmosphere exchange at the upper slope and depression sampling stations was relatively similar between campaigns with the stronger source activity in the wet season resulting from greater emissions from the lower slopes and hollows. Similarly, soil respiration was slightly greater during the wet season than dry season campaign with daily net CO 2 fluxes ranging from 0.2 to 6.0 g CO 2 -C m −2 d −1 and 0.4 to 5.1 g CO 2 -C m −2 d −1 . In both campaigns this exchange was greater on the foot-and toeslope than on the backslope with greatest emissions from the depression sampling stations and lowest emissions from the upper slopes. Underpinning these gas exchanges and capturing the temporal variability of the longer-term measurements (Table 1), air and soil temperatures were greater in the wet than the dry season. During the wet season, daily air temperature at the surface, and soil temperature at 5 and 10 cm depth ranged from 7.4 to 27.6 • C, 4.2 to 15.4 • C, and 4.2 to 13.8 • C, respectively. During the dry season, daily air temperature at the surface and soil temperature at 5 and 10 cm depth for the dataset ranged from 5.2 to 24.1 • C, 4 to 13.2 • C, and 4.6 to 11.3 • C, respectively. However, on average there was little spatial variability with respect to temperature among the sampling stations within each campaign. Soil O 2 concentrations ranged from anoxic to atmospheric and as with the longer-term measurements (Table 1) were lower in the wet (Figure 3b) than dry (Figure 4b) season campaign. On average, evidence of widespread anoxia was only apparent for the hollow and some of the lower slope sampling stations during the wet season campaign (Figure 3b). On the final day of both the wet and dry season campaigns soil CH 4 concentrations were greater than that of the atmosphere and ranged from 2.2 ppm to 6.1% and 2.4 ppm to 0.1%, respectively. Water-table depth below the surface was greater during the dry than wet season campaign with daily measurements ranging from 1.4 to 20.0 cm and −5.1 to 20 cm. On average and, as with soil O 2 concentration, evidence of inundation near the surface was only apparent for the hollow and some of the lower slope sampling stations during the wet season (Figure 3c) but not dry season campaign (Figure 4c). As with the longer-term measurements of soil water content (Table 1), WFPS was relatively similar between the wet and dry season campaigns with daily measurements ranging from 53.1 to 100.0% and 58.7 to 100.0%, respectively. On average, WFPS was greatest at the hollow sampling stations in both the wet (Figure 3d) and dry season (Figure 4d) campaigns. Sampling station means, as plotted in Figures 3 and 4, and standard deviations by morphological group and season associated with these patterns are reported in Table S2 of the supplementary material. deviations by morphological group and season associated with these patterns are reported in Table  S2 of the supplementary material.

Relationships between Soil-Atmosphere Gas Exchanges and Environmental Conditions
Clear temporal relationships between the measured variables were not found over the short period of the wet and dry season campaigns. Spatially, mean net CH 4 flux during the wet season was strongly negatively correlated with O 2 concentration and more weakly correlated with WFPS, and water-table depth below the surface (Table 5). In the dry season, the correlation between net CH 4 flux and these variables was similar, with a slightly weaker correlation with O 2 concentration and slightly stronger correlations with WFPS and water-table depth. In both campaigns relatively strong co-correlations were also found amongst O 2 concentration, water-table depth and WFPS. Excluding sampling stations with mean net fluxes greater than 0.5 mg CH 4 -C m −2 d −1 , the observed relationship between net CH 4

Relationships between Soil-Atmosphere Gas Exchanges and Environmental Conditions
Clear temporal relationships between the measured variables were not found over the short period of the wet and dry season campaigns. Spatially, mean net CH4 flux during the wet season was strongly negatively correlated with O2 concentration and more weakly correlated with WFPS, and water-table depth below the surface (Table 5). In the dry season, the correlation between net CH4 flux and these variables was similar, with a slightly weaker correlation with O2 concentration and slightly stronger correlations with WFPS and water-table depth. In both campaigns relatively strong co-correlations were also found amongst O2 concentration, water-table depth and WFPS. Excluding sampling stations with mean net fluxes greater than 0.5 mg CH4-C m −2 d −1 , the observed relationship between net CH4 flux and O2 concentration held for both the wet (n = 12, ρ = −0.67) and dry (n = 19, ρ = −0.74) season campaigns, whilst relationships with water-table depth and WFPS were insignificant.  Table 5. Spearman's rank correlation coefficients for spatial relationships between sampling station mean net CH 4 flux and net CO 2 flux, soil O 2 concentration, water-table depth, and WFPS for the wet (n = 24) and dry (n = 24) season intensive campaigns. Following the coefficients, * indicates p < 0.05 and ** indicates p < 0.01.

Campaign
Net

Gross Rates of Production and Consumption
Incubations of soil from 0-5 cm depth acted as both net sources and sinks for CH 4 with a mean flux of 10.9 (27.01) ng CH 4 -C g dry soil −1 h −1 ( Figure 5). In contrast, incubations of soil from 5-15 cm principally acted as net sinks with a mean flux of −1.1 (1.39) ng CH 4 -C g dry soil −1 h −1 . Significant production of 13 CH 4 in 11 out of 24 of the incubations of soils from 0-5 cm precluded reliable estimation of gross rates of CH 4 production and consumption, whilst, a strong positive linear correlation (p < 0.01, r 2 = 0.98) between observed and modelled CH 4 concentrations was apparent for all remaining observations. Excluding those that invalidated the model assumptions, the mean net CH 4 flux for the remaining incubations of soil from 0-5 cm depth was −0.6 (2.87) ng CH 4 -C g dry soil −1 h −1 . Net CO 2 fluxes associated with these incubations and those of soils from 5-15 cm were respectively 22.5 (20.16) and 5.3 (6.6) µg CO 2 -C g dry soil −1 h −1 . Mean estimated rates of gross production and consumption were respectively 4.8 (4.32) and 5.3 (5.8) ng CH 4 -C g dry soil −1 h −1 for soils from 0-5 cm depth and 1.7 (1.48) and 2.8 (2.74) ng CH 4 -C g dry soil −1 h −1 for soils from 5-15 cm depth. Across incubations from both depths, these gross rates of production and consumption were positively linearly correlated (p < 0.01, r 2 = 0.79). Similarly, gross production rate was strongly correlated with the total rate of carbon mineralisation (p < 0.01, r 2 = 0.65) but not the methanogenic fraction of carbon mineralisation (p > 0.1, r 2 = 0.02). flux of 10.9 (27.01) ng CH4-C g dry soil h ( Figure 5). In contrast, incubations of soil from 5-15 cm principally acted as net sinks with a mean flux of −1.1 (1.39) ng CH4-C g dry soil −1 h −1 . Significant production of 13 CH4 in 11 out of 24 of the incubations of soils from 0-5 cm precluded reliable estimation of gross rates of CH4 production and consumption, whilst, a strong positive linear correlation (p < 0.01, r 2 = 0.98) between observed and modelled CH4 concentrations was apparent for all remaining observations. Excluding those that invalidated the model assumptions, the mean net CH4 flux for the remaining incubations of soil from 0-5 cm depth was −0.6 (2.87) ng CH4-C g dry soil −1 h −1 . Net CO2 fluxes associated with these incubations and those of soils from 5-15 cm were respectively 22.5 (20.16) and 5.3 (6.6) μg CO2-C g dry soil −1 h −1 . Mean estimated rates of gross production and consumption were respectively 4.8 (4.32) and 5.3 (5.8) ng CH4-C g dry soil −1 h −1 for soils from 0-5 cm depth and 1.7 (1.48) and 2.8 (2.74) ng CH4-C g dry soil −1 h −1 for soils from 5-15 cm depth. Across incubations from both depths, these gross rates of production and consumption were positively linearly correlated (p < 0.01, r 2 = 0.79). Similarly, gross production rate was strongly correlated with the total rate of carbon mineralisation (p < 0.01, r 2 = 0.65) but not the methanogenic fraction of carbon mineralisation (p > 0.1, r 2 = 0.02).

Figure 5.
Observed rates of net CH 4 flux and modelled rates of gross CH 4 production and gross CH 4 consumption by incubation depth. Box lower hinge, middle, and upper hinge respectively indicate the 25% quartile, median, and 75% quartile, and points indicate the rates associated with individual incubations. Incubations excluded from determination of gross production and consumption rates, because significant 13 CH 4 production invalidated the modelling approach, are indicated in red. Subsequently, no modelled production or consumption rates are plotted from these observations. Incubations retained for determination of gross process rates are indicated in blue.

Discussion
Extension of the longer-term monthly measurements of net CH 4 flux from a period of one to two and half years corroborated the previous indication [19]; i.e., that the soils in this study site are a net source of atmospheric CH 4 that peaks during the wet season (Table 1 and Figure 2) in response to increases in temperature and decreases in soil O 2 concentration ( Table 2). The fact that temporal variations in monthly mean net CH 4 fluxes within the longer-term plots were best explained by mean monthly soil temperatures presumably reflects the fact that coincidence of warmer and wetter conditions ( Figure S1) influences factors, such as microbial metabolic activity and plant phenology, with positive influences on the availability of the niches and substrates that support CH 4 production. Weaker negative correlations between these fluxes and mean monthly soil O 2 concentrations may result from two factors. First, temporal variations in the soil conditions, such as the supply of substrates from plants, governing the CH 4 cycle which co-correlate with soil temperature are likely complex and beyond the scope of the data presented. Second, as soil O 2 concentrations were only measured at three of the five sampling stations in each longer-term plot, the apparently important role of hotspots in these soils may not have been fully captured by these measurements. Along with the observation by Veber et al. [20], that peatlands in the Colombian paramo also act as sources of atmospheric CH 4 with median rates of emission up to 55 mg CH 4 -C m −2 d −1 , these data suggest that wet and humid paramo and puna ecosystems function as regional hotspots for CH 4 emission when compared to sink activity reported for other Andean and Western Amazonian upland environments. For example, studies in Ecuador [79] and Peru [69] have demonstrated that montane and premontane forests on the eastern flank of the Andes are net sinks for atmospheric CH 4 on the order of −0.2 to −1.5 mg CH 4 -C m −2 d −1 . Similarly, Palm et al. [80] illustrate that upland soils under a variety of land-uses in the Peruvian Amazon, with average ecosystem exchange rates between −0.7 and 0.4 mg CH 4 -C m −2 d −1 , principally act as a CH 4 sink with source activity only observed in sites of high-input cropping agriculture. Indeed, the need to include these montane environments in regional budgetary considerations is highlighted by the fact that emissions across the landscape considered here are comparable to those for wet upland soils, reaching up to 18 mg CH 4 -C m −2 d −1 as reviewed by Spahni et al. [13], and that wet season emissions from the longer-term wet footslope plot were of the same order of magnitude as diffusive fluxes, with means in the range of 110 to 470 mg CH 4 -C m −2 d −1 , reported for flooded forests and palm swamps in the Western Amazon [81,82]. To assess the contribution made by the soils of humid puna and paramo ecosystems, respectively covering approximately 230,000 and 41,000 km 2 of the Andes [22], at the regional scale requires a better understanding of how representative limited examples of source behaviour such as that reported here are across their wider extent.
Considerable spatial differences ( Figure 2) in soil-atmosphere CH 4 exchange are apparent both between and within measurement plots for these longer-term data suggesting that up-scaling the activity of such landscapes requires appropriate land-cover estimates at meso-and microtopographical scales [18,83,84]. For example, the strong source behaviour of the wet footslope was driven by emission hotspots at some of the sampling locations within this plot. Indeed, instances of strong emissions, weak emissions, and uptake of CH 4 during the same sampling period within a plot are not uncommon, suggesting that the constraints on soil-atmosphere CH 4 exchange vary considerably at scales of tens of metres. It is reasonable to presume that these variations result from the influence of micro-topography on the flow of water, substrates and nutrients through this landscape [31]. However, the fact that source activity is less prevalent in the dry footslope, which we may expect to be more prone to inundation, than on the summit or backslope suggests that traditional models used to explain soil-atmosphere CH 4 exchange may not be appropriate across this landscape, and that inclusion of upland soils in up-scaling efforts is likely critical due to their greater areal coverage.
The drivers of the spatial variability observed in the longer-term measurements are elucidated by the more detailed measurements made during the intensive seasonal campaigns. The spatial and seasonal patterns of soil-atmosphere CH 4 exchange captured by these campaigns (Figures 3 and 4a) are broadly consistent with those shown by the longer-term measurements (Figure 2). It is clear from the site characterisations that hotspots of emission and differences in net exchange between the longer-term wet and dry footslope plots result from the occurrence of two different peat forming morphological landforms on the foot-and toeslope of this area. Namely, moss-accumulating hollows that are prone to inundation and act as strong emitters of CH 4 , or better-drained depressions that tend to consume atmospheric CH 4 . Furthermore, characterisation of the conditions and behaviour of the upper and lower portions of the backslope indicate that emissions from longer-term summit and backslope plots originate from upland organo-mineral soils (Table 4) and that the potential for CH 4 emissions from these soils is greater than might be expected based on topographic position alone (Table 3).
Transitions between CH 4 uptake and emission across the landscape, even after exclusion of emission hotspots associated with inundation, in both the wet and dry season intensive campaigns were best explained by the availability of O 2 (Table 5). These spatial relationships indicate the role of decreases from atmospheric concentrations of O 2 in promoting the development of anoxic zones suitable for methanogenesis and, subsequently, CH 4 production even when the bulk soil remained relatively oxic (Figures 3 and 4b). This finding is in line with our first hypothesis (H1) that bulk soil O 2 concentration would better reflect constraints on methanogenesis across this landscape than water-table depth or WFPS. In this respect, these patterns ultimately reflect the influence of variations in soil water and aerobic biological activity on the transport and depletion of soil O 2 [52]. The fact that emissions were observed even in oxic soils supports evidence for the ubiquity of simultaneous methanogenic and methanotrophic activity in heterogeneous soil environments [61]. Indeed, the presence of above atmospheric soil CH 4 concentrations at the end of both campaigns, irrespective of a sampling locations role as a sink or source of CH 4 , highlights the combined role of methanogens and high-capacity-low-affinity and low-capacity-high affinity methanotrophs in the cycling of CH 4 in these soils.
Strong net CH 4 emissions in oxic incubations of soil from 0-5 cm but not 5-15 cm depth ( Figure 5) indicate that, as in other wet organic rich soils [12] in the tropics, the direction of field soil-atmosphere exchange is principally driven by variations in the activity of CH 4 cycling microbial communities close to the surface. Such a situation supports the inference that field net CH 4 flux observations are driven by the influence of soil O 2 concentration measured at 10 cm depth through controls on the balance between methanogenesis and methanotrophic processes in the superficial oxic soils rather than methanogenic activity in deeper, potentially anoxic horizons. Whilst we were unable to deconvolve gross process rates for incubations exhibiting the highest rates of emission, positive correlations between gross production and consumption support the indication that methanotrophic communities typically have the capacity to track increases in CH 4 production [59,61,64]. Indeed, coincidence of positive soil-atmosphere CH 4 gradients and net uptake in the field indicates the potential of high-capacity-low-affinity methanotrophic communities to attenuate emissions in these soils is high [12,45]. Contrary to our second hypothesis (H2), we found that the increase in total C mineralisation, rather than the methanogenic fraction of C mineralisation, best explained increases in gross CH 4 production in our incubations. Such relationships between the methanogenic fraction of C mineralisation and gross production rates have been interpreted as reflecting the availability of C to methanogens within highly reduced microsites absent of more energetically favourable energy sources, such as nitrate, ferric iron, and sulphate [59,61,64]. The alternative scenario observed here suggests that substrate competitions are not a major constraint on methanogenesis across the soils of this landscape. This could be because high labile soil C availability and relatively stable soil moisture conditions in the superficial organic soils lead to the rapid and widespread reduction of terminal electron acceptors such as ferric iron that might otherwise suppress methanogenesis [61,85,86].

Conclusions
Here we report soil-atmosphere CH 4 exchanges from a humid puna ecosystem in the Southeastern Peruvian Andes, and highlight the need to consider the soils of these environments as sources of CH 4 in regional atmospheric budgets. We show that CH 4 emissions occur not only from peat soils on relatively flat ground but also from organo-mineral soils on ridges and steep slopes. Seasonal variations in net CH 4 fluxes reflect the influence temperature and rainfall on below-ground O 2 availability, and the balance between methanogenic and methanotrophic activity. Production of CH 4 in superficial soils is widespread, and spatial variability in net CH 4 fluxes is principally constrained by variations in soil O 2 concentration. Due to the size of the tussocks formed by the dominant grass species in this ecosystem our measurements excluded their direct influence on gas exchanges. However, the abundant nature of rushes and sedges suggests that future studies should consider the role and implications