Explore the vast range of titles available.

We quantify an elevated occurrence of abrupt changes in ocean environmental conditions under human-induced climate forcing using Earth system model output through a novel analysis method that compares the temporal evolution of the forcings applied with the development of local ocean state changes for temperature, oxygen concentration, and carbonate ion concentration. Through a multi-centennial Earth system model experiment, we show that such an increase is not fully reversible after excess greenhouse gas emissions go back to zero. The increase in occurrence of regional abrupt changes in marine environmental conditions has not yet been accounted for adequately in climate impact analyses that usually associate ecosystem shifts large-scale variability or extreme events. Estimates for remaining greenhouse gas emission targets need thus to be more conservative.

The globally integrated sea–air anthropogenic carbon dioxide (CO2) flux from 1990 to 2009 is determined from models and data-based approaches as part of the Regional Carbon Cycle Assessment and Processes (RECCAP) project. Numerical methods include ocean inverse models, atmospheric inverse models, and ocean general circulation models with parameterized biogeochemistry (OBGCMs). The median value of different approaches shows good agreement in average uptake. The best estimate of anthropogenic CO2 uptake for the time period based on a compilation of approaches is −2.0 Pg C yr−1. The interannual variability in the sea–air flux is largely driven by large-scale climate re-organizations and is estimated at 0.2 Pg C yr−1 for the two decades with some systematic differences between approaches. The largest differences between approaches are seen in the decadal trends. The trends range from −0.13 (Pg C yr−1) decade−1 to −0.50 (Pg C yr−1) decade−1 for the two decades under investigation. The OBGCMs and the data-based sea–air CO2 flux estimates show appreciably smaller decadal trends than estimates based on changes in carbon inventory suggesting that methods capable of resolving shorter timescales are showing a slowing of the rate of ocean CO2 uptake. RECCAP model outputs for five decades show similar differences in trends between approaches.

The recently developed Norwegian Earth System Model (NorESM) is employed for simulations contributing to the CMIP5 (Coupled Model Intercomparison Project phase 5) experiments and the fifth assessment report of the Intergovernmental Panel on Climate Change (IPCC-AR5). In this manuscript, we focus on evaluating the ocean and land carbon cycle components of the NorESM, based on the preindustrial control and historical simulations. Many of the observed large scale ocean biogeochemical features are reproduced satisfactorily by the NorESM. When compared to the climatological estimates from the World Ocean Atlas (WOA), the model simulated temperature, salinity, oxygen, and phosphate distributions agree reasonably well in both the surface layer and deep water structure. However, the model simulates a relatively strong overturning circulation strength that leads to noticeable model-data bias, especially within the North Atlantic Deep Water (NADW). This strong overturning circulation slightly distorts the structure of the biogeochemical tracers at depth. Advancements in simulating the oceanic mixed layer depth with respect to the previous generation model particularly improve the surface tracer distribution as well as the upper ocean biogeochemical processes, particularly in the Southern Ocean. Consequently, near-surface ocean processes such as biological production and air-sea gas exchange, are in good agreement with climatological observations. The NorESM adopts the same terrestrial model as the Community Earth System Model (CESM1). It reproduces the general pattern of land-vegetation gross primary productivity (GPP) when compared to the observationally based values derived from the FLUXNET network of eddy covariance towers. While the model simulates well the vegetation carbon pool, the soil carbon pool is smaller by a factor of three relative to the observational based estimates. The simulated annual mean terrestrial GPP and total respiration are slightly larger than observed, but the difference between the global GPP and respiration is comparable. Model-data bias in GPP is mainly simulated in the tropics (overestimation) and in high latitudes (underestimation). Within the NorESM framework, both the ocean and terrestrial carbon cycle models simulate a steady increase in carbon uptake from the preindustrial period to the present-day. The land carbon uptake is noticeably smaller than the observations, which is attributed to the strong nitrogen limitation formulated by the land model.

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of the El Niño–Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in CMIP5: (1) the simulated global burned area being less than half of that of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past 2 decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere mid-latitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet–dry conditions, indicating the need to improve fuel wetness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest a post-processing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products.

