The sensitivity f th energy budget and hydrological cycle to CO 2 and solar forcing

Introduction Conclusions References


Introduction
Human activities primarily affect the climate system in two competing ways: greenhouse gases warm the atmosphere by partly absorbing longwave radiation, while aerosols have a predominantly cooling effect by scattering incoming shortwave radiation (Wild et al., 2008).Both forcing agents alter the energy budget of the Earth (Kiehl and Trenberth, 1997;Trenberth et al., 2009), which triggers responses through complex feedback mechanisms in order to reach a new equilibrium state.Among all these mechanisms, the ones modifying the processes leading to precipitation formation are of particular interest because human societies as well as ecosystems will likely be affected by changing precipitation patterns.Precipitation, and its energy equivalent, latent heat, are variables that belong to both the energy budget and hydrological cycle (e.g.Bosilovich et al., 2008), hence the need to analyze them jointly.
Recent studies have highlighted the dependency of precipitation changes on different emission scenarios, i.e. on the different forcing agents (Shiogama et al., 2010a;Lambert and Allen, 2009).It is widely accepted that global mean precipitation change per unit temperature change is more sensitive to changes in solar radiation than to changes in CO 2 concentrations (Allen and Ingram, 2002;Gillett et al., 2004;Andrews et al., 2009;Bala et al., 2010).This difference in precipitation response has been investigated from the perspective of the fast (weeks to months and caused by the forcing agent directly) and slow (years to centuries and caused by changes in surface temperature) responses (Lambert and Faull, 2007;Bala et al., 2008Bala et al., , 2010;;Andrews et al., 2010;Cao et al., 2011).These studies found that the fast or forcing-dependent response of precipitation is different for CO 2 and solar forcings but that the slow or feedback-dependent response does not depend strongly on the nature of the forcing agent.Increasing CO 2 concentrations lead to small but rapid increases in tropospheric temperatures while surface temperature remains unchanged initially.Atmospheric stability is increased, reducing convection and therefore precipitation in the first days to months (Andrews et al., 2010).In addition, Andrews et al. (2011) and Cao et al. (2012)  showed that on time scales of a few days and over land, increasing CO 2 concentrations lead to a reduction in plant transpiration, which also contributes to precipitation reduction.While increasing solar radiation leads to an increased absorption of shortwave radiation in the atmosphere, also inducing a weak reduction in precipitation on short time scales, this direct effect of the forcing agents is almost negligible compared to that of CO 2 (O'Gorman et al., 2012).Many studies have quantified the climate responses in simulations where the forcing is increased instantaneously (e.g.Bala et al., 2010).While much can be learned from those, there is also currently a need to understand transient climate change as this is what is occurring in the real world.The aim of this study is therefore to quantify the transient response of the energy budget and the hydrological cycle to different forcing agents, globally and zonally.A particular focus is given to the assumption of linear additivity of the forcing's responses, i.e. whether the response to individual forcings can be added to estimate the response to the combined forcing.Meehl et al. (2004) found this assumption to be valid for the global mean response to the twentieth-century forcings.
However, Meehl et al. (2003) found that the model's temperature response to solar forcing is amplified when it is combined with anthropogenic forcings, which they interpreted as a non-linear response of the climate system to solar forcing.Results presented in Jonko et al. (2012) further suggest that global mean surface temperature neither exactly doubles from 2 × CO 2 to 4 × CO 2 nor from 4 × CO 2 to 8 × CO 2 .Testing the linear additivity of the global mean temperature change is particularly relevant for impact studies based on pattern scaling.This technique relies on the fact that the scaled pattern of change of surface variables in any scenario is approximately constant over time (Santer and Wigley, 1990).The resulting main advantage is that Global or Regional Climate Models (GCMS or RCMs) do not need to actually simulate all scenarios since the information about the pattern of change is approximately the same for all scenarios (Giorgi, 2008).Pattern scaling has been widely used for temperature and precipitation (e.g.Mitchell and Hulme, 1999;Meehl et al., 2007;Watterson, 2008) but recently Shiogama et al. (2010a) showed that the technique is less reliable for precipitation patterns, which Figures are more sensitive to the nature of the forcing agent.Besides the assumption that the pattern of change remains the same under low and high forcings and under forcings caused by different agents, another assumption needs to be made about the scaling factor.Global mean temperature change is usually used as the scaling factor and is estimated with energy balance models (Ruosteenoja et al., 2007).However, in those models, the global mean temperature response does not depend on the nature of the forcing but only on its total radiative perturbation unless efficacies for individual forcings are individually prescribed.In the way forcing agents are implemented in these models, the assumption that the temperature response to forcings is linear is made except if the feedback parameter is parametrized as a function of the forcing.The assumption of linear additivity was tested to some extent in Mitchell (2003) and will be tested systematically in this study for a broad range of variables using transient solar and CO 2 forcing increase simulations in a GCM.
A description of the climate model experiments is provided in Sect. 2. Section 3 summarizes the results and discusses the two key aspects considered, the assumption of the linear additivity of the responses to different forcing agents or forcing magnitudes, and the differences in the energy budget and hydrological cycle responses depending on the forcing agent.The conclusions are presented in Sect. 4.

A set of idealized transient simulations is performed with the NCAR Community Climate
System Model version 3.5 (CCSM3.5)(Collins et al., 2006;Gent et al., 2010).The finite volume dynamical core of this fully coupled ocean and atmosphere model has a spatial resolution of 1.9 • in latitude and 2.5 • in longitude, with 26 levels in the vertical.to doubled values (from 355 ppm up to 710 ppm and labeled "2x" hereafter).The next scenario is a 2 % yr −1 increase in CO 2 (from 355 to 1420 ppm and labeled "4x" hereafter).Then, to investigate the response of the energy budget and hydrological cycle to changes in solar radiation, the solar constant is transiently increased relative to the control value to reach a radiative forcing that corresponds theoretically to a doubling or quadrupling of CO 2 .According to the best estimate provided by Myhre et al. (1998), the radiative forcing of a doubling in CO 2 concentrations would be about 3.7 W m −2 .
Correspondingly, in the solar forcing simulations the solar constant is increased by the best estimate radiative forcings ∆F corresponding to each CO 2 scenario (i.e.3.7 or 7.4 W m −2 ) multiplied by 4 (to take account of the geometry of the Earth) and divided by 0.7 (the planetary albedo being approximately 0.3).The two solar scenarios are labeled "37" and "74".Finally a scenario with a 1 % yr −1 increase in CO 2 and an increase in solar forcing that reaches a 3.7 W m −2 radiative forcing at the end of the simulation is run (labeled "372x").In addition, a 25 yr simulation starting in year 400 for each scenario is performed with instantaneous increase of the forcing agents (2 × CO 2 , 4 × CO 2 , 3.7 W m −2 , 7.4 W m −2 and 2 × CO 2 combined with 3.7 W m −2 ) to infer the adjusted forcing.The adjusted forcing is different form the radiative forcing in that it includes the rapid adjustments occuring within a few days in the troposphere and land-surface (Forster et al., 2013).The method by Gregory et al. (2004)  As mentioned before, the scenarios are designed to test the assumption of linear additivity in the response to solar and CO 2 forcings.If the assumption is valid, scaling the global mean response at the end of the simulation in any variable in the 2x scenario by the adjusted forcing in 2x should be equal to the scaled responses in the 4x scenario.Similarly, scaling the global mean responses in the 37 scenario by its adjusted forcing should give similar results as the scaled responses in the 74 scenario.Finally, adding the responses of the 2x and 37 scenario and scaling them by the sum of the adjusted forcing in 37 and 2x should be equal to the scaled responses in the 372x scenario.A two-sample Kolmogorov-Smirnov test with significance level 0.05 is used to assess whether the distributions of the five ensemble members of two scenarios (the scaled 2x and scaled 4x for example) are statistically significantly different.

Linear additivity of the responses
A widely-made assumption in climate science is that the responses to individual forcings are linearly additive.Techniques that rely on this assumption are for example detection and attribution (Barnett et al., 2005) or pattern scaling (Mitchell, 2003).Pattern scaling approaches are based on the fact that the response pattern in one scenario is very similar to the pattern in another scenario.For the idealized simulations performed in this study, all response patterns of variables from the energy budget and hydrological cycle are highly correlated.As an example, the response pattern of total precipitation averaged over the last 30 yr of each scenario has a correlation of ρ = 0.93 between 2x and 4x and of ρ = 0.91 between 37 and 74.Even between the different forcing agents, correlation coefficients are high (ρ = 0.94 between 2x and 37 and ρ = 0.92 between 4x and 74).However, pattern scaling techniques further make the assumption that the responses to different forcing agents or to forcing agents of different intensity add linearly.overestimation of the changes in precipitation with this technique.Here we test the linear additivity assumption for changes in global mean temperature and changes in variables of the energy budget and hydrological cycle in a GCM.The objective of this section is not to make any judgment on the validity of the mentioned techniques.Depending on the purpose, it can be then decided whether the errors introduced by the assumption of linear additivity are negligible or not.First, we use the instantaneous forcing increase simulations described in the previous section and the method proposed in Gregory et al. (2004) to estimate the adjusted forcing in each scenario.There are however other ways to estimate the radiative forcing, as described in Hansen et al. (2002) or Shine et al. (2003).The adjusted forcing is 3.77 W m −2 for the 2x scenario, 8.15 W m −2 for 4x, 2.81 W m −2 for 37, 5.61 W m −2 for 74 and 6.84 W m −2 for 372x.The adjusted forcing in the 2x scenario agrees well with the theoretical value from Myhre et al. (1998).However, the adjusted forcing in the 4x scenario is larger than expected, and the ratio of adjusted forcings between 4x and 2x is 2.16.For the solar scenarios, the ratio of adjusted forcings is indeed 2 but the adjusted forcings are smaller than expected.Hansen et al. (2005) and Schmidt et al. (2012) showed that the efficacy of solar forcing is smaller than that of CO 2 forcing.The net shortwave flux is indeed as prescribed, 3.7 and 7.4 W m −2 (see Fig. 2b), however the increase in net longwave flux, due to the black-body response of the warmer surface (see Fig. 2a), compensates the incoming solar energy, and the net energy flux at TOA is small.Changes in energy fluxes in the scenarios will be further discussed in the next section.The fact that the adjusted forcings do not add linearly in the CO 2 scenarios and that the adjusted forcings in solar simulations are lower than expected from Myhre et al. (1998) are important since such results cannot be obtained from simple energy balance models.In these models, all forcing agents are added into a total radiative forcing value, and by construction, there is no forcing-dependency of the temperature response (except if it is explicitly parametrized).
The linear additivity assumption is tested for variables of the energy budget and hydrological cycle and results are presented in  1 indicate that the linear additivity assumption is not valid for most variables and in general, the responses to high forcings are larger than expected.Differences on the order of 10 % might be considered acceptable for some applications given that model or observational uncertainties are potentially much larger.As an example, Mitchell (2003) concludes that pattern scaling is generally accurate with significant errors not exceeding 2.8 % of the global mean temperature change and 11 % of the global mean precipitation change.Results obtained with the CCSM3.5 model show much larger differences.One could argue that these errors are avoidable as computational capacity is becoming less of a limiting factor for any given model, but on the other hand models are getting more expensive, such that the wall-clock time spent for a single simulation has remained remarkably constant over the past decades.However, other sources of errors are less easy to get rid of, for example structural or observational uncertainties.
As mentioned above, the system is obviously not in equilibrium 30 yr after the transient forcing increase and this is one reason for the non-validity of the linear additivity assumption.Even though for global mean temperature and in simple models, the ratio of warming per unit forcing is roughly constant even for the transient case (Gregory and Forster, 2008;Knutti et al., 2008), that is unlikely to hold for other variables and local Introduction

Conclusions References
Tables Figures

Back Close
Full changes.Long model simulations show that it takes hundreds of years to reach equilibrium (Gregory et al., 2004;Stouffer, 2004).Even after 300 yr, scenarios with CO 2 forcing (see Fig. 1a) show that the system is not in equilibrium as surface temperature continues to rise.However, there is currently a need for assessments of transient climate change since this is what is happening in the real world, the climate system being never (for long, at least) in equilibrium.In the recently published Representative Concentration Pathways (RCP) scenarios, stabilization of the forcings do not occur before 2100 and some pathways are even bell-shaped, with stabilization occuring after 300 yr (Moss et al., 2010).The fact that the responses of most variables does not scale linearly with the forcing in transient climate changes has important implications for the scaling of climate change patterns based on simple energy balance models.

Temperature response and forcing-feedback decomposition
For the rest of the study, the responses are scaled by their respective adjusted forcings only in figures concerning the linear additivity assumption, i.e.Figs.1b-d and 5.In the other cases, the raw responses will be considered, for several reasons.First, scaling the responses would make the assumption that each variable at each grid point scales linearly with the adjusted forcing.While scaling the responses of diagnostic variables might be justified, other quantities such as the zonal mean profile of specific humidity or residence time of water vapor in the atmosphere cannot necessarily be scaled with the adjusted forcing.Further, the differences in the raw responses between CO 2 and solar scenarios that will be discussed below are large and scaling the responses would not change the conclusions and/or make those differences even larger.
The larger temperature changes at the surface in the CO 2 compared to solar scenarios seen in Fig. 1a are caused by the larger adjusted forcing values mentioned above.The scaled global mean surface temperature anomalies are almost the same for the CO 2 and solar cases, although slightly larger in the solar cases due to the fact that solar forcing is more efficient at warming the surface.As seen in Fig. 1a for the 4x scenario, the temperature increase appears to be strong during the first 70 yr but after CO 2 402 Introduction

Conclusions References
Tables Figures

Back Close
Full concentrations remain constant, a decrease occurs, followed by a slow increase.This "overshoot" comes from the northern high latitudes (see region definitions below), the temperature increase is monotonic in the other sub-regions (not shown).An explanation for this behavior might be that when the Meridional Overturning Circulation (MOC) is reduced as it occurs in warming scenarios, the ocean transiently takes up a lot of energy which is compensated by a temporary decrease in atmospheric temperature (Knutti and Stocker, 2000).It is however difficult to judge whether this mechanism is a robust response of the climate system or just a feature of the CCSM3.5 model.Motivated by the fact that precipitation in CCSM3.5 is expected to increase in high latitudes and near the equator but decrease in mid-latitudes (see Schaller et al., 2011, Fig. 3), anomalies are also computed for three sub-regions: the tropics (15 • S-15 • N), the mid-latitudes ( 50• S-15 • S and 15 • N-50 • N) and the high latitudes ( 90• S-50 • S and 50 • N-90 • N).The surface temperature anomalies of the last 30 yr of each simulation scaled by their respective adjusted forcings are shown in Fig. 1b-d for the three subregions.According to the Kolmogorov-Smirnov test, all but one pair of scenarios are significantly different from each other, indicating that the non-validity of the linear additivity assumption does not depend on the region considered in solar scenarios but that it arises from higher latitudes in CO 2 scenarios.Figure 1b-d further shows that the scaled annual mean temperature responses to CO 2 and solar forcing in different regions are significantly different, although the idea underlying the radiative forcing concept would suggest that they are similar.The scaled surface temperature response to solar forcing is larger at low latitudes while it is larger at high latitudes in the CO 2 scenarios.This results from the fact that solar forcing acts primarily at low latitudes and, although also present in solar scenarios, polar amplification is stronger in CO 2 scenarios.
As a further step to investigate the non-linearity of the responses, the method described by Andrews (2009) is used to separate the total response into a part caused by the forcing directly and a part resulting from the change in surface temperature, i.e. the feedback-dependent part.We repeat the same analysis as in Andrews (2009) for the Introduction

Conclusions References
Tables Figures

Back Close
Full five ensemble members of all CO 2 and solar scenarios and for the energy budget fluxes as well as precipitation.The anomaly N i in a given flux i is given by: where F i is the part of this anomaly caused by the forcing directly, α i is the feedback parameter for the given energy flux i , and ∆T is the change in surface temperature.
By regressing the annual mean values from year 71 to 100 of N i against the annual mean values of ∆T , one can determine the feedback parameters α i for each of the five ensemble members separately in the five different scenarios.Once α i is determined, the forcing F i and feedback parts α i ∆T can be calculated for each of the five ensemble members, again separately, in the five scenarios.Andrews (2009) assumes that the feedback parameters are forcing-independent.We find this assumption to be valid in our CO 2 and solar forcing scenarios (not shown).For most variables however, the feedback parameters appear not to be well constrained in CCSM3.5.Reasons might be that the inter-annual variability is large or that the climate change signal is weak for those variables.The feedback parameter for precipitation however is well constrained and repeating the linear additivity test shows that the non-linearity arises from the feedback part in both CO 2 and solar forcing scenarios.The forcing part adds linearly but the feedback part is larger by 36 % in 4x compared to 2x and by 20 % in 74 compared to 37.This is expected since temperature also behaves non-linearly as shown previously.
Further, the regression method applied to CCSM3.5 confirms results from previous studies showing that the forcing contribution is negative in CO 2 scenarios but close to zero in the solar case (Bala et al., 2010;Andrews and Forster, 2010;Cao et al., 2011).The feedback part is positive for CO 2 and solar forcings but is larger in the CO 2 case since it follows the larger temperature increase in these unscaled scenarios.For the rest of the study the focus will be on differences in the processes, and therefore the comparison between CO 2 and solar scenarios of the same intensity.

Changes in the energy budget
The net shortwave (SW) flux change at TOA in the solar scenarios roughly corresponds to the imposed solar forcings as shown in Fig. 2b.This is not necessarily the case in all models, but in CCSM3.5 the SW clear-sky flux change is very large (around 11 W m −2 in 74) and is balanced by a relatively large increase in SW cloud radiative forcing (SWCRF) of more than 3 W m −2 (see Fig. 3h).The net SW flux change in the CO 2 scenarios is much smaller than the SW clear-sky flux change and cloud radiative forcing.In the solar scenarios, the longwave (LW) fluxes become more negative, indicating that as a consequence of the heating of the surface, the Earth emits more LW radiation to space (see Fig. 2a).In contrast, in the CO 2 scenarios, the surface also becomes warmer during the simulation but the increasing CO 2 concentrations trap the upward LW radiation, hence a positive anomaly at TOA.For TOA LW, the differences between high and low forcing scenarios are much larger for solar than for CO 2 scenarios, where the values are almost the same.This is due to the strongly negative LW cloud radiative forcing in 4x (see Fig. 3g).
The response of the system at TOA to solar forcing is large in both the SW and LW compared to the relatively small responses at TOA for CO 2 forcing, although one might expect that solar forcing would act mostly in the SW range and CO 2 in the LW range (Meehl et al., 2003;Bala et al., 2010).Still, the net TOA flux anomalies are in the same range for both forcings and would be almost equal if they were scaled by their respective adjusted forcing values (see Fig. 2c).Looking at the net energy flux changes at the surface (turbulent and radiative) in Fig. 2i, the picture is very similar indicating that the atmosphere only takes up a small part of the increased energy in the system.This is consistent with observations indicating that most of the excess energy is taken up by the oceans (Levitus et al., 2001).
Considering the SW and LW fluxes at the surface again highlights the different processes occuring in the CO 2 and solar scenarios.by the surface in the CO 2 scenarios compared to the control simulation, the opposite is true in the solar scenarios (see Fig. 2e).Besides the obvious fact that there is less net SW radiation at TOA in CO 2 compared to solar scenarios, this outcome can be further explained by an increase in low clouds in CO 2 scenarios (see Fig. 3f) rather than by changes in surface albedo.The time series of the surface net LW flux anomalies reflect surface warming and show the large increase in back radiation in the CO 2 scenarios (see Fig. 2d).In those scenarios, changes in LW back radiation are caused by increases in CO 2 , water vapor and low clouds, while in the solar scenarios, only the increasing greenhouse effect of water vapor and low clouds is seen.In addition, changes in surface temperature are larger in CO 2 scenarios, which in itself causes a larger LW back radiation, and consequently larger increases in water vapor.
Changes in global annual mean precipitation can be understood either from an atmospheric (Mitchell et al., 1987;Allen and Ingram, 2002) or a surface energy budget perspective (Boer, 1993;Wild et al., 2008;Andrews et al., 2009).Both interpretations have been shown to be valid and from a surface energy budget point of view, changes in the global hydrological cycle are primarily driven by changes in the energy available at the surface for the turbulent energy fluxes.The increase in net surface radiative flux, i.e. the sum of the changes in net surface LW ∆LW surf and in net surface SW ∆SW surf (left hand side of Eq. ( 2) and Fig. 2f), indicates that there is potentially more energy available at the surface for evaporation in the solar than in the corresponding CO 2 scenarios, mostly due to the negative surface SW anomaly in the CO 2 scenarios.In this case, a scaling of the variables would make the differences between CO 2 and solar scenarios even larger.sensible heat is strongly dependent on the forcing agent, and has a lower magnitude in solar compared to CO 2 scenarios (see Fig. 2h).SH decreases due to a decrease in the air-sea temperature difference at the surface caused by an increased opacity in the LW of the atmosphere associated with the increase in specific humidity (see Fig. 4b and d) (Stephens and Ellis, 2008;O'Gorman et al., 2012).The SH decrease is less pronounced in solar scenarios due to the fact that solar forcing primarily acts at warming the surface while CO 2 forcing induces a warming of the whole troposphere.
The SH decrease is however dominated by the ocean response: over land, SH slightly increases (around 1 W m −2 ) due to surface warming, and again, increases more in solar scenarios due to the strong SW forcing.

Changes in the hydrological cycle
Understanding how clouds will change in the future remains a major challenge since climate models have difficulties in simulating low clouds in particular (Stephens, 2005;Zhang et al., 2005).Total cloud cover increases in all five scenarios with no strong dependency on the forcing agent (not shown) but most CMIP3 and CMIP5 models actually show a decrease in total cloud cover (Bender, 2011;Andrews et al., 2012).Gregory and Webb (2008) as well as Andrews and Forster (2008) showed that fast cloud adjustments caused by an increase in CO 2 forcing are responsible for most of the spread among models.After removing these fast adjustments, cloud feedbacks appear more consistent in the different models.Low cloud cover increases more in solar scenarios as shown in Fig. 3f.These clouds likely have a cooling effect because they reflect SW radiation, which would act as a negative feedback to the climate system.The changes in high-level cloud cover shown in Fig. 3d are very small and they should be interpreted carefully.Still, increasing high cloud cover could lead to more warming because high clouds may trap upward LW radiation, which would imply a positive feedback in CO 2 scenarios and a negative feedback in solar scenarios due to the slight increase and decrease in high cloud cover, respectively, in these scenarios.Introduction

Conclusions References
Tables Figures

Back Close
Full Vertical profiles of temperature anomalies in 4x and 74 are shown in Fig. 4a and c respectively.The warming of the troposphere and cooling of the stratosphere is characteristic for CO 2 -induced climate change while the stratosphere is warming overall in solar scenarios.Given the Clausius-Clapeyron relationship, atmospheric water vapor change will closely follow temperature change, at a rate of around 7 % K −1 for temperatures typically found at the surface, but by a rate of around 15 % K −1 at temperatures encountered in the upper troposphere.In the lower troposphere, the temperature increase is stronger in the unscaled CO 2 scenarios and consequently, specific humidity increases more relative to solar scenarios, as shown in Fig. 4b and d.On the other hand, latent heat flux increases more in solar scenarios, hence a larger precipitation increase compared to CO 2 scenarios (see Fig. 3a).Mitchell et al. (1987) showed that changes in precipitation due to warming occur at a smaller rate than changes in water vapor because the former are limited by the heat balance of the atmosphere.Even though both specific humidity and precipitation increase with warming, there is a difference between their respective increase rates, which must be compensated by a weakening of the atmospheric circulation (Held and Soden, 2006).In global climate model projections, Held and Soden (2006) found that the change in convective mass flux decreases as a consequence of the fact that precipitation increases only by about 2 % K −1 while water vapor increases by 7 % K −1 .Here, the precipitation response in CO 2 scenarios is more muted compared to the response in water vapor than in solar scenarios.A weakening of the atmospheric circulation is equivalent to an increase in residence time of water vapor in the atmosphere, i.e. the ratio between total precipitable water and global precipitation rate, as shown by Douville et al. (2002) and Bosilovich et al. (2005).The residence time in the control simulation is around nine days and it increases in all scenarios by up to almost two days (see Fig. 3i).However, the residence time of water in the atmosphere is larger in CO 2 scenarios compared to solar scenarios of the same intensity, implying a weaker atmospheric circulation in the CO 2 scenarios.Since residence time is defined as the ratio of two variables, scaling the responses would lead to the same result.Introduction

Conclusions References
Tables Figures

Back Close
Full To obtain additional information on the processes responsible for changes in precipitation, changes in large-scale and convective precipitation are considered separately.The unscaled time series of global mean anomalies in large-scale and convective precipitation are shown in Fig. 3b-c and the scaled regional anomalies for large-scale and convective precipitation are shown in Fig. 5.The linear additivity assumption is also tested in the three sub-regions in a similar way as in Fig. 1b-d.For convective precipitation, the linear additivity assumption is almost never valid (see Fig. 5a-c) and in the high latitudes, the direction of change is even inverted between 37 and 74.For large-scale precipitation, the assumption of linear additivity is valid in mid-latitudes for all scenarios but also in high latitudes for the CO 2 cases (see Fig. 5d-f).The invalidity of the linear-additivity assumption therefore cannot easily be attributed to specific regions.However, large-scale precipitation processes seem to be better represented by the assumption than convective precipitation.
Focusing again on the differences between CO 2 and solar scenarios, the changes in convective precipitation seem to follow the changes in surface temperature shown in Fig. 1b-d.Convective precipitation increases more in the tropics and mid-latitudes in solar scenarios due to the stronger warming while in the high latitudes, the convective precipitation response is larger in CO 2 scenarios due to the stronger polar amplification in these scenarios (see Figs. 5c and 1d).According to Held and Soden (2006), the decrease in convective mass flux should be smaller in solar simulations due to a smaller difference between the increase in moisture and precipitation.This would be in line with the stronger increase in convective precipitation in solar compared to CO 2 scenarios in the tropics, while the changes in the other sub-regions are more comparable (see Fig. 5a-c).However, in the CCSM3.5 model, the stronger hydrological sensitivity to solar forcing compared to CO 2 forcing per unit of surface warming is primarily caused by a larger increase in large-scale precipitation (see Fig. 3b). Figure 5d-f further shows that the main difference between CO 2 and solar scenarios arises in the mid-latitudes, where large-scale precipitation decreases in CO 2 scenarios but increases slightly in solar scenarios.Introduction

Conclusions References
Tables Figures

Back Close
Full The responses in large-scale precipitation are complex to interpret as they not only depend on energetic and thermodynamic constraints but also on changes in atmospheric circulation in the system.Therefore, changes in the pole-to-equator surface air temperature gradient are calculated to investigate changes in atmospheric circulation.We use an adapted version of the meridional temperature gradient (MTG) index defined by Gitelman et al. (1997) to quantify differences between the scenarios.Using the sub-regions defined above, we calculate a MTG index in each hemisphere as the difference between the respective high latitudes and the tropics.With increasing global mean temperature and as already observed, Gitelman et al. (1997) showed that the MTG decreases, leading to weaker mid-latitudes eddies.In our simulations, the pole-to-equator gradient in the Southern Hemisphere decreases for all scenarios, in particular in the CO 2 cases (see Fig. 6b).In the Northern Hemisphere, the MTG index also decreases except for 37 where it slightly increases, which, according to Gitelman et al. (1997), would imply stronger mid-latitudes eddies (see Fig. 6a).
Finally, energy budgets for the 4x and 74 simulations are calculated in both high latitude boxes (90 • S to 50 • S and 50 • N to 90 • N) to quantify the changes in poleward energy transport through the atmosphere.The poleward atmospheric heat transport (AT) is calculated as the difference between the net energy flux at top of atmosphere and the net energy flux at surface in each box.In the last 30 yr of constant forcing the atmosphere is in quasi-equilibrium and thus we assume zero heat storage in the high latitude boxes to calculate the meridional atmospheric energy transport (AT) as in Rugenstein et al. ( 2013). Figure 7 shows that the net TOA energy flux anomaly is much larger in 4x for both boxes and interestingly, it is even negative in the northern high latitudes in 74.While the net surface energy flux anomalies are more similar between CO 2 and solar simulation, the separation in LW, LH and SH appears more forcing-dependent.
Note that the anomalies shown are not scaled but scaling them with the adjusted forcing would make the anomalies in 74 increase, or decrease if negative, relative to the anomalies in 4x.As a consequence of the difference between both scenarios in the net TOA energy flux anomalies, the increase in poleward energy transport through the Figures

Back Close
Full atmosphere is significantly larger in solar scenarios (AT in Fig. 7).Since a solar forcing increase acts primarily in the tropics, the energy received has to be redistributed.In contrast, a CO 2 forcing is homogeneously distributed over the globe and produces a stronger polar amplification, causing the temperature difference between poles and equator to decrease.Therefore, the climate system maintains a stronger atmospheric circulation in solar compared to CO 2 scenarios, but in both scenarios the atmospheric circulation is weaker compared to the control simulation.Given that moisture increases in both scenarios compared to the control simulation, the same or a larger transport of energy, as Fig. 7 shows, can still be achieved by a weaker atmospheric circulation (Trenberth, 2011).Although large-scale precipitation depends not only on the poleward atmospheric energy transport, the presented results are well in line with the fact that large-scale precipitation increases more in solar scenarios compared to CO 2 scenarios at higher latitudes (see Fig. 5e-f).

Conclusions
In order to better understand changes in the energy budget and hydrological cycle temperature or total radiative forcing.The linear additivity assumption is also tested for surface temperature, large-scale and convective precipitation in the tropics, midlatitudes and high latitudes and appears to be not valid in general, regardless of the sub-region considered.
The best estimate of the radiative forcing of a doubling of CO 2 provided by Myhre et al. ( 1998) is often used in the literature to construct a solar simulation of the same intensity.This is what is done here but the temperature response to solar forcing is found to be significantly smaller than CO 2 forcing, as shown in previous studies (e.g.Hansen et al., 2005;Schmidt et al., 2012).On one hand, this is due to the fact that the efficacy of solar forcing is smaller than the efficacy of CO 2 forcing, but on the other hand, the adjusted forcing of a quadrupling of CO 2 appeared to be more than twice the best estimate for a CO 2 doubling provided by Myhre et al. (1998).This outcome indicates limitations in the traditional definition of radiative forcing.It is therefore central to consider the efficacy of a forcing agent (Hansen et al., 2005), but also alternative techniques for calculating forcings (e.g.Hansen et al., 2002;Shine et al., 2003).
Following Andrews (2009), the forcing and feedback-dependent terms are calculated and the non-linear response in precipitation is shown to be caused by the feedbackdependent term.To provide new insights on reasons for the weaker increase in precipitation in CO 2 compared to solar forcing scenarios, the physical processes altering the components of the energy budget and hydrological cycle in the CCSM3.5 model are analyzed.When the raw modelled responses are considered and no scaling performed, surface warming appears weaker in solar compared to CO 2 scenarios.Given the Clausius-Clapeyron relationship, the water vapor increase in the lower troposphere follows the temperature increase and is therefore also weaker in solar scenarios.At the same time, the energy available at the surface for evaporation, and consequently precipitation as both are equal in the global annual mean, is larger in the unscaled solar scenarios.For the unscaled CO 2 scenarios, the opposite is true: the water vapor increase is relatively strong but the precipitation increase is relatively weak compared to solar scenarios.The more muted the water vapor and precipitation increases are, the Figures weaker the atmospheric circulation will become.This is shown by the longer residence times in CO 2 scenarios compared to solar scenarios.Finally, these conclusions are supported by the fact that the decrease in meridional temperature gradient expected with global warming is larger in the CO 2 case, indicating that less energy is available for mid-latitude eddies.Energy budgets over high latitude boxes in each hemisphere confirm that the poleward energy transport in the atmosphere increases more in solar scenarios.These results are in line with the stronger increase in large-scale precipitation in solar scenarios, while the increase in convective precipitation is more similar between CO 2 and solar forcing scenarios.Scaling the responses in each scenario by their corresponding adjusted forcing as is done to test the linear additivity assumption leads to the same conclusions.
It is important to stress that the presented results are based on one global climate model and cannot claim universal validity.Still, it would be useful to compare them with the same scenarios performed with different GCMs to assess whether the described processes are robust.For the IPCC reports, the focus is on comparing the models running the SRES and RCP scenarios where the composition of the forcings is complex.Pendergrass and Hartmann (2012) showed that in simulations used for the IPCC Fourth Assessment Report (IPCC, 2007) the black carbon forcing was not specified and was therefore different in each model, making model inter-comparison with such simulations difficult to interpret.Understanding how individual models respond to a given forcing is only possible with idealized simulations as presented in this study.We suggest that this could represent a systematic way to inter-compare and to some extent evaluate the models, and might provide additional information to classical model evaluations against observed climatological mean fields.The ability to reproduce the observed climate is a priority in model inter-comparisons, but since observational errors remain substantial, in particular over oceans for precipitation (Liu et al., 2012), inter-comparing models in simple and idealized simulations might represent a complementary tool.Simulations with a one percent per year increase in CO 2 are available in the CMIP archive and one could imagine that the modeling groups could perform Introduction

Conclusions References
Tables Figures

First, a 600
-year present-day control simulation is run with the CCSM3.5 model and equilibrium is reached after around 150 yr.The simulations with CO 2 and solar forcings are ensemble members that branch out at years 400, 430, 450, 470 and 500 of the control simulation to sample different initial conditions.Each case (scenario hereafter) Discussion Paper | Discussion Paper | Discussion Paper | therefore consists of five initial condition ensemble members to robustly quantify the model internal variability.The first scenario is a 1 % yr −1 transient increase of CO 2 up is used to estimate the adjusted forcing as the intercept in a linear regression of changes in net radiation at Top Of Atmosphere (TOA) against the changes in surface temperature.The end of the transient phase of all five scenarios is reached during year 69 relative to the branching point of the scenarios, and for the rest of the study, we consider annual mean anomalies of the years 71 to 100 compared to the average over 100 yr in the control simulation.All fluxes are defined as positive downward and consequently, positive anomalies represent a gain of energy for the surface, and negative ones, a loss.Discussion Paper | Discussion Paper | Discussion Paper | Shiogama et al. (2010b) identified limitations of this assumption in the form of potential Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | While less SW radiation is absorbed 405 Discussion Paper | Discussion Paper | Discussion Paper | 2) The change in net radiative flux at the surface is partitioned into changes in latent and sensible heat flux (∆LH and ∆SH respectively) and some of this excess energy will be taken up by the ocean ∆NET surf .The major part of the net radiative flux at the surface is used to increase the latent heat flux as shown in Fig. 2g.The decrease in 406 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | in a warmer climate, idealized simulations are performed with the NCAR CCSM3.5 model.The scenarios are also designed to test the assumption of linear additivity of the response to different forcing agents and magnitudes in a fully coupled climate model for transient climate change.The responses of most variables of the energy budget and hydrological cycle do not add linearly in the 30 yr after stabilization in the global mean.The fact that most climate variables as simulated by CCSM3.5 do not respond linearly to forcings can be relevant for detection and attribution and pattern scaling techniques.Depending on the application, the errors introduced by assuming linear additivity when it does not apply might be considered negligible or not.In any case, these results cannot be captured properly by models of lower complexity, which are often used to inform policy makers or for impact studies, and are implicit when characterizing the overall magnitude of climate change or a target for stabilization in terms of global mean Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | similar experiments with solar forcing only, and even with aerosol forcing.Then robust physical responses across the models could be identified and used as a complementary piece of information.Discussion Paper | Discussion Paper | Discussion Paper | Held, I. M. and Soden, B. J.: Robust responses of the hydrological cycle to global warming, J. Climate, 19, 5686-5699, 2006.408, 409 IPCC: Climate change 2007: The physical science basis, in: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M.Discussion Paper | Discussion Paper | Discussion Paper | Stouffer, R. J.: Time scales of climate response, J. Climate, 17, 209-217, doi:10.1175/1520-0442(2004)017<0209:TSOCR>2.0.CO;2, 2004.402 Trenberth, K. E.: Changes in precipitation with climate change, Clim.Res.
Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 1 .Fig. 2 .Fig. 3 .Fig. 4 .Fig. 5 .
Fig. 1.(a) Time series of the global annual mean anomaly of surface air temperature (TAS) [K] for all five scenarios.Averages over five initial condition ensemble members for the 100 yr simulations and over 3 members for the 300 yr simulations are shown.Shading extends from the lowest to highest value obtained by any of the five ensemble members.Regional surface air temperature anomalies scaled by the adjusted forcing [K Wm −2 ] in the scenarios are shown for (b) the tropics, (c) the mid-latitudes and (d) the high latitudes (see text for definition of the sub-regions).Asterisks are displayed whenever the pairs of scenarios (2x and 4x, 37 and 74, 372x and 2x+37) are significantly different from each other (two-sample Kolmogorov-Smirnov test, significance level 0.05) and therefore indicate when the linear additivity assumption is not valid.Red crosses are outliers.

Fig. 6 .
Fig. 6.Annual mean meridional surface air temperature gradient calculated as the difference between 15 • S-15 • N and 90 • N-50 • N or 50 • S-90 • S for (a) the Northern (NH) or (b) the Southern Hemisphere (SH) for all five scenarios.Averages over the five ensemble members are shown.

Fig. 7 .
Fig. 7. Annual mean anomalies of the Earth's energy budget components for (a) the 4x and (b) the 74 simulation.The energy budget anomalies are calculated for two boxes (90 • S to 50 • S and 50 • N to 90 • N) and the change in poleward atmospheric energy transport (AT) is also indicated.Values are in W m −2 and energy fluxes are defined as positive downward and northward.

Table 1 .
The response in a given variable Figures in a low scenario scaled by the adjusted forcing in this scenario is compared with the response in the scaled high forcing scenario.A two-sample Kolmogorov-Smirnov test is used to assess if the five ensemble members of the compared scenarios are significantly different.In Table1, the differences (response of the scaled high forcing scenarios minus scaled response of the low forcing scenarios) are shown as percentages of the reference value, i.e. the scaled responses to the low forcing scenarios.Positive and negative values therefore mean that for a given variable, the high forcing scenario responses are significantly larger and smaller, respectively, than expected from the linear additivity assumption.No values are shown for variables where the linear additivity assumption is valid, i.e. differences are not significantly different.The values shown represent non-linearities arising from long-term feedbacks.Overall, the results shown in Table . 395 Watterson, I. G.: Calculation of probability density functions for temperature and precipitation change under global warming, J. Geophys.Res.-Atmos., 113, D12106, doi:10.1029/2007JD009254,2008.396 Wild, M., Grieser, J., and Schaer, C.: Combined surface solar brightening and increasing greenhouse effect support recent intensification of the global land-based hydrological cycle, Geophys.Res.Lett., 35, L17706, doi:10.1029/2008GL034842,2008.395, 406 Zhang, M. H., Lin, W. Y., Klein, S. A., Bacmeister, J. T., Bony, S., Cederwall, R. T., Del Genio, A. D., Hack, J. J., Loeb, N. G., Lohmann, U., Minnis, P., Musat, I., Pincus, R., Stier, P., Suarez, M. J., Webb, M. J., Wu, J. B., Xie, S. C., Yao, M. S., and Zhang, J. H.: Comparing

Table 1 .
Results from the linear additivity test on chosen variables of the energy budget and hydrological cycle scaled by the adjusted forcing in the scenario.Positive values indicate that the high forcing scenario has a larger response than expected doubling/addition of the low forcing scenarios.Only statistically significant differences are shown, the linear additivity assumption is valid where no value is given.