1. Introduction
Since the middle of the twentieth-century orography has been known to have a substantial impact on large-scale circulation and stationary waves (Charney and Eliassen 1949; Bolin 1950). Orography subsequently affects regional climates around the world, as reflected in the climatology of surface temperatures (Seager et al. 2002) and precipitation (Broccoli and Manabe 1992; Wills and Schneider 2015). Traditionally the North American Cordillera (or Rocky Mountains) and the Tibetan Plateau and Himalaya have been considered the most important topographic features influencing atmospheric flow (e.g., Held et al. 2002). However, recent work suggests that the smaller Mongolian plateau to the north of the Himalaya is the most important region of Northern Hemisphere (NH) orography during boreal winter, due to its location within the latitude of strong near-surface zonal winds (White et al. 2017, 2018).
Earth’s atmospheric stationary waves are manifest in the planetary-scale wavelike structures in climatological sea level pressure or geopotential height. In the NH the dominant features of the wintertime stationary waves are so prominent they have been given names: the Aleutian low, the Icelandic low, and the Siberian high. The forcing of stationary waves is mainly by orography (Charney and Eliassen 1949; Bolin 1950) and zonal asymmetries in atmospheric diabatic heating associated with land–sea temperature contrasts (e.g., Smagorinsky 1953; Manabe and Terpstra 1974; Hoskins and Karoly 1981; Valdes and Hoskins 1991). In addition, orography can also influence diabatic heating, affecting the stationary waves (Chang 2009). Surface heat fluxes associated with ocean currents can also play a modifying role (Garfinkel et al. 2020).
From the 1980s to 2010 numerous studies examined the relative contributions of “orographic” and “thermal” forcing to the observed stationary waves, attempting to quantify the role of orography. This research, collectively using linear and nonlinear stationary wave models, as well as full general circulation models (GCMs), finds that orography produces anywhere between 30% and 66% of the observed stationary wave amplitude (Held 1983; Chen and Trenberth 1988a,b; Nigam et al. 1988; Valdes and Hoskins 1989), with recent estimates toward the lower end of this range (Ting et al. 2001; Held et al. 2002; Chang 2009). The large range of results can be at least partially understood from differences in the background flow into which the orography was placed, including differences in low-level winds (Valdes and Hoskins 1989; Held and Ting 1990; Ringler and Cook 1995; Held et al. 2002), as well as model differences in atmospheric dissipation strength and imposed drag (Valdes and Hoskins 1989). Additionally, as suggested by Held et al. (2002) and shown by Chang (2009), orography impacts the patterns of atmospheric diabatic heating, and thus a separation of the stationary waves into contributions from “orography” and “thermal” forcing is not a well-posed problem. Indeed, using an intermediate complexity GCM Garfinkel et al. (2020) show that such strong nonlinearities exist between different forcing mechanisms of the stationary waves, that, over the Pacific, the linear addition of the sources gives a response of different sign to the full nonlinear result. Garfinkel et al. (2020) also show that zonal asymmetries in ocean heat fluxes play a small, but not insignificant, role through their influence on the distribution of sea surface temperatures (SSTs). Using a comprehensive GCM Brayshaw et al. (2008, 2009) explored the influence of idealized landmasses (e.g., rectangular continents), orography (e.g., Gaussian shaped mountains), and SST gradients. In this study we focus only on the influence of present-day orography versus land–sea contrast.
In comparison with stationary waves, less attention has been paid to the impact of orography on the axisymmetric (zonal mean) flow. Held et al. (2002) notes that the structure of the tropospheric zonal mean zonal wind
Multiple mechanisms exist for how orography impacts the stationary waves and axisymmetric flow. Orography acts as a large-scale obstruction to the flow, creating stationary waves (Charney and Eliassen 1949; Bolin 1950); these stationary waves can also influence the axisymmetric flow through wave–mean-flow interactions (Eliassen and Palm 1961; Charney and Drazin 1961; Andrews and McIntyre 1976; Edmon et al. 1980). Other impacts of orography on the circulation that are at least partially resolved in models include orographic pressure torque, created by gradients in pressure between the western and eastern flanks of mountains (e.g., Wahr and Oort 1984), and differential heating (both radiative and latent) caused by orography (e.g., Ringler and Cook 1995; Wills and Schneider 2018). Orography also affects circulation through sub-grid-scale impacts that must be parameterized at the resolution of past and current GCMs (e.g., Wallace et al. 1983; Boer et al. 1984). Sandu et al. (2016) give a thorough description of the sub-grid-scale effects of orography and their typical representation in atmospheric models. In short, parameterized physics typically include: the forcing of gravity waves, which break close to the surface as well as in the upper atmosphere affecting the circulation (e.g., Boer et al. 1984; McFarlane 1987; Miller et al. 1989; Shaw et al. 2009) and the deceleration of the near-surface flow by turbulent or frictional small-scale processes, including sub-grid-scale blocking of flow (e.g., Webster et al. 2003; Pithan et al. 2016; Sandu et al. 2019). Parameterizing the sub-grid-scale effects of orographic affects both atmospheric
We perform a comprehensive series of experiments using NCAR’s Whole Atmosphere Community Climate Model, version 6 (WACCM6), in which global mountains are removed, with and without the simultaneous removal of sub-grid-scale orography. Setting the sub-grid-scale orography to zero removes the effects of unresolved orographic drag parameterized in WACCM6, including gravity waves, nonlinear surface drag, and turbulent form drag; hereafter, we refer to this as the response to “parameterized orographic drag,” or, for brevity, “parameterized orography.” We quantify the role of both resolved and parameterized orography on the wintertime (December–February in the NH) axisymmetric flow and stationary waves. We briefly discuss the sensitivity to the treatment of sea surface temperatures in section 4. Through analysis of additional mountain/no-mountain experiments in simulations in which Earth is rotated backward, we highlight the strong sensitivity of the orographic response to rotation direction, showing that the current distribution of orography on our planet relative to the coastlines seems to produce a particularly strong orographic response.
2. Models and experiments
We use the WACCM6, a dynamical atmosphere model with 70 vertical levels, including a well-resolved stratosphere (Gettelman et al. 2019). In addition to high vertical resolution in the upper atmosphere, the WACCM has sophisticated parameterizations of the drag from unresolved orography (Gettelman et al. 2019). Effects of unresolved orography are parameterized by two separate schemes: a turbulent orographic form drag (TOFD) scheme, which parameterizes the boundary layer form drag (Beljaars et al. 2004), and a scheme, which parameterizes the effects of gravity wave drag on both the near-surface and upper atmosphere, and incorporates near-surface nonlinear drag processes (Scinocca and McFarlane 2000). These parameterizations, although sophisticated, remain estimates of the true impact of unresolved orography. In addition to gravity waves forced by sub-grid-scale orography, the WACCM6 includes parameterization of frontal and convectively generated gravity waves (Gettelman et al. 2019; Richter et al. 2010).
For our main simulations a horizontal resolution of 0.9° × 1.25° (hereafter ~1°) is used. All simulations have atmospheric constituent values of the present day (circa year 2000), and the SSTs are fixed to a seasonally varying observed climatology from the late twentieth century (Hurrell et al. 2008). The WACCM6 model reproduces tropospheric and stratospheric circulations well, including stratospheric variability in the form of sudden stratospheric warmings (Gettelman et al. 2019). Throughout this paper we further evaluate the WACCM by comparing it with the ERA-Interim data (ERA-I; Dee et al. 2011).
To investigate the role of orography in shaping the zonal mean flow and stationary waves four WACCM6 simulations are performed: CTL, no mountain (NM), CTLsmooth, and NMrough. The CTL simulation has present-day topography, while the NM simulation has all topography flattened to sea level. In addition to flattening the resolved orography in the NM experiment, we reduce all boundary condition variables associated with sub-grid-scale orographic variability to 0 (the value over the ocean), effectively switching off all parameterized responses to unresolved orography, including the turbulent orographic form drag, orographic gravity wave drag, and the nonlinear near-surface drag. The CTLsmooth experiment has smooth mountains, i.e., we keep the resolved orography of CTL, but use the sub-grid-scale orographic variability of NM; the NMrough simulation is the opposite to CTLsmooth, i.e., zero resolved orography but the same sub-grid-scale orographic variability as CTL. The simulations are 20 years in length (following a 1-yr spinup period). All responses analyzed in this study are large relative to the interannual variability, and so 20 years is found to be sufficient to obtain a robust signal (not shown). Significance testing confirms that this is true even in the stratosphere, where interannual variability is relatively large. We run sensitivity experiments at 1.9° × 2.5° (hereafter ~2°) to test the sensitivity of our results to horizontal resolution. We also run sensitivity experiments using an earlier version of the model, the WACCM4, which uses different parameterizations of the unresolved orography. Our conclusions are unchanged by these changes in resolution or model version (results presented in the online supplemental material).
Orographic forcing can also influence SSTs through air–sea heat fluxes (Okajima and Xie 2007; Chang 2009); this may affect the impact of orography on diabatic heating, although historically this effect has been considered to be small (Held 1983). Additionally, changes in SSTs may alter circulation patterns and thus the winds impinging on orography; this can affect the atmospheric response to orography. We repeat some of our key simulations with the atmosphere model coupled to a slab ocean, which allows the SSTs to adjust to changes in air–sea heat fluxes; this has little impact on our conclusions. These results are not shown, but are discussed briefly in section 4. Mountains are also known to affect ocean circulation (Kitoh 2002; Sinha et al. 2012) through impacts on wind stress and freshwater fluxes (Warren 1983; Emile-Geay 2003), which can further affect surface heat fluxes. Here we focus on the dynamical atmospheric response to global mountain ranges and do not consider dynamical ocean changes.
3. Results
In section 3a we present the effects of orography on the wintertime axisymmetric flow, including a separation of the impacts of resolved orography and parameterized orographic drag, and compare the relative importance of different mechanisms. Following this, in section 3b we present the impacts of orography on wintertime stationary waves.
a. Zonal mean response to orography
Orography substantially decelerates NH wintertime zonal mean zonal wind

