1 Introduction

Tropical peatlands are a globally important carbon stock of 87–350 GtC (Gumbricht et al. 2017; Page et al. 2011) and are crucial for mitigating climate change (Murdiyarso et al. 2013). In addition, tropical peatlands are among the most efficient terrestrial ecosystems for carbon sequestration (Dommain et al. 2011; Page et al. 2004), as continuous input of organic material from lowland tropical evergreen vegetation combined with anaerobic soil conditions lead to a build-up of soil organic matter (SOM) over time (Jauhiainen et al. 2012). In natural conditions, the portion of SOM that is decomposed and emitted as carbon dioxide (CO2) or methane (CH4) is usually outweighed by the continuous input of fresh litter and roots (Jauhiainen et al. 2005; Hergoualc’h and Verchot 2011; Hoyos-Santillan et al. 2015). Models such as the Holocene Peat Model—HPM (Frolking et al. 2010; Kurnianto et al. 2015) use this balance to predict long-term peat accumulation via vegetation-specific characteristics such as decomposition speed and primary production.

Peru is estimated to harbor one of the largest extents of tropical peatlands in the world (Gumbricht et al. 2017; Draper et al. 2014; Page et al. 2011). Mauritia flexuosa-dominated palm swamp forests are the dominant peatland ecosystem type in the Peruvian Amazon (Draper et al. 2014). Interest in the carbon pools of the Peruvian Amazon peatlands has increased in recent years. Since Lähteenoja et al. (2009a) explored their extent, research further expanded into other fields such as palaeoecology (Roucoux et al. 2013), C stocks estimates (Draper et al. 2014), M. flexuosa management (Virapongse et al. 2017), the characterization of degradation (Hergoualc’h et al. 2017), and greenhouse gas (GHG) fluxes (Teh et al. 2017). The challenges that Peruvian peatlands face are substantially different from those in Southeast Asia. Therein Southeast Asia, peatlands are under great pressure from agricultural expansion, artificial drainage, and fires, which result in considerable GHG emissions (Gaveau et al. 2014; Hergoualc’h and Verchot 2014). In contrast, anthropogenic degradation of peatland in the Peruvian Amazon is mostly related to recurrent harvesting of M. flexuosa palms from natural stands without drainage or fire. The fruits from M. flexuosa palms (locally referred to as “Aguaje”) and palm weevils (Rhynchophorus palmarum) that grow inside dead palms are highly demanded products in the regional market and are important sources of vitamins and proteins for rural communities (Pacheco Santos 2005). Even though more sustainable (climbing) techniques exist (Horn et al. 2012), fruit harvesting continues to involve cutting down entire palms. The effect of this practice substantially alters the composition and structure of the forest (Hergoualc’h et al. 2017); however, its impact on long-term peat accumulation remains unstudied.

Information about driving factors of GHG fluxes is crucial to accurately model and predict long-term changes in the soil carbon pool and exchanges of GHG with the atmosphere (Ryan and Law 2005). Hereafter and throughout the text, soil CO2, CH4, and nitrous oxide (N2O) fluxes are collectively referred to as GHG fluxes. Factors that drive GHG flux dynamics include, among others, soil substrate quality, nutrient availability, and aeration. Soil water-filled pore space (WFPS) is commonly used as an indicator for soil aeration. For a wide range of soils, the WFPS has been reported to affect the emissions of CO2 (e.g., Husen et al. 2014; Howard and Howard 1993), CH4 (e.g., Verchot et al. 2000; Del Grosso et al. 2000), and N2O (e.g., van Lent et al. 2015; Davidson et al. 2000). Maximum respiration rates typically occur around 60% WFPS (Linn and Doran 1984); below this level, the microbial activity is limited by water, and above it, oxygen deficiency hampers microbial respiration (Moyano et al. 2013). N2O emissions predominantly arise from nitrification and denitrification (Davidson et al. 2000), and maximum N2O emission reported by Van Lent et al. (2015) was around 60% WFPS for a wide range of tropical mineral soils. Methanogenesis is the anaerobic microbial decomposition of organic material, which occurs in waterlogged soils or in anaerobic microaggregates, whereas methanotrophy takes place in parts of the soil where oxygen is available. Both processes concur in soils, and—among other factors—are modulated by soil structure and WFPS% (Smith et al. 2003; Teh et al. 2005). In waterlogged soils at 100% WFPS, methanogenesis is likely the dominant process, whereas at WFPS < 100% CH4 is solely produced in anoxic microsites and part of the CH4 is oxidized.

