Search

Live fuel moisture content (LFMC) plays a critical role in wildfire dynamics, but little is known about responses of LFMC to multivariate climate change, e.g., warming temperature, CO2 fertilization, and altered precipitation patterns, leading to a limited prediction ability of future wildfire risks. Here, we use a hydrodynamic demographic vegetation model to estimate LFMC dynamics of chaparral shrubs, a dominant vegetation type in fire-prone southern California. We parameterize the model based on observed shrub allometry and hydraulic traits and evaluate the model's accuracy through comparisons between observed and simulated LFMC of three plant functional types (PFTs) under current climate conditions. Moreover, we estimate the number of days per year of LFMC below 79% (which is a critical threshold for wildfire danger rating of southern California chaparral shrubs) from 1960 to 2099 for each PFT and compare the number of days below the threshold for medium and high greenhouse gas emission scenarios (RCP4.5 and 8.5). We find that climate change could lead to more days per year (5.2%-14.8% increase) with LFMC below 79% between the historical (1960-1999) and future (2080-2099) periods, implying an increase in wildfire danger for chaparral shrubs in southern California. Under the high greenhouse gas emission scenario during the dry season, we find that the future LFMC reductions mainly result from a warming temperature, which leads to 9.1%-18.6% reduction in LFMC. Lower precipitation in the spring leads to a 6.3%-8.1% reduction in LFMC. The combined impacts of warming and precipitation change on fire season length are equal to the additive impacts of warming and precipitation change individually. Our results show that the CO2 fertilization will mitigate fire risk by causing a 3.5%-4.8% increase in LFMC. Our results suggest that multivariate climate change could cause a significant net reduction in LFMC and thus exacerbate future wildfire danger in chaparral shrub systems.

Water deficit in the atmosphere and soil are two key interactive factors that constrain transpiration and vegetation productivity. It is not clear which of these two factors is more important for the water and carbon flux response to drought stress in ecosystems. In this study, field data and numerical modeling were used to isolate their impact on evapotranspiration (ET) and gross primary productivity (GPP) at a tropical forest site in Barro Colorado Island (BCI), Panama, focusing on their response to the drought induced by the El Niño event of 2015–2016. Numerical simulations were performed using a plant hydrodynamic scheme (HYDRO) and a heuristic approach that ignores stomatal sensitivity to leaf water potential in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). The sensitivity of canopy conductance (Gs) to vapor pressure deficit (VPD) obtained from eddy-covariance fluxes and measured sap flux shows that, at both ecosystem and plant scale, soil water stress is more important in limiting Gs than VPD at BCI during the El Niño event. The model simulations confirmed the importance of water stress limitation on Gs, but overestimated the VPD impact on Gs compared to that estimated from the observations. We also found that the predicted soil moisture is less sensitive to the diversity of plant hydraulic traits than ET and GPP. During the dry season at BCI, seasonal ET, especially soil evaporation at VPD > 0.42 kPa, simulated using HYDRO and ELM, were too strong and will require alternative parameterizations.

Tropical forest degradation from logging, fire, and fragmentation not only alters carbon stocks and carbon fluxes, but also impacts physical land surface properties such as albedo and roughness length. Such impacts are poorly quantified to date due to difficulties in accessing and maintaining observational infrastructures, as well as the lack of proper modeling tools for capturing the interactions among biophysical properties, ecosystem demography, canopy structure, and biogeochemical cycling in tropical forests. As a first step to address these limitations, we implemented a selective logging module into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) by mimicking the ecological, biophysical, and biogeochemical processes following a logging event. The model can specify the timing and aerial extent of logging events, splitting the logged forest patch into disturbed and intact patches; determine the survivorship of cohorts in the disturbed patch; and modifying the biomass and necromass (total mass of coarse woody debris and litter) pools following logging. We parameterized the logging module to reproduce a selective logging experiment at the Tapajós National Forest in Brazil and benchmarked model outputs against available field measurements. Our results suggest that the model permits the coexistence of early and late successional functional types and realistically characterizes the seasonality of water and carbon fluxes and stocks, the forest structure and composition, and the ecosystem succession following disturbance. However, the current version of FATES overestimates water stress in the dry season and therefore fails to capture seasonal variation in latent and sensible heat fluxes. Moreover, we observed a bias towards low stem density and leaf area when compared to observations, suggesting that improvements are needed in both carbon allocation and establishment of trees. The effects of logging were assessed by different logging scenarios to represent reduced impact and conventional logging practices, both with high and low logging intensities. The model simulations suggest that in comparison to old-growth forests the logged forests rapidly recover water and energy fluxes in 1 to 3 years. In contrast, the recovery times for carbon stocks, forest structure, and composition are more than 30 years depending on logging practices and intensity. This study lays the foundation to simulate land use change and forest degradation in FATES, which will be an effective tool to directly represent forest management practices and regeneration in the context of Earth system models.

Results from the fully and biogeochemically coupled simulations in which CO2 increases at a rate of 1%yr-1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon-concentration and carbon-climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77±0.37 ° C EgC-1 and is similar to that found in CMIP5 models (1.63±0.48 °C EgC-1) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

Plant functional traits determine vegetation responses to environmental variation, but variation in trait values is large, even within a single site. Likewise, uncertainty in how these traits map to Earth system feedbacks is large. We use a vegetation demographic model (VDM), the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to explore parameter sensitivity of model predictions, and comparison to observations, at a tropical forest site: Barro Colorado Island in Panama. We define a single 12-dimensional distribution of plant trait variation, derived primarily from observations in Panama, and define plant functional types (PFTs) as random draws from this distribution. We compare several model ensembles, where individual ensemble members vary only in the plant traits that define PFTs, and separate ensembles differ from each other based on either model structural assumptions or non-trait, ecosystem-level parameters, which include (a) the number of competing PFTs present in any simulation and (b) parameters that govern disturbance and height-based light competition. While singlePFT simulations are roughly consistent with observations of productivity at Barro Colorado Island, increasing the number of competing PFTs strongly shifts model predictions towards higher productivity and biomass forests. Different ecosystem variables show greater sensitivity than others to the number of competing PFTs, with the predictions that are most dominated by large trees, such as biomass, being the most sensitive. Changing disturbance and height-sorting parameters, i.e., the rules of competitive trait filtering, shifts regimes of dominance or coexistence between early-and late-successional PFTs in the model. Increases to the extent or severity of disturbance, or to the degree of determinism in height-based light competition, all act to shift the community towards early-successional PFTs. In turn, these shifts in competitive outcomes alter predictions of ecosystem states and fluxes, with more early-successional-dominated forests having lower biomass. It is thus crucial to differentiate between plant traits, which are under competitive pressure in VDMs, from those model parameters that are not and to better understand the relationships between these two types of model parameters to quantify sources of uncertainty in VDMs.

