Explore the vast range of titles available.

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. Plain Language Summary Land surface models (LSMs) are the part of climate models that simulate processes happening at the Earth's surface. These include reflection of the sunlight, evaporation from ecosystems, and the amount of carbon from human emissions that the land takes up. LSMs also need to simulate how human management of the land surface changes the climate both directly (e.g., via the effect on evaporation) and in the long term (via changing the amount of carbon stored in wood and soil). Not surprisingly, trying to make a single mathematical representation of all of these different parts of the Earth system is difficult. Here we discuss themes that repeatedly affect all teams developing LSMs: how to manage the increasing number of complicated model components, how to represent the high degree of variability of the land surface, and how to predict how the properties of the surface (particularly those of plant communities) will change. These are large problems, with no obvious easy solutions. We hope to spark discussion and investment into their resolution, concomitant with the increasing importance of LSMs as our best tools for translating possible trajectories of climate change into impacts on humans, ecosystems, food and water supplies, and river systems. Key Points Land surface models have grown in complexity, and new methods of managing this complexity are required for scientific understanding New methods are also needed to represent, classify, and benchmark models across the multidimensional heterogeneity of the land surface A further challenge is to constrain model parameters in ways that are consistent with allowing long‐term ecological dynamics to occur

Permafrost soils contain enormous amounts of organic carbon whose stability is contingent on remaining frozen. With future warming, these soils may release carbon to the atmosphere and act as a positive feedback to climate change. Significant uncertainty remains on the postthaw carbon dynamics of permafrost-affected ecosystems, in particular since most of the carbon resides at depth where decomposition dynamics may differ from surface soils, and since nitrogen mineralized by decomposition may enhance plant growth. Here we show, using a carbon–nitrogen model that includes permafrost processes forced in an unmitigated warming scenario, that the future carbon balance of the permafrost region is highly sensitive to the decomposability of deeper carbon, with the net balance ranging from 21 Pg C to 164 Pg C losses by 2300. Increased soil nitrogen mineralization reduces nutrient limitations, but the impact of deep nitrogen on the carbon budget is small due to enhanced nitrogen availability from warming surface soils and seasonal asynchrony between deeper nitrogen availability and plant nitrogen demands. Although nitrogen dynamics are highly uncertain, the future carbon balance of this region is projected to hinge more on the rate and extent of permafrost thaw and soil decomposition than on enhanced nitrogen availability for vegetation growth resulting from permafrost thaw.

Soil carbon release remains a highly uncertain climate feedback. Research now shows that the temperature control on carbon turnover is more sensitive in cold climates, supporting projections of a strong carbon–climate feedback from northern soils. 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) 1 , 2 , 3 . To assess the climatological temperature sensitivity of soil carbon, we calculate apparent soil carbon turnover times 4 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 5 , 6 and demonstrate a method for ESMs to capture this emergent behaviour.

The authors analyze global climate model predictions of soil temperature [from the Coupled Model Intercomparison Project phase 5 (CMIP5) database] to assess the models’ representation of current-climate soil thermal dynamics and their predictions of permafrost thaw during the twenty-first century. The authors compare the models’ predictions with observations of active layer thickness, air temperature, and soil temperature and with theoretically expected relationships between active layer thickness and air temperature annual mean- and seasonal-cycle amplitude. Models show a wide range of current permafrost areas, active layer statistics (cumulative distributions, correlations with mean annual air temperature, and amplitude of seasonal air temperature cycle), and ability to accurately model the coupling between soil and air temperatures at high latitudes. Many of the between-model differences can be traced to differences in the coupling between either near-surface air and shallow soil temperatures or shallow and deeper (1 m) soil temperatures, which in turn reflect differences in snow physics and soil hydrology. The models are compared with observational datasets to benchmark several aspects of the permafrost-relevant physics of the models. The CMIP5 models following multiple representative concentration pathways (RCP) show a wide range of predictions for permafrost loss: 2%–66% for RCP2.6, 15%–87% for RCP4.5, and 30%–99% for RCP8.5. Normalizing the amount of permafrost loss by the amount of high-latitude warming in the RCP4.5 scenario, the models predict an absolute loss of 1.6 ± 0.7 million km² permafrost per 1°C high-latitude warming, or a fractional loss of 6%–29% °C−1.

Permafrost soils contain enormous amounts of organic carbon, which could act as a positive feedback to global climate change due to enhanced respiration rates with warming. We have used a terrestrial ecosystem model that includes permafrost carbon dynamics, inhibition of respiration in frozen soil layers, vertical mixing of soil carbon from surface to permafrost layers, and CH4 emissions from flooded areas, and which better matches new circumpolar inventories of soil carbon stocks, to explore the potential for carbon-climate feedbacks at high latitudes. Contrary to model results for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4), when permafrost processes are included, terrestrial ecosystems north of 60°N could shift from being a sink to a source of CO2 by the end of the 21st century when forced by a Special Report on Emissions Scenarios (SRES) A2 climate change scenario. Between 1860 and 2100, the model response to combined CO2 fertilization and climate change changes from a sink of 68 Pg to a 27 + -7 Pg sink to 4 + -18 Pg source, depending on the processes and parameter values used. The integrated change in carbon due to climate change shifts from near zero, which is within the range of previous model estimates, to a climate-induced loss of carbon by ecosystems in the range of 25 + -3 to 85 + -16 Pg C, depending on processes included in the model, with a best estimate of a 62 + -7 Pg C loss. Methane emissions from high-latitude regions are calculated to increase from 34 Tg CH4/y to 41–70 Tg CH4/y, with increases due to CO2 fertilization, permafrost thaw, and warming-induced increased CH4 flux densities partially offset by a reduction in wetland extent.