GHG production is also influenced by the quality and quantity of fresh roots and litter that enter the soil (Updegraff et al. 1995). For instance, in the DAYCENT model, the carbon-to nitrogen-ratio (C/N ratio) of different plant parts is an input for SOM turnover rates, and consecutively influences C mineralization, which are in turn controlled by factors such as soil moisture, temperature, and texture (Parton et al. 1993, Del Grosso et al. 2001). In addition, Oktarita et al. (2017) found increased soil N2O when soil C/N ratios decreased. We hypothesize that palm harvesting leads to changes in vegetation composition, and in turn to alteration of litter inputs by secondary regrowth and the C/N ratio of the soil substrate. We expect this change in substrate to result in increased decomposition in areas where palms are harvested as compared to undisturbed conditions, which would be reflected in higher GHG emissions in those disturbed sites.

In saturated conditions, aerating roots potentially play an important role for soil GHG fluxes. Plants need to adapt themselves to cope with anoxic conditions (Bruhn et al. 2012). Such adaptations are generally aimed at increasing the influx of oxygen to the root zone through adventitious roots, lenticels, and enlarged aerenchymous tissues (Haase and Rätsch 2010). This acquired oxygen is used for a variety of essential plant functions, such as root respiration (Colmer and Voesenek 2009). A thickened exodermis around the aerenchyma acts like a barrier and only allows for oxygen to leak near and around the root tip (van Noordwijk et al. 1998). Granville (1969) described the presence of such tissues in the aerating pneumatophores of M. flexuosa palms. These pneumatophores may also conduct N2O or CH4 produced in the anoxic soil layers to the atmosphere, as has been shown for aerenchymous tissues in the stalks of rice (see, e.g., Minoda and Kimura 1994) as well as for tree stem lenticels in temperate and tropical forested wetlands (Gauci et al. 2010; Pangala et al. 2013). In this way, CH4 is directly emitted from the soil, and avoids potential oxidation in upper aerobic soil layers. Similarly, dissolved N2O and CH4 can be taken up by roots and subsequently emitted during root or leaf respiration (Pihlatie et al. 2005; Gauci et al. 2010). The existence and the contribution of pneumatophores to total soil GHG emissions have not been reported for palm swamp forests in the Amazon.

Current practices of palm harvesting in the Peruvian Amazon may lead to changes in soil moisture, alteration of the quality of substrate inputs, and decreases in aerating pneumatophores. However, to date, there is no assessment on how these changes could affect GHG emissions. In this context, we set as a first objective for this study to test the response of GHG production to varying levels of WFPS, across a gradient of forest degradation involving substrate alteration. Secondly, we aimed to evaluate the contribution of aerating pneumatophores to surface GHG fluxes in undisturbed and degraded palm swamp sites.

2 Method

We conducted two experiments. The first one was to determine the response of soil microbial GHG emissions to variations in soil moisture and was performed in vitro with samples without live roots. In the second experiment, we estimated the potential of root pneumatophores to conduct GHGs by relating pneumatophore densities to corresponding in situ soil GHG fluxes. The effect of degradation was assessed by including three sites differing in level of degradation. We acknowledge (and further discuss) the limits of our experimental design, which include no field-level replication for each degradation level.

2.1 Study site

The study was conducted southwest of the city of Iquitos, in the province of Loreto in the Peruvian Amazon. The area exhibits a tropical humid climate with an average annual precipitation of 3087 mm and a weak dry season between June and September (Marengo 1998). Mean annual temperature is 27 °C, with average daily minima around 20–22 °C and maxima around 29–31 °C (Marengo 1998). Humidity is at 80–90% year-round.

The palm swamp peatlands were located near Lake Quistococha (S 3° 49.75000′ W 73° 19.11333′). Peat deposits up to 5 m thick have been reported, with the 390–400 cm layer radiocarbon dated at 2335 ± 15 cal. BP (Lähteenoja et al. 2009a). Permanently waterlogged palm swamp originated around 1000 years ago, while the M. flexuosa-dominated vegetation cover observed today was established around 600 years ago (Roucoux et al. 2013). The flooding regime plays an important role for the vegetation development during these time scales. There is currently no man-made drainage in the area and the water table rarely goes deeper than 20 cm below soil surface level (Kelly et al. 2014). The peatlands occasionally flood; Roucoux et al. (2013) mention flooding events in 1998 (30 cm) and 2012 (100 cm). We observed another 100-cm flooding during the 2015 El Niño (van Lent, unpublished data). The area likely receives nutrients during these flooding events, as well as during the annual Amazon River flood pulses. In the region, both minerotrophic (nutrient rich) and ombrotrophic (nutrient poor) peatlands have been reported (Lähteenoja et al. 2009b). Teh et al. (2017) characterized the Quistococha peatland in transition between minerotrophic and ombrotrophic conditions, while and Lähteenoja et al. (2009b) characterized it as minerotrophic.

