Search

Abstract: Terrestrial evapotranspiration is the second‐largest component of the land water cycle, linking the water, energy, and carbon cycles and influencing the productivity and health of ecosystems. The dynamics of ET across a spectrum of spatiotemporal scales and their controls remain an active focus of research across different science disciplines. Here, we provide an overview of the current state of ET science across in situ measurements, partitioning of ET, and remote sensing, and discuss how different approaches complement one another based on their advantages and shortcomings. We aim to facilitate collaboration among a cross‐disciplinary group of ET scientists to overcome the challenges identified in this paper and ultimately advance our integrated understanding of ET.

Many parts of the world face severe and prolonged drought conditions, stressing the sustainability of water resources and agriculture. Transitioning to water-efficient crops is one strategy that can help adapt to water scarcity. An emerging drought-tolerant crop of interest is chia (Salvia hispanica). Yet, no study has compared its large-scale water use dynamics to those of widely established crops across the globe. Here, we use satellite data over multiple years to assess the water use efficiency of chia, alfalfa, corn, and soybean globally. Results show that chia consumed 13-38% less water than alfalfa, corn, and soy and assimilated 14-20% more carbon per amount of water used. Substituting 10% of Southwest United States alfalfa cultivation with chia would save 184 million liters of water per growing season, equivalent to the annual water consumption of 1,300 households. Future research shall explore the economic, societal, and environmental ramifications of substituting alfalfa with chia in dry areas worldwide. These insights can guide decision-makers in promoting sustainable agriculture and water resource management.

A significant challenge in global change research is understanding how vegetation interacts with the environment to influence ecosystem gross primary productivity (GPP) through carbon assimilation. One emerging objective is to consistently predict GPP fluctuations worldwide by establishing a robust scaling relationship between GPP measured at flux towers and satellite spectral reflectance data. However, a major hurdle in achieving this goal is the discrepancy in spatial resolution between early satellite measurements and eddy flux measurements. By using a large set of growing season data covering 100 site-years in North and Central America, we explored the potential of transforming incident and reflected shortwave (Rg) and photosynthetically active radiation (PAR) measurements into a broadband normalized difference vegetation index (NDVI) and near-infrared (NIR) reflectance of vegetation (NIRv) which simultaneously explains the GPP variability. We found that the broadband NDVI and NIRv derived from Rg and PAR measurements at the daily time scale were highly correlated with Planet Fusion, Landsat-8/9, and Sentinel-2 narrowband NDVI and NIRv across a wide range of climate and ecological gradients. The differences between satellite and broadband NDVI and NIRv were found to be significantly associated with soil background variations, phenological stages, water stress and signal saturation of broadband NIR reflectance at high biomass. The seasonal variability of broadband NDVI and NIRv remarkably captured the seasonality of vegetation phenology, evaporative fraction, GPP and rainfall in different ecosystems. Although saturation of GPP at high NDVI was evident, a linear relationship between broadband NIRv times incident PAR versus GPP indicated the effectiveness of NIRv-based approach to capture the hidden light use efficiency impacts on GPP. Our study concludes that inexpensive measurement of Rg and PAR components can provide reliable information on NDVI, NIRv, and GPP uninterruptedly. This enhances the sensing capability of flux tower sites without requiring additional spectrometer measurements. The proposed in-situ vegetation indices make a compelling case on using radiation signals for handshaking between ecosystem-scale measurements and remote sensing observables relevant to carbon uptake.

The temperature sensitivity of ecosystem respiration regulates how the terrestrial carbon sink responds to a warming climate but has been difficult to constrain observationally beyond the plot scale. Here we use observations of atmospheric CO2 concentrations from a network of towers together with carbon flux estimates from state-of-the-art terrestrial biosphere models to characterize the temperature sensitivity of ecosystem respiration, as represented by the Arrhenius activation energy, over various North American biomes. We infer activation energies of 0.43 eV for North America and 0.38 eV to 0.53 eV for major biomes therein, which are substantially below those reported for plot-scale studies (approximately 0.65 eV). This discrepancy suggests that sparse plot-scale observations do not capture the spatial-scale dependence and biome specificity of the temperature sensitivity. We further show that adjusting the apparent temperature sensitivity in model estimates markedly improves their ability to represent observed atmospheric CO2 variability. This study provides observationally constrained estimates of the temperature sensitivity of ecosystem respiration directly at the biome scale and reveals that temperature sensitivities at this scale are lower than those based on earlier plot-scale studies. These findings call for additional work to assess the resilience of large-scale carbon sinks to warming.

