Search

The goal of sustainable coffee production requires multiple functions from agroforestry systems. Many are difficult to quantify and data are lacking, hampering the choice of shade tree species and agronomic management. Process-based modelling may help quantify ecosystem services and disservices. We introduce and apply coffee agroforestry model CAF2021 (https://doi.org/10.5281/zenodo.5725277). The model allows for complex systems with up to three shade tree species. It simulates coffee yield, timber and fruit production by shade trees, soil loss in erosion, C-sequestration, N-fixation, -emission and -leaching. To calibrate the model, we used multivariate data from 32 different treatments applied in two long-term coffee agroforestry experiments in Costa Rica and Nicaragua. Without any further calibration, the model was then applied to agroforestry systems on 89 farms in Costa Rica and 79 in Guatemala where yields had been reported previously in farmer interviews. Despite wide variation in environmental and agronomic conditions, the model explained 36% of yield variation in Costa Rica but only 15% in Guatemala. Model analysis quantified trade-offs between yield and other ecosystem services as a function of fertilisation and shading.

The impact of atmospheric reactive nitrogen (Nr) deposition on carbon (C) sequestration in soils and biomass of unfertilized, natural, semi-natural and forest ecosystems has been much debated. Many previous results of this dC∕dN response were based on changes in carbon stocks from periodical soil and ecosystem inventories, associated with estimates of Nr deposition obtained from large-scale chemical transport models. This study and a companion paper (Flechard et al., 2020) strive to reduce uncertainties of N effects on C sequestration by linking multi-annual gross and net ecosystem productivity estimates from 40 eddy covariance flux towers across Europe to local measurement-based estimates of dry and wet Nr deposition from a dedicated collocated monitoring network. To identify possible ecological drivers and processes affecting the interplay between C and Nr inputs and losses, these data were also combined with in situ flux measurements of NO, N2O and CH4 fluxes; soil NO−3 leaching sampling; and results of soil incubation experiments for N and greenhouse gas (GHG) emissions, as well as surveys of available data from online databases and from the literature, together with forest ecosystem (BASFOR) modelling. Multi-year averages of net ecosystem productivity (NEP) in forests ranged from −70 to 826 g C m−2 yr−1 at total wet + dry inorganic Nr deposition rates (Ndep) of 0.3 to 4.3 g N m−2 yr−1 and from −4 to 361 g C m−2 yr−1 at Ndep rates of 0.1 to 3.1 g N m−2 yr−1 in short semi-natural vegetation (moorlands, wetlands and unfertilized extensively managed grasslands). The GHG budgets of the forests were strongly dominated by CO2 exchange, while CH4 and N2O exchange comprised a larger proportion of the GHG balance in short semi-natural vegetation. Uncertainties in elemental budgets were much larger for nitrogen than carbon, especially at sites with elevated Ndep where Nr leaching losses were also very large, and compounded by the lack of reliable data on organic nitrogen and N2 losses by denitrification. Nitrogen losses in the form of NO, N2O and especially NO−3 were on average 27 % (range 6 %–54 %) of Ndep at sites with Ndep 3 g N m−2 yr−1. Such large levels of Nr loss likely indicate that different stages of N saturation occurred at a number of sites. The joint analysis of the C and N budgets provided further hints that N saturation could be detected in altered patterns of forest growth. Net ecosystem productivity increased with Nr deposition up to 2–2.5 g N m−2 yr−1, with large scatter associated with a wide range in carbon sequestration efficiency (CSE, defined as the NEP ∕ GPP ratio). At elevated Ndep levels (> 2.5 g N m−2 yr−1), where inorganic Nr losses were also increasingly large, NEP levelled off and then decreased. The apparent increase in NEP at low to intermediate Ndep levels was partly the result of geographical cross-correlations between Ndep and climate, indicating that the actual mean dC∕dN response at individual sites was significantly lower than would be suggested by a simple, straightforward regression of NEP vs. Ndep. Financial support. This research has been supported by the Acknowledgements. The authors gratefully acknowledge financial support by the European Commission through the two FP6 integrated projects CarboEurope Integrated Project (project no. GOCE-CT-2003-505572) and NitroEurope Integrated Project (project no. 017841), the FP7 ECLAIRE project (grant agreement no. 282910), and the ABBA COST Action ES0804. We are also thankful for funding from the French GIP-ECOFOR consortium under the F-ORE-T forest observation and experimentation network, as well as from the MDM-2017-0714 Spanish grant. We are grateful to Christian Bernhofer, Robert Clement, Han Dolman, Axel Don, Eric Dufrêne, Damiano Gianelle, Ruediger Grote, Anders Lindroth, John Moncrieff, Dario Papale, Corinna Rebmann and Alex Vermeulen for the data they provided, as well as to Klau-dia Ziemblińska for her comments on the paper. Computer time for EMEP model runs was supported by the Research Council of Norway through the NOTUR project EMEP (NN2890K). Finalization of the paper was supported by the UK Natural Environment Research Council award number NE/R016429/1 as part of the UK-SCAPE programme delivering national capability. We also wish to thank two anonymous referees for their constructive criticism of the paper.