Transient Climate Response Research Articles

Machine learning of cloud types in satellite observations and climate models

Abstract. Uncertainty in cloud feedbacks in climate models is a major limitation in projections of future climate. Therefore, evaluation and improvement of cloud simulation are essential to ensure the accuracy of climate models. We analyse cloud biases and cloud change with respect to global mean near-surface temperature (GMST) in climate models relative to satellite observations and relate them to equilibrium climate sensitivity, transient climate response and cloud feedback. For this purpose, we develop a supervised deep convolutional artificial neural network for determination of cloud types from low-resolution (2.5∘×2.5∘) daily mean top-of-atmosphere shortwave and longwave radiation fields, corresponding to the World Meteorological Organization (WMO) cloud genera recorded by human observers in the Global Telecommunication System (GTS). We train this network on top-of-atmosphere radiation retrieved by the Clouds and the Earth’s Radiant Energy System (CERES) and GTS and apply it to the Coupled Model Intercomparison Project Phase 5 and 6 (CMIP5 and CMIP6) model output and the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) and the Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2) reanalyses. We compare the cloud types between models and satellite observations. We link biases to climate sensitivity and identify a negative linear relationship between the root mean square error of cloud type occurrence derived from the neural network and model equilibrium climate sensitivity (ECS), transient climate response (TCR) and cloud feedback. This statistical relationship in the model ensemble favours models with higher ECS, TCR and cloud feedback. However, this relationship could be due to the relatively small size of the ensemble used or decoupling between present-day biases and future projected cloud change. Using the abrupt-4×CO2 CMIP5 and CMIP6 experiments, we show that models simulating decreasing stratiform and increasing cumuliform clouds tend to have higher ECS than models simulating increasing stratiform and decreasing cumuliform clouds, and this could also partially explain the association between the model cloud type occurrence error and model ECS.

Much of zero emissions commitment occurs before reaching net zero emissions

We explore the response of the Earth’s coupled climate and carbon system to an idealized sequential addition and removal of CO2 to the atmosphere, following a symmetric and continuous emissions pathway, in contrast to the discontinuous emissions pathways that have largely informed our understanding of the climate response to net zero and net negative emissions to date. We find, using both an Earth system model and an ensemble of simple climate model realizations, that warming during the emissions reduction and negative emissions phases is defined by a combination of a proportionality of warming to cumulative emissions characterized by the transient climate response to emissions (TCRE), and a deviation from that proportionality that is governed by the zero emissions commitment (ZEC). About half of the ZEC is realized before reaching zero emissions, and the ZEC thus also controls the timing between peak cumulative CO2 emissions and peak temperature, such that peak temperature may occur before peak cumulative emissions if ZEC is negative, underscoring the importance of ZEC in climate policies aimed to limit peak warming. Thus we argue that ZEC is better defined as the committed warming relative to the expected TCRE proportionality, rather than as the additional committed warming that will occur after reaching net zero CO2 emissions. Once established, the combined TCRE and ZEC relationship holds almost to complete removal of prior cumulative CO2 emissions. As cumulative CO2 emissions approach zero through negative CO2 emissions, CO2 concentrations drop below preindustrial values, while residual long-term climate change continues, governed by multicentennial dynamical processes.

Inter‐Model Differences in Future Summer Onset Over the Northern High Latitudes

AbstractClimate models consistently project a significant advance of summer onset over the Northern high latitude lands in the late 21st century. However, large inter‐model differences exist, in particular, over Northwest Russia, Bering Sea region and Northeastern Canada. This study quantifies the influences of background global warming rate and regional snow‐albedo feedback on the inter‐model spreads in summer onset projections. Linear regression analysis reveals that models' transient climate response (TCR) values explain around 39%–50% of regional inter‐model variances, indicating the dominant contribution of models' global warming rate. When removing TCR influences, the summer onset and snow cover exhibit significant inter‐model correlations over Bering Sea region and Northeastern Canada, explaining 25%–30% of the inter‐model variances. These results suggest that global climate sensitivity as well as regional snow‐albedo feedback need to be considered for reliable projection of future summer season timing at regional scales.

