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 m-2, comprised of 1.81 (±0.09) W m-2 from CO2, 1.08 (± 0.21) W m-2 from other well-mixed greenhouse gases, -1.01 (± 0.23) W m-2 from aerosols and -0.09 (±0.13) W m-2 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 m-2. The majority of the remaining 0.21 W m-2 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 m-2, 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.
Introduction
Effective radiative forcing (ERF) has gained acceptance as the most useful measure of defining the impact on Earth's energy imbalance to a radiative perturbation . These perturbations can be anthropogenic or natural in origin and include changes in greenhouse gas concentrations, aerosol burdens, land use characteristics, solar activity and volcanic eruptions. Since the start of the industrial era until the present day, anthropogenic forcing has typically been increasing and has been the dominant component of the total forcing on the Earth system except for brief periods following large volcanic eruptions . The main constituents of anthropogenic ERF are a positive forcing from greenhouse gases and a partially offsetting negative forcing from aerosols. While greenhouse gas forcing is reasonably well-known, aerosol forcing is more uncertain due to the spatial variation of aerosols, their short atmospheric lifetime and their complex interactions with clouds .
ERF is useful because equilibrium temperatures are more closely related to surface warming in the forcing–feedback relationship of Earth's atmosphere:
ΔN=F-λΔT,
where ΔN, F, λ and ΔT are the top-of-atmosphere (TOA) energy imbalance, (effective) radiative forcing, climate feedback parameter and change in global-mean surface air temperature respectively. showed that using ERF rather than radiative forcing (RF) reduces the need for forcing-specific efficacy values (the temperature response per unit forcing), first introduced by as an observation that different values of λ better predicted ΔT for different forcing agents under RF. Conversely, evaluating ERF is less straightforward than RF, requiring climate model integrations, and numerous different methods of calculating ERF exist with their own benefits and drawbacks .
The difference between ERF and RF is that ERF includes all tropospheric and land surface adjustments, whereas RF only includes the adjustment due to stratospheric temperature change . Adjustments are often termed “rapid” ; however, there is no formal separation of adjustments and climate feedbacks based on timescale alone . It is conceptually more appropriate to divide adjustments as those changes in state that occur purely as a result of the action of a forcing agent from slow feedbacks that occur as a result of a change in global-mean surface temperature. The instantaneous radiative forcing (IRF) is the initial perturbation to Earth's radiation budget and unlike the RF and ERF does not include adjustments. By analysing atmosphere-only climate simulations using fixed climatological sea surface temperatures (SSTs) and sea ice distributions, surface-temperature-driven feedbacks are largely suppressed except for a small contribution from land surface warming or cooling , allowing for adjustments to be diagnosed from atmospheric state changes . This provides insight into the mechanisms contributing to the effective radiative forcing. For example, the ERF of black carbon is half of the impact estimated from its IRF as a consequence of its strong atmospheric absorption and adjustments arising from how it perturbs tropospheric heating rates, affecting the distribution of tropospheric temperatures, water vapour and clouds .
The experimental protocol for determining (effective) radiative forcing in models has been extended since
Phase 5 of the Coupled Model Intercomparison Project (CMIP5). CMIP5 included experiments for present-day (year 2000) all-aerosol and sulfate-only forcing (; CMIP5 experiment labels sstClimAerosol and sstClimSulfate) and 4×CO2 forcing sstClim4xCO2; with respect to a pre-industrial baseline with climatological SSTs and sea ice distributions (sstClim). A handful of IRF outputs from quadrupled CO2 experiments were also obtained. For CMIP6, the Radiative Forcing Model Intercomparison Project RFMIP; provides a number of present-day time-slice and historical-to-future transient experiments designed to evaluate the ERF in climate models for different forcing agents, providing insight into why climate models respond the way they do to particular forcings. This is important when diagnosing climate feedbacks , given the role of forcing in Earth's energy budget as in Eq. (), and knowledge of forcing is required for attribution of historical temperature change and evaluating non-CO2 contributions to remaining carbon budgets and in future scenario projections . Effective radiative forcings derived from models can be used to validate assumptions derived from other lines of evidence, particularly for aerosol forcing, as is done by the Intergovernmental Panel on Climate Change (IPCC) in their periodic Assessment Reports.
Models and experimental protocol
We use results from 17 state-of-the-art atmospheric general circulation models (GCMs) and Earth system models (ESMs) contributing to Tier 1 of RFMIP (Table 1) as part of CMIP6 . In addition, GISS-E2-1-G provided two physics variants, r1i1p1f1 and r1i1p3f1 (hereafter “p1” and “p3”), with aerosol treatments that are different enough to justify treating the variants as separate models, bringing the total to 18. Models with diagnostics available on the Earth System Grid Federation (ESGF) up until 13 May 2020 have been analysed. Each model is run in atmosphere-only mode using pre-industrial climatologies of sea surface temperatures (SSTs) and sea ice distributions from at least 30 years of the same model's corresponding coupled pre-industrial control run (piControl; ). RFMIP's Tier 1 calls for 30-year time-slice experiments forced with 4× pre-industrial CO2 concentrations (RFMIP name piClim-4xCO2), all present-day anthropogenic forcers (piClim-anthro), present-day well-mixed greenhouse gases (piClim-ghg), present-day aerosols (piClim-aer) and present-day land use (piClim-lu) in this fixed-SST configuration.
All forcing components that are not perturbed in a particular experiment remain at pre-industrial (year 1850) values, and “present day” is defined as year 2014 conditions. A 30-year experiment with pre-industrial conditions, piClim-control, is also performed as a reference case, and all results presented in this paper are with reference to piClim-control, accounting for the possibility that models may have a non-zero pre-industrial TOA flux imbalance. Results from the 4×CO2 experiment are also rescaled to the ratio of 2014 to 1850 CO2 concentrations of approximately 1.4× pre-industrial by a factor of 0.2266, being the ratio of RF from 1.4×CO2 to 4×CO2 from the formula. This is performed to isolate an estimate of the CO2-only contribution to the present-day forcing and is based on year 1850 and year 2014 CO2 concentrations of 284.32 and 397.55 ppm respectively along with the 1850 concentrations of 808.25 ppb for CH4 and 273.02 ppb for N2O. Except where explicitly stated, we present results from this experiment as 1.4×CO2.
Contributing climate models to RFMIP-ERF Tier 1. The adjustment time is based on approximately how long stratospheric temperatures take to equilibriate in the 4×CO2 experiment (Fig. ). ISCCP (International Satellite Cloud Climatology Project) simulator diagnostics are indicated where existent.
ModelAtmospheric resolution (long × lat)AdjustmentModelISCCPReferencetimescaleyearssimulator(years)ACCESS-CM21.875∘×1.25∘, 85 levels to 85 km130CanESM52.81∘×2.81∘, 49 levels to 1 hPa150allCESM21.25∘×0.9∘, 32 levels to 2.25 hPa130allCNRM-CM6-11.4∘×1.4∘, 91 levels to 0.01 hPa530CO2,ghg, aer, anthroCNRM-ESM2-11.4∘×1.4∘, 91 levels to 0.01 hPa1530allEC-Earth30.7∘×0.7∘, 91 levels to 0.01 hPa130GFDL-CM41.25∘×1∘, 33 levels to 1 hPa130allGFDL-ESM41.25∘×1∘, 49 levels to 1 hPa130GISS-E2-1-Ga2.5∘×2∘, 40 levels to 0.1 hPa531/41bHadGEM3-GC31-LL1.875∘×1.25∘, 85 levels to 85 km130allIPSL-CM6A-LR2.5∘×1.27∘, 79 levels to 80 km1030allMIROC61.4∘×1.4∘, 81 levels up to 0.004 hPa130aerMPI-ESM1-2-LR1.875∘×1.875∘, 47 levels up to 0.01 hPa131MRI-ESM2-01.125∘×1.125∘, 80 levels to 0.01 hPa130allNorESM2-LM2.5∘×1.875∘, 32 levels to 3 hPa130NorESM2-MM1.25∘×0.9375∘, 32 levels to 3 hPa130UKESM1-0-LL1.875∘×1.25∘, 85 levels to 85 km345all
a GISS-E2-1-G produced two physics variants for piClim-control and piClim-aer; physics_version=1 (p1) includes aerosol and ozone specified by pre-computed transient fields , and physics_version = 3 (p3) includes aerosol–cloud interactions. Both physics versions are analysed in this paper and treated as separate models.
b 41 years for r1i1p3f1.
The experiments and results presented in this study follow on from the assessment of ERF and adjustments in 11 models contributing to the Precipitation Driver and Response Model Intercomparison Project PDRMIP; see in . In idealised experiments of 2×CO2 concentrations, 3×CH4 concentrations, 10× black carbon (BC) emissions or burdens, 5×SO4 emissions or burdens and a 2 % solar constant increase were analysed from CMIP5-era and interim models. Only the 4×CO2 experiment has a similar experiment for comparison in , whereas the RFMIP protocol focuses more on combinations of anthropogenic forcers. In addition, extended model diagnostics allow us to determine cloud responses and aerosol forcing in more detail in this study.
Effective radiative forcing
Using climatological SSTs allows for ERF to be diagnosed as the difference of top-of-atmosphere net radiative flux between a given forcing experiment and a pre-industrial control simulation . Using 30-year time slices generally results in standard absolute errors of less than 0.1 W m-2. Although inter-annual variability affects the diagnosed ERF using this climatological-SST method, the standard error in the estimates obtained is much smaller than using a fully coupled ocean–atmosphere model with a Gregory regression , and as such fewer model years are needed to diagnose ERF. Two advantages of this is that it reduces the computational burden for modelling centres and can also be used to diagnose forcings of the order of 0.1 W m-2. For this reason, the climatological-SST method is implemented to derive forcing in RFMIP, and ERF in this paper (without qualifier) is taken to mean this.
The climatological-SST method of deriving ERF includes the TOA flux changes resulting from land surface warming or cooling as part of the ERF. Conceptually, any land surface temperature change as a response to forcing should be excluded in the same way that SST changes are , but prescribing land surface temperatures is difficult in GCMs, and this has not been performed in RFMIP. In essence, the goal is to completely isolate the forcing from any surface temperature change (ΔT) or feedbacks (λ) in Eq. (). We test several methods to correct for adjustments to attempt to isolate forcing at ΔT=0 (also performed in ):
Effective radiative forcing (ERF). ERF is reserved to mean the TOA flux difference between a perturbed and control simulation, with climatological SSTs and sea ice distributions and no correction for land surface temperature change, as in , , , and .
Effective radiative forcing using a Gregory regression (ERF_reg). ERF_reg is calculated from each model's CMIP abrupt-4xCO2 experiment by regressing the annual temperature anomaly compared to the same model's pre-industrial control (piControl) against the annual TOA energy imbalance anomaly ΔN in Eq. () and finding the intercept at ΔT=0, as in . This is done for the first 20 years of model output to avoid the changing value of λ over time present in many models ; using the full 150 years tends to underestimate the forcing (denoted ERF_reg150; Table S1 in the Supplement). It is only possible to determine ERF_reg for 4×CO2, as coupled abrupt forcing experiments are not performed for other forcing agents as part of CMIP6.
Stratospherically adjusted radiative forcing (RF). All tropospheric and surface adjustments, calculated using radiative kernels (Sect. ), are subtracted from the ERF, leaving just the stratospheric temperature adjustment to the IRF. The RF is included for historical comparison, although it is usually calculated using an offline method such as fixed dynamical heating . It should be noted that the stratospheric adjustment is included in all definitions of ERF.
Land-surface-corrected effective radiative forcing (ERF_ts). Land surface temperature change adjustment is subtracted from the climatological-SST ERF using the surface temperature radiative kernel.
Tropospherically corrected effective radiative forcing (ERF_trop). In addition to land surface warming a proportion of tropospheric temperature and water vapour change is subtracted from the ERF using radiative kernels, by assuming a fixed lapse rate in the troposphere based on the land surface temperature change. The remaining tropospheric temperature change when the constant lapse rate is subtracted is treated as the tropospheric temperature adjustment. The water vapour correction from the land surface warming is taken as the fraction of the adjustment from the constant lapse rate to the total tropospheric temperature adjustment. The surface albedo change is also removed, whereas no cloud adjustment is included justified by cloud adjustments to a large extent depending on heating or cooling in the troposphere . This was known as ERF_kernel in .
Feedback-corrected effective radiative forcing (ERF_λ). An amount corresponding to the global-average near-surface air temperature (GSAT) warming multiplied by the model's climate feedback parameter from its corresponding CMIP abrupt-4xCO2 run is subtracted from the fixed-SST ERF. The same value of λ from abrupt-4xCO2 is applied to the GSAT change in all experiments. This method was first investigated by and is known as ERF_fSST_ΔTland in .
Effective radiative forcing from each Tier 1 time-slice RFMIP experiment for each model (W m-2). Also shown is the 4×CO2 ERF scaled to 2014 concentrations (as 1.4×CO2) and the residual forcing (anthropogenic, WMGHGs, aerosol and land use). WMGHGs: well-mixed greenhouse gases. Note that not all models performed all experiments.
Comparison of radiative forcing (RF, which by definition includes stratospheric temperature adjustment), effective radiative forcing with tropospheric correction (ERF_trop), effective radiative forcing with land surface kernel correction (ERF_ts), feedback-corrected ERF (ERF_λ) and fixed-SST ERF. For CO2 forcing, ERF from a Gregory regression (ERF_reg) from each model's corresponding abrupt-4xCO2 CMIP simulation is also given. The ERF is compared with the IRF and adjustments (RA) for each of the present-day RFMIP-ERF time-slice experiments (1.4×CO2 is shown instead of 4×CO2 for better comparison with other forcing agents). Individual models are numbered.
Table shows the ERF diagnosed from each forcing and each model using the climatological-SST method, and Fig. shows the ERF, diagnosed IRF and adjustments from each RFMIP Tier 1 experiment. Values for the different methods for calculating forcing are given in Tables S1–S5. Instantaneous forcing (IRF) is calculated as the difference of the ERF and the sum of adjustments, with an exception being land use forcing where IRF is calculated directly from the surface albedo kernel. In keeping with the definitions of ERF and adjustments, IRF is defined at the TOA in this study. Adjustment calculations are explained in detail in Sect. .
For ease of comparison we show 1.4×CO2 instead of 4×CO2, with the scaling to present-day concentrations assumed to apply to ERF, IRF and all adjustments proportionally. Figure also shows the ERF_reg (for 4×CO2), ERF_ts, ERF_λ and RF. In general, the methods that correct for land surface temperature change (ERF_ts, ERF_trop and ERF_λ) result in forcings that are slightly stronger than non-corrected ERF, with ERF_trop giving the largest forcing change (8.84 W m-2 for 4×CO2 in ERF_trop compared to 7.98 W m-2 for ERF). For CO2, ERF_reg results in a similar mean estimate of ERF to the fixed-SST method. Excluding CNRM-ESM2-1 for reasons described in the next section, the 4×CO2 ERF_reg is 8.09 W m-2 compared to 7.99 W m-2 for ERF.
Forcing adjustmentsNon-cloud adjustments
Adjustments to the radiative forcing describe flux changes resulting from changing atmospheric or surface state, in response to a forcing, but unrelated to the change in globally averaged surface temperature (thus decoupling them from climate feedbacks; ). Adjustments to non-cloud changes in this study are calculated using radiative kernels . The difference in an atmospheric state variable x (air temperature, surface temperature, specific humidity or surface albedo) between a forcing perturbation (pert) and piClim-control (base) is multiplied by the kernel Kx to derive the adjustment Ax:
Ax=Kx(xpert-xbase).
The radiative kernel describes the change in TOA fluxes for a unit change in state for x∈{T,Ts,q,α} where T is atmospheric air temperature, Ts is surface temperature, q is water vapour and α is surface albedo. KT and Kq are four dimensional (month, pressure level, latitude and longitude), and KTs and Kα are three dimensional (month, latitude and longitude). Kernels are produced for both long-wave and short-wave radiation changes. Typical unit changes are 1 K for temperature, the change in specific humidity that maintains constant relative humidity for a temperature increase of 1 K for water vapour and 1 % additive for surface albedo. For the division of temperature into stratospheric and tropospheric components, the WMO (World Meteorological Organization) definition of the lapse-rate tropopause is used from each model's piClim-control run, using geopotential height as an approximation of geometric height on model pressure levels.
The water vapour kernel describes the change in TOA flux for a perturbation that maintains relative humidity for a temperature increase of 1 K, the effect being that specific humidity increases. The assumption therefore is that relative humidity is approximately constant between perturbation and control runs, which is found to be true in coupled experiments . Note that the difference in states is taken for the logarithm of water vapour concentration in Eq. (). More details on the application of the kernel method can be found in Supplementary Material.
In this paper we use radiative kernels derived from the atmospheric component of the HadGEM3-GC31-LL model (HadGEM3-GA7.1), interpolated to the 19 standard CMIP6 pressure levels . With the exception of stratospheric temperature adjustments to greenhouse gas forcing, structural differences introduced by using different kernels are well within 0.1 W m-2, and the HadGEM3-GA7.1 kernel is representative of the population of radiative kernels commonly used in the literature for tropospheric and surface adjustments (Fig. S1 in the Supplement); we use this particular kernel for its improved stratospheric resolution as outlined in .
Stratospheric adjustments to greenhouse-gas-driven experiments are expected to equilibriate within a few model months . We find that the time to reach equilibrium varies between models for a 4×CO2 forcing. Figure shows the time taken for the stratospheric temperature adjustment, and hence stratospheric temperatures, to adjust to a 4×CO2 forcing. In CNRM-ESM2-1, concentrations of CO2 are relaxed towards the 4× pre-industrial level below 560 hPa and allowed to propagate throughout the atmosphere, therefore taking around 15 years to reach an approximate uniform atmospheric concentration. A similar specification is implemented in the abrupt-4xCO2 run of CNRM-ESM2-1, causing ERF_reg to be biased low (Fig. ). This highlights one advantage of the fixed-SST based methods over the Gregory regression, as these “spin-up” years can simply be discarded with a fixed-SST measure of ERF. The CO2 treatment in CNRM-ESM2-1 is in contrast to the physical climate model from the same group (CNRM-CM6-1). However, even in some physical models, we find that the time to reach equilibrium varies between models and may be up to 10 years (e.g. in IPSL-CM6A-LR; Fig. ). For this reason, we discard the first few years of model output where the stratosphere is still adjusting to a forcing for the 4×CO2, well-mixed greenhouse gas (WMGHG) and anthropogenic forcing experiments (Table 1). We find this issue is not present in the aerosol or land use experiments. It is important to emphasise that our stratospheric adjustment is calculated in a different way to the usual RF method, which uses an offline radiative transfer method. It may therefore be the case that differences are due to a change in tropopause height in greenhouse-gas-driven experiments .
Transient response of the stratospheric temperature adjustment to a 4×CO2 forcing. The small spike in year 6 in CanESM5 is due to an unseasonably low tropical tropopause in July of year 6, resulting in much of the temperature adjustment at the 100 hPa level to be counted in the stratosphere.
Cloud adjustments
The radiative effect of clouds depends on their coverage (both within layer and total), ice water content, liquid water content, droplet effective radius and ice particle habit. Cloud properties vary extensively from model to model, and unlike pressure level diagnostics of temperature and humidity, cloud diagnostics are not output on 19 standard pressure levels in CMIP. A number of different approaches have therefore been used to estimate cloud adjustments, depending on availability of diagnostics and model-specific setup, and we can exploit methods originally designed for cloud feedback calculations for calculating adjustments. Where cloud adjustments can be calculated with more than one method, we take the mean of each available method. In some models and experiments, cloud adjustments cannot be calculated, and no estimate is made.
ISCCP simulator kernel
The ISCCP simulator provides a joint 7×7 histogram of cloud visible-wavelength optical depth (τ) and cloud top pressure (CTP). These outputs can be multiplied by the ISCCP simulator kernel to estimate the impact of cloud changes on top-of-atmosphere fluxes. Ten models included ISCCP simulator diagnostics within their RFMIP output (Table 1).
The ISCCP simulator kernel reports all flux changes resulting from clouds. For CO2, WMGHG and land use forcings, it is assumed that cloud droplet effective radius does not change (except for the land use experiment in NorESM2-LM as discussed in Sect. , but this model did not include ISCCP simulator diagnostics), and therefore in these experiments the SW (short-wave) flux changes from the ISCCP simulator kernel are treated as the cloud adjustment. For aerosol and total anthropogenic forcing this is usually not the case, as most models include aerosol–radiation interactions (significant in the SW), with ice particle behaviour also changing in the MRI-ESM2-0, MIROC6 and CESM2 models, which affects LW (long-wave) fluxes. NorESM2-LM also includes the effects of mineral dust and BC on heterogeneous ice nucleation . Following we treat the cloud-albedo response to aerosols as part of the IRF, and the ISCCP simulator kernel is unable to separate this effect from any adjustment. We assume that any LW effect from aerosol–cloud interactions is small except in those models that include aerosol effects on ice clouds.
Approximate partial radiative perturbation with liquid water path adjustment
The approximate partial radiative perturbation (APRP; Sect. ) method uses standard climate model diagnostics to estimate the components of SW ERF attributed to cloud fraction change and all-sky and clear-sky scattering and absorption. With no changes in aerosol forcing, the changes in cloud absorption, cloud scattering and cloud amount calculated from APRP can be taken to be the SW cloud adjustment. We use this estimate for CO2, WMGHG and land use forcing.
For aerosol forcing, the effect of cloud amount changes calculated by APRP (ACLT) is an adjustment, but the cloud scattering is a combination of radiative forcing due to aerosol–cloud interactions (RFaci), treated as part of the IRF, and adjustments due to cloud liquid water path (LWP) changes ALWP;. For the LWP adjustment we use a relationship obtained in in which LWP adjustment (W m-2) scales linearly with the vertically integrated in-cloud liquid water path (kg m-2):
ALWP=-100037.6clwvipert-clivipertcltpert/100-clwvibase-clivibasecltbase/100,
where clwvi, clivi and clt are the CMIP6 variable labels for total cloud water path, ice water path and total cloud fraction in percent. We then isolate the RFaci as
RFaci=ERFaci-ALWP-ACLT,
with ERFaci, the effective radiative forcing due to aerosol–cloud interactions, calculated from APRP (Sect. ).
For anthropogenic total forcing, the RFaci calculated in Eq. () from the aerosol forcing experiment is subtracted from the total derived cloud change under APRP, which includes contributions from greenhouse gases and land use as well as RFaci. For models not including ice cloud nucleation, the LW cloud adjustment for aerosols is estimated from the change in cloud radiative effect (CRE; difference between all-sky and clear-sky fluxes). For other experiments this results in a biased estimate of cloud adjustment due to masking of LW adjustments.
A direct estimate of cloud radiative effect can be obtained by substituting model cloud fields into an offline radiative transfer model. We perform these offline calculations using the SOCRATES (Suite of Community Radiative Transfer codes based on Edwards and Slingo) radiative transfer code . This is produced by substituting fields of three-dimensional cloud fraction, cloud water content and cloud ice content from each model and experiment into a climatology for the year 2014 provided by ERA5 (ECMWF Reanalysis) . Taking the cloud fields in each experiment minus those from the control gives ALWP+ACLT in each model. As only monthly-mean diagnostics are available from models in general, we only attempt this in the LW, which is assumed to be less biased than the SW . The monthly-mean cloud fraction, ice water content and liquid water content variables in all experiments are scaled by a model-dependent factor that ranges between 0.68 and 1.5 to ensure that TOA LW outgoing flux is approximately 240.2 W m-2 in the control experiment, in line with TOA observations .
Kernel masking
In the land use experiment, IRF is directly estimated from the surface albedo kernel such that IRF = Aα. As there are no other unknowns in the kernel decomposition, cloud adjustments can be calculated using the difference between all-sky and clear-sky fluxes such that
Ac=(ERF-ERFclr)-(Aα-Aαclr)-∑i∈{T,Ts,q}(Ai-Aiclr),
where the clr superscript refers to fluxes calculated with clear-sky radiative kernels.
Multi-model results
Figure shows the contribution to the total adjustment in each experiment from land surface temperature, tropospheric temperature, stratospheric temperature, water vapour, surface albedo and clouds. No corrections for tropospheric or land surface warming as discussed in Sect. have been performed for these results.
Adjustments broken down by mechanism in each of the present-day RFMIP-ERF time-slice experiments. Black and grey numbers indicate individual models; coloured bars indicate the multi-model mean.
Figure shows the effect on TOA radiative flux arising from cloud responses from the ISCCP simulator for each experiment from models that provided these diagnostics (Table 1). In this figure, histogram boxes not marked with a cross are where 75 % or more of the models agree on the sign of the cloud fraction or radiative flux change, following .
Global-mean change in ISCCP-simulated cloud fraction (clisccp) in CTP-τ space (first column) and consequential changes in SW (second column), LW (third column) and net (fourth column) radiation when convoluted with the ISCCP cloud kernel. Grey crosses show where less than 75 % of models agree on the sign. Figure shows the multi-model mean cloud fraction and radiative effect.
For 1.4×CO2 the change in cloud fraction, as well as the radiative fluxes, is scaled down from the 4×CO2 experiment using .
Carbon dioxide
The multi-model mean ERF from a quadrupling of CO2 is 7.98 W m-2 (±0.38 W m-2; all ranges given as 1 standard deviation). A point of comparison for ERF is , who computed estimates of and parametric fits for radiative forcing accounting for masking by clouds and stratospheric temperature equilibration, using a tropical and mid-latitude profile to represent the global mean. The implied ERF from RFMIP models for 2×CO2 is 3.81 (±0.18) W m-2 when scaling down the 4×CO2 results using the formula, comparable to a radiative forcing of 3.80 W m-2 for a doubling of CO2 in . Both estimates are slightly higher than the best estimate of 3.71 W m-2 from the IPCC's Fifth Assessment Report (AR5; ). Scaling down the 4×CO2 forcing using , our derived multi-model mean for 1.4×CO2 is 1.81 (±0.09) W m-2. As shown in Fig. and discussed in Sect. , ERF is approximately equal to RF for CO2, and we apply the Etminan formula to ERF.
The 4×CO2 ERF from 17 CMIP6 models is larger, but not significantly so (p value of 0.13 using Welch's t test), than the 4×CO2 ERF from 13 CMIP5 models of 7.53 (±0.89) W m-2. In addition, CMIP6 models are notable for their smaller spread in CO2 ERF than CMIP5 models (Fig. ). show that ERF_reg150 for 4×CO2 also increases in CMIP6 compared to CMIP5 and attribute 20 % of the increase in multi-model mean effective climate sensitivity (ECS) in CMIP6 to this.
We note that a long-standing problem in GCMs has been on the diversity in the forcing of CO2, which may result both from model broadband radiation parameterisation error in the IRF component and differences in base state climatology between models. The reduction in spread of CO2 forcing in CMIP6 may be indicative that model radiation parameterisations are improving, for example as documented in HadGEM3-GC31-LL and UKESM1-0-LL , but it could also be from a convergence in model base states, including clouds.
Aerosol and 4×CO2 effective radiative forcing from CMIP5 sstClim4xCO2 and sstClimAerosol experiments and CMIP6 RFMIP experiments. Numbers at the bottom of each plot give the number of participating models.
The breakdown of ERF into adjustments is shown in Table with the corresponding 4×CO2 values in Table S6. Stratospheric temperature adjustment dominates for CO2-driven simulations, which is well-known . Tropospheric adjustments approximately sum to zero such that the overall adjustment approximately equals the stratospheric adjustment, and RF is a good approximation to ERF . Nevertheless, individual tropospheric adjustments are non-zero and significant. A warming land surface and troposphere leads to a negative adjustment (more outgoing LW radiation to space) that is partially offset by increased tropospheric water vapour (analogous to the water vapour feedback). Cloud adjustments are overall positive, dominated by a reduction in mid-troposphere clouds driven by tropospheric warming, leading to a positive SW radiative effect (Fig. ). The LW effect is small in comparison so that the SW effect dominates the net cloud adjustment.
ERF, IRF and adjustments (W m-2) by component from 1.4×CO2, scaled down from the 4×CO2 RFMIP experiment. Adj.: adjustment; ts: surface temperature; ta_tr: tropospheric temperature; ta_st: stratospheric temperature; hus: water vapour; albedo: surface albedo; cl: clouds.
The spatial pattern of adjustments is shown in Fig. . In Figs. 6 to 8, 12 and 13, cloud changes are only shown from the ISCCP simulator kernels in panels (g–i) and are not the means of all participating models, whereas ERF and non-cloud adjustments in panels (a–f) are multi-model means. Hatched areas are defined where less than 75 % of models agree on the sign of the change. Stratospheric cooling is spatially uniform and results in a positive adjustment of +0.61 W m-2, i.e. around one-third of the total ERF. Tropospheric temperature adjustments are globally negative and robust. Cloud changes show several robust spatial patterns, including positive changes over land in Eurasia and North America.
Multi-model mean spatial patterns of (a) effective radiative forcing, (b–f) adjustments and (g–i) cloud contributions to ERF for 1.4×CO2. Hatched regions are where less than 75 % of models agree on the sign of the change.
As for Table but for 1850–2014 well-mixed greenhouse gas forcing.
The ERF from all well-mixed greenhouse gases is evaluated to be 2.89 (±0.19 W m-2) for 1850–2014, implying a contribution of 1.08 (±0.21) W m-2 from non-CO2 WMGHGs (uncertainties in quadrature, and this definition excludes changes in ozone). Tier 1 of RFMIP does not contain additional granularity to break down non-CO2 forcing by species; however dedicated experiments to derive ERF from methane, nitrous oxide and halocarbons separately are part of the protocol for the Aerosol and Chemistry Model Intercomparison Project AerChemMIP;.
There is also a substantial adjustment arising from WMGHG forcing, and again this is mostly driven by stratospheric cooling implied by the observation that ERF and RF are approximately equal. This confirms PDRMIP (Precipitation Driver and Response Model Intercomparison Project) model behaviour for CO2 and CH4 forcing , which found that tropospheric and land adjustments, while individually significant, approximately sum to zero, leaving just the stratospheric temperature adjustment. Unlike in , who found that the stratospheric temperature adjustment to methane was approximately zero, we find a larger stratospheric temperature adjustment for WMGHGs compared to CO2 implying a positive non-CO2 WMGHG stratospheric adjustment, although this cannot be attributed to individual gases.
The multi-model mean non-CO2 WMGHG ERF of 1.08 W m-2 is close to the 1850–2014 RF of 1.09 W m-2 made up of CH4 (0.55 W m-2) plus N2O (0.17 W m-2) from , plus halocarbons (0.37 W m-2) using relationships from .
As for CO2-only forcing, the total adjustment approximately equals the stratospheric temperature adjustment, implying that tropospheric and surface adjustments approximately cancel (Table ) so that the spread in their sum is smaller than for each component individually. For the ISCCP simulator cloud adjustments, a similar pattern can be seen from all WMGHGs to CO2-only forcing, with a larger reduction in mid-troposphere cloud fraction leading to a greater positive SW adjustment which dominates the net adjustment.
The spread in ERF and stratospheric temperature adjustments is larger for WMGHG than for CO2 forcing alone. One factor may be the inclusion or exclusion of stratospheric chemistry, which affects ozone formation. The effect can be seen by comparing Earth system (ESM) and physical models from the same group: the UKESM1-0-LL ESM (model 18) to the HadGEM3-GC31-LL physical model (model 11) and CNRM-ESM2-1 (model 5) to CNRM-CM6-1 (model 4). The physical models show ERFs around 0.2 W m-2 greater than the ESMs, a greater IRF and a smaller stratospheric temperature adjustment. Additionally, for UKESM1-0-LL, large and compensating ERFs from CH4 (+0.93 W m-2) and halocarbons (-0.33 W m-2), resulting from interactive chemistry, bring the total WMGHG ERF closer to the no-chemistry ERFs total from HadGEM3-GC31-LL .
The spatial patterns are overall similar to the CO2 experiment (Fig. ) with a larger magnitude.
As Fig. but for present-day WMGHG forcing.
As for Table but for 1850–2014 aerosol forcing.
No.ModelERFIRFAdj.tsta_trta_sthusalbedocl1ACCESS-CM2-1.090.070.110.01-0.00-0.032CanESM5-0.85-0.51-0.340.02-0.16-0.100.180.01-0.283CESM2-1.37-1.430.06-0.02-0.00-0.070.10-0.010.054CNRM-CM6-1-1.15-1.190.040.070.15-0.00-0.08-0.06-0.035CNRM-ESM2-1-0.74-0.750.010.060.10-0.01-0.06-0.02-0.066EC-Earth3-0.80-0.66-0.140.06-0.020.01-0.02-0.04-0.147GFDL-CM4-0.73-0.56-0.170.05-0.03-0.020.07-0.07-0.178GFDL-ESM4-0.70-0.37-0.330.04-0.15-0.010.10-0.05-0.269GISS-E2-1-G p1-1.32-0.46-0.860.060.22-0.03-0.09-0.09-0.9310GISS-E2-1-G p3-0.93-1.000.070.050.13-0.03-0.06-0.00-0.0111HadGEM3-GC31-LL-1.10-1.04-0.060.050.100.01-0.05-0.04-0.1412IPSL-CM6A-LR-0.63-0.60-0.030.050.06-0.110.01-0.00-0.0313MIROC6-1.04-1.130.100.060.170.01-0.10-0.03-0.0215MRI-ESM2-0-1.21-0.46-0.740.04-0.24-0.000.17-0.02-0.6816NorESM2-LM-1.21-1.09-0.110.000.03-0.050.02-0.04-0.0817NorESM2-MM-1.26-1.10-0.160.030.01-0.020.05-0.04-0.1918UKESM1-0-LL-1.11-0.97-0.140.060.010.010.01-0.05-0.18Mean-1.01-0.83-0.180.050.03-0.020.01-0.03-0.20SD0.230.310.270.020.120.040.080.030.25AerosolsForcing and adjustments
Present-day aerosol ERF is -1.01 (±0.23) W m-2 from 17 models. The full range of aerosol ERF estimates for 2014 versus 1850 is -0.63 to -1.37 W m-2. This is a narrower range of ERF than similar experiments performed with CMIP5 models for year 1850 and year 2000 forcings , particularly in relation to the lower (more negative) bound of aerosol forcing. Based on the 2000–1850 estimate of -1.17 (±0.30) W m-2 from , aerosol forcing in CMIP6 models is less negative than in CMIP5, but this difference again is not significant (p value of 0.15). Some of this multi-model mean difference is likely due to lower emissions of aerosol precursors in 2014 relative to 2000 along with updated historical estimates for CMIP6 , although it is not clear that this explains the reduction in model spread in CMIP6. It should also be borne in mind that our range does not include the E3SM (Energy Exascale Earth System Model) model, which diagnosed aerosol forcing to be -1.65 W m-2 for 2005–2014 from a pair of parallel all-forcing and pre-industrial aerosol forcing atmosphere-only runs Fig. 25 in. This highlights the likelihood that the inclusion of more models submitting results to RFMIP would extend the CMIP6 range of aerosol forcing, but the same may also have been true in CMIP5, where only a subset of models performed the sstClimAerosol experiment.
Atmospheric adjustments are small in magnitude in the aerosol forcing experiment but large enough such that there is a noticeable difference between ERF and RF (Fig. ; Table S3). The small non-cloud adjustments in most models show that the aerosol forcing is dominated by scattering aerosols (sulfate, organics and, for a limited number of models, nitrates) rather than black carbon . Additionally, in two of the four models that provide the single-forcing BC experiment in AerChemMIP (Aerosol Chemistry Model Intercomparison Project; CNRM-ESM2-1 and UKESM1-0-LL), the overall adjustment is small , in contrast to findings in PDRMIP models . In MRI-ESM2-0 (model 15) there are strong tropospheric temperature and cloud changes to black carbon forcing resulting in a negative adjustment overall .
For aerosol forcing, the aerosol–cloud interactions dominate, with an increase in cloud optical depth at all cloud heights. As cloud droplet effective radius decreases, cloud albedo, and hence optical depth, increases. This also implies that absorbing aerosols play only a minor role in most models, as BC induces strong adjustments that cause a general increase in cloud height in PDRMIP models from an increasing tropospheric stability (; Fig. S2). There is no evidence of this in the RFMIP aerosol forcing experiment, although some models do also include aerosol–cloud interactions from BC, and the effect may be due to the BC forcing being a smaller fraction of the total aerosol forcing than sulfate . Figure S2 shows ISCCP simulator results for the five PDRMIP experiments from the CMIP5-era HadGEM2-ES model, where it can be seen that the aerosol forcing experiment is qualitatively more similar to the 5×SO4 forcing experiment than the 10×BC experiment in PDRMIP. The increase in cloud albedo leads to a strong negative SW radiative effect that is partially compensated by LW effects (note that the ISCCP simulator kernel does not distinguish RFaci from adjustments).
Unlike for WMGHGs, aerosol forcing adjustments are dominated by cloud effects with only small non-cloud components (Table ). For aerosol forcing, all model years are used, as the stratospheric temperature adjustment is negligible. The spread in values of cloud adjustments is large and spans positive and negative values. This reconfirms that atmospheric processes in response to aerosol forcing remains one of the largest uncertainties in climate models. There is also a spread in tropospheric temperature and water vapour adjustments with multi-model means near zero, suggesting that some models respond to aerosols with substantial atmospheric warming or cooling.
For many regions, particularly southern Asia and the eastern Pacific, the aerosol ERF is driven by large and negative cloud changes (Fig. ). The small adjustment overall and increase in cloud optical depth for all ISCCP cloud categories suggest this is driven by an increase in cloud condensation nuclei leading to a more negative RFaci. There are some regions such as the Sahara in which a positive ERF arises and is not easily explained by any adjustment component. This may be a reduction in mineral dust loading and increase in BC loading, leading to a positive forcing (e.g. as seen in NorESM2-LM, Fig. S3).
As Fig. but for present-day aerosol forcing.
The total derived cloud adjustment for aerosols is -0.20 W m-2, derived of -0.04 W m-2 from SW cloud liquid water path adjustment, -0.13 W m-2 from SW cloud fraction change and -0.03 W m-2 from cloud changes in the LW (Table S7).
Relationship to climate sensitivity
The increase in the upper bound, and in the overall spread, of ECS in the CMIP6 model population compared to CMIP5 is well-documented . Figure shows the relationships between ECS and transient climate response (TCR) and aerosol ERF in CMIP6, taking ECS and TCR from each model's abrupt-4xCO2 and 1pctCO2 CMIP runs respectively. There are weak and non-significant positive correlations between ECS and aerosol forcing (r=0.12) and between TCR and aerosol forcing (r=0.26). This suggests that, as a population, models with high sensitivity are not tuning present-day aerosol forcing to be strong in order to reproduce observed warming
MPI-ESM1-2 is the only documented exception. MIROC6 did tune the aerosol forcing to better correspond to the AR5 best estimate but explicitly did not tune for surface temperature.
: it would be expected that these correlations would be negative if this was the case . In CMIP5 models, aerosol forcing was stronger in models with higher ECS and TCR, but not significantly so , although significance emerges if one considers only models which include an aerosol indirect effect . In CMIP3 there was a strong and significant negative correlation between climate sensitivity and aerosol forcing . It may be the case that aerosol forcing over some historical periods is stronger in CMIP6 than in CMIP5, as despite higher climate sensitivity, CMIP6 models warm less than CMIP5 models and observations up until 2000 .
Relationship between (a) ECS and (b) TCR and aerosol ERF in the CMIP6 model ensemble. MPI-ESM1.2-LR (model 14) did not produce the piClim-aer experiment.
Decomposition of aerosol forcing into aerosol–radiation and aerosol–cloud effects
The approximate partial radiative perturbation (APRP) method can be used to decompose short-wave (SW) aerosol forcing into aerosol–radiation interactions (ERFari), aerosol–cloud interactions (ERFaci) and the surface albedo adjustment .
In Sect. we compare other methods to estimate ERFari and ERFaci. ERFari is the component of aerosol forcing that arises from the direct radiative effect of aerosol absorption and scattering (RFari) plus any adjustments (formerly known as the semi-direct effect) arising from perturbations in tropospheric heating rates, humidity and their consequential effects on where clouds form . These adjustments tend to be strong for black carbon but weak for scattering aerosol . ERFaci is composed of any changes in cloud albedo resulting from aerosols acting as cloud condensation nuclei and changing cloud droplet effective radius (RFaci, formerly the first indirect or Twomey effect; ) plus adjustments relating to cloud lifetime and precipitation efficiency that changes liquid water path and cloud fraction (formerly second indirect or Albrecht effect; ). RFaci tends to be strong for sulfate aerosol, but several models also include cloud interactions to other aerosol species, and four models (CESM2, MIROC6, MRI-ESM2-0 and NorESM2-LM) include aerosol interaction on ice clouds. The direct plus Twomey effects (RFari+aci) are treated as the IRF component of aerosol forcing, with the remaining components of ERFari+aci as adjustments .
There is no equivalent long-wave (LW) method to APRP, so we take the approach of and use the cloud radiative effect to decompose LW ERF into ERFari and ERFaci. The advantages of these techniques are that they only require standard CMIP output, and all participating models can provide estimates. Results are displayed in Table and shown in Fig. . In Table S8 the equivalent SW ERFari for clear-sky conditions are shown.
Components of the aerosol forcing diagnosed from the approximate partial radiative perturbation (for SW aerosol components) and from the cloud radiative effect (for LW components). Black diamonds represent multi-model means; black bars show 1 standard deviation.
Contribution of the components of effective radiative forcing from the present-day aerosol time-slice RFMIP experiment. ARI: aerosol–radiation interaction; ACI: aerosol–cloud interaction; scat: scattering; abs: absorption; amt: cloud amount; sum: summation.
SW LW Net ARI ACI ModelscatabssumscatabsamtsumARI + ACIARIACIARI + ACIARIACIARI + ACIACCESS-CM2-0.790.30-0.48-0.80-0.01-0.12-0.93-1.420.260.040.30-0.22-0.89-1.11CanESM5-0.600.52-0.08-0.880.06-0.14-0.96-1.040.06-0.08-0.02-0.02-1.04-1.06CESM2-0.260.350.09-1.700.03-0.01-1.68-1.580.050.100.160.15-1.57-1.43CNRM-CM6-1-0.610.20-0.41-0.77-0.050.05-0.77-1.170.15-0.050.09-0.26-0.82-1.08CNRM-ESM2-1-0.420.18-0.24-0.52-0.04-0.03-0.59-0.830.10-0.020.08-0.14-0.61-0.75EC-Earth3-0.750.32-0.43-0.34-0.04-0.08-0.46-0.890.06-0.020.04-0.37-0.48-0.85GFDL-CM4-0.650.41-0.24-0.540.00-0.09-0.62-0.860.13-0.060.07-0.11-0.68-0.80GFDL-ESM4-0.650.53-0.12-0.590.01-0.06-0.63-0.750.06-0.17-0.11-0.06-0.80-0.86GISS-E2-1-G p1-0.910.19-0.720.060.01-0.94-0.87-1.580.190.100.29-0.52-0.77-1.29GISS-E2-1-G p3-0.970.18-0.79-0.25-0.01-0.06-0.32-1.110.18-0.050.13-0.61-0.37-0.97HadGEM3-GC31-LL-0.770.30-0.47-0.75-0.01-0.07-0.83-1.300.20-0.000.19-0.28-0.83-1.11IPSL-CM6A-LR-0.600.17-0.43-0.26-0.010.06-0.21-0.650.06-0.07-0.01-0.37-0.28-0.65MIROC6-0.480.09-0.39-1.04-0.06-0.01-1.12-1.510.180.380.57-0.21-0.73-0.94MRI-ESM2-0-0.700.19-0.51-1.71-0.09-0.38-2.18-2.680.051.471.52-0.46-0.70-1.16NorESM2-LM-0.420.21-0.21-1.07-0.00-0.10-1.17-1.380.070.140.21-0.14-1.03-1.17NorESM2-MM-0.400.22-0.18-1.100.02-0.17-1.25-1.430.150.040.20-0.03-1.20-1.23UKESM1-0-LL-0.710.34-0.37-0.82-0.01-0.10-0.93-1.300.18-0.010.17-0.19-0.94-1.13Mean-0.630.28-0.35-0.77-0.01-0.13-0.91-1.260.130.100.23-0.23-0.81-1.04SD0.180.120.220.460.040.220.470.460.060.360.360.190.300.20
The total ERFari+aci from the APRP method is -1.04 (±0.20) W m-2, agreeing very well with the ERF estimate of -1.01 (±0.23) W m-2. ERFari+aci is approximately 22 % from ERFari and 78 % from ERFaci and is comprised of an SW contribution of -1.26 W m-2 offset by a LW contribution of +0.23 W m-2. The model spread in both the SW and LW individual components is larger than for the net forcing. This is driven by the four models that include ice cloud interactions that show positive LW ERFaci offset by strong negative SW ERFaci. MRI-ESM2.0 in particular has a very large positive LW ERFaci of +1.47 W m-2, which comes from ice cloud nucleation by black carbon aerosols with temperature below -38∘C in high-level clouds in the tropics . For the SW component the ERFari–ERFaci split is approximately 28 % to 72 %.
Multi-model mean SW ERFari is -0.35 W m-2, comprised of an absorption of +0.28 W m-2 offset by scattering of -0.63 W m-2. The SW ERFaci is -0.91 W m-2, made up of scattering (-0.77 W m-2), absorption (-0.01 W m-2) and cloud fraction change (-0.13 W m-2).
Comparison of ERFari and ERFaci methods
Eight models also archived radiation diagnostics from aerosol-free radiation calls (the double-call method), as recommended by , which allows for separation into ERFari and ERFaci. This can be compared with the APRP estimates in the SW and cloud radiative effect for the LW. Figure shows different methods of estimating ERFaci and ERFari from the aerosol forcing experiment. For ERFaci in both the SW and LW, different methods provide similar estimates. For ERFari, the APRP and double-call methods sometimes disagree on the sign for SW forcing, but this component is relatively small compared to the SW ERFaci where estimates are generally more consistent between APRP and the double call. Similarly in the LW, the CRE and double-call methods produce similar results for ERFaci with larger relative differences for the smaller ERFari component. The double-call method is considered to be quite reliable, but the good agreement between the APRP or CRE and double-call methods suggest that these simpler tools are useful tool to diagnose ERFari and ERFaci from climate models, which is advantageous due to there being no requirement for specialised model diagnostics.
Land use change
Land use ERF is small and not significant at -0.09 (±0.13) W m-2. Forcing and adjustments are difficult to distinguish from zero, and it is unlikely that this forcing played a large role historically for global-mean impacts. In 13 of the 14 models that ran this experiment, land use ERF is negative, and the multi-model mean and standard deviation are affected by a relatively large positive forcing in the NorESM2-LM model. In Fig. S3g we show that this is due to cloud adjustments in this model. This is a consequence of interactive isoprene and monoterpene specified from the land surface changes, causing a reduction in organic matter, reducing cloud condensation nuclei and increasing SW cloud adjustment (unlike for the aerosol forcing experiment, the Twomey effect in response to a land use forcing is treated as an adjustment and not a forcing, because anthropogenic aerosol emissions are not perturbed). In other models, where the ERF is small and negative, it should also be borne in mind that internal variability may make it more difficult to isolate the forcing signal from the noise in free-running simulations , although the multi-model mean is likely to be more robust than individual model results. This experiment was partly motivated by a large land use forcing of -0.4 W m-2 in the CMIP5 HadGEM2-ES model , which showed a large change in regional dust loading that contributed to this forcing. Our multi-model mean ERF of -0.09 W m-2 (-0.12 W m-2 if NorESM2-LM is excluded) agrees well with an observational-constrained analysis from CMIP5 models of -0.11 W m-2 and is within the likely range of the AR5 assessment of -0.15 (-0.05 to -0.25) W m-2.
The radiative forcing from land use change is driven by the resulting change in surface albedo. For example, deforestation for agricultural use converts relatively dark forest cover to brighter cropland, exerting a negative forcing . The surface albedo kernel-derived flux change is taken to be the IRF. It is not a perfect measure, as it includes changes in snow and ice cover over land and any biophysical response, as both changes in land surface temperatures and surface properties can affect snow cover. However, the land surface temperature change is very small in the land use experiment, evidenced by the small land surface temperature adjustment in Fig. . In Fig. S4 we show changes in aerosol optical depth at 550 nm for models that provided this diagnostic. There is diversity in the model aerosol loadings to land use forcing that does not appear to explain the diversity in land use ERF between models. In particular, CanESM5 has a strong aerosol optical depth increase to land use change but a relatively weak ERF of -0.08 W m-2. Changes in surface properties such as how snow cover settles over different land types and the biophysical response are not easy to discern from model output. Again, all available model years are used because stratospheric temperature adjustment does not play a large role.
Comparison of methods to estimate ERFaci and ERFari from the aerosol experiment. CRE/resid. (residual) is the LW cloud radiative effect for ERFaci and the difference of LW ERF and CRE for ERFari. Not all methods are available in all models.
The spatial pattern of land use forcing and adjustments (Fig. ) is generally not significant in many parts of the world due to the small size of the forcing. The exception to this is water vapour and SW cloud adjustments over the Amazon; deforestation from the pre-industrial era to the present day is likely to have reduced evapotranspiration from vegetation, reducing tropospheric humidity and low-level cloud cover. These spatial patterns are also coincident with a decrease in organic carbon loading in NorESM2-LM (Fig. S3).
As Fig. but for present-day land use forcing. Panel (f) represents the IRF in this experiment.
Anthropogenic total
The total anthropogenic ERF for 1850–2014 stands at 2.00 (±0.23) W m-2. The inter-model spread is larger, both in relative and absolute terms, in the total anthropogenic forcing than it is for any of its individual components, suggesting that individual models respond very differently to the same combinations of forcing. In the absence of non-linearities between forcing components, the residual ERF of +0.21 W m-2 from the land use, aerosol and WMGHG components compared to the total anthropogenic would mostly be comprised of ozone forcing, although the sum of individual forcings does not necessarily equal the total forcing in some models . As for the aerosol forcing experiment, there is no significant correlation between total anthropogenic forcing and ECS or TCR.
The total anthropogenic forcing shows the offsetting influences of the greenhouse gas and aerosol forcing components on the ERF, IRF and adjustments. The total anthropogenic ISCCP simulator cloud changes are also a combination of the WMGHG and aerosol contributions, with the net effect being dominated by aerosol. For non-cloud adjustments the combination of strong positive adjustments from the greenhouse gas forcing with a small negative adjustment from the aerosol forcing results in an adjustment that is of comparable magnitude to the IRF. The exception is the GISS-E2-1-G model r1i1p1f1 variant that has a very strong negative cloud adjustment, driven by a large increase in cloud fraction in the aerosol experiment (Table ). Additionally, stratospheric adjustment is stronger for total anthropogenic forcing than for WMGHGs alone, suggesting a role for ozone forcing in contributing to this adjustment.
As for Table but for 1850–2014 anthropogenic forcing.
The pattern of anthropogenic forcing is spatially inhomogeneous; it is positive where aerosol forcing is weak and negative where localised aerosol–cloud effects dominate (Fig. ). The influence of WMGHG forcing on temperature and water vapour adjustments, and of aerosol forcing on the cloud response, is evident. For all forcings, but particularly for land use, aerosol and total anthropogenic, many of the forcing and adjustment terms do not show robust signals regionally. This indicates that adjustments are best considered as global-mean quantities that affect the globally resolved forcing–feedback framework as in Eq. ().
As Fig. but for present-day anthropogenic forcing.
Conclusions
Effective radiative forcing is the driving process behind long-term changes in global-mean surface temperature. As ERF is now preferred to RF, climate models are the best tools we have to determine the heating impacts of various species on the Earth atmosphere system.
From CMIP5 to CMIP6, both CO2 and aerosol forcing have become more consistent across the population of participating models. This has helped to address the concern from CMIP5 that forcing was poorly characterised in CMIP5 models and inconsistently determined . Multi-model mean CO2 and all-WMGHG ERF estimates agree very well with RF estimates from using a line-by-line radiative transfer model (Sect. 5.1 and 5.2). A comprehensive review of aerosol forcing placed the 16 %–84 % uncertainty range in present-day aerosol ERF at -1.60 to -0.65 W m-2. Results from CMIP6 models show a relatively tight spread of -1.37 to -0.63 W m-2 for the full range. Although 17 models is a reasonable sample size of the CMIP6 population, more models may submit forcing results to CMIP6 that would widen this range (and indeed, we would encourage modelling groups to do so). One example is E3SM which did not perform the RFMIP aerosol forcing experiment but where it would be likely that the 1850–2014 aerosol forcing would be more negative than -1.37 W m-2Fig. 25 in. While the increase in 4×CO2 forcing compared to CMIP5 may explain some of the increase in climate sensitivity in CMIP6 models , the model range of present-day aerosol forcing does not, particularly for the upper bound.
We determine a multi-model mean anthropogenic ERF of 2.00 (±0.23) W m-2 for 1850–2014. This is less than the anthropogenic ERF in AR5 for 1850–2011 of 2.24 W m-2 (; although this figure has a wide uncertainty range), and extrapolating trends forward would suggest an anthropogenic ERF of around 2.4 W m-2 from AR5 for 1850–2014. The two main reasons for this difference are a stronger negative aerosol forcing in CMIP6 compared to the AR5 assessment (-1.01 W m-2 in CMIP6 for 1850–2014 versus -0.72 W m-2 in AR5 for 1850–2011) and a weaker ozone forcing (+0.21 W m-2 versus +0.31 W m-2) if residual anthropogenic forcing is attributed to ozone.
Forcing adjustments produce insight into the atmospheric mechanisms that contribute to ERF. Warming of the troposphere results in a negative adjustment due to the increase in outgoing LW radiation, and increasing water vapour counteracts this effect partially by its role as a greenhouse gas. All models agree on tropospheric warming and moistening for WMGHG and all anthropogenic forcing. These tropospheric adjustments are small for aerosol forcing, but models do not agree on the sign of the change.
The instantaneous radiative forcing and cloud adjustments are generally the largest sources of inter-model spread in the forcing component in climate models. Since IRF is not directly calculated in this study, some of this spread may be from residuals in the kernel decomposition, and the true spread in IRF may be smaller than reported here. One strand of RFMIP will include benchmarking of GCM radiative transfer against line-by-line codes. Radiative transfer is a well-grounded theoretical problem where the diversity in line-by-line codes is small , so this component of inter-model diversity has a measurable yardstick for improvement. Cloud responses are more difficult to constrain and exhibit a wide range of behaviour to both greenhouse gas and aerosol forcing. However, progress is beginning to be made. For greenhouse gas forcing, techniques from the climate feedback literature that have observational parallels, such as analysing cloud-controlling factors , can be applied to adjustments. Use of the ISCCP simulator diagnostics with the ISCCP cloud kernel, another method conceptualised by climate feedback investigations , allows for cloud adjustments to be calculated, directly facilitating better inter-model comparison. For aerosol forcing, observational methods exist to determine RFari and RFaci using satellite and reanalysis data . Ultimately, reducing uncertainty in effective radiative forcing will reduce uncertainty in climate projections due to the central role of forcing in driving Earth's global-mean temperature response.
Data availability
RFMIP model data used in this study are freely available from the CMIP6 repository on the Earth System Grid Federation nodes (https://esgf-node.llnl.gov/search/cmip6/, ). The HadGEM3-GA7.1 kernels are available at https://doi.org/10.5281/zenodo.3594673.
The supplement related to this article is available online at: https://doi.org/10.5194/acp-20-9591-2020-supplement.
Author contributions
CJS co-ordinated the project, analysed the data and led the writing of the paper. RJK produced adjustment calculations using radiative kernels. GM provided adjustment calculations using offline radiation simulations. KA provided model results of effective radiative forcing and effective climate sensitivity. WC provided analysis of different definitions of effective radiative forcing. RP and PMF oversaw the RFMIP-ERF project in coordination with the World Climate Research Programme under the CMIP6 protocol. All other authors ran models and published results on the Earth System Grid Foundation, without which this work would not have been possible. All authors contributed to the writing and review.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
We thank Ed Gryspeerdt, Mark Zelinka, Roland Séférian, and the participants and hosts of the Tri-MIP-athlon consortium meetings in Reading in 2018 and Princeton in 2019 for fruitful discussions. We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP6. We thank the climate modelling groups for producing and making available their model output, the Earth System Grid Federation (ESGF) for archiving the data and providing access, and the multiple funding agencies who support CMIP6 and ESGF.
Financial support
Christopher J. Smith was supported by a NERC-IIASA Collaborative Research Fellowship (no. NE/T009381/1). Christopher J. Smith and Piers M. Forster were supported by the European Union's Horizon 2020 Research and Innovation programme (CONSTRAIN project, grant agreement no. 820829). Ryan J. Kramer is supported by an appointment to the NASA Postdoctoral Program at NASA Goddard Space Flight Center. Robert Pincus was supported by the US Department of Energy's Office of Biological and Environmental Research from Lawrence Berkeley National Laboratory (grant no. 7457436). Timothy Andrews, Fiona M. O'Connor, Eddy Robertson and Andy Wiltshire were supported by the Met Office Hadley Centre Climate Programme funded by BEIS and Defra and the Newton Fund through the Met Office Climate Science for Service Partnership Brazil (CSSP Brazil). Fiona M. O'Connor also acknowledges support from the European Union's Horizon 2020 Research and Innovation programme (CRESCENDO project, grant agreement no. 641816). Alf Kirkevåg and Dirk Olivié were supported by the Research Council of Norway (grant nos. 229771, 285003 and 285013), by Notur/NorStore (grant nos. NN2345K and NS2345K) and through the European Union's Horizon 2020 Research and Innovation programme (grant no. 280060). Hideo Shiogama was supported by TOUGOU (MEXT, Japan). The MIROC6 simulations were performed using the Earth Simulator at JAMSTEC and the NEC SX at NIES. The CMIP6 project at IPSL used the HPC resources of TGCC (project gencmip6, grant nos. 2017-R0040110492 and 2018-R0040110492) provided by GENCI (Grand Équipement National de Calcul Intensif). The CESM project is supported primarily by the National Science Foundation. Part of this material is based upon work supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation (cooperative agreement no. 1852977). The simulations with EC-Earth3 were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at the National Supercomputer Centre (NSC) partially funded by the Swedish Research Council (grant agreement no. 2016-07213). Martin Dix and Chloe Mackallah were supported by funding from the Earth Systems and Climate Change Hub of the Australian government's National Environmental Science Program. ACCESS modelling was undertaken with the assistance of resources from the National Computational Infrastructure (NCI Australia), an NCRIS-enabled capability supported by the Australian government.
Review statement
This paper was edited by Hailong Wang and reviewed by two anonymous referees.
References
Albrecht, B. A.: Aerosols, cloud microphysics, and fractional cloudiness,
Science, 245, 1227–1231, 1989.Allen, R., Amiri-Farahani, A., Lamarque, J.-F., Smith, C., Shindell, D.,
Hassan, T., and Chung, C.: npj Climate and Atmospheric Science, 2, 16,
10.1038/s41612-019-0073-9, 2019.Andrews, T., Gregory, J. M., Webb, M. J., and Taylor, K. E.: Forcing, feedbacks
and climate sensitivity in CMIP5 coupled atmosphere-ocean climate models,
Geophys. Res. Lett., 39, L09712, 10.1029/2012GL051607, 2012.Andrews, T., Betts, R. A., Booth, B. B. B., Jones, C. D., and Jones, G. S.:
Effective radiative forcing from historical land use change, Clim. Dynam.,
48, 3489–3505, 10.1007/s00382-016-3280-7, 2017.Andrews, T., Andrews, M. B., Bodas-Salcedo, A., Jones, G. S., Kuhlbrodt, T.,
Manners, J., Menary, M. B., Ridley, J., Ringer, M. A., Sellar, A. A., Senior,
C. A., and Tang, Y.: Forcings, Feedbacks, and Climate Sensitivity in
HadGEM3-GC3.1 and UKESM1, J. Adv. Model. Earth Syst., 11, 4377–4394,
10.1029/2019MS001866, 2019.Armour, K.: Energy budget constraints on climate sensitivity in light of
inconstant climate feedbacks, Nat. Clim. Change, 7, 331–335,
10.1038/nclimate3278, 2017.Bellouin, N., Quaas, J., Morcrette, J.-J., and Boucher, O.: Estimates of aerosol radiative forcing from the MACC re-analysis, Atmos. Chem. Phys., 13, 2045–2062, 10.5194/acp-13-2045-2013, 2013.Bellouin, N., Davies, W., Shine, K. P., Quaas, J., Mülmenstädt, J., Forster, P. M., Smith, C., Lee, L., Regayre, L., Brasseur, G., Sudarchikova, N., Bouarar, I., Boucher, O., and Myhre, G.: Radiative forcing of climate change from the Copernicus reanalysis of atmospheric composition, Earth Syst. Sci. Data, 12, 1649–1677, 10.5194/essd-12-1649-2020, 2020a.Bellouin, N., Quaas, J., Gryspeerdt, E., Kinne, S., Stier, P., Watson-Parris,
D., Boucher, O., Carslaw, K., Christensen, M., Daniau, A.-L., Dufresne,
J.-L., Feingold, G., Fiedler, S., Forster, P., Gettelman, A., Haywood, J.,
Lohmann, U., Malavelle, F., Mauritsen, T., McCoy, D., Myhre, G.,
Mülmenstädt, J., Neubauer, D., Possner, A., Rugenstein, M., Sato, Y.,
Schulz, M., Schwartz, S., Sourdeval, O., Storelvmo, T., Toll, V., Winker, D.,
and Stevens, B.: Bounding global aerosol radiative forcing of climate change,
Rev. Geophys., 58, e2019RG000660, 10.1029/2019RG000660, 2020b.Betts, R.: Offset of the potential carbon sink from boreal forestation by
decreases in surface albedo, Nature, 408, 187–190, 10.1038/35041545,
2000.Betts, R. A.: Biogeophysical impacts of land use on present-day climate:
near-surface temperature change and radiative forcing, Atmos. Sci. Lett., 2,
39–51, 10.1006/asle.2001.0037, 2001.
Bi, D., Dix, M., Marsland, S., O’Farrell, S., Sullivan, A., Bodman, R., Law,
R., Harman, I., Srbinovsky, J., Rashid, H., Dobrohotoff, P., Mackallah, C.,
Yan, H., Hirst, A., Savita, A., Dias, F. B., Woodhouse, M., Fiedler, R., and
Heerdegen, A.: Configuration and spin-up of ACCESS-CM2, the new generation
Australian Community Climate and Earth System Simulator Coupled Model,
Journal of Southern Hemisphere Earth Systems Science, submitted, 2020.Block, K. and Mauritsen, T.: Forcing and feedback in the MPI-ESM-LR coupled
model under abruptly quadrupled CO2, J. Adv. Model. Earth Syst., 5,
696–691, 10.1002/jame.20041, 2013.
Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster,
P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh,
S., Sherwood, S., Stevens, B., and Zhang, X.: Clouds and Aerosols, in:
Climate Change 2013: The Physical Science Basis. Contribution of Working
Group I to the Fifth Assessment Report of the Intergovernmental Panel on
Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M.,
Allen, S., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P.,
571–658, Cambridge University Press, Cambridge, United Kingdom and New York,
NY, USA, 2013.Boucher, O., Servonnat, J., Albright, A. L., Aumont, O., Balkanski, Y.,
Bastrikov, V., Bekki, S., Bonnet, R., Bony, S., Bopp, L., Braconnot, P.,
Brockmann, P., Cadule, P., Caubel, A., Cheruy, F., Cozic, A., Cugnet, D.,
D'Andrea, F., Davini, P., de Lavergne, C., Denvil, S., Deshayes, J.,
Devilliers, M., Ducharne, A., Dufresne, J.-L., Dupont, E., Éthé, C.,
Fairhead, L., Falletti, L., Foujols, M.-A., Gardoll, S., Gastineau, G.,
Ghattas, J., Grandpeix, J.-Y., Guenet, B., Guez, L., Guilyardi, E.,
Guimberteau, M., Hauglustaine, D., Hourdin, F., Idelkadi, A., Joussaume, S.,
Kageyama, M., Khodri, M., Krinner, G., Lebas, N., Levavasseur, G., Lévy,
C., Li, L., Lott, F., Lurton, T., Luyssaert, S., Madec, G., Madeleine, J.-B.,
Maignan, F., Marchand, M., Marti, O., Mellul, L., Meurdesoif, Y., Mignot, J.,
Musat, I., Ottlé, C., Peylin, P., Planton, Y., Polcher, J., Rio, C.,
Rousset, C., Sepulchre, P., Sima, A., Swingedouw, D., Thiéblemont, R.,
Traoré, A.-K., Vancoppenolle, M., Vial, J., Vialard, J., Viovy, N., and
Vuichard, N.: Presentation and evaluation of the IPSL-CM6A-LR climate
model, J. Adv. Model. Earth Syst., 12, e2019MS002010. 10.1029/2019MS002010, 2020.Chung, E.-S. and Soden, B. J.: An assessment of methods for computing radiative
forcing in climate models, Environ Res. Lett., 10, 074004,
10.1088/1748-9326/10/7/074004, 2015.Chylek, P., Vogelsang, T. J., Klett, J. D., Hengartner, N., Higdon, D., Lesins,
G., and Dubey, M. K.: Indirect Aerosol Effect Increases CMIP5 Models’
Projected Arctic Warming, J. Climate, 29, 1417–1428,
10.1175/JCLI-D-15-0362.1, 2016.Collins, W. J., Lamarque, J.-F., Schulz, M., Boucher, O., Eyring, V., Hegglin, M. I., Maycock, A., Myhre, G., Prather, M., Shindell, D., and Smith, S. J.: AerChemMIP: quantifying the effects of chemistry and aerosols in CMIP6, Geosci. Model Dev., 10, 585–607, 10.5194/gmd-10-585-2017, 2017.Copernicus Climate Change Service: ERA5: Fifth generation of ECMWF
atmospheric reanalyses of the global climate,
available at: https://cds.climate.copernicus.eu/cdsapp#!/home (last access: 28 October 2019), 2017.Danabasoglu, G., Lamarque, J.-F., Bacmeister, J., Bailey, D. A., DuVivier,
A. K., Edwards, J., Emmons, L. K., Fasullo, J., Garcia, R., Gettelman, A.,
Hannay, C., Holland, M. M., Large, W. G., Lauritzen, P. H., Lawrence, D. M.,
Lenaerts, J. T. M., Lindsay, K., Lipscomb, W. H., Mills, M. J., Neale, R.,
Oleson, K. W., Otto-Bliesner, B., Phillips, A. S., Sacks, W., Tilmes, S., van
Kampenhout, L., Vertenstein, M., Bertini, A., Dennis, J., Deser, C., Fischer,
C., Fox-Kemper, B., Kay, J. E., Kinnison, D., Kushner, P. J., Larson, V. E.,
Long, M. C., Mickelson, S., Moore, J. K., Nienhouse, E., Polvani, L., Rasch,
P. J., and Strand, W. G.: The Community Earth System Model Version 2 (CESM2),
J. Adv. Model. Earth Syst., 12, e2019MS001916, 10.1029/2019MS001916,
2020.
Dunne et al.: The GFDL Earth System Model version 4.1 (GFDL-ESM4.1): Model
description and simulation characteristics, J. Adv. Model. Earth Syst., in
prep., 2020.Edwards, J. M. and Slingo, A.: Studies with a flexible new radiation code. I:
Choosing a configuration for a large-scale model, Q. J. Roy. Meteor. Soc.,
122, 689–719, 10.1002/qj.49712253107, 1996.Etminan, M., Myhre, G., Highwood, E. J., and Shine, K. P.: Radiative forcing of
carbon dioxide, methane, and nitrous oxide: A significant revision of the
methane radiative forcing, Geophys. Res. Lett., 43, 12614–12623,
10.1002/2016GL071930, 2016.Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958, 10.5194/gmd-9-1937-2016, 2016.Flynn, C. M. and Mauritsen, T.: On the climate sensitivity and historical warming evolution in recent coupled model ensembles, Atmos. Chem. Phys., 20, 7829–7842, 10.5194/acp-20-7829-2020, 2020.Forster, P., Richardson, T., Maycock, A., Smith, C., Samset, B., Myhre, G.,
Andrews, T., Pincus, R., and Schulz, M.: Recommendations for diagnosing
effective radiative forcing from climate models from CMIP6, J. Geophys.
Res., 121, 12460–12475, 10.1002/2016JD025320, 2016.
Forster, P., Maycock, A., McKenna, C., and Smith, C.: Latest climate models
confirm need for urgent mitigation, Nat. Clim. Change, 10, 7–10, 2020.Forster, P. M., Andrews, T., Good, P., Gregory, J. M., Jackson, L. S., and
Zelinka, M.: Evaluating adjusted forcing and model spread for historical and
future scenarios in the CMIP5 generation of climate models, J. Geophys.
Res.-Atmos., 118, 1139–1150, 10.1002/jgrd.50174, 2013.Forster, P. M. d. F. and Shine, K. P.: Radiative forcing and temperature trends
from stratospheric ozone changes, J. Geophys. Res.-Atmos., 102,
10841–10855, 10.1029/96JD03510, 1997.Ghan, S. J.: Technical Note: Estimating aerosol effects on cloud radiative forcing, Atmos. Chem. Phys., 13, 9971–9974, 10.5194/acp-13-9971-2013, 2013.Gidden, M. J., Riahi, K., Smith, S. J., Fujimori, S., Luderer, G., Kriegler, E., van Vuuren, D. P., van den Berg, M., Feng, L., Klein, D., Calvin, K., Doelman, J. C., Frank, S., Fricko, O., Harmsen, M., Hasegawa, T., Havlik, P., Hilaire, J., Hoesly, R., Horing, J., Popp, A., Stehfest, E., and Takahashi, K.: Global emissions pathways under different socioeconomic scenarios for use in CMIP6: a dataset of harmonized emissions trajectories through the end of the century, Geosci. Model Dev., 12, 1443–1475, 10.5194/gmd-12-1443-2019, 2019.Golaz, J.-C., Caldwell, P. M., Van Roekel, L. P., Petersen, M. R., Tang, Q.,
Wolfe, J. D., Abeshu, G., Anantharaj, V., Asay-Davis, X. S., Bader, D. C.,
Baldwin, S. A., Bisht, G., Bogenschutz, P. A., Branstetter, M., Brunke,
M. A., Brus, S. R., Burrows, S. M., Cameron-Smith, P. J., Donahue, A. S.,
Deakin, M., Easter, R. C., Evans, K. J., Feng, Y., Flanner, M., Foucar,
J. G., Fyke, J. G., Griffin, B. M., Hannay, C., Harrop, B. E., Hoffman,
M. J., Hunke, E. C., Jacob, R. L., Jacobsen, D. W., Jeffery, N., Jones,
P. W., Keen, N. D., Klein, S. A., Larson, V. E., Leung, L. R., Li, H.-Y.,
Lin, W., Lipscomb, W. H., Ma, P.-L., Mahajan, S., Maltrud, M. E., Mametjanov,
A., McClean, J. L., McCoy, R. B., Neale, R. B., Price, S. F., Qian, Y.,
Rasch, P. J., Reeves Eyre, J. E. J., Riley, W. J., Ringler, T. D., Roberts,
A. F., Roesler, E. L., Salinger, A. G., Shaheen, Z., Shi, X., Singh, B.,
Tang, J., Taylor, M. A., Thornton, P. E., Turner, A. K., Veneziani, M., Wan,
H., Wang, H., Wang, S., Williams, D. N., Wolfram, P. J., Worley, P. H., Xie,
S., Yang, Y., Yoon, J.-H., Zelinka, M. D., Zender, C. S., Zeng, X., Zhang,
C., Zhang, K., Zhang, Y., Zheng, X., Zhou, T., and Zhu, Q.: The DOE E3SM
Coupled Model Version 1: Overview and Evaluation at Standard Resolution,
J. Adv. Model. Earth Syst., 11, 2089–2129,
10.1029/2018MS001603, 2019.Gregory, J., Ingram, W., Palmer, M., Jones, G., Stott, P., Thorpe, R., Lowe,
J., Johns, T., and Williams, K.: A new method for diagnosing radiative
forcing and climate sensitivity, Geophys. Res. Lett., 31, L03205,
10.1029/2003GL018747, 2004.Gryspeerdt, E., Mülmenstädt, J., Gettelman, A., Malavelle, F. F., Morrison, H., Neubauer, D., Partridge, D. G., Stier, P., Takemura, T., Wang, H., Wang, M., and Zhang, K.: Surprising similarities in model and observational aerosol radiative forcing estimates, Atmos. Chem. Phys., 20, 613–623, 10.5194/acp-20-613-2020, 2020.Hansen, J., Sato, M., Ruedy, R., Nazarenko, L., Lacis, A., Schmidt, G. A.,
Russell, G., Aleinov, I., Bauer, M., Bauer, S., Bell, N., Cairns, B., Canuto,
V., Chandler, M., Cheng, Y., Del Genio, A., Faluvegi, G., Fleming, E.,
Friend, A., Hall, T., Jackman, C., Kelley, M., Kiang, N., Koch, D., Lean, J.,
Lerner, J., Lo, K., Menon, S., Miller, R., Minnis, P., Novakov, T., Oinas,
V., Perlwitz, J., Perlwitz, J., Rind, D., Romanou, A., Shindell, D., Stone,
P., Sun, S., Tausnev, N., Thresher, D., Wielicki, B., Wong, T., Yao, M., and
Zhang, S.: Efficacy of climate forcings, J. Geophys. Res.-Atmos., 110,
d18104, 10.1029/2005JD005776, 2005.Haustein, K., Allen, M., Forster, P., Otto, F., Mitchell, D., Matthews, H., and
Frame, D.: A real-time Global Warming Index, Sci. Rep., 7, 15417, 10.1038/s41598-017-14828-5,
2017.Held, I. M. and Shell, K. M.: Using Relative Humidity as a State Variable in
Climate Feedback Analysis, J. Climate, 25, 2578–2582,
10.1175/JCLI-D-11-00721.1, 2012.Held, I. M. and Soden, B. J.: Water vapor feedback and global warming, Annu.
Rev. Energ. Env., 25, 441–475, 10.1146/annurev.energy.25.1.441, 2000.Held, I. M., Guo, H., Adcroft, A., Dunne, J. P., Horowitz, L. W., Krasting, J.,
Shevliakova, E., Winton, M., Zhao, M., Bushuk, M., Wittenberg, A. T., Wyman,
B., Xiang, B., Zhang, R., Anderson, W., Balaji, V., Donner, L., Dunne, K.,
Durachta, J., Gauthier, P. P. G., Ginoux, P., Golaz, J.-C., Griffies, S. M.,
Hallberg, R., Harris, L., Harrison, M., Hurlin, W., John, J., Lin, P., Lin,
S.-J., Malyshev, S., Menzel, R., Milly, P. C. D., Ming, Y., Naik, V.,
Paynter, D., Paulot, F., Rammaswamy, V., Reichl, B., Robinson, T., Rosati,
A., Seman, C., Silvers, L. G., Underwood, S., and Zadeh, N.: Structure and
Performance of GFDL's CM4.0 Climate Model, J. Adv. Model. Earth Syst., 11, 3691–3727,
10.1029/2019MS001829, 2019.Hoesly, R. M., Smith, S. J., Feng, L., Klimont, Z., Janssens-Maenhout, G., Pitkanen, T., Seibert, J. J., Vu, L., Andres, R. J., Bolt, R. M., Bond, T. C., Dawidowski, L., Kholod, N., Kurokawa, J.-I., Li, M., Liu, L., Lu, Z., Moura, M. C. P., O'Rourke, P. R., and Zhang, Q.: Historical (1750–2014) anthropogenic emissions of reactive gases and aerosols from the Community Emissions Data System (CEDS), Geosci. Model Dev., 11, 369–408, 10.5194/gmd-11-369-2018, 2018.Huang, Y.: On the longwave climate feedbacks, J. Climate, 26, 7603–7610,
10.1175/JCLI-D-13-00025.1, 2013.Johnson, B. T., Haywood, J. M., and Hawcroft, M. K.: Are Changes in Atmospheric
Circulation Important for Black Carbon Aerosol Impacts on Clouds,
Precipitation, and Radiation?, J. Geophys. Res.-Atmos., 124, 7930–7950,
10.1029/2019JD030568,
2019.Kamae, Y. and Watanabe, M.: On the robustness of tropospheric adjustment in
CMIP5 models, Geophys. Res. Lett., 39, L23808, 10.1029/2012GL054275,
2012.Kelley, M., Schmidt, G. A., Nazarenko, L., Miller, R. L., Bauer, S. E., Ruedy,
R., Russell, G. L., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Cesana,
G., Cheng, Y., Clune, T. L., Cook, B., Cruz, C. A., Genio, A. D. D.,
Elsaesser, G. S., Faluvegi, G., Kiang, N. Y., Kim, D., Lacis, A. A.,
Leboissetier, A., LeGrande, A. N., Lo, K. K., Marshall, J. C., McDermid, S.,
Matthews, E. E., Mezuman, K., Murray, L. T., Oinas, V., Orbe, C.,
Garcia-Pando, C. P., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A.,
Shindell, D. T., Sun, S., Tausnev, N., Tsigaridis, K., Tselioudis, G., Weng,
E., Wu, J., and Yao, M.-S.: GISS-E2.1: Configurations and Climatology, J.
Adv. Model. Earth Syst., 12, e2019MS002025, 10.1029/2019MS002025, 2020.Kiehl, J. T.: Twentieth century climate model response and climate sensitivity,
Geophys. Res. Lett., 34, L22710, 10.1029/2007GL031383, 2007.Kirkevåg, A., Grini, A., Olivié, D., Seland, Ø., Alterskjær, K., Hummel, M., Karset, I. H. H., Lewinschal, A., Liu, X., Makkonen, R., Bethke, I., Griesfeller, J., Schulz, M., and Iversen, T.: A production-tagged aerosol module for Earth system models, OsloAero5.3 – extensions and updates for CAM5.3-Oslo, Geosci. Model Dev., 11, 3945–3982, 10.5194/gmd-11-3945-2018, 2018.Klein, S. and Jakob, C.: Validation and Sensitivities of Frontal Clouds
Simulated by the ECMWF Model, Mon. Weather Rev., 127, 2514–2531,
10.1175/1520-0493(1999)127<2514:VASOFC>2.0.CO;2, 1999.Klein, S., Hall, A., Norris, J., and Pincus, R.: Low-Cloud Feedbacks from
Cloud-Controlling Factors: A Review., Surv. Geophys., 38, 1307–1329, 10.1007/s10712-017-9433-3, 2017.Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.: Historical (1850–2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039, 10.5194/acp-10-7017-2010, 2010.Lejeune, Q., Davin, E. L., Duveiller, G., Crezee, B., Meier, R., Cescatti, A., and Seneviratne, S. I.: Biases in the albedo sensitivity to deforestation in CMIP5 models and their impacts on the associated historical Radiative Forcing, Earth Syst. Dynam. Discuss., 10.5194/esd-2019-94, in review, 2020.Loeb, N. G., Doelling, D. R., Wang, H., Su, W., Nguyen, C., Corbett, J. G.,
Liang, L., Mitrescu, C., Rose, F. G., and Kato, S.: Clouds and the Earth’s
Radiant Energy System (CERES) Energy Balanced and Filled (EBAF)
Top-of-Atmosphere (TOA) Edition-4.0 Data Product, J. Climate, 31, 895–918,
10.1175/JCLI-D-17-0208.1, 2018.Mauritsen, T., Bader, J., Becker, T., Behrens, J., Bittner, M., Brokopf, R.,
Brovkin, V., Claussen, M., Crueger, T., Esch, M., Fast, I., Fiedler, S.,
Fläschner, D., Gayler, V., Giorgetta, M., Goll, D. S., Haak, H.,
Hagemann, S., Hedemann, C., Hohenegger, C., Ilyina, T., Jahns, T.,
Jimenéz-de-la Cuesta, D., Jungclaus, J., Kleinen, T., Kloster, S., Kracher,
D., Kinne, S., Kleberg, D., Lasslop, G., Kornblueh, L., Marotzke, J., Matei,
D., Meraner, K., Mikolajewicz, U., Modali, K., Möbis, B., Müller, W. A.,
Nabel, J. E. M. S., Nam, C. C. W., Notz, D., Nyawira, S.-S., Paulsen, H.,
Peters, K., Pincus, R., Pohlmann, H., Pongratz, J., Popp, M., Raddatz, T. J.,
Rast, S., Redler, R., Reick, C. H., Rohrschneider, T., Schemann, V., Schmidt,
H., Schnur, R., Schulzweida, U., Six, K. D., Stein, L., Stemmler, I.,
Stevens, B., von Storch, J.-S., Tian, F., Voigt, A., Vrese, P., Wieners,
K.-H., Wilkenskjeld, S., Winkler, A., and Roeckner, E.: Developments in the
MPI-M Earth System Model version 1.2 (MPI-ESM1.2) and Its Response to
Increasing CO2, J. Adv. Model. Earth Syst., 11, 998–1038,
10.1029/2018MS001400, 2019.Meinshausen, M., Vogel, E., Nauels, A., Lorbacher, K., Meinshausen, N., Etheridge, D. M., Fraser, P. J., Montzka, S. A., Rayner, P. J., Trudinger, C. M., Krummel, P. B., Beyerle, U., Canadell, J. G., Daniel, J. S., Enting, I. G., Law, R. M., Lunder, C. R., O'Doherty, S., Prinn, R. G., Reimann, S., Rubino, M., Velders, G. J. M., Vollmer, M. K., Wang, R. H. J., and Weiss, R.: Historical greenhouse gas concentrations for climate modelling (CMIP6), Geosci. Model Dev., 10, 2057–2116, 10.5194/gmd-10-2057-2017, 2017.Mülmenstädt, J., Gryspeerdt, E., Salzmann, M., Ma, P.-L., Dipu, S., and Quaas, J.: Separating radiative forcing by aerosol-cloud interactions and rapid cloud adjustments in the ECHAM-HAMMOZ aerosol-climate model using the method of partial radiative perturbations, Atmos. Chem. Phys., 19, 15415–15429, 10.5194/acp-19-15415-2019, 2019.
Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang,
J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock,
A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and Natural
Radiative Forcing, in: Climate Change 2013: The Physical Science Basis.
Contribution of Working Group I to the Fifth Assessment Report of the
Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D.,
Plattner, G.-K., Tignor, M., Allen, S., Boschung, J., Nauels, A., Xia, Y.,
Bex, V., and Midgley, P., Cambridge University Press,
Cambridge, United Kingdom and New York, NY, USA, 659–740, 2013.
Myhre, G., Forster, P., Samset, B., Hodnebrog, Ø., Sillmann, J.,
Aalbergsjø, S., Andrews, T., Boucher, O., Faluvegi, G., Fläschner, D.,
, Iversen, T., Kasoar, M., Kharin, V., Kirkevå g, A., Lamarque, J.-F.,
Olivié, D., Richardson, T. B., Shindell, D., Shine, K. P., Stjern, C. W.,
Takemura, T., Voulgarakis, A., and Zwiers, F.: PDRMIP: A precipitation driver
and response model intercomparison project-Protocol and preliminary
results, B. Am. Meteorol. Soc., 98, 1185–1198, 2017.O'Connor, F. M., Abraham, N. L., Dalvi, M., Folberth, G., Griffiths, P., Hardacre, C., Johnson, B. T., Kahana, R., Keeble, J., Kim, B., Morgenstern, O., Mulcahy, J. P., Richardson, M. G., Robertson, E., Seo, J., Shim, S., Teixeira, J. C., Turnock, S., Williams, J., Wiltshire, A., and Zeng, G.: Assessment of pre-industrial to present-day anthropogenic climate forcing in UKESM1, Atmos. Chem. Phys. Discuss., 10.5194/acp-2019-1152, in review, 2020.Oshima, N., Yukimoto, S., Deushi, M., Koshiro, T., Kawai, H., Tanaka, T. Y.,
and Yoshida, K.: Global and Arctic Radiative Forcing of Anthropogenic Gases
and Aerosols in MRI-ESM2.0, Prog. Earth Planet. Sci., 7, 38, 10.1186/s40645-020-00348-w, 2020.Pendergrass, A. G., Conley, A., and Vitt, F. M.: Surface and top-of-atmosphere radiative feedback kernels for CESM-CAM5, Earth Syst. Sci. Data, 10, 317–324, 10.5194/essd-10-317-2018, 2018.Pincus, R., Mlawer, E. J., Oreopoulos, L., Ackerman, A. S., Baek, S., Brath,
M., Buehler, S. A., Cady-Pereira, K. E., Cole, J. N. S., Dufresne, J.-L.,
Kelley, M., Li, J., Manners, J., Paynter, D. J., Roehrig, R., Sekiguchi, M.,
and Schwarzkopf, D. M.: Radiative flux and forcing parameterization error in
aerosol-free clear skies, Geophys. Res. Lett., 42, 5485–5492,
10.1002/2015GL064291, 2015.Pincus, R., Forster, P. M., and Stevens, B.: The Radiative Forcing Model Intercomparison Project (RFMIP): experimental protocol for CMIP6, Geosci. Model Dev., 9, 3447–3460, 10.5194/gmd-9-3447-2016, 2016.Richardson, T., Forster, P., Maycock, A., Smith, C., Wood, T., Andrews, T.,
Boucher, O., Faluvegi, G., Fläschner, D., Hodnebrog, O., Kasoar, M.,
Kirkevag, A., Lamarque, J.-F., Mülmenstädt, J., Olivié, D.,
Samset, B., Shawki, D., Shindell, D., Takemura, T., and Voulgarakis, A.:
Efficacy of climate forcings in PDRMIP models, J. Geophys. Res., 124, 12824–12844,
10.1029/2019JD030581, 2019.Santer, B. D., Wehner, M. F., Wigley, T. M. L., Sausen, R., Meehl, G. A.,
Taylor, K. E., Ammann, C., Arblaster, J., Washington, W. M., Boyle, J. S.,
and Brüggemann, W.: Contributions of Anthropogenic and Natural Forcing to
Recent Tropopause Height Changes, Science, 301, 479–483,
10.1126/science.1084123, 2003.Seland, Ø., Bentsen, M., Seland Graff, L., Olivié, D., Toniazzo, T., Gjermundsen, A., Debernard, J. B., Gupta, A. K., He, Y., Kirkevåg, A., Schwinger, J., Tjiputra, J., Schancke Aas, K., Bethke, I., Fan, Y., Griesfeller, J., Grini, A., Guo, C., Ilicak, M., Hafsahl Karset, I. H., Landgren, O., Liakka, J., Onsum Moseid, K., Nummelin, A., Spensberger, C., Tang, H., Zhang, Z., Heinze, C., Iverson, T., and Schulz, M.: The Norwegian Earth System Model, NorESM2 – Evaluation of theCMIP6 DECK and historical simulations, Geosci. Model Dev. Discuss., 10.5194/gmd-2019-378, in review, 2020.Séférian, R., Nabat, P., Michou, M., Saint-Martin, D., Voldoire, A., Colin,
J., Decharme, B., Delire, C., Berthet, S., Chevallier, M., Sénési, S.,
Franchisteguy, L., Vial, J., Mallet, M., Joetzjer, E., Geoffroy, O.,
Guérémy, J.-F., Moine, M.-P., Msadek, R., Ribes, A., Rocher, M., Roehrig,
R., Salas-y Mélia, D., Sanchez, E., Terray, L., Valcke, S., Waldman, R.,
Aumont, O., Bopp, L., Deshayes, J., Éthé, C., and Madec, G.: Evaluation of
CNRM Earth-System model, CNRM-ESM 2-1: role of Earth system processes in
present-day and future climate, J. Adv. Model. Earth Syst., 11, 4182–4227,
10.1029/2019MS001791, 2019.Sellar, A. A., Jones, C. G., Mulcahy, J., Tang, Y., Yool, A., Wiltshire, A.,
O'Connor, F. M., Stringer, M., Hill, R., Palmieri, J., Woodward, S., de Mora,
L., Kuhlbrodt, T., Rumbold, S., Kelley, D. I., Ellis, R., Johnson, C. E.,
Walton, J., Abraham, N. L., Andrews, M. B., Andrews, T., Archibald, A. T.,
Berthou, S., Burke, E., Blockley, E., Carslaw, K., Dalvi, M., Edwards, J.,
Folberth, G. A., Gedney, N., Griffiths, P. T., Harper, A. B., Hendry, M. A.,
Hewitt, A. J., Johnson, B., Jones, A., Jones, C. D., Keeble, J., Liddicoat,
S., Morgenstern, O., Parker, R. J., Predoi, V., Robertson, E., Siahaan, A.,
Smith, R. S., Swaminathan, R., Woodhouse, M. T., Zeng, G., and Zerroukat, M.:
UKESM1: Description and evaluation of the UK Earth System Model, J. Adv.
Model. Earth Syst., 11, 4513–4558, 10.1029/2019MS001739, 2019.Shell, K., Kiehl, J., and Shields, C.: Using the Radiative Kernel Technique to
Calculate Climate Feedbacks in NCAR's Community Atmospheric Model, J.
Climate, 21, 2269–2282, 10.1175/2007JCLI2044.1, 2008.Sherwood, S. C., Bony, S., Boucher, O., Bretherton, C., Forster, P. M.,
Gregory, J. M., and Stevens, B.: Adjustments in the Forcing-Feedback
Framework for Understanding Climate Change, B. Am. Meteorol. Soc., 96,
217–228, 10.1175/BAMS-D-13-00167.1, 2015.Shine, K. P., Cook, J., Highwood, E. J., and Joshi, M. M.: An alternative to
radiative forcing for estimating the relative importance of climate change
mechanisms, Geophys. Res. Lett., 30, 10.1029/2003GL018141, 2047, 2003.Smith, C.: HadGEM3-GA7.1 radiative kernels, Zenodo, 10.5281/zenodo.3594673 (last access: 13 May 2020), 2019.Smith, C. J., Forster, P. M., Allen, M., Leach, N., Millar, R. J., Passerello, G. A., and Regayre, L. A.: FAIR v1.3: a simple emissions-based impulse response and carbon cycle model, Geosci. Model Dev., 11, 2273–2297, 10.5194/gmd-11-2273-2018, 2018a.Smith, C. J., Kramer, R. J., Myhre, G., Forster, P. M., Soden, B. J., Andrews,
T., Boucher, O., Faluvegi, G., Fläschner, D., Hodnebrog, O., Kasoar, M.,
Kharin, V., Kirkevåg, A., Lamarque, J.-F., Mülmenstädt, J.,
Olivié, D., Richardson, T., Samset, B. H., Shindell, D., Stier, P.,
Takemura, T., Voulgarakis, A., and Watson-Parris, D.: Understanding Rapid
Adjustments to Diverse Forcing Agents, Geophys. Res. Lett., 45,
12023–12031, 10.1029/2018GL079826, 2018b.Smith, C. J., Kramer, R. J., and Sima, A.: The HadGEM3-GA7.1 radiative kernel: the importance of a well-resolved stratosphere, Earth Syst. Sci. Data Discuss., 10.5194/essd-2019-254, accepted, 2020.Soden, B., Held, I., Colman, R., Shell, K., Kiehl, J., and Shields, C.:
Quantifying Climate Feedbacks Using Radiative Kernels, J. Climate, 21,
3504–3520, 10.1175/2007JCLI2110.1, 2008.Soden, B. J., Collins, W. D., and Feldman, D. R.: Reducing uncertainties in
climate models, Science, 361, 326–327, 10.1126/science.aau1864, 2018.Stjern, C. W., Samset, B. H., Myhre, G., Forster, P. M., Hodnebrog, Ø.,
Andrews, T., Boucher, O., Faluvegi, G., Iversen, T., Kasoar, M., Kharin, V.,
Kirkevåg, A., Lamarque, J.-F., Olivié, D., Richardson, T., Shawki, D.,
Shindell, D., Smith, C. J., Takemura, T., and Voulgarakis, A.: Rapid
Adjustments Cause Weak Surface Temperature Response to Increased Black Carbon
Concentrations, J. Geophys. Res.-Atmos., 122, 11462–11481,
10.1002/2017JD027326, 2017.Stouffer, R. J., Eyring, V., Meehl, G. A., Bony, S., Senior, C., Stevens, B.,
and Taylor, K. E.: CMIP5 Scientific Gaps and Recommendations for CMIP6,
B. Am. Meteorol. Soc., 98, 95–105,
10.1175/BAMS-D-15-00013.1, 2017.Swart, N. C., Cole, J. N. S., Kharin, V. V., Lazare, M., Scinocca, J. F., Gillett, N. P., Anstey, J., Arora, V., Christian, J. R., Hanna, S., Jiao, Y., Lee, W. G., Majaess, F., Saenko, O. A., Seiler, C., Seinen, C., Shao, A., Sigmond, M., Solheim, L., von Salzen, K., Yang, D., and Winter, B.: The Canadian Earth System Model version 5 (CanESM5.0.3), Geosci. Model Dev., 12, 4823–4873, 10.5194/gmd-12-4823-2019, 2019.Tang, T., Shindell, D., Faluvegi, G., Myhre, G., Olivié, D., Voulgarakis,
A., Kasoar, M., Andrews, T., Boucher, O., Forster, P., Hodnebrog, O.,
Iversen, T., Kirkevåg, A., Lamarque, J.-F., Richardson, T., Samset, B.,
Stjern, C., Takemura, T., and Smith, C.: Comparison of Effective Radiative
Forcing Calculations Using Multiple Methods, Drivers, and Models, J. Geophys.
Res.-Atmos., 124, 4382–4394, 10.1029/2018JD030188,
2019.Tatebe, H., Ogura, T., Nitta, T., Komuro, Y., Ogochi, K., Takemura, T., Sudo, K., Sekiguchi, M., Abe, M., Saito, F., Chikira, M., Watanabe, S., Mori, M., Hirota, N., Kawatani, Y., Mochizuki, T., Yoshimura, K., Takata, K., O'ishi, R., Yamazaki, D., Suzuki, T., Kurogi, M., Kataoka, T., Watanabe, M., and Kimoto, M.: Description and basic evaluation of simulated mean state, internal variability, and climate sensitivity in MIROC6, Geosci. Model Dev., 12, 2727–2765, 10.5194/gmd-12-2727-2019, 2019.
Taylor, K., Crucifix, M., Braconnot, P., Hewitt, C., Doutriaux, C., Broccoli,
A., Mitchell, J., and Webb, M.: Estimating shortwave radiative forcing and
response in climate models, J. Climate, 20, 2530–2543, 2007.Thornhill, G. D., Collins, W. J., Kramer, R. J., Olivié, D., O'Connor, F., Abraham, N. L., Bauer, S. E., Deushi, M., Emmons, L., Forster, P., Horowitz, L., Johnson, B., Keeble, J., Lamarque, J.-F., Michou, M., Mills, M., Mulcahy, J., Myhre, G., Nabat, P., Naik, V., Oshima, N., Schulz, M., Smith, C., Takemura, T., Tilmes, S., Wu, T., Zeng, G., and Zhang, J.: Effective Radiative forcing from emissions of reactive gases and aerosols – a multimodel comparison, Atmos. Chem. Phys. Discuss., 10.5194/acp-2019-1205, in review, 2020.Tokarska, K. B., Gillett, N. P., Arora, V. K., Lee, W. G., and Zickfeld, K.:
The influence of non-CO2 forcings on cumulative carbon emissions budgets,
Environ Res. Lett., 13, 034039, 10.1088/1748-9326/aaafdd, 2018.
Twomey, S.: The influence of pollution on the shortwave albedo of clouds, J.
Atmos. Sci., 34, 1149–1152, 1977.Vial, J., Dufresne, J.-L., and Bony, S.: On the interpretation of inter-model
spread in CMIP5 climate sensitivity estimates, Clim. Dynam., 41,
3339–3362, 10.1007/s00382-013-1725-9, 2013.Voldoire, A., Saint-Martin, D., Sénési, S., Decharme, B., Alias, A.,
Chevallier, M., Colin, J., Guérémy, J.-F., Michou, M., Moine, M.-P., Nabat,
P., Roehrig, R., Salas y Mélia, D., Séférian, R., Valcke, S., Beau, I.,
Belamari, S., Berthet, S., Cassou, C., Cattiaux, J., Deshayes, J., Douville,
H., Ethé, C., Franchistéguy, L., Geoffroy, O., Lévy, C., Madec, G.,
Meurdesoif, Y., Msadek, R., Ribes, A., Sanchez-Gomez, E., Terray, L., and
Waldman, R.: Evaluation of CMIP6 DECK Experiments With CNRM-CM6-1, J. Adv.
Model. Earth Syst., 11, 2177–2213, 10.1029/2019MS001683, 2019.Webb, M., Senior, C., Bony, S., and Morcrette, J.-J.: Combining ERBE and ISCCP
data to assess clouds in the Hadley Centre, ECMWF and LMD atmospheric climate
models, Clim. Dynam., 17, 905–922, 10.1007/s003820100157, 2001.Williams, K. D., Copsey, D., Blockley, E. W., Bodas-Salcedo, A., Calvert, D.,
Comer, R., Davis, P., Graham, T., Hewitt, H. T., Hill, R., Hyder, P., Ineson,
S., Johns, T. C., Keen, A. B., Lee, R. W., Megann, A., Milton, S. F., Rae, J.
G. L., Roberts, M. J., Scaife, A. A., Schiemann, R., Storkey, D., Thorpe, L.,
Watterson, I. G., Walters, D. N., West, A., Wood, R. A., Woollings, T., and
Xavier, P. K.: The Met Office Global Coupled Model 3.0 and 3.1 (GC3.0 and
GC3.1) Configurations, J. Adv. Model. Earth Syst., 10, 357–380,
10.1002/2017MS001115, 2018.World Climate Research Programme (WCRP): Coupled Model Intercomparison Project, Phase 6, Earth System Grid Federation, Lawrence Livermore National Laboratory, available at: https://esgf-node.llnl.gov/search/cmip6/, last access: 13 May 2020.Wyser, K., van Noije, T., Yang, S., von Hardenberg, J., O'Donnell, D., and Döscher, R.: On the increased climate sensitivity in the EC-Earth model from CMIP5 to CMIP6, Geosci. Model Dev. Discuss., 10.5194/gmd-2019-282, in review, 2019.Yukimoto, S., Kawai, H., Koshiro, T., Oshima, N., Yoshida, K., Urakawa, S.,
Tsujino, H., Deushi, M., Tanaka, T., Hosaka, M., Yabu, S., Yoshimura, H.,
Shindo, E., Mizuta, R., Obata, A., Adachi, Y., and Ishii, M.: The
Meteorological Research Institute Earth System Model Version 2.0, MRI-ESM2.0:
Description and Basic Evaluation of the Physical Component, J. Meteorol. Soc.
Jpn. Ser. II, 97, 931–965, 10.2151/jmsj.2019-051, 2019.Zelinka, M., Klein, S., and Hartmann, D.: Computing and Partitioning Cloud
Feedbacks Using Cloud Property Histograms, Part I: Cloud Radiative Kernels,
J. Climate, 25, 3715–3735, 10.1175/JCLI-D-11-00248.1, 2012.
Zelinka, M. D., Andrews, T., Forster, P. M., and Taylor, K. E.: Quantifying
components of aerosol-cloud-radiation interactions in climate models, J.
Geophys. Res.-Atmos., 119, 7599–7615, 10.1002/2014JD021710,
2014.Zelinka, M. D., Myers, T. A., McCoy, D. T., Po-Chedley, S., Caldwell, P. M.,
Ceppi, P., Klein, S. A., and Taylor, K. E.: Causes of Higher Climate
Sensitivity in CMIP6 Models, Geophys. Res. Lett., 47,
e2019GL085782, 10.1029/2019GL085782, 2020.