Land surface models (LSMs) are a vital tool for understanding, projecting, and predicting the dynamics of the land surface and its role within the Earth system, under global change. Driven by the need to address a set of key questions, LSMs have grown in complexity from simplified representations of land surface biophysics to encompass a broad set of interrelated processes spanning the disciplines of biophysics, biogeochemistry, hydrology, ecosystem ecology, community ecology, human management, and societal impacts. This vast scope and complexity, while warranted by the problems LSMs are designed to solve, has led to enormous challenges in understanding and attributing differences between LSM predictions. Meanwhile, the wide range of spatial scales that govern land surface heterogeneity, and the broad spectrum of timescales in land surface dynamics, create challenges in tractably representing processes in LSMs. We identify three “grand challenges” in the development and use of LSMs, based around these issues: managing process complexity, representing land surface heterogeneity, and understanding parametric dynamics across the broad set of problems asked of LSMs in a changing world. In this review, we discuss progress that has been made, as well as promising directions forward, for each of these challenges.

There is large uncertainty whether Amazon forests will remain a carbon sink as atmospheric CO2 increases. Hence, we simulated an old-growth tropical forest using six versions of four terrestrial models differing in scale of vegetation structure and representation of biogeochemical (BGC) cycling, all driven with CO2 forcing from the preindustrial period to 2100. The models were benchmarked against tree inventory and eddy covariance data from a Brazilian site for present-day predictions. All models predicted positive vegetation growth that outpaced mortality, leading to continual increases in present-day biomass accumulation. Notably, the two vegetation demographic models (VDMs) (ED2 and ELM-FATES) always predicted positive stem diameter growth in all size classes. The field data, however, indicated that a quarter of canopy trees didn't grow over the 15-year period, and while high interannual variation existed, biomass change was near neutral. With a doubling of CO2, three of the four models predicted an appreciable biomass sink (0.77 to 1.24 Mg ha−1 year−1). ELMv1-ECA, the only model used here that includes phosphorus constraints, predicted the lowest biomass sink relative to initial biomass stocks (+21%), lower than the other BGC model, CLM5 (+48%). Models projections differed primarily through variations in nutrient constraints, then carbon allocation, initial biomass, and density-dependent mortality. The VDM's performance was similar or better than the BGC models run in carbon-only mode, suggesting that nutrient competition in VDMs will improve predictions. We demonstrate that VDMs are comparable to nondemographic (i.e., “big-leaf”) models but also include finer scale demography and competition that can be evaluated against field observations.

Despite the fundamental importance of soil temperature for Earth's carbon and energy budgets, ecosystem functioning, and agricultural production, studies of climate change impacts on soil processes have mainly relied on air temperatures, assuming they are accurate proxies for soil temperatures. We evaluated changes in soil temperature, moisture, and air temperature predicted over the 21st century from 14 Earth system models. The model ensemble predicted a global mean soil warming of 2.3 ± 0.7 and 4.5 ± 1.1 °C at 100-cm depth by the end of the 21st century for RCPs 4.5 and 8.5, respectively. Soils at 100 cm warmed at almost exactly the same rate as near-surface (~1 cm) soils. Globally, soil warming was slightly slower than air warming above it, and this difference increased over the 21st century. Regionally, soil warming kept pace with air warming in tropical and arid regions but lagged air warming in colder regions. Thus, air warming is not necessarily a good proxy for soil warming in cold regions where snow and ice impede the direct transfer of sensible heat from the atmosphere to soil. Despite this effect, high-latitude soils were still projected to warm faster than elsewhere, albeit at slower rates than surface air above them. When compared with observations, the models were able to capture soil thermal dynamics in most biomes, but some failed to recreate thermal properties in permafrost regions. Particularly in cold regions, using soil warming rather than air warming projections may improve predictions of temperature-sensitive soil processes.

Key uncertainties in terrestrial carbon cycle projections revolve around the timing, direction, and magnitude of the carbon cycle feedback to climate change. This is especially true in carbon-rich Arctic ecosystems, where permafrost soils contain roughly one third of the world's soil carbon stocks, which are likely vulnerable to loss. Using an ensemble of soil biogeochemical models that reflect recent changes in the conceptual understanding of factors responsible for soil carbon persistence, we quantify potential soil carbon responses under two representative climate change scenarios. Our results illustrate that models disagree on the sign and magnitude of global soil changes through 2100, with disagreements primarily driven by divergent responses of Arctic systems. These results largely reflect different assumptions about the nature of soil carbon persistence and vulnerabilities, underscoring the challenges associated with setting allowable greenhouse gas emission targets that will limit global warming to 1.5°C.

Vegetation tolerance to drought depends on an array of site-specific environmental and plant physiological factors. This tolerance is poorly understood for many forest types despite its importance for predicting and managing vegetation stress. We analyzed the relationships between precipitation variability and forest die-off in California's Sierra Nevada and introduce a new measure of drought tolerance that emphasizes plant access to subsurface moisture buffers. We applied this metric to California's severe 2012-2015 drought, and show that it predicted the patterns of tree mortality. We then examined future climate scenarios, and found that the probability of droughts that lead to widespread die-off increases threefold by the end of the 21st century. Our analysis shows that tree mortality in the Sierra Nevada will likely accelerate in the coming decades and that forests in the Central and Northern Sierra Nevada that largely escaped mortality in 2012-2015 are vulnerable to die-off.

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.