Rising atmospheric CO₂ will make Earth warmer, and many studies have inferred that this warming will cause droughts to become more widespread and severe. However, rising atmospheric CO₂ also modifies stomatal conductance and plant water use, processes that are often are overlooked in impact analysis. We find that plant physiological responses to CO₂ reduce predictions of future drought stress, and that this reduction is captured by using plant-centric rather than atmosphere-centric metrics from Earth system models (ESMs). The atmosphere-centric Palmer Drought Severity Index predicts future increases in drought stress for more than 70% of global land area. This area drops to 37% with the use of precipitation minus evapotranspiration (P-E), a measure that represents the water flux available to downstream ecosystems and humans. The two metrics yield consistent estimates of increasing stress in regions where precipitation decreases are more robust (southern North America, northeastern South America, and southern Europe). The metrics produce diverging estimates elsewhere, with P-E predicting decreasing stress across temperate Asia and central Africa. The differing sensitivity of drought metrics to radiative and physiological aspects of increasing CO₂ partly explains the divergent estimates of future drought reported in recent studies. Further, use of ESM output in offline models may double-count plant feedbacks on relative humidity and other surface variables, leading to overestimates of future stress. The use of drought metrics that account for the response of plant transpiration to changing CO₂, including direct use of P-E and soil moisture from ESMs, is needed to reduce uncertainties in future assessment.

Earth’s high latitudes are projected to experience warmer and wetter summers in the future but ramifications for soil thermal processes and permafrost thaw are poorly understood. Here we present 2750 end of summer thaw depths representing a range of vegetation characteristics in Interior Alaska measured over a 5 year period. This included the top and third wettest summers in the 91-year record and three summers with precipitation close to mean historical values. Increased rainfall led to deeper thaw across all sites with an increase of 0.7 ± 0.1 cm of thaw per cm of additional rain. Disturbed and wetland sites were the most vulnerable to rain-induced thaw with ~1 cm of surface thaw per additional 1 cm of rain. Permafrost in tussock tundra, mixed forest, and conifer forest was less sensitive to rain-induced thaw. A simple energy budget model yields seasonal thaw values smaller than the linear regression of our measurements but provides a first-order estimate of the role of rain-driven sensible heat fluxes in high-latitude terrestrial permafrost. This study demonstrates substantial permafrost thaw from the projected increasing summer precipitation across most of the Arctic region.

Carbon cycle feedbacks represent large uncertainties in climate change projections, and the response of soil carbon to climate change contributes the greatest uncertainty to this. Future changes in soil carbon depend on changes in litter and root inputs from plants and especially on reductions in the turnover time of soil carbon ( τ s ) with warming. An approximation to the latter term for the top one metre of soil (Δ C s,τ ) can be diagnosed from projections made with the CMIP6 and CMIP5 Earth System Models (ESMs), and is found to span a large range even at 2 °C of global warming (−196 ± 117 PgC). Here, we present a constraint on Δ C s,τ , which makes use of current heterotrophic respiration and the spatial variability of τ s inferred from observations. This spatial emergent constraint allows us to halve the uncertainty in Δ C s,τ at 2 °C to −232 ± 52 PgC. The fate of the carbon locked away in soil is uncertain, and there are vast differences between models. Here the authors apply observational, spatio-temporal constraints on carbon turnover projections and find that uncertainty in estimations of carbon dynamics are reduced by 50%.

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. Key Point The ILAMB benchmarking system broadly compares models to observational data sets and provides a synthesis of overall performance

The permafrost zone is expected to be a substantial carbon source to the atmosphere, yet large-scale models currently only simulate gradual changes in seasonally thawed soil. Abrupt thaw will probably occur in <20% of the permafrost zone but could affect half of permafrost carbon through collapsing ground, rapid erosion and landslides. Here, we synthesize the best available information and develop inventory models to simulate abrupt thaw impacts on permafrost carbon balance. Emissions across 2.5 million km2 of abrupt thaw could provide a similar climate feedback as gradual thaw emissions from the entire 18 million km2 permafrost region under the warming projection of Representative Concentration Pathway 8.5. While models forecast that gradual thaw may lead to net ecosystem carbon uptake under projections of Representative Concentration Pathway 4.5, abrupt thaw emissions are likely to offset this potential carbon sink. Active hillslope erosional features will occupy 3% of abrupt thaw terrain by 2300 but emit one-third of abrupt thaw carbon losses. Thaw lakes and wetlands are methane hot spots but their carbon release is partially offset by slowly regrowing vegetation. After considering abrupt thaw stabilization, lake drainage and soil carbon uptake by vegetation regrowth, we conclude that models considering only gradual permafrost thaw are substantially underestimating carbon emissions from thawing permafrost.Analyses of inventory models under two climate change projection scenarios suggest that carbon emissions from abrupt thaw of permafrost through ground collapse, erosion and landslides could contribute significantly to the overall permafrost carbon balance.