Terrestrial open water evaporation is difficult to measure both in situ and remotely yet is critical for understanding changes in reservoirs, lakes, and inland seas from human management and climatically altered hydrological cycling. Multiple satellite missions and data systems (e.g., ECOSTRESS, OpenET) now operationally produce evapotranspiration (ET), but the open water evaporation data produced over millions of water bodies are algorithmically produced differently than the main ET data and are often overlooked in evaluation. Here, we evaluated the open water evaporation algorithm, AquaSEBS, used by ECOSTRESS and OpenET against 19 in situ open water evaporation sites from around the world using MODIS and Landsat data, making this one of the largest open water evaporation validations to date. Overall, our remotely sensed open water evaporation retrieval captured some variability and magnitude in the in situ data when controlling for high wind events (instantaneous: r2 = 0.71; bias = 13% of mean; RMSE = 38% of mean). Much of the instantaneous uncertainty was due to high wind events (u > mean daily 7.5 m·s-1) when the open water evaporation process shifts from radiatively-controlled to atmospherically-controlled; not accounting for high wind events decreases instantaneous accuracy significantly (r2 = 0.47; bias = 36% of mean; RMSE = 62% of mean). However, this sensitivity minimizes with temporal integration (e.g., daily RMSE = 1.2-1.5 mm·day-1). To benchmark AquaSEBS, we ran a suite of 11 machine learning models, but found that they did not significantly improve on the process-based formulation of AquaSEBS suggesting that the remaining error is from a combination of the in situ evaporation measurements, forcing data, and/or scaling mismatch; the machine learning models were able to predict error well in and of itself (r2 = 0.74). Our results provide confidence in the remotely sensed open water evaporation data, though not without uncertainty, and a foundation by which current and future missions may build such operational data.

Drought is often thought to reduce ecosystem photosynthesis. However, theory suggests there is potential for increased photosynthesis during meteorological drought, especially in energy-limited ecosystems. Here, we examine the response of photosynthesis (gross primary productivity, GPP) to meteorological drought across the water-energy limitation spectrum. We find a consistent increase in eddy covariance GPP during spring drought in energy-limited ecosystems (83% of the energy-limited sites). Half of spring GPP sensitivity to precipitation was predicted solely from the wetness index (R2 = 0.47, p 30° N). We then compare these results to terrestrial biosphere model outputs and remote sensing products. In contrast to trends detected in eddy covariance data, model mean GPP always declined under spring precipitation deficits after controlling for air temperature and light availability. While remote sensing products captured the observed negative spring GPP sensitivity in energy-limited ecosystems, terrestrial biosphere models proved insufficiently sensitive to spring precipitation deficits.

Water availability plays a critical role in shaping terrestrial ecosystems, particularly in low- and mid-latitude regions. The sensitivity of vegetation growth to precipitation strongly regulates global vegetation dynamics and their responses to drought, yet sensitivity changes in response to climate change remain poorly understood. Here we use long-term satellite observations combined with a dynamic statistical learning approach to examine changes in the sensitivity of vegetation greenness to precipitation over the past four decades. We observe a robust increase in precipitation sensitivity (0.624% yr-1) for drylands, and a decrease (-0.618% yr-1) for wet regions. Using model simulations, we show that the contrasting trends between dry and wet regions are caused by elevated atmospheric CO2 (eCO2). eCO2 universally decreases the precipitation sensitivity by reducing leaf-level transpiration, particularly in wet regions. However, in drylands, this leaf-level transpiration reduction is overridden at the canopy scale by a large proportional increase in leaf area. The increased sensitivity for global drylands implies a potential decrease in ecosystem stability and greater impacts of droughts in these vulnerable ecosystems under continued global change.

Model Intercomparison Projects (MIPs) are fundamental to our understanding of how the land surface responds to changes in climate. However, MIPs are challenging to conduct, requiring the organization of multiple, decentralized modeling teams throughout the world running common protocols. We explored centralizing these models on a single supercomputing system. We ran nine offline terrestrial biosphere models through the Terrestrial Biosphere Model Farm: CABLE, CENTURY, HyLand, ISAM, JULES, LPJ-GUESS, ORCHIDEE, SiB-3, and SiB-CASA. All models were wrapped in a software framework driven with common forcing data, spin-up, and run protocols specified by the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) for years 1901-2100. We ran more than a dozen model experiments. We identify three major benefits and three major challenges. The benefits include: (a) processing multiple models through a MIP is relatively straightforward, (b) MIP protocols are run consistently across models, which may reduce some model output variability, and (c) unique multimodel experiments can provide novel output for analysis. The challenges are: (a) technological demand is large, particularly for data and output storage and transfer; (b) model versions lag those from the core model development teams; and (c) there is still a need for intellectual input from the core model development teams for insight into model results. A merger with the open-source, cloud-based Predictive Ecosystem Analyzer (PEcAn) ecoinformatics system may be a path forward to overcoming these challenges.

