Monitoring of carbon-water fluxes at Eurasian meteorological stations using random forest and remote sensing.

Simulating the carbon-water fluxes at more widely distributed meteorological stations based on the sparsely and unevenly distributed eddy covariance flux stations is needed to accurately understand the carbon-water cycle of terrestrial ecosystems. We established a new framework consisting of machine learning, determination coefficient (R2), Euclidean distance, and remote sensing (RS), to simulate the daily net ecosystem carbon dioxide exchange (NEE) and water flux (WF) of the Eurasian meteorological stations using a random forest model or/and RS. The daily NEE and WF datasets with RS-based information (NEE-RS and WF-RS) for 3774 and 4427 meteorological stations during 2002-2020 were produced, respectively. And the daily NEE and WF datasets without RS-based information (NEE-WRS and WF-WRS) for 4667 and 6763 meteorological stations during 1983-2018 were generated, respectively. For each meteorological station, the carbon-water fluxes meet accuracy requirements and have quasi-observational properties. These four carbon-water flux datasets have great potential to improve the assessments of the ecosystem carbon-water dynamics.

Leaf-level coordination principles propagate to the ecosystem scale.

Fundamental axes of variation in plant traits result from trade-offs between costs and benefits of resource-use strategies at the leaf scale. However, it is unclear whether similar trade-offs propagate to the ecosystem level. Here, we test whether trait correlation patterns predicted by three well-known leaf- and plant-level coordination theories - the leaf economics spectrum, the global spectrum of plant form and function, and the least-cost hypothesis - are also observed between community mean traits and ecosystem processes. We combined ecosystem functional properties from FLUXNET sites, vegetation properties, and community mean plant traits into three corresponding principal component analyses. We find that the leaf economics spectrum (90 sites), the global spectrum of plant form and function (89 sites), and the least-cost hypothesis (82 sites) all propagate at the ecosystem level. However, we also find evidence of additional scale-emergent properties. Evaluating the coordination of ecosystem functional properties may aid the development of more realistic global dynamic vegetation models with critical empirical data, reducing the uncertainty of climate change projections.

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.

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.

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.

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'.

Nutrients and water availability constrain the seasonality of vegetation activity in a Mediterranean ecosystem.

Anthropogenic nitrogen (N) deposition and resulting differences in ecosystem N and phosphorus (P) ratios are expected to impact photosynthetic capacity, that is, maximum gross primary productivity (GPPmax ). However, the interplay between N and P availability with other critical resources on seasonal dynamics of ecosystem productivity remains largely unknown. In a Mediterranean tree-grass ecosystem, we established three landscape-level (24 ha) nutrient addition treatments: N addition (NT), N and P addition (NPT), and a control site (CT). We analyzed the response of ecosystem to altered nutrient stoichiometry using eddy covariance fluxes measurements, satellite observations, and digital repeat photography. A set of metrics, including phenological transition dates (PTDs; timing of green-up and dry-down), slopes during green-up and dry-down period, and seasonal amplitude, were extracted from time series of GPPmax and used to represent the seasonality of vegetation activity. The seasonal amplitude of GPPmax was higher for NT and NPT than CT, which was attributed to changes in structure and physiology induced by fertilization. PTDs were mainly driven by rainfall and exhibited no significant differences among treatments during the green-up period. Yet, both fertilized sites senesced earlier during the dry-down period (17-19 days), which was more pronounced in the NT due to larger evapotranspiration and water usage. Fertilization also resulted in a faster increase in GPPmax during the green-up period and a sharper decline in GPPmax during the dry-down period, with less prominent decline response in NPT. Overall, we demonstrated seasonality of vegetation activity was altered after fertilization and the importance of nutrient-water interaction in such water-limited ecosystems. With the projected warming-drying trend, the positive effects of N fertilization induced by N deposition on GPPmax may be counteracted by an earlier and faster dry-down in particular in areas where the N:P ratio increases, with potential impact on the carbon cycle of water-limited ecosystems.

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.

Potential and limitations of inferring ecosystem photosynthetic capacity from leaf functional traits.

The aim of this study was to systematically analyze the potential and limitations of using plant functional trait observations from global databases versus in situ data to improve our understanding of vegetation impacts on ecosystem functional properties (EFPs). Using ecosystem photosynthetic capacity as an example, we first provide an objective approach to derive robust EFP estimates from gross primary productivity (GPP) obtained from eddy covariance flux measurements. Second, we investigate the impact of synchronizing EFPs and plant functional traits in time and space to evaluate their relationships, and the extent to which we can benefit from global plant trait databases to explain the variability of ecosystem photosynthetic capacity. Finally, we identify a set of plant functional traits controlling ecosystem photosynthetic capacity at selected sites. Suitable estimates of the ecosystem photosynthetic capacity can be derived from light response curve of GPP responding to radiation (photosynthetically active radiation or absorbed photosynthetically active radiation). Although the effect of climate is minimized in these calculations, the estimates indicate substantial interannual variation of the photosynthetic capacity, even after removing site-years with confounding factors like disturbance such as fire events. The relationships between foliar nitrogen concentration and ecosystem photosynthetic capacity are tighter when both of the measurements are synchronized in space and time. When using multiple plant traits simultaneously as predictors for ecosystem photosynthetic capacity variation, the combination of leaf carbon to nitrogen ratio with leaf phosphorus content explains the variance of ecosystem photosynthetic capacity best (adjusted R2 = 0.55). Overall, this study provides an objective approach to identify links between leaf level traits and canopy level processes and highlights the relevance of the dynamic nature of ecosystems. Synchronizing measurements of eddy covariance fluxes and plant traits in time and space is shown to be highly relevant to better understand the importance of intra- and interspecific trait variation on ecosystem functioning.

Ground-based optical measurements at European flux sites: a review of methods, instruments and current controversies.

This paper reviews the currently available optical sensors, their limitations and opportunities for deployment at Eddy Covariance (EC) sites in Europe. This review is based on the results obtained from an online survey designed and disseminated by the Co-cooperation in Science and Technology (COST) Action ESO903-"Spectral Sampling Tools for Vegetation Biophysical Parameters and Flux Measurements in Europe" that provided a complete view on spectral sampling activities carried out within the different research teams in European countries. The results have highlighted that a wide variety of optical sensors are in use at flux sites across Europe, and responses further demonstrated that users were not always fully aware of the key issues underpinning repeatability and the reproducibility of their spectral measurements. The key findings of this survey point towards the need for greater awareness of the need for standardisation and development of a common protocol of optical sampling at the European EC sites.