Articles | Volume 17, issue 21
Research article
10 Nov 2020
Research article |  | 10 Nov 2020

Stem and soil nitrous oxide fluxes from rainforest and cacao agroforest on highly weathered soils in the Congo Basin

Najeeb Al-Amin Iddris, Marife D. Corre, Martin Yemefack, Oliver van Straaten, and Edzo Veldkamp

Although tree stems act as conduits for greenhouse gases (GHGs) produced in the soil, the magnitudes of tree contributions to total (soil + stem) nitrous oxide (N2O) emissions from tropical rainforests on heavily weathered soils remain unknown. Moreover, soil GHG fluxes are largely understudied in African rainforests, and the effects of land-use change on these gases are identified as an important research gap in the global GHG budget. In this study, we quantified the changes in stem and soil N2O fluxes with forest conversion to cacao agroforestry. Stem and soil N2O fluxes were measured monthly for a year (2017–2018) in four replicate plots per land use at three sites across central and southern Cameroon. Tree stems consistently emitted N2O throughout the measurement period and were positively correlated with soil N2O fluxes. 15N-isotope tracing from soil mineral N to stem-emitted 15N2O and correlations between temporal patterns of stem N2O emissions, soil–air N2O concentration, soil N2O emissions and vapour pressure deficit suggest that N2O emitted by the stems originated predominantly from N2O produced in the soil. Forest conversion to extensively managed, mature (>20 years old) cacao agroforestry had no effect on stem and soil N2O fluxes. The annual total N2O emissions were 1.55 ± 0.20 kg N ha−1 yr−1 from the forest and 1.15 ± 0.10 kg N ha−1 yr−1 from cacao agroforestry, with tree N2O emissions contributing 11 % to 38 % for forests and 8 % to 15 % for cacao agroforestry. These substantial contributions of tree stems to total N2O emissions highlight the importance of including tree-mediated fluxes in ecosystem GHG budgets. Taking into account that our study sites' biophysical characteristics represented two-thirds of the humid rainforests in the Congo Basin, we estimated a total N2O source strength for this region of 0.18 ± 0.05 Tg N2O-N yr−1.

1 Introduction

The trace gas nitrous oxide (N2O) has become the main stratospheric ozone-depleting substance produced by human activities (Ravishankara et al., 2009) and is after carbon dioxide and methane (CH4) the most important anthropogenic greenhouse gas (GHG) (Denman et al., 2007). Humid tropical soils are considered one of the most important global N2O sources (Denman et al., 2007; Werner et al., 2007a), with tropical rainforests alone estimated to contribute between 0.9 and 4.5 Tg N2O-N yr−1 to the global N2O source of about 16 Tg N2O-N yr−1 (Bouwman et al., 1995; Breuer et al., 2000; Werner et al., 2007a). However, ground-based, bottom-up N2O emission estimates appear to be in stark contrast to the high emissions estimated from top-down approaches such as modelling and global N2O atmospheric inversions (Huang et al., 2008; Thompson et al., 2014). Nevertheless, there exists considerable uncertainty in both approaches (Davidson and Kanter, 2014), especially for the tropics (Valentini et al., 2014). Recent studies suggest two possible reasons for large uncertainties in bottom-up approaches: “missing” emission pathways such as trees (Welch et al., 2019) and a strong geographic bias of measured N2O fluxes from tropical forests.

Most of the studies on soil N2O fluxes from tropical ecosystems were conducted in South and Central America (Davidson and Verchot, 2000; Matson et al., 2017; Neill et al., 2005; Wolf et al., 2011), tropical Asia (Hassler et al., 2017; Purbopuspito et al., 2006; Veldkamp et al., 2008; Verchot et al., 2006; Werner et al., 2006), and Australia (Breuer et al., 2000; Kiese et al., 2003). Africa remains the continent with the fewest published field studies on soil N2O fluxes from the tropical forest biome. After the pioneering work by Serca et al. (1994), very few field studies have been conducted, most of which were either not replicated with independent plots or only with short measurement campaigns (Castaldi et al., 2013; Gütlein et al., 2018; Wanyama et al., 2018; Werner et al., 2007b). The remaining studies were based on laboratory incubations, which cannot be translated to actual field conditions. Consequently, field-based studies with sufficient spatial and temporal coverage are critical for improving the highly uncertain N2O sink and source estimates for Africa (Kim et al., 2016b; Valentini et al., 2014).

The Congo Basin is the second largest intact tropical forest in the world and constitutes one of the most important carbon (C) and biodiversity reservoirs globally. Behind the DR Congo, Cameroon is the second highest deforested country in the Congo Basin with about 75 % of its forest being subject to pressure from other land uses including agroforestry (Dkamela, 2010). Conversion of forests to traditional cacao agroforestry (CAF) systems has been well documented in Cameroon (Saj et al., 2013; Sonwa et al., 2007; Zapfack et al., 2002). Presently, an estimated 400 000 ha is under CAF on small family farms of approximately 1 to 3 ha (Kotto et al., 2002; Saj et al., 2013). These CAF systems are commonly established under the shade of the forests' remnant trees and are characterized by the absence of fertilizer inputs and low yields of up to 1 t cacao beans ha−1 (Saj et al., 2013).

Changes in land use have been found to affect soil N2O emissions due to changes in soil N availability, vegetation and management practices such as N fertilization (Corre et al., 2006; Davidson and Verchot, 2000; Groffman et al., 2000; Hassler et al., 2017; Veldkamp et al., 2020). In particular, unfertilized agroforestry and agricultural systems have been found to have N2O fluxes comparable to those from the reference forests (Hassler et al., 2017), whereas N-fertilized systems tend to have higher N2O fluxes than the previous forest due to elevated soil mineral N following fertilization (Verchot et al., 2006). This is in line with postulations of the conceptual hole-in-the-pipe (HIP) model, which suggests that the magnitude of N2O emissions from the soil is largely controlled first by soil N availability and second by soil water content (Davidson et al., 2000). A systematic comparison between a reference land use and a converted system for quantifying land-use change effects on GHG fluxes is virtually lacking for the Congo Basin and thus an important knowledge gap in the GHG budget of Africa (Valentini et al., 2014).

Tree stems have been found to act as conduits for soil N2O in wetlands, mangroves and well-drained forests (Kreuzwieser et al., 2003; Rusch and Rennenberg, 1998; Welch et al., 2019), facilitating the transport from the soil, where N2O is produced or consumed by microbial nitrification and denitrification processes, to the atmosphere. Findings of strong declines in N2O emissions with increasing stem height (Barba et al., 2019; Díaz-Pinés et al., 2016; Rusch and Rennenberg, 1998; Wen et al., 2017) suggest that N2O is mainly emitted through the stems and less likely through the leaves. Trees adapted to wetlands and mangroves have aerenchyma systems through which N2O can be transported from the soil into the tree by both gas diffusion and transpiration stream, with exchange to the atmosphere predominantly through the stem lenticels (Rusch and Rennenberg, 1998; Wen et al., 2017). However, for trees on well-drained soils, a different transport mechanism appears to be dominant: transpiration drives the xylem sap flow in which dissolved N2O is transported from the soil to the tree and emitted to the atmosphere through the stem surface and stomata (Machacova et al., 2013; Wen et al., 2017). Recent evidence shows that trees can also act as N2O sinks (Barba et al., 2019; Machacova et al., 2017), highlighting the need for further research of the stem N2O flux magnitudes and their mechanisms.