Idealised and hindcast simulations performed with the stand-alone ocean carbon-cycle configuration of the Norwegian Earth System Model (NorESM-OC) are described and evaluated. We present simulation results of three different model configurations (two different model versions at different grid resolutions) using two different atmospheric forcing data sets. Model version NorESM-OC1 corresponds to the version that is included in the NorESM-ME1 fully coupled model, which participated in CMIP5. The main update between NorESM-OC1 and NorESM-OC1.2 is the addition of two new options for the treatment of sinking particles. We find that using a constant sinking speed, which has been the standard in NorESM's ocean carbon cycle module HAMOCC (HAMburg Ocean Carbon Cycle model), does not transport enough particulate organic carbon (POC) into the deep ocean below approximately 2000 m depth. The two newly implemented parameterisations, a particle aggregation scheme with prognostic sinking speed, and a simpler scheme that uses a linear increase in the sinking speed with depth, provide better agreement with observed POC fluxes. Additionally, reduced deep ocean biases of oxygen and remineralised phosphate indicate a better performance of the new parameterisations. For model version 1.2, a re-tuning of the ecosystem parameterisation has been performed, which (i) reduces previously too high primary production at high latitudes, (ii) consequently improves model results for surface nutrients, and (iii) reduces alkalinity and dissolved inorganic carbon biases at low latitudes. We use hindcast simulations with prescribed observed and constant (pre-industrial) atmospheric CO2 concentrations to derive the past and contemporary ocean carbon sink. For the period 1990-1999 we find an average ocean carbon uptake ranging from 2.01 to 2.58 Pg C yr1 depending on model version, grid resolution, and atmospheric forcing data set.

Chemical state analyses of the atmosphere based on data assimilation may be degraded by inconsistent covariances of background and observation errors. An efficient method to calculate consistency diagnostics for background and observation errors in observation space is applied to analyses of the four‐dimensional variational stratospheric chemistry data assimilation system SACADA (Synoptic Analysis of Chemical Constituents by Advanced Data Assimilation). A background error covariance model for the assimilation of Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) ozone retrievals is set up and optimized. It is shown that a significant improvement of the assimilation system performance is attained through the use of this covariance model compared to a simple covariance formulation, which assumes background errors to be a fixed fraction of the field value. The forecast skill, measured by the distance between the model forecast and MIPAS observations, is shown to improve. Further, an evaluation of analyses with independent data from the Halogen Observation Experiment (HALOE), the Stratospheric Aerosol and Gas Experiment II (SAGE II), and ozone sondes reveals that the standard deviation of ozone analyses with respect to these instruments is reduced throughout the middle stratosphere. Compared to the impact of background error variances on analysis quality, it is found that the precise specification of spatial background error correlations appears to be less critical if observations are spatially and temporally dense. Results indicate that ozone forecast errors of a state of the art stratospheric chemistry assimilation system are of the same order of magnitude as MIPAS observation errors.

A novel stratospheric chemical data assimilation system has been developed and applied to Environmental Satellite Michelson Interferometer for Passive Atmospheric Sounding (ENVISAT/MIPAS) data, aiming to combine the sophistication of the four‐dimensional variational (4D‐var) technique with flow‐dependent covariance modeling and also to improve numerical performance. The system is tailored for operational stratospheric chemistry state monitoring. The atmospheric model of the assimilation system includes a state‐of‐the‐art stratospheric chemistry transport module along with its adjoint and the German weather service's global meteorological forecast model, providing meteorological parameters. Both models share the same grid and same advection time step, to ensure dynamic consistency without spatial and temporal interpolation errors. A notable numerical efficiency gain is obtained through an icosahedral grid. As a novel feature in stratospheric variational data assimilation a special focus was placed on an optimal spatial exploitation of satellite data by dynamic formulation of the forecast error covariance matrix, providing potential vorticity controlled anisotropic and inhomogeneous influence radii. In this first part of the study the design and numerical features of the data assimilation system is presented, along with analyses of two case studies and a posteriori validation. Assimilated data include retrievals of O3, CH4, N2O, NO2, HNO3, and water vapor. The analyses are compared with independent observations provided by Stratospheric Aerosol and Gas Experiment II (SAGE II) and Halogen Occultation Experiment (HALOE) retrievals. It was found that there are marked improvements for both analyses and assimilation based forecasts when compared with control model runs without any data ingestion.