Through an Ice Sheet, Darkly

Through an Ice Sheet, Darkly

Impact of negative and positive CO2 emissions on global warming metrics using an ensemble of Earth system model simulations

Abstract. The benefits of implementing negative emission technologies in the global warming response to cumulative carbon emissions until the year 2420 are assessed following the shared socioeconomic pathway (SSP) 1-2.6, the sustainable development scenario, with a comprehensive set of intermediate-complexity Earth system model integrations. Model integrations include 86 different model realisations covering a wide range of plausible climate states. The global warming response is assessed in terms of two key climate metrics: the effective transient climate response to cumulative CO2 emissions (eTCRE), measuring the surface warming response to cumulative carbon emissions and associated non-CO2 forcing, and the effective zero emissions commitment (eZEC), measuring the extent of any continued warming after net-zero CO2 emissions are reached. The transient climate response to cumulative CO2 emissions (TCRE) is estimated as 2.2 K EgC−1 (median value) with a 10 %–90 % range of 1.75 to 3.13 K EgC−1 in 2100, approximated from the eTCRE by removing the contribution of non-CO2 forcing. During the positive emission phase, the eTCRE decreases from 2.71 (2.0 to 3.65) to 2.61 (1.91 to 3.62) K EgC−1 due to a weakening in the dependence of radiative forcing on atmospheric carbon, which is partly opposed by an increasing fraction of the radiative forcing warming the surface as the ocean stratifies. During the net negative and zero emission phases, a progressive reduction in the eTCRE to 2.0 (1.39 to 2.96) K EgC−1 is driven by the reducing airborne fraction as atmospheric CO2 is drawn down mainly by the ocean. The model uncertainty in the slopes of warming versus cumulative CO2 emissions varies from being controlled by the radiative feedback parameter during positive emissions to being affected by carbon-cycle parameters during net negative emissions, consistent with the drivers of uncertainty diagnosed from the coefficient of variation of the contributions in the eTCRE framework. The continued warming after CO2 emissions cease and remain at zero gives a model mean eZEC of −0.03 K after 25 years, which decreases in time to −0.21 K at 90 years after emissions cease. However, there is a spread in the ensemble with a temperature overshoot occurring in 20 % of the ensemble members at 25 years after cessation of emissions. If net negative emissions are included, there is a reduction in atmospheric CO2 and there is a decrease in temperature overshoot so that the eZEC is positive in only 5 % of the ensemble members. Hence, incorporating negative emissions enhances the ability to meet climate targets and avoid risk of continued warming after net zero is reached.

Future Climate Change Under SSP Emission Scenarios With GISS‐E2.1

AbstractThis paper presents the response to anthropogenic forcing in the GISS‐E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non‐interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One‐Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1‐G) and HYCOM (E2.1‐H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1‐2.6, SSP2‐4.5, SSP4‐6.0, and SSP5‐8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high‐mitigation scenarios SSP1‐1.9 and SSP1‐2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1‐H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5‐8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1‐G (−1.21 × 106 km2/°C) and E2.1‐H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5‐8.5 in the E2.1‐G model, comparable to the reduction in the corresponding CMIP5 GISS‐E2 RCP8.5 simulation. Both low‐end scenarios SSP1‐1.9 and SSP1‐2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century).

A defense of usable climate mitigation science: how science can contribute to social movements

