Small global-mean cooling due to volcanic radiative forcing

In both the observational record and atmosphere-ocean general circulation model (AOGCM) simulations of the last ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}150 years, short-lived negative radiative forcing due to volcanic aerosol, following explosive eruptions, causes sudden global-mean cooling of up to ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}0.3 K. This is about five times smaller than expected from the transient climate response parameter (TCRP, K of global-mean surface air temperature change per W m−2 of radiative forcing increase) evaluated under atmospheric CO2 concentration increasing at 1 % yr−1. Using the step model (Good et al. in Geophys Res Lett 38:L01703, 2011. doi:10.1029/2010GL045208), we confirm the previous finding (Held et al. in J Clim 23:2418–2427, 2010. doi:10.1175/2009JCLI3466.1) that the main reason for the discrepancy is the damping of the response to short-lived forcing by the thermal inertia of the upper ocean. Although the step model includes this effect, it still overestimates the volcanic cooling simulated by AOGCMs by about 60 %. We show that this remaining discrepancy can be explained by the magnitude of the volcanic forcing, which may be smaller in AOGCMs (by 30 % for the HadCM3 AOGCM) than in off-line calculations that do not account for rapid cloud adjustment, and the climate sensitivity parameter, which may be smaller than for increasing CO2 (40 % smaller than for 4 × CO2 in HadCM3).