Atmospheric humidity and soil moisture in the Amazon forest are tightly coupled to the region's water balance, or the difference between two moisture fluxes, evapotranspiration minus precipitation (ET-P). However, large and poorly characterized uncertainties in both fluxes, and in their difference, make it challenging to evaluate spatiotemporal variations of water balance and its dependence on ET or P. Here, we show that satellite observations of the HDO/H2O ratio of water vapor are sensitive to spatiotemporal variations of ET-P over the Amazon. When calibrated by basin-scale and mass-balance estimates of ET-P derived from terrestrial water storage and river discharge measurements, the isotopic data demonstrate that rainfall controls wet Amazon water balance variability, but ET becomes important in regulating water balance and its variability in the dry Amazon. Changes in the drivers of ET, such as above ground biomass, could therefore have a larger impact on soil moisture and humidity in the dry (southern and eastern) Amazon relative to the wet Amazon.

Climate change, elevating atmosphere CO2 (eCO2) and increased nitrogen deposition (iNDEP) are altering the biogeochemical interactions between plants, microbes and soils, which further modify plant leaf carbon‑nitrogen (C:N) stoichiometry and their carbon assimilation capability. Many field experiments have observed large sensitivity of leaf C:N ratio to eCO2 and iNDEP. However, the large-scale pattern of this sensitivity is still unclear, because eCO2 and iNDEP drive leaf C:N ratio toward opposite directions, which are further compounded by the complex processes of nitrogen acquisition and plant-and-microbial nitrogen competition. Here, we attempt to map the leaf C:N ratio spatial variation in the past 5 decades in China with a combination of data-driven model and process-based modeling. These two approaches showed consistent results. Over different regions, we found that leaf C:N ratio had significant but uneven changes between 2 time periods (1960-1989 and 1990-2015): a 5% ± 8% increase for temperate grasslands in northern China, a 3% ± 6% increase for boreal grasslands in western China, and by contrast, a 7% ± 6% decrease for temperate forests in southern China, and a 3% ± 5% decrease for boreal forests in northeastern China. Additionally, the structural equation models indicated that the leaf C:N change was sensitive to ΔNDEP, ΔCO2 and ΔMAT rather than ΔMAP and ecosystem types. Process-based modeling suggested that iNDEP was the main source of soil mineral nitrogen change, dominating leaf C:N ratio change in most areas in China, while eCO2 led to leaf C:N ratio increase in low iNDEP area. This study also indicates that the long-term leaf C:N ratio acclimation was dominated by climate constraint, especially temperature, but was constrained by soil N availability over decade scale.

Accurate estimation of gross primary production (GPP), the amount of carbon absorbed by plants via photosynthesis, is of great importance for understanding ecosystem functions, carbon cycling, and climate-carbon feedbacks. Remote sensing has been widely used to quantify GPP at regional to global scales. However, polar-orbiting satellites (e.g., Landsat, Sentinel, Terra, Aqua, Suomi NPP, JPSS, OCO-2) lack the capability to examine the diurnal cycles of GPP because they observe the Earth's surface at the same time of day. The Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS), launched in June 2018, observes the land surface temperature (LST) at different times of day with high spatial resolution (70 m × 70 m) from the International Space Station (ISS). Here, we made use of ECOSTRESS data to predict instantaneous GPP with high spatial resolution for different times of day using a data-driven approach based on machine learning. The predictive GPP model used instantaneous ECOSTRESS LST observations along with the daily enhanced vegetation index (EVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS), land cover type from the National Land Cover Database (NCLD), and instantaneous meteorological data from the ERA5 reanalysis dataset. Our model estimated instantaneous GPP across 56 flux tower sites fairly well (R2 = 0.88, Root Mean Squared Error (RMSE) = 2.42 μmol CO2 m−2 s−1). The instantaneous GPP estimates driven by ECOSTRESS LST captured the diurnal variations of tower GPP for different biomes. We then produced multiple high resolution ECOSTRESS GPP maps for the central and northern California. We found distinct changes in GPP at different times of day (e.g., higher in late morning, peak around noon, approaching zero at dusk), and clear differences in productivity across landscapes (e.g., savannas, croplands, grasslands, and forests) for different times of day. ECOSTRESS GPP also captured the seasonal variations in the diurnal cycling of photosynthesis. This study demonstrates the feasibility of using ECOSTRESS data for producing instantaneous GPP (i.e., GPP for the acquisition time of the ECOSTRESS data) for different times of day. The ECOSTRESS GPP can shed light on how plant photosynthesis and water use vary over the course of the diurnal cycle and inform agricultural management and future improvement of terrestrial biosphere/land surface models.

