Search

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.

Large datasets of greenhouse gas and energy surface-atmosphere fluxes measured with the eddy-covariance technique (e.g., FLUXNET2015, AmeriFlux BASE) are widely used to benchmark models and remote-sensing products. This study addresses one of the major challenges facing model-data integration: To what spatial extent do flux measurements taken at individual eddy-covariance sites reflect model- or satellite-based grid cells? We evaluate flux footprints—the temporally dynamic source areas that contribute to measured fluxes—and the representativeness of these footprints for target areas (e.g., within 250–3000 m radii around flux towers) that are often used in flux-data synthesis and modeling studies. We examine the land-cover composition and vegetation characteristics, represented here by the Enhanced Vegetation Index (EVI), in the flux footprints and target areas across 214 AmeriFlux sites, and evaluate potential biases as a consequence of the footprint-to-target-area mismatch. Monthly 80% footprint climatologies vary across sites and through time ranging four orders of magnitude from 103 to 107 m2 due to the measurement heights, underlying vegetation- and ground-surface characteristics, wind directions, and turbulent state of the atmosphere. Few eddy-covariance sites are located in a truly homogeneous landscape. Thus, the common model-data integration approaches that use a fixed-extent target area across sites introduce biases on the order of 4%–20% for EVI and 6%–20% for the dominant land cover percentage. These biases are site-specific functions of measurement heights, target area extents, and land-surface characteristics. We advocate that flux datasets need to be used with footprint awareness, especially in research and applications that benchmark against models and data products with explicit spatial information. We propose a simple representativeness index based on our evaluations that can be used as a guide to identify site-periods suitable for specific applications and to provide general guidance for data use.

Leaf-level gas exchange data support the mechanistic understanding of plant fluxes of carbon and water. These fluxes inform our understanding of ecosystem function, are an important constraint on parameterization of terrestrial biosphere models, are necessary to understand the response of plants to global environmental change, and are integral to efforts to improve crop production. Collection of these data using gas analyzers can be both technically challenging and time consuming, and individual studies generally focus on a small range of species, restricted time periods, or limited geographic regions. The high value of these data is exemplified by the many publications that reuse and synthesize gas exchange data, however the lack of metadata and data reporting conventions make full and efficient use of these data difficult. Here we propose a reporting format for leaf-level gas exchange data and metadata to provide guidance to data contributors on how to store data in repositories to maximize their discoverability, facilitate their efficient reuse, and add value to individual datasets. For data users, the reporting format will better allow data repositories to optimize data search and extraction, and more readily integrate similar data into harmonized synthesis products. The reporting format specifies data table variable naming and unit conventions, as well as metadata characterizing experimental conditions and protocols. For common data types that were the focus of this initial version of the reporting format, i.e., survey measurements, dark respiration, carbon dioxide and light response curves, and parameters derived from those measurements, we took a further step of defining required additional data and metadata that would maximize the potential reuse of those data types. To aid data contributors and the development of data ingest tools by data repositories we provided a translation table comparing the outputs of common gas exchange instruments. Extensive consultation with data collectors, data users, instrument manufacturers, and data scientists was undertaken in order to ensure that the reporting format met community needs. The reporting format presented here is intended to form a foundation for future development that will incorporate additional data types and variables as gas exchange systems and measurement approaches advance in the future. The reporting format is published in the U.S. Department of Energy's ESS-DIVE data repository, with documentation and future development efforts being maintained in a version control system.

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.

Large spatial-scale effects of climate extremes on gross primary production (GPP), the largest terrestrial carbon flux, are highly uncertain even as these extremes increase in frequency and extent. Here we report the impacts of spring warming and summer drought in 2012 on GPP across the contiguous United States (CONUS) using estimates from four GPP models: Vegetation Photosynthesis Model (VPM), MOD17A2H V006, Carnegie-Ames-Stanford Approach, and Simple Biosphere/Carnegie-Ames-Stanford Approach. VPM simulations are driven by Moderate Resolution Imaging Spectroradiometer, North American Regional Reanalysis climate data, and C3 and C4 cropland maps from the United States Department of Agriculture Cropland Data Layer data set. Across 25 eddy covariance flux tower sites, GPP estimates from VPM (GPPVPM) showed better accuracy in terms of cross-site variability and interannual variability (R2 = 0.84 and 0.46, respectively) when compared to MOD17 GPP. We further assessed the spatial and temporal (seasonal) consistency between GPP products and the Global Ozone Monitoring Experiment-2 solar-induced chlorophyll fluorescence over CONUS during 2008–2014. The results suggested that GPPVPM agrees best with solar-induced chlorophyll fluorescence across space and time, capturing seasonal dynamics and interannual variations. Anomaly analyses showed that increased GPP during the spring compensated for the reduced GPP during the summer, resulting in near-neutral changes in annual GPP for the CONUS. This study demonstrates the importance of assessing the impacts of different types and timing of climate extremes on GPP and the need to improve light use efficiency models by incorporating C3 and C4 plant functional types.

