ALBERT

Description: In the last decades terrestrial ecosystems have reabsorbed on average more than one-quarter of anthropogenic emissions (Le Quéré et al., 2018). However, this large carbon sink is modulated by climate and is therefore highly variable in time and space. The magnitude and temporal changes in the sensitivity of terrestrial CO2 fluxes to climate drivers are key factors to determine future atmospheric CO2 concentration and climate trajectories. In the literature, there is so far a strong focus on the climatic controls of daily and long-term variability, while less is known about the key drivers at a seasonal timescale and about their variation over time (Wohlfahrt et al., 2008). This latter temporal scale is relevant to assess which climatic drivers dominate the seasonality of the fluxes and to understand which factors limit the CO2 exchange during the course of the year. Here, we investigate the global sensitivity of net terrestrial CO2 fluxes, derived from atmospheric inversion, to three key climate drivers (i.e. global radiation and temperature from WFDEI and soil water content from ERA-Interim) from weekly to seasonal temporal scales, in order to explore the short-term interdependence between climate and the terrestrial carbon budget. We observed that the CO2 exchange is controlled by temperature during the carbon uptake period over most of the land surface (from 55 % to 52 % of the total surface), while radiation is the most widespread dominant climate driver during the carbon release period (from 64 % to 70 % of the total surface). As expected, soil water content plays a key role in arid regions of the Southern Hemisphere during both the carbon uptake and the carbon release period. Looking at the decadal trend of these sensitivities (1985–2016) we observed that the importance of radiation as a driver is increasing over time, while we observed a decrease in sensitivity to temperature in Eurasia. Overall, we show that flux temporal variation due to a specific driver has been dominated by the temporal changes in ecosystem sensitivity (i.e. the response of ecosystem to climate) rather than to the temporal variability of the climate driver itself over the last decades. Ultimately, this analysis shows that the ecosystem response to climate is significantly changing both in space and in time, with potential repercussion on the future terrestrial CO2 sink and therefore on the role that land may play in climate trajectories.

Description: Wetlands are a major source of methane (CH4) and contribute between 30 and 40% to the total CH4 emissions. Wetland CH4 emissions depend on temperature, water table depth, and both the quantity and quality of organic matter. Global warming will affect these three drivers of methanogenesis, raising questions about the feedbacks between natural methane production and climate change. Until present the large-scale response of wetland CH4 emissions to climate has been investigated with land-surface models that have produced contrasting results. Here, we produce a novel global estimate of wetland methane emissions based on atmospheric inverse modeling of CH4 fluxes and observed temperature and precipitation. Our data-driven model suggests that by 2100, current emissions may increase by 50% to 80%, which is within the range of 50% and 150% reported in previous studies. This finding highlights the importance of limiting global warming below 2°C to avoid substantial climate feedbacks driven by methane emissions from natural wetlands.

Description: Sun-induced chlorophyll fluorescence (SIF) retrieved from satellite spectrometers can be a highly valuable proxy for photosynthesis. The SIF signal is very small and notoriously difficult to measure, requiring sub-nanometre spectral-resolution measurements, which to date are only available from atmospheric spectrometers sampling at low spatial resolution. For example, the widely used SIF dataset derived from the GOME-2 mission is typically provided in 0.5∘ composites. This paper presents a new SIF dataset based on GOME-2 satellite observations with an enhanced spatial resolution of 0.05∘ and an 8 d time step covering the period 2007–2018. It leverages on a proven methodology that relies on using a light-use efficiency (LUE) modelling approach to establish a semi-empirical relationship between SIF and various explanatory variables derived from remote sensing at higher spatial resolution. An optimal set of explanatory variables is selected based on an independent validation with OCO-2 SIF observations, which are only sparsely available but have a high accuracy and spatial resolution. After bias correction, the resulting downscaled SIF data show high spatio-temporal agreement with the first SIF retrievals from the new TROPOMI mission, opening the path towards establishing a surrogate archive for this promising new dataset. We foresee this new SIF dataset becoming a valuable asset for Earth system science in general and for monitoring vegetation productivity in particular. The dataset is available at https://doi.org/10.2905/21935FFC-B797-4BEE-94DA-8FEC85B3F9E1 (Duveiller et al., 2019).