Abstract: Large uncertainties in North American terrestrial carbon fluxes hinder regional climate projections. Terrestrial biosphere models (TBMs), the essential tools for understanding continental‐scale carbon cycle, diverge on whether temperate forests or croplands dominate carbon uptake in North America. Evidence from novel photosynthetic proxies, such as those based on chlorophyll fluorescence, has cast doubt on the “weak cropland, strong forest” carbon uptake patterns simulated by most TBMs. However, no systematic evaluation of TBMs has yet been attempted to pin down space‐time patterns that are most consistent with regional CO2 observational constraints. Here, we leverage atmospheric CO2 observations and satellite‐observed photosynthetic proxies to understand emergent space‐time patterns in North American carbon fluxes from a large suite of TBMs and data‐driven models. To do so, we evaluate how well the atmospheric signals resulting from carbon flux estimates reproduce the space‐time variability in atmospheric CO2, as is observed by a network of continuous‐monitoring towers over North America. Models with gross or net carbon fluxes that are consistent with the observed CO2 variability share a salient feature of growing‐season carbon uptake in Midwest US croplands. Conversely, the remaining models place most growing‐season uptake in boreal or temperate forests. Differences in model explanatory power depend mainly on the simulated annual cycles of cropland uptake—especially, the timing of peak uptake—rather than the distribution of annual mean fluxes across biomes. Our results suggest that improved model representation of cropland phenology is crucial to robust, policy‐relevant estimation of North American carbon exchange.

The 2017–2027 National Academies' Decadal Survey, Thriving on Our Changing Planet, recommended Surface Biology and Geology (SBG) as a “Designated Targeted Observable” (DO). The SBG DO is based on the need for capabilities to acquire global, high spatial resolution, visible to shortwave infrared (VSWIR; 380–2500 nm; ~30 m pixel resolution) hyperspectral (imaging spectroscopy) and multispectral midwave and thermal infrared (MWIR: 3–5 μm; TIR: 8–12 μm; ~60 m pixel resolution) measurements with sub-monthly temporal revisits over terrestrial, freshwater, and coastal marine habitats. To address the various mission design needs, an SBG Algorithms Working Group of multidisciplinary researchers has been formed to review and evaluate the algorithms applicable to the SBG DO across a wide range of Earth science disciplines, including terrestrial and aquatic ecology, atmospheric science, geology, and hydrology. Here, we summarize current state-of-the-practice VSWIR and TIR algorithms that use airborne or orbital spectral imaging observations to address the SBG DO priorities identified by the Decadal Survey: (i) terrestrial vegetation physiology, functional traits, and health; (ii) inland and coastal aquatic ecosystems physiology, functional traits, and health; (iii) snow and ice accumulation, melting, and albedo; (iv) active surface composition (eruptions, landslides, evolving landscapes, hazard risks); (v) effects of changing land use on surface energy, water, momentum, and carbon fluxes; and (vi) managing agriculture, natural habitats, water use/quality, and urban development. We review existing algorithms in the following categories: snow/ice, aquatic environments, geology, and terrestrial vegetation, and summarize the community-state-of-practice in each category. This effort synthesizes the findings of more than 130 scientists.

In an era of rapid global change, our ability to understand and predict Earth's natural systems is lagging behind our ability to monitor and measure changes in the biosphere. Bottlenecks to informing models with observations have reduced our capacity to fully exploit the growing volume and variety of available data. Here, we take a critical look at the information infrastructure that connects ecosystem modeling and measurement efforts, and propose a roadmap to community cyberinfrastructure development that can reduce the divisions between empirical research and modeling and accelerate the pace of discovery. A new era of data-model integration requires investment in accessible, scalable, and transparent tools that integrate the expertise of the whole community, including both modelers and empiricists. This roadmap focuses on five key opportunities for community tools: the underlying foundations of community cyberinfrastructure; data ingest; calibration of models to data; model-data benchmarking; and data assimilation and ecological forecasting. This community-driven approach is a key to meeting the pressing needs of science and society in the 21st century.

