We present top-of-atmosphere and surface radiative kernels based on the atmospheric component (GA7.1) of the HadGEM3 general circulation model developed by the UK Met Office. We show that the utility of radiative kernels for forcing adjustments in idealised CO2 perturbation experiments is greatest where there is sufficiently high resolution in the stratosphere in both the target climate model and the radiative kernel. This is because stratospheric cooling to a CO2 perturbation continues to increase with height, and low-resolution or low-top kernels or climate model output are unable to fully resolve the full stratospheric temperature adjustment. In the sixth phase of the Coupled Model Intercomparison Project (CMIP6), standard atmospheric model data are available up to 1 hPa on 19 pressure levels, which is a substantial advantage compared to CMIP5. We show in the IPSL-CM6A-LR model where a full set of climate diagnostics are available that the HadGEM3-GA7.1 kernel exhibits linear behaviour and the residual error term is small, as well as from a survey of kernels available in the literature that in general low-top radiative kernels underestimate the stratospheric temperature response. The HadGEM3-GA7.1 radiative kernels are available at https://doi.org/10.5281/zenodo.3594673 .
Introduction
Radiative kernels describe how a small change in an atmospheric state variable affects the Earth's energy balance . They allow an analysis of climate feedbacks or forcing adjustments from standardised climate model diagnostics such as those from the Coupled Model Intercomparison Projects (CMIPs). The use of radiative kernels is efficient, removing the need for time- and memory-consuming calculations of climate feedbacks online through partial radiative perturbation calculations or offline using a stand-alone version of the model radiative transfer code .
A radiative kernel KX is in effect a four-dimensional (time, height, latitude, longitude) array describing how radiation fluxes R change with an atmospheric state variable X:
KX(t,z,y,x)=∂R∂X(t,z,y,x).
Although strictly not a partial differential equation, this statement provides a concise written form and is used by others e.g.. R may be upwelling, downwelling or net, shortwave or longwave radiation changes at any atmospheric level. Most commonly net top-of-atmosphere (TOA), surface and tropopause-level fluxes are of greatest interest. X here represents atmospheric temperature (Ta), surface (skin) temperature (Ts), water vapour (q) and surface albedo (α). Although other (non-cloud) variables may also be relevant, the majority of adjustments are expected to be captured under this framework . For determining adjustments to a radiative forcing AX, the kernel KX is multiplied by the change in atmospheric state variable ΔX between two integrations of a climate model such that
AX=KXΔX.ΔX is calculated as the difference of two atmosphere-only climate model integrations using climatological sea-surface temperatures and sea-ice distributions, one of which is driven by a forcing perturbation (e.g. a quadrupling of CO2) and the other a control. For temperature and albedo the adjustment is linear with ΔX and logarithmic for water vapour (, and Supplement of , describe how the adjustment to water vapour is applied in practice). For determining climate feedbacks λX, the perturbation is normalised by the change in global mean near-surface air temperature T such that
λX=KX∂X∂T.
The individual contributions from each feedback component λX contribute the total climate feedback λ=λTa+λTs+λq+λα+λc where “c” represents cloud feedback in the forcing-feedback representation of the Earth's energy budget ΔN=F-λΔT. Here, ΔN is the Earth's energy imbalance and F the effective radiative forcing. Likewise, the effective radiative forcing can be decomposed into
F=Fi+ATa+ATs+Aq+Aα+Ac,
with Fi being the instantaneous radiative forcing (IRF).
Usage of radiative kernels assumes that radiative perturbations change linearly with changes in atmospheric state. Where perturbations are small, linearity is an appropriate assumption both for feedbacks and adjustments .
Cloud adjustments and feedbacks cannot be determined directly using atmospheric state kernels. They may be diagnosed using the cloud kernel based on International Satellite Cloud Climatology Project (ISCCP) simulator diagnostics or from the residual of all-sky and clear-sky radiative kernels . For adjustments this calculation is
Ac=(F-Fclr)-(Fi-Ficlr)-∑X∈{Ta,Ts,q,α}(AX-AXclr),
where the “clr” superscript represents fluxes calculated in the absence of clouds. In Eq. (), the instantaneous radiative forcing must be known or estimated. This method is commonly used, requiring the production of all-sky and clear-sky kernel sets to calculate all-sky and clear-sky adjustments.
This paper introduces the top-of-atmosphere and surface radiative kernels from the HadGEM3-GA7.1 model. HadGEM3-GA7.1 has 85 vertical levels up to 85 km vertical height (about 0.005 hPa). The construction of the HadGEM3-GA7.1 kernel was motivated by the observation that adjustments to a doubling of CO2 in models participating in the Precipitation Driver and Response Model Intercomparison Project (PDRMIP; ) were around 0.3 W m-2 larger using the ECMWF-Oslo kernel than other kernels used in the same study Supplement Fig. 3. The ECMWF-Oslo kernel was built from ECMWF-Interim reanalysis data which have 60 vertical levels up to 0.1 hPa. Most other kernels used in were derived from climate models with a lower model top, and it is possible that low-top kernels were underestimating the stratospheric adjustment in .
For stratospheric temperature adjustments we compare the HadGEM3-GA7.1 kernel to other kernels in the literature using available 4×CO2 results from climate models contributing to the Radiative Forcing Model Intercomparison Project (RFMIP). We find that, in general, kernels based upon climate models with a high stratospheric resolution are better at resolving the stratospheric adjustment to a CO2 forcing.
Construction of the HadGEM3-GA7.1 radiative kernel
One year of a pre-industrial, atmosphere-only (i.e. with climatological sea-surface temperatures and sea-ice distributions) integration of the HadGEM3-GA7.1 general circulation model was run. HadGEM3-GA7.1 is the atmospheric component of the HadGEM3-GC3.1 physical model and UKESM1.0 Earth system model that represents the UK research community's contribution to CMIP6. The model was run at LL (N96) resolution with a latitude–longitude grid of 1.25∘ by 1.875∘ and 85 vertical levels extending up to 85 km (approximately 0.005 hPa) and a native model time step of 20 min. A pre-industrial base climatology for the kernels was chosen because the first identified use case for this kernel set was the RFMIP-ERF Tier 1 single-forcing experiments . These experiments compare positive (e.g. greenhouse gas) and negative (e.g. aerosol) forcing perturbations with respect to a pre-industrial control baseline.
Model diagnostics of air temperature, specific humidity, surface (skin) temperature, surface albedo (ratio of broadband upwelling to downwelling shortwave surface radiation), model level pressure, surface pressure, cloud fraction, cloud water content, cloud ice content, effective (time-averaged) solar zenith angle and grid box daylight fraction every two model hours were saved. Two-hourly sampling is considered to give an appropriate representation of the diurnal cycle of clouds and reduce biases due to variations in solar geometry with longitude compared to longer time steps, while keeping computational demands to a minimum. These model outputs were transplanted into an offline version of the SOCRATES radiative transfer code (version 17.03; ) and top-of-atmosphere and surface radiative fluxes calculated for each 2 h time step in both the shortwave and longwave spectra, for all sky and clear sky. SOCRATES is a broadband radiation code that uses six bands in the shortwave and nine bands in the longwave and is the same radiation scheme used in the HadGEM3-GC3.1 and UKESM1.0 climate models. Aerosols were neglected and greenhouse gases, including the prescribed CMIP6 monthly climatology in ozone concentrations, were set to their pre-industrial (1850) values. Following the protocol for RFMIP , sea-surface temperatures and sea-ice distributions from 50 years of the HadGEM3-GC3.1 coupled model were used to build the climatology . The kernel is therefore dependent on two aspects of the HadGEM3-GA7.1 model: the pre-industrial background climatology (including clouds) and the broadband version of the radiation code.
To build the kernel, each vertical level of the model on each 2 h time step was perturbed separately, firstly by 1 K for air temperature and secondly by a perturbation in specific humidity that maintains relative humidity for an increase in 1 K (without actually changing the layer temperature). The surface temperature and surface albedo were also perturbed by 1 K and 1 % (additive) individually each time step. For each perturbation, surface and TOA fluxes are again saved for clear sky and all sky in the shortwave (SW) and longwave (LW), and the difference compared to the control simulation gives the radiative kernel for each model level or surface. Building the kernels took in total approximately 3 months of computing time on 24 processors on the University of Leeds “cluj” Linux cluster. Running the base HadGEM3 model at higher (MM) resolution, with 5 times as many grid points as LL, was not determined to be necessary, as kernels are usually ultimately regridded to the resolutions of other CMIP6 models which are approximately as coarse as LL (e.g. Table 1 in ).
Following this, the air temperature and water vapour kernel outputs were normalised by multiplying by 10 000 Pa/pthick, where pthick [Pa] is the thickness of each level in pressure coordinates. This allows the 85-level native model kernel to be reduced down to the 19-level standard CMIP6 pressure levels by providing a weighted average contribution to each pressure level. The kernels are further averaged by month. In the 19-level format they can be used with standard “Amon” model output from any CMIP6 model, which is one of the key advantages of radiative kernels.
Kernel profilesTop-of-atmosphere kernels
Figure shows the TOA radiative kernels for HadGEM3-GA7.1 for clear sky and all sky. The air temperature, all-sky kernel (Fig. a) shows a peak in cooling in the tropical upper troposphere, showing the importance of this region for changes in radiative balance. There are also substantial contributions to the TOA radiation balance from the lower troposphere in the mid-latitudes. For clear sky (Fig. b) there is more latitude-height homogeneity in the troposphere, showing the impact of removing clouds. A key feature of the air temperature kernels is the increasing strength of the LW outgoing radiation with increasing stratospheric height. The temperature kernel is negative throughout the atmosphere, in keeping with the fact that an increase in temperature results in additional Planck emission of LW radiation to space.
Water vapour kernels (Fig. c, d) also show a peak in the upper tropical troposphere, which is opposite in sign to the negative temperature adjustment owing to the fact that water vapour is a significant greenhouse gas. In contrast to the temperature kernel, the water vapour kernel is very insensitive in the dry upper stratosphere.
The impact of cloud masking is more easily seen for the surface temperature kernels (Fig. e, f) and surface albedo kernels (Fig. g, h).
Top-of-atmosphere radiative kernels from HadGEM3-GA7.1.
Surface kernels
Surface kernels are most useful for determining precipitation adjustments and feedbacks , where the precipitation adjustment is proportional to the atmospheric absorption, calculated as the difference in TOA and surface adjustments. Figure shows the surface radiative kernels for HadGEM3-GA7.1 for clear sky and all sky. Both the air temperature (Fig. a, b) and water vapour (Fig. c, d) kernels are more sensitive for perturbations close to the surface than higher in the atmosphere (note non-linear colour scales). Cloud masking for the surface temperature kernel has less of an effect for surface fluxes than for TOA fluxes (Fig. e, f), whereas the surface albedo kernel shows quite a similar spatial pattern (Fig. g, h) to its the TOA counterpart.
Surface radiative kernels from HadGEM3-GA7.1.
Air temperature radiative kernels available in the literature, truncated at 300 hPa to show the stratospheric temperature contribution. Blank areas are above the top level of the kernel.
Comparison to other kernels for stratospheric temperature
Figure shows the air temperature kernel for the stratosphere and upper troposphere for a selection of kernels available in the literature (Fig. ). In all cases, radiative kernels have been interpolated from their native vertical resolution (except for CCSM4, which is available only on the standard 17 CMIP5 pressure levels) to the 19 CMIP6 pressure levels for consistency with CMIP6 model output. For our calculations of stratospheric temperature adjustment, where kernels do not extend up to the 1 hPa top level of CMIP6 model output, kernels have been extended upwards using the value from the highest level where data do exist, but in Fig. missing data have been masked out. This extending upwards of the top level has been applied previously in adjustment calculations where the top level of the climate model is higher than the top level of the kernel (e.g. in ). However, extending the top level of a radiative kernel upwards cannot make up for the fact that more radiation is emitted to space from the upper stratosphere for each additional K of temperature change. For kernels built from underlying atmospheric profiles where the top of the profile is not sufficiently high or with too coarse a resolution in the stratosphere, this additional upper stratospheric cooling is missed. In Fig. , it can be seen that the kernels based on a high-top atmospheric model with a high number of native model levels – ECHAM6, ECMWF-Oslo, ECMWF-RRTM and HadGEM3-GA7.1 – have a marked increase in both the magnitude and the rate of negative LW outgoing flux at the 5 and 1 hPa levels.
Radiative kernels considered in this study.
Native modelTopThird levelBase modelvertical levelslevel (hPa)(hPa)ReferenceBMRC178.7553.63CCSM4171030CESM303.6414.36ECHAM5191050.39ECHAM6470.00990.11ECMWF-RRTM2416ECMWF-Oslo600.110.5GFDL253.3253.63HadGEM2382.9913.02HadGEM3-GA7.1850.0050.03this study
The consequences for a greenhouse-gas-induced stratospheric cooling are such that the additional stratospheric adjustment from greater cooling high in the stratosphere is not accounted for with either kernels or models that are truncated too low. Figure shows the atmospheric temperature anomalies simulated in atmosphere-only simulations from CMIP6 models participating in RFMIP-ERF Tier 1 experiments for a 30-year time-slice simulation where CO2 concentrations are quadrupled relative to a pre-industrial control (piClim-4×CO2). Stratospheric cooling continues to increase above 10 hPa in all models where data are available. A similar pattern of stratospheric cooling occurs in the piClim-ghg experiment which evaluates forcing from present-day greenhouse gases (not shown). Standard CMIP6 diagnostics call for model output on 19 pressure levels: 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10, 5 and 1 hPa, whereas in CMIP5 the standard set of 17 pressure levels did not include 5 and 1 hPa. Therefore, CMIP5 models were missing important additional stratospheric cooling where kernels were used for adjustment calculations.
Atmospheric temperature differences for RFMIP models for the piClim-4×CO2 experiment minus piClim-control experiment.
The truncation of stratospheric height in “low-top” radiative kernels (those with a top level lower in altitude than 1 hPa) has substantial consequences for adjustments to a greenhouse gas forcing. Figure shows the stratospheric temperature adjustment for 4×CO2 in 13 models contributing to RFMIP. A simplified tropopause definition is used here, borrowed from , of a linear-in-latitude ramp from 100 hPa at the Equator to 300 hPa at the poles. There is a spread of around 1 W m-2 in calculated stratospheric temperature adjustment for each model using the full range of kernels, which is about 13 % of the effective radiative forcing (ERF) for a quadrupling of CO2 from these models . It can be seen in Fig. that the kernel estimates for stratospheric adjustment to CO2 forcing are clustered into two groups and one high outlier for most models. The one model where kernel estimates are not clearly separated into high and low clusters is the GFDL-CM4 model, for which data are missing at 1 hPa. The low-top radiative kernels, with the exception of GFDL and ECHAM5, produce substantially lower estimates of the stratospheric temperature adjustment than the “high-top” kernels (HadGEM3-GA7.1, ECHAM6 and ECMWF-Oslo). The high outlier, ECMWF-RRTM, has a large flux change at the 1 hPa level. We show in Sect. that adjustments calculated using the HadGEM3-GA7.1 kernel in the IPSL-CM6A-LR model for a quadrupled CO2 experiment provide small residuals (i.e. the adjustments are appropriately captured), suggesting that assuming there are no compensating errors, low-top kernels would underestimate the stratospheric temperature response and produce larger residuals.
Stratospheric temperature adjustments calculated from RFMIP piClim-4×CO2 experiments using all kernels available in this study.
To investigate the source of the differences between kernels in more detail, Fig. shows the vertically integrated kernels from 1 hPa, 10 hPa, the full stratosphere (, definition) and the full atmosphere. As previously, kernels that do not include the uppermost layers have been extended based on the highest layer for which data are reported. Figure a gives the temperature adjustment from a uniform 1 K increase throughout the atmosphere. Here, there is little variation between kernels, and much of the total adjustment comes from the troposphere where the bulk of the atmosphere is present. Notably, the whole atmosphere adjustment to a 1 K temperature increase is of the order -3 W m-2, approximating (if slightly underestimating) the Planck feedback (-3.2 W m-2 K-1).
Figure b shows the layer contributions when each kernel is applied to the atmospheric profile from the IPSL-CM6A-LR model for piClim-4×CO2. In contrast to the isothermal case, the top 1 hPa, top 10 hPa and stratospheric temperature adjustments show substantial diversity between kernels. The large contributions from these layers follow from Fig. , where the 1 and 10 hPa layers are of the order 10 and 20 K cooler in the 4×CO2 run than in the control. The HadGEM3-GA7.1, GFDL, ECMWF-Oslo, ECMWF-RRTM and ECHAM6 kernels have a larger adjustment from the top 10 hPa than the other kernels. Despite being nominally a low-top kernel, the GFDL kernel has a similar magnitude and gradient of cooling between 10 and 5 hPa to the high-top kernels (Fig. ) and produces a similar result of total stratospheric adjustment to HadGEM3-GA7.1. The ECHAM5 kernel has more cooling around the 100 hPa level than any other kernel (Fig. ) and also has stratospheric adjustments that are similar to HadGEM3-GA7.1.
The differences between models for the top 10 hPa propagate through to the tropopause level and whole atmosphere. Differences in tropospheric adjustments between kernels (whole atmosphere minus stratosphere) are small, showing that choice of atmospheric base state and radiative transfer is not critical for tropospheric temperature adjustments. This analysis gives further confidence that the choice of radiative kernel is not that important in climate feedback studies (assuming state changes are sufficiently small to remain linear; ) as stratospheric temperature differences are small when differencing coupled atmosphere–ocean simulations and kernels show similar behaviour in the troposphere.
Layer contributions to the TOA temperature adjustment in each kernel considered in this study. From top to bottom the four bars show 1 hPa to TOA, 10 hPa to TOA, tropopause to TOA (i.e. stratospheric adjustment) and surface to TOA. (a) The effect of a uniform 1 K increase in atmospheric temperature on TOA fluxes, as diagnosed directly from the radiative kernel. (b) The kernels convoluted with atmospheric profiles from the IPSL-CM6A-LR model under an atmosphere-only quadrupled CO2 run (RFMIP piClim-4×CO2 minus piClim-control experiment).
IPSL-CM6A-LR double-call results for 4×CO2 experiments. IRFs are given in watts per square metre (W m-2).
Base climatologySecond callIRF LWIRF SWIRF NetIRF LW CSIRF SW CSIRF Net CSPre-industrial4×CO23.660.834.495.020.465.484×CO2Pre-industrial4.940.815.756.260.466.72MeanMean4.300.825.125.640.466.10
IPSL-CM6A-LR forcing and adjustments for the 4×CO2 experiment using the HadGEM3-GA7.1 kernel. Fluxes are given in W m-2. ATa strat. and ATa trop. are stratospheric and tropospheric temperature adjustments.
ERFIRFATa strat.ATa trop.ATsAqAαAc (Eq. )Ac (ISCCP kernel)ϵclrϵallLW5.334.302.74-1.38-0.490.52-0.66-0.750.280.38SW2.680.820.110.181.601.81-0.02-0.23Net8.015.122.74-1.38-0.490.630.180.941.060.260.15Linearity of the HadGEM3-GA7.1 kernel
This section shows the linear behaviour of the HadGEM3-GA7.1 kernel used with IPSL-CM6A-LR climate model output. IPSL-CM6A-LR is chosen as all required diagnostics are available in this model to evaluate linearity, including estimates of the IRF from double calls, and ISCCP simulator diagnostics for clouds. This model is also representative of the RFMIP population, with 4×CO2 ERF and adjustments close to the multi-model average .
Linearity is evaluated by the size of the residual term, ϵ, which is any TOA flux changes not explained by instantaneous radiative forcing or kernel-calculated adjustments. A guideline of linearity for the kernel method is that the residual should be within 10 % of the ERF. We take two different approaches to calculate the residual. The first assumes that we have knowledge of the cloud adjustment term Ac (e.g. from the ISCCP simulator or NASA A-Train kernels convoluted with ISCCP simulated cloud output from climate models ), and knowledge of the IRF (Fi) from a double call in the online model. This is rare in practice, as double calls and ISCCP cloud diagnostics are not routinely archived on the Earth System Grid Federation (ESGF), the main distributed data source for CMIP6 output (although model participation is improving for ISCCP clouds). From this “perfect information” method we can calculate an all-sky residual ϵall as
ϵall=F-Fi-ATa-ATs-Aq-Aα-Ac.
The second method does not require knowledge of the cloud adjustments but does require knowledge of the clear-sky IRF, and in this case the clear-sky residual term ϵclr can be calculated as
ϵclr=Fclr-Ficlr-ATaclr-ATsclr-Aqclr-Aαclr.
This clear-sky residual definition is more common in the literature , although often the clear-sky IRF is estimated rather than calculated directly. However, in some circumstances, clear-sky and all-sky IRF are known to be identically zero (e.g. in the LW spectrum to a change in the solar constant; ). In these cases, Eq. () can be used with Fi=Ficlr=0 to determine Ac, which is then plugged into Eq. () to calculate ϵall.
IPSL-CM6A-LR used two sets of double calls. In the RFMIP piClim-4×CO2 experiment (30-year time-slice atmosphere-only run with quadrupled CO2) the second radiation call saw a pre-industrial CO2 concentration. In the piClim-control experiment (pre-industrial atmosphere-only run) the second radiation call saw 4×CO2. The resulting IRF depends on the direction of the double call and is related to the underlying 4×CO2 or pre-industrial climatology. The 4×CO2 climatology that sees a pre-industrial second radiation call results in an IRF that is 1.26 W m-2 greater than the pre-industrial climatology that sees a 4×CO2 second radiation call (Table ). We take the mean of the two simulations to be the IRF.
Table shows ERF, IRF, adjustments and residuals using the HadGEM3-GA7.1 radiative kernel with the IPSL-CM6A-LR model output and cloud adjustments from the ISCCP simulator kernel. We use ISCCP kernel values in the calculation of ϵall and also compare the ISCCP values to the cloud-masking estimate of cloud adjustment from Eq. (). For LW forcing the residuals are 0.28 W m-2 for ϵclr and 0.38 W m-2 for ϵall. Residuals are present possibly due to a slight breakdown in the linearity assumption for a forcing as large as 4×CO2; however, the residuals are comfortably within the 10 % linearity guideline. SW residuals are also within 10 % of the ERF, with ϵclr being particularly small. For the net fluxes, forcings add but residuals partly cancel, such that ϵclr and ϵall are 3.2 % and 1.9 % of the ERF respectively.
IPSL-CM6A-LR forcing and adjustments for the 4×CO2 experiment using the HadGEM3-GA7.1 kernel and ISCCP cloud kernel compared to the multi-model mean and standard deviation (1σ) in from 11 CMIP5 models. Fluxes are given in W m-2.
Our results in Table can be compared with the results of for 11 CMIP5 models. The instantaneous forcing and tropospheric adjustments from IPSL-CM6A-LR with the HadGEM3-GA7.1 kernel are in line with the CMIP5 forcing and adjustments except for the stratospheric temperature adjustment, which is outside the 2σ range from CMIP5 models (Table ). As discussed in , the 4×CO2 ERF in available CMIP6 models (7.98 W m-2) is (non-significantly) greater than in corresponding CMIP5 sstClim4×CO2 experiments (7.53 W m-2), which is also the case for IPSL-CM6A-LR in our comparison to . IPSL-CM6A-LR is near the centre of the CMIP6 range for stratospheric temperature adjustment and so is a representative model of this ensemble. This could suggest stratospheric temperature adjustment increase as one driver of the increase in ERF between CMIP5 and CMIP6 models, although as most CMIP5 output is only on 17 model levels up to 10 hPa a formal comparison is difficult.
The stratospheric temperature adjustment is the only adjustment estimate that varies significantly between radiative kernels . If a kernel that did not resolve stratospheric temperature adjustment well was used instead of HadGEM3-GA7.1, this adjustment (ATa strat. in Table ) would be smaller, and the overall residuals for LW and net responses to 4×CO2 would be larger. From Fig. it can be seen that some kernels produce a stratospheric temperature adjustment around 0.7 W m-2 lower than the HadGEM3-GA7.1 kernel, which would lead to residuals of the order 1 W m-2 using these kernels, or more than 10 % of the ERF.
Data availability
The HadGEM3-GA7.1 radiative kernels are available at 10.5281/zenodo.3594673.
Conclusions
This paper serves two purposes – it introduces the radiative kernel based on the high-top HadGEM3-GA7.1 general circulation model, and it compares estimates of the stratospheric temperature adjustment obtained with a variety of different radiative kernels for quadrupled CO2 experiments. The HadGEM3-GA7.1 kernel is the first to our knowledge that has been produced using a CMIP6-era model, with a focus on the 19 pressure level diagnostics available in CMIP6 output, although other kernels in the literature have used high-top atmospheric profiles . Radiative kernels are produced for both top-of-atmosphere and surface fluxes and are available on the native 85-level hybrid height grid in addition to the 19 CMIP6 pressure levels.
We show that there is a significant diversity, of about 1 W m-2 or 13 % of the ERF for a quadrupling of CO2, for estimates of stratospheric temperature adjustments to CO2 depending on the radiative kernel used to derive the estimate. As tropospheric and land surface adjustments vary little between kernels to a variety of different forcing agents , these differences in stratospheric temperature adjustments lead to differing estimates of the total adjustment, and also of the IRF if it is calculated as a residual .
Climate feedbacks are little affected by the choice of kernel, due to the fact that stratospheric temperatures readjust quickly to an imposed forcing in coupled model simulations .
While only one model (IPSL-CM6A-LR) archived IRF from a double call and a rigorous multi-model test is not possible, we show that the HadGEM3-GA7.1 kernel diagnoses IRF and adjustments with a small residual owing to the increased stratospheric resolution available compared to many CMIP3- and CMIP5-era kernels. We suggest that radiative kernels with a higher stratospheric resolution and model top are better able to fully capture stratospheric adjustments to CO2 forcing in general and generate smaller residuals. This effect has become more prominent with the additional 5 and 1 hPa model levels archived as standard in processed CMIP6 model output compared to CMIP5. Archiving instantaneous radiative forcing from more models would be beneficial to further test the linearity assumption of the radiative kernel method.
Author contributions
CJS produced the HadGEM3-GA7.1 radiative kernel and led the writing of the manuscript. RJK provided calculations of stratospheric adjustment to 4×CO2 for all kernels considered in this paper. AS provided double-call results from the IPSL-CM6A-LR model.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
This work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk, last access: 4 September 2020).
Financial support
Christopher J. Smith was supported by a NERC/IIASA Collaborative Research Fellowship (NE/T009381/1) and the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 820829 (CONSTRAIN project). Ryan J. Kramer is supported by an appointment to the NASA Postdoctoral Program at the NASA Goddard Space Flight Center.
Review statement
This paper was edited by David Carlson and reviewed by Yi Huang and one anonymous referee.
ReferencesAndrews, 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 Sy., 11, 4377–4394,
10.1029/2019MS001866, 2019.Block, K. and Mauritsen, T.: Forcing and feedback in the MPI-ESM-LR coupled
model under abruptly quadrupled CO2, J. Adv. Model. Earth Sy., 5, 676–691,
10.1002/jame.20041, 2013.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, 2015a.Chung, E.-S. and Soden, B. J.: An Assessment of Direct Radiative Forcing,
Radiative Adjustments, and Radiative Feedbacks in Coupled Ocean–Atmosphere
Models, J. Climate, 28, 4152–4170, 10.1175/JCLI-D-14-00436.1,
2015b.Colman, R. A. and McAvaney, B. J.: On tropospheric adjustment to forcing and
climate feedbacks, Clim. Dynam., 36, 1649, 10.1007/s00382-011-1067-4,
2011.Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P.,
Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P.,
Bechtold, P., Beljaars, A. C. M., van de Berg, I., Biblot, J., Bormann, N.,
Delsol, C., Dragani, R., Fuentes, M., Greer, A. J., Haimberger, L., Healy, S.
B., Hersbach, H., Holm, E. V., Isaksen, L., Kallberg, P., Kohler, M.,
Matricardi, M., McNally, A. P., Mong-Sanz, B. M., Morcette, J.-J., Park,
B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thepaut, J. N., and Vitart,
F.: The ERA-Interim reanalysis: Configuration and performance of the data
assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597, 10.1002/qj.828, 2011.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.Huang, Y.: On the Longwave Climate Feedbacks, J. Climate, 26, 7603–7610,
10.1175/JCLI-D-13-00025.1, 2013.Huang, Y., Xia, Y., and Tan, X.: On the pattern of CO2 radiative forcing and
poleward energy transport, J. Geophys. Res.-Atmos., 122,
10578–10593, 10.1002/2017JD027221, 2017.Jonko, A. K., Shell, K. M., Sanderson, B. M., and Danabasoglu, G.: Climate
Feedbacks in CCSM3 under Changing CO2 Forcing. Part I: Adapting the Linear
Radiative Kernel Technique to Feedback Calculations for a Broad Range of
Forcings, J. Climate, 25, 5260–5272, 10.1175/JCLI-D-11-00524.1, 2012.
Manners, J., Edwards, J. M., Hill, P., and Thelen, J.: SOCRATES (Suite Of
Community RAdiative Transfer codes based on Edwards and Slingo) technical
guide, Tech. rep., Met Office, UK, 2015.Mulcahy, J. P., Jones, C., Sellar, A., Johnson, B., Boutle, I. A., Jones, A.,
Andrews, T., Rumbold, S. T., Mollard, J., Bellouin, N., Johnson, C. E.,
Williams, K. D., Grosvenor, D. P., and McCoy, D. T.: Improved Aerosol
Processes and Effective Radiative Forcing in HadGEM3 and UKESM1, J. Adv.
Model. Earth Sy., 10, 2786–2805, 10.1029/2018MS001464, 2018.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., Kirkevag, A., Lamarque, J., Olivié, D.,
Richardson, T., Shindell, D., Shine, K., Stjern, C., 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, 10.1175/BAMS-D-16-0019.1, 2017.Myhre, G., Kramer, R. J., Smith, C. J., Hodnebrog, O., Forster, P., Soden,
B. J., Samset, B. H., Stjern, C. W., Andrews, T., Boucher, O., Faluvegi, G.,
Fläschner, D., Kasoar, M., Kirkevåg, A., Lamarque, J.-F., Olivié,
D., Richardson, T., Shindell, D., Stier, P., Takemura, T., Voulgarakis, A.,
and Watson-Parris, D.: Quantifying the Importance of Rapid Adjustments for
Global Precipitation Changes, Geophys. Res. Lett., 45, 11399–11405,
10.1029/2018GL079474, 2018.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., 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.Previdi, M.: Radiative feedbacks on global precipitation, Environ. Res. Lett., 5,
025211, 10.1088/1748-9326/5/2/025211, 2010.Sanderson, B. M. and Shell, K. M.: Model-Specific Radiative Kernels for
Calculating Cloud and Noncloud Climate Feedbacks, J. Climate, 25, 7607–7624,
10.1175/JCLI-D-11-00726.1, 2012.Shell, K. M., Kiehl, J. T., and Shields, C. A.: 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.Smith, C.: HadGEM3-GA7.1 radiative kernels, Zenodo, 10.5281/zenodo.3594673, 2019.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, 2018.Smith, C. J., Kramer, R. J., Myhre, G., Alterskjær, K., Collins, W., Sima, A., Boucher, O., Dufresne, J.-L., Nabat, P., Michou, M., Yukimoto, S., Cole, J., Paynter, D., Shiogama, H., O'Connor, F. M., Robertson, E., Wiltshire, A., Andrews, T., Hannay, C., Miller, R., Nazarenko, L., Kirkevåg, A., Olivié, D., Fiedler, S., Lewinschal, A., Mackallah, C., Dix, M., Pincus, R., and Forster, P. M.: Effective radiative forcing and adjustments in CMIP6 models, Atmos. Chem. Phys., 20, 9591–9618, 10.5194/acp-20-9591-2020, 2020.Soden, B. J., Held, I. M., Colman, R., Shell, K. M., Kiehl, J. T., and Shields,
C. A.: 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.
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.Wetherald, R. T. and Manabe, S.: Cloud Feedback Processes in a General
Circulation Model, J. Atmos. Sci., 45, 1397–1416,
10.1175/1520-0469(1988)045<1397:CFPIAG>2.0.CO;2, 1988.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 Sy., 10, 357–380,
10.1002/2017MS001115, 2018.Yue, Q., Kahn, B. H., Fetzer, E. J., Schreier, M., Wong, S., Chen, X., and
Huang, X.: Observation-Based Longwave Cloud Radiative Kernels Derived from
the A-Train, J. Climate, 29, 2023–2040,
10.1175/JCLI-D-15-0257.1, 2016.Zelinka, M. D., Klein, S. A., and Hartmann, D. L.: 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.Zhang, M. and Huang, Y.: Radiative Forcing of Quadrupling CO2, J.
Climate, 27, 2496–2508, 10.1175/JCLI-D-13-00535.1, 2014.