• Energy Research
• natural sciences close

Abstract. 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 NO3- 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 NO3- were on average 27 % (range 6 %–54 %) of Ndep at sites with Ndep < 1 g N m−2 yr−1 versus 65 % (range 35 %–85 %) for 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.

Severe drought events are known to cause important reductions of gross primary productivity ( GPP ) in forest ecosystems. However, it is still unclear whether this reduction originates from stomatal closure (Stomatal Origin Limitation) and/or non-stomatal limitations (Non-SOL). In this study, we investigated the impact of edaphic drought in 2018 on GPP and its origin (SOL, NSOL) using a dataset of 10 European forest ecosystem flux towers. In all stations where GPP reductions were observed during the drought, these were largely explained by declines in the maximum apparent canopy scale carboxylation rate V CMAX,APP (NSOL) when the soil relative extractable water content dropped below around 0.4. Concurrently, we found that the stomatal slope parameter ( G 1 , related to SOL) of the Medlyn et al . unified optimization model linking vegetation conductance and GPP remained relatively constant. These results strengthen the increasing evidence that NSOL should be included in stomatal conductance/photosynthesis models to faithfully simulate both GPP and water fluxes in forest ecosystems during severe drought. This article is part of the theme issue ‘Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale’.

AbstractThe European CARBOEUROPE/FLUXNET monitoring sites, spatial remote sensing observations via the EOS‐MODIS sensor and ecosystem modelling provide independent and complementary views on the effect of the 2003 heatwave on the European biosphere's productivity and carbon balance. In our analysis, these data streams consistently demonstrate a strong negative anomaly of the primary productivity during the summer of 2003. FLUXNET eddy‐covariance data indicate that the drop in productivity was not primarily caused by high temperatures (‘heat stress’) but rather by limitation of water (drought stress) and that, contrary to the classical expectation about a heat wave, not only gross primary productivity but also ecosystem respiration declined by up to more than to 80 gC m−2 month−1. Anomalies of carbon and water fluxes were strongly correlated. While there are large between‐site differences in water‐use efficiency (WUE, 1–6 kg C kg−1 H2O) here defined as gross carbon uptake divided by evapotranspiration (WUE=GPP/ET), the year‐to‐year changes in WUE were small (<1 g kg−1) and quite similar for most sites (i.e. WUE decreased during the year of the heatwave). Remote sensing data from MODIS and AVHRR both indicate a strong negative anomaly of the fraction of absorbed photosynthetically active radiation in summer 2003, at more than five standard deviations of the previous years. The spatial differentiation of this anomaly follows climatic and land‐use patterns: Largest anomalies occur in the centre of the meteorological anomaly (central Western Europe) and in areas dominated by crops or grassland. A preliminary model intercomparison along a gradient from data‐oriented models to process‐oriented models indicates that all approaches are similarly describing the spatial pattern of ecosystem sensitivity to the climatic 2003 event with major exceptions in the Alps and parts of Eastern Europe, but differed with respect to their interannual variability.