Explore the vast range of titles available.

The nitrogen cycle and its effect on carbon uptake in the terrestrial biosphere is a recent progression in earth system models. As with any new component of a model, it is important to understand the behaviour, strengths, and limitations of the various process representations. Here we assess and compare five land surface models with nitrogen cycles that are used as the terrestrial components of some of the earth system models in CMIP6. The land surface models were run offline with a common spin-up and forcing protocol. We use a historical control simulation and two perturbations to assess the model nitrogen-related performances: a simulation with atmospheric carbon dioxide increased by 200 ppm and one with nitrogen deposition increased by 50 kgN ha−1 yr−1. There is generally greater variability in productivity response between models to increased nitrogen than to carbon dioxide. Across the five models the response to carbon dioxide globally was 5 % to 20 % and the response to nitrogen was 2 % to 24 %. The models are not evenly distributed within the ensemble range, with two of the models having low productivity response to nitrogen and another one with low response to elevated atmospheric carbon dioxide, compared to the other models. In all five models individual grid cells tend to exhibit bimodality, with either a strong response to increased nitrogen or atmospheric carbon dioxide but rarely to both to an equal extent. However, this local effect does not scale to either the regional or global level. The global and tropical responses are generally more accurately modelled than boreal, tundra, or other high-latitude areas compared to observations. These results are due to divergent choices in the representation of key nitrogen cycle processes. They show the need for more observational studies to enhance understanding of nitrogen cycle processes, especially nitrogen-use efficiency and biological nitrogen fixation.

The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W/sq. m, comprised of 1.81 (±0.09) W/sq. m from CO2, 1.08 (± 0.21) W/sq. m from other well-mixed greenhouse gases, −1.01 (± 0.23) W/sq. m from aerosols and −0.09 (±0.13) W/sq. m from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W/sq. m. The majority of the remaining 0.21 W/sq. m is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W/sq. m, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

We document the development of the first version of the U.K. Earth System Model UKESM1. The model represents a major advance on its predecessor HadGEM2‐ES, with enhancements to all component models and new feedback mechanisms. These include a new core physical model with a well‐resolved stratosphere; terrestrial biogeochemistry with coupled carbon and nitrogen cycles and enhanced land management; tropospheric‐stratospheric chemistry allowing the holistic simulation of radiative forcing from ozone, methane, and nitrous oxide; two‐moment, five‐species, modal aerosol; and ocean biogeochemistry with two‐way coupling to the carbon cycle and atmospheric aerosols. The complexity of coupling between the ocean, land, and atmosphere physical climate and biogeochemical cycles in UKESM1 is unprecedented for an Earth system model. We describe in detail the process by which the coupled model was developed and tuned to achieve acceptable performance in key physical and Earth system quantities and discuss the challenges involved in mitigating biases in a model with complex connections between its components. Overall, the model performs well, with a stable pre‐industrial state and good agreement with observations in the latter period of its historical simulations. However, global mean surface temperature exhibits stronger‐than‐observed cooling from 1950 to 1970, followed by rapid warming from 1980 to 2014. Metrics from idealized simulations show a high climate sensitivity relative to previous generations of models: Equilibrium climate sensitivity is 5.4 K, transient climate response ranges from 2.68 to 2.85 K, and transient climate response to cumulative emissions is 2.49 to 2.66 K TtC−1. Plain Language Summary We describe the development and behavior of UKESM1, a novel climate model that includes improved representations of processes in the atmosphere, ocean, and on land. These processes are inter‐related: For example, dust is produced on the land and blown up into the atmosphere where it affects the amount of sunlight falling on Earth. Dust can also be dissolved in the ocean, where it affects marine life. This in turn changes both the amount of carbon dioxide absorbed by the ocean and the material emitted from the surface into the atmosphere, which has an affect on the formation of clouds. UKESM1 includes many processes and interactions such as these, giving it a high level of complexity. Ensuring realistic process behavior is a major challenge in the development of our model, and we have carefully tested this. UKESM1 performs well, correctly exhibiting stable results from a continuous pre‐industrial simulation (used to provide a reference for future experiments) and showing good agreement with observations toward the end of its historical simulations. Results for some properties—including the degree to which average surface temperature changes with increased amounts of carbon dioxide in the atmosphere—are examined in detail. Key Points UKESM1 represents a major advance over its predecessor HadGEM2‐ES, both in the complexity of its components and its internal coupling The complex coupling presents challenges to the model development; we document the tuning process employed to obtain acceptable performance UKESM1 performs well, having a stable pre‐industrial state and showing good agreement with observations in a wide variety of contexts