The present study is part of a long-term experiment monitoring soil GHG fluxes along the previously mentioned gradient of forest degradation. The experimental design comprised an intact site (“I,” S 03° 49.949′ W 073° 18.851′), a moderately degraded site (“mD,” S 03° 50.364′ W 073° 19.501′), and a heavily degraded site (“hD,” S 03° 48.539′ W 073° 18.428′). The sites are part of the same peatland complex and are 1.3–1.7 km distant to the Itaya River, one of the anastomosing channels of the Amazon river (Fig. 1).

Fig. 1
figure 1

Map of the study location near Iquitos, Peru. The locations along the Itaya River of the intact (I), moderately (mD), and highly (hd) degraded sites are indicated with black dots. Both degraded sites are close to settlements (white lines represent streets). The black line around the Quistococha lake (dark gray) delimits the border of the reserve

The I site was located within the Quistococha regional reserve (369 ha), an officially protected area since it was registered as a national touristic park in 1984 (Resolución Suprema No 223-84 ITI/tur). The other two sites were adjacent to the reserve and are utilized by local villagers for extraction of M. flexuosa fruits, collection of palm weevils (R. palmarum), and timber harvesting. The common, unsustainable practice is to cut down entire palms for fruit and weevil collection (Horn et al. 2012). Recurrent harvesting has resulted in a reduced canopy closure at the mD site, which is close to a village founded in 2011. The hD site is close to a village that expanded in 2014, according to satellite images, and has few standing palm trees left.

2.2 Soil and litterfall properties

Peat depth was measured with a soil probe at six locations per site (Bhomia et al. unpublished). Eighteen samples per site were collected from the top 5 cm using a metal ring (radius = 4.5 cm). Samples were dried until constant weight at 60 °C and bulk density was determined from the dry mass per volume of the ring. Six samples per site (three from hollows and three from hummocks) were ground, homogenized, and sent for determination of total carbon (C) and nitrogen (N) by the induction furnace method (Costech EA C-N Analyzer); Calcium (Ca), Potassium (K), Magnesium (Mg), Sodium (Na) by the ammonium acetate method (Ross and Ketterings 2011); and Copper (Cu), Manganese (Mn), Zinc (Zn), and Phosphorus (P) by the Mehlich 3 method (Ziadi and Tran 2007). All analyses were conducted by the University of Hawaii-Hilo. The Ca and Mg content is useful for determining the nutrient status of the sites. A Ca/Mg ratio < 6 classifies as a nutrient poor, rain-fed ombrotrophic peatland; while a ratio > 6 is indicative of inputs from river water and classifies as a minerotrophic peatland (Lähteenoja et al. 2009b).

Vegetation litter was collected from 27 litter traps (total surface = 4.32 m2) per site in September 2016. Although flowering and fruiting can be seasonal, litterfall rates are dominated by leaf litter and samples were representative for our sites. Aguila-Pasquel et al. (2014) reports litterfall rates near our study location and showed that leaf litter C represents ~ 70% ± 2% throughout the year, and 77% in September alone. The litter was weighed and a ~ 100 g subsample was dried in an oven at 40 °C until constant weight. Litterfall C and N concentration was analyzed as described above.

2.3 In vitro incubation

The objective of this experiment was to evaluate the effect of WFPS on GHG flux rates along the gradient of forest degradation. For this purpose, two 500 g soil samples were excavated from the top 30 cm in three plots at ~ 150 m distance from each other within each site. Large roots were manually removed on-site and the remaining soil was air-dried at approximately 25–30 °C for 3 to 5 days. Next, smaller roots and aggregates were removed by gently crushing and sieving the soil (mesh size 2 mm). Then the six samples were mixed to obtain one homogenous sample per site and stored at 7 °C for a week until the start of the incubation.