Much of modern climate science is motivated by the problem of human-caused climate change or its potential solutions, and aims to be “usable” for relevant stakeholders. Sobel (2021) argues in this issue that the expectation for improved climate projections to drive mitigation of greenhouse gas emissions “can now be understood [as] naive about the role of politics, and the power of entrenched interests to inhibit climate action.” While he criticizes this “linear model” of the science–policy interface; he does not elaborate on alternative avenues for scientific advances to spur mitigation. Instead, he encourages physical climate scientists who wish to produce usable results to orient their work towards informing adaptation. He argues that, relative to mitigation science, adaptation science is more likely to be used by stakeholders because the remaining scientific uncertainties are larger and the social barriers to implementation are lower. We join Sobel in calling on physical climate scientists to reflect upon the pathways through which their research improves societal welfare. However, we argue that Sobel’s argument overlooks an important theory of change, namely that mitigation science is politically usable through the non-linear dynamics of social mobilization. Social theories of policy change suggest that organized groups play an outsized role in setting the policy agenda. Grassroots activism on climate, however, has historically been hindered by the abstract nature of climate change, paling in comparison to the lobbying and misinformation campaigns funded by the vested fossil fuel interests. We describe how two recent advances in mitigation science, the Transient Climate Response to cumulative Emissions (TCRE) and Extreme Event Attribution (EEA), have provided social movements with information that allows them to re-frame the climate change problem in a way that attributes blame for the problem, motivates collective action across a diverse coalition of stakeholders, and could plausibly compel policymakers to prioritize the issue in the coming years. Given the utmost importance of mitigation in preventing climate change at the source, we thus advocate for a broader agenda of usable climate research that includes co-production of both mitigation and adaptation science.

Does Model Calibration Reduce Uncertainty in Climate Projections?

Abstract Uncertainty in climate projections is large as shown by the likely uncertainty ranges in equilibrium climate sensitivity (ECS) of 2.5–4 K and in the transient climate response (TCR) of 1.4–2.2 K. Uncertainty in model projections could arise from the way in which unresolved processes are represented, the parameter values used, or the targets for model calibration. We show that, in two climate model ensembles that were objectively calibrated to minimize differences from observed large-scale atmospheric climatology, uncertainties in ECS and TCR are about 2–6 times smaller than in the CMIP5 or CMIP6 multimodel ensemble. We also find that projected uncertainties in surface temperature, precipitation, and annual extremes are relatively small. Residual uncertainty largely arises from unconstrained sea ice feedbacks. The more than 20-year-old HadAM3 standard model configuration simulates observed hemispheric-scale observations and preindustrial surface temperatures about as well as the median CMIP5 and CMIP6 ensembles while the optimized configurations simulate these better than almost all the CMIP5 and CMIP6 models. Hemispheric-scale observations and preindustrial temperatures are not systematically better simulated in CMIP6 than in CMIP5 although the CMIP6 ensemble seems to better simulate patterns of large-scale observations than the CMIP5 ensemble and the optimized HadAM3 configurations. Our results suggest that most CMIP models could be improved in their simulation of large-scale observations by systematic calibration. However, the uncertainty in climate projections (for a given scenario) likely largely arises from the choice of parameterization schemes for unresolved processes (“structural uncertainty”), with different tuning targets another possible contributor. Significance Statement Climate models represent unresolved phenomena controlled by uncertain parameters. Changes in these parameters impact how well a climate model simulates current climate and its climate projections. Multiple calibrations of a single climate model, using an objective method, to large-scale atmospheric observations are performed. These models produce very similar climate projections at both global and regional scales. An analysis that combines uncertainties in observations with simulated sensitivity to observations and climate response also has small uncertainty showing that, for this model, current observations constrain climate projections. Recently developed climate models have a broad range of abilities to simulate large-scale climate with only some improvement in their ability to simulate this despite a decade of model development.

CMIP6 Simulations With the CMCC Earth System Model (CMCC‐ESM2)