Description: Clouds play an important role in the climate system: (1) cooling Earth by reflecting incoming sunlight to space and (2) warming Earth by reducing thermal energy loss to space. Cloud radiative effects are especially important in polar regions and have the potential to significantly alter the impact of sea ice decline on the surface radiation budget. Using CERES (Clouds and the Earth's Radiant Energy System) data and 32 CMIP5 (Coupled Model Intercomparison Project) climate models, we quantify the influence of polar clouds on the radiative impact of polar sea ice variability. Our results show that the cloud short-wave cooling effect strongly influences the impact of sea ice variability on the surface radiation budget and does so in a counter-intuitive manner over the polar seas: years with less sea ice and a larger net surface radiative flux show a more negative cloud radiative effect. Our results indicate that 66±2% of this change in the net cloud radiative effect is due to the reduction in surface albedo and that the remaining 34±1 % is due to an increase in cloud cover and optical thickness. The overall cloud radiative damping effect is 56±2 % over the Antarctic and 47±3 % over the Arctic. Thus, present-day cloud properties significantly reduce the net radiative impact of sea ice loss on the Arctic and Antarctic surface radiation budgets. As a result, climate models must accurately represent present-day polar cloud properties in order to capture the surface radiation budget impact of polar sea ice loss and thus the surface albedo feedback.

Description: The comprehensive analysis of carbon stocks and fluxes of managed European forests is a prerequisite to quantify their role in biomass production and climate change mitigation. We applied the Carbon Budget Model (CBM) to 26 European countries, parameterized with country information on the historical forest age structure, management practices, harvest regimes and the main natural disturbances. We modeled the C stocks for the five forest pools plus harvested wood products (HWPs) and the fluxes among these pools from 2000 to 2030. The aim is to quantify, using a consistent modeling framework for all 26 countries, the main C fluxes as affected by land-use changes, natural disturbances and forest management and to assess the impact of specific harvest and afforestation scenarios after 2012 on the mitigation potential of the EU forest sector. Substitution effects and the possible impacts of climate are not included in this analysis. Results show that for the historical period from 2000 to 2012 the net primary productivity (NPP) of the forest pools at the EU level is on average equal to 639 Tg C yr−1. The losses are dominated by heterotrophic respiration (409 Tg C yr−1) and removals (110 Tg C yr−1), with direct fire emissions being only 1 Tg C yr−1, leading to a net carbon stock change (i.e., sink) of 110 Tg C yr−1. Fellings also transferred 28 Tg C yr−1 of harvest residues from biomass to dead organic matter pools. The average annual net sector exchange (NSE) of the forest system, i.e., the carbon stock changes in the forest pools including HWP, equals a sink of 122 Tg C yr−1 (i.e., about 19 % of the NPP) for the historical period, and in 2030 it reaches 126, 101 and 151 Tg C yr−1, assuming constant, increasing (+20 %) and decreasing (−20 %) scenarios, respectively, of both harvest and afforestation rates compared to the historical period. Under the constant harvest rate scenario, our findings show an incipient aging process of the forests existing in 1990: although NPP increases (+7 %), heterotrophic respiration increases at a greater rate (+13 %) and this leads to a decrease in the sink in the forest pools (−6 %) in 2030 compared to the historical period. By comparing the evolution of the biomass as a function of the NPP (i.e., the turnover time) for each country, we highlighted at least three groups of countries and turnover times. This means that, contrary to the assumptions proposed by other authors, this relationship cannot be assumed as a constant for all the EU countries, but specific conditions, such as the harvest rate, the current age structure and the forest composition, may contribute to the country-specific evolution of biomass stocks. The detailed picture of the C fluxes condensed in this study, and their evolution under different harvest scenarios, may represent both a benchmark for similar studies and a basis for broader analyses (e.g., including substitution effects of wood) on the mitigation potential of the EU forest sector.

Description: The terrestrial carbon fluxes show the largest variability among the components of the global carbon cycle and drive most of the temporal variations in the growth rate of atmospheric CO2 (Le Quéré 2014). Understanding the environmental controls and trends of the terrestrial carbon budget is therefore essential to predict the future trajectories of the CO2 airborne fraction and atmospheric concentrations. In the present work, patterns and controls of the inter-annual variability (IAV) of carbon Net Ecosystem Exchange (NEE) have been analysed using three different data-streams: ecosystem level observations from the FLUXNET database (La Thuille and 2015 releases), the MPI-MTE bottom-up product resulting from the global up-scaling of site-level fluxes, and the Jena CarboScope Inversion, a top-down estimate of surface fluxes obtained from observed CO2 concentrations and an atmospheric transport model. Consistencies and discrepancies in the temporal and spatial patterns and in the climatic and physiological controls of IAV were investigated between the three data sources. The global average of IAV at FLUXNET sites (~ 120 gC m−2 y−1), quantified as the standard deviation of annual NEE, was observed to peak in arid ecosystems and to be almost six times larger than the values calculated from the two global products (15 and 20 gC m−2 y−1 for MPI-MTE and Jena inversion, respectively). The two data-driven global products show that most of the temporal variability observed in the last three decades is due to yearly anomalies, whereas the temporal trends explain only about 15 % of the variability in the MPI-MTE product and 20 % in the Jena Inversion product. Both at site level and at global scale, the IAV of NEE is driven by the gross primary productivity and in particular by the cumulative carbon flux during the months when land acts as a sink. Altogether these results offer a broad view on the magnitude, spatial patterns and environmental drivers of IAV from a variety of data sources, that can be instrumental to improve our understanding of the terrestrial carbon budget and to validate the predictions of land surface models.

Description: Extreme climatic events, such as droughts and heat stress induce anomalies in ecosystem-atmosphere CO2 fluxes, such as gross primary production (GPP) and ecosystem respiration (Reco), and, hence, can change the net ecosystem carbon balance. However, despite our increasing understanding of the underlying mechanisms, the magnitudes of the impacts of different types of extremes on GPP and Reco within and between ecosystems remain poorly predicted. Here we aim to identify the major factors controlling the amplitude of extreme event impacts on GPP, Reco, and the resulting net ecosystem production (NEP). We focus on the impacts of heat and drought and their combination. We identified hydrometeorological extreme events in consistently downscaled water availability and temperature measurements over a 30 year time period. We then used FLUXNET eddy-covariance flux measurements to estimate the CO2 flux anomalies during these extreme events across dominant vegetation types and climate zones. Overall, our results indicate that short-term heat extremes increased respiration more strongly than they down-regulated GPP, resulting in a moderate reduction of the ecosystem’s carbon sink potential. In the absence of heat stress, droughts tended to have smaller and similarly dampening effects on both GPP and Reco, and, hence, often resulted in neutral NEP responses. The combination of drought and heat typically led to a strong decrease in GPP, whereas heat and drought impacts on respiration partially offset each other. Taken together, compound heat and drought events led to the strongest C sink reduction compared to any single-factor extreme. A key insight of this paper, however, is that duration matters most: for heat stress during droughts, the magnitude of impacts systematically increased with duration, whereas under heat stress without drought, the response of Reco over time turned from an initial increase to a down-regulation after about two weeks. This confirms earlier theories that not only the magnitude but also the duration of an extreme event determines its impact. Our study corroborates the results of several local site-level case studies, but as a novelty generalizes these findings at the global scale. Specifically, we find that the different response functions of the two antipodal land-atmosphere fluxes GPP and Reco can also result in increasing NEP during certain extreme conditions. Apparently counterintuitive findings of this kind bear great potential for scrutinizing the mechanisms implemented in state-of-the-art terrestrial biosphere models and provide a benchmark for future model development and testing.

Description: The terrestrial carbon fluxes show the largest variability among the components of the global carbon cycle and drive most of the temporal variations in the growth rate of atmospheric CO2. Understanding the environmental controls and trends of the terrestrial carbon budget is therefore essential to predict the future trajectories of the CO2 airborne fraction and atmospheric concentrations. In the present work, patterns and controls of the inter-annual variability (IAV) of carbon net ecosystem exchange (NEE) have been analysed using three different data streams: ecosystem-level observations from the FLUXNET database (La Thuile and 2015 releases), the MPI-MTE (model tree ensemble) bottom–up product resulting from the global upscaling of site-level fluxes, and the Jena CarboScope Inversion, a top–down estimate of surface fluxes obtained from observed CO2 concentrations and an atmospheric transport model. Consistencies and discrepancies in the temporal and spatial patterns and in the climatic and physiological controls of IAV were investigated between the three data sources. Results show that the global average of IAV at FLUXNET sites, quantified as the standard deviation of annual NEE, peaks in arid ecosystems and amounts to  ∼  120 gC m−2 y−1, almost 6 times more than the values calculated from the two global products (15 and 20 gC m−2 y−1 for MPI-MTE and the Jena Inversion, respectively). Most of the temporal variability observed in the last three decades of the MPI-MTE and Jena Inversion products is due to yearly anomalies, whereas the temporal trends explain only about 15 and 20 % of the variability, respectively. Both at the site level and on a global scale, the IAV of NEE is driven by the gross primary productivity and in particular by the cumulative carbon flux during the months when land acts as a sink. Altogether these results offer a broad view on the magnitude, spatial patterns and environmental drivers of IAV from a variety of data sources that can be instrumental to improve our understanding of the terrestrial carbon budget and to validate the predictions of land surface models.

Description: Terrestrial ecosystem respiration (Re) is the major source of CO2 release and constitutes the second largest carbon flux between the biosphere and atmosphere. Therefore, climate-driven changes of Re may greatly impact on future atmospheric CO2 concentration. The aim of this study was to derive an air temperature threshold for identifying the driving climate forces of the respiratory process in terrestrial ecosystems within different temperature zones. For this purpose, a global dataset of 647 site-years of ecosystem flux data collected at 152 sites has been examined. Our analysis revealed an ecosystem threshold of mean annual air temperature (MAT) of 11 ± 2.3 °C. In ecosystems with the MAT below this threshold, the maximum Re rates were primarily dependent on temperature and respiration was mainly a temperature-driven process. On the contrary, in ecosystems with the MAT greater than 11 ± 2.3 °C, in addition to temperature, other driving forces, such as water availability and surface heat flux, became significant drivers of the maximum Re rates and respiration was a multi-factor-driven process. The information derived from this study highlight the key role of temperature as main controlling factor of the maximum Re rates on a large fraction of the terrestrial biosphere, while other driving forces reduce the maximum Re rates and temperature sensitivity of the respiratory process. These findings are particularly relevant under the current scenario of rapid global warming, given that the potential climate-induced changes in ecosystem respiration may lead to substantial anomalies in the seasonality and magnitude of the terrestrial carbon budget.

Description: Accurate multi-point monitoring systems are required to derive atmospheric measurements of greenhouse gas concentrations both for the calculation of surface fluxes with inversion transport models and for the estimation of non-turbulent components of the mass balance equation (i.e. advection and storage fluxes) at eddy covariance sites. When a single analyser is used to monitor multiple sampling points, the deployment of buffer volumes (BV) along sampling lines can reduce the uncertainty due to the discrete temporal sampling of the signal. In order to optimize the use of buffer volumes we explored various setups by simulating their effect on time series of high-frequency CO2 concentration collected at three Fluxnet sites. Besides, we proposed a novel scheme to calculate half hourly weighted averages from discrete point samples, accounting for the probabilistic fraction of the signal generated in the averaging period. Results show that the use of BV with the new averaging scheme reduces the mean absolute errors (MAE) up to 80 % compared to set-up without BV and up to 60 % compared to the case with BV and a standard, non-weighted averaging scheme. The MAE of CO2 concentration measurements was observed to depend on the variability of the concentration field and on the size of BV, which therefore have to be carefully dimensioned. The optimal volume size depends on two main features of the instrumental setup: the number of measurement points and the time needed to sample at one point (i.e. line purging plus sampling time). A linear and consistent relationship was observed at all sites between the sampling frequency, which summarizes the two features mentioned above, and the renewal frequency associated to the volume. Ultimately, this empirical relationship can be applied to estimate the optimal volume size according to the technical specifications of the sampling system.