The global atmospheric energy transport analysed by a wavelengthbased scale separation
 ^{1}Department of Mathematics and Statistics, Arctic University of Norway, Tromsø, Norway
 ^{2}Department of Physics and Technology, Arctic University of Norway, Tromsø, Norway
 ^{3}Norwegian Meteorological Institute, Tromsø, Norway
 ^{4}Department of Earth Sciences, Uppsala University, Uppsala, Sweden
 ^{5}Centre of Natural Hazards and Disaster Science (CNDS), Uppsala, Sweden
 ^{6}Department of Meteorology, Stockholm University, Stockholm, Sweden
 ^{7}Bolin Centre for Climate Research, Stockholm, Sweden
 ^{1}Department of Mathematics and Statistics, Arctic University of Norway, Tromsø, Norway
 ^{2}Department of Physics and Technology, Arctic University of Norway, Tromsø, Norway
 ^{3}Norwegian Meteorological Institute, Tromsø, Norway
 ^{4}Department of Earth Sciences, Uppsala University, Uppsala, Sweden
 ^{5}Centre of Natural Hazards and Disaster Science (CNDS), Uppsala, Sweden
 ^{6}Department of Meteorology, Stockholm University, Stockholm, Sweden
 ^{7}Bolin Centre for Climate Research, Stockholm, Sweden
Correspondence: Patrick Johannes Stoll (patrick.stoll@uit.no)
Hide author detailsCorrespondence: Patrick Johannes Stoll (patrick.stoll@uit.no)
The atmosphere transports energy polewards by circulation cells and eddies. To the present day, there has been a knowledge gap regarding the preferred spatial scales and physical mechanisms of eddy energy transport. To fill the gap, we separate the meridional atmospheric energy transport in the ERA5 reanalysis by spatial scales and into quasistationary and transient flow patterns and latent and dry components.
Baroclinic instability is the major instability mechanism in the transient synoptic scales and is responsible for forming cyclones, anticyclones, and smallscale Rossby waves. At the planetary scales, circulation patterns are often induced by other mechanisms such as flow interaction with orography and land–sea heating contrasts. However, a separation between circulation patterns at the synoptic and planetary scales has yet to be established. We find that both baroclinically induced and transient energy transport is predominantly associated with eddies at wavelengths between 2000 and 8000 km. The maxima in both types of transport occur at wavelengths around 5000 km, in good agreement with linear baroclinic theory. Since these results are independent of latitude, we adapt the scale separation of the energy transport to be based on the wavelength instead of the previously used wavenumber. We define the synoptic transport by the wavelength band between 2000 and 8000 km.
We analyse the annual and seasonal mean in the energy transport components and their interannual variability. The scaleseparated transport components are fairly similar in both hemispheres. Transport by synoptic waves is the largest contributor to extratropical energy and moisture transport, mainly of a transient character, and is influenced little by seasonality. In contrast, transport by planetary waves depends highly on the season and has two distinct characteristics. (1) In the extratropical winter, planetary waves are important due to a large transport of dry energy. This planetary transport features the largest interannual variability of all components and is mainly quasistationary in the Northern Hemisphere but transient in its southern counterpart. (2) In the subtropical summer, quasistationary planetary waves are the most important transport component, mainly due to moisture transport, presumably associated with monsoons. In contrast to transport by planetary and synoptic waves, only a negligible amount of energy is transported by mesoscale eddies (< 2000 km).
 Article
(2687 KB) 
Supplement
(2374 KB)  BibTeX
 EndNote
Atmospheric motions reduce the thermal contrast created by differential solar heating between high and low latitudes (Hadley, 1735). Hence the atmospheric circulation transports large amounts of energy polewards (Oort and Peixóto, 1983) and has thereby a fundamental role in controlling the local weather and climate (e.g. Holton and Hakim, 2013; Vallis, 2017). The energy is primarily transported in the form of warm air (dry energy: comprising mainly enthalpy and potential energy) and water vapour (latent energy) which releases energy when condensating before precipitation (Peixoto and Oort, 1992).
The characteristics of the energy transport differ among latitudinal zones (Trenberth and Stepaniak, 2003a). In the tropics and subtropics, where the Coriolis effect is small, energy is predominantly transported by a zonally symmetric meridional overturning circulation, known as the Hadley cell (Hadley, 1735), and monsoon systems organised by quasistationary cells (Fig. 1a). In the extratropics, eddies transport the majority of the energy further poleward. Atmospheric eddies exist on a large range of scales, from planetary Rossby waves (Rossby, 1939) and transient synoptic Rossby waves, which vertically interact with synoptic cyclones (Bjerknes, 1919), to mesoscale disturbances, such as polar lows (Businger and Reed, 1989). Conventionally, the eddy transport is separated into a quasistationary component and a transient component (Fig. 1a), with the former representing monthlymean eddies and the latter fastervarying deviations from this mean (Oort and Peixóto, 1983). Transient eddies dominate the extratropical energy transport in both hemispheres (Fig. 1a), whereas quasistationary waves transport a considerable amount of energy in the higher latitudes of the Northern Hemisphere (NH) and the subtropics of both hemispheres. However, they are less relevant in the extratropical Southern Hemisphere (SH).
Recent studies demonstrate the usefulness of separating the energy transport by spatial scales into planetary and synoptic components. For instance, eddies at these two scales impact the Arctic differently (Baggett and Lee, 2015; Graversen and Burtu, 2016). The scale separation of the meridional energy transport by a zonal Fourier decomposition became popular in recent years. It was applied to study the mechanisms and impacts of energy transport in the Arctic region (Papritz and DunnSigouin, 2020; Graversen et al., 2021; Rydsaa et al., 2021; Hofsteenge et al., 2022) and the NH midlatitudes (Lembo et al., 2019). Further, numerous studies investigated the atmospheric dynamics in terms of Rossby waves at different spatial scales (e.g. Wirth et al., 2018; Röthlisberger et al., 2019).
So far, the energy transport across all latitudes has only been separated by a fixed wavenumber as shown in Fig. 1b and presented by Graversen and Burtu (2016). The separation between planetary and synoptic scales has often been established between waves 3 and 4 (Baggett and Lee, 2015; Heiskanen et al., 2020). For investigation of the scale of the transport at a specific latitude or a small zonal band, such a separation by wavenumber is appropriate. However, the wavelength associated with a given wavenumber is latitude dependent, which needs to be accounted for when defining the wavenumber separating spatial scales (Heiskanen et al., 2020). For instance, wave 4 corresponds to a wavelength of 8200 km at 35^{∘} but only to 2600 km at 75^{∘}. These two wavelengths can be interpreted to represent different spatial scales. Such a partitioning by wavenumber thus has two caveats: (i) towards the poles, all eddy transport converges to the planetary scale (Fig. 1b). (ii) In the subtropics, wavenumbers 1–3 appear to miss parts of planetary transport, as can be inferred from quasistationary eddies in the subtropical SH (Fig. 1a) being considerably larger than the planetary transport captured by wavenumbers 1–3 (Fig. 1b).
This study proposes a revised partitioning by spatial scales based on wavelength to circumvent the abovementioned problems associated with separation by wavenumber. With this, we partition the transport into planetary, synoptic, and mesoscale components to understand their role in transporting energy poleward. In addition, we apply the conventional decomposition of the transport into stationary and transient parts (Fig. 1a), as well as dry and latent components (Peixoto and Oort, 1992).
The new scaleseparation method is employed to analyse meridional energy transport by the global atmospheric circulation. Previous studies found that transport by quasistationary waves is important in the NH but almost negligible in the SH (e.g. Peixoto and Oort, 1992). Such quasistationary transport is often associated with the planetary scale, which appears to imply that planetary transport is irrelevant in the SH (e.g. Trenberth and Stepaniak, 2003a). However, this is not the case since planetary waves are strongly represented in the SH as shown in Fig. 1b. In contrast, transport by transient eddies is often associated with baroclinic eddies at the synoptic scale (e.g. Trenberth and Stepaniak, 2003a). However, transport at other scales could be transient as well. In this study, we underscore that the separation of transport into a quasistationary and transient contribution differs from a separation into planetary and synoptic scales. However, a correlation exists between the two, especially for the NH.
The following summarises the main research questions posed in this study:

At what scales does the atmosphere transport energy, and how does the scaleseparation method compare to the conventional separation method?

What characterises the atmospheric energy transport in different meridional bands?
These questions are investigated in Sects. 3 and 4. The main results are then summarised and discussed in Sect. 5. However, first, we present the utilised data and methods.
2.1 ERA5 reanalysis
The atmospheric energy transport for the period 1979–2021 is analysed with ERA5 (Hersbach et al., 2020), the fifth atmospheric reanalysis from the European Centre for MediumRange Weather Forecasts (ECMWF). ERA5 provides hourly fields at a spectral truncation of T639, equivalent to a grid spacing of 30 km, and with 137 vertical hybridsigma levels. For this study, fields of air temperature, specific humidity, geopotential height and horizontal wind components at all hybrid levels, and surface pressure and surface topography are used on a 0.25^{∘} × 0.25^{∘} horizontal grid spacing. Since reanalysis data are prone to include massflux inconsistencies (Trenberth, 1991), a barotropic massflux correction is applied to the wind field before calculating the energy transport (Graversen, 2006).
2.2 Zonal mean versus zonal integral
The zonal integral of the ERA5 energy transport (Fig. S1a) confirms that found in previous studies using datasets with lower spatial resolution (Peixoto and Oort, 1992; Trenberth and Caron, 2001; Graversen and Burtu, 2016). For instance, the zonally integrated poleward transport peaks at 4.8 × 10^{15} W in the NH and 5.6 × 10^{15} W in the SH at a latitude of 41^{∘} in both hemispheres. By computing the zonal integral of the energy transport, which depends on the length of the longitude circle, the transport becomes small at high latitudes since the longitudes converge (Fig. S1a). However, the local transport, expressed by the zonal mean, is considerable also in the polar regions (Fig. S1b). Hence, to compare the local importance of the atmospheric energy transport across all latitudes, we take a zonalmean perspective which provides the transport through an atmospheric column with 1 m width. The latitude of maximum zonalmean transport, namely the location at which the largest local meridional transport occurs, is at a latitude of 45^{∘} in both hemispheres (Fig. S1b), different from the maximum zonally integrated transport peaking at 41^{∘}.
For the calculation of convergence of energy transport, a 2^{∘} runningmean filter is used along the meridional dimension before computing the meridional derivatives.
2.3 Decomposition of the energy transport
The atmospheric energy transport and its components are characterised by a large daytoday variability (Swanson and Pierrehumbert, 1997; Messori and Czaja, 2013). Here we focus on the annualmean and seasonalmean energy transport, firstly to investigate the timemean behaviour of the atmospheric circulation and secondly to compare the newly introduced separation of the eddy transport based on spatial scales – described below – with the conventional separation into quasistationary and transient transport introduced by Oort and Peixóto (1983) and commonly used in the literature (e.g. Trenberth and Stepaniak, 2003a, b; Kaspi and Schneider, 2013). The conventional separation obtains quasistationary eddies from monthlymean fields; hence a comparison is only possible from a timemean perspective. Quasistationary eddies could be derived differently, for example, by application of a temporal lowpass filter. This would enable the comparison of the separation method by scale and by quasistationary and transient components without taking time averages. However, it would prevent a direct comparison with the literature and require a new ad hoc analysis, which is beyond the scope of this study.
We define the vertically integrated, zonalmean, timemean, meridional transport of energy as
where v is the meridional wind, E is the atmospheric energy, p is pressure, p_{s} is surface pressure, g is the gravitational constant, $\stackrel{\mathrm{\u203e}}{\cdot}$ denotes a monthly time average, and [⋅] denotes a zonal average. We decompose the total atmospheric energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, in three ways that can be applied in succession:
Firstly, the total atmospheric energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, is partitioned into latent energy transport, $\stackrel{\mathrm{\u0303}}{vQ}$, and dry energy transport, $\stackrel{\mathrm{\u0303}}{vD}$, following Oort and Peixóto (1983). This is achieved by separating the total atmospheric energy, E, in Eq. (1) into its latent component, Q, and dry component, D, which comprises the enthalpy, c_{p}T; potential energy, gz; and kinetic energy, $\frac{{\mathit{v}}^{\mathrm{2}}}{\mathrm{2}}$, as
with q being specific humidity, L being the latent heat release by condensation, c_{p} being specific heat capacity at constant pressure, T being temperature, and v being the horizontal wind vector.
Secondly, the total energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, expressing a timemean transport, is separated into transport by transient eddies, ${\stackrel{\mathrm{\u0303}}{vE}}^{\mathrm{tr}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}\left[\stackrel{\mathrm{\u203e}}{{v}^{\prime}{E}^{\prime}}\right]\frac{\mathrm{d}p}{g}$; quasistationary eddies, ${\stackrel{\mathrm{\u0303}}{vE}}^{\text{qs}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}\left[{\stackrel{\mathrm{\u203e}}{v}}^{*}{\stackrel{\mathrm{\u203e}}{E}}^{*}\right]\frac{\mathrm{d}p}{g}$; and the mean meridional circulation, ${\stackrel{\mathrm{\u0303}}{vE}}^{\mathrm{md}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}\left[\stackrel{\mathrm{\u203e}}{v}\right]\left[\stackrel{\mathrm{\u203e}}{E}\right]\frac{\mathrm{d}p}{g}$, by partitioning the zonalmean, timemean transport, $\left[\stackrel{\mathrm{\u203e}}{vE}\right]$, in Eq. (1) according to Oort and Peixóto (1983) and Peixoto and Oort (1992) by
with ${\cdot}^{\prime}$ denoting anomalies from the time mean and ${\cdot}^{*}$ denoting zonal anomalies from the zonal mean. The first equality separates the timemean transport into the transport by the timemean fields, $\left[\stackrel{\mathrm{\u203e}}{v}\stackrel{\mathrm{\u203e}}{E}\right]$, and by fastervarying fluctuations, $\left[\stackrel{\mathrm{\u203e}}{{v}^{\prime}{E}^{\prime}}\right]$. The second equality partitions the timemean term into a zonally symmetric transport component, $\left[\stackrel{\mathrm{\u203e}}{v}\right]\left[\stackrel{\mathrm{\u203e}}{E}\right]$, and a contribution by zonal waves, $\left[{\stackrel{\mathrm{\u203e}}{v}}^{*}{\stackrel{\mathrm{\u203e}}{E}}^{*}\right]$. The duration of the time mean, $\stackrel{\mathrm{\u203e}}{\cdot}$, for which eddies are considered quasistationary, is arbitrary. However, conventionally the mean is taken over 1month periods (e.g. Oort and Peixóto, 1983). Hereby, the quasistationary eddy transport, $\left[{\stackrel{\mathrm{\u203e}}{v}}^{*}{\stackrel{\mathrm{\u203e}}{E}}^{*}\right]$, provides the transport given by the monthlymean eddies, and the transient transport, $\left[\stackrel{\mathrm{\u203e}}{{v}^{\prime}{E}^{\prime}}\right]$, is the faster varying component. The annualmean energy transport partitioned in this conventional manner is depicted in Fig. 1a. The transport of latent energy, $\stackrel{\mathrm{\u0303}}{vQ}$, and dry energy, $\stackrel{\mathrm{\u0303}}{vD}$, is partitioned similarly.
Thirdly, the total energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, is separated into transport by eddies at the planetary scale, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$; the synoptic scale, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}$; and the mesoscale, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{ms}}$, as well as the meridional circulation, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}$. For the partitioning into spatial scales, a Fourier decomposition is applied along each latitude, ϕ; at each vertical model level; and at every time step for E and v^{1} following Graversen and Burtu (2016). The Fourier series of E yields
where x is the zonal coordinate, d=2πacos (ϕ); a=6371 km (Earth's radius); and the Fourier coefficients, ${a}_{n}^{E}$ and ${b}_{n}^{E}$, are given as
where n is zonal wavenumber. The Fourier coefficients for the meridional wind, v, ${a}_{n}^{v}$ and ${b}_{n}^{v}$, are derived in a similar manner.
Based on this Fourier decomposition, the timemean, zonalmean energy transport can be written as
Note that the cross terms ${a}_{n}^{v}$ and ${a}_{m}^{E}$ with n≠m, and similarly ${b}_{n}^{v}$ and ${b}_{m}^{E}$, vanish when computing a zonal average since the Fourier coefficients feature an orthogonal basis. In the following, we show how this decomposition relates to the decomposition by Oort and Peixóto (1983) expressed in Eq. (3).
The first term in Eq. (6) can be approximated by the first term on the righthand side (RHS) of Eq. (3):
since $\stackrel{\mathrm{\u203e}}{\left[{v}^{\prime}\right]\left[{E}^{\prime}\right]}$, the time mean of the transient part of the meridional circulation, is small. Hence, the timemean transport by the meridional overturning circulation can also be written as ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}\frac{\mathrm{1}}{\mathrm{4}}\stackrel{\mathrm{\u203e}}{{a}_{\mathrm{0}}^{v}{a}_{\mathrm{0}}^{E}}\frac{\mathrm{d}p}{g}$.
The second term in Eq. (6) is a timemean transport by eddies, $\left[\stackrel{\mathrm{\u203e}}{{v}^{*}{E}^{*}}\right]$, and can be separated into stationary and transient eddy transport:
The stationary eddy part is the second term on the RHS of Eq. (3) and can also be expressed by
The transient eddy part, $\left[\stackrel{\mathrm{\u203e}}{{{v}^{\prime}}^{*}{{E}^{\prime}}^{*}}\right]$, is the difference between the total transient transport, $\left[\stackrel{\mathrm{\u203e}}{{v}^{\prime}{E}^{\prime}}\right]$ (term 3 on RHS of Eq. 3), and the transient part of the meridional circulation $\stackrel{\mathrm{\u203e}}{\left[{v}^{\prime}\right]\left[{E}^{\prime}\right]}$, which is again small in the time mean, hence
Having shown the relation of the terms in the Fourier decomposition with the transient and quasistationary terms in the conventional decomposition, we now introduce the scale separation of the eddy transport. The instantaneous zonalmean transport by wave n is given by
To separate the energy transport between two scales based on a predefined wavelength, λ_{sep}, the latitudedependent (ϕ) wavenumber of separation, n_{sep}, is computed:
The wavenumber of separation between planetary and synoptic eddies, ${n}_{\mathrm{p}/\mathrm{s}}$, is computed at each latitude from the wavelength ${\mathit{\lambda}}_{\mathrm{p}/\mathrm{s}}=\mathrm{8000}$ km (lower black solid line in Fig. 2). Synoptic and mesoscale eddies are separated at a wavelength of ${\mathit{\lambda}}_{\mathrm{s}/\mathrm{m}}=\mathrm{2000}$ km. In Sect. 3, we argue for the usage of these two wavelengths for scaleseparating the energy transport. The wavenumbers of separation, as defined here, are real numbers. To ensure a continuous separation (Fig. S2b), instead of abrupt transitions (Fig. S2a), the transport by planetary, synoptic, and mesoscale eddies is defined by
with ⌊⋅⌋ denoting rounding to the integer part and ⌈⋅⌉ rounding to the least integer. This provides the latitudedependent separation depicted in Fig. 2 between the black curves at 2000 and 8000 km.
The partitioning of Eqs. (13)–(15) is applied to calculate the timemean transport by planetary eddies, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}[\stackrel{\mathrm{\u203e}}{{v}^{*}{E}^{*}}{]}_{\mathrm{p}}\frac{\mathrm{d}p}{g}$; synoptic eddies, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}[\stackrel{\mathrm{\u203e}}{{v}^{*}{E}^{*}}{]}_{\mathrm{sy}}\frac{\mathrm{d}p}{g}$; and mesoscale eddies, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{ms}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}[\stackrel{\mathrm{\u203e}}{{v}^{*}{E}^{*}}{]}_{\mathrm{ms}}\frac{\mathrm{d}p}{g}$. It is in a similar manner applied to calculate the transport by quasistationary eddies at different scales, for instance the transport by quasistationary planetary eddies, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\text{qs}}={\int}_{\mathrm{0}}^{{p}_{\mathrm{s}}}[{\stackrel{\mathrm{\u203e}}{v}}^{*}{\stackrel{\mathrm{\u203e}}{E}}^{*}{]}_{\mathrm{p}}\frac{\mathrm{d}p}{g}$. The transport by transient eddies at the different scales is simply obtained from the difference between the timemean transport and the quasistationary transport. For instance, the transport by transient planetary eddies is calculated as ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\mathrm{tr}}={\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}{\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\text{qs}}$.
An underlying interpretation is that at a longitude circle with for example 20 000 km extent (at a latitude of about 60^{∘}), wave 1 includes the transport at a scale of and larger than the longitude cycle, wave 2 is between 10 000 and 20 000 km, and so forth. Hereby, wave 3 with a wavelength of 6700 km provides the transport between its wavelength and 10 000 km and is therefore partitioned between the synoptic and planetary scales. Since the wavenumber of separation between the planetary and synoptic scale is ${n}_{\mathrm{p}/\mathrm{s}}=\mathrm{2.5}$ at a latitude of 60^{∘}, this wave is equally partitioned between the two scales (second term in Eq. 13 and first term in Eq. 14).
This continuous partitioning of the transport by the first wave smaller than 8000 km leads to the synoptic scale including less transport than the application of a strict separation (Fig. S2). Hence, waves at a wavelength around and larger than 8000 km are associated with the planetary scale, whereas the synoptic scale constitutes waves strictly smaller than 8000 km. Spectra depicting the separation into the scale components at different latitudes are provided in Fig. S3, and a detailed illustration of the separation is provided in Sect. S3 in the Supplement. It should be noted that Wolf and Wirth (2017) also applied a continuous scale separation for the latitudedependent detection of Rossby wave packages. The same procedure of scale separation is applied to the latent and dry energy transport.
Generally, the scale separation based on the Fourier decomposition is nonlocal, implying that the whole latitude circle influences the obtained eddies (Heiskanen et al., 2020). Therefore the Fourier decomposition is useful if the transport across the circle is governed by similar eddy scales, which we can observe from meteorological weather maps along zonal bands. Arguably, the zonal and meridional scales of atmospheric eddies match: from the investigation of meteorological weather maps, we know (i) that synopticscale cyclones are circular to the first order and (ii) that the meridional extent, i.e. the amplitude, of Rossby waves, appears to roughly match the distance between a trough and a ridge in the zonal direction. Hence, as many other studies (e.g. Graversen and Burtu, 2016; Lembo et al., 2019), we interpret the zonal wavenumber of an eddy, which is associated with a zonal wavelength at a given latitude, as its spatial scale.
Orlanski (1975) defines the mesoscale to be smaller than 2000 km, which is commonly accepted and hence used in this study as the lower threshold for the synoptic scale. However, a widely agreed separation between the planetary and synoptic scales is not established. Below, we argue for placing the separation between the planetary and synoptic scales at a wavelength of 8000 km, based on the argument that baroclinic instability is the major instability and hence the physical mechanism of wave generation at the synoptic scale. The exact value is to some degree arbitrary. As previously presented, this partitioning is done continuously such that transport at and around the wavelength of 8000 km is associated with the planetary scale, whereas synoptic waves are strictly smaller than this threshold.
3.1 Theoretical baroclinic argument
Baroclinic instability is recognised as the single dominant instability of the synoptic scale (Markowski and Richardson, 2011, p. 4). Thus, the synoptic scale should include most of the energy transports associated with eddies developing by baroclinic instability. Synoptic eddies are typically identified as cyclones and anticyclones in sealevel pressure, as well as shorter Rossby waves (Frederiksen, 1978) often referred to as transient synoptic Rossby waves (e.g. Röthlisberger et al., 2019; Ali et al., 2021). The cyclones and anticyclones interact vertically and baroclinically with the synoptic Rossby waves, generally recognised by uppertropospheric oscillations of the jet stream. The theoretical scale (wavelength) of baroclinic eddies is given by 3.9 times the Rossby deformation radius, L_{d}, and hence estimated to be 4000 km by Vallis (2017, p. 354) and 4800 km by Stoll et al. (2021) using more exact values obtained from extratropical cyclones. Note that a lowpressure (highpressure) system spans half a wavelength and has accordingly a typical diameter of around 2000 km.
A wavelength band between 2000–8000 km appears appropriate to capture the majority of the transport associated with baroclinically induced synoptic eddies for three reasons. (i) Cyclones and anticyclones feature some variability in their size but with a typical diameter between 1000 and smaller than 4000 km. (ii) Short, synoptic Rossby waves are considered to occur at zonal wavelengths within this band, with meridional amplitudes of around half a wavelength, as further discussed in Sect. 3.4. (iii) The nonlocal Fourier decomposition of the energy transport in situations of localised synoptic cyclones captures considerable transport at neighbouring waves to the cyclone (Heiskanen et al., 2020, Fig. 3), which may be considered a weakness of the Fourier decomposition. However, the band of wavenumbers around the localised eddy appears appropriate to capture the eddy.
3.2 Statistical baroclinic argument
Figure 3 depicts the anomalous daily energy transport by different waves at a lag time of 4 d after enhanced meridional temperature gradients (90th percentile). Different lag times (>0) provide similar results (Fig. S4), and the presented one of 4 d is in good agreement with Fig. 1a of Baggett and Lee (2015). The separation between waves featuring enhanced and decreased transport as a response to the increased thermal contrast is at around 8000 km for all latitudes in the extratropics. Waves whose energy transport is most closely associated with the increased temperature gradient have a wavelength between 3000 and 8000 km with a maximum transport anomaly for wavelengths around 5000–6000 km. Since baroclinic instability is induced by a horizontal temperature gradient, this provides strong statistical evidence that baroclinic instability is an important development mechanism for waves at wavelengths between 2000 and 8000 km, which is here defined as the synoptic scale. In the following, we refer to the anomalous energy transport 4 d after situations of enhanced meridional temperature gradients as a “baroclinically induced transport anomaly”. The strength of the baroclinically induced transport anomaly by the synoptic waves (Fig. 3 at a wavelength between 2000 and 8000 km) is around a third of the timemean energy transport of these waves (Fig. 2). The mean meridional temperature contrast will also induce baroclinic waves transporting energy, which is not captured by this methodology.
The latitudinal independence of the wavelength of the baroclinically induced transport anomaly provides a major argument for using a wavelength threshold to separate between the synoptic and planetary scales rather than a wavenumber threshold. It is surprising that the scale of maximum baroclinically induced transport anomaly is independent of the latitude since the Rossby deformation radius, ${L}_{\mathrm{d}}=\frac{NH}{f}$, estimating the size of baroclinic eddies (Vallis, 2017), depends inversely on the Coriolis parameter, f, increasing with latitude and depends linearly on the depth of the troposphere, H, and the tropospheric static stability, N, which rather decrease with latitude (Stoll et al., 2021). Hence, baroclinic eddies would be expected to be smaller at higher latitudes. A hypothesis for the latitude independence of the baroclinic eddies is as follows. Most extratropical cyclones originate from the midlatitudes, where the meridional temperature contrast is the largest. The size of a cyclone is set during the genesis stage when the fastestgrowing mode is prevailing. Many cyclones propagate to higher latitudes along the diagonal axis of the storm tracks (Shaw et al., 2016) and may preserve their original size.
We note that processes other than baroclinic instability may contribute to the formation of the heredefined synoptic eddies, such as heating contrasts. A future study will further investigate the drivers of eddy transport at different scales.
3.3 Transient and quasistationary eddies
In the following, we show that the partitioning at a wavelength of 8000 km approximately captures our intuitive understanding that synoptic cyclones and synoptic Rossby wave packets are transient by nature, whereas most quasistationary eddies are at planetary scales since they may be constrained by largescale orography and semistationary thermal forcing, such as heating contrasts between ocean and land (Vallis, 2017).
Indeed, the spectral decomposition of the annualmean energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, at different latitudes reveals that most eddies smaller than 8000 km are of a transient nature, whereas most of the quasistationary transport is at scales larger than 8000 km (Fig. 2). This is in good agreement with Dell'Aquila et al. (2005), who found in the average of the extratropical NH, that wave 3, here defined as planetary scale, features a typical period of around a month, whereas smaller eddies, here associated with the synoptic scale, are characterised by weekly or daily periods.
For the transient energy transport, the wavenumber of maximum transport is 6 or 7 in the subtropics and decreases towards the poles such that the corresponding wavelength of maximum transient transport is around 5000 to 6000 km for all latitudes. Also in the moisture transport, the transient component reaches its maximum around 5000 km (Fig. S5), which is in good agreement with Lee et al. (2019). This wavelength corresponds well to the preferred scale of baroclinic eddies (4000–5000 km: Vallis, 2017; Stoll et al., 2021) and of baroclinically induced transport anomalies presented in the previous section. The latitudinal independence of the wavelength of maximum transient transport supports separating the planetary and synoptic scale by using wavelength instead of wavenumber.
The preferred wavenumber and wavelength of the quasistationary transport, ${\stackrel{\mathrm{\u0303}}{vE}}^{\text{qs}}$, are latitude dependent (Fig. 2a). The wavelength of maximum quasistationary transport is larger than 8000 km everywhere besides close to the North Pole where the Fourier decomposition is not applicable and around 40^{∘} S where the quasistationary transport is almost negligible. A short analysis of the preferred scales of the quasistationary transport is provided in Sect. S5.
The separation of synoptic and planetary waves at a wavelength of 8000 km captures our intuitive understanding that (i) most transient transport occurs at the synoptic scale, whereas (ii) most quasistationary transport is at the planetary scale. The inverse is valid for the former but not the latter: planetary waves at a wavelength larger than 8000 km can be both transient and quasistationary (Fig. 2). In the latitude band between 70 and 30^{∘} S and poleward of 70^{∘} N, zones mainly characterised by open ocean, the transient component within the planetary transport is larger than the quasistationary component (Fig. 4). At latitudes with more land, the quasistationary component is responsible for 50 %–80 % of the planetary transport, meaning that also at these latitudes a considerable proportion of the planetary transport is transient. Also the moisture transport, $\stackrel{\mathrm{\u0303}}{vQ}$, comprises a planetary component that features both transient and quasistationary contributions, similar to the total energy transport, $\stackrel{\mathrm{\u0303}}{vE}$ (Fig. S6).
In contrast to planetary transport, the inverse of (i) above, that synoptic transport is mainly (70 %–100 %) of a transient nature, is valid at almost all latitudes (Fig. 4). This coincides with the transient character of synoptic cyclones and Rossby waves of short wavelengths. The small quasistationary contribution (0 %–30 %) to the synoptic transport is attributed to preferred spatial locations for synoptic activity. For instance, the NH Atlantic sector features more cyclonic activity than other longitudes, resulting in increased quasistationary transport in the time mean. This can be inferred from Rydsaa et al. (2021), who show a large timemean synoptic transport in the Atlantic sector for strong latent transport events in winter at 70^{∘} N. However, in a zonalmean perspective, the quasistationary contribution to the synoptic transport is small (< 30 %) compared to its transient part. Hence, for simplicity's sake, the synoptic transport is not separated into a transient and quasistationary contribution in the remainder of this study. An exception for the transient character of synoptic transport is over Antarctica, which is characterised by katabatic flows advecting cold air towards lower latitudes, occurring at preferred locations of drainage, hence including a large stationary component.
We conclude that the conventional decomposition of the transport into stationary and transient parts (Oort and Peixóto, 1983) is to some degree related but not equivalent to a separation into planetary and synoptic waves.
3.4 Comparison to previous studies
Previous studies performing a wave decomposition qualitatively agree on the separation at a wavelength of 8000 km:
Firstly, Baggett and Lee (2015) perform a Fourier decomposition of energy characteristics for the entire NH and demonstrate different life cycle behaviour of long (planetary) and short (synoptic) waves for a separation between waves 3 and 4. Ali et al. (2021) identify recurrent synopticscale transient Rossby wave packets in the midlatitudes at wavenumbers between 4 and 15. Wirth et al. (2018) note that planetaryscale Rossby waves are typically characterised by zonal wavenumbers 1–3, whereas higher wavenumbers characterise synopticscale Rossby waves.
Secondly, their separation into planetary and synoptic eddies by wavenumber is similar to the hereapplied partitioning by wavelength at a latitude of 53 ^{∘}, which is close to the location of maximum eddy activity (Fig. 1). For a hemispherewide separation into planetary and synoptic eddies with a single wavenumber threshold as utilised by Baggett and Lee (2015) and Ali et al. (2021), it appears crucial to capture the partitioning around this latitude.
Thirdly, at 70^{∘} N, a wavelength of 8000 km is associated with a separation number n_{sy}=1.7, meaning that 70 % of wave 2 is associated with the planetary scale and 30 % is associated with the synoptic scale. Accordingly, the lead–lag regression of the Arctic temperature on the latent transport presented in Fig. 6 of Graversen and Burtu (2016) shows a clear difference between wave 1, leading to heating of the Arctic, and waves 3 and larger, having a baroclinic signal. Wave 2 appears to share characteristics of both wave groups, in agreement with the hereapplied partitioning of wave 2 at 70^{∘} N. Also Heiskanen et al. (2020) find a good attribution of the meridional energy transport by an idealised cyclone to the synoptic scale if the separation between the planetary and synoptic scales is performed at rather small wavenumbers, even though they do not test a separation between waves 2 and 3.
In conclusion, separating at wavelengths of 8000 km is consistent with previous interpretations of planetaryscale and synopticscale transport. However, as mentioned earlier, a sharp threshold likely does not exist. Therefore, the separation wavelength of 8000 km is compared to wavelengths of 10 000 and 6000 km in the Supplement (Fig. S7). The main results of this study are not affected by the exact choice of the separation wavelength, as shortly discussed in the Supplement.
In this section, we analyse the atmospheric energy transport from ERA5 by utilising the scaleseparation method. The applied latitudinal bands used in this study are provided in Table 1.
4.1 Annualmean transport
4.1.1 Overview
The meridional atmospheric energy transport features considerable similarities in the two hemispheres in most components and for most climate zones. Hence, to simplify the comparison of the two hemispheres, we display the poleward transport of both hemispheres at different latitudes on a common x axis (Fig. 5a–c).
The annualmean, zonalmean poleward energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, is seamless in the sense that it resembles half of a sine curve for both hemispheres (black lines in Fig. 5a) as noted by Trenberth and Stepaniak (2003b). However, in the subtropics and midlatitudes, the energy transport is ∼ 15 % larger in the SH than in the NH. The difference is approximately balanced by more oceanic transport in the NH (Trenberth and Caron, 2001).
The energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, leads to an annualmean divergence of around 40 W m^{−2} in the tropics and subtropics and convergence poleward of a latitude of 40^{∘} of up to 100 W m^{−2} in the polar regions (Fig. 5d) in good agreement with Trenberth and Stepaniak (2003a). The smaller divergence close to the Equator is due to oceanic currents transporting heat to the subtropics, creating a cold tongue of seasurface temperatures over the equatorial Pacific (Trenberth and Stepaniak, 2003b).
Despite the total energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, appearing seamless, different transport mechanisms are important in the different climate zones, as reflected by considerable variations across latitudes in the moisture and dry transports (Fig. 5b and c) and the scale components of the transport.
The total annualmean moisture transport, $\stackrel{\mathrm{\u0303}}{vQ}$, of both hemispheres features equatorward extremes at around 10^{∘} (Fig. 5b), poleward maxima at around 40^{∘}, and decaying tails towards the poles. This leads to moisture divergence in the nonequatorial tropics and subtropics, whereas moisture converges in the equatorial regions and extratropics (Fig. 5e). The moisture transport is generally stronger in the Southern Hemisphere than Northern Hemisphere due to more evaporation on the water surfaces of the SH. Further, some moisture is transported across the Equator from the SH to the NH, leading to the highest convergence of moisture at around 7^{∘} N, which is in the annual mean the approximate location of the intertropical convergence zone (ITCZ).
The dry transport, $\stackrel{\mathrm{\u0303}}{vD}$, is plateaulike between a latitude of 10–65^{∘} (Fig. 5c) and decays towards the Equator and poles. Hence, it is mainly responsible for divergence in the equatorial regions and convergence in the polar regions (Fig. 5f). However, the scale components have different roles in the different climate zones.
4.1.2 Synoptic transport
Synopticscale waves dominate the energy transport in the midlatitudes of both hemispheres (Fig. 5a). However, synoptic transport is around 40 % stronger in the SH than in the NH. This is in broad agreement with the finding of Peixoto and Oort (1992) that transient transport is more relevant in the SH, but it should be kept in mind that some of the transient transport occurs at the planetary scale (Fig. 4). Despite a broad spectrum of possible waves, most energy is transported by eddies at wavelengths between 2000 and 8000 km, demonstrating the dominant role of synoptic waves for the midlatitudinal energy transport (Fig. S7). This provides a different perspective from previous studies finding that transient eddies are responsible for the majority of midlatitudinal energy transport (e.g. Peixoto and Oort, 1992), since planetary waves may also be transient. The peak transport by synoptic waves is at a latitude of around 45^{∘}, the location of the climatological storm tracks (Hoskins and Hodges, 2002, 2005), and coincides with the zone of largest total energy transport. Hence, synoptic waves transport energy from the subtropics (< 40^{∘}) to the extratropics (Fig. 5d).
Extratropical synoptic waves transport approximately twothirds of their energy in dry form and onethird in latent form (Fig. 5b, c). The latent contribution is slightly higher in the midlatitudes close to the subtropics and lower in the polar regions, where the cold atmosphere cannot hold much water vapour. Despite transporting more dry than moist energy, synoptic eddies are responsible for most extratropical moisture transport. Especially in the polar regions, almost all moisture transport is by synoptic waves, whereas dry transport in the polar regions occurs at both synoptic and planetary scales.
4.1.3 Planetary transport
The planetary energy transport is similar in both hemispheres, different from quasistationary transport which is mainly relevant in the NH (Figs. 1a, 4, 5a). The latter is in agreement with Trenberth and Stepaniak (2003a), underscoring that quasistationary transport is a primary factor in the extratropical NH. The latter study associates this quasistationary transport with the planetary scale, albeit without proof – a hypothesis confirmed by this study (Fig. 4). A new finding that could partly be inferred from Fig. S3 of Lembo et al. (2019) is the nearsymmetric structure of the planetary energy transport in both hemispheres. This symmetry could not have been anticipated by considering only quasistationary transport since the planetary transport in the SH is mainly transient (Fig. 4), in agreement with Mo (1986).
In both hemispheres, the planetary transport is similar in the subtropics and lower midlatitudes. In the higher midlatitudes and polar region of the SH, it is approximately 20 % weaker than in the corresponding NH regions. Hence, eddies at similar spatial scales are transporting the energy in both hemispheres (see also Fig. 2), likely due to similar physical mechanisms in both hemispheres forming the energytransporting eddies.
Generally, two patterns of planetary transport are identified:

In the subtropics, most planetary transport is associated with quasistationary moisture transport, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{p}}^{\text{qs}}$. These waves presumably represent subtropical highpressure and monsoon systems, which can prevail for several weeks. The planetary moisture transport leads to humidity divergence in the tropics and convergences in the sub and extratropics (Fig. 5e).

In the extratropics, planetary eddies mainly transport dry energy and only little moisture. In the polar regions, this planetary energy transport is almost as important as synoptic transport. The peak in energy transport by planetary waves, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, is at the polar boundary at a latitude of around 60^{∘} in both hemispheres (Fig. 5a, c), hence further poleward than the peak in synoptic transport at 45^{∘}. Thus, the planetary dry transport, ${\stackrel{\mathrm{\u0303}}{vD}}_{\mathrm{p}}$, leads to energy divergence in the subtropics and midlatitudes and convergence in the polar regions (Fig. 5f). These planetary waves are mainly transient (Fig. 4: 70 %) in the midlatitudinal SH, whereas the quasistationary component dominates (60 %) in the corresponding NH region, which agrees with Peixoto and Oort (1992). However, in the high latitudes of the SH, a considerable fraction of the planetary transport is quasistationary.
4.1.4 Meridional overturning circulation
The meridional overturning circulation, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}$, and its role in transporting energy poleward have long been known (e.g. Hadley, 1735; Lorenz, 1967; Peixoto and Oort, 1992). The Hadley circulation in the tropics dominates the energy transport in that region. The total energy transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}$, by the Hadley circulation is small compared to the thermally direct dry component, ${\stackrel{\mathrm{\u0303}}{vD}}_{\mathrm{md}}$, and to the thermally indirect latent component, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{md}}$, since the latter largely compensates the former (Fig. 5a–c). Further, the annualmean transport by the Hadley circulation is small since transport into the tropical winter hemisphere originates from the equatorial summer hemisphere, leading to compensation of summer and winter cells (see Sect. 4.2).
In the midlatitudes, the meridional circulation, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}$, features the thermally indirect Ferrel cell with a peak around 53^{∘}. The Ferrel cell is almost twice as strong in the SH than NH. The transport by the Ferrel cell is mainly in the form of dry energy, which is to a small extent compensated by thermally direct moisture transport. In the NH, the eddydriven Ferrel cell (Vallis, 2017), spans the whole extratropics, including the Arctic. In the SH, a thermally direct polar cell is evident in the dry transport, ${\stackrel{\mathrm{\u0303}}{vD}}_{\mathrm{md}}$, which is primarily driven by katabatic flow from Antarctica as noted by Trenberth and Stepaniak (2003a). Different to previous studies by for example Peixoto and Oort (1992) a NH polar cell is not evident in the hereutilised ERA5 dataset, neither in the annual mean nor in the summer or winter season (Sect. 4.2). In the Arctic, energy is transported by synoptic and planetary eddies favouring the formation of a Ferrel cell, whereas zonal symmetric katabatic flows, as observed in the Antarctic, do not develop due to the lack of a large ice dome centred over the pole.
4.1.5 Mesoscale transport
In contrast to synoptic and planetary waves, mesoscale waves, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{ms}}$, at scales smaller than 2000 km, are only responsible for a negligible part of the energy and moisture transport at all latitudes (Figs. 2, S7), in accordance with Graversen and Burtu (2016), who show that wavenumbers 0–10 are responsible for the majority of the energy transport at all latitudes. Atmospheric models may have larger challenges in reproducing mesoscale eddies as compared to eddies at larger scales. However, ERA5 at a horizontal grid spacing equivalent to 30 km can reproduce mesoscale cyclones, such as polar lows (Stoll et al., 2021), and hence the negligible importance of mesoscale eddies for the total energy transport appears trustworthy. Due to its negligible role, we include the mesoscale into the synoptic transport for the remainder of this study such that the sum of all scale components yields the total transport.
4.2 Seasonal transport
Some transport patterns become more apparent when seasons are analysed separately. The NH summer and the SH winter are defined as June to August, and the NH winter and the SH summer are December to February. In spring and autumn the energy transport is mainly similar to the annualmean transport (Figs. 5, S8).
Due to seasonal variations in the thermal contrast between the Equator and poles, more energy is transported poleward in the winter than in the summer hemisphere (Fig. 6a, d) in good agreement with previous studies (e.g. Peixoto and Oort, 1992). The seasonality in the transport is larger in the NH than the SH, as also noted by Trenberth and Stepaniak (2003a). This is affected by a larger annual temperature cycle in the NH due to its large continents having a smaller heat capacity than the oceans in the SH. Still, the two hemispheres share many characteristics in the atmospheric energy transport regarding seasonality.
The location separating northward and southward total transport, the energy flux equator (Adam et al., 2016), is at around a latitude of 10^{∘} in the summer hemisphere (Fig. 6d). This is close to the ascending branch of the Hadley circulation (Vallis, 2017, p. 514). In the tropics, the meridional overturning circulation, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}$, is most important (Fig. 6, blue). It transports energy from the summer tropics to the winter subtropics, whereas moisture, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{md}}$, is transported in the opposite direction.
In the extratropics of both hemispheres, the synoptic transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}$, is little influenced by seasonality, in broad agreement with results from wavenumber separated transport of the NH midlatitudes in Lembo et al. (2019) and with transient transport in Trenberth and Stepaniak (2003a). In contrast, the planetary transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, has a strong seasonality, being highly relevant in winter but almost absent in summer. Hence, the summer is dominated by synoptic transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}$. Differently, in winter, both planetary and synoptic waves are highly relevant for the energy transport. In the NH winter, planetary waves, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, contribute to more transport than do synoptic waves, mainly by their quasistationary component, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\text{qs}}$. In the SH winter, planetary transport is also important but mostly in its transient form, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\mathrm{tr}}$, whereas its stationary part is small as noted previously (Oort and Peixóto, 1983).
In the subtropics of the summer hemisphere, quasistationary planetary waves, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\text{qs}}$, are among the largest contributors to poleward energy transport (Fig. 6d). These waves are transporting energy mainly in the form of moisture, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{p}}^{\text{qs}}$ (Fig. 6e), presumably reflecting the summer monsoon and longlasting subtropical highpressure systems.
4.3 Interannual variability
In this study, the interannual variability in the energy transport is computed by the standard deviation of the annualmean transport. The total energy transport, $\stackrel{\mathrm{\u0303}}{vE}$, is only varying by a few percent between years for all latitudes (Fig. 7a). However, the interannual variability in the individual components is up to 3 times larger than the variability in the total transport, especially in the extratropics (Fig. 7a). Hence, large transport in one component is typically compensated by smaller transport in another, as noted by Trenberth and Stepaniak (2003a) and Lembo et al. (2019).
In contrast to the total energy transport, the variability of the total moisture transport, $\stackrel{\mathrm{\u0303}}{vQ}$, is larger than the variability of its individual scale components (Fig. 7b). Hence, the moisture transport components are not compensating each other in the same manner as the total energy transport components. Instead, the compensation of the components of the total energy transport is in the form of dry energy (Fig. 7c).
We hypothesise that the different covariability of the scale contributions for the total energy and moisture transport is due to their different underlying mechanisms. Preliminary results point towards the annualmean energy transport being induced by the meridional temperature gradient in the manner of a diffusion process with a globally almost constant diffusion coefficient. Hence, large transport in one component reduces the temperature gradient, leading to less transport in another. Differently, moisture is a tracer of the atmospheric circulation and, therefore, not described by a diffusion process such that the components do not compensate similarly to the total energy transport.
The tropics feature large interannual variability in the moisture transport, $\stackrel{\mathrm{\u0303}}{vQ}$, and approximately equally large variability in the dry transport, $\stackrel{\mathrm{\u0303}}{vD}$ (Fig. 7b, c). The large total variability in the moist and dry energy transports is almost entirely due to variability in the meridional overturning. This is not surprising since the meridional overturning is responsible for most of the moist and dry energy transport in the tropics (Figs. 5b, c, 6b, c, e, f). However, the variability in total transport, $\stackrel{\mathrm{\u0303}}{vE}$, is only a fifth of the moisture transport variability (Fig. 7a, b). Hence, a Hadley circulation causing a largerthanusual thermally direct dry energy transport in a single year also causes a largerthanusual thermally indirect latent energy transport in that year. This explains that the total energy transported by the Hadley circulation is more or less constant between years.
In the extratropics, the planetary transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, is the component exhibiting the largest interannual variability (Fig. 7a). This is remarkable since the annualmean planetary transport is mainly smaller (Fig. 5) than the synoptic transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}$. Consequently, the planetary transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, varies from year to year by approximately 10 % in the midlatitudes and 15 %–20 % in the polar regions. In contrast, the synoptic transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}$, varies only by around 5 %. The large planetary variability is mainly attributed to its quasistationary component, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\text{qs}}$, whereas its transient component, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\mathrm{tr}}$, is much less variable (not shown). This demonstrates the “quasiness” of the quasistationary transport, being the transport component that varies most from year to year. This presumably reflects the role of (inter)annual modes of climate variability in favouring different quasistationary circulations and planetary transports in different years.
The extratropical moisture transport, $\stackrel{\mathrm{\u0303}}{vQ}$, is varying approximately equally for planetary, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{p}}$, and synoptic waves, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{sy}}$ (Fig. 7b). Since synoptic waves are responsible for most of the moisture transport in the extratropics (Fig. 5b), the variability fraction is considerably higher for the planetary than the synoptic moisture transport.
In this study, we analyse the global atmospheric circulation by separating the meridional energy transport of the years 1979–2021 in the ERA5 reanalysis by spatial scales, by moist and dry components, and by quasistationary and transient parts. For separating the energy transport by scales for all latitudes, we apply a new approach by using the wavelength instead of the wavenumber utilised previously (e.g Graversen and Burtu, 2016).
We demonstrate that separating transport by synoptic and planetary eddies at a wavelength of 8000 km reflects the physically grounded distinction between baroclinically induced and other eddies. Moreover, we show that the wavelength of 8000 km is reasonable for separating at all latitudes. This separation wavelength approximately agrees with the conventional separation between transient and quasistationary eddies, as most eddy transport smaller than 8000 km is transient, whereas most quasistationary transport occurs at the planetary scale larger than 8000 km. Despite the latter, a considerable amount of the planetary energy transport is transient, especially in the extratropical SH.
It should be noted that the scale separation is implemented continuously so that the largest eddy with a wavelength smaller than 8000 km is partitioned between the planetary and synoptic scales. This implies that the planetary transport includes waves of around 8000 km in size, whereas synoptic transport is by waves strictly smaller than 8000 km.
Eddies with a wavelength of almost 8000 km may appear as large for being part of the synoptic scale. However, most baroclinically induced and most transient energy transport organises at wavelengths around 5000 km at all latitudes (Figs. 2, 3), agreeing well with the predicted length by drybaroclinic theory (Vallis, 2017). Despite the maximum at around 5000 km, the baroclinically induced and transient energy transport occurs in a wavelength band between approximately 2000 and 8000 km; hence separating at around 5000 km would be misleading.
Synoptic eddies are perceived as low and highpressure systems near the surface and smallscale Rossby waves that are often strongest around the tropopause. In the baroclinic development, nearsurface pressure systems vertically interact with the upperlevel Rossby waves. One synoptic wave includes both a low and a highpressure system. Hence, a 2000–8000 km range implies that synoptic cyclones and anticyclones feature a typical diameter of between 1000 and smaller than 4000 km and that the typical distance between the cores of two independent (anti)cyclones is between 2000 and 8000 km. Both appear reasonable from comparison with weather maps. Even in situations of cyclone clusters, arguably situations where the distance between cyclones is unusually low, the spacing between cyclonic centres is around 2000 km (e.g. Weijenborg and Spengler, 2020, Fig. 2d–f), the lower end of the heredefined synoptic scale. Note that the cyclones of that cluster share a common front and hence are not independent. The distance of the cyclones is in the northeastward direction, whereas the closest cyclones in the zonal direction are at larger distances of around 5000 km. Hence, our methodology relying on a zonal Fourier decomposition appears accurate for such situations.
A single, sharp spatial separation between the synoptic and planetary scales may not exist in the real world. However, a spatial decomposition is helpful for better understanding the atmospheric circulation and its impact on regional climate (Baggett and Lee, 2015; Graversen and Burtu, 2016; Röthlisberger et al., 2019). Here, we provide further arguments for the usefulness of separating the energy transport by spatial scale.
Different from the classical separation into quasistationary and transient energy transport, the spatial scales of the transport are quite similar in both hemispheres, pointing towards similarities in the contribution of active physical mechanisms. The most pronounced difference between the two separation methods is that planetary transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, is broadly similar in both hemispheres, whereas quasistationary transport, ${\stackrel{\mathrm{\u0303}}{vE}}^{\text{qs}}$, is mainly important in the NH.
In the annual mean, most energy and moisture in the extratropics are transported by synoptic eddies, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{sy}}$. It is remarkable that in the large range of atmospheric eddies, those at scales in the band between 2000 and 8000 km are responsible for most of the meridional energy transport for the whole extratropics. This points towards the importance of baroclinic instability for inducing the extratropical energy transport as long hypothesised (Holton and Hakim, 2013). However, other mechanisms may contribute to the formation of eddies at the synoptic scale. Hence, a future study is planned to investigate the causes and effects of energy transport at different scales.
The synoptic energy transport reveals to be influenced little by the season. In contrast, extratropical planetary transport, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}$, is of major importance only in the winter season mainly by transporting dry energy, ${\stackrel{\mathrm{\u0303}}{vD}}_{\mathrm{p}}$. In winter around the Arctic boundary, quasistationary planetary waves, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{p}}^{\text{qs}}$, dominate the energy transport. Such quasistationary planetary waves around the polar boundaries of both hemispheres feature the transport component with the largest interannual variability globally.
Other known characteristics of the atmospheric circulation are reproduced in this study, such as the dominance of the Hadley circulation in the tropics for transporting energy, ${\stackrel{\mathrm{\u0303}}{vE}}_{\mathrm{md}}$, to the winter hemisphere and transporting moisture, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{md}}$, in the opposite direction. In the subtropical summer, quasistationary planetary waves transporting moisture, ${\stackrel{\mathrm{\u0303}}{vQ}}_{\mathrm{p}}^{\text{qs}}$, are the largest contributor to poleward energy transport, which is associated with monsoon systems.
In this study, the atmospheric transport is analysed on an annualmean and seasonalmean basis, whereas on shorter timescales, the transport can be highly sporadic and displays large deviations from climatology (Swanson and Pierrehumbert, 1997; Messori and Czaja, 2013; Lembo et al., 2019). Some planetary transport events in the extratropics are for example equatorward. Hence, when the seasonalmean poleward planetary transport is close to zero, this can mean either that the planetary transport is generally small or that poleward and equatorward transport events are balancing each other. A followup study will investigate the intraseasonal distribution of the different transport components. Another option for investigating the spatiotemporal scale of eddies is the usage of Hayashi spectra as in Dell'Aquila et al. (2005) that performs a Fourier decomposition in both space and time.
The computed decomposition of the energy transport based on ERA5 and the code for the analysis are willingly provided on request.
The supplement related to this article is available online at: https://doi.org/10.5194/wcd412023supplement.
RGG calculated the Fourierdecomposed energy transport. PJS performed further data processing and visualisation. PJS, RGG, and GM contributed to interpreting the results and writing the manuscript.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Thanks go to ECMWF for providing access to data from the ERA5 reanalysis. The data were partly processed with the supercomputer Fram and stored at NIRD, both provided by the Norwegian Research Infrastructure Services (NRIS) Sigma2 AS (project nos. NN9348K and NS9063K, respectively). We also thank two anonymous reviewers and the editors for their constructive comments, which considerably improved the paper.
This research has been supported by the Research Council of Norway (NFR) under the projects “The role of the atmospheric energy transport in recent Arctic climate change” (no. 280727) and “Stability of the Arctic climate” (no. 314570).
This paper was edited by Camille Li and reviewed by two anonymous referees.
Adam, O., Bischoff, T., and Schneider, T.: Seasonal and Interannual Variations of the Energy Flux Equator and ITCZ. Part I: Zonally Averaged ITCZ Position, J. Climate, 29, 3219–3230, https://doi.org/10.1175/JCLID150512.1, 2016. a
Ali, S. M., Martius, O., and Röthlisberger, M.: Recurrent Rossby wave packets modulate the persistence of dry and wet spells across the globe, Geophys. Res. Lett., 48, e2020GL091452, https://doi.org/10.1029/2020GL091452, 2021. a, b, c
Baggett, C. and Lee, S.: Arctic warming induced by tropically forced tapping of available potential energy and the role of the planetaryscale waves, J. Atmos. Sci., 72, 1562–1568, 2015. a, b, c, d, e, f
Bjerknes, J.: On the structure of moving cyclones, Mon. Weather Rev., 47, 95–99, https://doi.org/10.1175/15200493(1919)47<95:OTSOMC>2.0.CO;2, 1919. a
Businger, S. and Reed, R. J.: Cyclogenesis in cold air masses, Weather Forecast., 4, 133–156, 1989. a
Dell'Aquila, A., Lucarini, V., Ruti, P. M., and Calmanti, S.: Hayashi spectra of the northern hemisphere midlatitude atmospheric variability in the NCEP–NCAR and ECMWF reanalyses, Clim. Dynam., 25, 639–652, 2005. a, b
Frederiksen, J.: Instability of planetary waves and zonal flows in twolayer models on a sphere, Q. J. Roy. Meteor. Soc., 104, 841–872, 1978. a
Graversen, R. G.: Do changes in the midlatitude circulation have any impact on the Arctic surface air temperature trend?, J. Climate, 19, 5422–5438, 2006. a
Graversen, R. G. and Burtu, M.: Arctic amplification enhanced by latent energy transport of atmospheric planetary waves, Q. J. Roy. Meteor. Soc., 142, 2046–2054, 2016. a, b, c, d, e, f, g, h, i, j
Graversen, R. G., Stoll, P. J., and Rydsaa, J. H.: Hvorfor er Arktis det området i verden med raskest oppvarming?, Naturen, 145, 160–167, 2021. a
Hadley, G.: VI. Concerning the cause of the general tradewinds, Philos. T. Roy. Soc. Lond., 39, 58–62, 1735. a, b, c
Heiskanen, T., Graversen, R. G., Rydsaa, J. H., and Isachsen, P. E.: Comparing wavelet and Fourier perspectives on the decomposition of meridional energy transport into synoptic and planetary components, Q. J. Roy. Meteor. Soc., 146, 2717–2730, 2020. a, b, c, d, e
Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., MuñozSabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, 2020. a
Hofsteenge, M. G., Graversen, R. G., Rydsaa, J. H., and Rey, Z.: The impact of atmospheric Rossby waves and cyclones on the Arctic sea ice variability, Clim. Dynam., 59, 579–594, https://doi.org/10.1007/s0038202206145z, 2022. a
Holton, J. and Hakim, G.: An Introduction to Dynamic Meteorology, vol. 5, Academic Press, New York, USA, ISBN 9780123848666, 2013. a, b
Hoskins, B. J. and Hodges, K. I.: New perspectives on the Northern Hemisphere winter storm tracks, J. Atmos. Sci., 59, 1041–1061, 2002. a
Hoskins, B. J. and Hodges, K. I.: A new perspective on Southern Hemisphere storm tracks, J. Climate, 18, 4108–4129, 2005. a
Kaspi, Y. and Schneider, T.: The role of stationary eddies in shaping midlatitude storm tracks, J. Atmos. Sci., 70, 2596–2613, 2013. a
Lee, S., Woods, C., and Caballero, R.: Relation between Arctic moisture flux and tropical temperature biases in CMIP5 simulations and its fingerprint in RCP8.5 projections, Geophys. Rese. Lett., 46, 1088–1096, 2019. a
Lembo, V., Messori, G., Graversen, R., and Lucarini, V.: Spectral Decomposition and Extremes of Atmospheric Meridional Energy Transport in the Northern Hemisphere Midlatitudes, Geophys. Res. Lett., 46, 7602–7613, https://doi.org/10.1029/2019GL082105, 2019. a, b, c, d, e, f, g
Lorenz, E.: The nature and theory of the general circulation of the atmosphere, World Meteorological Organization, No. 218, TP. 115, https://library.wmo.int/index.php?lvl=notice_display&id=5571#.YtEtPYTMJro (last access: 1 December 2022), 1967. a
Markowski, P. and Richardson, Y.: Mesoscale meteorology in midlatitudes, vol. 2, John Wiley & Sons, Chichester, UK, https://doi.org/10.1002/9780470682104, 2011. a
Messori, G. and Czaja, A.: On the sporadic nature of meridional heat transport by transient eddies, Q. J. Roy. Meteor. Soc., 139, 999–1008, 2013. a, b
Mo, K. C.: QuasiStationary States in the Southern Hemisphere, Mon. Weather Rev., 114, 808–823, https://doi.org/10.1175/15200493(1986)114<0808:QSSITS>2.0.CO;2, 1986. a
Oort, A. H. and Peixóto, J. P.: Global angular momentum and energy balance requirements from observations, Adv. Geophys., 25, 355–490, 1983. a, b, c, d, e, f, g, h, i, j
Orlanski, I.: A rational subdivision of scales for atmospheric processes, B. Am. Meteorol. Soc., 56, 527–530, 1975. a
Papritz, L. and DunnSigouin, E.: What configuration of the atmospheric circulation drives extreme net and total moisture transport into the Arctic, Geophys. Res. Lett., 47, e2020GL089769, https://doi.org/10.1029/2020GL089769, 2020. a
Peixoto, J. P. and Oort, A. H.: Physics of climate, American Institute of Physics, ISBN 9780883187128, 1992. a, b, c, d, e, f, g, h, i, j, k
Rossby, C.G.: Relation between variations in the intensity of the zonal circulation of the atmosphere and the displacements of the semipermanent centers of action, J. Mar. Res., 2, 38–55, 1939. a
Röthlisberger, M., Frossard, L., Bosart, L. F., Keyser, D., and Martius, O.: Recurrent synopticscale Rossby wave patterns and their effect on the persistence of cold and hot spells, J. Climate, 32, 3207–3226, 2019. a, b, c
Rydsaa, J. H., Graversen, R., Heiskanen, T. I. H., and Stoll, P.: Changes in atmospheric latent energy transport into the Arctic: Planetary versus synoptic scales, Q. J. Roy. Meteor. Soc., 147, 2281–2292, https://doi.org/10.1002/qj.4022, 2021. a, b
Shaw, T., Baldwin, M., Barnes, E. A., Caballero, R., Garfinkel, C., Hwang, Y.T., Li, C., O'gorman, P., Rivière, G., Simpson, I. R., and Voigt, A.: Storm track processes and the opposing influences of climate change, Nat. Geosci., 9, 656–664, 2016. a
Stoll, P. J., Spengler, T., Terpstra, A., and Graversen, R. G.: Polar lows – moistbaroclinic cyclones developing in four different vertical wind shear environments, Weather Clim. Dynam., 2, 19–36, https://doi.org/10.5194/wcd2192021, 2021. a, b, c, d
Swanson, K. L. and Pierrehumbert, R. T.: Lowertropospheric heat transport in the Pacific storm track, J. Atmos. Sci., 54, 1533–1543, 1997. a, b
Trenberth, K. E.: Climate diagnostics from global analyses: Conservation of mass in ECMWF analyses, J. Climate, 4, 707–722, 1991. a
Trenberth, K. E. and Caron, J. M.: Estimates of meridional atmosphere and ocean heat transports, J. Climate, 14, 3433–3443, 2001. a, b
Trenberth, K. E. and Stepaniak, D. P.: Covariability of components of poleward atmospheric energy transports on seasonal and interannual timescales, J. Climate, 16, 3691–3705, 2003a. a, b, c, d, e, f, g, h, i, j
Trenberth, K. E. and Stepaniak, D. P.: Seamless poleward atmospheric energy transports and implications for the Hadley circulation, J. Climate, 16, 3706–3722, 2003b. a, b, c
Vallis, G. K.: Atmospheric and oceanic fluid dynamics, Cambridge University Press, https://doi.org/10.1017/9781107588417, 2017. a, b, c, d, e, f, g, h
Weijenborg, C. and Spengler, T.: Diabatic heating as a pathway for cyclone clustering encompassing the extreme storm Dagmar, Geophys. Res. Lett., 47, e2019GL085777, https://doi.org/10.1029/2019GL085777, 2020. a
Wirth, V., Riemer, M., Chang, E. K., and Martius, O.: Rossby wave packets on the midlatitude waveguide – A review, Mon. Weather Rev., 146, 1965–2001, 2018. a, b
Wolf, G. and Wirth, V.: Diagnosing the Horizontal Propagation of Rossby Wave Packets along the Midlatitude Waveguide, Mon. Weather Rev., 145, 3247–3264, https://doi.org/10.1175/MWRD160355.1, 2017. a
To be precise, the Fourier decomposition is applied on the mass flux, $v\phantom{\rule{0.125em}{0ex}}\mathrm{d}p/g$, since the thickness, dz, times density, ρ, providing the mass thickness of the model layers, $\mathrm{d}p/g=\mathit{\rho}\mathrm{d}z$, is nonconstant in hybridsigma coordinates.