AbstractThis article introduces the second generation CMCC Earth System Model (CMCC‐ESM2) that extends a number of marine and terrestrial biogeochemical processes with respect to its CMIP5 predecessor. In particular, land biogeochemistry was extended to a wider set of carbon pools and plant functional types, along with a prognostic representation of the nitrogen cycle. The marine ecosystem representation was reshaped toward an intermediate complexity of lower trophic level interactions, including an interactive benthic compartment and a new formulation of heterotrophic bacterial population. Details are provided on the model setup and implementation for the different experiments performed as contribution to the sixth phase of the Coupled Model Intercomparison Project. CMCC‐ESM2 shows an equilibrium climate sensitivity of 3.57°C and a transient climate response of 1.97°C which are close to the CMIP5 and CMIP6 multi‐model averages. The evaluation of the coupled climate‐carbon response in the historical period against available observational datasets show a consistent representation of both physical and biogeochemical quantities. However, the land carbon sink is found to be weaker than the current global carbon estimates and the simulated marine primary production is slightly below the satellite‐based average over recent decades. Future projections coherently show a prominent global warming over the northern hemisphere with intensified precipitations at high latitudes. The expected ranges of variability for oceanic pH and oxygen, as well as land carbon and nitrogen soil storage, compare favorably with those assessed from other CMIP6 models.

Transient Influence of the Reduction of Deepwater Formation on Ocean Heat Uptake and Heat Budgets in the Global Climate System

AbstractThe formation and spreading of dense deepwater in the polar regions play a key role in one of the most important climate systems, namely ocean meridional overturning circulation, and the deepwater formation is projected to decrease under the global warming. However, the impact of the reduced deepwater formation on the climate system has not been explored in detail. Here, we performed a series of numerical experiments with a climate model where the downward water mass transport through the bottom boundary layer is artificially reduced to quantitatively evaluate its impacts on the transient ocean and climate responses. It is demonstrated that changes in deepwater formation have non‐negligible impacts on not only ocean heat content but also the Earth's radiation budget at the top of the atmosphere: reduction in deepwater formation in high‐latitude oceans causes warming of bottom water, cooling of the ocean surface, and a subsequent decrease in outgoing longwave radiation.

The fractional energy balance equation for climate projections through 2100

Abstract. We produce climate projections through the 21st century using the fractional energy balance equation (FEBE): a generalization of the standard energy balance equation (EBE). The FEBE can be derived from Budyko–Sellers models or phenomenologically through the application of the scaling symmetry to energy storage processes, easily implemented by changing the integer order of the storage (derivative) term in the EBE to a fractional value. The FEBE is defined by three parameters: a fundamental shape parameter, a timescale and an amplitude, corresponding to, respectively, the scaling exponent h, the relaxation time τ and the equilibrium climate sensitivity (ECS). Two additional parameters were needed for the forcing: an aerosol recalibration factor α to account for the large aerosol uncertainty and a volcanic intermittency correction exponent ν. A Bayesian framework based on historical temperatures and natural and anthropogenic forcing series was used for parameter estimation. Significantly, the error model was not ad hoc but rather predicted by the model itself: the internal variability response to white noise internal forcing. The 90 % credible interval (CI) of the exponent and relaxation time were h=[0.33, 0.44] (median = 0.38) and τ=[2.4, 7.0] (median = 4.7) years compared to the usual EBE h=1, and literature values of τ typically in the range 2–8 years. Aerosol forcings were too strong, requiring a decrease by an average factor α=[0.2, 1.0] (median = 0.6); the volcanic intermittency correction exponent was ν=[0.15, 0.41] (median = 0.28) compared to standard values α=ν=1. The overpowered aerosols support a revision of the global modern (2005) aerosol forcing 90 % CI to a narrower range [−1.0, −0.2] W m−2. The key parameter ECS in comparison to IPCC AR5 (and to the CMIP6 MME), the 90 % CI range is reduced from [1.5, 4.5] K ([2.0, 5.5] K) to [1.6, 2.4] K ([1.5, 2.2] K), with median value lowered from 3.0 K (3.7 K) to 2.0 K (1.8 K). Similarly we found for the transient climate response (TCR), the 90 % CI range shrinks from [1.0, 2.5] K ([1.2, 2.8] K) to [1.2, 1.8] K ([1.1, 1.6] K) and the median estimate decreases from 1.8 K (2.0 K) to 1.5 K (1.4 K). As often seen in other observational-based studies, the FEBE values for climate sensitivities are therefore somewhat lower but still consistent with those in IPCC AR5 and the CMIP6 MME. Using these parameters, we made projections to 2100 using both the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP) scenarios, and compared them to the corresponding CMIP5 and CMIP6 multi-model ensembles (MMEs). The FEBE historical reconstructions (1880–2020) closely follow observations, notably during the 1998–2014 slowdown (“hiatus”). We also reproduce the internal variability with the FEBE and statistically validate this against centennial-scale temperature observations. Overall, the FEBE projections were 10 %–15 % lower but due to their smaller uncertainties, their 90 % CIs lie completely within the GCM 90 % CIs. This agreement means that the FEBE validates the MME, and vice versa.