Introduction
The global-mean surface air temperature, expressed as the difference T from an unperturbed steady state, is widely used as an indicator of the magnitude of global climate change, both in observations and in simulations of the past and future.Changes in T occur as a result of unforced (internally generated) variability of the climate system on all timescales, and in response to radiative forcing of the climate system.
Anthropogenic radiative forcing, mainly due to wellmixed greenhouse gases and tropospheric aerosols, has increased monotonically and rather smoothly during the "historical" period i.e. since the mid-nineteenth century, the period for which we have instrumental estimates of T (Fig. 1a), and is projected to continue to rise during the present century at a rate which depends on the emissions scenario.For example, the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change considered a set of scenarios under which the nominal radiative forcing at 2100 (relative to pre-industrial, regarded as the unperturbed steady state) ranges between 2.6 and 8.5 W m −2 (e.g.Fig. 12.4 of Collins et al. 2013).
When integrated with historical changes in radiative forcing agents, coupled atmosphere-ocean general circulation models (AOGCMs) show an ensemble-mean historical warming trend due to anthropogenic forcing that is very similar to the observed, as many studies have demonstrated (recently assessed by Bindoff et al. 2013).For example, in Fig. 2a we compare the ensemble-mean T from Abstract In both the observational record and atmosphere-ocean general circulation model (AOGCM) simulations of the last ∼150 years, short-lived negative radiative forcing due to volcanic aerosol, following explosive eruptions, causes sudden global-mean cooling of up to ∼0.3 K.This is about five times smaller than expected from the transient climate response parameter (TCRP, K of globalmean surface air temperature change per W m −2 of radiative forcing increase) evaluated under atmospheric CO 2 concentration increasing at 1 % yr −1 .Using the step model (Good et al. in Geophys Res Lett 38:L01703, 2011.doi:1 0.1029/2010GL045208), we confirm the previous finding (Held et al. in J Clim 23:2418-2427, 2010. doi:10.1175/2009JCLI3466.1)that the main reason for the discrepancy is the damping of the response to short-lived forcing by the thermal inertia of the upper ocean.Although the step model includes this effect, it still overestimates the volcanic cooling simulated by AOGCMs by about 60 %.We show that this remaining discrepancy can be explained by the magnitude of the volcanic forcing, which may be smaller in AOGCMs (by 30 % for the HadCM3 AOGCM) than in offline calculations that do not account for rapid cloud adjustment, and the climate sensitivity parameter, which may be smaller than for increasing CO 2 (40 % smaller than for 4 × CO 2 in HadCM3).
the "historical" simulations of 16 AOGCMs of the Coupled Model Intercomparison Project Phase 5 (CMIP5, Table 1, black line) with an observational estimate (HadCRUT4, green line, Morice et al. 2012).We use one integration of each AOGCM, and the historical T in each is the difference from its parallel control experiment with constant preindustrial atmospheric composition.
The historical simulations also include natural forcing, due to variability of solar irradiance and to aerosol injected into the stratosphere by explosive volcanic eruptions (henceforth referred to as "volcanic aerosol").For a few years following the eruption, volcanic aerosol causes a net negative radiative forcing ("volcanic forcing") of the climate system, by reflection of sunlight (shortwave radiation), partly offset by absorption of outgoing longwave radiation by the volcanic aerosol (Oman et al. 2005;Forster and Taylor 2006).The timeseries of volcanic forcing for the historical period from the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (Fig. 1a; Myhre et al. 2013) shows a maximum magnitude of −3.6 W m −2 for the Krakatau eruption of 1883.The Pinatubo eruption of 1992 was the next largest, and there have been no such large events since.Other estimates of time-dependent volcanic F(t) are similar (e.g.Figs. 2, 5 of Forster et al. 2013).
The CMIP5 ensemble-mean T shows a sudden globalmean cooling of 0.1-0.3K caused by the negative forcing from each historical volcanic eruption (Fig. 2a).The observed T timeseries has larger interannual variability, because unforced interannual variability is independent in    each model; the interannual standard deviation of T would be reduced by a factor of √ 16 in the mean if it were the same in all the individual models.Against this background, it is not simple to evaluate the cooling due to each volcano in reality.A thorough evaluation of the AOGCMs compared with observations involves removing the influence of and considering the relationship with unforced modes of variability such as El Niño (Thompson et al. 2009;Driscoll et al. 2012;Ding et al. 2014;Maher et al. 2015).In this work, our interest is the factors affecting the magnitude of cooling as simulated by the AOGCMs.
With time-dependent forcing F(t) that increases at a roughly constant rate, experiments with AOGCMs show that F(t) = ρT (t) is a fairly good approximation on time- scales from about 10 years to several decades (Raper et al. 2002;Gregory and Forster 2008;Gregory et al. 2015), where the climate resistance ρ (W m −2 K −1 ) is a modeldependent property of the climate system.The physical interpretation of this simple model is that the radiative forcing F is balanced by heat loss αT to space and heat uptake κT by the ocean, which holds the great majority of the heat capacity of the climate system (Levitus et al. 2001;Church et al. 2013).The climate feedback parameter α and the ocean heat uptake efficiency κ are both positive, and ρ = α + κ.In this picture, T is a surface skin temperature, with negligible thermal inertia, determined by the Earth energy balance F = N + αT, where N is the net down- ward radiative heat flux at the top of the atmosphere, and N = κT if we neglect heat storage other than in the ocean.We call F = ρT the "zero-layer model" (Bouttes et al. 2013), because it does not include any finite heat capacity.
We call 1/ρ the "transient climate response parameter" (TCRP, K W −1 m 2 , Gregory et al. 2015).It is the increase in global-mean temperature per unit increase in radiative forcing during time-dependent climate change.The standard benchmark for the predicted AOGCM response to anthropogenic forcing is the transient climate response (TCR, Cubasch et al. 2001), evaluated under the idealised 1pctCO2 scenario, in which the atmospheric CO 2 concentration increases at 1 % yr −1 .The TCR is defined as T after 70 years, the time of 2 × CO 2 i.e. double the initial concentration.The TCR and TCRP are related by TCR = F 2× × TCRP = F 2× /ρ, where F 2× is the radia- tive forcing of 2 × CO 2 . (We use use the term "TCRP" rather than "transient climate sensitivity", which has also been suggested e.g.Held et al. 2010, to avoid confusion with the equilibrium climate sensitivity, and following an analogy with the relationship between the equilibrium climate sensitivity in K and the climate sensitivity parameter in K W −1 m 2 .)Gregory and Forster (2008) and Held et al. (2010) pointed out that in observations and simulations of historical climate change the response of T to volcanic forcing is much smaller than would be expected from ρ calculated for CO 2 forcing from idealised scenarios of CO 2 increase.We demonstrate this by comparing the AOGCM ensemblemean historical T with the ensemble mean of estimates derived from the AR5 historical forcing according to the zero-layer model F = ρT (Fig. 2a, black and blue lines) using ρ for each model from its own 1pctCO2 experiment (Table 1).The time-profile of anthropogenic warming is reasonably well-reproduced by the zero-layer model, with somewhat overestimated magnitude-this could be because the TCRP is larger at the higher CO 2 concentration at which it is evaluated under 1pctCO2 (Gregory et al. 2015).However, the short-lived pronounced volcanic cooling is hugely exaggerated by the zero-layer model.
As a simple measure of this effect, for each of the six major eruptions named in Fig. 2, we compute the "volcanic cooling" ∆T = min(T (t), T (t + 1))− −mean(T (t − 2), T (t − 1)), where t is the year of the erup- tion, for both the zero-layer and the AOGCM T. This quantity measures the maximum cooling caused by the volcanic For all of them a historical experiment for 1860-2010 was available, with time-dependent volcanic aerosol, solar irradiance, and anthropogenic changes in atmospheric composition.The column headed "Nat?" indicates the AOGCMs (marked "Y") for which a historical-Nat experiment for 1860-2010 was used, with natural forcings only (volcanic and solar, not anthropogenic).The climate feedback parameter α, ocean heat uptake efficiency κ and climate resistance κ are in W m −2 K −1 and evaluated from the time-means of years 61-80 in 1pctCO2 experiments using the 4 × CO 2 forcing derived from years 1-20 of abrupt4xCO2 experiments (the same as in Table 1 of Gregory et al. 2015) forcing with respect to the years immediately before.
According to a one-parameter regression (requiring an intercept of zero) of the ensemble-mean ∆T from the zero-layer model against the AOGCM ensemble-mean ∆T, the zero- layer model overestimates the cooling by a factor of about five.Held et al. (2010, their Fig. 3) show that the overestimate by the zero-layer model can be explained by its neglect of the heat capacity C u of the upper ocean, which is important for episodic forcing (like volcanoes), but not for gradual multidecadal forcing change (like 1pctCO2).They demonstrate this by using the one-layer model (so called by Geoffroy et al. 2013, and"upper-layer model" by Gregory et al. 2015) where γ T is the rate of heat loss from the upper ocean to the deep ocean beneath, which is treated as an infinite heatsink.In the limit C u → 0, the upper-layer model becomes the zero-layer model F = ρT, if we identify γ with κ so that ρ = α + γ.
If a forcing F is imposed instantaneously at t = 0 and held constant, T = (F/ρ)(1 − e −ρt/C u ) in the upper-layer model.From CMIP5 abrupt4xCO2 experiments, in which CO 2 is quadrupled at the start and subsequently held constant, Geoffroy et al. (2013) found that C u has the heat (1) capacity of a few tens of metres of water, and the response timescale C u /ρ is 4.1 ± 1.0 years (τ f in their Table 4); Held et al. (2010) assume 4 years.For times which are much shorter than this, C u dT /dt = Fe −ρt/C u ≃ F and ρT ≪ F. i.e. most of F is absorbed by the upper-ocean heat capacity (of course volcanic forcing is negative, and "absorption" means loss of heat in this case) rather than by heat loss to the deep ocean or through climate feedback.For volcanic forcing which is large in magnitude for only a year or two, the TCRP and the zero-layer model are therefore not applicable.
Considering the Earth energy balance F = N + αT, we see that the prediction of T in response to volcanic forcing is affected by uncertainties in the forcing F and the feedback α, as well as in the ocean heat uptake N. Two issues in particular have been identified in previous studies of the AOGCM simulation of volcanic response (e.g.Wigley et al. 2005;Boer et al. 2007;Bender et al. 2010), namely that volcanic F is not exactly known in AOGCMs because it is not usually diagnosed, and that α for volcanic forc- ing might not be the same as for CO 2 .As an alternative to analysis of the AOGCM experiments, in this work we estimate the AOGCM response to episodic volcanic forcing using the "step" model (described in the next section), with which we compare the influences on the AOGCM response of N, F and α in Sects.2, 3 and 4 respectively, In Sect. 5 we draw conclusions about the TCRP as applied to volcanic forcing, including comparison with the results of Merlis et al. (2014).