Wintertime (DJF in the NH and JJA in the SH) zonal mean zonal wind
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

Wintertime (DJF in the NH and JJA in the SH) zonal mean zonal wind
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
Wintertime (DJF in the NH and JJA in the SH) zonal mean zonal wind
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
Figure 2 shows NH

(a)–(c) DJF
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

(a)–(c) DJF
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
(a)–(c) DJF
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
In the free troposphere there is little to no contribution from gravity wave drag in the CTL simulation (Fig. 2d), as expected, and a deceleration from resolved EP convergence (Fig. 2g), consistent with previous research (Edmon et al. 1980; Li et al. 2011). In the stratosphere there is generally deceleration from both orographic gravity wave forcing and resolved waves, with stronger forcing from the resolved waves (cf. Figs. 2d,g). In the NM simulation (second row), there is no orographic gravity wave drag, and the deceleration from EP flux convergence is reduced in the stratosphere. In the high-latitude stratosphere, there is divergence of EP fluxes in the NM simulation, suggesting a local wave source and acceleration of the zonal flow by wave activity; this is consistent with observations of the Southern Hemisphere (e.g., Hartmann et al. 1984; Randel and Lait 1991).
Figure 2c shows that orography greatly slows down
Other impacts of parameterized orographic drag, including the drag from the TOFD parameterization and nonlinear processes near the surface, likely also play a significant role in the low-level

As in Fig. 2, but for (top) CTL − NM (repeated from bottom row of Fig. 2), (middle) the response to resolved orography (CTLsmooth − NM), and (bottom) the response to parameterized orography (CTL − CTL smooth).
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