Biased Estimates of Equilibrium Climate Sensitivity and Transient Climate Response Derived From Historical CMIP6 Simulations

AbstractThis study assesses the effective climate sensitivity (EffCS) and transient climate response (TCR) derived from global energy budget constraints within historical simulations of eight CMIP6 global climate models (GCMs). These calculations are enabled by use of the Radiative Forcing Model Intercomparison Project (RFMIP) simulations, which permit accurate quantification of the radiative forcing. Long‐term historical energy budget constraints generally underestimate EffCS from CO2 quadrupling and TCR from CO2 ramping, owing to changes in radiative feedbacks and changes in ocean heat uptake efficiency. Atmospheric GCMs forced by observed warming patterns produce lower values of EffCS that are more in line with those inferred from observed historical energy budget changes. The differences in the EffCS estimates from historical energy budget constraints of models and observations are traced to discrepancies between modeled and observed historical surface warming patterns.

Ocean Model Formulation Influences Transient Climate Response

AbstractThe transient climate response (TCR) is 20% higher in the Alfred Wegener Institute Climate Model (AWI‐CM) compared to the Max Planck Institute Earth System Model (MPI‐ESM) whereas the equilibrium climate sensitivity (ECS) is by up to 10% higher in AWI‐CM. These results are largely independent of the two considered model resolutions for each model. The two coupled CMIP6 models share the same atmosphere‐land component ECHAM6.3 developed at the Max Planck Institute for Meteorology (MPI‐M). However, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice‐ocean model developed at MPI‐M and the FESOM sea ice‐ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is related to ocean heat uptake in response to greenhouse gas forcing. Specifically, AWI‐CM simulations show stronger surface heating than MPI‐ESM simulations while the latter accumulate more heat in the deeper ocean. The vertically integrated ocean heat content is increasing slower in AWI‐CM model configurations compared to MPI‐ESM model configurations in the high latitudes. Weaker vertical mixing in AWI‐CM model configurations compared to MPI‐ESM model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell and Ross Gyres and the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI‐CM model configurations and the presence of a warming hole in MPI‐ESM model configurations. All these differences occur largely independent of the considered model resolutions.

Abrupt climate changes in the last two deglaciations simulated with different Northern ice sheet discharge and insolation

There were significant differences between the last two deglaciations, particularly in Atlantic Meridional Overturning Circulation (AMOC) and Antarctic warming in the deglaciations and the following interglacials. Here, we present transient simulations of deglaciation using a coupled atmosphere–ocean general circulation model for the last two deglaciations focusing on the impact of ice sheet discharge on climate changes associated with the AMOC in the first part, and the sensitivity studies using a Northern Hemisphere ice sheet model in the second part. We show that a set of abrupt climate changes of the last deglaciation, including Bolling–Allerod warming, the Younger Dryas, and onset of the Holocene were simulated with gradual changes of both ice sheet discharge and radiative forcing. On the other hand, penultimate deglaciation, with the abrupt climate change only at the beginning of the last interglacial was simulated when the ice sheet discharge was greater than in the last deglaciation by a factor of 1.5. The results, together with Northern Hemisphere ice sheet model experiments suggest the importance of the transient climate and AMOC responses to the different orbital forcing conditions of the last two deglaciations, through the mechanisms of mass loss of the Northern Hemisphere ice sheet and meltwater influx to the ocean.

