NEE estimates 2006–2019 over Europe from a pre-operational ensemble-inversion system

Abstract. 3-hourly Net Ecosystem Exchange (NEE) is estimated at spatial scales of 0.25 degrees over the European continent, based on the pre-operational inverse modelling framework CarboScope Regional (CSR) for the years 2006 to 2019. To assess the uncertainty originating from the choice of a-priori flux models and observational data, ensembles of inversions were produced using three terrestrial ecosystem flux models, two ocean flux models, and three sets of atmospheric stations. We find that the station set ensemble accounts for 61 % of the total spread of the annually aggregated fluxes over the full domain when varying all these elements, while the biosphere and ocean ensembles resulted in much smaller contributions to the spread of 28 % and 11 %, respectively. These percentages differ over the specific regions of Europe, based on the availability of atmospheric data. For example, the spread of the biosphere ensemble is prone to be larger in regions that are less constrained by CO₂ measurements. We further investigate the unprecedented increase of temperature and simultaneous reduction of Soil Water Content (SWC) observed in 2018 and 2019. We find that NEE estimates during these two years suggest an impact of drought occurrences represented by the reduction of Net Primary Productivity (NPP), which in turn lead to less CO₂ uptake across Europe in 2018 and 2019, resulting in anomalies up to 0.13 and 0.07 PgC yr^-1 above the climatological mean, respectively. Annual temperature anomalies also exceeded the climatological mean by 0.46 °C in 2018 and by 0.56 °C in 2019, while standardized-precipitation-evaporation-index (SPEI) anomalies declined to −0.20 and −0.05 SPEI units below the climatological mean in both 2018 and 2019, respectively. Therefore, the biogenic fluxes showed a weaker sink of CO₂ in both 2018 and 2019 (−0.22±0.05 and −0.28±0.06 PgC yr^-1, respectively) in comparison with the mean −0.36±0.07 PgC yr^-1 calculated over the full analysed period (i.e., fourteen years). These translate into a continental-wide reduction of the annual sink by 39 % and 22 %, respectively, larger than the typical year-to-year standard deviation of 19 % observed over the full period.