Land surface models range in complexity of terrestrial evapotranspiration, yet it is unknown how model complexity translates to accuracy of modeled evapotranspiration estimates. Here, we use the International Land Model Benchmarking system to assess ET estimates from three models of varying complexity driven by the same forcing datasets: an earth system model, a terrestrial biosphere model, and a stand-alone ET model. The performance assessment includes both temporal and spatial evaluation, and different plant functional types across China. Our results indicate that the most complex model, an earth system model, performed best against the benchmarking datasets and metrics. Terrestrial biosphere model performed best in simulating inter-annual variability of ET, while earth system model performed best in simulating the seasonal cycle. The more complex models (earth system model and terrestrial biosphere model) perform better in forest, shrub and crop ecosystems, while the simpler model (stand-alone ET model) perform better in grass ecosystems. Our study demonstrates the impact of model complexity on ET estimates and highlights directions for future ET model improvements.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Terrestrial vegetation removes CO2 from the atmosphere; an important climate regulation service that slows global warming. This 119 Pg C per annum transfer of CO2 into plants-gross primary productivity (GPP)-is the largest land carbon flux globally. While understanding past and anticipated future GPP changes is necessary to support carbon management, the factors driving long-term changes in GPP are largely unknown. Here we show that 1901 to 2010 changes in GPP have been dominated by anthropogenic activity. Our dual constraint attribution approach provides three insights into the spatiotemporal patterns of GPP change. First, anthropogenic controls on GPP change have increased from 57% (1901 decade) to 94% (2001 decade) of the vegetated land surface. Second, CO2 fertilization and nitro  gen deposition are the most important drivers of change, 19.8 and 11.1 Pg C per annum (2001 decade) respectively, especially in the tropics and industrialized areas since the 1970's. Third, changes in climate have functioned as fertilization to enhance GPP (1.4 Pg C per annum in the 2001 decade). These findings suggest that, from a land carbon balance perspective, the Anthropocene began over 100 years ago and that global change drivers have allowed GPP uptake to keep pace with anthropogenic emissions.

