Global net land carbon uptake or net biome production (NBP) has increased during recent decades
. Whether its temporal variability and autocorrelation have changed during this period, however, ...remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink
. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO
concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.
Full text
Available for:
GEOZS, IJS, IMTLJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZAGLJ
Earth System Models typically use static responses to temperature to calculate photosynthesis and respiration, but experimental evidence suggests that many plants acclimate to prevailing ...temperatures. We incorporated representations of photosynthetic and leaf respiratory temperature acclimation into the Community Land Model, the terrestrial component of the Community Earth System Model. These processes increased terrestrial carbon pools by 20 Pg C (22%) at the end of the 21st century under a business‐as‐usual (Representative Concentration Pathway 8.5) climate scenario. Including the less certain estimates of stem and root respiration acclimation increased terrestrial carbon pools by an additional 17 Pg C (~40% overall increase). High latitudes gained the most carbon with acclimation, and tropical carbon pools increased least. However, results from both of these regions remain uncertain; few relevant data exist for tropical and boreal plants or for extreme temperatures. Constraining these uncertainties will produce more realistic estimates of land carbon feedbacks throughout the 21st century.
Key Points
Most global land surface models do not account for plant temperature acclimation
With acclimation, future carbon gain was larger in high latitudes and smaller in the tropics
Studies determining acclimation in tropical and boreal plants are needed
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Dry deposition of ozone is an important sink of ozone in near‐surface air. When dry deposition occurs through plant stomata, ozone can injure the plant, altering water and carbon cycling and reducing ...crop yields. Quantifying both stomatal and nonstomatal uptake accurately is relevant for understanding ozone's impact on human health as an air pollutant and on climate as a potent short‐lived greenhouse gas and primary control on the removal of several reactive greenhouse gases and air pollutants. Robust ozone dry deposition estimates require knowledge of the relative importance of individual deposition pathways, but spatiotemporal variability in nonstomatal deposition is poorly understood. Here we integrate understanding of ozone deposition processes by synthesizing research from fields such as atmospheric chemistry, ecology, and meteorology. We critically review methods for measurements and modeling, highlighting the empiricism that underpins modeling and thus the interpretation of observations. Our unprecedented synthesis of knowledge on deposition pathways, particularly soil and leaf cuticles, reveals process understanding not yet included in widely used models. If coordinated with short‐term field intensives, laboratory studies, and mechanistic modeling, measurements from a few long‐term sites would bridge the molecular to ecosystem scales necessary to establish the relative importance of individual deposition pathways and the extent to which they vary in space and time. Our recommended approaches seek to close knowledge gaps that currently limit quantifying the impact of ozone dry deposition on air quality, ecosystems, and climate.
Plain Language Summary
The removal of tropospheric ozone at Earth's surface (often called dry deposition) is important for our understanding of air pollution, ecosystem health, and climate. Several processes contribute to dry deposition of ozone. While we have basic knowledge of these processes, we lack the ability to robustly estimate changes in ozone dry deposition through time and from one place to another. Here we review ozone deposition processes, measurements, and modeling and propose steps necessary to close gaps in understanding. A major conclusion revealed by our review is that most deposition processes can be fairly well described from a theoretical standpoint, but the relative importance of the various processes remains uncertain. We suggest that progress can be made by establishing multiyear measurements of ozone dry deposition at a limited set of sites around the world and coordinating these measurements with laboratory and field experiments that can be integrated with theory through carefully designed modeling studies.
Key Points
Ozone dry deposition through pathways other than plant stomata is critical for describing the total terrestrial ozone sink
Process‐level knowledge of ozone deposition pathways is missing from the models used to quantify deposition impacts on the Earth system
Long‐term ozone flux and related measurements are key for establishing relative importance of individual pathways
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Earth system models (ESMs) rely on the calculation of canopy conductance in land surface models (LSMs) to quantify the partitioning of land surface energy, water, and CO2 fluxes. This is achieved by ...scaling stomatal conductance, gw, determined from physiological models developed for leaves. Traditionally, models for gw have been semi‐empirical, combining physiological functions with empirically determined calibration constants. More recently, optimization theory has been applied to model gw in LSMs under the premise that it has a stronger grounding in physiological theory and might ultimately lead to improved predictive accuracy. However, this premise has not been thoroughly tested. Using original field data from contrasting forest systems, we compare a widely used empirical type and a more recently developed optimization‐type gw model, termed BB and MED, respectively. Overall, we find no difference between the two models when used to simulate gw from photosynthesis data, or leaf gas exchange from a coupled photosynthesis‐conductance model, or gross primary productivity and evapotranspiration for a FLUXNET tower site with the CLM5 community LSM. Field measurements reveal that the key fitted parameters for BB and MED, g1B and g1M, exhibit strong species specificity in magnitude and sensitivity to CO2, and CLM5 simulations reveal that failure to include this sensitivity can result in significant overestimates of evapotranspiration for high‐CO2 scenarios. Further, we show that g1B and g1M can be determined from mean ci/ca (ratio of leaf intercellular to ambient CO2 concentration). Applying this relationship with ci/ca values derived from a leaf δ13C database, we obtain a global distribution of g1B and g1M, and these values correlate significantly with mean annual precipitation. This provides a new methodology for global parameterization of the BB and MED models in LSMs, tied directly to leaf physiology but unconstrained by spatial boundaries separating designated biomes or plant functional types.
We show that current optimization and empirical methods to determine stomatal conductance in Earth system models (ESMs) yield similar results and exhibit similar species‐specific sensitivities to environmental variables, including CO2. Failure to account for different stomatal sensitivities across species or vegetation types will lead to significant errors in ESM simulations, particularly for terrestrial water fluxes. The key fitted parameter calibrating optimization or empirical stomatal conductance models to different global vegetation types may be determined from mean annual precipitation, as shown here for the empirical Ball‐Berry slope, overcoming constraints imposed by inflexibilities in current global calibration methods.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
The divergence among Earth system models in the terrestrial carbon cycle has prompted interest in how to reduce uncertainty. Previous studies have identified model structural uncertainty arising from ...process parameterizations and parameter values. The current study highlights the importance of climate forcing in generating carbon cycle uncertainty. We use simulations in which three models (Community Land Model version 4 (CLM4), CLM4.5, CLM5) with substantially different carbon cycles are forced with two climate reconstructions (CRUNCEPv7, Global Soil Wetness Project 3 version 1 (GSWP3v1)) to examine the contributions of model structure and climate to uncertainty in the carbon cycle over the period 1850–2014. Climate uncertainty for global annual net biome production exceeds one third of total uncertainty (defined as the sum of climate and model structure uncertainty) in the first half of the twentieth century, but declines after the 1950s. Global annual gross primary productivity, net primary productivity, heterotrophic respiration, and vegetation and soil carbon stocks have substantial climate uncertainty (relative to total uncertainty) throughout the simulation period. Climate forcing contributes more than one half of total uncertainty for these carbon cycle fluxes and stocks throughout boreal North America and Eurasia, some midlatitude regions, and in eastern Amazonia and western equatorial Africa during the decade 2000–2009. Comparison with observationally based data sets of the carbon cycle using model benchmarking methods provides insight into strengths and deficiencies among models and climate forcings, but we caution against overreliance on benchmarking to discriminate among models. The conceptualization of uncertainty arising from this study implies embracing multiple feasible model simulations rather than focusing on which model or simulation is best.
Key Points
Two climate reconstructions commonly used to force terrestrial biosphere models over the industrial era differ in climate
Climate forcing produces uncertainty in the carbon cycle, and climate uncertainty exceeds model structural uncertainty in many regions
Models differ in their sensitivity to climate forcing, especially in boreal North America and Eurasia
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory ...and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm/yr (6.56×10^4 cu.km/yr) to 617.1 mm/yr (6.87×10^4 cu.km/yr). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm/sq.yr with a significance level of p<0.05 and 0.38 mm yr−2 with a significance level of p<0.05, respectively). In contrast, the ensemble mean of the LSMs showed no statistically significant change (0.23 mm/sq.yr, p>0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm/yr. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET.
Atmospheric vapor pressure deficit (VPD) is a critical variable in determining plant photosynthesis. Synthesis of four global climate datasets reveals a sharp increase of VPD after the late 1990s. In ...response, the vegetation greening trend indicated by a satellite-derived vegetation index (GIMMS3g), which was evident before the late 1990s, was subsequently stalled or reversed. Terrestrial gross primary production derived from two satellite-based models (revised EC-LUE and MODIS) exhibits persistent and widespread decreases after the late 1990s due to increased VPD, which offset the positive CO
fertilization effect. Six Earth system models have consistently projected continuous increases of VPD throughout the current century. Our results highlight that the impacts of VPD on vegetation growth should be adequately considered to assess ecosystem responses to future climate conditions.
Uncertainty in model initial states produces uncertainty in climate simulations because of unforced variability internal to the climate system. Climate scientists use initial-condition ensembles to ...separate the forced signal of climate change from the unforced internal variability. Our analysis of an 11-member initial-condition ensemble from the Community Earth System Model Version 2 that spans the period 1850-2014 shows that a similar ensemble approach is needed to robustly assess trends in the terrestrial carbon cycle. Uncertainty in model initialization gives rise to internal variability that masks trends in carbon fluxes, and also creates spurious unforced trends, during the period 1960-2014 across North America, meaning that a single model realization can diverge from the observational record or from other models simply because of random behavior. The forced response is, however, evident in the ensemble mean and emerges from the noise of unforced variability at decadal timescales. Our results suggest that trends in the observational record must be interpreted with caution because of multiple possible histories that would have been observed if the sequence of internal variability had unfolded differently. Furthermore, internal variability produces irreducible uncertainty in the carbon cycle, leading to ambiguity in the magnitude and sign of carbon cycle trends, especially at small spatial scales and short timescales. The small spread in initial land carbon pools at 1850 suggests that internal climate variability arising from atmospheric and oceanic initialization, not the biogeochemical initialization, is the predominant cause of carbon cycle variability among ensemble members. Initial-condition ensembles with other Earth system models are needed to develop a multi-model understanding of internal variability in the terrestrial carbon cycle.
Earth system models are intended to make long‐term projections, but they can be evaluated at interannual and seasonal time scales. Although the Community Earth System Model (CESM2) showed ...improvements in a number of terrestrial carbon cycle benchmarks, relative to its predecessor, our analysis suggests that the interannual variability (IAV) in net terrestrial carbon fluxes did not show similar improvements. The model simulated low IAV of net ecosystem production (NEP), resulting in a weaker than observed sensitivity of the carbon cycle to climate variability. Low IAV in net fluxes likely resulted from low variability in gross primary productivity (GPP)—especially in the tropics—and a high covariation between GPP and ecosystem respiration. Although lower than observed, the IAV of NEP had significant climate sensitivities, with positive NEP anomalies associated with warmer and drier conditions in high latitudes, and with wetter and cooler conditions in mid and low latitudes. We identified two dominant modes of seasonal variability in carbon cycle flux anomalies in our fully coupled CESM2 simulations that are characterized by seasonal amplification and redistribution of ecosystem fluxes. Seasonal amplification of net and gross carbon fluxes showed climate sensitivities mirroring those of annual fluxes. Seasonal redistribution of carbon fluxes is initiated by springtime temperature anomalies, but subsequently negative feedbacks in soil moisture during the summer and fall result in net annual carbon losses from land. These modes of variability are also seen in satellite proxies of GPP, suggesting that CESM2 appropriately represents regional sensitivities of photosynthesis to climate variability on seasonal time scales.
Plain Language Summary
Earth system models that are intended to make climate change projections also represent the global exchange of carbon dioxide (CO2) between the atmosphere, oceans, and land. As such, the growth rate and variability of CO2 concentrations in the atmosphere provide a robust measurement to evaluate models. We looked at the interannual variability (IAV) of terrestrial carbon fluxes and their sensitivity to variations in temperature and water that were simulated by the Community Earth System Model and compared them to observations. We found that the model underestimates the IAV of net terrestrial carbon fluxes, especially in the tropics. We also identified two modes of variability that correspond to an increase in summer land carbon uptake (amplification) and a change in the seasonal timing (redistribution) of land carbon fluxes. Notably, seasonal redistribution was initialized by warmer springs that increased early‐season productivity, but subsequent water limitations in the summer and fall resulted in lower‐than‐average productivity over the growing season and net annual C losses to the atmosphere. Similar patterns of seasonal amplification and redistribution are seen in satellite observations of photosynthesis, suggesting that the model is realistically simulating characteristics of terrestrial ecosystems necessary for capturing carbon cycle‐climate feedbacks.
Key Points
Simulated net carbon fluxes show low interannual variability because of low tropical variability and high correlation in component fluxes
Seasonal amplification and redistribution of fluxes characterize dominant modes of carbon cycle variability, consistent with observations
Seasonal redistribution is characterized by warm springs, dry summers, and can result in a net carbon loss from land
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK