The global potential for increased storage of carbon on land.

Constraining the climate crisis requires urgent action to reduce anthropogenic emissions while simultaneously removing carbon dioxide from the atmosphere. Improved information about the maximum magnitude and spatial distribution of opportunities for additional land-based removals of CO2 is needed to guide on-the-ground decision-making about where to implement climate change mitigation strategies. Here, we present a globally consistent spatial dataset (approximately 500-m resolution) of current, potential, and unrealized potential carbon storage in woody plant biomass and soil organic matter. We also provide a framework for prioritizing actions related to the restoration, management, and maintenance of woody carbon stocks and associated soils. By comparing current to potential carbon storage, while excluding areas critical to food production and human habitation, we find 287 petagrams (PgC) of unrealized potential storage opportunity, of which 78% (224 PgC) is in biomass and 22% (63 PgC) is in soil. Improved management of existing forests may offer nearly three-fourths (206 PgC) of the total unrealized potential, with the majority (71%) concentrated in tropical ecosystems. However, climate change is a source of considerable uncertainty. While additional research is needed to understand the impact of natural disturbances and biophysical feedbacks, we project that the potential for additional carbon storage in woody biomass will increase (+17%) by 2050 despite projected decreases (−12%) in the tropics. Our results establish an absolute reference point and conceptual framework for national and jurisdictional prioritization of locations and actions to increase land-based carbon storage.

Global patterns of drought recovery.

Drought, a recurring phenomenon with major impacts on both human and natural systems, is the most widespread climatic extreme that negatively affects the land carbon sink. Although twentieth-century trends in drought regimes are ambiguous, across many regions more frequent and severe droughts are expected in the twenty-first century. Recovery time-how long an ecosystem requires to revert to its pre-drought functional state-is a critical metric of drought impact. Yet the factors influencing drought recovery and its spatiotemporal patterns at the global scale are largely unknown. Here we analyse three independent datasets of gross primary productivity and show that, across diverse ecosystems, drought recovery times are strongly associated with climate and carbon cycle dynamics, with biodiversity and CO2 fertilization as secondary factors. Our analysis also provides two key insights into the spatiotemporal patterns of drought recovery time: first, that recovery is longest in the tropics and high northern latitudes (both vulnerable areas of Earth's climate system) and second, that drought impacts (assessed using the area of ecosystems actively recovering and time to recovery) have increased over the twentieth century. If droughts become more frequent, as expected, the time between droughts may become shorter than drought recovery time, leading to permanently damaged ecosystems and widespread degradation of the land carbon sink.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

RCP8.5 tracks cumulative CO2 emissions.

Climate simulation-based scenarios are routinely used to characterize a range of plausible climate futures. Despite some recent progress on bending the emissions curve, RCP8.5, the most aggressive scenario in assumed fossil fuel use for global climate models, will continue to serve as a useful tool for quantifying physical climate risk, especially over near- to midterm policy-relevant time horizons. Not only are the emissions consistent with RCP8.5 in close agreement with historical total cumulative CO2 emissions (within 1%), but RCP8.5 is also the best match out to midcentury under current and stated policies with still highly plausible levels of CO2 emissions in 2100.

The terrestrial biosphere as a net source of greenhouse gases to the atmosphere.

The terrestrial biosphere can release or absorb the greenhouse gases, carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O), and therefore has an important role in regulating atmospheric composition and climate. Anthropogenic activities such as land-use change, agriculture and waste management have altered terrestrial biogenic greenhouse gas fluxes, and the resulting increases in methane and nitrous oxide emissions in particular can contribute to climate change. The terrestrial biogenic fluxes of individual greenhouse gases have been studied extensively, but the net biogenic greenhouse gas balance resulting from anthropogenic activities and its effect on the climate system remains uncertain. Here we use bottom-up (inventory, statistical extrapolation of local flux measurements, and process-based modelling) and top-down (atmospheric inversions) approaches to quantify the global net biogenic greenhouse gas balance between 1981 and 2010 resulting from anthropogenic activities and its effect on the climate system. We find that the cumulative warming capacity of concurrent biogenic methane and nitrous oxide emissions is a factor of about two larger than the cooling effect resulting from the global land carbon dioxide uptake from 2001 to 2010. This results in a net positive cumulative impact of the three greenhouse gases on the planetary energy budget, with a best estimate (in petagrams of CO2 equivalent per year) of 3.9 ± 3.8 (top down) and 5.4 ± 4.8 (bottom up) based on the GWP100 metric (global warming potential on a 100-year time horizon). Our findings suggest that a reduction in agricultural methane and nitrous oxide emissions, particularly in Southern Asia, may help mitigate climate change.

Patterns of post-drought recovery are strongly influenced by drought duration, frequency, post-drought wetness, and bioclimatic setting.

Understanding vegetation recovery after drought is critical for projecting vegetation dynamics in future climates. From 1997 to 2009, Australia experienced a long-lasting drought known as the Millennium Drought (MD), which led to widespread reductions in vegetation productivity. However, vegetation recovery post-drought and its determinants remain unclear. This study leverages remote sensing products from different sources-fraction of absorbed photosynthetically active radiation (FPAR), based on optical data, and canopy density, derived from microwave data-and random forest algorithms to assess drought recovery over Australian natural vegetation during a 20-year period centered on the MD. Post-drought recovery was prevalent across the continent, with 6 out of 10 drought events seeing full recovery within about 6 months. Canopy density was slower to recover than leaf area seen in FPAR. The probability of full recovery was most strongly controlled by drought return interval, post-drought hydrological condition, and drought length. Full recovery was seldom observed when drought events occurred at intervals of 3 months or less, and moderately dry (standardized water balance anomaly [SWBA] within [-1, -0.76]) post-drought conditions resulted in less complete recovery than wet (SWBA > 0.3) post-drought conditions. Press droughts, which are long term but not extreme, delayed recovery more than pulse droughts (short term but extreme) and led to a higher frequency of persistent decline. Following press droughts, the frequency of persistent decline differed little among biome types but peaked in semi-arid regions across aridity levels. Forests and savanna required the longest recovery times for press drought, while grasslands were the slowest to recover for pulse drought. This study provides quantitative thresholds that could be used to improve the modeling of ecosystem dynamics post-drought.

Ghosts of the past: how drought legacy effects shape forest functioning and carbon cycling.

Multi-year lags in tree drought recovery, termed 'drought legacy effects', are important for understanding the impacts of drought on forest ecosystems, including carbon (C) cycle feedbacks to climate change. Despite the ubiquity of lags in drought recovery, large uncertainties remain regarding the mechanistic basis of legacy effects and their importance for the C cycle. In this review, we identify the approaches used to study legacy effects, from tree rings to whole forests. We then discuss key knowledge gaps pertaining to the causes of legacy effects, and how the various mechanisms that may contribute these lags in drought recovery could have contrasting implications for the C cycle. Furthermore, we conduct a novel data synthesis and find that legacy effects differ drastically in both size and length across the US depending on if they are identified in tree rings versus gross primary productivity. Finally, we highlight promising approaches for future research to improve our capacity to model legacy effects and predict their impact on forest health. We emphasise that a holistic view of legacy effects - from tissues to whole forests - will advance our understanding of legacy effects and stimulate efforts to investigate drought recovery via experimental, observational and modelling approaches.

Global vegetation biomass production efficiency constrained by models and observations.

Plants use only a fraction of their photosynthetically derived carbon for biomass production (BP). The biomass production efficiency (BPE), defined as the ratio of BP to photosynthesis, and its variation across and within vegetation types is poorly understood, which hinders our capacity to accurately estimate carbon turnover times and carbon sinks. Here, we present a new global estimation of BPE obtained by combining field measurements from 113 sites with 14 carbon cycle models. Our best estimate of global BPE is 0.41 ± 0.05, excluding cropland. The largest BPE is found in boreal forests (0.48 ± 0.06) and the lowest in tropical forests (0.40 ± 0.04). Carbon cycle models overestimate BPE, although models with carbon-nitrogen interactions tend to be more realistic. Using observation-based estimates of global photosynthesis, we quantify the global BP of non-cropland ecosystems of 41 ± 6 Pg C/year. This flux is less than net primary production as it does not contain carbon allocated to symbionts, used for exudates or volatile carbon compound emissions to the atmosphere. Our study reveals a positive bias of 24 ± 11% in the model-estimated BP (10 of 14 models). When correcting models for this bias while leaving modeled carbon turnover times unchanged, we found that the global ecosystem carbon storage change during the last century is decreased by 67% (or 58 Pg C).

