Conservative N cycling despite high atmospheric deposition in early successional African tropical lowland forests

Across the tropics, the share of secondary versus primary forests is strongly increasing. The high rate of biomass accumulation during this secondary succession relies on the availability of essential nutrients, such as nitrogen (N). Nitrogen primarily limits many young secondary forests in the tropics. However, recent studies have shown that forests of the Congo basin are subject to high inputs of atmospheric N deposition, potentially alleviating this N limitation in early succession. To address this hypothesis, we assessed the N status along a successional gradient of secondary forests in the Congo basin. In a set-up of 18 plots implemented along six successional stages, we quantified year-round N deposition, N leaching, N2O emission and the N flux of litterfall and fine root assimilation. Additionally, we determined the N content and C:N stoichiometry for canopy leaves, fine roots, and litter, as well as δ15N of canopy leaves. We confirmed that these forests receive high amounts of atmospheric N deposition, with an increasing deposition as forest succession proceeds. Additionally, we noted lower C:N ratios, and higher N leaching losses, N2O emission, and foliar δ15N in older secondary forest (60 years). In contrast, higher foliar, litter and root C:N ratios, and lower foliar δ15N, N leaching, and N2O emission in young (< 20 years) secondary forest were observed. Altogether, we show that despite high N deposition, this early forest succession still shows conservative N cycling characteristics, which are likely indicating N limitation early on in secondary forest succession. As secondary succession advances, the N cycle gradually becomes more open.