As in Fig. 2, but for (top) CTL − NM (repeated from bottom row of Fig. 2), (middle) the response to resolved orography (CTLsmooth − NM), and (bottom) the response to parameterized orography (CTL − CTL smooth).
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
As in Fig. 2, but for (top) CTL − NM (repeated from bottom row of Fig. 2), (middle) the response to resolved orography (CTLsmooth − NM), and (bottom) the response to parameterized orography (CTL − CTL smooth).
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
Resolved orographic pressure torque contributes to the resolved response to orography. This torque is estimated for the CTL simulation using the surface pressure method of Wahr and Oort (1984): summing a value proportional to ps × 0.5(Hi+1 − Hi−1) over all grid boxes i along a circle of constant latitude, where ps is surface pressure and H is surface height. This gives values in quantitative agreement with observations from Tucker (1960) and Madden and Speth (1995) for the CTL simulation. Between 40° and 70°N, where orography has the largest impact on
Resolved orography versus parameterized orographic drag
We now decompose the total response to orography into that forced by the large-scale resolved orography (CTLsmooth − NM), and that forced by parameterized orographic drag (CTL − CTLsmooth). The decomposition can be seen in Fig. 3, where the first row repeats the CTL − NM fields shown in the bottom row of Fig. 2, the middle row shows the impact of resolved orography, and the bottom row that of parameterized orographic drag. The total response to parameterized orography is a combination of the parameterized gravity wave response (Fig. 3f), parameterized low-level drag (TOFD and nonlinear drag from the orographic gravity wave scheme; not shown), as well as subsequent changes in resolved wave forcing (Fig. 3i). Parameterized orography has a much greater influence on
Parameterized orography dominates both changes in orographic gravity wave drag (as expected) as well as changes in resolved wave forcing, at least away from the surface (cf. Figs. 3d–i). Resolved orography acts as a strong source of resolved wave activity at the surface around 40°–60°N, where strong surface westerlies interact with the major NH mountain ranges (White et al. 2017); this surface wave source can be seen in the upward EP flux arrows emanating from the bottom of Fig. 3h. Most of this wave activity, however, converges below 500 hPa, with relatively little able to propagate into the stratosphere to converge and decelerate stratospheric
b. Stationary wave response to orography
We analyze stationary waves through the climatological eddy geopotential height:

DJF NH eddy geopotential height Z′ (in m) at 45°N for (a) ERA-I, (b) CTL, (c) NM, and (d) CTL − NM.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

DJF NH eddy geopotential height Z′ (in m) at 45°N for (a) ERA-I, (b) CTL, (c) NM, and (d) CTL − NM.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
DJF NH eddy geopotential height Z′ (in m) at 45°N for (a) ERA-I, (b) CTL, (c) NM, and (d) CTL − NM.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
The spatial distribution of the response to all orography in Fig. 4d is broadly consistent with previous studies (e.g., Held 1983). Conversely, the spatial pattern of the NM stationary wave pattern, shown in Fig. 4c, shows some differences, including a higher wavenumber, than in some previous studies (e.g., Held 1983; Nigam et al. 1986). As will be discussed in section 4, the NM stationary wave pattern is very sensitive to small differences in tropical SSTs, and thus differences between NM simulations in different studies may be due to differences in tropical SST fields.
To study the horizontal structure of the stationary wave response to orography, Fig. 5 shows maps of Z at 1000, 300, and 30 hPa, with contours showing the full fields and colors showing the eddy component, Z′. The WACCM reproduces the ERA-I spatial pattern of stationary waves well (cf. Fig. 5 top row with the middle row). The midlatitude stationary wave pattern is much weaker in NM than in CTL throughout the troposphere and stratosphere. Near the surface, the Z′ pattern in CTL is predominantly due to the presence of orography: without mountains the near-surface stationary wave has little resemblance to the CTL stationary wave (cf. Figs. 5b,c). In the upper troposphere the NM eddy field (Fig. 5f) exhibits a spatial structure similar to that of CTL (Fig. 5e), but with muted amplitude. This suggests a constructive interference between the NM and orographically forced stationary waves, noted by Held (1983), Held et al. (2002), and Chang (2009). This can also be seen in the cross section of Z′ in Fig. 4, where the free troposphere response to orography (Fig. 4d) generally projects onto the NM stationary wave pattern (Fig. 4c), with a strong low centered around 150°E and highs at around 130° and 30°W. This constructive interference, which is sensitive to the tropical SSTs specified in the NM case (see section 4), helps augment the magnitude of the free troposphere stationary wave in the CTL simulation, particularly the low in the western Pacific (120°E–180°) and the high over the Atlantic and eastern Europe (60°–0°W). The NM stratospheric stationary wave is weak and has little resemblance to that in the CTL simulation (cf. Fig. 4c with Figs. 4d and 5g–i).

DJF NH geopotential height Z (in m), full fields (contours), and eddy component Z′ (colors) at (a)–(c) 1000, (d)–(f) 300, and (g)–(i) 30 hPa for (top) ERA-I, (middle) CTL, and (bottom) NM. Prominent highs and lows in the full surface fields are denoted with H and L, respectively, in (a)–(c). The contour interval for the full fields is 50, 200, and 250 m at 1000, 300, and 30 hPa, respectively.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