Quantifying non-CO2 contributions to remaining carbon budgets

The IPCC Special Report on 1.5 °C concluded that anthropogenic global warming is determined by cumulative anthropogenic CO2 emissions and the non-CO2 radiative forcing level in the decades prior to peak warming. We quantify this using CO2-forcing-equivalent (CO2-fe) emissions. We produce an observationally constrained estimate of the Transient Climate Response to cumulative carbon Emissions (TCRE), giving a 90% confidence interval of 0.26–0.78 °C/TtCO2, implying a remaining total CO2-fe budget from 2020 to 1.5 °C of 350–1040 GtCO2-fe, where non-CO2 forcing changes take up 50 to 300 GtCO2-fe. Using a central non-CO2 forcing estimate, the remaining CO2 budgets are 640, 545, 455 GtCO2 for a 33, 50 or 66% chance of limiting warming to 1.5 °C. We discuss the impact of GMST revisions and the contribution of non-CO2 mitigation to remaining budgets, determining that reporting budgets in CO2-fe for alternative definitions of GMST, displaying CO2 and non-CO2 contributions using a two-dimensional presentation, offers the most transparent approach.

Enhanced hydrological cycle increases ocean heat uptake and moderates transient climate change.

The large-scale moistening of the atmosphere in response to increasing greenhouse gases amplifies the existing patterns of precipitation minus evaporation (P-E) which, in turn, amplifies the spatial contrast in sea surface salinity (SSS). Through a series of transient CO2 doubling experiments, we demonstrate that surface salinification driven by the amplified dry conditions (P-E < 0), primarily in the subtropical ocean, accelerates ocean heat uptake. The salinification also drives the sequestration of upper-level heat into the deeper ocean, reducing the thermal stratification and increasing the heat uptake through a positive feedback. The change in Atlantic Meridional Overturning Circulation due to salinification plays a secondary role in heat uptake. Consistent with the heat uptake changes, the transient climate response would increase by approximately 0.4 K without this process. Observed multi-decadal changes in subsurface temperature and salinity resembles those simulated, indicating that anthropogenically-forced changes in salinity are likely enhancing the ocean heat uptake.

EC-Earth3-AerChem: a global climate model with interactive aerosols and atmospheric chemistry participating in CMIP6

Abstract. This paper documents the global climate model EC-Earth3-AerChem, one of the members of the EC-Earth3 family of models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6). EC-Earth3-AerChem has interactive aerosols and atmospheric chemistry and contributes to the Aerosols and Chemistry Model Intercomparison Project (AerChemMIP). In this paper, we give an overview of the model, describe in detail how it differs from the other EC-Earth3 configurations, and outline the new features compared with the previously documented version of the model (EC-Earth 2.4). We explain how the model was tuned and spun up under preindustrial conditions and characterize the model's general performance on the basis of a selection of coupled simulations conducted for CMIP6. The net energy imbalance at the top of the atmosphere in the preindustrial control simulation is on average −0.09 W m−2 with a standard deviation due to interannual variability of 0.25 W m−2, showing no significant drift. The global surface air temperature in the simulation is on average 14.08 ∘C with an interannual standard deviation of 0.17 ∘C, exhibiting a small drift of 0.015 ± 0.005 ∘C per century. The model's effective equilibrium climate sensitivity is estimated at 3.9 ∘C, and its transient climate response is estimated at 2.1 ∘C. The CMIP6 historical simulation displays spurious interdecadal variability in Northern Hemisphere temperatures, resulting in a large spread across ensemble members and a tendency to underestimate observed annual surface temperature anomalies from the early 20th century onwards. The observed warming of the Southern Hemisphere is well reproduced by the model. Compared with the ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis version 5 (ERA5), the surface air temperature climatology for 1995–2014 has an average bias of −0.86 ± 0.05 ∘C with a standard deviation across ensemble members of 0.35 ∘C in the Northern Hemisphere and 1.29 ± 0.02 ∘C with a corresponding standard deviation of 0.05 ∘C in the Southern Hemisphere. The Southern Hemisphere warm bias is largely caused by errors in shortwave cloud radiative effects over the Southern Ocean, a deficiency of many climate models. Changes in the emissions of near-term climate forcers (NTCFs) have significant effects on the global climate from the second half of the 20th century onwards. For the SSP3-7.0 Shared Socioeconomic Pathway, the model gives a global warming at the end of the 21st century (2091–2100) of 4.9 ∘C above the preindustrial mean. A 0.5 ∘C stronger warming is obtained for the AerChemMIP scenario with reduced emissions of NTCFs. With concurrent reductions of future methane concentrations, the warming is projected to be reduced by 0.5 ∘C.