The most important soil parameters found to influence tree stem N2O fluxes include soil water content (Machacova et al., 2016; Rusch and Rennenberg, 1998), soil N2O fluxes (Díaz-Pinés et al., 2016; Wen et al., 2017), soil temperature (Machacova et al., 2013) and soil–air N2O concentration within the rooting zone (Machacova et al., 2013; Wen et al., 2017). These studies also reported environmental parameters, such as air temperature and vapour pressure deficit, to drive stem N2O fluxes due to their influence on transpiration (O'Brien et al., 2004). For temperate forests on well-drained soil, annual stem N2O fluxes have been found to contribute up to 10 % of the ecosystem N2O emissions (Wen et al., 2017). However, until now, there have been no ground-based spatial extrapolation of the contribution of stem N2O emissions from tropical forests on well-drained soils. Hence, there is a need for concurrent quantifications of the contributions of stem and soil N2O fluxes so as to provide insights into the source strengths of N2O emissions from tropical African land uses and to improve estimates of N2O emissions from the region.

Our present study addresses these knowledge gaps by providing year-round measurements of stem and soil N2O fluxes from forests and converted CAF systems with spatially replicated plots in the Congo Basin as well as stem N2O fluxes of 23 tree species that have not been measured before. Our findings contribute to the much-needed improvement of GHG budget from this region. Our study aimed to (i) assess whether trees in tropical rainforests and CAF are important conduits of N2O, (ii) quantify changes in soil–atmosphere N2O fluxes with forest conversion to CAF, and (iii) determine the temporal and spatial controls of stem and soil N2O fluxes. We hypothesized that (i) stem and soil N2O fluxes from these extensively managed CAF systems (unfertilized and manual harvest) will be comparable to the natural forests, and (ii) the seasonal pattern of stem emissions will parallel that of soil N2O emissions, and both will have similar soil and climatic controlling factors.

2 Materials and methods

2.1 Study area and experimental design

Our study was conducted at three study sites located in southern and central Cameroon, where natural forests are predominantly converted to CAF (Sonwa et al., 2007). Sites in the southern region were located around the villages of Aloum (2.813 N, 10.719 E; 651 m above sea level, a.s.l.) and Biba Yezoum (3.158 N, 12.292 E; 674 m a.s.l.), and the third site was located around the village of Tomba (3.931 N, 12.430 E; 752 m a.s.l.) in the central region (Fig. B1 in Appendix B). The mean annual air temperature across the three sites is 23.5 C (, 2019), and the soil temperature ranged from 21.6 to 24.4 C during our measurement period from May 2017 to April 2018. The study sites span an annual precipitation from 1576 mm yr−1 in central Cameroon to 2064 mm yr−1 in the south of Cameroon (Table A1 in Appendix A;, 2019). Precipitation occurs in a bimodal pattern, with two dry seasons (< 120 mm monthly rainfall) occurring from July to August and December to February. All sites are situated on heavily weathered soils classified as Ferralsols (FAO classification; IUSS Working Group WRB, 2015). Geologically, Tomba and Biba Yezoum are underlain by middle to superior Precambrian basement rocks (metamorphic schists, phyllites and quartzites), whereas the Aloum site is situated on inferior Precambrian basement rocks (inferior gneiss and undifferentiated gneiss) (Gwanfogbe et al., 1983).

At each of the three sites, we studied two land-use systems: the reference forest and the converted CAF system. Additional information on vegetation and site characteristics is reported in Table A1. These CAF sites were established right after clearing the natural forests, where remnant forest trees were retained by farmers to provide shade for understorey cacao trees (Theobroma cacao). Cacao planting and localized weeding were all done manually using hand tools. Interviews of farm owners indicated that there had been no mineral fertilization at any of the CAF sites. The ages of the CAF since conversion varied between 22 and  45 years.

We selected four replicate plots (50 m × 50 m each with a minimum distance of 100 m between plots) per land-use type within each site (Fig. B1), totalling 24 plots that were all located on relatively flat topography. Within each plot, all stems including cacao trees with a diameter at breast height (DBH) ≥10 cm were identified and measured for DBH and height. We conducted N2O flux measurements and soil and meteorological parameters in the inner 40 m × 40 m area within each plot to minimize edge effects. To check that soil conditions were comparable between the reference forests and converted CAF, we compared a land-use-independent soil characteristic, i.e. clay content at 30–50 cm depth, between these land uses at each site. Since we did not find significant differences in clay contents between the forest and CAF at each site (Table 1), we inferred that land-use types within each site had comparable initial soil characteristics prior to conversion, and any differences in N2O fluxes and soil controlling factors can be attributed to land-use conversion.

Table 1Mean (± SE, n=4) soil biochemical characteristics in the top 50 cma depth in forest and cacao agroforestry (CAF) within each site in the Congo Basin, Cameroon. Means followed by different lowercase letters indicate significant differences between land-use types within each site, and different capital letters indicate significant differences among the three sites within a land-use type (ANOVA with Tukey's honest significance difference (HSD) test or Kruskal–Wallis ANOVA with a multiple-comparison extension test at P≤0.05).

a Values are the depth-weighted average, except for clay content (30–50 cm) and stocks of soil organic C and total N, which are the sum of the entire 50 cm depth. b ECEC: effective cation exchange capacity; Exch. bases: sum of exchangeable Ca, Mg, K and Na.

Download Print Version | Download XLSX

For measurements of stem N2O fluxes, we selected six cacao trees per replicate plot in the CAF and six trees representing the most dominant species within each replicate plot in the forest, based on their importance value index (IVI) (Table A2). The species IVI is a summation of the relative density, relative frequency and relative dominance of the tree species (Curtis and McIntosh, 1951). For a given species, the relative density refers to its total number of individuals in the four forest plots at each site, the relative frequency refers to its occurrence among the four forest plots and the relative dominance refers to its total basal area in the four forest plots, all expressed as percentages of all species. These 24 trees measured at each site (6 trees × 4 forest plots) included nine species at the Aloum site, seven species at the Biba Yezoum site and 10 species at the Tomba site (species are specified in Fig. 1; Table A2). The trees were measured for stem N2O fluxes at 1.3 m height above the ground at monthly intervals from May 2017 to April 2018. Furthermore, we assessed the influence of tree height on stem N2O fluxes by conducting additional measurements on 16 individual trees per land use in May 2018; these trees were included in the monthly measurements but were additionally measured at three stem heights (1.3, 2.6 and 3.9 m from the ground) per tree in the forest, and at two heights (1.3 and 2.6 m) per tree in the CAF due to the limited height of the cacao trees.

Figure 1Stem N2O fluxes from 22 tree species at three forest sites (Aloum, Biba Yezoum and Tomba) across central and south Cameroon in the Congo Basin. Boxes (25th, median and 75th percentile) and whiskers (1.5× interquartile range) are based on N2O fluxes measured monthly from May 2017 to April 2018 for each tree species, and the values in parentheses represent the number of trees measured per species. There were no differences in N2O fluxes among species (linear mixed-effect models with Tukey's HSD at P≥0.27).


For soil N2O flux measurements, we installed four permanent chamber bases per replicate plot which were randomly distributed within the inner 40 m × 40 m area. We conducted monthly measurements of soil N2O fluxes from May 2017 to April 2018 as well as meteorological and soil variables known to control N2O emission (see below).

2.2 Measurement of stem and soil N2O fluxes

We measured in situ stem N2O fluxes using stem chambers made from transparent polyethylene terephthalate foil, as described by Wen et al. (2017). About 1 month prior to measurement, we applied acetic acid-free silicone sealant strips (Otto Seal® S110, Hermann Otto GmbH, Fridolfing, Germany) about 1 cm wide and 20 cm apart around the surface of the tree stems (between 1.2 and 1.4 m height from the ground) that stayed permanently to ensure that all the stem chambers had airtight seals (Fig. B2). As many of the measured trees have buttresses (rendering stem chambers impossible to attach at low stem height, e.g. Fig. B2), we chose the measurements at an average of 1.3 m height (or between 1.2 and 1.4 m), congruent to the standard measurement of DBH. Since chamber installation is quick, chambers were newly installed on each sampling date, using the silicone sealant strips as a mark to ensure that the same 0.2 m length stem section was measured. We wrapped a piece of foil (cut approximately 50 cm longer than the measured stem circumference and fitted with a Luer lock sampling port) around each stem. Using a gas-powered heat gun, we “shrank” the top and bottom parts of the foil to fit closely onto the silicone strips, leaving 0.2 m length between the top and bottom silicone strips, which served as the chamber for collecting gas samples (Fig. B2). We then wrapped strips of polyethylene foam around the edges of the foil and adjusted the foam tightly using lashing straps equipped with ratchet tensioners (two straps at the top and two at the bottom). The lashing straps adjusted the flexible foam and the foil (on top of the silicone strips) to any irregularities on the bark and ensured an airtight fitting. After installation, we completely evacuated the air inside the stem chamber using a syringe fitted with a Luer lock one-way check valve. Afterwards, we used a manual hand pump to refill the stem chamber with a known volume of ambient outside air for correct calculation of stem N2O flux. A 25 mL air sample was taken with a syringe through the Luer lock sampling port immediately after refilling the stem chamber with ambient air and then again after 20, 40 and 60 min. Each air sample was immediately stored in pre-evacuated 12 mL Labco Exetainers with rubber septa (Labco Limited, Lampeter, UK), maintaining an overpressure.

In May 2018, we conducted a 15N tracing experiment at the Tomba site as a follow-on study to elucidate the source of stem N2O emissions. The tracing was conducted in three replicate plots per land use, where one tree was selected in each plot. Around each selected tree, 290 mg 15N (in the form of (15NH4)2SO4 with 98 % 15N) dissolved in 8 L distilled water was applied evenly onto the soil surface of 0.8 m2 around the tree using a watering can (equivalent to 10 mm of rain). The water-filled pore space (WFPS) in the top 5 cm depth was 49 ± 1 % and 52 ± 2 % for the forest and CAF, respectively, which were within the range of monthly average WFPS of these plots (Fig. 2i). Based on the monthly average soil mineral N concentrations in these plots, the applied 15N was only 20 % of the extant mineral N in the top 10 cm of soil (resulting in a starting enrichment of 17 % 15N), such that we only minimally changed the substrate, which could influence N2O flux, similar to that described by Corre et al. (2014). Stem and soil 15N2O fluxes were measured 1 d, 7 d and 14 d following 15N application, and on each sampling day gas samples were taken at 0, 30 and 60 min after chamber closure. The gas samples were stored in new pre-evacuated glass containers (100 mL) with rubber septa and transported to the University of Göttingen, Germany, for analysis. We also stored 15N2O standards in similar 100 mL glass containers, which were brought to Cameroon and back to Germany, to have the same storage duration as the gas samples in order to check for leakage; we found no difference in 15N2O with the original standard at our laboratory.

Figure 2Mean (± SE, n=4) stem N2O fluxes (a, b, c), soil N2O fluxes (d, e, f) and water-filled pore space (g, h, i) in Aloum site (a, d, g), Biba Yezoum site (b, e, h) and Tomba site (c, f, i) in the Congo Basin, Cameroon, measured monthly from May 2017 to April 2018; grey shadings mark the dry season.


We measured soil N2O fluxes using vented, static chambers made from polyvinyl chloride that were permanently inserted  0.02 m into the soil at least 1 month prior to the start of measurements, as described in our earlier studies (e.g., Corre et al., 2014; Koehler et al., 2009; Müller et al., 2015). On each sampling day, we covered the chamber bases with vented, static polyethylene hoods (0.04 m2 in area and  11 L total volume) equipped with Luer lock sampling ports. Soil N2O fluxes were then determined by taking four gas samples (25 mL each) at 2, 12, 22 and 32 min after chamber closure. The samples were taken with a syringe and immediately injected into pre-evacuated 12 mL Exetainers as described above.

Concurrent to the stem and soil N2O flux measurements, we sampled soil–air N2O concentrations at 50 cm depth from permanently installed stainless-steel probes (1 mm internal diameter) located  1 m from the measured trees. The stainless-steel probes were installed 1 month prior to the start of measurements. Luer locks were attached to the probes, and on each sampling day the probes were first cleared of any previous accumulation of N2O concentration by removing 5 mL air volume using a syringe and discarding it. We then took 25 mL gas samples and stored them in pre-evacuated 12 mL Exetainers as described above.

2.3 N2O analysis and flux rate calculation

The N2O concentrations in the gas samples were analysed using a gas chromatograph equipped with an electron capture detector, a make-up gas of 5 % CO2–95 % N2 (SRI 8610C, SRI Instruments Europe GmbH, Bad Honnef, Germany) and an autosampler (AS-210, SRI Instruments). 15N2O was analysed on an isotope ratio mass spectrometer (IRMS) (Finnigan DELTAplus XP, Thermo Electron Corporation, Bremen, Germany). We calculated N2O fluxes from the linear change in concentrations over time of chamber closure, and we adjusted the fluxes with air temperature and atmospheric pressure, measured at each replicate plot on each sampling day. We included zero and negative fluxes in our data analysis.

We up-scaled the measured stem N2O fluxes (considering trees ≥10 cm DBH) to annual values on a ground area in the following steps: (1) the relationship between stem N2O fluxes and stem heights was modelled from the 16 individual trees per land use (see above) that were measured at multiple heights, from which we observed decreases in stem N2O fluxes with increasing stem heights. A linear function was statistically the best fit characterizing these decreases in stem N2O fluxes with height. (2) Using this linear function and considering the stem surface area as a frustum with 20 cm increments, the tree-level N2O fluxes on each sampling day were calculated for the regularly measured six trees per plot. (3) The annual tree-level N2O fluxes from these regularly measured six trees per plot were calculated using a trapezoidal interpolation between the tree-level N2O fluxes (step 2) and measurement day intervals from May 2017 to April 2018. (4) The annual tree-level N2O fluxes were then extrapolated on a ground-area basis for each replicate plot as follows (Eq. 1):

(1) annual stem N 2 O flux ( kg N 2 O - N ha - 1 yr - 1 ) = X 1 24 / DBH 1 - 24 24 DBH n A ,

where X1–24 and DBH1–24 are the corresponding annual tree-level N2O flux (kg N2O-N yr−1 of each tree; step 3) and DBH (cm) of each of the 24 measured trees (6 trees × 4 plots) per land use at each site, DBHn is the individual tree DBH (cm) measured for all trees (with ≥10 cm DBH) present within the inner 40 m × 40 m area of each plot (Table A1), Σ is the sum of the annual N2O fluxes of all trees within each plot (kg N2O-N yr−1), and A is the plot area (0.16 ha).

For step 4 of the CAF plots, the annual stem N2O flux was the sum of the cacao and shade trees (Table A1); as these shade trees were remnants of the original forest, we used the average annual tree-level N2O flux of the measured trees in the corresponding paired forest plots multiplied by the actual DBH of the shade trees in the CAF plots. This spatial extrapolation based on trees' DBH of each plot was also supported by the fact that there were no significant differences in stem N2O fluxes among tree species (Fig. 1).

Annual soil N2O fluxes from each plot were calculated using the trapezoidal rule to interpolate the measured fluxes from May 2017 to April 2018, as employed in our earlier studies (e.g., Koehler et al., 2009; Veldkamp et al., 2013). Finally, the annual N2O fluxes from each replicate plot were represented by the sum of the stem and soil N2O fluxes.

2.4 Soil and meteorological variables

We measured soil temperature, WFPS and extractable mineral N in the top 5 cm of depth concurrent to stem and soil N2O flux measurements on each sampling day. The soil temperature was measured  1 m away from the soil chambers using a digital thermometer (GTH 175, Greisinger Electronic GmbH, Regenstauf, Germany). We determined soil WFPS and extractable mineral N by pooling soil samples from four sampling locations within 1 m from each soil chamber in each replicate plot. Gravimetric moisture content was determined by oven-drying the soils at 105 C for 24 h, and WFPS was calculated using a particle density of 2.65 g cm−3 for mineral soil and our measured soil bulk density (Table 1). Soil mineral N (NO3- and NH4+) was extracted in the field by putting a subsample of soil into a pre-weighed bottle containing 150 mL of 0.5 M K2SO4. The bottles were weighed and then shaken for 1 h, and the solution was filtered through pre-washed (with 0.5 M K2SO4) filter papers. The extracts were immediately frozen and later transported to the University of Göttingen, where NH4+ and NO3- concentrations were analysed using continuous-flow injection calorimetry (SEAL Analytical AA3, SEAL Analytical GmbH, Norderstedt, Germany) (described in detail by Hassler et al., 2015). The dry mass of soil extracted for mineral N was calculated using the measured gravimetric moisture content.

During each measurement day, we set up a portable weather station at each site to record relative humidity and air temperature over the course of each sampling day at 15 min intervals. We calculated vapour pressure deficit (VPD) as the difference between saturation vapour pressure (based on its established equation with air temperature) and actual vapour pressure (using saturation vapour pressure and relative humidity; Allen et al., 1998).

Soil biochemical characteristics were measured in April 2017 at all 24 plots. We collected soil samples from the top 50 cm depth, where changes in soil biochemical characteristics resulting from land-use changes have been shown to occur (van Straaten et al., 2015; Tchiofo Lontsi et al., 2019). In each plot, we collected 10 soil samples from the top 0–10 cm and five soil samples each from 10–30 and 30–50 cm depths; in total, we collected 480 soil samples from the 24 plots. The soil samples were air-dried, sieved (2 mm) and transported to the University of Göttingen, where they were dried again at 40 C before analysis. Soil pH was analysed from a 1 : 4 soil-to-distilled water ratio. Soil texture for each plot was determined using the pipette method after iron oxide and organic matter removal (Kroetsch and Wang, 2008). Effective cation exchange capacity (ECEC) and exchangeable cation concentrations (Ca, Mg, K, Na, Al, Fe, Mn) were determined by percolating the soil samples with unbuffered 1 M NH4Cl, and the extracts were analysed using an inductively coupled plasma atomic emission spectrometer (ICP-AES; iCAP 6300 Duo VIEW ICP spectrometer, Thermo Fisher Scientific GmbH, Dreieich, Germany). Soil subsamples were ground and analysed for total organic C and N using a CN analyser (vario EL cube; Elementar Analysis Systems GmbH, Hanau, Germany), and the soil 15N natural abundance signatures were determined using IRMS (DELTAplus; Finnigan MAT, Bremen, Germany). Soil organic carbon (SOC) and total N stocks were calculated for the top 50 cm in both land uses. We used the bulk density of the reference forest for calculating the SOC and total N stocks of the converted CAF in order to avoid overestimations of element stocks resulting from increases in soil bulk densities following land-use conversion (van Straaten et al., 2015; Veldkamp, 1994).

To evaluate the representativeness of our study area with the rest of the Congo Basin rainforest, we estimated the proportion of the Congo rainforest area that has similar biophysical conditions (elevation, precipitation ranges and soil type) as our study sites (Table A1). Using the FAO's Global Ecological Zone map for the humid tropics, we identified the areal coverage of (i) Ferralsols (FAO Harmonized World Soil Database; FAO/IIASA/ISRIC/ISS-CAS/JRC, 2012) with (ii) elevation  1000 m a.s.l. (SRTM digital elevation model; Jarvis et al., 2008) and (iii) a precipitation range between 1500 and 2100 mm yr−1 (WorldClim dataset; Hijmans et al., 2005) within the six Congo rainforest countries (Fig. B3). This analysis was conducted using QGIS version 3.6.3.

2.5 Statistical analyses

Statistical comparisons between land uses or among sites for stem and soil N2O fluxes were performed on the monthly measurements and not on the annual values as the latter are trapezoidal interpolations. As the six trees and four chambers per plot were considered subsamples representing each replicate plot, we conducted the statistical analysis using the means of the six trees and of the four chambers on each sampling day for each replicate plot (congruent to our previous studies, e.g., Hassler et al., 2017; Matson et al., 2017). We tested each parameter for normal distribution (Shapiro–Wilk's test) and homogeneity of variance (Levene's test), and we applied a logarithmic or square root transformation when these assumptions were not met. For the repeatedly measured parameters, i.e. stem and soil N2O fluxes and the accompanying soil variables (temperature, WFPS, NH4+ and NO3- concentrations), differences between land-use types for each site or differences among sites for each land-use type were tested using linear mixed-effect (LME) models with land use or site as the fixed effect and replicate plots and sampling days as random effects (Crawley, 2009). We assessed significant differences between land uses or sites using analysis of variance (ANOVA) with Tukey's HSD test.

We also analysed whether there were differences in stem N2O fluxes among tree species across four forest plots at each site as well as across the three sites. Similar LME analysis was carried out with tree species as the fixed effect, and the random effects were trees belonging to each species and sampling days; only, for this test, we used individual trees as the random effect because most of the tree species (selected based on their IVI; see Sect. 2.1) were not present in all plots, which is typical in species-diverse tropical forest. For soil biochemical characteristics that were measured once (Table 1), one-way ANOVA followed by a Tukey HSD test was used to assess the differences between land uses or sites for the variables with normal distribution and homogenous variance; otherwise, we applied Kruskal–Wallis ANOVA with a multiple-comparison extension test.

To determine the temporal controls of soil and meteorological variables (temperature, WFPS, NH4+ and NO3- concentrations, soil–air N2O concentration, VPD) on stem and soil N2O fluxes, we conducted Spearman's rank correlation tests using the means of the four replicate plots for each land use on each sampling day. For each land use, the correlation tests were conducted across sites and sampling days (n=33, from 3 sites × 11 monthly measurements). To determine the spatial controls of soil biochemical characteristics (which were measured once, Table 1) on stem and soil N2O fluxes, we used the plots' annual N2O emissions and tested with Spearman's rank correlation across land uses and sites (n=24, from three sites times two land uses times four replicate plots). The statistical significance for all the tests was set at P≤0.05. All statistical analyses were conducted using the open-source software R 3.5.2 (R Core Team, 2018).

3 Results

3.1 Stem N2O emissions

Stem N2O emissions differed neither between forest and CAF at each site (P= 0.15–0.76; Table 2) nor among the three sites for each land use (P= 0.16–0.78; Table 2). There were also no differences in stem N2O emissions among tree species in forest plots at each site as well as across the three sites (P= 0.06–0.39; Fig. 1). For the forests, stem N2O emissions exhibited a seasonal pattern with larger fluxes in the wet season than in the dry season at all sites (all P<0.01; Table A3; Fig. 2a, b, c). However, for the CAF, we observed seasonal differences only at the Aloum site (P< 0.01; Table A4; Fig. 2a). Contributions of annual stem N2O emissions reached up to one-third of the total (stem + soil) N2O emissions from the forests (Table 2).

Table 2Mean (± SE, n=4) stem and soil N2O emission as well as annual stem, soil and total (soil + stem) N2O fluxes from forest and cacao agroforestry (CAF) within each site in the Congo Basin, Cameroon. Means followed by different lowercase letters indicate significant differences between land-use types within each site, and different capital letters indicate significant differences among the three sites within a land-use type (linear mixed-effect models with Tukey's HSD at P≤0.05).

* Annual stem and soil N2O fluxes were not statistically tested for differences among sites or between land-use types since these annual values are trapezoidal extrapolations. Annual stem N2O emissions in parentheses are from cacao trees only.

Download Print Version | Download XLSX

From the 15N-tracing experiment, stem 15N-N2O emissions mirrored soil 15N-N2O emissions from both land uses (Fig. 3). About 1 d after 15N addition to the soil, substantial 15N-N2O was emitted from the stem as well as from the soil. This diminished within 2 weeks as the added 15N recycled within the soil N cycling processes, diluting the 15N signatures; nevertheless, the 15N signatures of stem- and soil-emitted N2O remained elevated above the natural abundance level (Fig. 3).

Figure 3Mean (± SE, n=3) 15N2O fluxes from stems (a, unit is per square metre of stem area) and soil (b, unit is per square metre of ground area) in the Congo Basin, Cameroon. In May 2018, 290 mg 15N (in the form of (15NH4)2SO4 with 98 % 15N) was dissolved in 8 L distilled water and sprayed within a 0.8 m2 area around each tree (equal to 10 mm rain), which was only 20 % of the extant mineral N in the top 10 cm of soil and 49 ± 1 % and 52±2 % water-filled pore space for the forest and CAF, respectively, comparable to the soil water content of the site (Fig. 2).


Across the study period, stem N2O emissions from the forests were positively correlated with air temperature, soil–air N2O concentrations and VPD (Table 3) and negatively correlated with WFPS and NH4+ contents (Table 3). The negative correlation of stem N2O emissions with WFPS was possibly spurious, as this correlation may have been driven by the autocorrelation between WFPS and air temperature (Spearman's ρ=-0.59, P< 0.01, n=33). In CAF, stem N2O emissions were only positively correlated with soil N2O emissions (Table 3).

Table 3Spearman correlation coefficients of stem N2O flux (µg N m−2 stem h−1) and soil N2O flux (µg N m−2 soil h−1) with air temperature (C), water-filled pore space (WFPS) (%, top 5 cm depth), extractable NH4+ (mg N kg−1, top 5 cm depth), soil–air N2O concentration (ppm N2O at 50 cm depth) and vapour pressure deficit (VPD) (kPa), using the monthly means of the four replicate plots per land use across the three sites from May 2017 to April 2018 (n=33).

a P≤0.05. b P≤0.01.

Download Print Version | Download XLSX

We detected no difference in WFPS between the forest and CAF (P= 0.15–0.28; Table 4) at any of the sites. For the CAF, we detected higher WFPS in the wet season compared to the dry season at two sites (P< 0.01; Table A4; Fig. 2g, h) whereas there was no seasonal difference in WFPS for the forests at any sites (P= 0.31–0.92; Table A3; Fig. 2g, h, i). At all three sites, the dominant form of mineral N was NH4+ (Table 4). There was generally no difference in soil NH4+ and NO3- between the wet and dry seasons (P= 0.12–0.93), except for the forests at two sites with larger values in the dry than wet season (P< 0.01; Tables S2, S3).

Table 4Mean (± SE, n=4) water-filled pore space (WFPS) and extractable mineral N in the top 5 cm of soil in forest and cacao agroforestry (CAF) within each site in Congo Basin, Cameroon, measured monthly from May 2017 to April 2018.

* Means followed by different lowercase letters indicate significant differences between land-use types within each site and different capital letters indicate significant differences among the three sites within a land-use type (linear mixed-effect models with Tukey's HSD at P≤0.05).

Download Print Version | Download XLSX

3.2 Soil N2O emissions

Soil N2O emissions did not differ between forest and CAF at any site (P= 0.06–0.86; Table 2). Similarly, no differences in soil N2O emissions were detected among sites for each land use (P= 0.26–0.44; Table 2). Soil N2O emissions exhibited consistent seasonal patterns with larger fluxes in the wet than dry season for both land uses (all P< 0.01; Tables S2, S3; Fig. 2d, e, f).

Over the measurement period, soil N2O emissions from the forests were positively correlated with soil–air N2O concentrations and negatively correlated with NH4+ contents (Table 3). In the CAF, soil N2O emissions were positively correlated with WFPS and soil–air N2O concentrations and negatively correlated with air temperatures (Table 3). We did not detect any correlation between annual total N2O fluxes and soil physical and biochemical characteristics. This was not surprising as the ranges of these soil characteristics were relatively small among sites, which reduces the likelihood that significant correlations will be detected.

3.3 Soil biochemical characteristics

Soil physical characteristics (clay content, bulk density) did not differ between forest and CAF at any of the sites (Table 1). Across sites, Biba Yezoum had lower clay content compared to the other sites for each land use (P< 0.01). Generally, the forest showed higher SOC and total N compared to the CAF (P< 0.01–0.05; Table 1).  Soil 15N natural abundance signatures, as an index of the long-term soil N availability, were generally similar between the forest and CAF except at the Aloum site (P< 0.01; Table 1). Soil C  N ratio, another proxy for the long-term soil N status, was higher in the forest than in the CAF at all sites (P< 0.01–0.05). Soil pH and exchangeable bases were lower in the forest compared to the CAF at all sites and the converse was true for exchangeable Al (P< 0.01–0.05; Table 1). Soil ECEC did not differ between the land uses at two sites (P< 0.01; Table 1) and all were low congruent to Ferralsol soils.

4 Discussion

4.1 Stem and soil N2O emissions from the forest

There has been no study on tree stem N2O emission from Africa, nor has any study been reported for the Congo Basin on soil N2O emission with year-round measurements and spatial replication. Stems consistently emitted N2O in both land uses (Table 2; Figs. 1, 2a, b, c), exemplifying that tropical trees on well-drained soils were important contributors of ecosystem N2O emission. So far, there are only two tree species of tropical lowland forest reported with measurements of stem N2O emissions (Welch et al., 2019). Our present study included 23 tree species and their comparable stem N2O emissions, at least from highly weathered Ferralsol soils, across sites over a year of measurements, which provided support to our spatial extrapolation based on DBH of trees at the sites. Mean stem N2O fluxes from our study were within the range of those reported for temperate forests (0.01–2.2 µg N m−2 stem h−1; Díaz-Pinés et al., 2016; Machacova et al., 2016; Wen et al., 2017) but substantially lower than the reported stem N2O emissions of 51–759 µg N m−2 stem h−1 for a humid forest in Panama (Welch et al., 2019). However, Welch et al. (2019) measured stem N2O emissions at a lower stem height (0.3 m) compared to our study (1.3 m), which may partly explain their much larger N2O emissions, as other studies reported that larger N2O emissions occur closer to the stem base of trees (Barba et al., 2019; Díaz-Pinés et al., 2016). Moreover, the consistently higher stem than soil N2O emissions found by Welch et al. (2019), which we did not observe in our study, may point to production of N2O within the stem (e.g., Lenhart et al., 2019). Nonetheless, such high stem N2O emissions as reported by Welch et al. (2019) have not been observed anywhere else under field conditions. We did not find an effect of tree diameter sizes on stem N2O fluxes at our study sites. This was due to the narrow range between the DBH of our measured trees (10–18 cm DBH for cacao trees and 10–30 cm DBH for the forest trees), which reflected the mean stem diameter of trees in our sites (Table A1). Future studies should incorporate trees of wide-ranging diameter size classes, if present at the site, as they may influence N2O flux estimates at the ecosystem scale.

Our annual soil N2O emissions from forests (Table 2) were lower than the reported global average for humid tropical forests (2.81 kg N ha−1 yr−1; summarized by Castaldi et al., 2013). In contrast, the N2O emissions from our forest soils were comparable to those reported for lowland forests on Ferralsol soils in Panama (0.35–1.07 kg N ha−1 yr−1; Matson et al., 2017) and lowland forests on Acrisol soils in Indonesia (0.9 and 1.0 kg N ha−1 yr−1; Hassler et al., 2017). These were possibly due to the soil N availability at our forest sites generally similar to that of forest sites in Panama and Indonesia, indicated by their comparable soil mineral N contents and soil 15N natural abundance signatures.

In comparison with studies from sub-Saharan Africa, annual soil N2O emissions from our forests were lower than the annual N2O emissions reported for the Mayombe forest in the Congo (2.9 kg N ha−1 yr−1; Serca et al., 1994), Kakamega mountain rainforest in Kenya (2.6 kg N ha−1 yr−1; Werner et al., 2007b) and Ankasa rainforest in Ghana (2.3 kg N ha−1 yr−1; Castaldi et al., 2013) but similar in magnitude to those reported for Mau Afromontane forest in Kenya (1.1 kg N ha−1 yr−1; Wanyama et al., 2018). Although these African sites have precipitation levels and highly weathered acidic soils similar to those of our study sites, the Kakamega rainforest in Kenya had higher SOC (7.9 %–20 %) and N contents (0.5 %–1.6 %) in the topsoil layer compared to our forest sites (2.8 %–4.7 % SOC, 0.2 %–0.4 % total N), which may explain its correspondingly higher soil N2O emissions. The study in Congo (Serca et al., 1994), however, was conducted only in a short campaign (2 rainy months and 1 dry month) with less sampling frequency and spatial replication, which may not be a good representation of the spatial and temporal dynamics of soil N2O fluxes to achieve an annual and large-scale estimate.

4.2 Source of tree stem N2O emissions and their contribution to total (stem + soil) N2O emissions

Emitted N2O from stems was found to originate predominantly from N2O produced in the soil, as shown by the 15N-tracing experiment (Fig. 3). Additionally, the positive correlations of stem N2O emissions with soil–air N2O concentrations and soil N2O emissions (Table 3) suggest that the seasonal variation in stem N2O emissions (Table A3; Fig. 2) was likely driven by the temporal dynamics of produced N2O in the soil, which partly supported our second hypothesis. While there have been suggestions of within-tree N2O production (e.g., Lenhart et al., 2019), our finding from the 15N tracing experiment, combined with the correlations of stem N2O emissions with VPD and air temperature, pointed to a transport mechanism of dissolved N2O in soil water by transpiration stream, which has been reported to be important for upland trees that do not have aerenchyma (Machacova et al., 2016; Welch et al., 2019; Wen et al., 2017).

The contributions of up-scaled stem N2O emissions from our studied forests to total (stem + soil) N2O emissions (Table 2) were higher than those reported for temperate forests (1 %–18 %; Díaz-Pinés et al., 2016; Machacova et al., 2016; Wen et al., 2017). Given the higher stem N2O emissions in the wet than dry seasons (Table A3), coupled with the fact that we consistently measured positive fluxes or net stem N2O emissions throughout our measurement period (Fig. 2), we conclude that tree stems in these well-drained Ferralsol soils were efficient conduits for releasing N2O from the soil. This has significant implications in particular during the rainy season as this pathway bypasses the chance for complete denitrification (N2O to N2 reduction) in the soil.

4.3 Factors controlling temporal variability of stem and soil N2O fluxes

The positive correlation of stem N2O emissions with VPD and air temperature in the forest suggests transport of N2O via sap flow, for which the sap flow had been shown to be stimulated with increasing VPD and air temperature (McJannet et al., 2007; O'Brien et al., 2004). Soil water containing dissolved N2O is transported through the xylem via the transpiration stream and eventually emitted from the stem surface to the atmosphere (Díaz-Pinés et al., 2016; Welch et al., 2019; Wen et al., 2017).

Soil moisture has been shown to strongly affect the seasonal variation in soil N2O emissions from tropical ecosystems, with increases in soil N2O emissions by predominantly the denitrification process at high WFPS (Corre et al., 2014; Koehler et al., 2009; Matson et al., 2017; Werner et al., 2006). The larger stem N2O emissions from the forest and soil N2O emissions from both land uses in the wet season than the dry season (Tables S2, S3) signified favourable soil N2O production during the wet season, which suggests that denitrification was the dominant N2O-producing process. However, the moderate WFPS across the year (Table 4) suggests that nitrification may also have contributed to N2O emissions, especially at Biba Yezoum (with lower rainfall and clay contents; Tables 1, S1) where the low WFPS (Table 4) likely favoured nitrification (Corre et al., 2014). For the forest, the negative correlation of the stem and soil N2O emissions with soil NH4+ (Tables 3, S2) may be indicative of a conservative soil N cycle at our forest sites, as supported by the dominance of soil NH4+ over NO3- (Table 2) and by the lower soil N2O emissions at our sites compared to NO3--dominated systems (Davidson et al., 2000). Although the soil mineral N content alone does not indicate the N-supplying capacity of the soil, the relative contents of NH4+ over NO3- can be a good indicator of whether the soil N cycling is conservative with low N2O losses or increasingly leaky (Corre et al., 2010, 2014).

4.4 Land-use change effects on soil N2O emissions

The annual soil N2O emissions from CAF (Table 2) were comparable with those reported for rubber agroforestry in Indonesia (0.6–1.2 kg N ha−1 yr−1; Hassler et al., 2017) and from multistrata agroforestry systems in Peru (0.6 kg N ha−1 yr−1; Palm et al., 2002). However, our soil N2O emissions from CAF were higher than those from an extensively managed home garden in Tanzania (0.35 N ha−1 yr−1; Gütlein et al., 2018). In a review, Kim et al. (2016a) reported mean annual N2O emission from agroforestry systems to be 7.7 kg N ha−1 yr−1. Most of the data used in their review were from intensively managed agroforestry systems with varied fertilizer inputs, which were absent in our extensively managed CAF systems. In line with this, our measured soil N2O emissions from the CAF were also lower than the emissions reported for 10–23-year-old CAF in Indonesia (3.1 kg N ha−1 yr−1; Veldkamp et al., 2008). Our measured N2O emissions provide the first estimates for traditional CAF systems in Africa, as these production systems were not represented in extrapolation of GHG budgets despite their extensive coverage in Africa.

Soil N2O emissions did not differ between forest and CAF systems, which supported our first hypothesis. This is possibly due to the presence of leguminous trees in both systems (Table A1), which can compensate for N export from harvest and other losses (Erickson et al., 2002; Veldkamp et al., 2008). Although studies have hinted at increased N2O emissions from managed systems that utilize leguminous trees as cover crops (Veldkamp et al., 2008), the similar abundance of leguminous trees between forest and CAF at our sites may have offset this effect (Table A1). Previous studies have indeed reported similar soil N2O fluxes between reference forests and unfertilized agroforestry systems (van Lent et al., 2015). Despite the general absence of heavy soil physical disturbance, cultivation and fertilization in these traditional CAF systems, some soil biochemical characteristics have decreased (Table 1); however, these did not translate into detectable differences in soil N2O emissions from those from forest.

4.5 Implications

The biophysical conditions of our forest sites were representative of approximately two-thirds of the rainforest area in the Congo Basin (1.137×106 km2; Fig. B3), considering the same Ferralsol soils, similar elevation ( 1000 m a.s.l.), and annual rainfall between 1500 and 2100 mm yr−1. Using the total (soil + stem) N2O emission from our forest sites (1.55 ± 0.20 N2O-N kg ha−1 yr−1; Table 2), our extrapolated emission for two-thirds of the Congo Basin was 0.18 ± 0.05 Tg N2O-N yr−1 (error estimate is the 95 % confidence interval). This accounted for 52 % of the earlier estimate of soil N2O emissions from tropical rainforests in Africa (0.34 Tg N2O-N yr−1; Werner et al., 2007a), or 25 % based on a more recent estimate (0.72 Tg N2O-N yr−1; Valentini et al., 2014). We acknowledge, however, that there are uncertainties in our extrapolation (as is the case of these cited estimates) because our up-scaling approach from plot to regional level did not account for the spatial variability of large-scale drivers of soil N2O emissions, such as soil texture, landforms and vegetation characteristics. These limitations of our estimate of N2O source strength for the Congo Basin rainforests call for further investigations in Africa to address the geographic bias of studies in the tropical region (e.g., Powers et al., 2011). The most important consideration in the bottom-up spatial extrapolation approach is to recognize at the outset that the design of the field quantification must reflect the landscape-scale drivers of the studied process; e.g. land-use types (reflecting management), soil texture (as a surrogate of parent material) and climate are landscape-scale controllers of soil N, C and GHG fluxes (e.g., Corre et al., 1999; Hassler et al., 2017; Silver et al., 2000; Veldkamp et al., 2008, 2013), whereas topography (reflecting soil types, moisture regimes, fertility) is the main driver within a landscape (e.g., Corre et al., 1996, 2002; Groffman and Tiedje, 1989; Pennock and Corre, 2001). Process-based models and geographic information system databases can be combined with field-based measurements for improved extrapolation.

Our year-round measurements of stem and soil N2O fluxes were the first detailed study carried out in the Congo Basin, with key implications on improved estimates of N2O budget for Africa. Our results revealed that trees on well-drained, highly weathered soils served as an important N2O emission pathway, with the potential to overlook up to 38 % of N2O emissions if trees are not considered in the ecosystem N2O budget. Our measured tree species spanned different life history strategies and functional traits (a mixture of pioneers, non-pioneer light demanders and shade tolerants; Table A2); the lack of species-specific differences suggests that our findings could be more widely generalizable across communities with different species compositions, at least from highly weathered soils. However, the narrow range of tree DBH classes of our measured trees may have important implications for stands of different successional stages or ages, as stem diameter size, wood density and other physiological characteristics may possibly influence stem N2O fluxes (Machacova et al., 2019; Welch et al., 2019). Also, the possibility of large N2O fluxes at the stem base near the ground (Barba et al., 2019; Welch et al., 2019), which we could not measure due to irregular surface of buttresses, warrants further investigation. All this combined may imply that our quantified stem N2O emissions result in a conservative estimate of the overall stem N2O budget from this important region. Forest conversion to traditional, mature (>20 years old) CAF systems had no effect on stem and soil N2O emissions because of similarities in soil moisture and soil texture, absence of fertilizer application, and comparable abundance of leguminous trees in both land uses, which can compensate for N export from harvest or other losses. Further multi-temporal and spatially replicated studies are needed to provide additional insights into the effect of forest conversion to other land uses on GHG fluxes from the African continent in order to improve GHG budget estimations for the region.

Appendix A

Table A1Vegetation and site characteristics of the study sites on highly weathered soils in the Congo Basin, Cameroon. All vegetation characteristics were determined from trees with ≥10 cm diameter at breast height in both forest and cacao agroforestry.

a For cacao agroforestry, the first values are for both cacao and remnant shade trees, and the second values in parentheses are for cacao trees only. b (2019).

Download Print Version | Download XLSX

Table A2Ecological and functional traits of the measured trees, selected from the most dominant tree species at each site, based on their importance value index (IVI is equal to relative density plus relative frequency plus relative dominance; Curtis and McIntosh, 1951). For a given species, the relative density refers to its total number of individuals in the four forest plots at each site, the relative frequency refers to its occurrence among the four forest plots and the relative dominance refers to its total basal area in the four forest plots, all expressed as percentages of all species.

a Each species was assigned to one of the three regeneration guilds defined by Hawthorne (1995): SB: shade-bearer; NPLD: non-pioneer light demander; P: pioneer. b Global Wood Density Database (Brown, 1997; Zanne et al., 2009).

Download Print Version | Download XLSX

Table A3Seasonal mean (± SE, n=4) water-filled pore space (WFPS), extractable mineral N (measured in the top 5 cm of soil) and nitrous oxide (N2O) fluxes in forests on highly weathered soils in the Congo Basin, Cameroon. Means followed by different lowercase letters indicate significant differences between seasons for each site (linear mixed-effect models with Tukey's HSD at P≤0.05).

Download Print Version | Download XLSX

Table A4Seasonal mean (± SE, n=4) water-filled pore space (WFPS), extractable mineral N (measured in the top 5 cm of soil) and nitrous oxide (N2O) fluxes in cacao agroforestry sites located on highly weathered soils in the Congo Basin, Cameroon. Means followed by different lowercase letters indicate significant differences between seasons for each site (linear mixed-effect models with Tukey's HSD at P≤0.05).

Download Print Version | Download XLSX

Appendix B

Figure B1Location of the study sites in Cameroon, showing the four replicate plots per land use (green for forests and orange for cacao agroforestry) at one site.

Figure B2Sampling set-up for stem nitrous oxide (N2O)-flux measurement at three stem heights in a rainforest in the Congo Basin, Cameroon.


Figure B3Map of the Congo Basin rainforest (green) spanning across the six major Congo Basin countries. The brown shaded area represents the proportion of the Congo rainforest with biophysical conditions similar to those in our study sites (Ferralsol soils, ≤1000 m elevation and 1500–2100 mm yr−1 precipitation).

Data availability

Data are available from the Göttingen Research Online repository: (Iddris et al., 2020).

Author contributions

EV and MDC conceived the research project; NAAI carried out fieldwork and analysed data; NAAI and OvS performed GIS analysis; NAAI and MDC interpreted data and wrote the paper; EV, OvS and MY revised the draft paper.

Competing interests

The authors declare that they have no conflict of interest.


We gratefully acknowledge our counterparts in Cameroon, the International Institute for Tropical Agriculture (IITA) for granting us access to and use of their storage facilities. We are especially grateful to our Cameroonian field assistants Leonel Boris Gadjui Youatou, Narcis Lekeng, Yannick Eyenga Alfred, Denis Djiyo and all the fieldworkers for their great support with field measurements as well as Raphael Manu for helping with the GIS work and Rodine Tchiofo Lontsi for many discussions on soil processes and Cameroonian settings. We also thank the village leaders and local plot owners for granting us access to their forest and cacao farms. We thank Andrea Bauer, Kerstin Langs, Martina Knaust and Lars Szwec for their assistance with laboratory analyses.

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant nos. VE 219/14-1 and STR 1375/1-1).

Review statement

This paper was edited by Lutz Merbold and reviewed by Yit Arn Teh and Vincent Gauci.


Allen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Determination of ET0, crop evapotranspiration, Guidel. Comput. Crop Water Requir. Irrig. Drain. Pap. 56, 300 pp., available at: of contens_files.pdf (last access: 3 May 2019), 1998. 

Barba, J., Poyatos, R., and Vargas, R.: Automated measurements of greenhouse gases fluxes from tree stems and soils: magnitudes, patterns and drivers, Sci. Rep.-UK, 9, 1–13,, 2019. 

Bouwman, A. F., Van Der Hoek, K. W., and Olivier, J. G. J.: Uncertainties in the global source distribution of nitrous oxide, J. Geophys. Res., 100, 2785–2800,, 1995. 

Breuer, L., Papen, H., and Butterbach-Bahl, K.: N2O emission from tropical forest soils of Australia, J. Geophys. Res.-Atmos., 105, 26353–26367,, 2000. 

Brown, S.: Estimating biomass and biomass change of tropical forests: a primer, UN FAO Forestry Paper 134, FAO, Rome, 1997. 

Castaldi, S., Bertolini, T., Valente, A., Chiti, T., and Valentini, R.: Nitrous oxide emissions from soil of an African rain forest in Ghana, Biogeosciences, 10, 4179–4187,, 2013. Cameroon climate, available at:, last access: 21 May 2019. 

Corre, M. D., van Kessel, C., and Pennock, D. J.: Landscape and Seasonal Patterns of Nitrous Oxide Emissions in a Semiarid Region, Soil Sci. Soc. Am. J., 60, 1806–1815,, 1996. 

Corre, M. D., Pennock, D. J., Van Kessel, C., and Elliott, D. K.: Estimation of annual nitrous oxide emissions from a transitional grassland-forest region in Saskatchewan, Canada, Biogeochemistry, 44, 29–49,, 1999. 

Corre, M. D., Schnabel, R. R., and Stout, W. L.: Spatial and seasonal variation of gross nitrogen transformations and microbial biomass in a Northeastern US grassland, Soil Biol. Biochem., 34, 445–457,, 2002. 

Corre, M. D., Dechert, G., and Veldkamp, E.: Soil nitrogen cycling following montane forest conversion in Central Sulawesi, Indonesia, Soil Sci. Soc. Am. J., 70, 359–366,, 2006. 

Corre, M. D., Veldkamp, E., Arnold, J., and Joseph Wright, S.: Impact of elevated N input on soil N cycling and losses in old-growth lowland and montane forests in Panama, Ecology, 91, 1715–1729,, 2010. 

Corre, M. D., Sueta, J. P., and Veldkamp, E.: Nitrogen-oxide emissions from tropical forest soils exposed to elevated nitrogen input strongly interact with rainfall quantity and seasonality, Biogeochemistry, 118, 103–120,, 2014. 

Crawley, M. J.: The R Book, John Wiley & Sons Ltd, Chichester, UK, 2009. 

Curtis, J. T. and McIntosh, R. P.: An Upland Forest Continuum in the Prairie-Forest Border Region of Wisconsin, Ecology, 32, 476–496,, 1951. 

Davidson, E. A. and Kanter, D.: Inventories and scenarios of nitrous oxide emissions, Environ. Res. Lett., 9, 105012,, 2014. 

Davidson, E. A. and Verchot, L. V.: Testing the hole-in-the-pipe model of nitric and nitrous oxide emissions from soils using the TRAGNET database, Global Biogeochem. Cy., 14, 1035–1043,, 2000. 

Davidson, E. A., Keller, M., Erickson, H. E., Verchot, L. V., and Veldkamp, E.: Testing a Conceptual Model of Soil Emissions of Nitrous and Nitric Oxides, Bioscience, 50, 667,[0667:tacmos];2, 2000. 

Denman, K. L., Brasseur, G., Chidthaisong, A., Ciais, P., Cox, P. M., Dickinson, R. E., Hauglustaine, D., Heinze, C., Holland, E., Jacob, D., Lohmann, U., Ramachandran, S., da Silva Dias, P. L., Wofsy, S. C., and Zhang, X.: Couplings Between Changes in the Climate System and Biogeochemistry, in Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, available at: (last access: 24 July 2019), 2007. 

Díaz-Pinés, E., Heras, P., Gasche, R., Rubio, A., Rennenberg, H., Butterbach-Bahl, K., and Kiese, R.: Nitrous oxide emissions from stems of ash (Fraxinus angustifolia Vahl) and European beech (Fagus sylvatica L.), Plant Soil, 398, 35–45,, 2016. 

Dkamela, G. P.: The context of REDD+ in Cameroon: Drivers, agents and institutions, Occasional., CIFOR, Bogor, Indonesia, 2010. 

Erickson, H. E., Davidson, E. A., and Keller, M.: Former land-use and tree species affect nitrogen oxide emissions from a tropical dry forest, Oecologia, 130, 297–308,, 2002. 

FAO/IIASA/ISRIC/ISS-CAS/JRC: Harmonized World Soil Database (version 1.2), FAO, Rome, Italy and IIASA, Laxenburg, Austria, available at: (last access: 13 September 2019), 2012. 

Groffman, P. M. and Tiedje, J. M.: Denitrification in north temperate forest soils: Spatial and temporal patterns at the landscape and seasonal scales, Soil Biol. Biochem., 21, 613–620,, 1989. 

Groffman, P. M., Brumme, R., Butterbach-Bahl, K., Dobbie, K. E., Mosier, A. R., Ojima, D., Papen, H., Parton, W. J., Smith, K. A., and Wagner-Riddle, C.: Evaluating annual nitrous oxide fluxes at the ecosystem scale, Global Biogeochem. Cy., 14, 1061–1070,, 2000. 

Gütlein, A., Gerschlauer, F., Kikoti, I., and Kiese, R.: Impacts of climate and land use on N2O and CH4 fluxes from tropical ecosystems in the Mt. Kilimanjaro region, Tanzania, Glob. Change Biol., 24, 1239–1255,, 2018. 

Gwanfogbe, M., Meligui, A., Moukam, J., and Nguoghia, J.: Geography of Cameroon, Macmillan Education Ltd, Hong Kong, 1983. 

Hassler, E., Corre, M. D., Tjoa, A., Damris, M., Utami, S. R., and Veldkamp, E.: Soil fertility controls soil–atmosphere carbon dioxide and methane fluxes in a tropical landscape converted from lowland forest to rubber and oil palm plantations, Biogeosciences, 12, 5831–5852,, 2015. 

Hassler, E., Corre, M. D., Kurniawan, S., and Veldkamp, E.: Soil nitrogen oxide fluxes from lowland forests converted to smallholder rubber and oil palm plantations in Sumatra, Indonesia, Biogeosciences, 14, 2781–2798,, 2017. 

Hawthorne, W. D.: Ecological profiles of Ghanaian forest trees, Tropical forestry papers 29, Oxford Forestry Institute, University of Oxford, Oxford, UK, 1995. 

Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., and Jarvis, A.: Very high resolution interpolated climate surfaces for global land areas, Int. J. Climatol., 25, 1965–1978,, 2005. 

Huang, J., Golombeck, A., Prinn, R. G., Weiss, R. F., Fraser, P. J., Simmonds, P., Dlugokencky, E. J., Hall, B., Elkins, J., Steele, L. P., Langenfelds, R. L., Krummel, P. B., Dutton, G., and Porter, L.: Estimation of regional emissions of nitrous oxide from 1997 to 2005 using multinetwork measurements, a chemical transport model, and an inverse method, J. Geophys. Res.-Atmos., 113, 1–19,, 2008. 

Iddris, N. A., Corre, M. D., Yemefack, M., van Straaten, O., and Veldkamp, E.: Replication Data for: Stem and soil nitrous oxide fluxes from rainforest and cacao agroforest on highly weathered soils in the Congo Basin, Göttingen Research Online/Data, V1,, 2020. 

IUSS Working Group WRB: World Reference Base for Soil Resources 2014, update 2015 International soil classification system for naming soils and creating legends for soil maps, World Soil Resources Reports No. 106, FAO, Rome, 2015. 

Jarvis, A., Reuter, H. I., Nelson, A., and Guevara, E.: Hole-filled SRTM for the globe, Version 4. CGIAR-CSI SRTM 90 m Database, Int. Cent. Trop. Agric. Cali, Columbia, available at: (last access: September 2017), 2008. 

Kiese, R., Hewett, B., Graham, A., and Butterbach-Bahl, K.: Seasonal variability of N2O emissions and CH4 uptake by tropical rainforest soils of Queensland, Australia, Global Biogeochem. Cy., 17, 1043,, 2003. 

Kim, D. G., Kirschbaum, M. U. F., and Beedy, T. L.: Carbon sequestration and net emissions of CH4 and N2O under agroforestry: Synthesizing available data and suggestions for future studies, Agr. Ecosyst. Environ., 226, 65–78,, 2016a. 

Kim, D.-G., Thomas, A. D., Pelster, D., Rosenstock, T. S., and Sanz-Cobena, A.: Greenhouse gas emissions from natural ecosystems and agricultural lands in sub-Saharan Africa: synthesis of available data and suggestions for further research, Biogeosciences, 13, 4789–4809,, 2016b. 

Koehler, B., Corre, M. D., Veldkamp, E., Wullaert, H., and Wright, S. J.: Immediate and long-term nitrogen oxide emissions from tropical forest soils exposed to elevated nitrogen input, Glob. Change Biol., 15, 2049–2066,, 2009. 

Kotto, J. S., Moukam, A., Njomgang, R., Tiki-Manga, T., Tonye, J., Diaw, C., Gockowski, J., Hauser, S., Weise, S. F., Nwaga, D., Zapfack, L., Palm, C. A., Woomer, P., Gillison, A., Bignell, D., and Tondoh, J.: Alternatives to slash-and-burn in Indonesia: summary report & synthesis of phase II in Cameroon, Nairobi, Kenya, 2002. 

Kreuzwieser, J., Buchholz, J., and Rennenberg, H.: Emission of Methane and Nitrous Oxide by Australian Mangrove Ecosystems, Plant Biol., 5, 423–431,, 2003. 

Kroetsch, D. and Wang, C.: Particle size distribution, in: Soil Sampling and Methods of Analysis, 2nd Edn., 713–725, 2008. 

Lenhart, K., Behrendt, T., Greiner, S., Steinkamp, J., Well, R., Giesemann, A., and Keppler, F.: Nitrous oxide effluxes from plants as a potentially important source to the atmosphere, New Phytol., 221, 1398–1408,, 2019. 

Machacova, K., Papen, H., Kreuzwieser, J., and Rennenberg, H.: Inundation strongly stimulates nitrous oxide emissions from stems of the upland tree Fagus sylvatica and the riparian tree Alnus glutinosa, Plant Soil, 364, 287–301,, 2013. 

Machacova, K., Bäck, J., Vanhatalo, A., Halmeenmäki, E., Kolari, P., Mammarella, I., Pumpanen, J., Acosta, M., Urban, O., and Pihlatie, M.: Pinus sylvestris as a missing source of nitrous oxide and methane in boreal forest, Sci. Rep., 6, 1–8,, 2016. 

Machacova, K., Maier, M., Svobodova, K., Lang, F., and Urban, O.: Cryptogamic stem covers may contribute to nitrous oxide consumption by mature beech trees, Sci. Rep.-UK, 7, 1–7,, 2017. 

Machacova, K., Vainio, E., Urban, O., and Pihlatie, M.: Seasonal dynamics of stem N2O exchange follow the physiological activity of boreal trees, Nat. Commun., 10, 1–13,, 2019. 

Matson, A. L., Corre, M. D., Langs, K., and Veldkamp, E.: Soil trace gas fluxes along orthogonal precipitation and soil fertility gradients in tropical lowland forests of Panama, Biogeosciences, 14, 3509–3524,, 2017. 

McJannet, D., Fitch, P., Disher, M., and Wallace, J.: Measurements of transpiration in four tropical rainforest types of north Queensland, Australia, Hydrol. Process., 21, 3549–3564,, 2007. 

Müller, A. K., Matson, A. L., Corre, M. D., and Veldkamp, E.: Soil N2O fluxes along an elevation gradient of tropical montane forests under experimental nitrogen and phosphorus addition, Front. Earth Sci., 3, 1–12,, 2015. 

Neill, C., Steudler, P. A., Garcia-Montiel, D. C., Melillo, J. M., Feigl, B. J., Piccolo, M. C., and Cerri, C. C.: Rates and controls of nitrous oxide and nitric oxide emissions following conversion of forest to pasture in Rondônia, Nutr. Cycl. Agroecosys., 71, 1–15,, 2005. 

O'Brien, J. J., Oberbauer, S. F., and Clark, D. B.: Whole tree xylem sap flow responses to multiple environmental variables in a wet tropical forest, Plant Cell Environ., 27, 551–567,, 2004. 

Palm, C. A., Alegre, J. C., Arevalo, L., Mutuo, P. K., Mosier, A. R., and Coe, R.: Nitrous oxide and methane fluxes in six different land use systems in the Peruvian Amazon, Global Biogeochem. Cy., 16, 1073,, 2002. 

Pennock, D. J. and Corre, M. D.: Development and application of landform segmentation procedures, Soil Till. Res., 58, 151–162,, 2001. 

Powers, J. S., Corre, M. D., Twine, T. E., and Veldkamp, E.: Geographic bias of field observations of soil carbon stocks with tropical land-use changes precludes spatial extrapolation, P. Natl. Acad. Sci. USA, 108, 6318–6322,, 2011. 

Purbopuspito, J., Veldkamp, E., Brumme, R., and Murdiyarso, D.: Trace gas fluxes and nitrogen cycling along an elevation sequence of tropical montane forests in Central Sulawesi, Indonesia, Global Biogeochem. Cy., 20, 1–11,, 2006. 

R Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2018. 

Ravishankara, A. R., Daniel, J. S., and Portmann, R. W.: Nitrous oxide (N2O): The dominant ozone-depleting substance emitted in the 21st century, Science, 326, 123–125,, 2009. 

Rusch, H. and Rennenberg, H.: Black alder (Alnus glutinosa (L.) Gaertn.) trees mediate methane and nitrous oxide emission from the soil to the atmosphere, Plant Soil, 201, 1–7,, 1998. 

Saj, S., Jagoret, P., and Todem Ngogue, H.: Carbon storage and density dynamics of associated trees in three contrasting Theobroma cacao agroforests of Central Cameroon, Agroforest. Syst., 87, 1309–1320,, 2013. 

Serca, D., Delmas, R., Jambert, C., and Labroue, L.: Emissions of nitrogen oxides from equatorial rain forest in central Africa, Tellus B, 46, 243–254,, 1994. 

Silver, W. L., Neff, J., McGroddy, M., Veldkamp, E., Keller, M., and Cosme, R.: Effects of Soil Texture on Belowground Carbon and Nutrient Storage in a Lowland Amazonian Forest Ecosystem, Ecosystems, 3, 193–209,, 2000. 

Sonwa, D. J., Nkongmeneck, B. A., Weise, S. F., Tchatat, M., Adesina, A. A., and Janssens, M. J. J.: Diversity of plants in cocoa agroforests in the humid forest zone of Southern Cameroon, Biodivers. Conserv., 16, 2385–2400,, 2007. 

Tchiofo Lontsi, R., Corre, M. D., van Straaten, O., and Veldkamp, E.: Changes in soil organic carbon and nutrient stocks in conventional selective logging versus reduced-impact logging in rainforests on highly weathered soils in Southern Cameroon, Forest Ecol. Manage., 451, 117522,, 2019. 

Thompson, R. L., Chevallier, F., Crotwell, A. M., Dutton, G., Langenfelds, R. L., Prinn, R. G., Weiss, R. F., Tohjima, Y., Nakazawa, T., Krummel, P. B., Steele, L. P., Fraser, P., O'Doherty, S., Ishijima, K., and Aoki, S.: Nitrous oxide emissions 1999 to 2009 from a global atmospheric inversion, Atmos. Chem. Phys., 14, 1801–1817,, 2014. 

Valentini, R., Arneth, A., Bombelli, A., Castaldi, S., Cazzolla Gatti, R., Chevallier, F., Ciais, P., Grieco, E., Hartmann, J., Henry, M., Houghton, R. A., Jung, M., Kutsch, W. L., Malhi, Y., Mayorga, E., Merbold, L., Murray-Tortarolo, G., Papale, D., Peylin, P., Poulter, B., Raymond, P. A., Santini, M., Sitch, S., Vaglio Laurin, G., van der Werf, G. R., Williams, C. A., and Scholes, R. J.: A full greenhouse gases budget of Africa: synthesis, uncertainties, and vulnerabilities, Biogeosciences, 11, 381–407,, 2014. 

van Lent, J., Hergoualc'h, K., and Verchot, L. V.: Reviews and syntheses: Soil N2O and NO emissions from land use and land-use change in the tropics and subtropics: a meta-analysis, Biogeosciences, 12, 7299–7313,, 2015. 

van Straaten, O., Corre, M. D., Wolf, K., Tchienkoua, M., Cuellar, E., Matthews, R. B., and Veldkamp, E.: Conversion of lowland tropical forests to tree cash crop plantations loses up to one-half of stored soil organic carbon, P. Natl. Acad. Sci. USA, 112, 9956–9960,, 2015. 

Veldkamp, E.: Organic Carbon Turnover in Three Tropical Soils under Pasture after Deforestation, Soil Sci. Soc. Am. J., 58, 175–180,, 1994. 

Veldkamp, E., Purbopuspito, J., Corre, M. D., Brumme, R., and Murdiyarso, D.: Land use change effects on trace gas fluxes in the forest margins of Central Sulawesi, Indonesia, J. Geophys. Res.-Biogeo., 113, 1–11,, 2008. 

Veldkamp, E., Koehler, B., and Corre, M. D.: Indications of nitrogen-limited methane uptake in tropical forest soils, Biogeosciences, 10, 5367–5379,, 2013. 

Veldkamp, E., Schmidt, M., Powers, J. S., and Corre, M. D.: Deforestation and reforestation impacts on soils in the tropics, Nat. Rev. Earth Environ., 1, 590–605,, 2020. 

Verchot, L. V., Hutabarat, L., Hairiah, K., and van Noordwijk, M.: Nitrogen availability and soil N2O emissions following conversion of forests to coffee in southern Sumatra, Global Biogeochem. Cy., 20, 1–12,, 2006. 

Wanyama, I., Pelster, D. E., Arias-Navarro, C., Butterbach-Bahl, K., Verchot, L. V., and Rufino, M. C.: Management intensity controls soil N2O fluxes in an Afromontane ecosystem, Sci. Total Environ., 624, 769–780,, 2018. 

Welch, B., Gauci, V., and Sayer, E. J.: Tree stem bases are sources of CH4 and N2O in a tropical forest on upland soil during the dry to wet season transition, Glob. Change Biol., 25, 361–372,, 2019.  

Wen, Y., Corre, M. D., Rachow, C., Chen, L., and Veldkamp, E.: Nitrous oxide emissions from stems of alder, beech and spruce in a temperate forest, Plant Soil, 420, 423–434,, 2017. 

Werner, C., Zheng, X., Tang, J., Xie, B., Liu, C., Kiese, R., and Butterbach-Bahl, K.: N2O, CH4 and CO2 emissions from seasonal tropical rainforests and a rubber plantation in Southwest China, Plant Soil, 289, 335–353,, 2006. 

Werner, C., Butterbach-Bahl, K., Haas, E., Hickler, T., and Kiese, R.: A global inventory of N2O emissions from tropical rainforest soils using a detailed biogeochemical model, Global Biogeochem. Cy., 21, GB3010,, 2007a. 

Werner, C., Kiese, R., and Butterbach-Bahl, K.: Soil-atmosphere exchange of N2O, CH4, and CO2 and controlling environmental factors for tropical rain forest sites in western Kenya, J. Geophys. Res., 112, D03308,, 2007b. 

Wolf, K., Veldkamp, E., Homeier, J., and Martinson, G. O.: Nitrogen availability links forest productivity, soil nitrous oxide and nitric oxide fluxes of a tropical montane forest in southern Ecuador, Global Biogeochem. Cy., 25, GB4009,, 2011. 

Zanne, A. E., Lopez-Gonzalez, G., Coomes, David A., Ilic, J., Jansen, S., Lewis, S. L., Miller, R. B., Swenson, N. G., Wiemann, M. C., and Chave, J.: Data from: Towards a worldwide wood economics spectrum, v5, Dryad, Dataset,, 2009. 

Zapfack, L., Engwald, S., Sonké, B., Achoundong, G., and Madong, B. A.: The impact of land conversion on plant biodiversity in the forest zone of Cameroon, Biodivers. Conserv., 11, 2047–2061,, 2002. 

Short summary
We quantified the changes in stem and soil nitrous oxide (N2O) fluxes with forest conversion to cacao agroforestry in the Congo Basin, Cameroon. All forest and cacao trees consistently emitted N2O, contributing 8–38 % of the total (soil and stem) emissions. Forest conversion to extensively managed (>–20 years old) cacao agroforestry had no effect on stem and soil N2O fluxes. Our results highlight the importance of including tree-mediated fluxes in the ecosystem-level N2O budget.
Final-revised paper