Exploring the sensitivity of Northern Hemisphere atmospheric circulation to different surface temperature forcing using a Statistical Dynamical Atmospheric Model

Climate and weather conditions in the mid-latitudes are strongly driven by the large-scale atmosphere circulation. Observational data indicates that important components of the large-scale circulation have changed in recent decades including the strength and the width of the Hadley cell, jets, storm tracks and planetary waves. 15 Here, we use a new statistical-dynamical atmosphere model (SDAM) to test the individual sensitivities of the large-scale atmospheric circulation to changes in the zonal temperature gradient, meridional temperature gradient and global-mean temperature. We analyse the Northern Hemisphere Hadley circulation, jet streams, storm tracks and planetary waves by systematically altering the zonal temperature asymmetry, the meridional temperature gradient, and the global mean temperature. Our results show that the strength of the Hadley cell, storm tracks and jet streams depend in terms of relative 20 changes almost linearly on both the global mean temperature and the meridional temperature gradient whereas the zonal temperature asymmetry has little or no influence. The magnitude of planetary waves is affected by all three temperature components, as expected from theoretical dynamical considerations. The width of the Hadley cell behaves nonlinearly with respect to all three temperature components in the SDAM. Moreover, some of these observed large-scale atmospheric changes are expected from dynamical equations and therefore an important part of model validation. 25