Tropical forests are expected to green up with increasing atmospheric CO2 concentrations, but primary productivity may be limited by soil nutrient availability. However, rarely have canopy-scale measurements been assessed against soil measurements in the tropics. Here, we sought to assess remotely sensed canopy greenness against steep soil nutrient gradients across 50 1-ha mature forest plots in Panama. Contrary to expectations, increases in in situ extractable soil phosphorus (P) and base cations (K, Mg) corresponded to declines in remotely sensed mean annual canopy greenness (r2 = 0.77-0.85; p

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

Terrestrial ecosystems carbon and water cycles are tightly coupled through photosynthesis and evapotranspiration processes. The ratios of carbon stored to carbon uptake and water loss to carbon gain are key ecophysiological indicators essential to assess the magnitude and response of the terrestrial plant to the changing climate. Here, we use estimates from 10 terrestrial ecosystem models to quantify the impacts of climate, atmospheric CO2 concentration, and nitrogen (N) deposition on water use efficiency (WUE), and carbon use efficiency (CUE). We find that across models, WUE increases over the 20th Century particularly due to CO2 fertilization and N deposition and compares favorably to experimental studies. Also, the results show a decrease in WUE with climate for the last 3 decades, in contrasts with up-scaled flux observations that demonstrate a constant WUE. Modeled WUE responds minimally to climate with modeled CUE exhibiting no clear trend across space and time. The divergence between simulated and observationally-constrained WUE and CUE is driven by modeled NPP and autotrophic respiration, nitrogen cycle, carbon allocation, and soil moisture dynamics in current ecosystem models. We suggest that carbon-modeling community needs to reexamine stomatal conductance schemes and the soil-vegetation interactions for more robust modeling of carbon and water cycles.

Clarifying how increased atmospheric CO2 concentration (eCO2) contributes to accelerated land carbon sequestration remains important since this process is the largest negative feedback in the coupled carbon–climate system. Here, we constrain the sensitivity of the terrestrial carbon sink to eCO2 over the temperate Northern Hemisphere for the past five decades, using 12 terrestrial ecosystem models and data from seven CO2 enrichment experiments. This constraint uses the heuristic finding that the northern temperate carbon sink sensitivity to eCO2 is linearly related to the site-scale sensitivity across the models. The emerging data-constrained eCO2 sensitivity is 0.64 ± 0.28 PgC yr−1 per hundred ppm of eCO2. Extrapolating worldwide, this northern temperate sensitivity projects the global terrestrial carbon sink to increase by 3.5 ± 1.9 PgC yr−1 for an increase in CO2 of 100 ppm. This value suggests that CO2 fertilization alone explains most of the observed increase in global land carbon sink since the 1960s. More CO2 enrichment experiments, particularly in boreal, arctic and tropical ecosystems, are required to explain further the responsible processes.

Elevated CO2 (eCO2) experiments provide critical information to quantify the effects of rising CO2 on vegetation1–6. Many eCO2 experiments suggest that nutrient limitations modulate the local magnitude of the eCO2 effect on plant biomass1,3,5, but the global extent of these limitations has not been empirically quantified, complicating projections of the capacity of plants to take up CO27,8. Here, we present a data-driven global quantification of the eCO2 effect on biomass based on 138 eCO2 experiments. The strength of CO2 fertilization is primarily driven by nitrogen (N) in ~65% of global vegetation and by phosphorus (P) in ~25% of global vegetation, with N- or P-limitation modulated by mycorrhizal association. Our approach suggests that CO2 levels expected by 2100 can potentially enhance plant biomass by 12 ± 3% above current values, equivalent to 59 ± 13 PgC. The global-scale response to eCO2 we derive from experiments is similar to past changes in greenness9 and biomass10 with rising CO2, suggesting that CO2 will continue to stimulate plant biomass in the future despite the constraining effect of soil nutrients. Our research reconciles conflicting evidence on CO2 fertilization across scales and provides an empirical estimate of the biomass sensitivity to eCO2 that may help to constrain climate projections.

Characterizing the terrestrial carbon, water and energy cycles depends strongly on a capacity to accurately reproduce the spatial and temporal dynamics of land surface evaporation. For this, and many other reasons, monitoring terrestrial evaporation across multiple space and time scales has been an area of focused research for many decades. Much of this activity has been supported by developments in satellite remote sensing, which have been leveraged to deliver new process insights, model development and methodological improvements. In this Special Issue, published contributions explored a range of research topics directed towards the enhanced estimation of terrestrial evaporation. Here we summarize these cutting-edge efforts and provide an overview of some of the state-of-the-art approaches for retrieving this key variable. Some perspectives on outstanding challenges, issues, and opportunities are also presented.

The ability to accurately predict ecosystem drought response and recovery is necessary to produce reliable forecasts of land carbon uptake and future climate. Using a suite of models from the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), we assessed modeled net primary productivity (NPP) response to, and recovery from, drought events against a benchmark derived from tree ring observations between 1948 and 2008 across forested regions of the US and Europe. We find short lag times (0-6 months) between climate anomalies and modeled NPP response. Although models accurately simulate the direction of drought legacy effects (i.e. NPP decreases), projected effects are approximately four times shorter and four times weaker than observations suggest. This discrepancy between observed and simulated vegetation recovery from drought reveals a potential critical model deficiency. Since productivity is a crucial component of the land carbon balance, models that underestimate drought recovery time could overestimate predictions of future land carbon sink strength and, consequently, underestimate forecasts of atmospheric CO2.

Vegetation impacts on ecosystem functioning are mediated by mycorrhizas, plant-fungal associations formed by most plant species. Ecosystems dominated by distinct mycorrhizal types differ strongly in their biogeochemistry. Quantitative analyses of mycorrhizal impacts on ecosystem functioning are hindered by the scarcity of information on mycorrhizal distributions. Here we present global, high-resolution maps of vegetation biomass distribution by dominant mycorrhizal associations. Arbuscular, ectomycorrhizal, and ericoid mycorrhizal vegetation store, respectively, 241 ± 15, 100 ± 17, and 7 ± 1.8 GT carbon in aboveground biomass, whereas non-mycorrhizal vegetation stores 29 ± 5.5 GT carbon. Soil carbon stocks in both topsoil and subsoil are positively related to the community-level biomass fraction of ectomycorrhizal plants, though the strength of this relationship varies across biomes. We show that human-induced transformations of Earth's ecosystems have reduced ectomycorrhizal vegetation, with potential ramifications to terrestrial carbon stocks. Our work provides a benchmark for spatially explicit and globally quantitative assessments of mycorrhizal impacts on ecosystem functioning and biogeochemical cycling.

