Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown.

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

Environmental factors contributing to variations in CO2 flux over a barley-rice double-cropping paddy field in the Korean Peninsula.

Understanding the CO2 flux over agricultural crop fields is critical because the temporal cycle is driven by both ecological environment and anthropogenic change. We analyzed the net ecosystem exchange of CO2 measured over a barley-rice double-cropping field using the eddy covariance method for 5 years. We conducted gap-filling based on u*-threshold criteria and partitioned the net ecosystem exchange into gross primary production and respiration. The relative importance analysis of solar radiation, temperature, soil heat flux, soil water content, and vapor deficit revealed that solar radiation and temperature were the dominant contributors to net ecosystem exchange. The annual variation in the net ecosystem exchange followed a bimodal pattern driven by CO2 uptake by both barley and rice, displaying two negative peaks in late April and mid-August. The elongation stages of the crops exhibited the highest flux. Gross primary production and respiration were closely related to solar radiation and nighttime temperature, respectively. The relative importance of the other environmental variables was affected by the cultivation season and irrigation water. In the period of rice cultivation, respiration was approximately 3 µmol m-2 s-1 higher during rice drainage than during the flooded period. The accumulated net ecosystem production was estimated to be 315 gC m-2 and 349 gC m-2 for the barley and rice growing periods, respectively, and 649 gC m-2 for the annual total. These values are comparable with the results of other studies on barley-rice double-cropping fields.

A remote sensing-based three-source energy balance model to improve global estimations of evapotranspiration in semi-arid tree-grass ecosystems.

It is well documented that energy balance and other remote sensing-based evapotranspiration (ET) models face greater uncertainty over water-limited tree-grass ecosystems (TGEs), representing nearly 1/6th of the global land surface. Their dual vegetation strata, the grass-dominated understory and tree-dominated overstory, make for distinct structural, physiological and phenological characteristics, which challenge models compared to more homogeneous and energy-limited ecosystems. Along with this, the contribution of grasses and trees to total transpiration (T), along with their different climatic drivers, is still largely unknown nor quantified in TGEs. This study proposes a thermal-based three-source energy balance (3SEB) model, accommodating an additional vegetation source within the well-known two-source energy balance (TSEB) model. The model was implemented at both tower and continental scales using eddy-covariance (EC) TGE sites, with variable tree canopy cover and rainfall (P) regimes and Meteosat Second Generation (MSG) images. 3SEB robustly simulated latent heat (LE) and related energy fluxes in all sites (Tower: LE RMSD ~60 W/m2 ; MSG: LE RMSD ~90 W/m2 ), improving over both TSEB and seasonally changing TSEB (TSEB-2S) models. In addition, 3SEB inherently partitions water fluxes between the tree, grass and soil sources. The modelled T correlated well with EC T estimates (r > .76), derived from a machine learning ET partitioning method. The T/ET was found positively related to both P and leaf area index, especially compared to the decomposed grass understory T/ET. However, trees and grasses had contrasting relations with respect to monthly P. These results demonstrate the importance in decomposing total ET into the different vegetation sources, as they have distinct climatic drivers, and hence, different relations to seasonal water availability. These promising results improved ET and energy flux estimations over complex TGEs, which may contribute to enhance global drought monitoring and understanding, and their responses to climate change feedbacks.

Heatwave breaks down the linearity between sun-induced fluorescence and gross primary production.

Sun-induced fluorescence in the far-red region (SIF) is increasingly used as a remote and proximal-sensing tool capable of tracking vegetation gross primary production (GPP). However, the use of SIF to probe changes in GPP is challenged during extreme climatic events, such as heatwaves. Here, we examined how the 2018 European heatwave (HW) affected the GPP-SIF relationship in evergreen broadleaved trees with a relatively invariant canopy structure. To do so, we combined canopy-scale SIF measurements, GPP estimated from an eddy covariance tower, and active pulse amplitude modulation fluorescence. The HW caused an inversion of the photosynthesis-fluorescence relationship at both the canopy and leaf scales. The highly nonlinear relationship was strongly shaped by nonphotochemical quenching (NPQ), that is, a dissipation mechanism to protect from the adverse effects of high light intensity. During the extreme heat stress, plants experienced a saturation of NPQ, causing a change in the allocation of energy dissipation pathways towards SIF. Our results show the complex modulation of the NPQ-SIF-GPP relationship at an extreme level of heat stress, which is not completely represented in state-of-the-art coupled radiative transfer and photosynthesis models.