DJF NH geopotential height Z (in m), full fields (contours), and eddy component Z′ (colors) at (a)–(c) 1000, (d)–(f) 300, and (g)–(i) 30 hPa for (top) ERA-I, (middle) CTL, and (bottom) NM. Prominent highs and lows in the full surface fields are denoted with H and L, respectively, in (a)–(c). The contour interval for the full fields is 50, 200, and 250 m at 1000, 300, and 30 hPa, respectively.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
DJF NH geopotential height Z (in m), full fields (contours), and eddy component Z′ (colors) at (a)–(c) 1000, (d)–(f) 300, and (g)–(i) 30 hPa for (top) ERA-I, (middle) CTL, and (bottom) NM. Prominent highs and lows in the full surface fields are denoted with H and L, respectively, in (a)–(c). The contour interval for the full fields is 50, 200, and 250 m at 1000, 300, and 30 hPa, respectively.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
Resolved orography versus parameterized orographic drag
The stationary wave response is now separated into the response to resolved and parameterized orography using the CTLsmooth simulation (Fig. 6). Near the surface, both resolved and parameterized orography play approximately equal roles in the total stationary wave response (cf. Figs. 6b,c). In the upper troposphere and stratosphere the response to resolved orography is slightly stronger than the response to parameterized orography (cf. Fig. 6e with Figs. 6f,h with Fig. 6i). The total response to orography is spatially much more similar to the response to resolved orography (cf. Fig. 6d with Figs. 6e,g with Fig. 6h), than to the parameterized orography (cf. Fig. 6d with Figs. 6f,g with Fig. 6i), suggesting the response to resolved orography is dominant. This conclusion holds when we use the NMrough experiment to separate resolved and parameterized orography, although the stationary wave response to parameterized orography is much stronger when resolved orography is present (i.e., CTL − CTLsmooth) than when resolved orography is absent (i.e., NMrough − NM).

As in Fig. 5, but for the (top) total orographic response (CTL − NM), (middle) resolved orography response (CTLsmooth − NM) and (bottom) response to parameterized orography (CTL − CTL smooth). The contour interval for the full fields is 50, 50, and 200 m at 1000, 300, and 30 hPa, respectively.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

As in Fig. 5, but for the (top) total orographic response (CTL − NM), (middle) resolved orography response (CTLsmooth − NM) and (bottom) response to parameterized orography (CTL − CTL smooth). The contour interval for the full fields is 50, 50, and 200 m at 1000, 300, and 30 hPa, respectively.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
As in Fig. 5, but for the (top) total orographic response (CTL − NM), (middle) resolved orography response (CTLsmooth − NM) and (bottom) response to parameterized orography (CTL − CTL smooth). The contour interval for the full fields is 50, 50, and 200 m at 1000, 300, and 30 hPa, respectively.
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
4. Discussion
a. A proposed mechanism for the impacts of orography on stratospheric
That parameterized orography impacts stratospheric

Schematic showing the influence of orography on DJF NH stratospheric
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1

Schematic showing the influence of orography on DJF NH stratospheric
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
Schematic showing the influence of orography on DJF NH stratospheric
Citation: Journal of the Atmospheric Sciences 78, 7; 10.1175/JAS-D-20-0300.1
It is notable that resolved wave convergence in the stratosphere is more impacted by the presence of parameterized orography than by the presence of resolved orography (see Fig. 3 and Fig. S3). Resolved orography acts as a surface source of resolved waves, producing a small increase of EP flux propagation into the stratosphere, which results in a small reduction of
b. Impact of SSTs
All of the simulations presented in this paper use fixed, present-day, SSTs. We investigated the sensitivity of our results to this choice of lower boundary condition by running simulations with the atmosphere model coupled to a slab-ocean model, thus allowing the SSTs to change in response to changing surface heat fluxes, but not considering ocean circulation changes. Through a series of experiments in which we limit changes in SST to the tropics and/or extratropics, we find that the extratropical SST changes caused by the presence of orography have little impact on the results and conclusions presented. Both the strength and the spatial pattern of stationary waves in the NM simulations, however, are remarkably sensitive to relatively small (~1 K) changes in the tropical SSTs. This is consistent with previous work finding a strong sensitivity of midlatitude circulation and stationary waves to tropical SSTs (e.g., Yin and Battisti 2001; Held et al. 2002; Inatsu et al. 2002), although our results suggest that in the real world orography reduces the sensitivity somewhat. Tropical SST changes induced by substantial shifts in ocean circulation due to the presence/absence of orography could affect our conclusions.
c. The importance of the location and orientation of orography
To further understand the results presented in this study, and test how sensitive the atmospheric response to orography is to the relative distribution and orientation of orography and land boundaries on present-day Earth, we also run simulations at ~2° horizontal resolution with a backward-rotating Earth, CTLBW and NMBW. We set Earth’s angular rotation rate ωBW = −1.0ω. These simulations have the same imposed SSTs as the standard simulations, including the effects of western boundary currents, thereby ignoring changes in air–sea fluxes and ocean circulation from the reversed rotation direction. Sensitivity experiments using slab oceans with zonal mean q fluxes (i.e., assuming no longitudinal dependence of oceanic heat transport) in both forward and backward simulations find that our conclusions are insensitive to this choice of SSTs. In these retrograde-rotating Earth simulations, there is easterly flow in the midlatitudes, and thus the location of the most prominent orography relative to the boundaries where air flows from ocean to land is changed substantially from the real world, as well as substantial changes in the orientation of mountain ranges relative to the impinging winds. The results are summarized in the supplementary material; see Figs. S4 and S5. We use these experiments to highlight the importance of the orography–land configuration on the effect of orography on the atmospheric circulation. For a study of the climate of a retrograde-rotating Earth, including a dynamic coupled ocean, see Mikolajewicz et al. (2018).
The strength of
d. Stationary waves and meridional heat transport
Last, we consider the impact of the substantially different stationary wave amplitudes on the atmospheric meridional heat transport in these different climates. Previous work has found that, for any given model, the maximum meridional heat transport is remarkably constant across a range of different climates, despite large differences in the partitioning of energy transport: e.g., stationary versus transient and/or moisture versus sensible heat transport (Donohoe and Battisti 2012; Donohoe et al. 2020). These studies looked at realistic climates ranging from the Last Glacial Maximum to CO2 forcing experiments. Our NM simulations, particularly NMBW, provide a stronger test of the invariance in maximum meridional heat transport. Table 1 summarizes the maximum meridional atmospheric heat transport (AHT) calculated using the difference between top of atmosphere radiative fluxes and surface heat fluxes (as in Donohoe et al. 2020). Despite markedly different climates and stationary wave strengths, the maximum AHT varies by less than 4%. Any changes in absorbed solar radiation from cloud or surface albedo changes must therefore be almost completely compensated by changes in outgoing longwave radiation. Interestingly, simply changing the version of the same model (WACCM4 and WACCM6) gives differences of ~12% for the CTL simulation, a much larger difference than removing mountains and rotating Earth backward in the same model. These experiments provide further evidence of the remarkable consistency of maximum AHT, for a given model, to a vast range of climates.
Maximum annual mean NH atmospheric meridional heat transport (AHT) in various simulations. Values in parentheses show differences from the CTL simulation. Standard error of the mean shown for the CTL is based on assuming independence of each year in the 20-yr simulation. Conclusions are very similar for the SH.