Introduction
Large-scale atmosphere dynamics including Hadley cells, jet streams, storm tracks and planetary waves play a key role in the general circulation of the atmosphere, determining climatic conditions worldwide. There is emerging evidence that these largescale dynamics have significantly changed in recent observations and/or might do so in climate change projections of the 21st 5 century.
Hadley cells are large-scale atmospheric circulations in the tropics and responsible for the transport of heat and moisture from the equator to the mid-latitudes (D'Agostino and Lionello, 2016). Recent literature suggests that the Hadley cells have widened in the past few decades due to ozone depletion in the Southern Hemisphere (Kang et al., 2011;Son et al., 2009) and due to increases in different warming agents as carbon dioxide and tropospheric ozone in the Northern Hemisphere (Allen et al., 10 2012). Further possible drivers of the Hadley cell dynamics are sea surface temperature variations, which may lead to tropical contraction (Allen et al., 2014), stratospheric cooling, global warming or changes in baroclinic eddy phase speeds (Chen et al., 2014). While most studies based on reanalysis products observe that together with this widening there is also a strengthening of the Hadley cells (Nguyen et al., 2013), modelling studies project a weakening over the 21 st century under a high emission scenario . This discrepancy may be explained by relatively short observational records, large natural variability 15 or model deficiencies (Allen et al., 2014).
Changes in the strength and the width of the Hadley cell circulation have strong implications for a variety of atmospheric phenomena such as jet streams, extratropical storms and planetary waves. Jet streams are upper-level fast currents of westerly winds that circulate around the hemisphere. The jet stream is caused by momentum conservation, when air of the Hadley cell moves poleward and by meridional contrasts in solar heating (Woollings, 20 2010). Therefore, changes of the Hadley cell directly influence the jet. Storm tracks are associated with fast-traveling Rossby waves that move from west to east, which drive the jet stream. (Totz et al., 2018). Under climate change, a weakening of the northern subtropical jet stream in winter is observed by most studies (Molnos et al., 2017;Rikus, 2015) though a few studies observe also a strengthening trend (Pena-Ortiz et al., 2013). The underlying drivers for those changes are still debated. One possible driver is the Arctic amplification leading to a decreasing meridional temperature gradient and a weakening of the jet 25 stream and storm tracks.
Storm tracks play a crucial role in modulating precipitation in the Earth system (Hawcroft et al., 2012;Lehmann and Coumou, 2015;Raible et al., 2007). Yin (2005) studied storm tracks under climate change using 15 coupled climate models and found that storm tracks shift poleward and intensify under climate change in winter In addition, O'Gorman, (2010) discovered that the storm track intensity is nonlinearly related to global warming. Lehmann and Coumou (2015) show that during the satellite 30 period (1979-2013) storm tracks have significantly weakened in boreal summer. They argue that the weakening might be related to a decreasing meridional temperature gradient. Consistently, CMIP5 models project a further decrease in summer storm track activity under a high-emission scenario until the end of the 21st century (Lehmann et al., 2014).
In contrast to the changes in these large-scale circulation components, a recent study by Barnes and Polvani (2013) indicates that there are no significant trends in the strength of long-term quasi-stationary planetary waves under climate change.
Planetary waves are disturbances of longer time-scales (beyond 15 days) and are manifested by a displacement of the circumpolar flow out of zonal symmetry, as is apparent at high latitude. They lead to a gradual undulation of the jet stream about latitude circles (Salby, 2012). These planetary waves strongly interact with storm tracks in the mid-latitudes and impact 5 their location and activity. (Branstator, 1994;Hoskins and Valdes, 1990) . Especially the barotropic component of the planetary waves has a strong influence on the storm tracks and can change the distribution of storms by steering the disturbances (Branstator, 1994).
The large-scale atmospheric circulation is mainly a result of large-scale temperature differences between the equator and the poles, and between the ocean and land. A further key driver is the angular momentum as shown by the study of Egger et al. 10 (2007) Different modelling studies have thus tried to examine the influences of different temperature sources on changes in mid-latitude circulation. Using an idealized general circulation model, Butler et al. (2010) separate the temperature effects by using different heating sources. Their main findings are that a warming tropical troposphere causes poleward shifts in the extratropical tropospheric storm tracks and a weakened stratospheric Brewer-Dobson Circulation, a polar stratosphere cooling leads also to a shift of the extratropical storm tracks and a warming at the polar surface results in a equatorward shift of the 15 storm tracks. With their approach Butler et al. can attribute which forcing has the most important influence on the shift of jet streams, storm tracks etc. In a study from 2011 Butler et al. presented an alternative perspective on the response of the midlatitude tropospheric circulation to zonal-mean tropical heating. The projection of the heating onto the isentropic surfaces at extratropical latitudes drive the poleward shift in wave generation at lower levels. In addition, the poleward shift in the heat fluxes within the troposphere and the diffusive nature of eddy fluxes of the polar vortex lead to a poleward shift in wave 20 breaking near the tropospause.
However, none of the above-mentioned studies analyzed the impact of zonal temperature asymmetry and meridional temperature gradient as well as global mean temperature by changing the temperature components directly without adding a heat source or increasing CO2. For that reason, changes in the Hadley cell, jet streams, planetary waves and storm tracks might not be directly assigned to individual temperature components. 25 In this study, we analyse the full temperature space with a novel systematic approach using a statistical-dynamical atmospheric model (Totz et al., 2018). This way we can scan for 'non-linearities' in the system (i.e. the HCE might be very sensitive to meridional temperature gradients only for a narrow range of meridional temperature gradient values, and outside of that range it is not sensitive). It is important to know such non-linearities as it could imply more abrupt changes under global warming.
The statistical-dynamical atmosphere model Aeolus 1.0 (Totz et al., 2018) is explained in section 2 in more detail. In section 30 3, we describe the data used for the experiments, the separation of the temperature components and the analysed dynamical variables. In section 4 we present the results and in section 5 compare them with the literature. In section 6 we provide a short discussion on the robustness and interpretation of our results. We conclude with a summary in section 6.
It is a 2.5dimensional model with the vertical dimension coarsely resolved and therefore belongs to the model class of 5 intermediate complexity atmosphere models (Claussen et al., 2002;Petoukhov et al., 2000). Aeolus is based on time-averaged (over short time scales) equations in which transient eddies are parameterized in terms of the large scale field (Peixoto and Oort, 1992;Saltzman, 1978). This means that instead of solving every eddy directly, only the ensemble mean eddy characteristics (in terms of heat, water vapour and momentum transport) are solved. The essential difference compared to more widely used general circulation models (GCMs) is thus the point of truncation in the frequency spectrum of atmospheric motion 10 (Saltzman, 1978). This different approach allows much coarser spatial and temporal discretization, making SDAMs computationally efficient, because the synoptic waves are parameterized in terms of the large scale wind field which is the basic idea of a statistical-dynamical method. The model includes the two-way dynamical interactions between parameterized storms and resolved westerlies. A full description can be found in Coumou et al. (2011).
The model has 5 vertical levels in the troposphere with the model top at 10000m altitude. Aeolus 1.0 has a "dummy" 15 stratosphere (i.e. its physics and dynamics are not resolved) to have a boundary condition at the top of the troposphere. In this experiment we excluded topographic influences and it is an atmosphere-only setup using prescribed sea level temperatures. Aeolus 1.0 contains a full hydrological cycle consisting of three-layer stratiform cloud plus convective cloud scheme as presented and validated in Eliseev et al. (2013). The convective plus 3-layer stratiform cloud scheme includes low-level, midlevel and upper-level stratiform clouds. The equation for humidity is a prognostic equation and described in Petoukhov et al. 20 (2000).
Parameter values for the dynamical core were taken from the calibration process, as described in Totz et al. (2018), which optimizes the model's representation of the tropospheric large-scale circulation.
For more information such as comparisons of Aeolus 1.0 with GCMs or the equations for planetary waves, zonal mean meridional and zonal mean velocities as well as the azonal wind velocities, we refer the reader to Totz et al. (2018). 25

Aeolus forcing parameters
The simulations are forced by climatological  winter mean (December-January-February (DJF)) data of sea surface temperature and specific humidity at the surface, using ERA-Interim data from the European Centre for Medium-Range Weather Forecasts (ECMWF) (Dee et al., 2011). First, the data is regridded to 3.75°× 3.75° (longitude × latitude). In this experiment, atmospheric temperatures extrapolated to sea level and humidity are prescribed to decouple the dynamics from diabatic heating and associated temperature changes. This way, the dynamical core equilibrates to the prescribed temperature and humidity patterns without any additional complicating factors. 5 In Aeolus stand-alone mode, the output is on the same time step as the input, i.e. seasonal means. However, internally, the time-stepping of the solver is one hour resolution. In Aeolus, the synoptic wind fields u' and v' are parameterized in terms of the large-scale long-term fields (i.e. the seasonal mean fields). Detailed information is provided in Coumou et al. (2011).
We then test how different changes in the temperature profiles affect different aspects of the circulation. We do thousands of individual simulations to disentangle and separately analyse the effect of global mean temperature, equator-to-pole temperature 10 gradient and east-west temperature differences.

Specifications of the surface temperature
In this and the following sections, the angle brackets denote time-averaged quantities, the overbar denotes zonal mean quantities, the prime indicates synoptic scale components (2 -6 days period), and the star indicates azonal components, i.e. 15 deviations from the zonal mean.
For the sensitivity analysis we vary three different temperature components: (1) the meridional temperature gradient , the zonal temperature asymmetry * (i.e. deviations from the zonal mean) and (3) the global mean temperature .
We change the temperature for each grid cell with respect to parameters for the three components in three steps. First, the parameter is used to vary the meridional temperature gradient by cooling/warming the poles 20 whereby and are respectively latitude and longitude, ( , ) is the original temperature, ( ) the temperature at the equator and 1 ( , ) is the altered temperature. In temperature 1 only the meridional temperature gradient is altered/ updated and the temperature change of the azonal component as well as the global mean temperature are done in the next steps. The parameter is varied between 0.75 and 1.1 (with steps of 0.025) to examine the behavior of the dynamical core under conditions between −25% to +10% of their present-day wintertime climatological values. 25 In the second step the global mean temperature is adapted to the originally global mean temperature and global mean temperature is varied by the parameter Whereby Mean( 1 ) is the global mean temperature of 1 and Mean( ) is the global mean temperature of .
In the third step, the parameter is used to alter the azonal temperature, which is added to the zonal mean temperature 1 .
This way = 1, = 1 and = 0 indicate present day conditions. The temperature perturbation Final is the final temperature based on all three temperature components and is used as model input to which the dynamical core 5 equilibrates. The temperature perturbations are applied at sea level and propagate to the upper levels based on the lapse rate equation. A schematic plot of the different temperature perturbations is shown in Fig. 1.
Also the parameter is varied between 0.75 and 1.1 (with steps of 0.025). These limits roughly correspond to expected temperature gradients during the last glacial maximum scenario and in a 2 × 2 scenario (Coumou et al., 2011). Large zonal temperature differences, i.e. large values of , imply strong temperature deviations between land masses and oceans. 10 Small values represent amplified warming of the poles, compared to the equator, and thus a reduced meridional temperature gradient.
The parameter is altered between −4K and +4K (with steps of 4K) relative to the climatological present-day (PD) temperature (1979 -2014). This range covers climate variability over the past million years and possible near future changes.
For each temperature component we determine its influence on the strength, and its width of the Hadley cell, as well as the 15 strength of zonal-mean jets, storm tracks and planetary waves in the Northern Hemisphere.
Since the changes for the zonal temperature asymmetry and the meridional temperature gradients are relative, it can be that in terms of absolute values the zonal temperature asymmetry is greater than that of the meridional temperature gradient.
To force the Aeolus dynamical core, we use perturbed surface temperature profiles derived from the ERA-Interim winter climatology as explained above. We performed 2025 simulations with a regular 3-dimensional parameter space using the 20 multi-run simulation environment SimEnv, which provides a tool to inspect the model's behaviour in the parameter space by discrete numerical sampling (Flechsig et al., 2013).

Dynamical variables
To obtain the strength of the jet stream for this analysis, we use seasonally (DJF) averaged zonal mean and vertical averaged 25 from 3000m to 9000m (corresponding to 200mb to 700mb) zonal wind 〈 〉 ̅̅̅̅ . For simplicity, we define the jet stream strength as the maximum of 〈 〉 ̅̅̅̅ between 10°N and 80°N of the vertical averaged wind velocity We define the strength of the Hadley cell as the maximum of the meridional mass stream function integrated between 200mb -700 mb as well as between 0° and 90° latitude and the width as the distance between the mass flux zero-crossings near 0° and 30° latitude.
As a measure of storm track activity we calculate the eddy kinetic energy ( ′ = 0.5( ′2 ± ′2 )) , whereby ′ and ′ are the 5 zonal and meridional synoptic wind velocity. We average the eddy kinetic energy over all five pressure levels and calculate the maximum between 10° N and 80° N to analyse the strength and shift of the storm track activity.
We calculate the strength of the planetary waves by averaging all positive values between 20°N−80°N of the azonal wind components 〈 * 〉 and 〈 * 〉.

10
We compare and analyze the zonal mean dynamical variables of eddy kinetic energy 〈 ′ 〉 ̅̅̅̅̅̅ (which captures storm track activity), zonal mean zonal wind velocity 〈 〉 ̅̅̅̅ and the vertical integral of the lower tropospheric integrated mass flux 〈 〉 as well as azonal wind velocities 〈 * 〉 and 〈 * 〉.

Strength and width of the Hadley cell
The integrated mass flux in the lower troposphere of the present-day modelled climatological NH winter values (black line in Fig. 2) captures well the shape of the red curve from ERA-Interim data. In particular, the maximum strength, defined as the minimum between the zero-crossings, is close to the ERA-Interim data. There exist bigger differences in the SH. This model 20 bias might be related to the missing Antarctic ice sheet, upper-tropospheric ozone, the constant lapse rate assumption, or fundamental limitations of the equations. Another reason could be due to the statistical nature of the eddy representation in the SDAM since , the summer HC is dominated by eddy momentum flux divergence (Schneider and Bordoni, 2008).
Zero-crossing refers to the point in the graph where the function ( ) crosses the = 0 line. The modelled Hadley cell's width, defined as the distance between the mass flux zero-crossings near 0° and 30° latitude, is smaller than in ERA-Interim. 25 For further analysis we plot the width (Fig. 3) and strength (Fig. 4)  is the difference between the present day global mean temperature and the altered global mean temperature.
In general, both a stronger meridional temperature gradient and a stronger zonal temperature contrast lead to a nonlinear broadening of the Hadley cell width: For meridional temperature gradients smaller than today, the width is smaller for larger global mean temperature. The Hadley cell expands as well for stronger meridional temperature gradient and even for a greater zonal temperature asymmetry, but with a smaller rate.
For meridional temperature gradients larger than today, the influence of global mean temperature and meridional temperature gradient are less significant, and the values of the zonal temperature asymmetry play a bigger role.In this case larger values of 5 the zonal temperature asymmetry lead to a larger width of the Hadley cell. However, depending on the global mean temperature an increase of the zonal temperature asymmetry can also lead to a decrease of the Hadley cell width, e.g. , ≥ 3, 0.95 ≥ ≥ 0.9 and ≥ 0.85.
The Hadley cell width shows larger changes in response to changes in the meridional temperature gradient than for changes in the zonal temperature asymmetry, indicating that the former has a stronger relative influence. 10 The Hadley cell strengthens with increasing meridional gradient and depends stronger on global mean temperature than on (Fig. 4).

Strength of the jet stream 15
The jet stream locations and strengths for both hemispheres are detectable as two distinct maxima in the zonal-mean zonal wind velocity in ERA-Interim (Fig. 5a). Aeolus reasonably reproduces the main jet stream features in terms of spatial position and magnitude (Fig. 5c). The modelled magnitude of the jet in the Northern Hemisphere is in better agreement with reanalysis data than in the Southern Hemisphere. This is likely related to the fact that Aeolus is not coupled to an ice model and thus effects from the Antarctic ice sheet are not considered. The model reasonably reproduces near-surface tropical easterlies ("trade 20 winds") at low latitude. For a decrease of 10 % of the temperature gradient the velocity weakens to 20 m/s. With an increase of 10 % of the meridional 25 temperature, the jet stream velocity is more than 30 m/s and with 10 per cent stronger in the Northern Hemisphere. This is also observed in Figure 6 where the jet stream strength is shown as function of , and , .
The strength of the jet stream is sensitive to the meridional temperature gradient and to the global mean temperature. The zonal temperature contrasts have little influence on the jet stream strength.

Strength of the storm track activity
The NH winter climatological (1979 -2014) storm track's activity in Aeolus (Fig. 7(c)) is similar to ERA-Interim data ( Fig.   7(a)).

Figures 7(b)-(d) show that storm track activity increases with increasing temperature gradient. 5
The strength of the storm track activity depends on all three components , and , (Fig. 8) in a way that the influence of dominates the influence on storm track activity. The increased global mean temperature leads to a general strengthening of the storm track activity (Fig. 8).

Wind velocity of the planetary waves 10
The wind velocity of the planetary waves is roughly as sensitive to as to , both in terms of 〈 * 〉 ( Fig. 9(a)) and in terms of 〈 * 〉 ( Fig. 9(b)). Both meridional and zonal wind directions exhibit the same relationship such that larger meridional and zonal temperature asymmetries lead to stronger winds. In addition, if the scaling factor of the meridional temperature gradient is smaller than the scaling factor of the zonal temperature asymmetry, the strength of the planetary waves increases faster with higher relative changes of the meridional temperature gradient then if both have similar relative changes 15 of magnitude. Moreover, if the relative change of the zonal temperature asymmetry is smaller than the relative changes of the meridional temperature gradient, the strength of the planetary waves increases faster with higher zonal temperature asymmetry then if both have a similar magnitude. This behaviour leads to a curved structure of the azonal winds.
The global mean temperature has a positive but only weak influence on the strength of planetary waves. 20

Discussion
For all investigated atmosphere variables we observe a strengthening for higher global mean temperature and higher absolute meridional temperature gradients and only a weak (strong) dependence on the zonal temperature asymmetry for storm tracks (planetary waves and Hadley cell width), which we discuss in comparison with results from literature in the following sections. 25 However, most previous studies have analysed only the combined effect of changes in several temperature components making a direct comparison difficult.

Strength and width of the Hadley cell
The strength of the Hadley cell depends strongly on the meridional temperature gradient with a stronger Hadley circulation for larger meridional gradient (Fig.4). Its strength is much less sensitive to global mean temperature. A rising global mean 5 temperature leads to a weakening of the Hadley cell explainable by the enhanced latent heat release under warmer conditions. The specific humidity increases faster with temperature than precipitation and therefore the strength weakens (Hill et al., 2015).
Finally, it is almost insensitive to zonal temperature asymmetries.
In addition, our analysis suggests that the Hadley cell width depends nonlinearly on all three temperature components (Fig. 3). In agreement to our results, Mitas and Clement (2005) detected a strengthening of the Hadley cell in their analyses using several reanalysis data sets, a rawsonde data set and a model data set. However, they found great differences between different 25 data sets. Lu et al. (2007) found a robust weakening and a poleward expansion of the Hadley circulation in response to increased GHG forcing in simulations of the 21st century climate taken from the A2 scenario of the IPCC AR4 project. Lu et al. (2008) analysed the change in the zonal mean atmospheric circulation under global warming in comparison with the response to El Niño forcing, by examining the CMIP5 model simulations. They used again the A2 scenario to simulate global warming. 30 Under global warming due to higher CO2 concentrations the Hadley cell weakens and expands northwards together with a poleward shift of the jet stream. Based on our results, we can assume that "El Niño-like" enhanced warming leads to a stronger zonal temperature asymmetry (and a higher global mean temperature) resulting in a stronger Hadley cell, whereas the CO2 concentration leads to a weaker meridional temperature gradient (and a higher global mean temperature) and as a consequence the Hadley cell weakens. This can also explain the widening of the Hadley cell, which we observe in our experiments as well: A decreased meridional temperature gradient for warmer global mean temperature than today can lead to a smaller width of the Hadley cell and vice versa. Seo et al. (2014) investigated possible drivers of the Hadley cell such as the meridional temperature gradient, gross static 5 stability and tropopause height using CMIP5 climate models. Consistent with our results, they found a robust dependence between meridional temperature gradient and the strength of the Hadley cell in winter: A decreased meridional temperature gradient leads to a weakening of the Hadley cell.
In addition, D' Agostino et al. (2017) analysed and compared the Hadley cell during the last glacial maximum to global warming scenarios (RCP4.5 and RCP8.5) with a focus on the dependence on subtropical stability, near-surface meridional 10 potential temperature gradient, and the tropical tropopause level. They concluded that the meridional temperature gradient is a major driver for Hadley cell changes.
However, in both studies the atmospheric composition in terms of anthropogenic aerosols is changed and hence not only the meridional temperature gradient changes but also the global mean temperature and the zonal temperature asymmetry. This makes it difficult to attribute changes in the Hadley cell to one temperature component only. 15 Levine and Schneider (2015) investigated the sensitivity of the Hadley cell width to a equator-to-pole temperature contrasts and global-mean temperatures using two idealized GCMs.Their results confirm our findings that a an increasing equator-topole temperature leads to a widening of the Hadley cell (Levine and Schneider, 2015).

Strength of the jet stream 20
We show that the strength of the jet stream decreases with decreasing absolute meridional temperature gradient ( Fig. 5 and   Fig. 6). This is in agreement with findings from Polvani and Kushner (2002) and Haigh et al. (2005). Polvani and Kushner used a simple general circulation model and showed that for sufficiently strong cooling of the polar winter stratosphere, jet streams 25 weaken and shift poleward. Haigh et al. (2005) analyse the weakening and shift of the subtropical jet using a multiple regression analysis of the NCEP-NCAR reanalysis zonal mean zonal wind velocity. Furthermore, they show with a simple general circulation model that the generic heating of the lower stratosphere tends to weaken the subtropical jets.
In most observational studies, a weakening of the jet is observed over the last decades like Archer and Caldeira (2008) using NCEP and ERA-40 reanalysis sets, Rikus (2015) using ERA-40 data and using ERA-Interim data . 30 However, Pena-Ortiz et al. (2013) found that trends in both strength and position of the jet strongly vary between different reanalysis products.

Strength of the storm track activity
In this study, we observe a strengthening of storm track activity under increased global-mean temperature (Fig. 8). 5 Our results are supported by findings from McCabe et al. (2001), who observe a strengthening of the storm track activity with higher global temperature induced by GHG forcing. This is also in agreement with Yin (2005), who investigated 15 coupled climate models and showed that storm tracks intensify under global warming. In addition, Chang et al. (2012) found that storm tracks in the upper troposphere increase in winter using 23 CMIP5 models (below 300mb they found a slight decrease) and the strength of storm track activity depends strongly on the meridional temperature gradient, which is consistent to our findings. 10 This latter result is intuitive as the prime role of storm tracks within the general circulation is to transport heat poleward, with a stronger temperature contrast leading to enhanced heat transport. It also directly follows from the equation of eddy kinetic energy, which in the first place depends on the meridional temperature gradient (Coumou et al., 2011). Harvey et al. (2013 observe similar results using CMIP5 data: Larger temperature differences in the equator-to-pole temperature at upper-and lower-tropospheric levels lead to stronger storm activity . 15 In reanalysis data also a strengthening of the storm tracks can be observed (Schneidereit et al., 2007;Wang et al., 2006), which is supposedly because of the rising global mean temperature. The zonal temperature asymmetry could not be responsible for a strengthening, since the zonal temperature asymmetry should be reduced in winter due to global warming. According to our results, this would lead to a weakening of the storm track activity.
O'Gorman and  examined the response of storm tracks to different climate conditions simulating an 20 aquaplanet and by changing the longwave optical thickness in the radiation scheme of the GCM (representing variations in greenhouse gas concentrations). They found that eddy kinetic energy has a maximum for a climate with the global-mean temperature similar to that of present-day-climate. Lower or higher global-mean temperatures lead to significantly smaller values. In addition, they observed that the eddy kinetic energy increases monotonically with the meridional insolation gradient (representing changes in, for example, high-latitude surface albedo). 25 Similarly, Pfahl et al. (2015) investigated the behaviour of extratropical cyclones under strongly varying climate conditions using idealized climate model simulations in an aquaplanet setup. They changed the meridional insolation gradient together with the longwave optical thickness with shortwave parameters held constant. They found that the maximum of eddy kinetic energy is reached at a global mean temperature slightly warmer than present-day climate.
These results are different to our results, where no such peak in Eddy kinetic energy can be observed. The different results 30 may be explained by the different techniques applied to simulate higher global mean temperature. In our study, we directly change the temperature, whereas Pfahl et al. change the longwave optical thickness with shortwave parameters held constant, which represents variations in longwave absorbers like carbon dioxide and water vapour. These changes could also change the meridional and zonal temperature asymmetry leading to different results.
Nevertheless, we also observe a strong positive dependence between temperature gradient and Eddy kinetic energy.

Wind velocity of planetary waves 5
In our analysis the strengthening of the planetary wave wind velocities depends on all temperature components (Fig. 9). Larger meridional and zonal temperature asymmetries as well as global mean temperatures lead to stronger winds.
Since azonal wind components emerge due to zonal temperature asymmetries, it is expected that higher zonal temperature differences lead to stronger azonal wind components. Stronger temperature gradients cause stronger meridional wind velocities, which are deflected by the Coriolis force and therefore also the zonal wind velocities are stronger. Those wind 10 velocities are slowed down or accelerated due to topography, land-ocean-contrast and hence also the azonal component of the wind velocity will be stronger. In addition, a higher global mean temperature leads to more available energy in the atmosphere and therefore to a larger azonal wind velocity.
Under climate change the global mean temperature increases whereas the meridional temperature is expected to decrease. Our results suggest that this will have contrary effects on the strength of planetary waves. 15 Thus, this could explain the results, which Barnes and Polvani (2013) found by analysing the planetary waves with wave number 1 -6 as well as wave numbers 1 -3. She concluded that there is no significant trend in terms of the strength of the planetary waves and thus Arctic amplification does not play a dominant role for changing the undulations of the jet stream.

20
The large-scale dynamical atmospheric changes due to variations of the temperature components, as presented here, coincide with those expected from dynamical principles, with the possible exception of the Hadley cell width changes. This therefore first of all serves as an important validation of our SDAM atmosphere model. To our knowledge, no other SDAM model exists that captures these dynamical interactions between the large-scale circulation components of tropical circulation, jets, storm tracks and planetary waves. 25 The fine sampling resolution of the parameter space reveals that most variables have a linear relationship to the investigated temperature components. In retrospect, a coarser resolution sampling of the parameter space using a GCM might have sufficed to detect the relations between temperature components and dynamical variables.
However, only through our many runs using the fast SDAM this linearity could be shown.

Conclusion
In this paper, we present a study on multiple fundamental components of the large-scale atmosphere dynamics to different surface temperature forcing with the statistical-dynamical Atmosphere model Aeolus 1.0. Due to the statistical-dynamical approach, Aeolus 1.0 is much faster than GCMs, which allows us to do 1000s of individual simulations and thus test the 5 sensitivity of the dynamical fields to different surface temperature changes. This way one can disentangle and separately analyse the effect of global mean temperature, equator-to-pole temperature gradient and east-west temperature differences.
Therefore, we are one of the first, who scan the full temperature phase space. This way we can scan for 'non-linearities' in the system (i.e. the Hadley cell edge might be very sensitive to meridional temperature gradients only for a narrow range of temperature gradient values, and outside of that range it is not sensitive) as done for example in Levine and Schneider (2015). 10 It is important to know such non-linearities as it could imply more abrupt changes under global warming.
The model's climatology generally reproduces the dynamical fields of ERA-Interim, especially in the Northern Hemisphere, which is the focus of our analysis. If possible, we compare our findings with results of the literature and conclude that most modelled changes are in line with theory and simulations. Exceptions are the planetary waves, whichare well explained by theoretic dynamical consideration and the width of the Hadley, which could be a model artefact or a real feature. Latter should 15 be tested with GCMs.
These results also serve as an important validation of the dynamical core of the Aeolus. We could show that Aeolus is to our knowledge the first model that captures the dynamical interactions expected from dynamical principles between the large-scale circulation components of tropical circulation, jets, storm tracks and planetary waves. In future work we would like to use the gained knowledge to simulate only specific temperature component configurations to further explore the dependence of the 20 different atmospheric large-scale circulations on near-surface temperature profiles.

Code and data availability
All original data was downloaded from public archives. Code and data are stored in PIK's long term archive, and are made available to interested parties on request.    .
All other parameters are set to climatology. With larger gradients the storm track activity gets stronger. 10