Land-use change affects both the quality and quantity of soil organic carbon (SOC) and leads to changes in ecosystem functions such as productivity and environmental regulation. Future changes in SOC are, however, highly uncertain owing to its heterogeneity and complexity. In this study, we analyzed the outputs of simulations of SOC stock by Earth system models (ESMs), most of which are participants in the Land-Use Model Intercomparison Project. Using a common protocol and the same forcing data, the ESMs simulated SOC distribution patterns and their changes during historical (1850-2014) and future (2015-2100) periods. Total SOC stock increased in many simulations over the historical period (30 ± 67 Pg C) and under future climate and land-use conditions (48 ± 32 Pg C for ssp126 and 49 ± 58 Pg C for ssp370). Land-use experiments indicated that changes in SOC attributable to land-use scenarios were modest at the global scale, in comparison with climatic and rising CO2 impacts, but they were notable in several regions. Future net soil carbon sequestration rates estimated by the ESMs were roughly 0.4‰ yr−1 (0.6 Pg C yr−1). Although there were considerable inter-model differences, the rates are still remarkable in terms of their potential for mitigation of global warming. The disparate results among ESMs imply that key parameters that control processes such as SOC residence time need to be better constrained and that more comprehensive representation of land management impacts on soils remain critical for understanding the long-term potential of soils to sequester carbon.

Scenarios that limit global warming to below 2 degrees Centigrade by 2100 assume significant land-use change to support large-scale carbon dioxide (CO2) removal from the atmosphere by afforestation/reforestation, avoided deforestation, and Biomass Energy with Carbon Capture and Storage (BECCS). The more ambitious mitigation scenarios require even greater land area for mitigation and/or earlier adoption of CO2 removal strategies. Here we show that additional land-use change to meet a 1.5 degrees Centigrade climate change target could result in net losses of carbon from the land. The effectiveness of BECCS strongly depends on several assumptions related to the choice of biomass, the fate of initial above ground biomass, and the fossil-fuel emissions offset in the energy system. Depending on these factors, carbon removed from the atmosphere through BECCS could easily be offset by losses due to land-use change. If BECCS involves replacing high-carbon content ecosystems with crops, then forest-based mitigation could be more efficient for atmospheric CO2 removal than BECCS.

We present the compatible CO₂ emissions from fossil fuel (FF) burning and industry, calculated from the historical and Shared Socioeconomic Pathway (SSP) experiments of nine Earth system models (ESMs) participating in phase 6 of the Coupled Model Intercomparison Project (CMIP6). The multimodel mean FF emissions match the historical record well and are close to the data-based estimate of cumulative emissions (394 ± 59 GtC vs 400 ± 20 GtC, respectively). Only two models fall inside the observed uncertainty range; while two exceed the upper bound, five fall slightly below the lower bound, due primarily to the plateau in CO₂ concentration in the 1940s. The ESMs’ diagnosed FF emission rates are consistent with those generated by the integrated assessment models (IAMs) from which the SSPs’ CO₂ concentration pathways were constructed; the simpler IAMs’ emissions lie within the ESMs’ spread for seven of the eight SSP experiments, the other being only marginally lower, providing confidence in the relationship between the IAMs’ FF emission rates and concentration pathways. The ESMs require fossil fuel emissions to reduce to zero and subsequently become negative in SSP1-1.9, SSP1-2.6, SSP4-3.4, and SSP5-3.4over. We also present the ocean and land carbon cycle responses of the ESMs in the historical and SSP scenarios. The models’ ocean carbon cycle responses are in close agreement, but there is considerable spread in their land carbon cycle responses. Land-use and land-cover change emissions have a strong influence over the magnitude of diagnosed fossil fuel emissions, with the suggestion of an inverse relationship between the two.

Despite recent research identifying a clear anthropogenic impact on glacier recession, the effect of recent climate change on glacier-related hazards is at present unclear. Here we present the first global spatio-temporal assessment of glacial lake outburst floods (GLOFs) focusing explicitly on lake drainage following moraine dam failure. These floods occur as mountain glaciers recede and downwaste. GLOFs can have an enormous impact on downstream communities and infrastructure. Our assessment of GLOFs associated with the rapid drainage of moraine-dammed lakes provides insights into the historical trends of GLOFs and their distributions under current and future global climate change. We observe a clear global increase in GLOF frequency and their regularity around 1930, which likely represents a lagged response to post-Little Ice Age warming. Notably, we also show that GLOF frequency and regularity – rather unexpectedly – have declined in recent decades even during a time of rapid glacier recession. Although previous studies have suggested that GLOFs will increase in response to climate warming and glacier recession, our global results demonstrate that this has not yet clearly happened. From an assessment of the timing of climate forcing, lag times in glacier recession, lake formation and moraine-dam failure, we predict increased GLOF frequencies during the next decades and into the 22nd century.