Uncovering the critical soil moisture thresholds of plant water stress for European ecosystems.

Understanding the critical soil moisture (SM) threshold (θcrit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θcrit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θcrit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θcrit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θcrit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θcrit . We further found that soil texture dominates the spatial variability of θcrit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θcrit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.

COSORE: A community database for continuous soil respiration and other soil-atmosphere greenhouse gas flux data.

Globally, soils store two to three times as much carbon as currently resides in the atmosphere, and it is critical to understand how soil greenhouse gas (GHG) emissions and uptake will respond to ongoing climate change. In particular, the soil-to-atmosphere CO2 flux, commonly though imprecisely termed soil respiration (RS ), is one of the largest carbon fluxes in the Earth system. An increasing number of high-frequency RS measurements (typically, from an automated system with hourly sampling) have been made over the last two decades; an increasing number of methane measurements are being made with such systems as well. Such high frequency data are an invaluable resource for understanding GHG fluxes, but lack a central database or repository. Here we describe the lightweight, open-source COSORE (COntinuous SOil REspiration) database and software, that focuses on automated, continuous and long-term GHG flux datasets, and is intended to serve as a community resource for earth sciences, climate change syntheses and model evaluation. Contributed datasets are mapped to a single, consistent standard, with metadata on contributors, geographic location, measurement conditions and ancillary data. The design emphasizes the importance of reproducibility, scientific transparency and open access to data. While being oriented towards continuously measured RS , the database design accommodates other soil-atmosphere measurements (e.g. ecosystem respiration, chamber-measured net ecosystem exchange, methane fluxes) as well as experimental treatments (heterotrophic only, etc.). We give brief examples of the types of analyses possible using this new community resource and describe its accompanying R software package.

Drought and heatwave impacts on semi-arid ecosystems' carbon fluxes along a precipitation gradient.

The inter-annual variability (IAV) of the terrestrial carbon cycle is tightly linked to the variability of semi-arid ecosystems. Thus, it is of utmost importance to understand what the main meteorological drivers for the IAV of such ecosystems are, and how they respond to extreme events such as droughts and heatwaves. To shed light onto these questions, we analyse the IAV of carbon fluxes, its relation with meteorological variables, and the impact of compound drought and heatwave on the carbon cycle of two similar ecosystems, along a precipitation gradient. A four-year long dataset from 2016 to 2019 was used for the FLUXNET sites ES-LMa and ES-Abr, located in central (39°56'25″ N 5°46'28″ W) and southeastern (38°42'6″ N 6°47'9″ W) Spain. We analyse the physiological impact of compound drought and heatwave on the dominant tree species, Quercus ilex. Our results show that the gross primary productivity of the wetter ecosystem was less sensitive to changes in soil water content, compared to the dryer site. Still, the wetter ecosystem was a source of CO2 each year, owing to large ecosystem respiration during summer; while the dry site turned into a CO2 sink during wet years. Overall, the impact of the summertime compound event on annual CO2 fluxes was marginal at both sites, compared to drought events during spring or autumn. This highlights that drought timing is crucial to determine the annual carbon fluxes in these semi-arid ecosystems. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Sensitivity of gross primary productivity to climatic drivers during the summer drought of 2018 in Europe.

In summer 2018, Europe experienced a record drought, but it remains unknown how the drought affected ecosystem carbon dynamics. Using observations from 34 eddy covariance sites in different biomes across Europe, we studied the sensitivity of gross primary productivity (GPP) to environmental drivers during the summer drought of 2018 versus the reference summer of 2016. We found a greater drought-induced decline of summer GPP in grasslands (-38%) than in forests (-10%), which coincided with reduced evapotranspiration and soil water content (SWC). As compared to the 'normal year' of 2016, GPP in different ecosystems exhibited more negative sensitivity to summer air temperature (Ta) but stronger positive sensitivity to SWC during summer drought in 2018, that is, a stronger reduction of GPP with soil moisture deficit. We found larger negative effects of Ta and vapour pressure deficit (VPD) but a lower positive effect of photosynthetic photon flux density on GPP in 2018 compared to 2016, which contributed to reduced summer GPP in 2018. Our results demonstrate that high temperature-induced increases in VPD and decreases in SWC aggravated drought impacts on GPP. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Altered energy partitioning across terrestrial ecosystems in the European drought year 2018.

Drought and heat events, such as the 2018 European drought, interact with the exchange of energy between the land surface and the atmosphere, potentially affecting albedo, sensible and latent heat fluxes, as well as CO2 exchange. Each of these quantities may aggravate or mitigate the drought, heat, their side effects on productivity, water scarcity and global warming. We used measurements of 56 eddy covariance sites across Europe to examine the response of fluxes to extreme drought prevailing most of the year 2018 and how the response differed across various ecosystem types (forests, grasslands, croplands and peatlands). Each component of the surface radiation and energy balance observed in 2018 was compared to available data per site during a reference period 2004-2017. Based on anomalies in precipitation and reference evapotranspiration, we classified 46 sites as drought affected. These received on average 9% more solar radiation and released 32% more sensible heat to the atmosphere compared to the mean of the reference period. In general, drought decreased net CO2 uptake by 17.8%, but did not significantly change net evapotranspiration. The response of these fluxes differed characteristically between ecosystems; in particular, the general increase in the evaporative index was strongest in peatlands and weakest in croplands. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Bigleaf-An R package for the calculation of physical and physiological ecosystem properties from eddy covariance data.

We present the R package bigleaf (version 0.6.5), an open source toolset for the derivation of meteorological, aerodynamic, and physiological ecosystem properties from eddy covariance (EC) flux observations and concurrent meteorological measurements. A 'big-leaf' framework, in which vegetation is represented as a single, uniform layer, is employed to infer bulk ecosystem characteristics top-down from the measured fluxes. Central to the package is the calculation of a bulk surface/canopy conductance (Gs/Gc) and a bulk aerodynamic conductance (Ga), with the latter including formulations for the turbulent and canopy boundary layer components. The derivation of physical land surface characteristics such as surface roughness parameters, wind profile, aerodynamic and radiometric surface temperature, surface vapor pressure deficit (VPD), potential evapotranspiration (ET), imposed and equilibrium ET, as well as vegetation-atmosphere decoupling coefficients, is described. The package further provides calculation routines for physiological ecosytem properties (stomatal slope parameters, stomatal sensitivity to VPD, bulk intercellular CO2 concentration, canopy photosynthetic capacity), energy balance characteristics (closure, biochemical energy), ancillary meteorological variables (psychrometric constant, saturation vapor pressure, air density, etc.), customary unit interconversions and data filtering. The target variables can be calculated with a different degree of complexity, depending on the amount of available site-specific information. The utilities of the package are demonstrated for three single-level (above-canopy) eddy covariance sites representing a temperate grassland, a temperate needle-leaf forest, and a Mediterranean evergreen broadleaf forest. The routines are further tested for a two-level EC site (tree and grass layer) located in a Mediterranean oak savanna. The limitations and the ecophysiological interpretation of the derived ecosystem properties are discussed and practical guidelines are given. The package provides the basis for a consistent, physically sound, and reproducible characterization of biometeorological conditions and ecosystem physiology, and is applicable to EC sites across vegetation types and climatic conditions with minimal ancillary data requirements.

Plant functional traits and canopy structure control the relationship between photosynthetic CO2 uptake and far-red sun-induced fluorescence in a Mediterranean grassland under different nutrient availability.

Sun-induced fluorescence (SIF) in the far-red region provides a new noninvasive measurement approach that has the potential to quantify dynamic changes in light-use efficiency and gross primary production (GPP). However, the mechanistic link between GPP and SIF is not completely understood. We analyzed the structural and functional factors controlling the emission of SIF at 760 nm (F760 ) in a Mediterranean grassland manipulated with nutrient addition of nitrogen (N), phosphorous (P) or nitrogen-phosphorous (NP). Using the soil-canopy observation of photosynthesis and energy (SCOPE) model, we investigated how nutrient-induced changes in canopy structure (i.e. changes in plant forms abundance that influence leaf inclination distribution function, LIDF) and functional traits (e.g. N content in dry mass of leaves, N%, Chlorophyll a+b concentration (Cab) and maximum carboxylation capacity (Vcmax )) affected the observed linear relationship between F760 and GPP. We conclude that the addition of nutrients imposed a change in the abundance of different plant forms and biochemistry of the canopy that controls F760 . Changes in canopy structure mainly control the GPP-F760 relationship, with a secondary effect of Cab and Vcmax . In order to exploit F760 data to model GPP at the global/regional scale, canopy structural variability, biodiversity and functional traits are important factors that have to be considered.