A pilot incubation experiment was carried out to determine the appropriate incubation time and quantity of air-dried soil needed. The first objective of the pilot was to obtain a linear increase in CO2 concentrations over time, while not exceeding the maximum CO2 standard gas concentration used to calibrate the gas chromatograph (4990 ppm CO2). The second objective was to evaluate the amount of water needed to bring the samples to WFPS values ranging from 20 to 100%.

In the final incubation experiment, triplicates of 15 g air-dried soil were placed in an Erlenmeyer flask (600 ml) and treated with 0, 10, 20, 30, 40, 50, and 70 ml of DI water. The 70-ml treatment intentionally surpassed the water holding capacity to mimic flooded field conditions. Soil and water were mixed with a mechanical shaker and left uncapped to settle for 48 h, to allow the microbial community to stabilize following the drying and rewetting process. Afterwards, GHG flux measurements were performed on four consecutive days. At each day, the flasks were closed and four air samples were taken from the headspace at 1-h intervals. At the onset, flasks were vented with ambient air for 10 s using a vacuum pump and then closed with a rubber stopper equipped with a sampling port. After closure, air samples (10 ml) were taken with a disposable syringe and stored in 10 ml evacuated glass vials with septa caps until further analysis by gas chromatography. Laboratory air temperature and pressure were recorded after the first headspace sample was taken.

Between sampling dates, the flasks were left uncapped at room temperature (~ 21 °C) and the wet soil slurry weight was kept constant through addition D.I. water. At the end of the experiment, samples were oven-dried at 60 °C until constant weight, followed by sample-specific bulk density determination. The WFPS was calculated following the formula by Linn and Doran (1984), assuming a particle density of 1.4 g cm−3 (Driessen and Rochimah 1976).

GHG fluxes from individual flasks at each sampling day were calculated by linear regression of GHG concentration against time, expressed per mass of oven-dried soil. C mineralization was calculated as the site-averaged CO2 production rate per unit C in the soil. WFPS values were calculated for each water addition treatment (n = 12) to test the response of GHG fluxes to WFPS per site.

2.4 Pneumatophore gaseous exchanges experiment

The objective of this experiment was to determine the potential of root pneumatophores to conduct CO2, CH4, or N2O from the soil to the atmosphere. As part of a long-term GHG monitoring study, nine chambers per site were installed in July 2014 (I), September 2014 (hD), and April 2015 (hD) at > 1.5 m distance from a tree or a palm in order to avoid the elevated areas surrounding trunks (hummocks). The lower parts (hollows) are more often flooded, and highest density of pneumatophores is observed in hollows (Granville 1969). These areas represent > 80% of the total surface area in all three sites (data not shown). CO2, CH4, and N2O fluxes were measured using the static closed chamber method, and measured monthly between August and October 2015. The chambers (25 cm height and 30 cm diameter) were pushed 3–5 cm in the soil. Their lids were equipped with a center port for gas sampling and a vent to equalize pressure inside the chamber with that outside. At the start of each measurement, chambers were vented, closed, and gas samples (30 ml) were taken from the enclosed headspace at t = 0, 10, 20, and 30 min using a 50-ml disposable syringe. Twenty milliliters of this sample were injected in pre-evacuated glass vials (10 ml) to store the samples under over-pressurized conditions. The vials were sealed with silicon to prevent leakage during transportation by air from Iquitos to Lima. Air pressure and air and soil temperatures were monitored concomitantly with gas flux measurements. Soil temperature was measured with a probe outside, but within 20 cm of the chamber to prevent soil disturbance.

The number of pneumatophores within each permanent chamber was counted in September 2015. The water table depth (WTD) and WFPS were measured simultaneously with GHG sampling. The WTD was obtained from PVC wells (10 cm diameter, 1.5 m height) installed within 50 cm of each GHG flux chamber. The WFPS was calculated from soil samples collected from the top 10 cm soil layer at each GHG flux chamber, using the same method as described for the in vitro incubation.

2.5 Gas analysis and flux calculation

Gas samples were analyzed at the CIFOR laboratory in Lima, Peru using a gas chromatograph (GC, Perkin Elmer, USA) within 1 week of sample collection on average. The GC was equipped with a 63Ni electronic capture detector (ECD) for N2O analysis and with a flame ionization detector (FID) with a methanizer for analysis of CH4 and CO2. The flux was computed by linear regression against time using the four sampling points. Samples were discarded following a visual quality check for leakage or departure from linearity. Leakage corresponded to a sample for which the concentration of the three GHGs was similar to atmospheric concentrations (except for the sample taken immediately after closure). Departure from linearity of the regression happens because the chamber creates an artifact by reducing the concentration gradient between the soil and the atmosphere (Collier et al. 2014). This usually happens for the last sample time-point.

2.6 Statistical analysis

Statistical analysis was performed using the software IBM SPSS Statistics for Windows 21.0 (IBM Corp. 2012) and statistical significance was set at an alpha level of 5%. Normality was tested using the Shapiro-Wilk test, and visual interpretation of Q-Q plots. Comparisons between sampling dates and replicates were done with ANOVA’s and post-hoc Bonferroni tests for CO2 in case of significant differences between groups. The in vitro CH4, N2O, and the soil nutrients were compared with the Kruskal-Wallis test and pairwise comparisons. The comparisons between sites should be used with care due to the low sample size. Regression models were constructed with the average WFPS and GHG flux per water addition treatment. The response of the fluxes to WFPS was modeled using cubic and quadratic polynomial functions that allow for a decrease in fluxes at high WFPS%. The r2 values were used to estimate which model fit best the observations. The cumulative in situ CH4 fluxes and pneumatophore densities were positively skewed and therefore log-transformed to improve normality. Pneumatophore density was correlated with cumulative CH4 fluxes.

3 Results

3.1 Soil and litterfall properties

Soil and litterfall properties are presented in Table 1. The peat at the hD site was shallower than at the other sites. Bulk densities overall were low which is typical of peat soils. The Ca/Mg ratios were similar among sites and were > 6, indicating minerotrophic conditions. Soil P, C, and C/N ratio at the hD site exhibited lower values than the values at the other sites. Litterfall N and C/N, respectively, increased and decreased with increasing level of degradation. Differences between hummocks and hollows were non-significant, except for Mn, which showed higher concentrations in hummocks.

Table 1 Soil (0–5 cm) and litterfall properties at the intact (I), moderately (mD), and heavily (hD) degraded sites. Averages are given with standard errors (n = 6, bulk density n = 18) and letters indicate significant differences between sites (n = 1). Nutrients considered are Carbon (C), Nitrogen (N), Calcium (Ca), Potassium (K), Magnesium (Mg), Sodium (Na), Copper (Cu), Manganese (Mn), Zinc (Zn), and Phosphorus (P)

3.2 Response of in vitro GHG fluxes from root-free soils to WFPS variation

Average CO2 production across WFPS treatments was consistently lower (p = 0.01) in soils from the I site (1.67 ± 0.18 μg C-CO2 g−1 d.w. h−1), compared to production in soils from the mD (2.34 ± 0.16) and hD sites (2.45 ± 0.16) (Table 2). The difference between sites was more pronounced when expressed as C mineralization: 3.9 ± 0.4, 5.1 ± 0.4, and 7.0 ± 0.5 μg C g−1 C h−1 for the I, mD, and hD sites, respectively. At all sites, the best model fitting the response of CO2 to WFPS was cubic polynomial with maxima between 60 and 70% WFPS (Fig. 2). The difference in CO2 production between sites was largest at the maximum flux, where soils from the mD and hD site produced 134 and 139% of that produced by soils from the I site, respectively. Soil CO2 production remained relatively stable during the incubation period. There were only a few cases where fluxes varied between sampling dates (Appendix).

Table 2 Relationships between water-filled pore space (WFPS) and carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) production from in vitro incubations of soils from the intact (I), moderately (mD), and heavily (hD) degraded sites. Units are:μg C-CO2, ng C-CH4, and ng N-N2O g−1 d.w. h−1. Letters indicate significant differences between sites
Fig. 2
figure 2

In vitro response of soil carbon dioxide (CO2) (a), methane (CH4) (b), and nitrous oxide (N2O) (c) fluxes to varying water-filled pore space (WFPS). Curves indicate statistical models for the intact (I), moderately (mD), and heavily (hD) degraded sites. Soils incubated under flooded conditions are indicated by shaded areas. Error bars present the standard error (n = 12)

Fluxes of CH4 from the hD site were significantly higher compared to fluxes from the I and mD sites, at all WFPS levels (p < 0.01; Table 2). Soils from the I and mD sites were CH4 sinks for all water addition treatments. The soil from the hD site was a net CH4 sink below 70% WFPS (− 0.15 ± 0.09 ng C-CH4 g−1 d.w. h−1), thereafter switching to a net source (6.5 ± 0.6 ng C-CH4 g−1 d.w. h−1). The response of CH4 fluxes to varying WFPS was best described by a quadratic function. For the I and mD sites, the response was very weak, whereas for the hD site CH4 fluxes showed a strong increase above 54% WFPS. Soils from the I and mD sites did not display significant differences in flux rate over time, whereas soil CH4 production from the hD site increased over time (r2 = 0.14, p < 0.01) (Appendix).