Vegetation plays an important role in regulating global carbon cycles and is a key component of the Earth system models (ESMs) that aim to project Earth's future climate. In the last decade, the vegetation component within ESMs has witnessed great progress from simple "big-leaf" approaches to demographically structured approaches, which have a better representation of plant size, canopy structure, and disturbances. These demographically structured vegetation models typically have a large number of input parameters, and sensitivity analysis is needed to quantify the impact of each parameter on the model outputs for a better understanding of model behavior. In this study, we conducted a comprehensive sensitivity analysis to diagnose the Community Land Model coupled to the Functionally Assembled Terrestrial Simulator, or CLM4.5(FATES). Specifically, we quantified the first- and second-order sensitivities of the model parameters to outputs that represent simulated growth and mortality as well as carbon fluxes and stocks for a tropical site with an extent of 1×1°. While the photosynthetic capacity parameter (Vc;max25) is found to be important for simulated carbon stocks and fluxes, we also show the importance of carbon storage and allometry parameters, which determine survival and growth strategies within the model. The parameter sensitivity changes with different sizes of trees and climate conditions. The results of this study highlight the importance of understanding the dynamics of the next generation of demographically enabled vegetation models within ESMs to improve model parameterization and structure for better model fidelity.

Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (

Climate change is expected to increase the frequency of flooding events and, thus, the risks of flood-related mortality and infrastructure damage. Global-scale assessments of future flooding from Earth system models based only on precipitation changes neglect important processes that occur within the land surface, particularly plant physiological responses to rising CO2. Higher CO2 can reduce stomatal conductance and transpiration, which may lead to increased soil moisture and runoff in some regions, promoting flooding even without changes in precipitation. Here we assess the relative impacts of plant physiological and radiative greenhouse effects on changes in daily runoff intensity over tropical continents using the Community Earth System Model. We find that extreme percentile rates increase significantly more than mean runoff in response to higher CO2. Plant physiological effects have a small impact on precipitation intensity but are a dominant driver of runoff intensification, contributing to one half of the 99th and one third of the 99.9th percentile runoff intensity changes.

The increasing complexity of Earth system models has inspired efforts to quantitatively assess model fidelity through rigorous comparison with best available measurements and observational data products. Earth system models exhibit a high degree of spread in predictions of land biogeochemistry, biogeophysics, and hydrology, which are sensitive to forcing from other model components. Based on insights from prior land model evaluation studies and community workshops, the authors developed an open source model benchmarking software package that generates graphical diagnostics and scores model performance in support of the International Land Model Benchmarking (ILAMB) project. Employing a suite of in situ, remote sensing, and reanalysis data sets, the ILAMB package performs comprehensive model assessment across a wide range of land variables and generates a hierarchical set of web pages containing statistical analyses and figures designed to provide the user insights into strengths and weaknesses of multiple models or model versions. Described here is the benchmarking philosophy and mathematical methodology embodied in the most recent implementation of the ILAMB package. Comparison methods unique to a few specific data sets are presented, and guidelines for configuring an ILAMB analysis and interpreting resulting model performance scores are discussed. ILAMB is being adopted by modeling teams and centers during model development and for model intercomparison projects, and community engagement is sought for extending evaluation metrics and adding new observational data sets to the benchmarking framework.

In the last few decades, temperatures in the Arctic have increased twice as much as the rest of the globe. As permafrost thaws in response to this warming, large amounts of soil organic matter may become vulnerable to decomposition. Microbial decomposition will release carbon (C) from permafrost soils, however, warmer conditions could also lead to enhanced plant growth and C uptake. Field and modeling studies show high uncertainty in soil and plant responses to climate change but there have been few studies that reconcile field and model data to understand differences and reduce uncertainty. Here, we evaluate gross primary productivity (GPP), ecosystem respiration (Reco), and net ecosystem C exchange (NEE) from eight years of experimental soil warming in moist acidic tundra against equivalent fluxes from the Community Land Model during simulations parameterized to reflect the field conditions associated with this manipulative field experiment. Over the eight-year experimental period, soil temperatures and thaw depths increased with warming in field observations and model simulations. However, the field and model results do not agree on warming effects on water table depth; warming created wetter soils in the field and drier soils in the models. In the field, initial increases in growing season GPP, Reco, and NEE to experimentally-induced permafrost thaw created a higher C sink capacity in the first years followed by a stronger C source in years six through eight. In contrast, both models predicted linear increases in GPP, Reco, and NEE with warming. The divergence of model results from field experiments reveals the role subsidence, hydrology, and nutrient cycling play in influencing the C flux responses to permafrost thaw, a complexity that the models are not structurally able to predict, and highlight challenges associated with projecting C cycle dynamics across the Arctic.

Tree mortality is a key driver of forest community composition and carbon dynamics. Strong winds associated with severe convective storms are dominant natural drivers of tree mortality in the Amazon. Why forests vary with respect to their vulnerability to wind events and how the predicted increase in storm events might affect forest ecosystems within the Amazon are not well understood. We found that windthrows are common in the Amazon region extending from northwest (Peru, Colombia, Venezuela, and west Brazil) to central Brazil, with the highest occurrence of windthrows in the northwest Amazon. More frequent winds, produced by more frequent severe convective systems, in combination with well-known processes that limit the anchoring of trees in the soil, help to explain the higher vulnerability of the northwest Amazon forests to winds. Projected increases in the frequency and intensity of convective storms in the Amazon have the potential to increase wind-related tree mortality. A forest demographic model calibrated for the northwestern and the central Amazon showed that northwestern forests are more resilient to increased wind-related tree mortality than forests in the central Amazon. Our study emphasizes the importance of including wind-related tree mortality in model simulations for reliable predictions of the future of tropical forests and their effects on the Earth' system.

Understanding how anthropogenic CO2 emissions will influence future precipitation is critical for sustainably managing ecosystems, particularly for drought-sensitive tropical forests. Although tropical precipitation change remains uncertain, nearly all models from the Coupled Model Intercomparison Project Phase 5 predict a strengthening zonal precipitation asymmetry by 2100, with relative increases over Asian and African tropical forests and decreases over South American forests. Here we show that the plant physiological response to increasing CO2 is a primary mechanism responsible for this pattern. Applying a simulation design in the Community Earth System Model in which CO2 increases are isolated over individual continents, we demonstrate that different circulation, moisture and stability changes arise over each continent due to declines in stomatal conductance and transpiration. The sum of local atmospheric responses over individual continents explains the pan-tropical precipitation asymmetry. Our analysis suggests that South American forests may be more vulnerable to rising CO2 than Asian or African forests.

Thaw and release of permafrost carbon (C) due to climate change is likely to offset increased vegetation C uptake in northern high-latitude (NHL) terrestrial ecosystems. Models project that this permafrost C feedback may act as a slow leak, in which case detection and attribution of the feedback may be difficult. The formation of talik, a subsurface layer of perennially thawed soil, can accelerate permafrost degradation and soil respiration, ultimately shifting the C balance of permafrost-affected ecosystems from long-term C sinks to long-term C sources. It is imperative to understand and characterize mechanistic links between talik, permafrost thaw, and respiration of deep soil C to detect and quantify the permafrost C feedback. Here, we use the Community Land Model (CLM) version 4.5, a permafrost and biogeochemistry model, in comparison to long-term deep borehole data along North American and Siberian transects, to investigate thaw-driven C sources in NHL ( > 55° N) from 2000 to 2300. Widespread talik at depth is projected across most of the NHL permafrost region (14 millionĝ€†km2) by 2300, 6.2 millionĝ€†km2 of which is projected to become a long-term C source, emitting 10 Pgĝ€†C by 2100, 50 Pgĝ€†C by 2200, and 120 Pgĝ€†C by 2300, with few signs of slowing. Roughly half of the projected C source region is in predominantly warm sub-Arctic permafrost following talik onset. This region emits only 20 Pgĝ€†C by 2300, but the CLM4.5 estimate may be biased low by not accounting for deep C in yedoma. Accelerated decomposition of deep soil C following talik onset shifts the ecosystem C balance away from surface dominant processes (photosynthesis and litter respiration), but sink-to-source transition dates are delayed by 20-200 years by high ecosystem productivity, such that talik peaks early ( ĝ1/4 2050s, although borehole data suggest sooner) and C source transition peaks late ( ĝ1/4 2150-2200). The remaining C source region in cold northern Arctic permafrost, which shifts to a net source early (late 21st century), emits 5 times more C (95 Pgĝ€†C) by 2300, and prior to talik formation due to the high decomposition rates of shallow, young C in organic-rich soils coupled with low productivity. Our results provide important clues signaling imminent talik onset and C source transition, including (1) late cold-season (January-February) soil warming at depth ( ĝ1/4 2 m), (2) increasing cold-season emissions (November-April), and (3) enhanced respiration of deep, old C in warm permafrost and young, shallow C in organic-rich cold permafrost soils. Our results suggest a mosaic of processes that govern carbon source-to-sink transitions at high latitudes and emphasize the urgency of monitoring soil thermal profiles, organic C age and content, cold-season CO2 emissions, and atmospheric 14CO2 as key indicators of the permafrost C feedback.

Vegetation phenology plays an important role in regulating land-atmosphere energy, water, and trace-gas exchanges. Changes in spring greenup (SG) have been documented in the past half-century in response to ongoing climate change. We use normalized difference vegetation index generated from NOAA's advanced very high resolution radiometer data in the Global Inventory Modeling and Monitoring Study project over the 1982–2005 period, coupled with climate reanalysis (Climate Research Unit-National Centers for Environmental Prediction) to investigate the SG responses to preseason climate change in northern temperate and boreal regions. We compared these observed responses to the simulated SG responses to preseason climate inferred from the Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) over 1982–2005. The observationally inferred SG suggests that there has been an advance of about 1 days per decade between 1982 and 2005 in the northern midlatitude to high latitude, with significant spatial heterogeneity. The spatial heterogeneity of the SG advance results from heterogeneity in the change of the preseason climate as well as varied vegetation responses to the preseason climate across biomes. The SG to preseason temperature sensitivity is highest in forests other than deciduous needleleaf forests, followed by temperate grasslands and woody savannas. The SG in deciduous needleleaf forests, open shrublands, and tundra is relatively insensitive to preseason temperature. Although the extent of regions where the SG is sensitive to preseason precipitation is smaller than the extent of regions where the SG is sensitive to preseason temperature, the biomes that are more sensitive to temperature are also more sensitive to precipitation, suggesting the interactive control of temperature and precipitation. In the mean, the CMIP5 ESMs reproduced the dominant latitudinal preseason climate trends and SG advances. However, large biases in individual ESMs for the preseason period, climate, and SG sensitivity imply needed model improvements to climate prediction and phenological process parameterizations.

The projected loss of soil carbon to the atmosphere resulting from climate change is a potentially large but highly uncertain feedback to warming. The magnitude of this feedback is poorly constrained by observations and theory, and is disparately represented in Earth system models (ESMs). To assess the climatological temperature sensitivity of soil carbon, we calculate apparent soil carbon turnover times that reflect long-term and broad-scale rates of decomposition. Here, we show that the climatological temperature control on carbon turnover in the top metre of global soils is more sensitive in cold climates than in warm climates and argue that it is critical to capture this emergent ecosystem property in global-scale models. We present a simplified model that explains the observed high cold-climate sensitivity using only the physical scaling of soil freeze-thaw state across climate gradients. Current ESMs fail to capture this pattern, except in an ESM that explicitly resolves vertical gradients in soil climate and carbon turnover. An observed weak tropical temperature sensitivity emerges in a different model that explicitly resolves mineralogical control on decomposition. These results support projections of strong carbon-climate feedbacks from northern soils and demonstrate a method for ESMs to capture this emergent behaviour.

Droughts in the western United States are expected to intensify with climate change. Thus, an adequate representation of ecosystem response to water stress in land models is critical for predicting carbon dynamics. The goal of this study was to evaluate the performance of the Community Land Model (CLM) version 4.5 against observations at an old-growth coniferous forest site in the Pacific Northwest region of the United States (Wind River AmeriFlux site), characterized by a Mediterranean climate that subjects trees to water stress each summer. CLM was driven by site-observed meteorology and calibrated primarily using parameter values observed at the site or at similar stands in the region. Key model adjustments included parameters controlling specific leaf area and stomatal conductance. Default values of these parameters led to significant underestimation of gross primary production, overestimation of evapotranspiration, and consequently overestimation of photosynthetic 13C discrimination, reflected in reduced 13°C 12°C ratios of carbon fluxes and pools. Adjustments in soil hydraulic parameters within CLM were also critical, preventing significant underestimation of soil water content and unrealistic soil moisture stress during summer. After calibration, CLM was able to simulate energy and carbon fluxes, leaf area index, biomass stocks, and carbon isotope ratios of carbon fluxes and pools in reasonable agreement with site observations. Overall, the calibrated CLM was able to simulate the observed response of canopy conductance to atmospheric vapor pressure deficit (VPD) and soil water content, reasonably capturing the impact of water stress on ecosystem functioning. Both simulations and observations indicate that stomatal response from water stress at Wind River was primarily driven by VPD and not soil moisture. The calibration of the Ball-Berry stomatal conductance slope (mbb) at Wind River aligned with findings from recent CLM experiments at sites characterized by the same plant functional type (needleleaf evergreen temperate forest), despite significant differences in stand composition and age and climatology, suggesting that CLM could benefit from a revised mbb value of 6, rather than the default value of 9, for this plant functional type. Conversely, Wind River required a unique calibration of the hydrology submodel to simulate soil moisture, suggesting that the default hydrology has a more limited applicability. This study demonstrates that carbon isotope data can be used to constrain stomatal conductance and intrinsic water use efficiency in CLM, as an alternative to eddy covariance flux measurements. It also demonstrates that carbon isotopes can expose structural weaknesses in the model and provide a key constraint that may guide future model development.

Increasing atmospheric methane (CH4) concentrations have contributed to approximately 20% of anthropogenic climate change. Despite the importance of CH4 as a greenhouse gas, its atmospheric growth rate and dynamics over the past two decades, which include a stabilization period (1999-2006), followed by renewed growth starting in 2007, remain poorly understood. We provide an updated estimate of CH4 emissions from wetlands, the largest natural global CH4 source, for 2000-2012 using an ensemble of biogeochemical models constrained with remote sensing surface inundation and inventory-based wetland area data. Between 2000-2012, boreal wetland CH4 emissions increased by 1.2 Tg yr-1 (-0.2-3.5 Tg yr-1), tropical emissions decreased by 0.9 Tg yr-1 (-3.2-1.1 Tg yr-1), yet globally, emissions remained unchanged at 184 22 Tg yr-1. Changing air temperature was responsible for increasing high-latitude emissions whereas declines in low-latitude wetland area decreased tropical emissions; both dynamics are consistent with features of predicted centennial-scale climate change impacts on wetland CH4 emissions. Despite uncertainties in wetland area mapping, our study shows that global wetland CH4 emissions have not contributed significantly to the period of renewed atmospheric CH4 growth, and is consistent with findings from studies that indicate some combination of increasing fossil fuel and agriculture-related CH4 emissions, and a decrease in the atmospheric oxidative sink.

Evapotranspiration (ET) is a critical term in the surface energy budget as well as the water cycle. There are few direct measurements of ET, and thus the magnitude and variability are poorly constrained at large spatial scales. Estimates of the annual cycle of ET over the Amazon are critical because they influence predictions of the seasonal cycle of carbon fluxes, as well as atmospheric dynamics and circulation. The authors estimate ET for the Amazon basin using a water budget approach by differencing rainfall, discharge, and time-varying storage from the Gravity Recovery and Climate Experiment. It is found that the climatological annual cycle of ET over the Amazon basin upstream of óbidos shows suppression of ET during the wet season and higher ET during the dry season, consistent with flux-tower-based observations in seasonally dry forests. They also find a statistically significant decrease in ET over the time period 2002-15 of -1.46 mm yr-1. The direct estimate of the seasonal cycle of ET is largely consistent with previous indirect estimates, including energy-budget-based approaches, an upscaled station-based estimate, and land surface model estimates, but suggests that suppression of ET during the wet season is underestimated by existing products.

Measurements of the stable carbon isotope ratio (δ13C) on annual tree rings offer new opportunities to evaluate mechanisms of variations in photosynthesis and stomatal conductance under changing CO2 and climate conditions, especially in conjunction with process-based biogeochemical model simulations. The isotopic discrimination is indicative of the ratio between the CO2 partial pressure in the intercellular cavities and the atmosphere (ci=ca) and of the ratio of assimilation to stomatal conductance, termed intrinsic wateruse efficiency (iWUE). We performed isotope-enabled simulations over the industrial period with the land biosphere module (CLM4.5) of the Community Earth System Model and the Land Surface Processes and Exchanges (LPX-Bern) dynamic global vegetation model. Results for C3 tree species show good agreement with a global compilation of δ13C measurements on leaves, though modeled 13C discrimination by C3 trees is smaller in arid regions than measured. A compilation of 76 tree-ring records, mainly from Europe, boreal Asia, and western North America, suggests on average small 20th century changes in isotopic discrimination and in ci=ca and an increase in iWUE of about 27% since 1900. LPXBern results match these century-scale reconstructions, supporting the idea that the physiology of stomata has evolved to optimize trade-offs between carbon gain by assimilation and water loss by transpiration. In contrast, CLM4.5 simulates an increase in discrimination and in turn a change in iWUE that is almost twice as large as that revealed by the tree-ring data. Factorial simulations show that these changes are mainly in response to rising atmospheric CO2. The results suggest that the downregulation of ci=ca and of photosynthesis by nitrogen limitation is possibly too strong in the standard setup of CLM4.5 or that there may be problems associated with the implementation of conductance, assimilation, and related adjustment processes on long-term environmental changes.

Land surface models are useful tools to quantify contemporary and future climate impact on terrestrial carbon cycle processes, provided they can be appropriately constrained and tested with observations. Stable carbon isotopes of CO2 offer the potential to improve model representation of the coupled carbon and water cycles because they are strongly influenced by stomatal function. Recently, a representation of stable carbon isotope discrimination was incorporated into the Community Land Model component of the Community Earth System Model. Here, we tested the model's capability to simulate whole-forest isotope discrimination in a subalpine conifer forest at Niwot Ridge, Colorado, USA. We distinguished between isotopic behavior in response to a decrease of δ13C within atmospheric CO2 (Suess effect) vs. photosynthetic discrimination (Δcanopy), by creating a site-customized atmospheric CO2 and δ13C of CO2 time series. We implemented a seasonally varying Vcmax model calibration that best matched site observations of net CO2 carbon exchange, latent heat exchange, and biomass. The model accurately simulated observed δ13C of needle and stem tissue, but underestimated the δ13C of bulk soil carbon by 1-2%. The model overestimated the multiyear (2006-2012) average Δcanopy relative to prior data-based estimates by 2-4%. The amplitude of the average seasonal cycle of Δcanopy (i.e., higher in spring/fall as compared to summer) was correctly modeled but only when using a revised, fully coupled An - gs (net assimilation rate, stomatal conductance) version of the model in contrast to the partially coupled An - gs version used in the default model. The model attributed most of the seasonal variation in discrimination to An, whereas interannual variation in simulated Δcanopy during the summer months was driven by stomatal response to vapor pressure deficit (VPD). The model simulated a 10% increase in both photosynthetic discrimination and water-use efficiency (WUE) since 1850 which is counter to established relationships between discrimination and WUE. The isotope observations used here to constrain CLM suggest (1) the model overestimated stomatal conductance and (2) the default CLM approach to representing nitrogen limitation (partially coupled model) was not capable of reproducing observed trends in discrimination. These findings demonstrate that isotope observations can provide important information related to stomatal function driven by environmental stress from VPD and nitrogen limitation. Future versions of CLM that incorporate carbon isotope discrimination are likely to benefit from explicit inclusion of mesophyll conductance.

Wetlands are the largest global natural methane (CH4/ source, and emissions between 50 and 70° N latitude contribute 10-30% to this source. Predictive capability of land models for northern wetland CH4 emissions is still low due to limited site measurements, strong spatial and temporal variability in emissions, and complex hydrological and biogeochemical dynamics. To explore this issue, we compare wetland CH4 emission predictions from the Community Land Model 4.5 (CLM4.5-BGC) with siteto regional-scale observations. A comparison of the CH4 fluxes with eddy flux data highlighted needed changes to the model's estimate of aerenchyma area, which we implemented and tested. The model modification substantially reduced biases in CH4 emissions when compared with CarbonTracker CH4 predictions. CLM4.5 CH4 emission predictions agree well with growing season (May-September) CarbonTracker Alaskan regional-level CH4 predictions and sitelevel observations. However, CLM4.5 underestimated CH4 emissions in the cold season (October-April). The monthly atmospheric CH4 mole fraction enhancements due to wetland emissions are also assessed using the Weather Research and Forecasting-Stochastic Time-Inverted Lagrangian Transport (WRF-STILT) model coupled with daily emissions from CLM4.5 and compared with aircraft CH4 mole fraction measurements from the Carbon in Arctic Reservoirs Vulnerability Experiment (CARVE) campaign. Both the tower and aircraft analyses confirm the underestimate of cold-season CH4 emissions by CLM4.5. The greatest uncertainties in predicting the seasonal CH4 cycle are from the wetland extent, coldseason CH4 production and CH4 transport processes. We recommend more cold-season experimental studies in highlatitude systems, which could improve the understanding and parameterization of ecosystem structure and function during this period. Predicted CH4 emissions remain uncertain, but we show here that benchmarking against observations across spatial scales can inform model structural and parameter improvements.

In many ecosystems, nitrogen is the most limiting nutrient for plant growth and productivity. However, current Earth System Models (ESMs) do not mechanistically represent functional nitrogen allocation for photosynthesis or the linkage between nitrogen uptake and root traits. The current version of CLM (4.5) links nitrogen availability and plant productivity via (1) an instantaneous downregulation of potential photosynthesis rates based on soil mineral nitrogen availability, and (2) apportionment of soil nitrogen between plants and competing nitrogen consumers assumed to be proportional to their relative N demands. However, plants do not photosynthesize at potential rates and then downregulate; instead photosynthesis rates are governed by nitrogen that has been allocated to the physiological processes underpinning photosynthesis. Furthermore, the role of plant roots in nutrient acquisition has also been largely ignored in ESMs. We therefore present a new plant nitrogen model for CLM4.5 with (1) improved representations of linkages between leaf nitrogen and plant productivity based on observed relationships in a global plant trait database and (2) plant nitrogen uptake based on root-scale Michaelis-Menten uptake kinetics. Our model improvements led to a global bias reduction in GPP, LAI, and biomass of 70%, 11%, and 49%, respectively. Furthermore, water use efficiency predictions were improved conceptually, qualitatively, and in magnitude. The new model's GPP responses to nitrogen deposition, CO2 fertilization, and climate also differed from the baseline model. The mechanistic representation of leaf-level nitrogen allocation and a theoretically consistent treatment of competition with belowground consumers led to overall improvements in global carbon cycling predictions.

Soil is a complex system where biotic (e.g., plant roots, micro-organisms) and abiotic (e.g., mineral surfaces) consumers compete for resources necessary for life (e.g., nitrogen, phosphorus). This competition is ecologically significant, since it regulates the dynamics of soil nutrients and controls aboveground plant productivity. Here we develop, calibrate and test a nutrient competition model that accounts for multiple soil nutrients interacting with multiple biotic and abiotic consumers. As applied here for tropical forests, the Nutrient COMpetition model (N-COM) includes three primary soil nutrients (NHC+4 , NO-3 and POx representing the sum of PO3-4 , HPO2-4 and H2PO-4/and five potential competitors (plant roots, decomposing microbes, nitrifiers, denitrifiers and mineral surfaces). The competition is formulated with a quasi-steady-state chemical equilibrium approximation to account for substrate (multiple substrates share one consumer) and consumer (multiple consumers compete for one substrate) effects. N-COM successfully reproduced observed soil heterotrophic respiration, N2O emissions, free phosphorus, sorbed phosphorus and NHC 4 pools at a tropical forest site (Tapajos). The overall model uncertainty was moderately well constrained. Our sensitivity analysis revealed that soil nutrient competition was primarily regulated by consumer-substrate affinity rather than environmental factors such as soil temperature or soil moisture. Our results also imply that under strong nutrient limitation, relative competitiveness depends strongly on the competitor functional traits (affinity and nutrient carrier enzyme abundance).We then applied the N-COM model to analyze field nitrogen and phosphorus perturbation experiments in two tropical forest sites (in Hawaii and Puerto Rico) not used in model development or calibration. Under soil inorganic nitrogen and phosphorus elevated conditions, the model accurately replicated the experimentally observed competition among nutrient consumers. Although we used as many observations as we could obtain, more nutrient addition experiments in tropical systems would greatly benefit model testing and calibration. In summary, the N-COM model provides an ecologically consistent representation of nutrient competition appropriate for land BGC models integrated in Earth System Models.

We present an approach to estimate the feedback from large-scale thawing of permafrost soils using a simplified, data-constrained model that combines three elements: soil carbon (C) maps and profiles to identify the distribution and type of C in permafrost soils; incubation experiments to quantify the rates of C lost after thaw; and models of soil thermal dynamics in response to climate warming. We call the approach the Permafrost Carbon Network Incubation-Panarctic Thermal scaling approach (PInc-PanTher). The approach assumes that C stocks do not decompose at all when frozen, but once thawed follow set decomposition trajectories as a function of soil temperature. The trajectories are determined according to a three-pool decomposition model fitted to incubation data using parameters specific to soil horizon types. We calculate litterfall C inputs required to maintain steady-state C balance for the current climate, and hold those inputs constant. Soil temperatures are taken from the soil thermal modules of ecosystem model simulations forced by a common set of future climate change anomalies under two warming scenarios over the period 2010 to 2100. Under a medium warming scenario (RCP4.5), the approach projects permafrost soil C losses of 12.2-33.4 Pg C; under a high warming scenario (RCP8.5), the approach projects C losses of 27.9-112.6 Pg C. Projected C losses are roughly linearly proportional to global temperature changes across the two scenarios. These results indicate a global sensitivity of frozen soil C to climate change (γ sensitivity) of -14 to -19 Pg C °C(-1) on a 100 year time scale. For CH4 emissions, our approach assumes a fixed saturated area and that increases in CH4 emissions are related to increased heterotrophic respiration in anoxic soil, yielding CH4 emission increases of 7% and 35% for the RCP4.5 and RCP8.5 scenarios, respectively, which add an additional greenhouse gas forcing of approximately 10-18%. The simplified approach presented here neglects many important processes that may amplify or mitigate C release from permafrost soils, but serves as a data-constrained estimate on the forced, large-scale permafrost C response to warming.

We present an approach to estimate the feedback from large-scale thawing of permafrost soils using a simplified, data-constrained model that combines three elements: soil carbon (C) maps and profiles to identify the distribution and type of C in permafrost soils; incubation experiments to quantify the rates of C lost after thaw; and models of soil thermal dynamics in response to climate warming. We call the approach the Permafrost Carbon Network Incubation-Panarctic Thermal scaling approach (PInc-PanTher). The approach assumes that C stocks do not decompose at all when frozen, but once thawed follow set decomposition trajectories as a function of soil temperature. The trajectories are determined according to a three-pool decomposition model fitted to incubation data using parameters specific to soil horizon types. We calculate litterfall C inputs required to maintain steady-state C balance for the current climate, and hold those inputs constant. Soil temperatures are taken from the soil thermal modules of ecosystem model simulations forced by a common set of future climate change anomalies under two warming scenarios over the period 2010 to 2100. Under a medium warming scenario (RCP4.5), the approach projects permafrost soil C losses of 12.2-33.4 Pg C; under a high warming scenario (RCP8.5), the approach projects C losses of 27.9-112.6 Pg C. Projected C losses are roughly linearly proportional to global temperature changes across the two scenarios. These results indicate a global sensitivity of frozen soil C to climate change (γ sensitivity) of -14 to -19 Pg C °C(-1) on a 100 year time scale. For CH4 emissions, our approach assumes a fixed saturated area and that increases in CH4 emissions are related to increased heterotrophic respiration in anoxic soil, yielding CH4 emission increases of 7% and 35% for the RCP4.5 and RCP8.5 scenarios, respectively, which add an additional greenhouse gas forcing of approximately 10-18%. The simplified approach presented here neglects many important processes that may amplify or mitigate C release from permafrost soils, but serves as a data-constrained estimate on the forced, large-scale permafrost C response to warming.

To better understand sources of uncertainty in projections of terrestrial carbon cycle feedbacks, we present an approach to separate the controls on modeled carbon changes. We separate carbon changes into four categories using a linearized, equilibrium approach: those arising from changed inputs (productivity-driven changes), and outputs (turnover-driven changes), of both the live and dead carbon pools. Using Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations for five models, we find that changes to the live pools are primarily explained by productivity-driven changes, with only one model showing large compensating changes to live carbon turnover times. For dead carbon pools, the situation is more complex as all models predict a large reduction in turnover times in response to increases in productivity. This response arises from the common representation of a broad spectrum of decomposition turnover times via a multi-pool approach, in which flux-weighted turnover times are faster than mass-weighted turnover times. This leads to a shift in the distribution of carbon among dead pools in response to changes in inputs, and therefore a transient but long-lived reduction in turnover times. Since this behavior, a reduction in inferred turnover times resulting from an increase in inputs, is superficially similar to priming processes, but occurring without the mechanisms responsible for priming, we call the phenomenon "false priming", and show that it masks much of the intrinsic changes to dead carbon turnover times as a result of changing climate. These patterns hold across the fully coupled, biogeochemically coupled, and radiatively coupled 1 % yr-1 increasing CO2 experiments. We disaggregate inter-model uncertainty in the globally integrated equilibrium carbon responses to initial turnover times, initial productivity, fractional changes in turnover, and fractional changes in productivity. For both the live and dead carbon pools, inter-model spread in carbon changes arising from initial conditions is dominated by model disagreement on turnover times, whereas inter-model spread in carbon changes from fractional changes to these terms is dominated by model disagreement on changes to productivity in response to both warming and CO2 fertilization. However, the lack of changing turnover time control on carbon responses, for both live and dead carbon pools, in response to the imposed forcings may arise from a common lack of process representation behind changing turnover times (e.g., allocation and mortality for live carbon; permafrost, microbial dynamics, and mineral stabilization for dead carbon), rather than a true estimate of the importance of these processes.

A significant fraction of anthropogenic CO2 emissions is assimilated by tropical forests and stored as biomass, slowing the accumulation of CO2 in the atmosphere. Because different plant tissues have different functional roles and turnover times, predictions of carbon balance of tropical forests depend on how earth system models (ESMs) represent the dynamic allocation of productivity to different tree compartments. This study shows that observed allocation of productivity, biomass, and turnover times of main tree compartments (leaves, wood, and roots) are not accurately represented in Coupled Model Intercomparison Project Phase 5 ESMs. In particular, observations indicate that biomass saturates with increasing productivity. In contrast, most models predict continuous increases in biomass with increases in productivity. This bias may lead to an over-prediction of carbon uptake in response to CO2 or climate-driven changes in productivity. Compartment-specific productivity and biomass are useful benchmarks to assess terrestrial ecosystem model performance. Improvements in the predicted allocation patterns and turnover times by ESMs will reduce uncertainties in climate predictions.

Recent efforts to narrow the spread in equilibrium climate sensitivity (ECS) across global climate models have focused on identifying observationally based constraints, which are rooted in empirical correlations between ECS and biases in the models' present-day climate. This study reexamines one such constraint identified from CMIP3 models: the linkage between ECS and net top-of-the-atmosphere radiation biases in the Southern Hemisphere (SH). As previously documented, the intermodel spread in the ECS of CMIP3 models is linked to present-day cloud and net radiation biases over the midlatitude Southern Ocean, where higher cloud fraction in the present-day climate is associated with larger values of ECS. However, in this study, no physical explanation is found to support this relationship. Furthermore, it is shown here that this relationship disappears in CMIP5 models and is unique to a subset of CMIP models characterized by unrealistically bright present-day clouds in the SH subtropics. In view of this evidence, Southern Ocean cloud and net radiation biases appear inappropriate for providing observationally based constraints on ECS. Instead of the Southern Ocean, this study points to the stratocumulus-to-cumulus transition regions of the SH subtropical oceans as key to explaining the intermodel spread in the ECS of both CMIP3 and CMIP5 models. In these regions, ECS is linked to present-day cloud and net radiation biases with a plausible physical mechanism: models with brighter subtropical clouds in the present-day climate show greater ECS because 1) subtropical clouds dissipate with increasing CO2 concentrations in many models and 2) the dissipation of brighter clouds contributes to greater solar warming of the surface.

Century-long observed gridded land precipitation datasets are a cornerstone of hydrometeorological research. But recent work has suggested that observed Northern Hemisphere midlatitude (NHML) land mean precipitation does not show evidence of an expected negative response to mid-twentieth-century aerosol forcing. Utilizing observed river discharges, the observed runoff is calculated and compared with observed land precipitation. The results show a near-zero twentieth-century trend in observed NHML land mean runoff, in contrast to the significant positive trend in observed NHML land mean precipitation. However, precipitation and runoff share common interannual and decadal variability. An obvious split, or breakpoint, is found in the NHML land mean runoff-precipitation relationship in the 1930s. Using runoff simulated by six land surface models (LSMs), which are driven by the observed precipitation dataset, such breakpoints are absent. These findings support previous hypotheses that inhomogeneities exist in the early-twentieth-century NHML land mean precipitation record. Adjusting the observed precipitation record according to the observed runoff record largely accounts for the departure of the observed precipitation response from that predicted given the real-world aerosol forcing estimate, more than halving the discrepancy from about 6 to around 2 W m-2. Consideration of complementary observed runoff adds support to the suggestion that NHML-wide early-twentieth-century precipitation observations are unsuitable for climate change studies. The agreement between precipitation and runoff over Europe, however, is excellent, supporting the use of whole-twentieth-century observed precipitation datasets here.

Approximately 1700 Pg of soil carbon (C) are stored in the northern circumpolar permafrost zone, more than twice as much C than in the atmosphere. The overall amount, rate, and form of C released to the atmosphere in a warmer world will influence the strength of the permafrost C feedback to climate change. We used a survey to quantify variability in the perception of the vulnerability of permafrost C to climate change. Experts were asked to provide quantitative estimates of permafrost change in response to four scenarios of warming. For the highest warming scenario (RCP 8.5), experts hypothesized that C release from permafrost zone soils could be 19-45 Pg C by 2040, 162-288 Pg C by 2100, and 381-616 Pg C by 2300 in CO2 equivalent using 100-year CH4 global warming potential (GWP). These values become 50 % larger using 20-year CH4 GWP, with a third to a half of expected climate forcing coming from CH4 even though CH4 was only 2.3 % of the expected C release. Experts projected that two-thirds of this release could be avoided under the lowest warming scenario (RCP 2.6). These results highlight the potential risk from permafrost thaw and serve to frame a hypothesis about the magnitude of this feedback to climate change. However, the level of emissions proposed here are unlikely to overshadow the impact of fossil fuel burning, which will continue to be the main source of C emissions and climate forcing. © 2013 The Author(s).