Ecological memory of daily carbon exchange across the globe and its importance in drylands.

How do antecedent (past) conditions influence land-carbon dynamics after those conditions no longer persist? In particular, quantifying such memory effects associated with the influence of past environmental (exogenous) and biological (endogenous) conditions is crucial for understanding and predicting the carbon cycle. Here we show, using data from 42 eddy covariance sites across six major biomes, that ecological memory-decomposed into environmental and biological memory components-of daily net carbon exchange (NEE) is critical for understanding the land-carbon metabolism, especially in drylands for which memory explains ~ 32% of the variation in NEE. The strong environmental memory in drylands was primarily driven by short- and long-term moisture status. Moreover, the strength of environmental memory scales with increasing water stress. This universal scaling relationship, emerging within and among major biomes, suggests a potential adaptive response to water limitation. Our findings underscore the necessity of considering ecological memory in experiments, observations and modelling.

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions.

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land-atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO2) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century-long (1901-2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project and two observation-based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 1015 g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr-1 with a median value of 51 Pg C yr-1 during 2001-2010. The largest uncertainty in SOC stocks exists in the 40-65°N latitude whereas the largest cross-model divergence in Rh are in the tropics. The modeled SOC change during 1901-2010 ranges from -70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble-estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO2 and nitrogen deposition over intact ecosystems increased SOC stocks-even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Toward verifying fossil fuel CO2 emissions with the CMAQ model: motivation, model description and initial simulation.

UNLABELLED: Motivated by the question of whether and how a state-of-the-art regional chemical transport model (CTM) can facilitate characterization of CO2 spatiotemporal variability and verify CO2 fossil-fuel emissions, we for the first time applied the Community Multiscale Air Quality (CMAQ) model to simulate CO2. This paper presents methods, input data, and initial results for CO2 simulation using CMAQ over the contiguous United States in October 2007. Modeling experiments have been performed to understand the roles of fossil-fuel emissions, biosphere-atmosphere exchange, and meteorology in regulating the spatial distribution of CO2 near the surface over the contiguous United States. Three sets of net ecosystem exchange (NEE) fluxes were used as input to assess the impact of uncertainty of NEE on CO2 concentrations simulated by CMAQ. Observational data from six tall tower sites across the country were used to evaluate model performance. In particular, at the Boulder Atmospheric Observatory (BAO), a tall tower site that receives urban emissions from Denver CO, the CMAQ model using hourly varying, high-resolution CO2 fossil-fuel emissions from the Vulcan inventory and Carbon Tracker optimized NEE reproduced the observed diurnal profile of CO2 reasonably well but with a low bias in the early morning. The spatial distribution of CO2 was found to correlate with NO(x), SO2, and CO, because of their similar fossil-fuel emission sources and common transport processes. These initial results from CMAQ demonstrate the potential of using a regional CTM to help interpret CO2 observations and understand CO2 variability in space and time. The ability to simulate a full suite of air pollutants in CMAQ will also facilitate investigations of their use as tracers for CO2 source attribution. This work serves as a proof of concept and the foundation for more comprehensive examinations of CO2 spatiotemporal variability and various uncertainties in the future. IMPLICATIONS: Atmospheric CO2 has long been modeled and studied on continental to global scales to understand the global carbon cycle. This work demonstrates the potential of modeling and studying CO2 variability at fine spatiotemporal scales with CMAQ, which has been applied extensively, to study traditionally regulated air pollutants. The abundant observational records of these air pollutants and successful experience in studying and reducing their emissions may be useful for verifying CO2 emissions. Although there remains much more to further investigate, this work opens up a discussion on whether and how to study CO2 as an air pollutant.

Observational and model evidence together support wide-spread exposure to noncompensable heat under continued global warming.

As our planet warms, a critical research question is when and where temperatures will exceed the limits of what the human body can tolerate. Past modeling efforts have investigated the 35°C wet-bulb threshold, proposed as a theoretical upper limit to survivability taking into account physiological and behavioral adaptation. Here, we conduct an extreme value theory analysis of weather station observations and climate model projections to investigate the emergence of an empirically supported heat compensability limit. We show that the hottest parts of the world already experience these heat extremes on a limited basis and that under moderate continued warming parts of every continent, except Antarctica, will see a rapid increase in their extent and frequency. To conclude, we discuss the consequences of the emergence of this noncompensable heat and the need for incorporating different critical thermal limits into heat adaptation planning.

The Terrestrial Biosphere Model Farm.

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.

How close are we to the temperature tipping point of the terrestrial biosphere?

The temperature dependence of global photosynthesis and respiration determine land carbon sink strength. While the land sink currently mitigates ~30% of anthropogenic carbon emissions, it is unclear whether this ecosystem service will persist and, more specifically, what hard temperature limits, if any, regulate carbon uptake. Here, we use the largest continuous carbon flux monitoring network to construct the first observationally derived temperature response curves for global land carbon uptake. We show that the mean temperature of the warmest quarter (3-month period) passed the thermal maximum for photosynthesis during the past decade. At higher temperatures, respiration rates continue to rise in contrast to sharply declining rates of photosynthesis. Under business-as-usual emissions, this divergence elicits a near halving of the land sink strength by as early as 2040.

Author Correction: Modeling suggests fossil fuel emissions have been driving increased land carbon uptake since the turn of the 20th Century.

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

Modeling suggests fossil fuel emissions have been driving increased land carbon uptake since the turn of the 20th Century.

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.

Land carbon models underestimate the severity and duration of drought's impact on plant productivity.

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.

Effects of mutual shading of tree crowns on prediction of photosynthetic light-use efficiency in a coastal Douglas-fir forest.

Gross primary production (GPP) is often expressed as the product of absorbed photosynthetically active radiation and the efficiency (epsilon) with which a plant community uses absorbed radiation in biomass production. Light-use efficiency is affected by environmental stresses, and varies diurnally and seasonally. Uncertainty about epsilon can be a serious limitation when modeling GPP. An important determinant of epsilon is the amount and type of solar radiation incident on a canopy, because an abundance of light can trigger a photo-protective reaction, diminishing GPP. The radiation regime in a forest canopy is determined by the predominant sky conditions and by mutual shading of tree crowns. Shading effects, producing shifts in the amount of incident direct and diffuse solar radiation, have been largely ignored, however, because they depend on forest structure and are difficult to measure. We describe a new approach for estimating changes in mutual canopy shading throughout the day and year based on a canopy structure model derived from light detection and ranging (LiDAR). Proportions of canopy shading were then combined with eddy covariance data to assess the explanatory power for variance in epsilon by regression tree analysis over half-hourly, daily and weekly time scales. The approach explained between 75 and 97% of variance in epsilon, representing an increase of between 5 and 16% compared with models driven solely by meteorological variables.

Accelerating rates of Arctic carbon cycling revealed by long-term atmospheric CO2 measurements.

The contemporary Arctic carbon balance is uncertain, and the potential for a permafrost carbon feedback of anywhere from 50 to 200 petagrams of carbon (Schuur et al., 2015) compromises accurate 21st-century global climate system projections. The 42-year record of atmospheric CO2 measurements at Barrow, Alaska (71.29 N, 156.79 W), reveals significant trends in regional land-surface CO2 anomalies (ΔCO2), indicating long-term changes in seasonal carbon uptake and respiration. Using a carbon balance model constrained by ΔCO2, we find a 13.4% decrease in mean carbon residence time (50% confidence range = 9.2 to 17.6%) in North Slope tundra ecosystems during the past four decades, suggesting a transition toward a boreal carbon cycling regime. Temperature dependencies of respiration and carbon uptake suggest that increases in cold season Arctic labile carbon release will likely continue to exceed increases in net growing season carbon uptake under continued warming trends.