Averaged N2O fluxes from the I site (36.68 ± 7.86 ng N-N2O g−1 d.w. h−1) were significantly higher than fluxes from the mD (2.80 ± 0.65) and hD sites (0.08 ± 0.06) for WFPS > 70% (p < 0.01), but not for WFPS < 70% (p = 0.6). N2O production in soils from the I site steeply increased above 70% WFPS and decreased in the flooded treatment. N2O production in soils from the mD and hD site did not respond significantly to variation in WFPS. Production rates reduced over time for the > 67% WFPS treatments of the I (p < 0.01) and mD sites (p < 0.01); soils displayed overall larger fluxes on the first day of incubation compared to other days.

3.3 Contribution of aerating roots to in situ GHG fluxes

Pneumatophore density was higher at the I site (5.6 ± 1.2 pneumatophores dm−2) compared to density at the mD (1.9 ± 0.6) and hD (0.3 ± 0.2) sites (p < 0.01). At the I site, all chambers had at least 0.7 pneumatophore dm−2, whereas at the mD and hD sites 22 and 67% of the chambers, respectively, did not have pneumatophores at all. CH4 fluxes varied strongly within sites, and cumulative CH4 fluxes were not significantly different between sites (p = 0.45). Among all GHG, only CH4 fluxes were significantly correlated to pneumatophore density (Fig. 3). The relationship was most robust for the I site (r2 = 0.5, p = 0.03). Regression lines for pneumatophore densities and CO2 or N2O were insignificant and had nearly horizontal slopes. Water table depth (WTD) varied significantly between sites, in the order I < hD < mD with respective averages of − 21.9 ± 1.6, − 6.5 ± 1.4, and − 2.0 ± 1.2 cm below the soil surface, respectively. Average WFPS followed the same trend as the WTD, with I = hD < mD, and averages of 60 ± 2, 69 ± 4, and 94 ± 2%, respectively.

Fig. 3
figure 3

Relationship between pneumatophore density (Pn) and methane (CH4) fluxes in the intact (I), moderately degraded (mD), and heavily degraded (hD) sites. Greenhouse gas chambers with no pneumatophores were excluded (n = 9, 7, and 3, respectively). Equations for the regression lines are: ln(CH4)all sites = 0.2* ± 0.1 × ln(Pn) + 6.5** ± 0.3 (r2 = 0.2) and ln(CH4)intact = 0.5* ± 0.2 × ln(Pn) + 5.5** ± 0.6 (r2 = 0.5), * = p < 0.05, ** p < 0.01

4 Discussion

The topsoil at all three sites displayed minerotrophic properties, in agreement with findings by Lähteenoja et al. (2009a) at the Quistococha reserve. The sites exhibited similar levels of cations and a similar minerotrophic state. Accordingly, the lower soil C content and C/N ratio at the hD site compared to those at the other sites may indicate a higher degree of humification (Tfaily et al. 2014). The lower P content at the hD site than at the other sites may also suggest a more advanced decomposition status of the peat (Jordan et al. 2007; Könönen et al. 2015). Recurrent harvesting of M. flexuosa palms leads to shifts in forest structure and composition (Hergoualc’h et al. 2017). Bhomia et al. (unpublished) found that the I site harbored a higher density of dicot trees (~ 1500 ha−1) than the degraded sites (~ 700 ha−1). Furthermore, the density of palms overall and of M. flexuosa palms in particular was much lower at the hD site (53 and 16 ha−1, respectively) than at the I and mD sites (~ 300 and 170 ha−1, respectively). The presence of a high number of individuals of Cecropia spp.—a pioneer species observed at the hD site—is indicative of its degraded status. In addition, the gaps formed upon harvesting are colonized by small (< 2 m) fast-growing herbaceous vegetation (personal observation). Such drastic differences in forest structure can be expected to lead to lower above and belowground litter inputs at the degraded sites as compared to inputs at the I site. Differences in forest composition may also induce significant changes in litter quality, as indicated by the difference in litterfall C/N ratio between sites (Table 1). Decreased C/N ratio of litterfall with increased degradation is in agreement with observations by Feldpausch et al. (2004) along a sequence of secondary forest regeneration in Central Amazonia. During early forest succession the biomass accumulates rapidly and consists mainly of easily decomposable litter and less of wood (Gehring et al. 2005; Guariguata and Ostertag 2001). More readily decomposable litter may on the long-term lead to reduced peat accumulation at the degraded sites (Frolking et al. 2001).