On Spurious Causality, CO2, and Global Temperature

Stips et al. (2016) use information flows (Liang (2008, 2014)) to establish causality from various forcings to global temperature. We show that the formulas being used hinge on a simplifying assumption that is nearly always rejected by the data. We propose the well-known forecast error variance decomposition based on a Vector Autoregression as an adequate measure of information flow, and find that most results in Stips et al. (2016) cannot be corroborated. Then, we discuss which modeling choices (e.g., the choice of CO2 series and assumptions about simultaneous relationships) may help in extracting credible estimates of causal flows and the transient climate response simply by looking at the joint dynamics of two climatic time series.

Estimating Remaining Carbon Budgets Using Temperature Responses Informed by CMIP6

A remaining carbon budget (RCB) estimates how much CO2 we can emit and still reach a specific temperature target. The RCB concept is attractive since it easily communicates to the public and policymakers, but RCBs are also subject to uncertainties. The expected warming levels for a given carbon budget has a wide uncertainty range, which increases with less ambitious targets, i.e., with higher CO2 emissions and temperatures. Leading causes of RCB uncertainty are the future non-CO2 emissions, Earth system feedbacks, and the spread in the climate sensitivity among climate models. The latter is investigated in this paper, using a simple carbon cycle model and emulators of the temperature responses of the Earth System Models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) ensemble. Driving 41 CMIP6 emulators with 127 different emission scenarios for the 21st century, we find almost perfect linear relationship between maximum global surface air temperature and cumulative carbon emissions, allowing unambiguous estimates of RCB for each CMIP6 model. The range of these estimates over the model ensemble is a measure of the uncertainty in the RCB arising from the range in climate sensitivity over this ensemble, and it is suggested that observational constraints imposed on the transient climate response in the model ensemble can reduce uncertainty in RCB estimates.

Energy Budget Constraints on the Time History of Aerosol Forcing and Climate Sensitivity

AbstractAn observationally constrained time series of historical aerosol effective radiative forcing (ERF) from 1750 to 2019 is developed in this study. We find that the time history of aerosol ERFs diagnosed in CMIP6 models exhibits considerable variation and explore how the time history of aerosol forcing influences the probability distributions of present‐day aerosol forcing and emergent metrics such as climate sensitivity. Using a simple energy balance model, trained on CMIP6 climate models and constrained by observed near‐surface warming and ocean heat uptake, we derive estimates for the historical aerosol forcing. We find 2005–2014 mean aerosol ERF to be −1.1 (−1.8 to −0.5) W m−2relative to 1750. Assuming recently published historical emissions from fossil fuel and industrial sectors and biomass burning emissions from SSP2‐4.5, aerosol ERF in 2019 is −0.9 (−1.5 to −0.4) W m−2. There is a modest recovery in aerosol forcing (+0.025 W m−2 decade−1) between 1980 and 2014. This analysis also gives a 5%–95% range of equilibrium climate sensitivity of 1.8°C –5.1°C (best estimate 3.1°C) with a transient climate response of 1.2°C –2.6°C (best estimate 1.8°C).