Introduction
The Democratic Republic of the Congo is home to more than 60% of the Congo basin's forest, representing the second largest contiguous block of tropical forest on the planet (de Wasseige et al. 2012). These forests are important for hosting biodiversity (Slik et al. 2015), and play an important role in the global carbon cycle. In particular, the carbon sink of the intact mature African forests has shown to be stable over the past decades while the Amazonian carbon sink was declining (Hubau et al. 2020). However, these forests and the ecosystem services they provide are subject to increasing deforestation. With an ongoing population expansion in Central Africa (Gerland et al. 2014), shifting cultivation will continue to decrease the extent of pristine forests and increasingly convert it to secondary forests in the next decades (Curtis et al. 2018;Tyukavina et al. 2018). In general, secondary tropical forests exhibit high rates of biomass accumulation (Martin et al. 2013;Poorter et al. 2016), but this productivity shows regional variation and peaks in the young central African forests (7.6 Mg above-ground biomass ha −1 yr −1 , Suarez et al. 2019). This high productivity is sustained by a high uptake of essential nutrients, especially early on during the secondary succession (Mcgrath et al. 2001). Nitrogen (N) limits primary production of at least 18% of terrestrial ecosystems (Du et al. 2020) and N limitation may reduce the projected increase in terrestrial C uptake due to increasing CO 2 fertilization by 24% (Meyerholt et al. 2020). Furthermore, it is well known from the Neotropics that N limitation is also pervasive in many young secondary forests (Davidson et al. 2007;Batterman et al. 2013). As such, assessing shifts in nutrient cycling underlying carbon accumulation along secondary forest succession, is vital to understand potential thresholds in the recovery of these ecosystems (Townsend et al. 2011;Sullivan et al. 2019). In practice, studies on N recuperation in tropical regions are scarce, geographically biased and suggest an inconsistent timeframe of N cycle recuperation (Davidson et al. 2007;Winbourne et al. 2018). Most of the research has focused on Neotropical and Southeast Asian forests, while African forests received less research attention (Bauters et al. 2018(Bauters et al. , 2019Kukla et al. 2019;Wang et al. 2021). Furthermore, existing studies focus only on either soils or plants and rarely take a whole ecosystem approach (Davidson et al. 2007;Winbourne et al. 2018;Figueiredo et al. 2019;Sullivan et al. 2019). The balance between the N inputs, outputs and plant demand ultimately determines the forest N cycle (Vitousek et al. 1986) and earlier work has shown that the N cycling properties of forests can be recovered during secondary succession, i.e. that secondary forests progressively shift from a conservative N cycle in to a more open N cycle in old growth forest (Davidson et al. 2007;Sullivan et al. 2019).
Nitrogen enters natural terrestrial ecosystems via three main pathways: mineral weathering, biological N fixation (BNF), and atmospheric deposition. On highly weathered tropical soils, the contribution of N derived from mineral weathering is minor (Vitousek et al. 1986;Porder et al. 2005). BNF represents an important source of ecosystem N and was used as the dominant argument to explain N richness of tropical forests (Jenny 1950). However, BNF is an energetically expensive process, which can be alternatively activated or downregulated by N-fixing species in response to the ecosystem N status (Hedin et al. 2009;Barron et al. 2011). As such, BNF rates are generally high in the young successional stages where N demand is high relative to N availability (Barron et al. 2011;Batterman et al. 2013;Bauters et al. 2016). Forest canopies modify the chemical composition of rainfall in two main ways (De Schrijver et al. 2004). First, the forest canopy structure acts as a trap for atmospheric particles, in a way that canopies that are more complex can lead to higher dry deposition inputs to the forest floor (Bauters et al. 2021a). Second, canopy exchange adds (canopy leaching) or removes (canopy uptake) substantial nutrients to/from throughfall deposition that finally arrives on the forest floor (Parker 1983;Bortolazzi et al. 2021). Additionally, central African forest ecosystems are subjected to substantial N loads through atmospheric N deposition as NO 3 − (14.3 kg ha −1 yr −1 ), NH 4 + (12.2 kg ha −1 yr −1 ) and dissolved organic N (DON, 26.6 kg ha −1 yr −1 ) species (Bauters et al. 2019). DON represents an important constituent of total dissolved N (TDN) but is rarely quantified (Bauters et al. 2018(Bauters et al. , 2019Van Langenhove et al. 2020).
In the Congo basin, atmospheric N deposition is mostly due to biomass burning (Bauters et al. 2018). Nitrogen is removed from forest ecosystems through biomass export (e.g. logging), biomass burning, gaseous emissions (NOx, N 2 O and N 2 ), as well as dissolved inorganic N (DIN), organic N (DON) and particulate organic N by means of leaching and subsequent river export (Davidson and Martinelli 2009;Powers and Marín-Spiotta 2017). The soil microbial activity regulates-via mineralization and (de)nitrification-the amount of gaseous N emission and N bioavailability. Consequently, gross rates of mineralization and (de)nitrification can affect, at least indirectly, the rate of forest regrowth (Schmidt et al. 2007). Quantification of gaseous N emissions and leaching losses of bioavailable N species are commonly used to describe recuperation of open N cycles in late successional stages (Davidson et al. 2007;Hedin et al. 2009;Figueiredo et al. 2019).
In recent years, our understanding of N cycling in the central Congo basin has increased substantially (Bauters et al. 2016(Bauters et al. , 2018Baumgartner et al. 2020;Gallarotti et al. 2021). From these studies we know that the forests of the Congo basin receive high atmospheric N deposition (Bauters et al. 2018) but, at the same time, show relatively low rates of N 2 O emission (Gallarotti et al. 2021). Especially the high N deposition loads further question if the N cycle also evolves from a conservative to a more open cycle as forest succession proceeds, as it is the case in the Neotropics (Davidson et al. 2007;Figueiredo et al. 2019). To address this, we quantified N cycle metrics along a secondary successional gradient of Afrotropical lowland forests. Particularly, we aimed to quantify the effect of secondary succession stage on N deposition, N leaching, N emission and biomass N pools and -fluxes to identify recuperation from a conservative (in young secondary forests) to an open of N cycle (in old secondary forests). For this, we set up long-term monitoring sites, comprised of 18 plots in total, implemented in triplicate along six successional stages.

Study area
The study was conducted in post-agricultural forests at different successional stages within and around the Yoko Forest Reserve. The Yoko Forest Reserve (N00°17′; E25°18′; mean elevation 435 m a.s.l.) is located approximately 29 to 39 km southeast of Kisangani in DR Congo (Fig. 1). The vegetation of the region is semi-deciduous rainforest and the climate falls within the AF type in the classification of Köppen. The annual rainfall volume ranges between 1418 and 1915 mm and the annual mean temperature is 24.2 °C (De Ridder et al. 2014). We collected the samples in a setup of 18, 40 by 40 m plots, which were set up along a secondary succession gradient that consists of six successional stages, i.e., three plots per successional stage, beginning with agricultural fields (AG) and ending with an old growth forest (OG). The AG sites consisted of cassava fields that were prepared (e.g., the forest on the site was cleared and burned) in May 2018 and planted in July 2018. The secondary forests were aged 5, 12, 20 and 60-years-old. On basis of several field visits, interviews with landowners and local knowledge by farmers, plots were selected if (1) we were able to find at least three plots of the same age and land-use 1 3 Vol:. (1234567890) history (i.e. cut only once); and when (2) plots did not show signs of disturbances (e.g. trees harvested) since abandonment (Makelele et al. 2021). Soils in the region are highly weathered Ferrasols, being poor in nutrients, with low pH and dominated by a sandy texture (Table 1).
Nitrogen deposition and leaching: sampling protocol, sample and data analysis The bulk precipitation, throughfall and lysimeter samples were collected weekly during one hydrological year (i.e. from July 2018 to August 2019). For the throughfall and bulk deposition, the three agricultural plots (open field deposition) were used for bulk deposition estimates, and the twelve secondary forest plots for the respective throughfall estimates (i.e. 5yrs, 12yrs, 20yrs and 60yrs). Consequently, the bulk deposition was considered 'wet deposition', while the throughfall loads were considered 'wet + dry deposition + canopy exchange'. The OG plots were not sampled for N deposition and leaching because of logistical constraints but these plots were already studied earlier (see Bauters et al. 2018Bauters et al. , 2019. The bulk precipitation and throughfall samples were collected using polyethylene (PE) funnels, supported by a wooden pole of 1.5 m height that was connected to a 5-L PE container with a PE tube. A nylon mesh was placed in the neck of the funnels to avoid contamination of the collected water by large particles. The containers were buried and covered by litter to avoid the growth of algae and to keep the samples cooler. In total, we used a network of 144 collectors (i.e. 120 collectors for 15 plots of throughfall in the forest plots The study was conducted in and around the Yoko Forest Reserve (dark yellow circle in the inset) in the central Congo basin. AG, 5yrs, 12yrs, 20yrs, 60yrs and OG indicate the agricultural sites and 5, 12, 20, 60 years-old secondary forests and old growth forest plots, respectively. In addition, slash-and-burn events around our plots during the sampling period are shown 1 3 Vol.: (0123456789) and 24 collectors of 3 plots of bulk precipitation in the agricultural plots). The setup of bulk and throughfall collectors consisted of a rectangle of 8 by 24 m with 8 m distance between each collector. Given that we could not manage to have wet deposition funnels at our study site, we used the bulk deposition for wet deposition (Hofhansl et al. 2011;Van Langenhove et al. 2020).
Soil solution at 80 cm depth was sampled with 4 lysimeters per plot, installed in the vicinity of the throughfall and bulk precipitation collectors. These'suction cup' lysimeters consisted of a PVC tube fitted with a porous ceramic cup (Eijkelkamp Soil and Water, Giesbeek, The Netherlands) and connected to a buried opaque 2-L glass bottle by a PE tube. Using a battery-powered vacuum pump (Prenart Equipment, Copenhagen, Denmark), we applied an underpressure of -500 hPa at each sampling occasion.
For the throughfall, open field rainfall, and lysimeter sampling, the weekly sampling protocol was identical. As such, the water volume of each collector was measured in the field and the materials (i.e. recipients, funnels, mesh) were replaced and subsequently rinsed with distilled water. A plot-level composite sample of each collector category was made, considering the volume of each individual sample. These samples were transported to the lab in Kisangani, filtered with a 0.45 µm nylon membrane filter and put in the freezer for storage within 24 h. Finally, the samples were shipped in batch to Ghent University in Belgium for chemical analysis. The NH 4 + concentration was determined colorimetrically by the salycilate-nitroprusside method on an autoanalyzer (AA3; Bran and Luebbe, Norderstedt, Germany), while NO 3 − was measured by Ion Chromatography (Thermo-Scientific, Pittsburgh, Pennsylvania, USA). Due to some technical problems, NH 4 + could not be quantified for ca. 40% of the samples, and hence data from these samples was not included in the analysis. Total dissolved nitrogen (TDN) was determined in the water samples with the persulfate oxidation method. For this, an oxidizing solution containing NaOH, H 3 BO 3 and K 2 S 2 O 8 was added to the sample (Yasuitamura et al. 2020), which was subsequently placed in an autoclave at 121 °C for 1 h in order to convert NH 4 + and dissolved organic N (DON) into NO 3 − . We calculated the DON concentration as the difference between TDN and the sum of NH 4 + and NO 3 − analyzed in a non-digested sample.
The water flux for bulk precipitation and throughfall was calculated by dividing the average water volume by the surface area of the funnel. We calculated the intercept evaporation at plot scale as the difference between rainfall and throughfall under the assumption that the stem flow was negligible (Zimmermann et al. 2013). Nitrogen deposition fluxes were calculated by multiplying the water volumes with the corresponding concentration of N species. Further, we used the Table 1 Vegetation structure and soil characteristics of the plots at the study site. Plots were characterized regarding their stem density, basal area (BA), aboveground biomass (AGB) and soil texture, bulk density (BD), soil organic carbon (SOC) and pH are given are. AG, 5yrs, 12yrs, 20yrs, 60yrs and OG stand for the agricultural site, 5, 12, 20, 60 years-old secondary forest and old growth forest, respectively. Means, standard deviations (sd) and Kruskal-wallis P-values are reported *Significant p-value are highlighted in bold, the superscript letters represent significant difference between successional stages filtering approach of the canopy budget model to estimate dry deposition and canopy exchange from the throughfall and bulk precipitation data (Ulrich 1983;Van Langenhove et al. 2020). Consequently, the total atmospheric deposition was obtained by summing the bulk and the dry deposition, and the canopy exchange by subtracting total atmospheric deposition from throughfall deposition. The N leaching flux at 80 cm depth was calculated by the chloride mass balance (CMB) method (De Schrijver et al. 2004). This method is based on the assumption of mass conservation between the input of atmospheric chloride and the chloride flux in the subsoil (Eriksson and Khunakasem 1969).

Soil N metrics
We sampled topsoil (0-10 cm depth) once at the end of the dry season, at five locations and combined them into one composite sample per plot. The samples were oven dried at 60 °C for 48 h, homogenized and analyzed. The N content and δ 15 N was determined using an elemental analyzer (automated nitrogen carbon analyzer; ANCA-SL, SerCon, UK), interfaced with an isotope ratios mass spectrometer (IRMS; 20-22, SerCon, UK). The soil texture was determined by the pipet method (Burt 2004). After removal of the organic matter with H 2 O 2 30%, soil particles were dispersed by Na 4 P 2 O 7 (aq). The total soil N stock was calculated as the product of soil N content, bulk density, and the soil upper 10 cm, and extrapolated to the plot-scale surface area. We extracted NH 4 + and NO 3 − with 1 M KCl. For each sample, 30 g soil was extracted with 60 ml 1 M KCl, shaken for 1 h and filtered (MN615; Macherey-Nagel, Darmstadt, Germany). The NH 4 + and NO 3 − concentrations in the KCl extracts were determined colorimetrically using an auto analyzer (AA3; Bran and Luebbe, Norderstedt, Germany).
Nitrous oxide (N 2 O) emissions were determined every two weeks during one hydrological year (e.g. from April 2019 to May 2020 because of delayed arrival of static chambers at the study site) using manual static chambers (Hutchinson and Mosier 1981). The chamber bases were installed permanently to avoid soil disturbance. We used chambers made of PVC (diameter = 0.3 m, height = 0.3 m), equipped with an airtight lid, thermocouples, sampling ports, and a vent tube to avoid pressure disturbances. On every sampling occasion, chambers were closed for 1 h and the headspace samples taken at 20 min interval using a 20 mL syringe and stored in pre-evacuated 12 mL containers (Exetainer, Labco, Lampeter, UK). At each sampling time point, temperature inside the chamber was measured using a thermocouple (Type T, Omega Engineering, Stanford, CT, USA). The samples were analyzed at ETH Zurich for N 2 O using gas chromatography (456-GC; Scion Instruments, Livingston, UK). We calculated the N 2 O fluxes according to Eq. 1.
where V is the volume of static chamber, P is the pressure, R is the gas constant (0.08206 [L atm K −1 mol −1 ]), S is the area, T is the temperature and ΔC/Δt is the slope of linear regression model representing the rate of concentration change in one hour (Hutchinson and Mosier 1981).

Plant N metrics
A tree species inventory was conducted in the 40 by 40 m plots. The diameter at breast height (DBH: 1.3 m along the stem), or 30 cm above any deformities in case it occurred at that point of measurement, of trees with diameter ≥ 10 cm was measured, marked, and identified at species level. Fully expanded, sun leaves were sampled from the canopy for all species that contributed to 85% of the cumulative plot basal area. For each selected species at least three individuals were sampled and for each individual ten leaves were collected and measured (Pérez-Harguindeguy et al. 2013), more details on leaf trait sampling and measurement are available in Makelele et al. (2021). The leaves were oven-dried at 70 °C for 48 h and ground before chemical analyses. Leaf C and N content and leaf δ 15 N were analyzed using an elemental analyzer (ANCA-SL, Ser-Con, UK), interfaced with an isotope ratio mass spectrometer (20-22, SerCon, UK). We further calculated the community-weighted mean (CWM) per plot for these measured leaf functional traits (N, C:N ratios) using the species composition and average species leaf traits, weighing with basal area. Litter traps (diameter = 0.68 m, height = 1.5 m, mesh size = 1 mm) were installed in parallel to the throughfall funnels. Consequently, as for throughfall, eight litter traps per plots were set up in a rectangular 8 by 24 m sub-plot and were sampled once a week from July 2018 to August 2019 in twelve plots of secondary forests (5yrs, 12yrs, 20yrs and 60yrs). As we were interested in the fine litterfall, branches of diameter > 2 cm were discarded when present in our traps (Marthews et al. 2014). Subsequently, litter samples were oven dried at 70 °C for 24 h, weighed and ground before chemical analysis. A subset of litter samples, selected to be representative of seasons, was analyzed for N content using an elemental analyzer (ANCA-SL, SerCon, UK). The litter N flux was calculated by multiplying litter mass by litter N content.
We used ingrowth cores to assess fine roots production from December 2019 to February 2021. The sampling was conducted with 3 months intervals according to an established protocol (Metcalfe et al. 2007;Marthews et al. 2014). We washed the fine roots, oven dried them for 24 h at 70 °C, and finally weighted and ground them before chemical analysis. The amount of root dry mass collected after three months represents the quarterly root production and subsequently the annual production is the sum of the quarterly production. A subset of root samples, selected to be representative of seasons, was analyzed for N content using an elemental analyzer (ANCA-SL, SerCon, UK). The contribution of fine root production to the plot N was calculated by multiplying the fine root mass production by the corresponding fine root N content.

Statistical analysis
We conducted one-way analysis of variance (ANOVA) to test the variation of the different components of the N cycle and metrics between successional stages. Additionally, we applied the Tuckey HSD post-hoc test. In case the underlying assumptions of ANOVA, mainly homoscedasticity and normality, were not met, we used the Kruskal-Wallis and subsequently its post-hoc Dunn's test to analyze the data. Specifically, Kruskal-Wallis and Dunn's test were performed to compare between successional stages the flux of NO 3 − and NH 4 + in both the leaching and budget, N 2 O flux and leaf C:N. We performed the statistical analyses using the R programming language 3.6.1 (R Development core 2018).

Nitrogen input and output fluxes
The overall amount of N that reached the forest floor via throughfall was affected by the forest stand age via the canopy complexity and composition (Fig. 2, Table S1). During the first twenty years of succession the forests show negative values of canopy exchange that is indicative of net canopy N uptake and conservation (Table S1; Bortolazzi et al. 2021). The net canopy N uptake varied between 2.5% and 44.1% of total atmospheric N deposition. Secondly, young secondary forests (up to 20 years) showed a reduced canopy efficiency and high plot to plot variation in trapping dry N particles (Table S1). The 60 years secondary forests showed a high efficiency for trapping atmospheric dry N particles and positive canopy exchange values indicative of net canopy leaching (Table S1). Consequently, the amount of N reaching the forest floor (throughfall) was significantly higher in the 60 years secondary forests (P = 0.044, P = 0.026, 1 3 Vol:. (1234567890) P = 0.005, P = 0.004 for ammonium, nitrate, DON and TDN, respectively) than in the young forests. Wet and dry deposition contributed proportionally to the throughfall deposition in the 60 years old forest (being respectively 39% and 36%) and the canopy leaching contributed 25%.
In general, we observed a decrease in the leaching fluxes from AG to the 20 years-old secondary forest (Fig. 2, Table S1) and a re-increase above 20 years of succession. In the 60 years secondary forests, the leached NO 3 − showed variable δ 15 N and a low δ 18 O signature, while δ 18 O was higher for the other successional stages ( Figure S3). The partial atmospheric deposition -leaching balance was negative in the AG sites and approached zero after the first decade of succession (Fig. 2). Moreover, N concentration in the leachate of the agricultural sites showed a consistent decreasing trend over the time and were considerably higher during the first half of the sampling period ( Figure S2). Additionally, N 2 O fluxes increased significantly along the successional stages (P = 0.0062). The agricultural field and 60 years-old forest showed the lowest and highest values, respectively (Fig. 3). The total soil N and the available N pool (NH 4 + and NO 3 − ) and soil δ 15 N did not show a consistent change with forest succession status ( Figure S4). ), dissolved organic carbon (DON) and total dissolved nitrogen (TDN) input, output and partial deposition-leaching balance along the studied secondary succession gradient; a bulk (i.e. wet) and total atmospheric (i.e. wet + dry) deposition (kg N ha −1 yr −1 ), b throughfall deposition (wet + dry + canopy leaching) deposition (kg N ha −1 yr −1 ), c leaching loss at 80 cm depth (kg N ha −1 yr −1 ), d partial deposition leaching budget (kg N ha −1 yr −1 ). AG, 5yrs, 12yrs, 20yrs and 60yrs indicate the agricultural sites and 5, 12, 20 and 60 years-old secondary forests. The bars and diamond shaped points represent successional stages mean; circle points are plot level data and error bars are standard deviations 1 3 Vol.: (0123456789) Tissue N The weighted foliar N and C:N ratio increased and decreased, respectively, along the successional gradient and thus showed higher and lower values in the late successional stages (Fig. 3, Table S4). Weighted foliar δ 15 N and foliar N content were maximal in the late successional stages (i.e., 60 years-old secondary forest and old growth forest), though the foliar δ 15 N also showed a decreasing trend early on the succession (Fig. 3). Overall, the N recycled internally through leaf litter and fine root fluxes did not significantly change with successional stage. However, the amount of N transferred by litter deposition increased with stand age and ranged between 119 ± 45 kg N ha −1 yr −1 and 178 ± 51 kg N ha −1 yr −1 (Fig. 3). The litter N concentration showed an increasing trend over time with the highest value in the late successional stages. The litter C:N ratios were lowest in the 60 years-old growth forests relatively to the young secondary forests. The N amount assimilated by the fine roots was ranged between 47 ± 5 kg N ha −1 yr −1 and 139 ± 47 kg N ha −1 yr −1 (Fig. 3). Fine root N concentration (%) and fine root C:N ratio showed higher and lower values in the late successional stages, respectively.

Hydrologic balance
Structural vegetation parameters, such as basal area and canopy height, consistently changed along the successional gradient at our study site (Makelele et al. 2021). But an effect of forest structure recovery on water fluxes was only observed in the first decade of succession. Interception was highest in the 5 yearsold forests, while leaching was maximal in the agricultural sites and minimal in the 5 and 12 years-old forests ( Figure S1). The relative interception values are between the range of values found previously at our study site (i.e. 6% and 17%, Bauters et al. 2019), but are relatively low compared to the global forest interception average (Magliano et al. 2019;Yue et al. 2021). Forest composition and structure or individual species' crown morphology plays an important role in water interception (Yang et al. 2019;Yue et al. 2021). As such, high interception earlier in the succession could be attributable to the crown shape of Musanga cecropioides, the umbrella tree, dominating young secondary forest of the Congo basin and closing the canopy within the first years of forest regrowth (Ross 1954). On the other hand, Zimmermann et al. (2013) attributed recovery of forest interception to vegetation structure restoration in the Brazilian Amazon. However, contrasting results have been found in Mexican mountain forests where interception after 20 years forest regrowth was 50% of that of mature forest due to slow recovery of leaf area index and epiphytic biomass (Muñoz-Villers et al. 2012). High leaching from the agricultural site is the result of reduced vegetation interception and evaporation (Bruijnzeel 2004;Zimmermann et al. 2013). The reduced leaching in the 5-year-old forest may result from increased interception and legacy of soil impermeability from the previous agriculture activities (Bruijnzeel 2004). Globally, the effect of forest disturbance on the hydrological cycle depends on site, hydro-climate and landscape characteristics (Sandström 1995). Here, we show that the effect of shifting cultivation on the forest ecosystem hydrological process is only observed during the first years of succession (see also , Zimmermann et al. 2013).

Fig. 3
Soil N 2 O emission (kg N ha −1 yr −1 ) along secondary forest succession in the Congo basin. AG, 5yrs, 12yrs, 20yrs, 60yrs and OG stand respectively for the agricultural sites, 5, 12, 20, 60 years-old secondary forests and old growth forest. The bars and diamond shaped points represent successional stages mean; circle points are plot level data and error bars are standard deviations Bulk and throughfall N deposition One of the major environmental consequences of N pollution on the planet is the increase in atmospheric N deposition (Stevens 2019). Recent research has shown that the Congo basin is a hotspot for N deposition (Bauters et al. 2018(Bauters et al. , 2019, even more than previously predicted (Reay et al. 2008). Here, we confirm that the Congo basin is subjected to a high bulk N deposition (i.e. open field deposition) of 11.6 ± 2.6 kg N ha −1 yr −1 , lower but comparable to 18.2 kg N ha −1 yr −1 that has been previously found (Bauters et al. 2018). We also corroborated that DON is contributing roughly 50% to TDN deposition. The amount of bulk N deposition for the central Congo basin is substantially higher than what has been found previously in Central and South America, Southeast Asia and West Africa (Williams et al. 1997;Oziegbe et al. 2011;Kohler et al. 2013;Van Langenhove et al. 2020). For example, for a study using the same field protocol as this study conducted in French Guiana, the amount of TDN and inorganic N (NH 4 + and NO 3 − ) was, respectively, two and four times lower than what we found (Van Langenhove et al. 2020). However, the relative proportion of DON measured in French Guiana was in the same range, indicating that omitting DON quantification underestimates the magnitude of N deposition in tropical forests. In contrast to industrial and emerging countries, where elevated N deposition is derived from industrial and intensive agriculture activities, elevated N deposition at our study site is sourced from biomass burning (Chen et al. 2010;Bauters et al. 2019). Chen et al. (2010) estimated that between 50 and 80% of the N deposited in the Congo basin, South American and Southeast Asian tropical forest is fire-derived. Indeed, Bauters et al. (2018) showed that TDN deposited on forests of the Congo basin was linked to the biomass burning intensities in the prevailing wind directions over the seasons. Further, they showed a positive correlation between DON deposition and the concentration of condensed aromatics in rainfall, which are organic compounds that are exclusively sourced by combustion.
Furthermore, in this study, we show how the canopy composition and complexity, through dry particle interception and canopy exchange, sensibly modify the amount of N that reaches the forest floor by means of throughfall. Vegetation structure attributes, such as canopy height and density, play a great role in trapping atmospheric particles (Lawrence et al. 2008;Runyan et al. 2013;Powers and Marín-Spiotta 2017). Indeed, more complex canopies have been shown more efficient in trapping atmospheric particles and nutrients (Runyan et al. 2013;Bauters et al. 2021a) such that the amount of nutrients deposited increases with the canopy complexity along secondary succession (Bauters et al. 2021a). Consequently, the change of forest composition and structure as the succession proceed (Bauters et al. 2021a;Makelele et al. 2021) results in increasing total N deposition as shown here (i.e. between 12 and 60 years-old forests). The amount of N deposited locally can be highly variable and in some cases dependent on the distance between the site and an active slash-and-burn field and on the landscape configuration (Kohler et al. 2013). Kohler et al (2013) showed that plots closed to an active slash-and-burn field receive more atmospheric N as N-NO 3 − in comparison to plots located deep in the forest and relatively far from slash-and-burn activities. The distance between each individual plot and an active slash-and-burn field (Fig. 1) is potentially responsible for the high variability of atmospheric N deposition in young secondary forest (Table S1) and the unexpected high N deposition in the 5 years secondary forest. On the other hand, our results suggest that young secondary forest substantially rely on canopy N uptake, likely to compensate N limitation. Indeed, the canopy may also retain a substantial amount of N via leaves stomata and cuticle absorption (Bortolazzi et al. 2021). In contrast, canopy N leaching in the 60 years secondary forest may be indicative of N saturation of aged tropical forests. Altogether, the forests of the Congo basin are receiving high amounts of N from biomass burning occurring in shifting cultivation systems in the forest and savannah biomes of the Congo basin. However, the forest's successional status affects the amount of N that reaches the forest floor, through canopy trapping and exchange. On this basis, we postulate that the main external sources of N for tropical forest of the Congo basin shift from BNF in the youngest forests (Bauters et al. 2016) to atmospheric N deposition in the older forests. As such the canopy composition and complexity play an important role in the recuperation of the N cycle after a shifting cultivation event.

Soil and leaching N fluxes
The δ 15 N and available and total soil N stocks (Figures S4 and S5) did not change with the forest successional status. On the one hand, the high δ 15 N value (approx. 7.5‰, Figure S5) is indicative of a system with long-term open N cycling (i.e. relatively high inputs and outputs that results in isotop fractionation over the long term), which includes most humid tropical forests globally (Martinelli et al. 1999;Houlton and Bai 2009;Craine et al. 2015). On the other hand, our finding corroborates with a number of studies from the Neotropics and South East Asia (e.g. Markewitz et al. 2004;Winbourne et al. 2018;Kukla et al. 2019). Relatively extensive agricultural activity such as the shifting cultivation at our study site tends to preserve nitrogen (Bauters et al. 2021b). The conservation of substantial N in the soil should facilitate rapid recuperation of many processes of the N cycle (Davidson et al. 2007;Winbourne et al. 2018). However, the projected population growth in the Democratic Republic of the Congo will likely result in a more intense (i.e. reduction of fallow length) shifting cultivation. The reduction of fallow lengths and the repeated slash-and-burn cycles deplete soil N and other nutrients and might ultimately affect the recuperation trajectory of N and vegetation in the traditional shifting cultivation systems (Zarin et al. 2005;Moonen et al. 2019;Bauters et al. 2021b).
The N leaching fluxes were strikingly high in the agricultural sites and mostly composed of NO 3 − . The leached NO 3 − showed elevated δ 18 O values for all the successional stages (comparable to that of NO 3 − in rainfall; Pardo et al. 2004), except for the 60 years-old secondary forests ( Figure S3). This imply that the leached NO 3 − is at least partly composed of NO 3 − directly derived from atmospheric deposition. The NO3-leaching from the 60-yearold-forest exhibited consistently depleted δ 18 O-signature compared to the younger forest types. Given that these lower d 18 O values are indicative for nitrification-derived NO 3 − (Pardo et al. 2004), we conclude that the in-situ mineralization and nitrification in the 60-year-old forest is increasingly contributing to the overall leaching flux. The effect of fire on N leaching in a shifting cultivation system is sensibly detectable in the six first months ( Figure S2) that follow field burning (see also, da Silva Neto et al. 2019 for the effect on soil properties). The high N leaching flux during the agricultural phase is thus most likely a combined result of: (1) enhanced soil water recharge ( Figure S1, see also Owuor et al. 2016); (2) preceding burning, which increases soil available nutrients by incorporating nutrient rich-ash from aboveground biomass burning in the soil along with the chemical alteration of soil organic matter and biomass (Murphy et al. 2006;Kukla et al. 2019), and (3) reduced vegetation N uptake, especially before crop planting. Further, we showed that the amount of N (24.5 ± 29 kg N ha −1 yr −1 ) leached from the 60-years-old was higher compared to the findings of Bauters et al. (2019) in mature mixed and monodominant lowland forests at our study site, where they reported soil-leaching fluxes of 11.5 ± 2.8 and 15.0 ± 5.5 kg N ha −1 yr −1 , respectively.

N cycle recuperation
A recuperation of the N cycle during secondary succession is reflected in most analyzed plant tissuebased N metrics. We showed an increasing tissue N content and decreasing tissue C:N ratios along the succession. The latter is consistent with increasing N availability (Davidson et al. 2007;Hietz et al. 2011). Foliar δ 15 N values of the central Congo basin (3.1‰ ± 0.8‰ to 6.8‰ ± 0.3‰) fall within the range of values for the Amazon forest (0.9‰ to 10.0‰, Ometto et al. 2006;Davidson et al. 2007). The canopy δ 15 N signal translates ecosystem N processes and N input-output budgets (Compton et al. 2007). Conceptually, increasing N availability during secondary succession due to BNF and atmospheric N deposition results in loss of δ 15 N depleted N species and thus δ 15 N enrichment (Amundson et al. 2003). Our foliar δ 15 N data show two trends in the recuperation process of the N cycle of the central Congo basin: (1) initial higher foliar δ 15 N as a result from biomass burning is reduced during the first 20 years of succession due to 15 N depletion resulting from BNF -a trend corroborated by several earlier studies (Compton et al. 2007;Wang et al. 2007;Broadbent et al. 2014), and (2) after 20 years into the succession, we show a reincrease of foliar δ 15 N likely due to a declining BNF (Bauters et al. 2016) gradual increasing N availability, leading to increasing gaseous (Fig. 4) and hydrological loss in the late successional stage, which enriches the remaining plant available N pool in 15 N (Amundson et al. 2003;Davidson et al. 2007 In line with this trend in foliar δ 15 N, we showed a decreasing N leaching in the young forest and a reincrease during the late successional stages (Fig. 4). In contrast to the findings of Bauters et al. (2019) in a mature mixed forest, we showed that mineral N (mostly NO 3 − ) was the dominant leaching N species, which is in accordance to a more open N cycle in the late successional stages of tropical forests (Hedin et al. 2009;Elser 2011;Brookshire et al. 2012). In addition, we showed that all the deposited N in the late successional stage is fully recycled in the system before leaching, while the leaching N in the younger forests are composed of a mix between the atmospheric deposition and soil N ( Figure S3).
The transition of a more conservative N cycle in the young secondary tropical forest to a more open N cycle in late successional tropical forest is associated with increasing N 2 O emissions (Davidson et al. 2007). We showed increasing soil emissions of N 2 O along the successional stages. These trends corroborate with findings from similar studies from the Neotropics (Davidson et al. 2007) and South East Asia (Wang et al. 2021) with N 2 O emissions increasing from 3 to 70 years and 12 to 54 years of succession, respectively. Low N 2 O emissions in a young secondary forest results from high N losses during the agricultural phase and is thus generally used as an indicator of a conservative N cycle of younger forests. The N 2 O emissions in the 60 years-old secondary forest was double the amount emitted in the lowland OG forests (Fig. 4, see also Gallarotti et al. 2021;Barthel et al. 2022) and above the global tropical forest N 2 O emissions mean (Werner et al. 2007), but it was similar to the emissions in a secondary forest in Cameroun (Verchot et al. 2020) and OG forests in Ghana (Castaldi et al. 2013) and Congo (Serca et al. 1994). The N 2 O emission peak in the 60 years-old secondary forest additionally indicates that the N cycle is recovering to its status in a mature forest, following shifting cultivation (Winbourne et al. 2018).
Overall, a recuperation of a typically open N cycle of a mature lowland tropical forest is, in many cases, Fig. 4 Variation of vegetation tissue N metrics along secondary succession including foliar, litter and root C:N ratios; foliar, litter and root N concentration (%), litter and root internal N recycling (kg N ha −1 yr −1 ) and foliar δ 15 N (‰); 5yrs, 12yrs, 20yrs, 60yrs and OG stand respectively for 5, 12, 20 and 60 years-old secondary forests and the old growth forest. The bars and diamond shaped points represent successional stages mean; circle points are plot level data and error bars are standard deviations associated with increasing N 2 O emissions, increasing foliar δ 15 N as well as foliar, litter and root N content and decreasing foliar, root and litter C:N ratios, in line with what earlier studies reported (Davidson et al. 2007;Sullivan et al. 2019). The N status of the 60 years-old forests was similar to that of old growth forests suggesting that the N cycle at our study site recovered between 20 and 60 years. In conclusion, and in accordance with the points developed above, our data depict such a recuperation process despite elevated atmospheric N deposition. As such, this also suggests that N deposition does not alleviate N limitation early-on in secondary successions in Central Africa.