The changes in litter quality and the less-decomposed SOM pool as outlined above could explain the differences in C mineralization we found between sites. This is in line with Nilsson and Bohlin (1993), who found that more decomposed peat soils produced more CO2 than less decomposed, more fibrous peat soils in temperate climate. In contrast, incubation studies by Jauhiainen et al. (2016) and Swails et al. (2017) showed that soils from degraded peatlands with secondary regrowth or cultivated with oil palm produced less CO2 than soils from intact peat swamp forests. However, the degraded peatlands of Indonesia were drained and burnt, while sites in our study were not. The effect of such drainage and fire practices generally accelerates the peat mineralization (Hergoualc’h and Verchot 2011, 2014). An intact hydrology at our study area might prohibit the easily-decomposable litter in the degraded sites to decompose. During relative dry periods (60–70% WFPS), the in situ CO2 fluxes in degraded peat swamp forests could show a stronger response to lower WFPS as compared to undisturbed peat swamp forests. Periods with less precipitation are therefore expected to result in increased CO2 fluxes in situ.

The response of CO2 production to changes in WFPS followed a cubic polynomial relationship, which is in agreement with results reported by Howard and Howard (1993) for a wide range of soils and by Husen et al. (2014) for peat soils under oil palm cultivation in Indonesia. The maximum CO2 production at 60–70% WFPS for all three sites is in accordance with common values around 60% WFPS found in the literature (Davidson et al. 1998; Linn and Doran 1984), but slightly higher than the 50% value reported by Husen et al. (2014).

In the incubation study, CH4 was solely produced in soil from the heavily degraded site at anoxic conditions (WFPS > 100%), as has been reported by Smith et al. (2003) and similar to 88% WFPS reported by Melling et al. (2005) for tropical peat swamp forest. The soil samples from the I and mD sites were likely not lacking in C substrate for the methanogens (Table 1), and the water addition treatment mimicking flooded conditions should have ensured anaerobic soil samples. Possibly, CH4 production was inhibited due to high nitrate and other denitrification products, which inhibits methanogenesis in rice paddy soils (Kluber and Conrad 1998a, b). The absence of N2O production in the hD site could indicate a lack of nitrate or other denitrification product in these soils. In addition, the soils from the I and mD produced large amounts of N2O at WFPS% > 70%, which could indicate the inhibition of methanogenesis due to the presence of nitrate or other denitrification products. Most likely, the time for recovery to conditions favoring methanogenesis in the soil samples was inadequate. We did see an effect of increased CH4 production toward the end of the incubation experiment (Appendix), which confirms this lag in recovery.

N2O was predominantly produced at WFPS greater than 70% for soils from the I and mD sites. Mathieu et al. (2006) used 15N stable isotope tracers to demonstrate that in water-saturated soils denitrification accounted for 85–90% of emitted N2O. Therefore, we expect that N2O in our peat soils predominantly originated from denitrification instead of nitrification. In addition, we observed a reduced production with increased incubation time that could indicate a depletion of substrate for denitrification, or a further reduction to dinitrogen (N2) (Appendix). In contrast to the soils from the I and mD sites, the soils from the heavily degraded site did not produce any significant amount of N2O. The large CH4 flux from the same hD soils indicates highly reducing conditions in the soil sample. Possibly denitrification in these soils predominantly resulted in N2, which is the dominant product of denitrification in wet, anaerobic conditions (Davidson et al. 2000), while at the same it is non-toxic for methanogenesis (Roy and Conrad 1999). The maximum N2O production rate occurred at 100 and 84% WFPS for soils from the I and mD sites (Fig. 2c), respectively, much above the 60% value presented by Van Lent et al. (2015) for a wide range of tropical mineral soils. This indicates that nitrification is lacking in these peat soils and N2O predominantly originates from denitrification instead.