The El Niño-Southern Oscillation exerts a large influence on global climate regimes and on the global carbon cycle. Although El Niño is known to be associated with a reduction of the global total land carbon sink, results based on prognostic models or measurements disagree over the relative contribution of photosynthesis to the reduced sink. Here, we provide an independent remote sensing-based analysis on the impact of the 2015-2016 El Niño on global photosynthesis using six global satellite-based photosynthesis products and a global solar-induced fluorescence (SIF) dataset. An ensemble of satellite-based photosynthesis products showed a negative anomaly of -0.7 ± 1.2 PgC in 2015, but a slight positive anomaly of 0.05 ± 0.89 PgC in 2016, which when combined with observations of the growth rate of atmospheric carbon dioxide concentrations suggests that the reduction of the land residual sink was likely dominated by photosynthesis in 2015 but by respiration in 2016. The six satellite-based products unanimously identified a major photosynthesis reduction of -1.1 ± 0.52 PgC from savannahs in 2015 and 2016, followed by a highly uncertain reduction of -0.22 ± 0.98 PgC from rainforests. Vegetation in the Northern Hemisphere enhanced photosynthesis before and after the peak El Niño, especially in grasslands (0.33 ± 0.13 PgC). The patterns of satellite-based photosynthesis ensemble mean were corroborated by SIF, except in rainforests and South America, where the anomalies of satellite-based photosynthesis products also diverged the most. We found the inter-model variation of photosynthesis estimates was strongly related to the discrepancy between moisture forcings for models. These results highlight the importance of considering multiple photosynthesis proxies when assessing responses to climatic anomalies.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'.

Ten years of terrestrial water storage anomalies from the Gravity Recovery and Climate Experiment (GRACE) were used to estimate high latitude snowfall accumulation using a mass balance approach. The estimates were used to assess two common gauge-undercatch correction factors (CFs): Legates climatology (CF-L) utilized in the Global Precipitation Climatology Project (GPCP), and Fuchs dynamic correction model (CF-F) used in the Global Precipitation Climatology Centre (GPCC) Monitoring product. The two CFs can be different by more than 50%. CF-L tended to exceed CF-F over northern Asia and Eurasia, while the opposite was observed over North America. Estimates of snowfall from GPCP, GPCC-L (GPCC corrected by CF-L), and GPCC-F (GPCC corrected by CF-F) were 62%, 64%, and 46% more than GPCC over northern Asia and Eurasia. GRACE-based estimate (49% more than GPCC) was the closest to GPCC-F. We found that as near surface air temperature decreases, the products increasingly underestimated the GRACE-based snowfall accumulation. Overall, GRACE showed that CFs are effective in improving GPCC estimates. Furthermore, our case studies and overall statistics suggest that CF-F is likely more effective than CF-L in most of the high latitude regions studied here. GPCP showed generally better skill than GPCC-L, which might be related to the use of satellite data or additional quality controls on gauge inputs to GPCP. This study suggests that GPCP can be improved if it employs CF-L instead of CF-F to correct for gauge undercatch. However, this implementation requires further studies, region-specific analysis, and operational considerations.

Terrestrial carbon cycle models have incorporated increasingly more processes as a means to achieve more-realistic representations of ecosystem carbon cycling. Despite this, there are large across-model variations in the simulation and projection of carbon cycling. Several model intercomparison projects (MIPs), for example, the fifth phase of the Coupled Model Intercomparison Project (CMIP5) (historical simulations), Trends in Net Land-Atmosphere Carbon Exchange (TRENDY), and Multiscale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), have sought to understand intermodel differences. In this study, the authors developed a suite of new techniques to conduct post-MIP analysis to gain insights into uncertainty sources across 25 models in the three MIPs. First, terrestrial carbon storage dynamics were characterized by a three-dimensional (3D) model output space with coordinates of carbon residence time, net primary productivity (NPP), and carbon storage potential. The latter represents the potential of an ecosystem to lose or gain carbon. This space can be used to measure how and why model output differs. Models with a nitrogen cycle generally exhibit lower annual NPP in comparison with other models, and mostly negative carbon storage potential. Second, a transient traceability framework was used to decompose any given carbon cycle model into traceable components and identify the sources of model differences. The carbon residence time (or NPP) was traced to baseline carbon residence time (or baseline NPP related to the maximum carbon input), environmental scalars, and climate forcing. Third, by applying a variance decomposition method, the authors show that the intermodel differences in carbon storage can be mainly attributed to the baseline carbon residence time and baseline NPP (>90% in the three MIPs). The three techniques developed in this study offer a novel approach to gain more insight from existing MIPs and can point out directions for future MIPs. Since this study is conducted at the global scale for an overview on intermodel differences, future studies should focus more on regional analysis to identify the sources of uncertainties and improve models at the specified mechanism level.