5. Summary
The impact of global orography on the climatological wintertime circulation is examined with the WACCM6 at 1° horizontal resolution, revisiting the seminal work of the 1980s and 1990s with a state-of-the-art climate model. We also quantify the relative contributions from resolved and parameterized orography. Through comparisons of simulations with and without orography, we find the following:
Orography substantially slows
Orography contributes at least as much as land–sea contrast and SSTs to the NH stationary waves. The observed features of the stationary waves in the lower troposphere—including the Aleutian low and Siberian high—are weak or absent in the absence of orography (Fig. 5, left column). In the free troposphere, the stationary waves attributable to land–sea contrast and orography are of approximately equal amplitude. In the stratosphere, stationary waves are predominantly driven by orography. These results are broadly consistent with previous work, but suggest a more important role for orography than the most recent estimates (e.g., Ting et al. 2001; Held et al. 2002; Chang 2009). The increased importance of orography in our study relative to previous work may be due to the role of parameterized orographic drag in forcing stationary waves (Fig. 6, bottom row), or due to differences in tropical SSTs, which we find to have a strong influence on the NM stationary waves.
Stationary waves depend on the positioning of orography relative to the land–sea boundaries. By repeating our CTL and no-mountain (NM) simulations with Earth rotating backward, we find that the orographic impact on both
Amalgamating these results, along with the sensitivity of the NM stationary waves to tropical SSTs discussed in the previous section, leads to the conclusion that producing a precise estimate of the role of orography in the forcing of the observed stationary waves is a formidable, and potentially impossible, task. The substantial impact of orography on the zonal mean flow means that stationary wave models, in which zonal mean flow is prescribed, cannot provide a full, or potentially even a meaningful, answer to this particular problem. By ignoring the impact of orography on the zonal flow
To conclude, orography plays a substantial role in producing the observed stationary waves, more so than most recent work suggests, perhaps due to the inclusion of relatively sophisticated parameterizations of sub-grid-scale orographic drag in the model used in this study. Orography also has considerable impact on the zonal mean zonal wind, explaining almost all of the difference in strength between the SH and NH zonal mean zonal wind throughout the free troposphere and stratosphere. The impact of orography on the stratosphere is dominated by the parameterized effects of unresolved orographic drag, although this manifests through a catalytic positive feedback mechanism that involves resolved processes. This study also highlights the important role of parameterized orographic drag in both stationary waves and troposphere–stratosphere coupling; insights gained by this study may be relevant for better understanding tropospheric–stratospheric connections, such as those associated with the quasi-biennial oscillation (QBO) or sudden stratospheric warmings.
Acknowledgments
RHW acknowledges support from the Natural Sciences and Engineering Research Council of Canada (NSERC) (RGPIN-2020-05783), the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Grant Agreement 797961, from the National Science Foundation Grant AGS-1665247, and from the Tamaki Foundation. DSB was supported by the Tamaki Foundation. Simulations were performed at the Department of Atmospheric Sciences at the University of Washington, and on the NCAR’s Cheyenne computer under Grant UWAS0064. The code to produce the figures shown in this paper can be found on RHW’s github account (rhwhite): https://github.com/rhwhite/SupportingInformation/blob/main/JASStationaryWaves2021.ipynb. Climatological output data from the simulations can be found here: https://doi.org/10.5683/SP2/NTEUHF; for daily data please contact the authors.
REFERENCES
Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen-Palm relation and the mean zonal acceleration. J. Atmos. Sci., 33, 2031–2048, https://doi.org/10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.
Baldwin, J., A. Atwood, G. Vecchi, and D. Battisti, 2021: Outsize influence of Central American orography on global climate. AGU Adv., in press.
Beljaars, A. C. M., A. R. Brown, and N. Wood, 2004: A new parametrization of turbulent orographic form drag. Quart. J. Roy. Meteor. Soc., 130, 1327–1347, https://doi.org/10.1256/qj.03.73.
Boer, G. J., N. A. McFarlane, R. Laprise, J. D. Henderson, and J.-P. Blanchet, 1984: The Canadian Climate Centre spectral atmospheric general circulation model. Atmos.–Ocean, 22, 397–429, https://doi.org/10.1080/07055900.1984.9649208.
Bolin, B., 1950: On the influence of the Earth’s orography on the general character of the westerlies. Tellus, 2, 184–195, https://doi.org/10.3402/tellusa.v2i3.8547.
Brayshaw, D. J., B. Hoskins, and M. Blackburn, 2008: The storm-track response to idealized SST perturbations in an aquaplanet GCM. J. Atmos. Sci., 65, 2842–2860, https://doi.org/10.1175/2008JAS2657.1.
Brayshaw, D. J., B. Hoskins, and M. Blackburn, 2009: The basic ingredients of the North Atlantic storm track. Part I: Land–sea contrast and orography. J. Atmos. Sci., 66, 2539–2558, https://doi.org/10.1175/2009JAS3078.1.
Broccoli, A. J., and S. Manabe, 1992: The effects of orography on midlatitude Northern Hemisphere dry climates. J. Climate, 5, 1181–1201, https://doi.org/10.1175/1520-0442(1992)005<1181:TEOOOM>2.0.CO;2.
Butler, A. H., J. P. Sjoberg, D. J. Seidel, and K. H. Rosenlof, 2017: A sudden stratospheric warming compendium. Earth Syst. Sci. Data, 9, 63–76, https://doi.org/10.5194/essd-9-63-2017.
Chang, E. K. M., 2009: Diabatic and orographic forcing of northern winter stationary waves and storm tracks. J. Climate, 22, 670–688, https://doi.org/10.1175/2008JCLI2403.1.
Charney, J. G., and A. Eliassen, 1949: A numerical method for predicting the perturbations of the middle latitude westerlies. Tellus, 1 (2), 38–54, https://doi.org/10.3402/tellusa.v1i2.8500.
Charney, J. G., and P. G. Drazin, 1961: Propagation of planetary-scale disturbances from the lower into the upper atmosphere. J. Geophys. Res., 66, 83–109, https://doi.org/10.1029/JZ066i001p00083.
Chen, S.-C., and K. E. Trenberth, 1988a: Forced planetary waves in the Northern Hemisphere winter: Wave-coupled orographic and thermal forcings. J. Atmos. Sci., 45, 682–704, https://doi.org/10.1175/1520-0469(1988)045<0682:FPWITN>2.0.CO;2.
Chen, S.-C., and K. E. Trenberth, 1988b: Orographically forced planetary waves in the Northern Hemisphere winter: Steady state model with wave-coupled lower boundary formulation. J. Atmos. Sci., 45, 657–681, https://doi.org/10.1175/1520-0469(1988)045<0657:OFPWIT>2.0.CO;2.
Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828.
Donohoe, A., and D. S. Battisti, 2012: What determines meridional heat transport in climate models? J. Climate, 25, 3832–3850, https://doi.org/10.1175/JCLI-D-11-00257.1.
Donohoe, A., K. C. Armour, G. H. Roe, D. S. Battisti, and L. Hahn, 2020: The partitioning of meridional heat transport from the Last Glacial Maximum to CO2 quadrupling in coupled climate models. J. Climate, 33, 4141–4165, https://doi.org/10.1175/JCLI-D-19-0797.1.
Edmon, H. J., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen-Palm cross sections for the troposphere. J. Atmos. Sci., 37, 2600–2616, https://doi.org/10.1175/1520-0469(1980)037<2600:EPCSFT>2.0.CO;2.
Eliassen, A., and E. Palm, 1961: On the transfer of energy in stationary mountain waves. Geofys. Publ., 22, 1–23.
Emile-Geay, J., 2003: Warren revisited: Atmospheric freshwater fluxes and “Why is no deep water formed in the North Pacific.” J. Geophys. Res., 108, 3178, https://doi.org/10.1029/2001JC001058.
Garcia, R. R., D. R. Marsh, D. E. Kinnison, B. A. Boville, and F. Sassi, 2007: Simulation of secular trends in the middle atmosphere, 1950–2003. J. Geophys. Res., 112, D09301, https://doi.org/10.1029/2006JD007485.
Garcia, R. R., A. K. Smith, D. E. Kinnison, Á. Cámara, and D. J. Murphy, 2017: Modification of the gravity wave parameterization in the Whole Atmosphere Community Climate Model: Motivation and results. J. Atmos. Sci., 74, 275–291, https://doi.org/10.1175/JAS-D-16-0104.1.
Garfinkel, C. I., I. White, E. P. Gerber, M. Jucker, and M. Erez, 2020: The building blocks of Northern Hemisphere wintertime stationary waves. J. Climate, 33, 5611–5633, https://doi.org/10.1175/JCLI-D-19-0181.1.
Gettelman, A., and Coauthors, 2019: The Whole Atmosphere Community Climate Model version 6 (WACCM6). J. Geophys. Res. Atmos., 124, 12 380–12 403, https://doi.org/10.1029/2019JD030943.
Hartmann, D. L., C. R. Mechoso, and K. Yamazaki, 1984: Observations of wave–mean flow interaction in the Southern Hemisphere. J. Atmos. Sci., 41, 351–362, https://doi.org/10.1175/1520-0469(1984)041<0351:OOWMFI>2.0.CO;2.
Held, I. M., 1983: Stationary and quasi-stationary eddies in the extratropical troposphere: Theory. Large-Scale Dynamical Processes in the Atmosphere, B. J. Hoskins and R. P. Pearce, Eds., Academic Press, 127–167.
Held, I. M., and M. Ting, 1990: Orographic versus thermal forcing of stationary waves: The importance of the mean low-level wind. J. Atmos. Sci., 47, 495–500, https://doi.org/10.1175/1520-0469(1990)047<0495:OVTFOS>2.0.CO;2.
Held, I. M., M. Ting, and H. Wang, 2002: Northern winter stationary waves: Theory and modeling. J. Climate, 15, 2125–2144, https://doi.org/10.1175/1520-0442(2002)015<2125:NWSWTA>2.0.CO;2.
Hoskins, B. J., and D. J. Karoly, 1981: The steady linear response of a spherical atmosphere to thermal and orographic forcing. J. Atmos. Sci., 38, 1179–1196, https://doi.org/10.1175/1520-0469(1981)038<1179:TSLROA>2.0.CO;2.
Hurrell, J. W., J. J. Hack, D. Shea, J. M. Caron, and J. Rosinski, 2008: A new sea surface temperature and sea ice boundary dataset for the Community Atmosphere Model. J. Climate, 21, 5145–5153, https://doi.org/10.1175/2008JCLI2292.1.
Inatsu, M., H. Mukougawa, and S.-P. Xie, 2002: Stationary eddy response to surface boundary forcing: Idealized GCM experiments. J. Atmos. Sci., 59, 1898–1915, https://doi.org/10.1175/1520-0469(2002)059<1898:SERTSB>2.0.CO;2.
Kitoh, A., 2002: Effects of large-scale mountains on surface climate—A coupled ocean-atmosphere general circulation model study. J. Meteor. Soc. Japan, 80, 1165–1181, https://doi.org/10.2151/jmsj.80.1165.
Li, Q., H.-F. Graf, and X. Cui, 2011: The role of stationary and transient planetary waves in the maintenance of stratospheric polar vortex regimes in Northern Hemisphere winter. Adv. Atmos. Sci., 28, 187–194, https://doi.org/10.1007/s00376-010-9163-7.
Lindvall, J., G. Svensson, and R. Caballero, 2017: The impact of changes in parameterizations of surface drag and vertical diffusion on the large-scale circulation in the Community Atmosphere Model (CAM5). Climate Dyn., 48, 3741–3758, https://doi.org/10.1007/s00382-016-3299-9.
Madden, R. A., and P. Speth, 1995: Estimates of atmospheric angular momentum, friction, and mountain torques during 1987–1988. J. Atmos. Sci., 52, 3681–3694, https://doi.org/10.1175/1520-0469(1995)052<3681:EOAAMF>2.0.CO;2.
Manabe, S., and T. B. Terpstra, 1974: The effects of mountains on the general circulation of the atmosphere as identified by numerical experiments. J. Atmos. Sci., 31, 3–42, https://doi.org/10.1175/1520-0469(1974)031<0003:TEOMOT>2.0.CO;2.
Marsh, D. R., M. J. Mills, D. E. Kinnison, J.-F. Lamarque, N. Calvo, and L. M. Polvani, 2013: Climate change from 1850 to 2005 simulated in CESM1(WACCM). J. Climate, 26, 7372–7391, https://doi.org/10.1175/JCLI-D-12-00558.1.
McFarlane, N. A., 1987: The effect of orographically excited gravity wave drag on the general circulation of the lower stratosphere and troposphere. J. Atmos. Sci., 44, 1775–1800, https://doi.org/10.1175/1520-0469(1987)044<1775:TEOOEG>2.0.CO;2.
Mikolajewicz, U., and Coauthors, 2018: The climate of a retrograde rotating Earth. Earth Syst. Dyn., 9, 1191–1215, https://doi.org/10.5194/esd-9-1191-2018.
Miller, M. J., T. N. Palmer, and R. Swinbank, 1989: Parametrization and influence of subgridscale orography in general circulation and numerical weather prediction models. Meteor. Atmos. Phys., 40, 84–109, https://doi.org/10.1007/BF01027469.
Nakamura, N., J. Falk, and S. W. Lubis, 2020: Why are stratospheric sudden warmings sudden (and intermittent)? J. Atmos. Sci., 77, 943–964, https://doi.org/10.1175/JAS-D-19-0249.1.
Nigam, S., I. M. Held, and S. W. Lyons, 1986: Linear simulation of the stationary eddies in a general circulation model. Part I: The no-mountain model. J. Atmos. Sci., 43, 2944–2961, https://doi.org/10.1175/1520-0469(1986)043<2944:LSOTSE>2.0.CO;2.
Nigam, S., I. M. Held, and S. W. Lyons, 1988: Linear simulation of the stationary eddies in a GCM. Part II: The “mountain” model. J. Atmos. Sci., 45, 1433–1452, https://doi.org/10.1175/1520-0469(1988)045<1433:LSOTSE>2.0.CO;2.
Okajima, H., and S.-P. Xie, 2007: Orographic effects on the northwestern Pacific monsoon: Role of air-sea interaction. Geophys. Res. Lett., 34, L21708, https://doi.org/10.1029/2007GL032206.
Palmer, T. N., G. J. Shutts, and R. Swinbank, 1986: Alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parametrization. Quart. J. Roy. Meteor. Soc., 112, 1001–1039, https://doi.org/10.1002/qj.49711247406.
Pithan, F., T. G. Shepherd, G. Zappa, and I. Sandu, 2016: Climate model biases in jet streams, blocking and storm tracks resulting from missing orographic drag. Geophys. Res. Lett., 43, 7231–7240, https://doi.org/10.1002/2016GL069551.
Polichtchouk, I., T. G. Shepherd, and N. J. Byrne, 2018: Impact of parametrized nonorographic gravity wave drag on stratosphere-troposphere coupling in the Northern and Southern Hemispheres. Geophys. Res. Lett., 45, 8612–8618, https://doi.org/10.1029/2018GL078981.
Randel, W. J., and L. R. Lait, 1991: Dynamics of the 4-day wave in the Southern Hemisphere polar stratosphere. J. Atmos. Sci., 48, 2496–2508, https://doi.org/10.1175/1520-0469(1991)048<2496:DOTDWI>2.0.CO;2.
Richter, J. H., F. Sassi, and R. R. Garcia, 2010: Toward a physically based gravity wave source parameterization in a general circulation model. J. Atmos. Sci., 67, 136–156, https://doi.org/10.1175/2009JAS3112.1.
Ringler, T. D., and K. H. Cook, 1995: Orographically induced stationary waves: Dependence on latitude. J. Atmos. Sci., 52, 2548–2560, https://doi.org/10.1175/1520-0469(1995)052<2548:OISWDO>2.0.CO;2.
Sandu, I., P. Bechtold, A. Beljaars, A. Bozzo, F. Pithan, T. G. Shepherd, and A. Zadra, 2016: Impacts of parameterized orographic drag on the Northern Hemisphere winter circulation. J. Adv. Model. Earth Syst., 8, 196–211, https://doi.org/10.1002/2015MS000564.
Sandu, I., and Coauthors, 2019: Impacts of orography on large-scale atmospheric circulation. npj Climate Atmos. Sci., 2, 10, https://doi.org/10.1038/s41612-019-0065-9.
Scinocca, J. F., and N. A. McFarlane, 2000: The parametrization of drag induced by stratified flow over anisotropic orography. Quart. J. Roy. Meteor. Soc., 126, 2353–2393, https://doi.org/10.1002/qj.49712656802.
Seager, R., D. S. Battisti, J. Yin, N. Gordon, N. Naik, A. C. Clement, and M. A. Cane, 2002: Is the Gulf Stream responsible for Europe’s mild winters? Quart. J. Roy. Meteor. Soc., 128, 2563–2586, https://doi.org/10.1256/qj.01.128.
Shaw, T. A., M. Sigmond, T. G. Shepherd, and J. F. Scinocca, 2009: Sensitivity of simulated climate to conservation of momentum in gravity wave drag parameterization. J. Climate, 22, 2726–2742, https://doi.org/10.1175/2009JCLI2688.1.
Sigmond, M., and J. F. Scinocca, 2010: The influence of the basic state on the Northern Hemisphere circulation response to climate change. J. Climate, 23, 1434–1446, https://doi.org/10.1175/2009JCLI3167.1.
Sinha, B., A. T. Blaker, J. J.-M. Hirschi, S. Bonham, M. Brand, S. Josey, R. S. Smith, and J. Marotzke, 2012: Mountain ranges favour vigorous Atlantic meridional overturning. Geophys. Res. Lett., 39, L02705, https://doi.org/10.1029/2011GL050485.
Smagorinsky, J., 1953: The dynamical influence of large-scale heat sources and sinks on the quasi-stationary mean motions of the atmosphere. Quart. J. Roy. Meteor. Soc., 79, 342–366, https://doi.org/10.1002/QJ.49707934103.
Taguchi, M., and D. L. Hartmann, 2006: Increased occurrence of stratospheric sudden warmings during El Niño as simulated by WACCM. J. Climate, 19, 324–332, https://doi.org/10.1175/JCLI3655.1.
Ting, M., H. Wang, and L. Yu, 2001: Nonlinear stationary wave maintenance and seasonal cycle in the GFDL R30 GCM. J. Atmos. Sci., 58, 2331–2354, https://doi.org/10.1175/1520-0469(2001)058<2331:NSWMAS>2.0.CO;2.
Tucker, G. B., 1960: The atmospheric budget of angular momentum. Tellus, 12, 134–144, https://doi.org/10.3402/tellusa.v12i2.9385.
Valdes, P. J., and B. J. Hoskins, 1989: Linear stationary wave simulations of the time-mean climatological flow. J. Atmos. Sci., 46, 2509–2527, https://doi.org/10.1175/1520-0469(1989)046<2509:LSWSOT>2.0.CO;2.
Valdes, P. J., and B. J. Hoskins, 1991: Nonlinear orographically forced planetary waves. J. Atmos. Sci., 48, 2089–2106, https://doi.org/10.1175/1520-0469(1991)048<2089:NOFPW>2.0.CO;2.
van Niekerk, A., J. F. Scinocca, and T. G. Shepherd, 2017: The modulation of stationary waves, and their response to climate change, by parameterized orographic drag. J. Atmos. Sci., 74, 2557–2574, https://doi.org/10.1175/JAS-D-17-0085.1.
van Niekerk, A., I. Sandu, and S. B. Vosper, 2018: The circulation response to resolved versus parametrized orographic drag over complex mountain terrains. J. Adv. Model. Earth Syst., 10, 2527–2547, https://doi.org/10.1029/2018MS001417.
Wahr, J. M., and A. H. Oort, 1984: Friction- and mountain-torque estimates from global atmospheric data. J. Atmos. Sci., 41, 190–204, https://doi.org/10.1175/1520-0469(1984)041<0190:FAMTEF>2.0.CO;2.
Wallace, J. M., S. Tibaldi, and A. J. Simmons, 1983: Reduction of systematic forecast errors in the ECMWF model through the introduction of an envelope orography. Quart. J. Roy. Meteor. Soc., 109, 683–717, https://doi.org/10.1002/qj.49710946202.
Warren, B. A., 1983: Why is no deep water formed in the North Pacific? J. Mar. Res., 41, 327–347, https://doi.org/10.1357/002224083788520207.
Waugh, D. W., and L. M. Polvani, 2010: Stratospheric polar vortices. The Stratosphere: Dynamics, Transport, and Chemistry, L. M. Polvani, A. H. Sobel, and D. W. Waugh, Eds., Amer. Geophys. Union, 43–57.
Webster, S., A. R. Brown, D. R. Cameron, and C. P. Jones, 2003: Improvements to the representation of orography in the Met Office Unified Model. Quart. J. Roy. Meteor. Soc., 129, 1989–2010, https://doi.org/10.1256/qj.02.133.
White, R. H., D. S. Battisti, and G. H. Roe, 2017: Mongolian mountains matter most: Impacts of the latitude and height of Asian orography on Pacific wintertime atmospheric circulation. J. Climate, 30, 4065–4082, https://doi.org/10.1175/JCLI-D-16-0401.1.
White, R. H., D. S. Battisti, and A. Sheshadri, 2018: Orography and the boreal winter stratosphere: The importance of the Mongolian mountains. Geophys. Res. Lett., 45, 2088–2096, https://doi.org/10.1002/2018GL077098.
Wills, R. C., and T. Schneider, 2015: Stationary eddies and the zonal asymmetry of net precipitation and ocean freshwater forcing. J. Climate, 28, 5115–5133, https://doi.org/10.1175/JCLI-D-14-00573.1.
Wills, R. C., and T. Schneider, 2018: Mechanisms setting the strength of orographic Rossby waves across a wide range of climates in a moist idealized GCM. J. Climate, 31, 7679–7700, https://doi.org/10.1175/JCLI-D-17-0700.1.
Yin, J. H., and D. S. Battisti, 2001: The importance of tropical sea surface temperature patterns in simulations of Last Glacial Maximum climate. J. Climate, 14, 565–581, https://doi.org/10.1175/1520-0442(2001)014<0565:TIOTSS>2.0.CO;2.