Aerodynamic canopy height (ha) is the effective height of vegetation canopy for its influence on atmospheric fluxes and is a key parameter of surface-atmosphere coupling. However, methods to estimate ha from data are limited. This synthesis evaluates the applicability and robustness of the calculation of ha from eddy covariance momentum-flux data. At 69 forest sites, annual ha robustly predicted site-to-site and year-to-year differences in canopy heights (R2 = 0.88, 111 site-years). At 23 cropland/grassland sites, weekly ha successfully captured the dynamics of vegetation canopies over growing seasons (R2 > 0.70 in 74 site-years). Our results demonstrate the potential of flux-derived ha determination for tracking the seasonal, interannual, and/or decadal dynamics of vegetation canopies including growth, harvest, land use change, and disturbance. The large-scale and time-varying ha derived from flux networks worldwide provides a new benchmark for regional and global Earth system models and satellite remote sensing of canopy structure.

Light-use efficiency (LUE), which quantifies the plants' efficiency in utilizing solar radiation for photosynthetic carbon fixation, is an important factor for gross primary production estimation. Here we use satellite-based solar-induced chlorophyll fluorescence as a proxy for photosynthetically active radiation absorbed by chlorophyll (APARchl) and derive an estimation of the fraction of APARchl (fPARchl) from four remotely sensed vegetation indicators. By comparing maximum LUE estimated at different scales from 127 eddy flux sites, we found that the maximum daily LUE based on PAR absorption by canopy chlorophyll (εchlmax), unlike other expressions of LUE, tends to converge across biome types. The photosynthetic seasonality in tropical forests can also be tracked by the change of fPARchl, suggesting the corresponding (εchlmax) to have less seasonal variation. This spatio-temporal convergence of LUE derived from fPARchl can be used to build simple but robust gross primary production models and to better constrain process-based models.

Spatially and temporally continuous estimation of plant photosynthetic carbon fixation (or gross primary production, GPP) is crucial to our understanding of the global carbon cycle and the impact of climate change. Besides spatial, seasonal and interannual variations, GPP also exhibits strong diurnal variations. Satellite retrieved solar-induced chlorophyll fluorescence (SIF) provides a spatially continuous, but temporally discrete measurement of plant photosynthesis, and has the potential to be used to estimate GPP at global scale. However, it remains unclear whether the seasonal time series of SIF snapshots taken at a fixed time of the day can be used to infer daily total GPP variation at spatial and seasonal scales. In this study, we first used GPP estimates from 135 eddy covariance flux sites, covering a wide range of geographic locations and biome types, to investigate the relationship between the instantaneous GPP (GPPᵢₙₛₜ) and daily GPP (GPPdₐᵢₗy) on the seasonal course for different times of the day. Latitudinal and diurnal patterns were found to correspond to variations in photosynthetically active radiation (PAR) and light use efficiency (LUE), respectively. We then used the Soil-Canopy Observation Photosynthesis and Energy Balance (SCOPE) model and the FluxCom GPP product to investigate the instantaneous and daily SIF-GPP relationships at five flux tower sites along a latitudinal gradient and at a global scale for different biome types. The results showed that daily SIF had a stronger linear correlation with daily GPP than instantaneous SIF at the seasonal scale, with an instantaneous to daily SIF conversion factor following the latitudinal and seasonal pattern driven by PAR. Our study highlights the necessity to take the latitudinal and diurnal factors into consideration for SIF-GPP relationship analyses or for physiological phenology analyses based on SIF.

Elevated atmospheric CO2 concentrations are expected to enhance photosynthesis and reduce stomatal conductance, thus increasing plant water use efficiency. A recent study based on eddy covariance flux observations from Northern Hemisphere forests showed a large increase in inherent water use efficiency (IWUE). Here we used an updated version of the same data set and robust uncertainty quantification to revisit these contemporary IWUE trends. We tested the hypothesis that the observed IWUE increase could be attributed to interannual trends in plant functional traits, potentially triggered by environmental change. We found that IWUE increased by ~1.3% yr−1, which is less than previously reported but still larger than theoretical expectations. Numerical simulations with the Tethys-Chloris ecosystem model using temporally static plant functional traits cannot explain this increase. Simulations with plant functional trait plasticity, i.e., temporal changes in model parameters such as specific leaf area and maximum Rubisco capacity, match the observed trends in IWUE. Our results show that trends in plant functional traits, equal to 1.0% yr−1, can explain the observed IWUE trends. Thus, at decadal or longer time scales, trait plasticity could potentially influence forest water, carbon, and energy fluxes with profound implications for both the monitoring of temporal changes in plant functional traits and their representation in Earth system models.