The terrestrial carbon and water cycles are coupled through a multitude of connected processes among soil, roots, leaves, and the atmosphere. The strength and sensitivity of these couplings are not yet well known at the global scale, which contributes to uncertainty in predicting the terrestrial water and carbon budgets. We now have synchronous, global-scale satellite observations of critical terrestrial carbon and water cycle components: solar-induced chlorophyll fluorescence (SIF) and soil moisture. We used these observations within the framework of a global terrestrial biosphere model (Simplified Simple Biosphere Model version 2.0, SSiB2) to investigate carbon-water coupling processes. We updated SSiB2 to include a mechanistic representation of SIF and tested the sensitivity of model parameters to improve the simulation of both SIF and soil moisture with the ultimate objective of improving the first-order terrestrial carbon component, gross primary production. Although several vegetation parameters, such as leaf area index and the green leaf fraction, improved the simulated SIF, and several soil parameters, such as hydraulic conductivity, improved simulated soil moisture, their effects were mainly limited to their respective cycles. One root-mean-square error parameter emerged as the key coupler between the carbon and water cycles: the wilting point. Updates to the wilting point significantly improved the simulations for SIF and gross primary production although substantial mismatches with the satellite data still existed. This study demonstrates the value of synchronous global measurements of the terrestrial carbon and water cycles in improving the understanding of coupled carbon-water cycles.

NASA has launched the decade-long Arctic-Boreal Vulnerability Experiment (ABoVE). While the initial phases focus on field and airborne data collection, early integration with modeling activities is important to benefit future modeling syntheses. We compiled feedback from ecosystem modeling teams on key data needs, which encompass carbon biogeochemistry, vegetation, permafrost, hydrology, and disturbance dynamics. A suite of variables was identified as part of this activity with a critical requirement that they are collected concurrently and representatively over space and time. Individual projects in ABoVE may not capture all these needs, and thus there is both demand and opportunity for the augmentation of field observations, and synthesis of the observations that are collected, to ensure that science questions and integrated modeling activities are successfully implemented.

Water use efficiency (WUE), defined as the ratio of gross primary productivity and evapotranspiration at the ecosystem scale, is a critical variable linking the carbon and water cycles. Incorporating a dependency on vapor pressure deficit, apparent underlying WUE (uWUE) provides a better indicator of how terrestrial ecosystems respond to environmental changes than other WUE formulations. Here we used 20th century simulations from four terrestrial biosphere models to develop a novel variance decomposition method. With this method, we attributed variations in apparent uWUE to both the trend and interannual variation of environmental drivers. The secular increase in atmospheric CO2 explained a clear majority of total variation (66 ± 32%: mean ± one standard deviation), followed by positive trends in nitrogen deposition and climate, as well as a negative trend in land use change. In contrast, interannual variation was mostly driven by interannual climate variability. To analyze the mechanism of the CO2 effect, we partitioned the apparent uWUE into the transpiration ratio (transpiration over evapotranspiration) and potential uWUE. The relative increase in potential uWUE parallels that of CO2, but this direct CO2 effect was offset by 20 ± 4% by changes in ecosystem structure, that is, leaf area index for different vegetation types. However, the decrease in transpiration due to stomatal closure with rising CO2 was reduced by 84% by an increase in leaf area index, resulting in small changes in the transpiration ratio. CO2 concentration thus plays a dominant role in driving apparent uWUE variations over time, but its role differs for the two constituent components: potential uWUE and transpiration.

The Regional Hydrologic Extremes Assessment System (RHEAS) is a prototype software framework for hydrologic modeling and data assimilation that automates the deployment of water resources nowcasting and forecasting applications. A spatially-enabled database is a key component of the software that can ingest a suite of satellite and model datasets while facilitating the interfacing with Geographic Information System (GIS) applications. The datasets ingested are obtained from numerous space-borne sensors and represent multiple components of the water cycle. The object-oriented design of the software allows for modularity and extensibility, showcased here with the coupling of the core hydrologic model with a crop growth model. RHEAS can exploit multi-threading to scale with increasing number of processors, while the database allows delivery of data products and associated uncertainty through a variety of GIS platforms. A set of three example implementations of RHEAS in the United States and Kenya are described to demonstrate the different features of the system in real-world applications.