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.

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.

Joint optimization of land carbon uptake and albedo can help achieve moderate instantaneous and long-term cooling effects.

Both carbon dioxide uptake and albedo of the land surface affect global climate. However, climate change mitigation by increasing carbon uptake can cause a warming trade-off by decreasing albedo, with most research focusing on afforestation and its interaction with snow. Here, we present carbon uptake and albedo observations from 176 globally distributed flux stations. We demonstrate a gradual decline in maximum achievable annual albedo as carbon uptake increases, even within subgroups of non-forest and snow-free ecosystems. Based on a paired-site permutation approach, we quantify the likely impact of land use on carbon uptake and albedo. Shifting to the maximum attainable carbon uptake at each site would likely cause moderate net global warming for the first approximately 20 years, followed by a strong cooling effect. A balanced policy co-optimizing carbon uptake and albedo is possible that avoids warming on any timescale, but results in a weaker long-term cooling effect.

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.

Warming homogenizes apparent temperature sensitivity of ecosystem respiration.

Warming-induced carbon loss through terrestrial ecosystem respiration (Re) is likely getting stronger in high latitudes and cold regions because of the more rapid warming and higher temperature sensitivity of Re (Q 10). However, it is not known whether the spatial relationship between Q 10 and temperature also holds temporally under a future warmer climate. Here, we analyzed apparent Q 10 values derived from multiyear observations at 74 FLUXNET sites spanning diverse climates and biomes. We found warming-induced decline in Q 10 is stronger at colder regions than other locations, which is consistent with a meta-analysis of 54 field warming experiments across the globe. We predict future warming will shrink the global variability of Q 10 values to an average of 1.44 across the globe under a high emission trajectory (RCP 8.5) by the end of the century. Therefore, warming-induced carbon loss may be less than previously assumed because of Q 10 homogenization in a warming world.

Temperature thresholds of ecosystem respiration at a global scale.

Ecosystem respiration is a major component of the global terrestrial carbon cycle and is strongly influenced by temperature. The global extent of the temperature-ecosystem respiration relationship, however, has not been fully explored. Here, we test linear and threshold models of ecosystem respiration across 210 globally distributed eddy covariance sites over an extensive temperature range. We find thresholds to the global temperature-ecosystem respiration relationship at high and low air temperatures and mid soil temperatures, which represent transitions in the temperature dependence and sensitivity of ecosystem respiration. Annual ecosystem respiration rates show a markedly reduced temperature dependence and sensitivity compared to half-hourly rates, and a single mid-temperature threshold for both air and soil temperature. Our study indicates a distinction in the influence of environmental factors, including temperature, on ecosystem respiration between latitudinal and climate gradients at short (half-hourly) and long (annual) timescales. Such climatological differences in the temperature sensitivity of ecosystem respiration have important consequences for the terrestrial net carbon sink under ongoing climate change.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites.

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.

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

SoilTemp: A global database of near-surface temperature.

Current analyses and predictions of spatially explicit patterns and processes in ecology most often rely on climate data interpolated from standardized weather stations. This interpolated climate data represents long-term average thermal conditions at coarse spatial resolutions only. Hence, many climate-forcing factors that operate at fine spatiotemporal resolutions are overlooked. This is particularly important in relation to effects of observation height (e.g. vegetation, snow and soil characteristics) and in habitats varying in their exposure to radiation, moisture and wind (e.g. topography, radiative forcing or cold-air pooling). Since organisms living close to the ground relate more strongly to these microclimatic conditions than to free-air temperatures, microclimatic ground and near-surface data are needed to provide realistic forecasts of the fate of such organisms under anthropogenic climate change, as well as of the functioning of the ecosystems they live in. To fill this critical gap, we highlight a call for temperature time series submissions to SoilTemp, a geospatial database initiative compiling soil and near-surface temperature data from all over the world. Currently, this database contains time series from 7,538 temperature sensors from 51 countries across all key biomes. The database will pave the way toward an improved global understanding of microclimate and bridge the gap between the available climate data and the climate at fine spatiotemporal resolutions relevant to most organisms and ecosystem processes.

Publisher Correction: Stand age and species richness dampen interannual variation of ecosystem-level photosynthetic capacity.

In the version of this Article originally published, the wrong Supplementary Information pdf was uploaded, in which the figures did not correspond with those mentioned in the main text and the R code was not presented properly. This has now been replaced.

Solar-induced chlorophyll fluorescence exhibits a universal relationship with gross primary productivity across a wide variety of biomes.

In our recent study in Global Change Biology (Li et al., ), we examined the relationship between solar-induced chlorophyll fluorescence (SIF) measured from the Orbiting Carbon Observatory-2 (OCO-2) and gross primary productivity (GPP) derived from eddy covariance flux towers across the globe, and we discovered that there is a nearly universal relationship between SIF and GPP across a wide variety of biomes. This finding reveals the tremendous potential of SIF for accurately mapping terrestrial photosynthesis globally.

Solar-induced chlorophyll fluorescence is strongly correlated with terrestrial photosynthesis for a wide variety of biomes: First global analysis based on OCO-2 and flux tower observations.