Ocean heat uptake
The "step" model (Good et al. 2011) has the advantage that it avoids fitting any parametric form to the AOGCM results, unlike the upper-layer model and the upwelling-diffusion model of Wigley et al. (2005).Instead, it uses the response of the AOGCM itself to constant forcing instantaneously imposed (i.e. a "step-change").The step model relies on the assumption of linear systems theory (Good et al. 2015) that the response of the system depends linearly on forcing, so that the response to a sum of forcings equals the sum of responses to individual forcings.Therefore the response X(t) of any climate variable to a forcing scenario F(t) can be estimated as the sum of responses to a series of t annual forcing increments F(t) − F(t − 1), with F(0) = 0, accord- ing to where the sum is over years, and X s (t) is the response after t years to a constant forcing F s imposed at t = 0.
(2) For the sake of argument, let us idealise the forcing due to a single volcanic eruption as a step of F(<0) in year t = 1 and a step of −F in year t = 2 back to the initial level i.e. a pulse lasting for a year.According to Eq. 2, the response to the pulse with respect to the mean state is (only these two terms in the sum are non-zero).Although we do not use this equation in a form with continuous time, it is interesting to note that for small δt and in the limit δt → 0 the response to a delta-function pulse is the time-derivative of the step response.
Since X s (0) = 0 (before the forcing is switched on), X(1) = (F/F s )X s (1).Hence N in the year of the erup- tion is simply the first-year N s (1) in response to the step forcing, scaled by the ratio F/F s of the forcings, and N(1)/F = N s (1)/F s .In the ensemble mean of CMIP5 abrupt4xCO2 experiments for the models of Table 1, N s (1)/F s = 0.83 (estimated from the linear fit shown in Fig. 5 of Gregory et al. 2015).Thus the step model agrees with the upper-layer model in predicting N ≃ F during a short volcanic eruption.That is, a much larger proportion of (3) F is absorbed by the ocean than in the case of the gradually increasing 1pctCO2 forcing, for which N/F = κ/ρ ≃ 0.36 in the model mean at the time of 2 × CO 2 (Table 1).Running the step model with each AOGCM's abrupt4xCO2 T s (t) and the AR5 forcing timeseries (Fig. 1a), we estimate the ensemble-mean T(t) for the CMIP5 historical and "historicalNat" experiments (solid red lines in Fig. 2); the latter has the same natural forcing as the former, but no anthropogenic forcing.As Held et al. (2010) found with the upper-layer model, the volcanic cooling in the step model is much smaller than given by constant TCRP.However, it is still somewhat larger than in the AOGCMs.Computing ∆T as defined in Sect. 1 for the historical ensemble, we find that the step model overestimates the volcanic cooling by about 60 % on average.Moreover, the step-model simulations for historicalNat show a long-term negative trend not present in the AOGCM ensemble mean.
Since the CMIP5 historical experiments include solar forcing as well as volcanic, we have carried out a "hist-Vol" experiment, comprising an ensemble of four integrations with time-dependent historical volcanic aerosol 1860-2010 and no other forcing agents using the HadCM3 AOGCM (Gordon et al. 2000).This model was included in CMIP3 (the previous generation of the Coupled Model Intercomparison Project), performs well in comparison to many more recently developed AOGCMs (Reichler and Kim 2008) and is computationally relatively inexpensive by today's standards.(The experiments referred to in this paper are listed in Table 2.) Given ρ = 1.6 W m −2 K −1 for HadCM3 under CO 2 forcing (Gregory and Forster 2008), the zero-layer model predicts a maximum cooling of about 2 K during the largest events, whereas HadCM3 cools by 0.3 K at most relative to its control (Fig. 3).Using the HadCM3 abrupt4xCO2 T s (t) and the AR5 volcanic F(t), the step model gives a maximum cooling of about 0.6 K, and a long-term negative trend (Fig. 3), qualitatively similar to the CMIP5 historicalNat ensemble.
In summary, the results of the step model for CMIP5 and HadCM3 confirm that the zero-layer model's neglect of upper-ocean heat capacity is the main reason for the excessive cooling it predicts in response to volcanic forcing, but indicate that it is not the whole explanation.
3 Volcanic forcing CMIP5 does not include experiments which can be used to diagnose the volcanic forcing simulated by AOGCMs, unlike for CO 2 , and the lack of knowledge of F is an obstacle to analysis of the T response (Bender et al. 2010).Therefore we run a pair of experiments (sstPi and sst-PiHistVol, Table 2) to diagnose the historical volcanic Since the sea surface conditions are the same in the two experiments, global-mean T is almost the same (0.002 K cooler in sstPiHistVol; the effect of land temperature change is very small).Hence αT cancels out in the differ- ence in the energy balance F = N + αT between the two experiments, and F(t) is diagnosed as the difference in N(t) (Hansen et al. 2005;Held et al. 2010

; Andrews 2014).
There is statistical uncertainty in this estimate of F due to unforced variability in N, whose interannual standard deviation is 0.17 W m −2 in the HadCM3 atmosphere model with constant boundary conditions.Consequently we regard forcing with magnitude |F| < 1.65 × 0.17/ √ 2 = 0.22 W m −2 as insignificant (at the 10 % level) in the ensemble mean of the two integrations of sstPiHistVol.
The HadCM3 volcanic forcing (Fig. 1a) is positive in years when there is no volcanic aerosol (the majority of years, between major eruptions), because there is a nonzero concentration in the control.The alternative assumption, made in the AR5, that zero volcanic aerosol implies zero forcing, is incorrect for the real world, because the unperturbed natural state of the climate system includes the effects of occasional eruptions.The permanent cessation of volcanic eruptions would produce a climatic warming, so must imply a positive forcing (Gregory 2010;Gregory et al. 2013).In the long term, the small positive forcing in years with no volcanic aerosol is balanced by the small number of years with large negative forcing.On the other hand, if an AOGCM control integration does not include volcanic forcing, F = 0 for zero volcanic aerosol, and the negative time-mean volcanic forcing in the historical experiment produces a negative trend in ocean heat content.
Apart from this offset, the HadCM3 volcanic forcing timeseries is similar to the AR5 (Fig. 1a) and they are strongly correlated (0.95), as expected because they were derived by different methods from the same volcanic aerosol timeseries (Sato et al. 1993, updated).However, regression of annual means shows that the HadCM3 forcing is only 77 % of the AR5 (Fig. 1b).Because global-mean T is not exactly zero in sstPiHistVol (see above), the magnitude of F will be underestimated, but this effect is only 3 % (not shown).Regression of F against global-mean volcanic aerosol optical depth (AOD) gives a slope of −24.6 ± 0.2 for AR5 (blue line in Fig. 4), consistent with the AR5 formula of −25 W m −2 per unit AOD (Table 8.SM.8 in the supplementary material of Myhre et al. 2013), following the results of Hansen et al. (2005, their equation 2a) from the GISS AOGCM for Pinatubo.For HadCM3, the slope is only −19.0 ± 0.5 W m −2 (solid black line).
We have examined this relationship also in a pair of experiments carried out by Andrews (2014)  Atmosphere Intercomparison Model Project Phase II (Taylor et al. 2000;Hurrell et al. 2008).One experiment includes all forcings agents, both anthropogenic and natural, and the other anthropogenic only.By the same argument as for sstPiHistVol, the difference in N between these experiments is the natural radiative forcing.To obtain the volcanic contribution, we subtract the AR5 estimate of solar forcing.The AMIP period includes two major volcanic eruptions (El Chichon and Pinatubo), and the regression slope of volcanic F against AOD is −17.0 ± 1.0 W m −2 (red line in Fig. 4), weaker than in HadCM3 and further from the AR5 formula.
We propose that the difference from the AR5 formula is caused by rapid tropospheric adjustment (Myhre et al. 2013;Sherwood et al. 2015) in our AGCMs in response to volcanic aerosol.By including such adjustment, the AGCM diagnosis gives an effective radiative forcing.The AR5 formula is based on results from a version of the GISS AOGCM (Hansen et al. 2005) in which tropospheric adjustment to volcanic aerosol was apparently smaller.We apply the method of approximate partial radiative perturbation (APRP, Taylor et al. 2007) to estimate the shortwave effect of rapid cloud adjustment in HadCM3-A.It is a positive quantity, and arises from a reduction of cloud fraction and the planetary albedo when volcanic aerosol is imposed.With this positive adjustment subtracted, the HadCM3 forcing has a regression slope of −26.6 ± 0.5 W m −2 against the AOD, much closer to the AR5 value (Fig. 4).Thus shortwave cloud adjustment would be sufficient to explain the difference, but there may also be longwave cloud adjustment, which APRP cannot be used to evaluate.
Rapid adjustment in shortwave cloud radiative effect (CRE) has been previously been noted in response to CO 2 forcing in AGCMs (Gregory and Webb 2008) and AOGCMs (Andrews et al. 2012;Zelinka et al. 2013), due to reduction in cloudiness.In HadCM3, there is a positive shortwave cloud adjustment of 1.7 W m −2 included in the net 4 × CO 2 forcing of 7.7 W m −2 .It is interesting and worthy of further investigation that both negative shortwave volcanic and positive longwave CO 2 forcing produce a positive shortwave cloud forcing adjustment.
In view of the differences between the AR5 and sstPi-HistVol F(t), we rerun the step model for HadCM3 hist-Vol.First, we use the AR5 F(t) adjusted to have the same small time-mean (−0.04 W m −2 ) as the HadCM3 sstPi-HistVol F(t); this adjustment removes the long-term cooling trend (Fig. 3, dotted red line), which is due to the more negative time-mean (−0.24W m −2 ) of the AR5 F(t).In fact before the Krakatau eruption there is now a warming trend, because zero volcanic aerosol implies positive forcing.Second, using the sstPiHistVol F(t) instead of the AR5, the step-model estimate of the HadCM3 volcanic cooling becomes smaller in magnitude (dashed red line), and thus more similar to the AOGCM.
Rerunning the step model for the CMIP5 historicalNat ensemble using the HadCM3-A sstPiHistVol F(t) removes the long-term cooling trend and reduces the magnitude of the volcanic cooling (Fig. 2b), just as in HadCM3.This suggests that rapid adjustments may occur in CMIP5 AOGCMs that reduce the magnitude of volcanic forcing, and that this set of CMIP5 AOGCMs is near to steady state with time-mean historical volcanic forcing (like HadCM3), rather than with zero volcanic forcing (like CMIP3 models; Gregory 2010).
We found above that N s /F = 0.83 in the CMIP5 abrupt4xCO2 ensemble in the first year of forcing.In the historicalNat ensemble mean, for years with significant volcanic forcing (F < −0.22 W m −2 in the AR5 F time- series, the threshold obtained from HadCM3, Sect.3), the slope of the regression of N(t) against HadCM3 volcanic F(t) is 0.89 ± 0.08 (red line in Fig. 5), consistent with the expected value.This is further, although circumstantial, evidence that the volcanic forcing in the CMIP5 ensemble mean, as in HadCM3, is less than the AR5 formula.As similar conclusion has been reached by Larson and Portmann (2016), who use an inverse application of the step model to estimate AOGCM forcings.Their method assumes that the same climate feedback applies to volcanic forcing and CO 2 , which is the issue we examine in the next section.