Anthropogenic climate change is projected to lead to ocean warming, acidification, deoxygenation, reductions in near-surface nutrients, and changes to primary production, all of which are expected to affect marine ecosystems. Here we assess projections of these drivers of environmental change over the twenty-first century from Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) that were forced under the CMIP6 Shared Socioeconomic Pathways (SSPs). Projections are compared to those from the previous generation (CMIP5) forced under the Representative Concentration Pathways (RCPs). A total of 10 CMIP5 and 13 CMIP6 models are used in the two multi-model ensembles. Under the high-emission scenario SSP5-8.5, the multi-model global mean change (2080–2099 mean values relative to 1870–1899) ± the inter-model SD in sea surface temperature, surface pH, subsurface (100–600 m) oxygen concentration, euphotic (0–100 m) nitrate concentration, and depth-integrated primary production is +3.47±0.78 ∘C, -0.44±0.005, -13.27±5.28, -1.06±0.45 mmol m−3 and -2.99±9.11 %, respectively. Under the low-emission, high-mitigation scenario SSP1-2.6, the corresponding global changes are +1.42±0.32 ∘C, -0.16±0.002, -6.36±2.92, -0.52±0.23 mmol m−3, and -0.56±4.12 %. Projected exposure of the marine ecosystem to these drivers of ocean change depends largely on the extent of future emissions, consistent with previous studies. The ESMs in CMIP6 generally project greater warming, acidification, deoxygenation, and nitrate reductions but lesser primary production declines than those from CMIP5 under comparable radiative forcing. The increased projected ocean warming results from a general increase in the climate sensitivity of CMIP6 models relative to those of CMIP5. This enhanced warming increases upper-ocean stratification in CMIP6 projections, which contributes to greater reductions in upper-ocean nitrate and subsurface oxygen ventilation. The greater surface acidification in CMIP6 is primarily a consequence of the SSPs having higher associated atmospheric CO2 concentrations than their RCP analogues for the same radiative forcing. We find no consistent reduction in inter-model uncertainties, and even an increase in net primary production inter-model uncertainties in CMIP6, as compared to CMIP5.

The second version of the coupled Norwegian Earth System Model (NorESM2) is presented and evaluated. NorESM2 is based on the second version of the Community Earth System Model (CESM2) and shares with CESM2 the computer code infrastructure and many Earth system model components. However, NorESM2 employs entirely different ocean and ocean biogeochemistry models. The atmosphere component of NorESM2 (CAM-Nor) includes a different module for aerosol physics and chemistry, including interactions with cloud and radiation; additionally, CAM-Nor includes improvements in the formulation of local dry and moist energy conservation, in local and global angular momentum conservation, and in the computations for deep convection and air–sea fluxes. The surface components of NorESM2 have minor changes in the albedo calculations and to land and sea-ice models.We present results from simulations with NorESM2 that were carried out for the sixth phase of the Coupled Model Intercomparison Project (CMIP6). Two versions of the model are used: one with lower (∼ 2∘) atmosphere–land resolution and one with medium (∼ 1∘) atmosphere–land resolution. The stability of the pre-industrial climate and the sensitivity of the model to abrupt and gradual quadrupling of CO2 are assessed, along with the ability of the model to simulate the historical climate under the CMIP6 forcings. Compared to observations and reanalyses, NorESM2 represents an improvement over previous versions of NorESM in most aspects. NorESM2 appears less sensitive to greenhouse gas forcing than its predecessors, with an estimated equilibrium climate sensitivity of 2.5 K in both resolutions on a 150-year time frame; however, this estimate increases with the time window and the climate sensitivity at equilibration is much higher. We also consider the model response to future scenarios as defined by selected Shared Socioeconomic Pathways (SSPs) from the Scenario Model Intercomparison Project defined under CMIP6. Under the four scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5), the warming in the period 2090–2099 compared to 1850–1879 reaches 1.3, 2.2, 3.0, and 3.9 K in NorESM2-LM, and 1.3, 2.1, 3.1, and 3.9 K in NorESM-MM, robustly similar in both resolutions. NorESM2-LM shows a rather satisfactory evolution of recent sea-ice area. In NorESM2-LM, an ice-free Arctic Ocean is only avoided in the SSP1-2.6 scenario.

The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is -0.36 to 0.29 ∘C, with a model ensemble mean of -0.07 ∘C, median of -0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.