Solar-induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite-observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory-2 (OCO-2) provides the first opportunity to examine the SIF-GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO-2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO-2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm (R2  = 0.72, p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R2  = 0.57-0.79, p < 0.0001) except evergreen broadleaf forests (R2  = 0.16, p < 0.05) at the daily timescale. A higher slope was found for C4 grasslands and croplands than for C3 ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome-specific SIF-GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO-2 SIF generally had a better performance for predicting GPP than satellite-derived vegetation indices and a light use efficiency model. The universal SIF-GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO-2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.

Stand age and species richness dampen interannual variation of ecosystem-level photosynthetic capacity.

The total uptake of carbon dioxide by ecosystems via photosynthesis (gross primary productivity, GPP) is the largest flux in the global carbon cycle. A key ecosystem functional property determining GPP is the photosynthetic capacity at light saturation (GPPsat), and its interannual variability (IAV) is propagated to the net land-atmosphere exchange of CO2. Given the importance of understanding the IAV in CO2 fluxes for improving the predictability of the global carbon cycle, we have tested a range of alternative hypotheses to identify potential drivers of the magnitude of IAV in GPPsat in forest ecosystems. Our results show that while the IAV in GPPsat within sites is closely related to air temperature and soil water availability fluctuations, the magnitude of IAV in GPPsat is related to stand age and biodiversity (R2 = 0.55, P < 0.0001). We find that the IAV of GPPsat is greatly reduced in older and more diverse forests, and is higher in younger forests with few dominant species. Older and more diverse forests seem to dampen the effect of climate variability on the carbon cycle irrespective of forest type. Preserving old forests and their diversity would therefore be beneficial in reducing the effect of climate variability on Earth's forest ecosystems.

Stand age and species richness dampen interannual variation of ecosystem-level photosynthetic capacity.

The total uptake of carbon dioxide by ecosystems via photosynthesis (gross primary productivity, GPP) is the largest flux in the global carbon cycle. A key ecosystem functional property determining GPP is the photosynthetic capacity at light saturation (GPPsat), and its interannual variability (IAV) is propagated to the net land-atmosphere exchange of CO2. Given the importance of understanding the IAV in CO2 fluxes for improving the predictability of the global carbon cycle, we have tested a range of alternative hypotheses to identify potential drivers of the magnitude of IAV in GPPsat in forest ecosystems. Our results show that while the IAV in GPPsat within sites is closely related to air temperature and soil water availability fluctuations, the magnitude of IAV in GPPsat is related to stand age and biodiversity (R2 = 0.55, P < 0.0001). We find that the IAV of GPPsat is greatly reduced in older and more diverse forests, and is higher in younger forests with few dominant species. Older and more diverse forests seem to dampen the effect of climate variability on the carbon cycle irrespective of forest type. Preserving old forests and their diversity would therefore be beneficial in reducing the effect of climate variability on Earth's forest ecosystems.

Do the energy fluxes and surface conductance of boreal coniferous forests in Europe scale with leaf area?

Earth observing systems are now routinely used to infer leaf area index (LAI) given its significance in spatial aggregation of land surface fluxes. Whether LAI is an appropriate scaling parameter for daytime growing season energy budget, surface conductance (Gs ), water- and light-use efficiency and surface-atmosphere coupling of European boreal coniferous forests was explored using eddy-covariance (EC) energy and CO2 fluxes. The observed scaling relations were then explained using a biophysical multilayer soil-vegetation-atmosphere transfer model as well as by a bulk Gs representation. The LAI variations significantly alter radiation regime, within-canopy microclimate, sink/source distributions of CO2 , H2 O and heat, and forest floor fluxes. The contribution of forest floor to ecosystem-scale energy exchange is shown to decrease asymptotically with increased LAI, as expected. Compared with other energy budget components, dry-canopy evapotranspiration (ET) was reasonably 'conservative' over the studied LAI range 0.5-7.0 m2 m-2 . Both ET and Gs experienced a minimum in the LAI range 1-2 m2 m-2 caused by opposing nonproportional response of stomatally controlled transpiration and 'free' forest floor evaporation to changes in canopy density. The young forests had strongest coupling with the atmosphere while stomatal control of energy partitioning was strongest in relatively sparse (LAI ~2 m2 m-2 ) pine stands growing on mineral soils. The data analysis and model results suggest that LAI may be an effective scaling parameter for net radiation and its partitioning but only in sparse stands (LAI <3 m2 m-2 ). This finding emphasizes the significance of stand-replacing disturbances on the controls of surface energy exchange. In denser forests, any LAI dependency varies with physiological traits such as light-saturated water-use efficiency. The results suggest that incorporating species traits and site conditions are necessary when LAI is used in upscaling energy exchanges of boreal coniferous forests.

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.

Satellite chlorophyll fluorescence measurements reveal large-scale decoupling of photosynthesis and greenness dynamics in boreal evergreen forests.