The Paris agreement aims to hold global warming to well below 2 ∘C and to pursue efforts to limit it to 1.5 ∘C relative to the pre-industrial period. Recent estimates based on population growth and intended carbon emissions from participant countries suggest global warming may exceed this ambitious target. Here we present glacier volume projections for the end of this century, under a range of high-end climate change scenarios, defined as exceeding +2 ∘C global average warming relative to the pre-industrial period. Glacier volume is modelled by developing an elevation-dependent mass balance model for the Joint UK Land Environment Simulator (JULES). To do this, we modify JULES to include glaciated and unglaciated surfaces that can exist at multiple heights within a single grid box. Present-day mass balance is calibrated by tuning albedo, wind speed, precipitation, and temperature lapse rates to obtain the best agreement with observed mass balance profiles. JULES is forced with an ensemble of six Coupled Model Intercomparison Project Phase 5 (CMIP5) models, which were downscaled using the high-resolution HadGEM3-A atmosphere-only global climate model. The CMIP5 models use the RCP8.5 climate change scenario and were selected on the criteria of passing 2 ∘C global average warming during this century. The ensemble mean volume loss at the end of the century plus or minus 1 standard deviation is -64±5 % for all glaciers excluding those on the peripheral of the Antarctic ice sheet. The uncertainty in the multi-model mean is rather small and caused by the sensitivity of HadGEM3-A to the boundary conditions supplied by the CMIP5 models. The regions which lose more than 75 % of their initial volume by the end of the century are Alaska, western Canada and the US, Iceland, Scandinavia, the Russian Arctic, central Europe, Caucasus, high-mountain Asia, low latitudes, southern Andes, and New Zealand. The ensemble mean ice loss expressed in sea level equivalent contribution is 215.2±21.3 mm. The largest contributors to sea level rise are Alaska (44.6±1.1 mm), Arctic Canada north and south (34.9±3.0 mm), the Russian Arctic (33.3±4.8 mm), Greenland (20.1±4.4), high-mountain Asia (combined central Asia, South Asia east and west), (18.0±0.8 mm), southern Andes (14.4±0.1 mm), and Svalbard (17.0±4.6 mm). Including parametric uncertainty in the calibrated mass balance parameters gives an upper bound global volume loss of 281.1 mm of sea level equivalent by the end of the century. Such large ice losses will have inevitable consequences for sea level rise and for water supply in glacier-fed river systems.

Future climate change and increasing atmospheric CO2 are expected to cause major changes in vegetation structure and function over large fractions of the global land surface. Seven global vegetation models are used to analyze possible responses to future climate simulated by a range of general circulation models run under all four representative concentration pathway scenarios of changing concentrations of greenhouse gases. All 110 simulations predict an increase in global vegetation carbon to 2100, but with substantial variation between vegetation models. For example, at 4 °C of global land surface warming (510–758 ppm of CO2), vegetation carbon increases by 52–477 Pg C (224 Pg C mean), mainly due to CO2 fertilization of photosynthesis. Simulations agree on large regional increases across much of the boreal forest, western Amazonia, central Africa, western China, and southeast Asia, with reductions across southwestern North America, central South America, southern Mediterranean areas, southwestern Africa, and southwestern Australia. Four vegetation models display discontinuities across 4 °C of warming, indicating global thresholds in the balance of positive and negative influences on productivity and biomass. In contrast to previous global vegetation model studies, we emphasize the importance of uncertainties in projected changes in carbon residence times. We find, when all seven models are considered for one representative concentration pathway × general circulation model combination, such uncertainties explain 30% more variation in modeled vegetation carbon change than responses of net primary productivity alone, increasing to 151% for non-HYBRID4 models. A change in research priorities away from production and toward structural dynamics and demographic processes is recommended.

Achieving the long-term temperature goal of the Paris Agreement requires forest-based mitigation. Collective progress towards this goal will be assessed by the Paris Agreement’s Global stocktake. At present, there is a discrepancy of about 4 GtCO2 yr−1 in global anthropogenic net land-use emissions between global models (reflected in IPCC assessment reports) and aggregated national GHG inventories (under the UNFCCC). We show that a substantial part of this discrepancy (about 3.2 GtCO2 yr−1) can be explained by conceptual differences in anthropogenic forest sink estimation, related to the representation of environmental change impacts and the areas considered as managed. For a more credible tracking of collective progress under the Global stocktake, these conceptual differences between models and inventories need to be reconciled. We implement a new method of disaggregation of global land model results that allows greater comparability with GHG inventories. This provides a deeper understanding of model–inventory differences, allowing more transparent analysis of forest-based mitigation and facilitating a more accurate Global stocktake.