Climate feedback
The magnitude of cooling due to volcanic forcing is affected by the climate feedback parameter according to T = (F − N)/α.Therefore a third possible contribution to the smaller cooling than expected from the TCRP is that α for volcanic forcing is larger than for CO 2 forcing.The assumption that α is the same for both kinds of forcing is a motivation for evaluating climate feedbacks from observations and simulations of the response to the Pinatubo and other large eruptions (Soden et al. 2002;Forster and Collins 2004;Wigley et al. 2005;Merlis et al. 2014).There is evidence in support of the assumption, but there are uncertainties arising from forcing and unforced variability (Bender et al. 2010).
To evaluate α for volcanic forcing in HadCM3, we carried out a ten-year "step volcano" experiment, called "abruptPin" (analogous to abrupt4xCO2, Table 2), with constant volcanic aerosol as it was in January 1992, approximately the peak of the Pinatubo eruption.An advantage of this method is that we do not need to know the forcing in advance; since F = N − αT, regression of annual-mean N versus T gives both the forcing as the N-intercept and −α as the slope (Gregory et al. 2004).The N-T relationship is not quite linear (Fig. 6), which indicates that α is not constant, as has been found under constant 4 × CO 2 in most AOGCMs (Winton et al. 2010;Andrews et al. 2012;Gregory et al. 2015).Unlike CO 2 , volcanic forcing is short-lived, and the feedbacks in response to a "step volcano" sustained for many years are not relevant to historical eruptions.We therefore use only the first five years in the regression, although the results are not greatly different for ten years, since non-linearity is not pronounced on this timescale.The experiment comprised an ensemble of 14, to obtain adequate signal/noise.
Before considering the feedback, we check that the forcing is consistent with the results of Sect. 3. From the regression, F = −2.7 ± 0.1 W m −2 (Fig. 6b, cf.+7.7 W m −2 for abrupt4xCO2 in HadCM3).The global-mean AOD is 0.148, so the forcing is −18.4 ± 0.1 W m −2 per unit AOD, nearly the same as in sstPiHistVol.We presume it is slightly different because the geographical distribution of volcanic aerosol is variable in the historical timeseries and has a small effect on the forcing.APRP shows that there is a positive shortwave cloud adjustment of 0.9 ± 0.1 W m −2 to the forcing (from the intercept of this term in Fig. 6b), whereas the forcing adjustment from surface shortwave absorption is negligible.Excluding the adjustment increases the forcing to −25 W m −2 per unit AOD, in agreement with the AR5 formula (as in Sect.3).
We also evaluated the forcing of abruptPin following the method of sstPiHistVol (Sect.3), as the perturbation to N in experiment sstPiPin (Table 2), which has the constant volcanic aerosol of abruptPin in HadCM3-A with control sea surface conditions.By this method we obtain the same value of forcing (to one decimal place) as from abruptPin, and the value lies on the relationship between forcing and AOD found in sstPiHistVol (Fig. 4).Although in different ways, both methods estimate F as N under the influence of volcanic aerosol and in the absence of climate change; they are therefore expected nearly to agree, and have previously been found to do so (e.g.Hansen et al. 2005;Andrews et al. 2012).
The slope of the regression for abruptPin (Fig. 6) gives α = 1.81 ± 0.22 W m −2 K −1 , significantly larger (by about 40 %) than its value of 1.25 ± 0.04 W m −2 K −1 from the first 20 years of abrupt4xCO2 in HadCM3 (Andrews et al. 2015).(Non-linearity is statistically insignificant during these years; regression using only the first five gives α = 1.30 ± 0.15 W m −2 K −1 .)We note that α in abruptPin is consistent with α = 1.8 W m −2 K −1 esti- mated for infinitesimal changes in CO 2 relative to the control; in HadCM3, α is larger for smaller step-increases in CO 2 concentration (Good et al. 2011(Good et al. , 2012;;Gregory et al. 2015).
The ratio N s (1)/F s = −2.16± 0.12 ÷ −2.71 ± 0.13 = 0.80 ± 0.06 in abruptPin (shown in Fig. 6b, standard errors) agrees with the regression slope of N versus F in histVol (0.82 ± 0.10 , blue line in Fig. 5), in accordance with the step model (Sect.2).For abrupt4xCO2 F s in HadCM3, the ratio N s (1)/F s = 6.38 ± 0.16 ÷ 7.73 ± 0.14 = 0.83 ± 0.03 i.e. the fraction of forcing taken up is slightly greater but not significantly different for positive forcing.This is contrary to expectation that negative forcing would cause the ocean to lose heat more readily, by destabilising the ."SW" is shortwave, "LW" longwave, "CRE" cloud radiative effect, "APRP" approximate partial radiative perturbation.In (a) the clear-sky effect is the change in the radiative flux diagnosed with clouds ignored (referred to as Method II by Cess et al. 1993) and the CRE is the difference between the all-sky (i.e. with the cloud fraction simulated by the GCM) and the clear-sky fluxes.In (b) the APRP method is used to decompose the change in SW radiative flux under all-sky conditions into contributions from changes in cloud-free air, cloud and surface albedo vertical temperature profile (e.g.Stouffer and Manabe 1999;Bouttes et al. 2013;Merlis et al. 2014).Clearly this is not an important effect in HadCM3 for this magnitude of forcing.
We expect that the HadCM3 T(t) in histVol should be consistent with the feedbacks shown in abruptPin.To test this, we use abruptPin (instead of abrupt4xCO2) for T s (t) with F(t) from HadCM3-A sstPiHistVol in the step model.Since abruptPin was only 10 years long, we set T s (t) = T s (10) in Eq. 2 for all t > 10, which means that the cooling caused by a volcanic eruption abruptly vanishes 10 years later.By this time it is actually quite small anyway, because in response to any constant forcing the differences |T s (t) − T s (t − 1)| diminish with time (Fig. 6), so the response to a pulse forcing (Eq. 3) likewise decreases.This is consistent with our expectation that the influence of a brief forcing will gradually be forgotten as time passes and the system returns asymptotically to its unperturbed state.
With this combination of inputs, the step model gives the closest of its estimates to the HadCM3 histVol T(t) (blue line in Fig. 3).This confirms that the smaller climate sensitivity to volcanic forcing during the historical period than for elevated CO 2 is also a reason for the overestimation of volcanic cooling by the TCRP in HadCM3, although less important than ocean heat uptake and volcanic forcing.We are not able to test this possibility for other AOGCMs, since similar experiments to abruptPin have not been carried out.

Conclusions
The zero-layer model of time-dependent climate change T = F/ρ gives a fairly accurate reproduction of changes in global-mean surface air temperature T(t) as observed during the CMIP5 "historical" period (since the latter part of the 19th century) and as simulated by CMIP5 AOGCMs in response to the smoothly varying anthropogenic part of the forcing F(t).The transient climate response parameter used in the zero-layer model (TCRP, the increase in T per unit increase in F, in K W −1 m 2 ) is evaluated from idealised climate-change experiments with CO 2 increasing at 1 % yr −1 .The TCRP is the reciprocal of the climate resistance ρ (W m −2 K −1 ), which is the sum of the climate feed- back parameter α and the ocean heat uptake efficiency κ.
The zero-layer model ovestimates, by a factor of about five, the sudden cooling in the AOGCMs caused by the short-lived large negative forcing from volcanic aerosol following explosive volcanic eruptions.This could be due to errors in any of the factors involved, namely F, α and κ.With reference to the two-layer model of the Earth energy balance (Gregory 2000;Held et al. 2010;Geoffroy et al. 2013;Gregory et al. 2015).Held et al. (2010) attributed the overestimate to the zero layer model's neglect of the relatively small upper-layer heat capacity, which is unimportant for warming on multidecadal timescales forced by increasing CO 2 , but dominates the response of T to an impulsive forcing.Consequently κ does not correctly account for the ocean heat loss.
We confirm that most of the error is due to this effect by using the step model (Good et al. 2011) to estimate the AOGCM response to historical volcanic forcing from the AOGCM response to a step-change in CO 2 forcing.In the step model and the AOGCMs, ocean heat content change is more than 80 % of the volcanic forcing during the first year, because the upper ocean readily gives up heat, whereas heat uptake is less than 40 % of CO 2 forcing gradually increasing over decades, because it is limited by the less efficient thermal coupling to the deeper ocean.The difference is not due to the sign of the forcing (negative volcanic, positive CO 2 ), but its timescale.However, although much closer than the zero-layer model, the step model ovestimates the cooling simulated by AOGCMs by about 60 %, and we explain this remaining discrepancy in terms of forcing and feedback.
The AR5 formula (−25 W m −2 per unit AOD) overestimates the magnitude of the volcanic forcing in HadCM3 (−19 W m −2 per unit AOD) by about 30 %, and in HadGEM2 (−17 W m −2 per unit AOD) by about 50 %.We have shown that this can be explained in HadCM3 by a rapid positive shortwave cloud adjustment which reduces the magnitude of the negative volcanic forcing.This may be a model-specific result, but the step-model simulations suggest that the CMIP5 ensemble-mean volcanic forcing is more similar to that of HadCM3 than to the AR5 estimate, which was based on the GISS AOGCM (Hansen et al. 2005).
We find that the climate feedback parameter α is about 40 % greater (climate sensitivity parameter smaller) for volcanic forcing than for 4 × CO 2 in HadCM3.This could be related to the different natures of the forcing, or it might be that α depends on CO 2 concentration (Jonko et al. 2012;Meraner et al. 2013;Caballero and Huber 2013;Gregory et al. 2015), for which the value applicable in the historical period is nearer to the pre-industrial CO 2 concentration than to 4 × CO 2 .Because volcanic perturbations to climate in the historical period are short-lived and comparable in magnitude to unforced interannual variability, a large ensemble is needed to evaluate α with this precision, and consequently we cannot determine whether the CMIP5 models also show this effect.
It would be useful for investigation of volcanic forcing and feedback in CMIP6 if ensemble experiments were conducted with historical volcanic aerosol as the only forcing agent in each AOGCM, to diagnose the climate response, and with the corresponding AGCMs with prescribed sea surface conditions, to diagnose the radiative forcing.This would reveal whether other models also exhibit a cloud adjustment and a lower climate sensitivity for volcanic forcing, and offer the opportunity for analysis of the processes involved.Overestimated historical volcanic forcing, and overestimated climate sensitivity to such forcing, are possible explanations for the need to scale down the forcing in simple climate models in order to reproduce AOGCM results for volcanic cooling (Meinshausen et al. 2011;Lewis and Curry 2015).
An alternative to the use of a large ensemble to investigate the climate response to volcanic eruptions of the size that is typical of the historical period might be to improve the signal/noise by multiplying historical volcanic aerosol in an AOGCM by a large factor.However, this procedure may give inappropriate values for the feedback and forcing.For example, Jones et al. (2005) simulated a "supervolcano" of roughly the size of the Toba eruption 72 ka ago, which was two orders of magnitude greater than Pinatubo, using HadCM3 with Pinatubo volcanic aerosol multiplied by 100.The peak forcing was about −60 W m −2 , only 20 times greater than Pinatubo, and the climate feedback parameter α ≃ 4 W m −2 K −1 , more than twice the value (half the sensitivity) that we found for Pinatubo in HadCM3.The mechanisms for this are not a subject of the present work, but we note that α for volcanic aerosol of ten times the magnitude of Pinatubo exceeds 2 W m −2 K −1 in HadGEM2 and MPI-ESM1.1 (a revised version of MPI-ESM of Giorgetta et al. 2013) as well.
If α for volcanic forcing is larger than for CO 2 (its effi- cacy is less than unity, in the terms of Hansen et al. 2005), the method of Forster and Taylor (2006), used to evaluate the forcings in historical CMIP5 experiments by Forster et al. (2013), will underestimate the magnitude of the volcanic forcing, because it assumes that α is the same for all forcing agents.If α is not the same, the response to histori- cal volcanic eruptions cannot be used to estimate the effective or equilibrium climate sensitivity that applies to future CO 2 -forced climate change.Merlis et al. (2014) considered the possibility of using the climate response to volcanic (Pinatubo) forcing to place constraints on the TCRP for CO 2 increase in the GFDL-CM2.1 AOGCM.As we have seen, the volcanic cooling itself is not correctly predicted by the TCRP, but they follow an alternative and novel approach in their analysis by using the upper-layer model (C u dT /dt = F − ρT, Eq. 1) in a time-integral form where the time-integral runs from t = 0 when the vol- canic forcing begins, to t = τ, sufficiently long after the eruption that the recovery of T from the volcanic cooling is complete (they assume 15 years).Because the upper ocean heat content is the same at the end as at the beginning, its heat capacity C u is irrelevant in this integral form, which therefore agrees with the zero-layer model F = ρT.
The duration of the volcanic forcing is small compared with τ, so a reasonably accurate picture is that the forc- ing and the consequent sudden drop in T take place almost instantaneously at the start, and for most of the time-integration F is zero, while T is recovering back towards zero.This means that the time-integral method mainly measures the climate feedback parameter and the TCRP which apply to the relaxation towards the steady state following the perturbation, rather than during the perturbation itself, which happens very quickly.Merlis et al. (2014) find that the TCRP from the time-integral method is 5-15 % smaller than the TCRP for gradual CO 2 -forced warming in their AOGCM.It could be that this difference is partly due to a larger α in response to volcanic forcing, as we have found.
In summary, we conclude that the zero-layer model and the TCRP are not applicable to the rapid cooling due to volcanic forcing, mostly because of the importance of the upper-ocean heat capacity on short timescales.However, simple models may also overestimate the the volcanic cooling simulated by AOGCMs because the climate sensitivity parameter is smaller for volcanic forcing, and the volcanic effective radiative forcing is reduced by rapid adjustment.

Fig. 1 a
Fig. 1 a Timeseries of historical annual-mean radiative forcing F(t) assessed by Myhre et al. (2013) (AR5, Sect. 1) and for volcanic aerosol alone diagnosed from the HadCM3-A sstPiHistVol experiment (Sect.3).The first years of the six named major volcanic eruptions are indicated by the vertical lines.b Comparison of the AR5 estimate of historical volcanic forcing with the HadCM3-A sstPiHistVol estimate (Sect.3)

Fig. 2
Fig. 2 Global-mean surface air temperature T simulated by the ensemble mean of CMIP5 AOGCMs, compared with the ensemble means of estimates made from the historical forcing timeseries of Myhre et al. (2013) using constant TCRP (1/ρ) (Sect.1) and using the step model (Sect.2), (a) historical experiment (anthropogenic and natural forcings), relative to the time-mean of 1961-1990, also showing the HadCRUT4 observational T (Morice et al. 2012), (b) histori-calNat experiment (natural forcings only), relative to control, also showing results with forcing diagnosed from the HadCM3-A sstPi-HistVol experiment (Sect.3)

Fig. 4 Fig. 5
Fig. 4 Relationship between volcanic aerosol optical depth and volcanic forcing assessed by Myhre et al. (2013) (AR5) and diagnosed from HadCM3-A experiments with historical time-dependent volcanic aerosol (sstPiHistVol, Sect.3), constant volcanic aerosol approximately as at the peak of the Pinatubo eruption (sstPiPin, Sect.4), and HadGEM2-A AMIP experiments (Sect.3).Both the AOD and the forcing are shown as differences from the state of the model in the control experiment

Fig. 6
Fig. 6 Relationship between means for years 1-10 (years 9 and 10 are too close to be distinguished in the plot) of the change in globalmean downward radiative fluxes at the top of the atmosphere and the change in global-mean surface air temperature T in the ensemble mean of HadCM3 abruptPin integrations.Time runs from right to left, because the climate is cooling down.The lines show linear regressions of the radiative fluxes against T for years 1-5, whose slopes give the components of climate feedback (W m −2 K −1 ) shown in the key.A positive parameter indicates a negative feedback on climate change, because αT opposes F if α > 0.The net forcing F and N for the first year are marked in blue in (b)."SW" is shortwave, "LW" longwave, "CRE" cloud radiative effect, "APRP" approximate partial radiative perturbation.In (a) the clear-sky effect is the change in the radiative flux diagnosed with clouds ignored (referred to as Method II byCess et al. 1993) and the CRE is the difference between the all-sky (i.e. with the cloud fraction simulated by the GCM) and the clear-sky fluxes.In (b) the APRP method is used to decompose the change in SW radiative flux under all-sky conditions into contributions from changes in cloud-free air, cloud and surface albedo

Table 1
CMIP5 AOGCMs used in this work Global-mean surface air temperature change T with respect to control simulated by the ensemble of HadCM3 histVol experiments (Sect.2), and estimated by the step model using various combinations of forcing and response (Sects.2, 3 and 4).For the dotted red line, the AR5 forcing was adjusted by adding a constant so that it had the same time-mean as the HadCM3 forcing.The HadCM3 ensemble mean is shown by the thick black line and the envelope (from maximum to minimum in each year) of the four integrations by the grey shading

Table 2
Experiments carried out with the HadCM3 AOGCM, indicating in which section of the paper each is first discussed and the number of integrations in cases of an ensemble A (the atmosphere general circulation model, or AGCM, component of HadCM3), prescribing constant climatological sea surface boundary conditions (temperature and sea-ice) from the coupled HadCM3 control experiment.Experiment sstPi has volcanic aerosol prescribed at the constant level of the HadCM3 control, which is typical of the long-term mean; F = 0 by definition if this is regarded as the unperturbed state.Experiment sst-PiHistVol has time-varying historical volcanic aerosol.