Mid-to-high latitude forests play an important role in the terrestrial carbon cycle, but the representation of photosynthesis in boreal forests by current modelling and observational methods is still challenging. In particular, the applicability of existing satellite-based proxies of greenness to indicate photosynthetic activity is hindered by small annual changes in green biomass of the often evergreen tree population and by the confounding effects of background materials such as snow. As an alternative, satellite measurements of sun-induced chlorophyll fluorescence (SIF) can be used as a direct proxy of photosynthetic activity. In this study, the start and end of the photosynthetically active season of the main boreal forests are analysed using spaceborne SIF measurements retrieved from the GOME-2 instrument and compared to that of green biomass, proxied by vegetation indices including the Enhanced Vegetation Index (EVI) derived from MODIS data. We find that photosynthesis and greenness show a similar seasonality in deciduous forests. In high-latitude evergreen needleleaf forests, however, the length of the photosynthetically active period indicated by SIF is up to 6 weeks longer than the green biomass changing period proxied by EVI, with SIF showing a start-of-season of approximately 1 month earlier than EVI. On average, the photosynthetic spring recovery as signalled by SIF occurs as soon as air temperatures exceed the freezing point (2-3 °C) and when the snow on the ground has not yet completely melted. These findings are supported by model data of gross primary production and a number of other studies which evaluated in situ observations of CO2 fluxes, meteorology and the physiological state of the needles. Our results demonstrate the sensitivity of space-based SIF measurements to light-use efficiency of boreal forests and their potential for an unbiased detection of photosynthetic activity even under the challenging conditions interposed by evergreen boreal ecosystems.

The uncertain climate footprint of wetlands under human pressure.

Significant climate risks are associated with a positive carbon-temperature feedback in northern latitude carbon-rich ecosystems, making an accurate analysis of human impacts on the net greenhouse gas balance of wetlands a priority. Here, we provide a coherent assessment of the climate footprint of a network of wetland sites based on simultaneous and quasi-continuous ecosystem observations of CO2 and CH4 fluxes. Experimental areas are located both in natural and in managed wetlands and cover a wide range of climatic regions, ecosystem types, and management practices. Based on direct observations we predict that sustained CH4 emissions in natural ecosystems are in the long term (i.e., several centuries) typically offset by CO2 uptake, although with large spatiotemporal variability. Using a space-for-time analogy across ecological and climatic gradients, we represent the chronosequence from natural to managed conditions to quantify the "cost" of CH4 emissions for the benefit of net carbon sequestration. With a sustained pulse-response radiative forcing model, we found a significant increase in atmospheric forcing due to land management, in particular for wetland converted to cropland. Our results quantify the role of human activities on the climate footprint of northern wetlands and call for development of active mitigation strategies for managed wetlands and new guidelines of the Intergovernmental Panel on Climate Change (IPCC) accounting for both sustained CH4 emissions and cumulative CO2 exchange.

Joint control of terrestrial gross primary productivity by plant phenology and physiology.

Terrestrial gross primary productivity (GPP) varies greatly over time and space. A better understanding of this variability is necessary for more accurate predictions of the future climate-carbon cycle feedback. Recent studies have suggested that variability in GPP is driven by a broad range of biotic and abiotic factors operating mainly through changes in vegetation phenology and physiological processes. However, it is still unclear how plant phenology and physiology can be integrated to explain the spatiotemporal variability of terrestrial GPP. Based on analyses of eddy-covariance and satellite-derived data, we decomposed annual terrestrial GPP into the length of the CO2 uptake period (CUP) and the seasonal maximal capacity of CO2 uptake (GPPmax). The product of CUP and GPPmax explained >90% of the temporal GPP variability in most areas of North America during 2000-2010 and the spatial GPP variation among globally distributed eddy flux tower sites. It also explained GPP response to the European heatwave in 2003 (r(2) = 0.90) and GPP recovery after a fire disturbance in South Dakota (r(2) = 0.88). Additional analysis of the eddy-covariance flux data shows that the interbiome variation in annual GPP is better explained by that in GPPmax than CUP. These findings indicate that terrestrial GPP is jointly controlled by ecosystem-level plant phenology and photosynthetic capacity, and greater understanding of GPPmax and CUP responses to environmental and biological variations will, thus, improve predictions of GPP over time and space.

Differentiating moss from higher plants is critical in studying the carbon cycle of the boreal biome.

The satellite-derived normalized difference vegetation index (NDVI), which is used for estimating gross primary production (GPP), often includes contributions from both mosses and vascular plants in boreal ecosystems. For the same NDVI, moss can generate only about one-third of the GPP that vascular plants can because of its much lower photosynthetic capacity. Here, based on eddy covariance measurements, we show that the difference in photosynthetic capacity between these two plant functional types has never been explicitly included when estimating regional GPP in the boreal region, resulting in a substantial overestimation. The magnitude of this overestimation could have important implications regarding a change from a current carbon sink to a carbon source in the boreal region. Moss abundance, associated with ecosystem disturbances, needs to be mapped and incorporated into GPP estimates in order to adequately assess the role of the boreal region in the global carbon cycle.