The absolute flux rates from the incubated soils are of limited use; however, as presented above, we were mainly interested in the relative flux rates that are useful for interpretation and modeling of field-based fluxes. Nevertheless, the average peat CO2 and CH4 flux rates measured in vitro were in the same order of magnitude as in situ rates from the literature. A hypothetical 30-cm peat profile with the overall average bulk density of 0.1 g cm−3 and the average maximum in vitro production rate of 2.7 μg C-CO2 g−1 d.w. h−1 (Fig. 2a, 60–70% WFPS) would result in CO2 surface fluxes of 81 mg C-CO2 m−2 h−1. This compares well to the average value of 79 mg C-CO2 m−2 h−1 reported by a review of soil surface CO2 from tropical peat swamp forests in Southeast Asia (Hergoualc’h and Verchot 2011). The 30-cm depth is expected to represent the soil surface flux since deep peat layers contribute little to total GHG efflux (Moore and Dalva 1997; Jauhiainen et al. 2016). The average maximum in vitro CH4 production rate of the water-saturated soils (WFPS > 95%) from the hD site was 5.9 ng C-CH4 g−1 d.w. h−1 or 0.18 mg C-CH4 m−2 h−1 which is similar to the rate of 0.33 mg C-CH4 m−2 h−1 reported by Hergoualc’h and Verchot (2014), but an order of magnitude lower than the average of 1.50 mg C-CH4 m−2 h−1 reported from similar M. flexuosa palm swamp forests (Teh et al. 2017). In contrast to CO2 and CH4, N2O fluxes measured in vitro were much higher than rates measured in situ. The average maximum N2O production rate from in vitro incubated soils from the I and mD sites were 105 and 10 ng N-N2O g−1 d.w. h−1 (respectively, 84 and 100% WFPS), or 3.2 and 0.3 mg N-N2O m−2 h−1, respectively. For the soils from the intact site, this was an order of magnitude above the average of 0.74 mg N-N2O m−2 h−1 reported by Van Lent et al. (2015) for tropical peat soils. Differences between fluxes measured in situ and in vitro may originate from soil disturbance during sampling or incubation preparation, differences in active peat layers producing or consuming N2O and environmental conditions, or lack of a plant sink for mineralized N (Butterbach-Bahl et al. 2013). Lastly, an artifact could have been created due to an up-build of GHGs in the soil and water suspension between sampling days, which then diffused into the headspace after venting at the onset of the measurements. However, we took special care to minimize this effect by leaving the flasks open between sampling days.

Cumulative in situ CH4 emissions were correlated with pneumatophore density (Fig. 3), indicating a potential role for aerating roots to conduct CH4 to the atmosphere in these peat swamp forests. Recently, vegetation-based CH4 fluxes have been synthesized to represent 5–22% of the total global CH4 budget, of which 58–78% is conducted by a variety of plant structures and 22–42% is thought to be produced by plants themselves (Carmichael et al. 2014; Schlesinger and Bernhardt 2013). Pangala et al. (2013) roughly estimated that CH4 emissions from woody tree stems alone represented 62–81% of the total ecosystem flux in an Indonesian peat swamp forest, while the contribution of pneumatophores was thought to be negligible (< 2%). The specific contribution of roots to conduct CH4 has been further studied in Indian (Purvaja et al. 2004) and Australian (Kreuzwieser et al. 2003) mangroves; however, to our knowledge, this study is the first to do so for a tropical palm-dominated swamp forest. Our setup was unable to quantify the contribution of these roots to the total ecosystem flux, but it highlights the need for additional research in order to elucidate the role of pneumatophores for CH4 emissions in these ecosystems, and to investigate to what extent this observation holds with different vegetation compositions and soil types throughout the Amazon. Long-term GHG monitoring studies should therefore include aerating roots as potential sources in their experimental design.

The body of knowledge on tropical M. flexuosa-dominated peatlands in the Peruvian Amazon steadily increased in recent years, even though research in this remote region has its practical limitations. This study should be seen as a first exploration of soil moisture and roots as controlling factors for GHG fluxes along a degradation gradient in this region. Despite no site replication, we were able to find enough variation of the underlying drivers of the processes of interest to detect significant relationship. To be able to make inferences about the effects of degradation on biogeochemical cycles and GHG emissions, we would need to undertake a more extensive and replicated experiment. However, this study does show the importance of various site-specific factors that drive the spatial and temporal variations of GHG emissions, and that recurrent harvesting of M. flexuosa could alter the GHG balance on the long-term. Sustainable management of the M. flexuosa-dominated peat swamp forests is needed and should aim at providing a continuous source of income for many families, while conserving the peat and its role for climate change mitigation at the same time.