Potential soil methane oxidation in naturally regenerated oak-dominated temperate deciduous forest stands responds to soil water status regardless of their age—an intact core incubation study
Annals of Forest Science volume 79, Article number: 29 (2022)
Potential CH4 oxidation in the top soil layer increased with decreasing soil water content in spring but was inhibited during severe summer drought in naturally-regenerated oak-dominated temperate deciduous forest stands regardless of their age. No direct effect of mineral nitrogen on soil CH4 oxidation was found. Soil CH4 oxidation in temperate forests could be reduced by extreme climatic events.
The oxidation of atmospheric methane (CH4) by methanotrophic bacteria in forest soils is an important but overlooked ecosystem service.
Our objective was to determine which factors drive variations in soil CH4 oxidation in oak-dominated temperate deciduous forest stands of different ages.
Soil samples were collected in 16 stands aged 20 to 143 years in periods of high and low soil water content (SWC). The potential rate of soil CH4 oxidation was measured by incubating the first five centimetres of intact soil cores at 20 °C.
SWC was the main driver accounting for variations in CH4 oxidation. In spring, a two-fold reduction in SWC greatly increased CH4 oxidation. But when the soil was dry in late summer, a further reduction in SWC led to a decrease in CH4 oxidation in the top soil layer. No direct effect of mineral nitrogen on soil CH4 oxidation was found.
With regard to soil CH4 oxidation, naturally regenerated forest stands contribute equally to climate change mitigation regardless of their age. Considering future climate scenarios for Europe, soil CH4 sink in temperate forests could be reduced, due to both an increase in the number of flooding episodes in spring and drier summers.
The atmospheric concentration of methane (CH4), which is the second most important greenhouse gas contributing to global warming, has increased from 0.75 ppmv to 1.85 ppmv in the past 200 years (Saunois et al. 2020; IPCC 2021). The main CH4 sinks are dependent on oxidation with hydroxyl radicals in the troposphere and biological oxidation in the soil by CH4-oxidising bacteria (Wuebbles 2002). Among the different types of land uses, forest soils show the greatest CH4 oxidation rate at 4.2 kg CH4 ha−1 year−1 while the oxidation rate of other types of soils is about 1.6 kg CH4 ha−1 year−1 (Dutaur and Verchot 2007).
One of the factors that has the greatest impact on CH4 oxidation is the soil’s air-filled porosity, in other words, the open space between soil particles that is filled with air (Epron et al. 2016; Fest et al. 2017). The air-filled porosity depends on both bulk density (BD) and volumetric soil water content (SWC). The diffusion of gases, including CH4 and dioxygen (O2), from the atmosphere into the soil increases with decreasing SWC (Smith et al. 2003), which stimulates CH4 oxidation. On the other hand, high SWC, resulting in low diffusion of gases into the soil, leads to lower methanotrophic activity and can even cause anaerobic conditions that favour the production of CH4 by methanogenic archaea (Veldkamp et al. 2013). A seasonal variation in CH4 oxidation has been observed in several European forests, with higher oxidation in summer and early fall when soils are drier than in winter and earlier spring when SWC is high (Guckland et al. 2009; Epron et al. 2016).
Mineral nitrogen (N) affects CH4 oxidation, either inhibiting or stimulating it depending on the form of the nitrogen and its concentrations, and on the methanotroph community (Carlsen et al. 1991; Hütsch 1998; Saari et al. 2004; Bodelier and Laanbroek 2004; Alam and Jia 2012). Some methanotrophic bacteria switch to ammonia oxidation when ammonia concentrations in the soil are high because ammonia acts as a competitive inhibitor of CH4 monooxygenases (O’Neill and Wilkinson 1977). Nitrite and nitrate can also inhibit CH4 oxidation by methanotrophic bacteria (Seghers et al. 2003; Mochizuki et al. 2012). However, mineral N can also stimulate methanotroph activity when nitrogen is a limiting factor for microbial growth (Bodelier and Laanbroek 2004). In forests, the rate of CH4 oxidation was found to increase with the amount of available nitrogen (Goldman et al. 1995) and when low amounts of nitrogen fertilisers were added (Rigler and Zechmeister-Boltenstern 1999; Veldkamp et al. 2013; Martinson et al. 2021).
Mineral N and SWC vary seasonally and can also change over time with the aging of forest stands. Indeed, mineral N is thought to increase with stand age due to cumulative atmospheric N deposition, increased N mineralisation and the decreased N requirements of trees (Hansen et al. 2007; Christiansen and Gundersen 2011; Rahman et al. 2017). Higher porosity in older stands compared to young stands at least partly explained higher CH4 oxidation in some previous studies (Peichl et al. 2010; Christiansen and Gundersen 2011; Bárcena et al. 2014; Zeng et al. 2014). An increase in soil organic matter over time after afforestation usually results in an increase in soil porosity (Boivin et al. 2009). Moreover, decompaction resulting from root growth may also increase soil porosity (Bodner et al. 2014). Older stands may also have higher evapotranspiration, which would result in lower SWC, therefore improving gas diffusion (Farley et al. 2005; Hiltbrunner et al. 2012). However, it is not certain that SWC decreases as forest stands age because lower evapotranspiration in older stands than in young stands may be due to lower tree densities resulting in a lower leaf area index. Furthermore, greater tree height may result in lower stomatal conductance (Alsheimer et al. 1998; Delzon and Loustau 2005). A few studies have found an increase in CH4 uptake measured in situ with the aging of forest stands linked to a decrease in SWC while others did not find any effect of stand age on CH4 oxidation related to SWC (Peichl et al. 2010; Christiansen and Gundersen 2011; Hiltbrunner et al. 2012; Fest et al. 2015). However, most previous studies on the effect of stand age on CH4 oxidation were conducted on afforested chronosequence established on abandoned agricultural lands, not on naturally regenerated forest stands.
Our objective in this study was to determine the influence of soil bulk density, water content and mineral N on spatial and temporal variations in soil CH4 oxidation in a 470-ha oak-dominated temperate deciduous forest. We measured the potential rate of CH4 oxidation in the top soil layer, SWC, BD and mineral N concentrations on intact soil cores collected in 16 neighbouring stands of sessile oak, aged 20 to 143 years. We hypothesised that (1) variations in CH4 oxidation would be related to differences within and among-stand in SWC, BD and mineral N content, as well as their seasonal variations and (2) CH4 oxidation would increase with stand age, reflecting changes in SWC, BD and mineral N content over time.
2 Materials and methods
2.1 Description of the site and the stands
The study was conducted in the state forest of Champenoux (48°43′ N, 6°21′ E, mean elevation 260 m, 470 ha) in the northeast of France. The soils are neoluvisols on a carbonated clay horizon covered by rather deep loam (Bréda et al. 1993). The climate is semi-continental. Annual mean temperature for the past 30 years is 11 °C and annual average precipitation 808 mm based on data from the nearest weather station (Essey-les-Nancy, 12 km, Météo-France 2020).
Sixteen forest stands dominated by sessile oak (Quercus petraea (Matt.) Liebl.) were selected (Table 1). Thirteen of these stands, previously included in a chronosequence study (Genet et al. 2010), were between 40 and 143 years old; they were studied in fall 2018 and spring 2019. In summer 2019, five of the previous stands were resampled (age 40–60 years) and three younger stands (20 to 26 years old) were added to ensure that the absence of age effect on CH4 oxidation observed in older stands (see the Section 3), was also verified for younger stands. Special care was taken to minimise variations in soil type and topography. All the stands were even-aged and naturally regenerated. In each stand, a 1200-m2 plot was established and the girth at breast height of all the trees in the plot with a girth of more than 10 cm was measured to calculate stand basal area.
2.2 Soil sampling
Soil samples were collected during three campaigns: in early fall 2018 (September 17–October 9) and late summer 2019 (August 27–September 5) during periods of low SWC; and in early spring 2019 (March 25–April 13) during a period of high SWC. In fall 2018 and spring 2019, samples were collected in the oldest 13 of the 16 stands while in summer 2019 samples were collected in the eight youngest stands. Five stands were therefore sampled during the three sampling campaigns.
In each 1200-m2 plot, a 30 × 40 m-grid was installed, then subdivided into 12 squares (10 × 10 m). After removal of the fresh litter, metal cylinders (8 cm in diameter) were used to take soil samples (0–5 cm deep) within each square. Two or three samples in each stand were collected every two or three days over a 14-day period so that sampling day would not be a confounding factor for the stand.
2.3 Potential rates of CH4 oxidation and CO2 production in the top soil layer
Twelve cylinders each containing a soil sample were mounted with two plastic lids, leaving a headspace over the soil sample which, together with the tubing and the analyser cell, occupied a volume of 0.84 L. The 12 cylinders were connected to an off-axis integrated cavity output spectroscopy gas analyser (GGA-24EP, Los Gatos Research, Mountain View, CA, USA) that measured the methane, carbon dioxide and water vapour concentrations in the headspace atmosphere. The incubation of the intact cores was done in a climatic chamber (Memmert, IPP260) where the temperature was kept at 20 °C during the measurements. The gas concentrations in the cylinder headspace were measured for 1 to 2 h, recorded every 5 s and stored on a datalogger (CR1000, Campbell Scientific Inc., Logan, UT, USA). After the measurement of one cylinder was completed, the system was purged and two gases with known concentrations of CH4 (5 and 0 ppm) and CO2 (600 and 0 ppm) were flushed through the analyser for 2 min to calibrate the system before measuring the next cylinder. Solenoid valves (LVM10R1-6B-1-Q, SMC, Tokyo, Japan) controlled with relay controllers (SDM-CD16AC, Campbell Scientific Inc.) connected to the data-logger were used to switch measurements from one cylinder to the next.
Potential rates of CH4 oxidation (in nmol m−2 s−1) and CO2 production (in μmol m−2 s−1) were calculated from the slopes of the variations in concentrations over time, after discarding the first 360 s required to fully equilibrate the system. For CO2, the slope of a linear function was calculated based on the next 240 s of measurements. We preferred a linear function for CO2 over a short time period because CO2 diffusion became limited when the concentrations in the cylinder headspace increased. For CH4, the initial slope was determined by fitting an exponential function on all the remaining data. We chose an exponential function for CH4 because oxidation was expected to decrease along with CH4 concentration in the cylinder headspace due to CH4 limitation.
2.4 Soil mineral nitrogen, water content, and bulk density
The weight of the cylinder with the soil core inside was measured before and after the incubation, and the fresh weight of the soil core was calculated as the difference in mass between the full cylinder and the empty cylinder. After incubation, the samples were sieved through a 2-mm mesh and stored at 4 °C before further analyses. Soil mineral N concentrations were measured on a 20-g subsample of sieved soil. The soil was suspended in 0.1 L of potassium chloride solution (1 M), agitated for one hour in a rotary shaker and filtered through Whatman filter papers. The filtrate was stored at 4 °C for 2 weeks at the most before being analysed with a spectrometric continuous flow analyser (San++, Skalar, 4823 AA Breda, The Netherlands). Two different wavelengths (540 and 660 nm) were used to measure respectively ammonium (NH4+) and nitrite (NO2−) + nitrate (NO3−) nitrogen mass concentrations in the solution. Total mineral N was calculated as the sum of NH4+ and NO2− + NO3−.
To calculate SWC and BD, a 20-g subsample of the fine soil fraction was oven-dried at 105 °C for 24 h immediately following sieving, weighed and used to estimate the dry mass of the soil cores, knowing their total fresh mass and the subsample gravimetric water content.
In summer 2019 only, N mineralisation and net nitrification were measured by incubating at 20 °C the remaining sieved soil for 6 weeks during which SWC was kept constant. NH4+ and NO2− + NO3− concentrations at the end of the incubation period were measured as described above. Net N mineralisation was calculated as the difference between the final mineral N content at the end of the incubation period and the initial mineral N content at the beginning of the incubation period divided by the incubation time. Net nitrification was calculated likewise using NO2− + NO3− nitrogen mass instead of mineral N.
Statistical analyses were performed with the R statistical software version 3.6.1 (R Core Team 2019). We performed one-way ANOVA to test the influence of forest stand on CH4 oxidation, SWC, BD, soil mineral nitrogen and CO2 production in fall 2018, spring 2019 and summer 2019, with 12 replicates per stands. We used linear mixed-effect models (‘lmerTest’ package, Kuznetsova et al. 2017), with stands as random effect, to estimate the effect of the season (fixed effect) on the same variables than above, using data of 13 stands (stand ID 4-16) to compare fall 2018 and spring 2019, and data of five stands (stand ID 4-8) to compare the three sampling dates. In the latter case, contrasts were used to test for seasonal differences when the overall model was significant (p < 0.05).
At the stand scale, we used non-linear models to study the relationship between stand age and either tree density or stand basal area. We used linear models to test the influence of stand age and stand basal area on CH4 oxidation, SWC, mineral N concentration, CO2 production, and BD separately for each sampling dates. We tested the combination of the two factors with or without interaction, and each factor independently, and selected the model for which all effects were significant. If more than one model would have met this requirement, we would have selected the one with the lowest Akaike’s Information Criterion (AIC). Linear regressions between CH4 oxidation and the other variables were tested for each sampling date and are shown when the model was significant (p < 0.05). We fitted an exponential function on the relationship between SWC and nitrification.
At the sample scale, we used structural equation models (SEM) to analyse causal pathways through which soil water content and soil mineral nitrogen have direct effects on CH4 oxidation and CO2 production, and the covariation of these two processes at each sampling campaign. The initial model also included soil bulk density, but it was not included in the final model due to the lack of significant standardised path coefficients at all dates. SEM analyses were conducted using the R package “Lavaan” (Rosseel 2012). Models were retained if the χ2 values of the fits were not significantly different from 0 (p > 0.05), the comparative fit index (cfi) were close to 1, and the root mean square error of approximation (rmsea) and the standardised root mean square residual (srmr) were both close to 0 (Kline 2016).
Stand density decreased sharply with increasing stand age, especially in young stands up to 60 years of age, from 19,367 trees ha−1 to 550 trees ha−1 (Table 1, Fig. 1A). Stand basal area (BA) increased with stand age in young stands, from 12.2 m2 ha−1 to a maximum of 43.7 m2 ha−1 and tended to decrease with age in older stands (Table 1, Fig. 1B). Stand 11 exhibited noticeable high BA.
Potential CH4 oxidation, SWC, mineral N concentrations, and CO2 production in the top soil layer differed significantly from season to season while BD did not (Fig. 2). Higher CH4 oxidation and lower SWC were measured in fall 2018 (respectively 0.66 nmol m−2 s−1 and 0.17 m3 m−3 on average) than in spring 2019 (respectively 0.41 nmol m−2 s−1 and 0.36 m3 m−3 on average; p < 0.001 for both differences) for the 13 stands aged from 40 to 143 years. When considering the five stands that were sampled at the three dates (aged from 40 to 60 years), the potential CH4 oxidation in the top soil layer was lower in summer 2019 than in fall 2018 (respectively 0.42 nmol m−2 s−1 and 0.60 nmol m−2 s−1 on average, p < 0.001), while SWC was lower in summer 2019 than in fall 2018 (respectively 0.14 m3 m−3 and 0.17 m3 m−3 on average, p < 0.001). For the 13 stands 40 to 143 years old, mineral N was higher in spring 2019 (21 mg N kg−1 of dry soil) than in fall 2018 (12 mg N kg−1 of dry soil). However, in the five stands aged from 40 to 60 years, no significant differences in mineral N content were observed between fall 2018 and summer 2019. Ammonium represented on average 89% of the mineral N in our soils (data not shown). CO2 production was higher in spring 2019 than in fall 2018 (respectively 8.4 μmol m−2 s−1 and 2.7 μmol m−2 s−1 on average in 13 stands, p < 0.001), and lower in summer 2019 than in fall 2018 (respectively 0.5 μmol m−2 s−1 and 3.1 μmol m−2 s−1 on average in 5 stands, p < 0.001).
Bulk density varied among stands from 0.8 kg dm−3 to 1.05 kg dm−3, but with no trend related to stand age or stand basal area (data not shown). A significant effect of stand on potential CH4 oxidation and SWC was observed at the three sampling dates (Table 2). In spring 2019, average stand CH4 oxidation in the top soil layer ranged from 0.20 to 0.98 nmol m−2 s− 1 and was positively related to both stand age and basal area (Table 3) but not to stand age alone (Fig. 1C). A similar trend, yet not significant (p < 0.1) was also observed in fall 2018, with a range of average CH4 oxidation in the top soil layer between 0.47 and 0.84 nmol m−2 s−1, but not in the dry summer 2019 (range 0.28 to 0.62 nmol m−2 s−1). SWC was negatively related to stand basal area at all dates, although it was only a trend (p < 0.1) in fall 2018 and summer 2019 (Table 3). Therefore, a significant negative relationship between mean CH4 oxidation and SWC in the top soil layer was observed in spring (Fig. 3A, R2 = 0.82, p < 0.001). It should be noticed that the relationship was strongly influenced by stand 11, which exhibited high rates of CH4 oxidation and low SWC, probably related to its high stand basal area mentioned above. Nevertheless, the relationship remained significant even with the removal of stand 11 (R2 = 0.49, p = 0.011, not shown).
Soil mineral N differed among stands in fall 2018 and summer 2019, but not in spring (Table 2), with noticeable high concentration in stand 11 in fall 2018 and stand 5 in summer 2019. Soil mineral N was positively related to both stand age and basal area in fall 2018 (Table 3) but not at the other seasons. Between-stand variations in potential CH4 oxidation in the top soil layer was significantly related to soil mineral N in summer 2019 (Fig. 3B, R2 = 0.55, p = 0.036). However, the negative relation was strongly influenced by stand 5 and it should be noticed that if we remove stand 5, the relation is no longer significant. The net nitrification decreased exponentially with SWC in summer 2019 while no SWC effect was observed for the mineralisation of N (Fig. 4).
Soil CO2 production, used in this study as a proxy for microbial activity in the top soil layer, differed among stands at the three sampling dates (Table 2), and was negatively related to the stand basal area at all dates, although it was only a trend (p < 0.1) in fall 2018 (Table 3). Between-stand variations in CH4 oxidation was positively related to CO2 production in summer 2019 (Fig. 3C, respectively R2 = 0.68, p = 0.012).
Structural equation model (SEM) analysis conducted at the sample scale (156 soil cores in fall 2018 and spring 2019, and 96 cores in summer 2019) revealed a direct positive effect of SWC on CO2 production at all sampling dates while the direct effect of SWC on potential CH4 oxidation differed between the three sampling dates: not significant in fall 2018, negative in spring 2019, and positive in summer 2019 (Fig. 5). The significant covariation between CH4 oxidation and CO2 production in the top soil layer in summer 2019 confirmed the positive relationship found at stand level mentioned above (Fig. 3C). The direct effect of SWC on soil mineral nitrogen was positive in fall 2018 and spring 2019, but negative in summer 2019. Soil mineral nitrogen had a direct effect on CO2 production in spring 2019 only and had no effect on CH4 oxidation at any sampling date. This suggests that the negative relationship between CH4 oxidation and mineral N in the top soil layer at stand level mentioned above (Fig. 3B) was indirect, possibly mediated by soil water content.
4.1 Variations in methane oxidation were related to variations in soil water content
In agreement with our hypothesis (1), the potential CH4 oxidation measured at 20 °C on intact soil cores was higher in late summer and early fall than in early spring for all stands due to low SWC, which favoured gas diffusion, confirming the relationships observed between net surface fluxes measured in situ and SWC (Smith et al. 2003; Guckland et al. 2009; Epron et al. 2016). The fact that among-stand variations in CH4 oxidation were negatively related to variations in SWC only in early spring suggests that gas diffusion was no longer a limiting factor for CH4 oxidation in the top soil layer in late summer and early fall. The highest rate of CH4 oxidation in spring was observed in stand 11, which also had the lower SWC.
The potential CH4 oxidation in the top soil layer was higher in fall 2018 than in summer 2019 despite a significantly lower SWC in summer 2019 compared to early fall 2018. The increase in gas diffusivity in drier soils explained the increase in soil CH4 uptake measured in situ during moderate summer drought or experimental throughfall reduction (Epron et al. 2016; Fest et al. 2017). However, very low SWC values could have subjected the methanotrophic bacteria to water stress in summer 2019 (van den Pol-van Dasselaar et al. 1998; Borken et al. 2006; von Fischer et al. 2009). We are aware that the potential CH4 oxidation, assessed by incubation of intact cores, is not necessarily representative of the net methane flux that would be measured at the soil surface in situ. Nevertheless, Borken et al. (2006) reported a decrease in ex situ potential CH4 oxidation during severe drought that was also reflected in the in situ seasonal pattern in soil CH4 uptake. The higher CO2 production measured at 20 °C in fall 2018 compared to summer 2019, the positive relationship between CH4 oxidation and CO2 production in summer 2019 and the results of the SEM analysis support the idea that severe summer drought was stressful for the microbial communities living in the top soil layer, including methanotrophic bacteria, at our site. At low SWC (low soil water potential), limited diffusion of soluble organic substrates and osmotic stress impair microbial activity (Skopp et al. 1990; Stark and Firestone 1995; Moyano et al. 2012). Because methanotrophic bacteria gain their energy from the oxidation of CH4, which diffuses better in air than in water, osmotic stress may be the main constraint in late summer and early fall. On the other hand, in early spring when soil water content is high and air-filled porosity low, gas diffusion was the limiting factor for CH4 oxidation.
4.2 Variations in soil mineral nitrogen did not account for the variations in methane oxidation
Though soil mineral N is another important factor that could affect CH4 oxidation since it is known to have both inhibitory and stimulating effects on methanotrophs (Goldman et al. 1995; Mochizuki et al. 2012), we did not observed any clear relation between mineral N content and potential CH4 oxidation. Higher mineral N content in spring than in summer and early fall could be explained by higher N mineralisation in early spring along with low demand from trees (El Zein et al. 2011) and weaker mineralisation at the end of summer and early fall due to unfavourable edaphic conditions following a period of high consumption of belowground resources before the end of the growing season (Ellis 1974).
The weak negative relationship between potential CH4 oxidation and mineral N content in summer 2019, when N content was the lowest and water stress the highest, suggests that unfavourable edaphic conditions in summer 2019 impacted both methanotrophic activity and nitrification. Nitrification is reduced during drought mainly because diffusion of NH4 is limited in dry soils (Stark and Firestone 1995; Schimel 2018). Interestingly, NH4 was always the dominant form of mineral nitrogen at our site; however, in the five stands that were measured both in spring and summer 2019, NH4 represented 95% of the total mineral N in summer, and only 87% in spring. We therefore postulate that this weakly negative relationship between CH4 oxidation and mineral N content resulted from a concomitant effect of soil drying on methanotrophic activity and nitrification in the top soil layer rather than from an inhibitory effect of mineral nitrogen on CH4 oxidation. The lack of inhibitory effect of mineral N was consistent with a much lower mineral N content than in studies showing inhibition where it is usually related to the addition of nitrogen fertiliser (Seghers et al. 2003; Mochizuki et al. 2012).
4.3 Stand age has little influence on methane oxidation in the top soil layer and its putative drivers
Stand age alone did not explain variations of potential CH4 oxidation between stands, but an additive model including both BA and stand age did, except during the 2019 summer drought. We hypothesised that both bulk density and soil water content would decrease with stand age due to the decompaction resulting from root growth (Bodner et al. 2014), the increase in soil organic matter (Boivin et al. 2009) and higher evapotranspiration in older stands than in younger stands (Barker et al. 2009). However, neither bulk density nor soil water content varied with stand age at our site. This lack of stand age effect on bulk density can be explained by the initial low bulk density in our young stands (between 0.8 and 1.1 kg dm−3) compared to previous study of Christiansen and Gundersen (2011) where bulk densities were higher (between 1.2 and 1.7 kg dm−3). These different densities could be explained by the land use: Christiansen and Gundersen (2011) focused on the afforestation of arable lands that had been compacted by agricultural machinery for many years, while our forest stands were naturally regenerated.
In contrast to a previous study on a tree plantation chronosequence on former agricultural land (Hiltbrunner et al. 2012), SWC was not related to stand age in our site. The reason is that, in contrast to forest plantations, all naturally-regenerated stands in our study had a closed canopy, even the youngest ones, due to an initial very high tree density, which decreased with time due to thinning and mortality. SWC was negatively related to the basal area (BA) of the stand. However, as already reported for other forest chronosequences (Roberts et al. 2001; Delzon and Loustau 2005), BA was not continuously increasing with stand age but reached a maximum value before 100 years and decreased thereafter.
Mineral N is thought to increase over time due to a decrease in N demand by trees, cumulative atmospheric N deposition or an increase in N mineralisation but, in contrast to our expectations, there was no effect of stand age on mineral N that could have, in turn, affected CH4 oxidation. An increase in mineral N with stand age has been observed in Danish forests (Hansen et al. 2007; Christiansen and Gundersen 2011; Rahman et al. 2017). This increase was related to a decrease in N demand by trees after canopy closure, following a high demand for N by the trees in their early growth stages when they were building up their N-rich canopy biomass compartments, and to an increase in N mineralisation with age. However, other studies have found a decrease in N mineralisation with stand age because N may be immobilised due to the accumulation of polyphenolic compounds (DeLuca et al. 2002). In contrast to the plantation chronosequence in the Danish study mentioned above, our young naturally-regenerated stands had a closed canopy. In addition, the development of an understorey composed of young hornbeams (30 to 100 cm tall) in our stands older than 86 years (stands 9 to 16) could have increased the N demand in these stands and counterbalanced the lower N demand of the older oak trees. It should also be noticed that nitrogen deposition in our study area (8 kg N ha−1 in 2019 in a nearby beech stand belonging to French permanent plot network for the monitoring of forest ecosystems, RENECOFOR) was lower than in the Danish forests (Hansen et al. 2007; Christiansen and Gundersen 2011).
4.4 Implication for the forest methane sink in the context of climate change and climate change mitigation
It has been suggested that young regrowth forests take up more CO2 than old-growth forests while the latter store more carbon than the former (Schulze 2000; Stephenson et al. 2014; Pugh et al. 2019). However, the CH4 sink function of forest soils has not been included in the debate about whether it is better to conserve old forests or harvest and regenerate them. We found that there was no effect of stand age on potential CH4 oxidation in the top soil layer in a naturally regenerated oak-dominated deciduous temperate forest. Old stands consume as much atmospheric methane as younger stands and contribute equally to climate change mitigation. However, differences in both soil CH4 uptake measured in situ and CH4 oxidation measured ex situ have been reported between forest stands of different ages when forest regeneration occurred after a disturbance, such as severe forest fires, due to the effects of the disturbance on soil properties affecting gas diffusivity (Fest et al. 2015). The current increase in the use of heavy machinery for timber harvesting could alter the uncompacted nature of the soil in the future, potentially reducing methane oxidation in young stands in the next rotation. Frey et al. (2011) and Epron et al. (2016) have indeed reported a negative effect of soil compaction on in situ soil CH4 consumption.
In the past 50 years, the strength and frequency of summer droughts have increased in central and southern Europe (Vicente-Serrano et al. 2014; Trnka et al. 2016). This trend could also become significant in western Europe in the coming years, as suggested by the severe summer droughts in the last 3 years (2018–2020, 45% less summer precipitation and 1.4° more for the average summer temperature compared to the 30-year average based on data from the nearest weather station, Météo-France 2020). In the years to come, the frequency and severity of such droughts is likely to increase (Dai 2011; Trenberth et al. 2014; Spinoni et al. 2018; Grillakis 2019). The decrease in potential CH4 oxidation that we observed in the top soil layer at the end of the very long dry summer in 2019 will presumably occur more frequently, thus reducing the CH4 sink in temperate forest soils. While several studies have reported higher potential CH4 oxidation in the top soil layer than deeper in the soil profile (e.g. Reay et al. 2005) and close agreements between ex situ measurements of potential CH4 oxidation and in situ measurement of soil CH4 uptake (Borken et al. 2006), CH4 oxidation can happen across the whole aerobic soil profile (Fest et al. 2015). We cannot exclude that methanotrophic activity migrated deep into the soil profile during the dry summer, particularly because the high air-filled porosity and reduced oxidation in the upper soil may have favoured CH4 diffusion in deeper layers. This might have at least partly counterbalanced the decreased activity in the top soil layer. As mentioned above, the potential CH4 oxidation, assessed by incubation of intact cores, is not necessarily representative of the net methane flux that would be measured at the soil surface in situ. Further studies on the vertical distribution of soil CH4 oxidation, and its seasonal variation, are needed.
In addition to severe summer drought, an increase in precipitation in late autumn, winter and early spring, and extreme precipitation events in summer despite drier conditions, would not only increase flood hazard but would also increase the frequency and duration of waterlogging episodes (Christensen and Christensen 2003; Dankers and Feyen 2009), reducing the soil CH4 sink of temperate forest soils prone to waterlogging. Indeed, low gas diffusion in waterlogged forest soils is known to result in low in situ net methane uptake (Megonigal and Guenther 2008; Gundersen et al. 2012; Epron et al. 2016). Moreover, a recent study by Gatica et al. (2020) has shown that methanotrophic bacteria are more sensitive to warmer temperatures than methanogenic archaea. Wetter and warmer springs and flash floods in summer could favour methanogenic activity and decrease methane oxidation.
In this context of global warming and related extreme climatic events, our results emphasise that further studies should investigate the effects of long, severe summer droughts and waterlogging events have on the composition and activity of soil microbial communities, and how the interactions between climate and microbial community affect atmospheric methane oxidation, and more generally soil functions and their related forest-ecosystem services.
Alam MS, Jia Z (2012) Inhibition of methane oxidation by nitrogenous fertilizers in a paddy soil. Front Microbiol 3:246. https://doi.org/10.3389/fmicb.2012.00246
Alsheimer M, Köstner B, Falge E, Tenhunen JD (1998) Temporal and spatial variation in transpiration of Norway spruce stands within a forested catchment of the Fichtelgebirge, Germany. Ann For Sci 55:103–123. https://doi.org/10.1051/forest:19980107
Bárcena TG, D’Imperio L, Gundersen P et al (2014) Conversion of cropland to forest increases soil CH4 oxidation and abundance of CH4 oxidizing bacteria with stand age. Appl Soil Ecol 79:49–58. https://doi.org/10.1016/j.apsoil.2014.03.004
Barker CA, Amiro BD, Kwon H et al (2009) Evapotranspiration in intermediate-aged and mature fens and upland black spruce boreal forests. Ecohydrol 2:462–471. https://doi.org/10.1002/eco.74
Bodelier PLE, Laanbroek HJ (2004) Nitrogen as a regulatory factor of methane oxidation in soils and sediments. FEMS Microbiol Ecol 47:265–277. https://doi.org/10.1016/S0168-6496(03)00304-0
Bodner G, Leitner D, Kaul H-P (2014) Coarse and fine root plants affect pore size distributions differently. Plant Soil 380:133–151. https://doi.org/10.1007/s11104-014-2079-8
Boivin P, Schäffer B, Sturny W (2009) Quantifying the relationship between soil organic carbon and soil physical properties using shrinkage modelling. Eur J Soil Sci 60:265–275. https://doi.org/10.1111/j.1365-2389.2008.01107.x
Borken W, Davidson EA, Savage K et al (2006) Effect of summer throughfall exclusion, summer drought, and winter snow cover on methane fluxes in a temperate forest soil. Soil Biol Biochem 38:1388–1395. https://doi.org/10.1016/j.soilbio.2005.10.011
Bréda N, Cochard H, Dreyer E, Granier A (1993) Water transfer in a mature oak stand ( Quercus petraea ): seasonal evolution and effects of a severe drought. Can J For Res 23:1136–1143. https://doi.org/10.1139/x93-144
Carlsen HN, Joergensen L, Degn H (1991) Inhibition by ammonia of methane utilization in Methylococcus capsulatus (Bath). Appl Microbiol Biotechnol 35. https://doi.org/10.1007/BF00180649
Christensen JH, Christensen OB (2003) Severe summertime flooding in Europe. Nature 421:805–806. https://doi.org/10.1038/421805a
Christiansen JR, Gundersen P (2011) Stand age and tree species affect N2O and CH4 exchange from afforested soils. Biogeosciences 8:2535–2546. https://doi.org/10.5194/bg-8-2535-2011
Dai A (2011) Drought under global warming: a review. WIREs Clim Change 2:45–65. https://doi.org/10.1002/wcc.81
Dankers R, Feyen L (2009) Flood hazard in Europe in an ensemble of regional climate scenarios. J Geophys Res 114:D16108. https://doi.org/10.1029/2008JD011523
DeLuca T, Nilsson M-C, Zackrisson O (2002) Nitrogen mineralization and phenol accumulation along a fire chronosequence in northern Sweden. Oecologia 133:206–214. https://doi.org/10.1007/s00442-002-1025-2
Delzon S, Loustau D (2005) Age-related decline in stand water use: sap flow and transpiration in a pine forest chronosequence. Agric For Meteorol 129:105–119. https://doi.org/10.1016/j.agrformet.2005.01.002
Dutaur L, Verchot LV (2007) A global inventory of the soil CH4 sink. Global Biogeochem Cycles 21:GB4013. https://doi.org/10.1029/2006GB002734
El Zein R, Maillard P, Bréda N et al (2011) Are seasonal fluctuations of carbon and nitrogen non-structural compounds coupled within the stem sapwood of 50-year-old sessile oak and beech trees? Tree Physiol. https://doi.org/10.1093/treephys/tpr074
Ellis RC (1974) The seasonal pattern of nitrogen and carbon mineralization in forest and pasture soils in southern Ontario. Can J Soil Sci 54:15–28. https://doi.org/10.4141/cjss74-003
Epron D, Plain C, Ndiaye F-K et al (2016) Effects of compaction by heavy machine traffic on soil fluxes of methane and carbon dioxide in a temperate broadleaved forest. Forest Ecol Manag 382:1–9. https://doi.org/10.1016/j.foreco.2016.09.037
Farley KA, Jobbágy EG, Jackson RB (2005) Effects of afforestation on water yield: a global synthesis with implications for policy. Glob Change Biol 11:1565–1576. https://doi.org/10.1111/j.1365-2486.2005.01011.x
Fest B, Hinko-Najera N, von Fischer JC et al (2017) Soil methane uptake increases under continuous throughfall reduction in a temperate evergreen, broadleaved eucalypt forest. Ecosystems 20:368–379. https://doi.org/10.1007/s10021-016-0030-y
Fest B, Wardlaw T, Livesley SJ et al (2015) Changes in soil moisture drive soil methane uptake along a fire regeneration chronosequence in a eucalypt forest landscape. Glob Change Biol 21:4250–4264. https://doi.org/10.1111/gcb.13003
Frey B, Niklaus PA, Kremer J et al (2011) Heavy-machinery traffic impacts methane emissions as well as methanogen abundance and community structure in oxic forest soils. Appl Environ Microbiol 77:6060–6068. https://doi.org/10.1128/aem.05206-11
Gatica G, Fernández ME, Juliarena MP, Gyenge J (2020) Environmental and anthropogenic drivers of soil methane fluxes in forests: Global patterns and among-biomes differences. Glob Change Biol 26:6604–6615. https://doi.org/10.1111/gcb.15331
Genet H, Bréda N, Dufrêne E (2010) Age-related variation in carbon allocation at tree and stand scales in beech (Fagus sylvatica L.) and sessile oak (Quercus petraea (Matt.) Liebl.) using a chronosequence approach. Tree Physiol 30:177–192
Goldman MB, Groffman PM, Pouyat RV et al (1995) CH4 uptake and N availability in forest soils along an urban to rural gradient. Soil Biol Biochem 27:281–286. https://doi.org/10.1016/0038-0717(94)00185-4
Grillakis MG (2019) Increase in severe and extreme soil moisture droughts for Europe under climate change. Sci Total Environ 660:1245–1255. https://doi.org/10.1016/j.scitotenv.2019.01.001
Guckland A, Flessa H, Prenzel J (2009) Controls of temporal and spatial variability of methane uptake in soils of a temperate deciduous forest with different abundance of European beech (Fagus sylvatica L.). Soil Biol Biochem 41:1659–1667
Gundersen P, Christiansen JR, Alberti G et al (2012) The response of methane and nitrous oxide fluxes to forest change in Europe. Biogeosciences 9:3999–4012. https://doi.org/10.5194/bg-9-3999-2012
Hansen K, Rosenqvist L, Vesterdal L, Gundersen P (2007) Nitrate leaching from three afforestation chronosequences on former arable land in Denmark. Glob Chang Biol 13:1250–1264. https://doi.org/10.1111/j.1365-2486.2007.01355.x
Hiltbrunner D, Zimmermann S, Karbin S et al (2012) Increasing soil methane sink along a 120-year afforestation chronosequence is driven by soil moisture. Glob Change Biol 18:3664–3671. https://doi.org/10.1111/j.1365-2486.2012.02798.x
Hütsch BW (1998) Methane oxidation in arable soil as inhibited by ammonium, nitrite, and organic manure with respect to soil pH. Biol Fertil Soils 28:27–35. https://doi.org/10.1007/s003740050459
IPCC (2021) Climate change 2021: the physical science basis. Contribution of working group i to the sixth assessment report of the intergovernmental panel on climate change. Cambridge University Press, Cambridge
Kline RB (2016) Principles and practice of structural equation modeling, 4th edn. The Guilford Press, New York
Kuznetsova A, Brockhoff PB, Christensen RHB (2017) lmerTest package: tests in linear mixed effects models. J Stat Soft 82:i13. https://doi.org/10.18637/jss.v082.i13
Martinson GO, Müller AK, Matson AL et al (2021) Nitrogen and phosphorus control soil methane uptake in tropical montane forests. J Geophys Res Biogeosci 126. https://doi.org/10.1029/2020JG005970
Megonigal JP, Guenther AB (2008) Methane emissions from upland forest soils and vegetation. Tree Physiol 28:491–498
Météo-France (2020) Normales et records climatologiques 1991-2020 à Nancy-Essey. Infoclimat.fr https://www.infoclimat.fr/climatologie/normales-records/1991-2020/nancy-essey/valeurs/07180.html. Accessed 19 Mar 2020
Mochizuki Y, Koba K, Yoh M (2012) Strong inhibitory effect of nitrate on atmospheric methane oxidation in forest soils. Soil Biol Biochem 50:164–166. https://doi.org/10.1016/j.soilbio.2012.03.013
Moyano FE, Vasilyeva N, Bouckaert L et al (2012) The moisture response of soil heterotrophic respiration: interaction with soil properties. Biogeosciences 9:1173–1182. https://doi.org/10.5194/bg-9-1173-2012
O’Neill JG, Wilkinson JF (1977) Oxidation of ammonia by methane-oxidizing bacteria and the effects of ammonia on methane oxidation. J Gen Microbiol 100:407–412
Peichl M, Arain MA, Ullah S, Moore TR (2010) Carbon dioxide, methane, and nitrous oxide exchanges in an age-sequence of temperate pine forests. Glob Chang Biol 16:2198–2212
Pugh TAM, Lindeskog M, Smith B et al (2019) Role of forest regrowth in global carbon sink dynamics. Proc Natl Acad Sci USA 116:4382–4387. https://doi.org/10.1073/pnas.1810512116
R Core Team (2019) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna
Rahman MM, Bárcena TG, Vesterdal L (2017) Tree species and time since afforestation drive soil C and N mineralization on former cropland. Geoderma 305:153–161. https://doi.org/10.1016/j.geoderma.2017.06.002
Reay D, Nedwell D, Mcnamara N, Ineson P (2005) Effect of tree species on methane and ammonium oxidation capacity in forest soils. Soil Biol Biochem 37:719–730. https://doi.org/10.1016/j.soilbio.2004.10.004
Rigler E, Zechmeister-Boltenstern S (1999) Oxidation of ethylene and methane in forest soils—effect of CO2 and mineral nitrogen. Geoderma 90:147–159. https://doi.org/10.1016/S0016-7061(98)00099-8
Roberts S, Vertessy R, Grayson R (2001) Transpiration from Eucalyptus sieberi (L. Johnson) forests of different age. Forest Ecol Manag 143:153–161
Rosseel Y (2012) lavaan: An R package for structural equation modeling. J Stat Soft 48:i02. https://doi.org/10.18637/jss.v048.i02
Saari A, Rinnan R, Martikainen PJ (2004) Methane oxidation in boreal forest soils: kinetics and sensitivity to pH and ammonium. Soil Biol Biochem 36:1037–1046
Saunois M, Stavert AR, Poulter B et al (2020) The global methane budget 2000–2017. Earth Syst Sci Data 12:1561–1623. https://doi.org/10.5194/essd-12-1561-2020
Schimel JP (2018) Life in dry soils: effects of drought on soil microbial communities and processes. Annu Rev Ecol Evol Syst 49:409–432. https://doi.org/10.1146/annurev-ecolsys-110617-062614
Schulze E-D (2000) Climate change: managing forests after Kyoto. Science 289:2058–2059. https://doi.org/10.1126/science.289.5487.2058
Seghers D, Top EM, Reheul D et al (2003) Long-term effects of mineral versus organic fertilizers on activity and structure of the methanotrophic community in agricultural soils. Environ Microbiol 5:867–877. https://doi.org/10.1046/j.1462-2920.2003.00477.x
Skopp J, Jawson MD, Doran JW (1990) Steady-state aerobic microbial activity as a function of soil water content. Soil Sci Soc Am J 54:1619–1625
Smith KA, Ball T, Conen F et al (2003) Exchange of greenhouse gases between soil and atmosphere: interactions of soil physical factors and biological processes. Eur J Soil Sci 54:779–791. https://doi.org/10.1046/j.1351-0754.2003.0567.x
Spinoni J, Vogt JV, Naumann G et al (2018) Will drought events become more frequent and severe in Europe? Int J Climatol 38:1718–1736. https://doi.org/10.1002/joc.5291
Stark JM, Firestone MK (1995) Mechanisms for soil moisture effects on activity of nitrifying bacteria. Appl Environ Microbiol 61:218–221
Stephenson NL, Das AJ, Condit R et al (2014) Rate of tree carbon accumulation increases continuously with tree size. Nature 507:90–93. https://doi.org/10.1038/nature12914
Trenberth KE, Dai A, van der Schrier G et al (2014) Global warming and changes in drought. Nat Clim Chang 4:17–22. https://doi.org/10.1038/nclimate2067
Trnka M, Balek J, Štěpánek P et al (2016) Drought trends over part of Central Europe between 1961 and 2014. Clim Res 70:143–160. https://doi.org/10.3354/cr01420
van den Pol-van Dasselaar A, van Beusichem ML, Oenema O (1998) Effects of soil moisture content and temperature on methane uptake by grasslands on sandy soils. Plant Soil 204:213–222. https://doi.org/10.1023/A:1004371309361
Veldkamp E, Koehler B, Corre MD (2013) Indications of nitrogen-limited methane uptake in tropical forest soils. Biogeosciences 10:5367–5379. https://doi.org/10.5194/bg-10-5367-2013
Vicente-Serrano SM, Lopez-Moreno J-I, Beguería S et al (2014) Evidence of increasing drought severity caused by temperature rise in southern Europe. Environ Res Lett 9:044001. https://doi.org/10.1088/1748-9326/9/4/044001
von Fischer JC, Butters G, Duchateau PC et al (2009) In situ measures of methanotroph activity in upland soils: A reaction-diffusion model and field observation of water stress. J Geophys Res Biogeosci 114:G01015. https://doi.org/10.1029/2008JG000731
Wuebbles D (2002) Atmospheric methane and global change. Earth Sci Rev 57:177–210. https://doi.org/10.1016/S0012-8252(01)00062-9
Zeng X, Zhang W, Cao J et al (2014) Changes in soil organic carbon, nitrogen, phosphorus, and bulk density after afforestation of the “Beijing–Tianjin Sandstorm Source Control” program in China. Catena 118:186–194. https://doi.org/10.1016/j.catena.2014.01.005
We are grateful to J.-C. Olry, V. Rousselet and F. Bordat for their technical assistance in the field. We also acknowledge Dr N. Breda (INRAE) and C. Lemoine (ONF), who provided us with the age of the stands. We also thank A. Reichard and S. Bienaimé from the “Biogéochimie des Ecosystèmes Forestiers” laboratory who managed the measurement of the mineral N concentrations.
The custom code generated during the current study is available from the corresponding author on reasonable request.
Nicolas Bras was supported by a grant from the French PIA project “Lorraine Université d’Excellence” (ANR-15-IDEX-04-LUE). UMR 1434 SILVA is supported by the French National Research Agency through the Laboratory of Excellence ARBRE (ANR-11-LABX-0002-01).
Ethic approval and consent to participate
The authors declare that the study was not conducted on endangered, vulnerable, or threatened species.
Consent for publication
All authors gave their informed consent to this publication and its content.
The authors declare that they have no conflict of interest.
Handling Editor: Alexia Stokes
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Bras, N., Plain, C. & Epron, D. Potential soil methane oxidation in naturally regenerated oak-dominated temperate deciduous forest stands responds to soil water status regardless of their age—an intact core incubation study. Annals of Forest Science 79, 29 (2022). https://doi.org/10.1186/s13595-022-01145-9
- Methane oxidation
- Naturally regenerated forest
- Air-filled porosity
- Mineral nitrogen
- Water stress