## 1. Introduction

The global atmospheric circulation transports momentum, energy, moisture, and chemical constituents both horizontally and vertically, and consequently impacts both global and regional climate. The conventional Eulerian-mean circulation, averaged at constant pressure or height, displays a three-cell pattern in each hemisphere, with a tropical Hadley cell, a midlatitude Ferrel cell, and a polar cell (e.g., Peixoto and Oort 1992). By contrast, to the extent that an isentropic surface can be regarded as a material surface in the absence of diabatic heating or diffusion, the mean circulation averaged at constant potential temperature can approximate the Lagrangian motion of an air parcel. In the isentropic coordinates, the mean circulation exhibits one equator-to-pole cell, with a poleward circulation in the upper troposphere and an equatorward return flow near the surface (e.g., Townsend and Johnson 1985; Tung 1986; Iwasaki 1989; Juckes et al. 1994; Held and Schneider 1999; Tanaka et al. 2004; Czaja and Marshall 2006; Pauluis et al. 2008, 2010). The upward branch of the isentropic circulation corresponds to diabatic heating in the tropics and the descending branch in the polar region is accompanied by radiative cooling.

The isentropic mean circulation is approximated in more familiar geometric coordinates by the transformed Eulerian-mean (TEM) residual circulation, resulting in an analogous single-cell circulation from equator to pole (e.g., Andrews and McIntyre 1976, 1978; Edmon et al. 1980; Tung 1986; Andrews et al. 1987; Iwasaki 1989; McIntosh and McDougall 1996; Juckes 2001; Plumb and Ferrari 2005; Pauluis et al. 2011). In the TEM formulation, the residual circulation is generally formulated as the Eulerian-mean circulation plus an eddy term that corresponds to a wave-driven Stokes drift in the small-amplitude limit. Held and Schneider (1999) showed that owing to nearly neutral static stability in the boundary layer, the equatorward flow in the conventional TEM framework of Andrews and McIntyre (1976) does not close at the ground, implying a very thin near-surface layer. This contrasts with the isentropic circulation that an equatorward flow exists in a finite-depth isentropic layer where the layer intersects with the ground (Held and Schneider 1999). Several alternative approximations of isentropic circulation were proposed, such as considering the meridional temperature gradient rather than the static stability in the original TEM formulation (Held and Schneider 1999) or partitioning the eddy fluxes into adiabatic and diabatic components (Plumb and Ferrari 2005). More recently, Pauluis et al. (2011) generalize the TEM formulation to nonmonotonic vertical coordinates by assuming a Gaussian joint distribution for the meridional wind and the state variable of interest.

Another important aspect of the TEM circulation is the corresponding eddy forcing of angular momentum transport. Under the quasigeostrophic approximation, the eddy forcing of angular momentum can be expressed as the well-known Eliassen–Palm (EP) flux divergence (e.g., Andrews and McIntyre 1976; Edmon et al. 1980), corresponding to the nonacceleration theorem (e.g., Charney and Drazin 1961). This relationship can be generalized to finite-amplitude eddies in the isentropic coordinates on the sphere, treating the isentropic mean circulation as the analog for the TEM residual circulation (e.g., Andrews 1983; Tung 1986; Iwasaki 1989, 1998). Further, using a coordinate-independent formulation of Andrews and McIntyre (1978), Plumb and Ferrari (2005) generalized the TEM theory for nonquasigeostrophic eddies in geometric coordinates, and within this framework, Kuo et al. (2005) analyzed the potential vorticity (PV) homogenization in a cylinder flow.

In this paper, we present a new hybrid isobaric–isentropic diagnostic of the mean circulation and corresponding eddy forcing using the mass above the isentropes at each latitude and time as the vertical coordinate. Given the fact that the mass above the isentropic surface decreases monotonically with increasing potential temperature (i.e., positive static stability), we can define a new vertical coordinate corresponding to the mass above the isentrope of interest, termed as “equivalent pressure,” analogous to the familiar concept of equivalent latitude (Butchart and Remsberg 1986). This is essentially equivalent to changing the vertical coordinate of the isentropic formulation from potential temperature to the isentropic zonal-mean pressure in Iwasaki (1989) and Iwasaki (1998). However, our hybrid diagnostic framework is built on the pressure coordinates and can be readily generalized from dry isentropes to other quasi-conservative quantities by construction. Further, the relationship between the TEM and isentropic circulation is clearer in our hybrid formulation, and the problem of conventional TEM streamfunction near the ground can be remedied, at least partially, by replacing the Eulerian static stability with isentropic static stability.

The paper is organized as follows. We first introduce the basic formulations of hybrid isobaric–isentropic diagnostics in section 2. The isentropic mass-weighted mean static stability and meridional wind are discussed in comparison with their isobaric counterparts in section 3. In section 4, the thermodynamic balance between the mean circulation and entropy sources and sinks is derived and the mean streamfunction is presented for the reanalysis data. The corresponding angular momentum budget is investigated in section 5. The conclusions and discussion are provided in section 6. Some detailed derivations are offered in the appendices.

## 2. Formulation of hybrid isobaric–isentropic diagnostics

### a. Basic formulation

We first introduce the formulation of the diagnostics using the mass above isentropes as the vertical coordinate. This is equivalent to changing the vertical coordinate of the isentropic formulation from the potential temperature (e.g., Andrews 1983; Tung 1986) to the isentropic zonal-mean pressure (e.g., Iwasaki 1989, 1998; Tanaka et al. 2004). However, our formulation is derived from the pressure coordinate by treating entropy as a quasi-conservative tracer, analogous to the modified Lagrangian-mean formulation of Nakamura (1995) using a tracer as the meridional coordinate. The formulation can be readily generalized to moist entropy or other quasi-conservative tracers by construction.

*σ*by a Heaviside step function

*p*−

_{s}*p*) for the massless underground:

*p*is the surface pressure, and the zonal mean

_{s}*p*at the latitude

*φ*and time

*t*above the ground.

*θ*= Θ at

*φ*and

*t*(shading in Fig. 1a) can be written as

*dλ*is the isentropic displacement in longitude and

_{θ}*θ*− Θ) ensures the mass above the isentropic surface Θ, and

*σ*for the surface boundary condition.

*m*(

*φ*, Θ,

*t*) can be used to define a new vertical coordinate

*p*—the pressure of the isentropic line after a reversible (mass- and entropy-conserved) zonalization (e.g., from Fig. 1a to Fig. 1b):

_{e}*p*can be termed as the “Θ-equivalent pressure,” analogous to the familiar concept of the equivalent latitude for the diagnostic of isentropic transport (Butchart and Remsberg 1986). At each latitude and time, this yields a one-to-one correspondence of Θ = Θ[

_{e}*p*(Θ)] by

_{e}*p*

_{↑}, which is used as the vertical coordinate in the isentropic formulation of Iwasaki (1998),

*p*is identical to

_{e}*p*

_{↑}if the isentropic surface does not intersect the ground. But

*p*can differ from

_{e}*p*

_{↑}when the isentropic surface intersects the ground, as

*p*

_{↑}does not include the step function

*σ*in the formulation.

*X*, we can define the mass-weighted integral above the isentropic surface as (shading in Fig. 1a)

*θ*= Θ contour:

*X*. In practice,

*X*) with respect to

*m*, and therefore it is defined by quantities in the pressure coordinates.

*σ*, the isobaric mean can be written as the mass-weighted zonal mean at the pressure surface

*p*=

*p*(Fig. 1b), and an eddy term can be defined as the deviation therefrom. If the pressure surface does not intersect the ground, they reduce to the regular zonal mean and eddy terms:

_{e}Hereafter we set the pressure level in the isobaric mean *p* = *p _{e}*). For the isentropic mean, this results from a reversible (mass- and entropy-conserved) vertical adjustment of the isentropic line toward a zonally symmetric basic state, and

*p*has a one-to-one correspondence to the mass above the isentrope Θ by Eq. (4). As such, we construct a hybrid isobaric–isentropic diagnostic using

_{e}*p*=

*p*as the vertical coordinate.

_{e}### b. Data and analysis method

In this study, we employ the National Centers for Environmental Prediction (NCEP)–U.S. Department of Energy (DOE) Reanalysis 2, an updated product of the NCEP–National Center for Atmospheric Research (NCAR) reanalysis (Kanamitsu et al. 2002). The data are available at the same resolution as the NCEP–NCAR reanalysis, with a 2.5° × 2.5° horizontal resolution and 17 pressure levels. We analyze 6-hourly data for the period of 1979–2011, and the 6-hourly data are expected to capture the contributions of extreme weather events to the mass fluxes better than daily data.

To facilitate the vertical integral or derivative with respect to *p* or *p _{e}*, the meteorological fields are interpolated linearly onto 52 evenly spaced vertical levels from 20 to 1040 hPa with an increment of 20 hPa. By comparing the geopotential height with the surface orography, we determine whether a grid point is above or below the ground (i.e., the value of

*σ*). Consistent with aforementioned formulation, the underground world is set as motionless (i.e.,

*u*=

*υ*=

*ω*= 0) and massless (i.e.,

*σ*= 0). In practice, to ensure a monotonic relationship between Θ and

*p*in Eq. (4), we set a very small number for the mass density below the ground as

_{e}*σ*= 10

^{−20}, making it easier in coding than dealing with a constantly moving lower boundary. In terms of entropy, to be consistent with the isentropic coordinates in which the underground mass density −(1/

*g*) ∂

*p*/∂

*θ*is zero, the underground potential temperature is interpolated downward by assuming a very large value of ∂

*θ*/∂

*p*.

Additionally, in computing the isentropic mass-weighted integral [i.e., Eqs. (2) and (5)] at a latitude and time, we first sort *θ* in an ascending order for all of the points. Then *m* and *X*) can be obtained for each value of Θ by aggregating all of the points with *θ* > Θ. The isentropic mean, defined by Eq. (6), is calculated by a finite-differencing scheme of *X*) with respect to *m*. Finally, the equivalent pressure [equivalently the mass above an isentrope by Eq. (4)] and the corresponding isentropic means are interpolated onto the same 52 pressure levels as the isobaric means to facilitate a hybrid isobaric–isentropic diagnostic.

We illustrate the usefulness of our diagnostic using the example of weather at 0000 UTC 28 December 2010. In Fig. 2a, the mean sea level pressure (MSLP) displays a notable extratropical storm at about 40°N, 60°W, known as the December 2010 North American blizzard. The longitude–pressure cross section of the storm at 40°N (Fig. 2b) displays a poleward (solid line) surface wind on the east of 60°W associated with warm air and an equatorward (dashed line) surface wind on the west associated with cold air. Similar poleward/warm-air and equatorward/cold-air correlations are evident elsewhere in the lower troposphere, as expected from the poleward heat transport by extratropical storms. In Figs. 2c–e, the isobaric means (black) and isentropic means (red) are compared. The isobaric zonal-mean meridional wind is characterized by an equatorward wind at 200 hPa and a poleward wind near the surface—a typical circulation pattern of the Ferrel cell. However, when averaging with respect to the isentropic surface and then transformed to the equivalent pressure, the mean meridional wind exhibits a poleward flow within 500–850 hPa and an equatorward flow below 850 hPa, where a fraction of pressure surface is underground, as indicated by *σ* < 1. This corresponds to the single overturning cell of the isentropic mean circulation. The two equivalent-pressure layers of 500–850 and 850–1000 hPa approximately correspond to the two isentropic layers of 290–310 and 270–290 K (Fig. 2e), and in Fig. 2b, the warmer layer portrays a Rossby wave pattern and the colder layer is seen to intersect the ground. Interestingly, the isentropic mean shows much colder temperature near the surface than the isobaric mean, which may be attributed to the cold temperature advection associated with the isentropic equatorward flow. In spite of the circulation diagnostic of a single day, the general pattern of mean meridional circulation in equivalent pressure agrees well with the mean circulation in the isentropic coordinate depicted in Held and Schneider (1999).

An example of isobaric and isentropic mass-weighted means at 40°N at 0000 UTC 28 Dec 2010. (a) MSLP (hPa) is shown in shading, and the dashed line denotes 40°N. (b) The longitude–pressure cross section at 40°N, where the shading denotes potential temperature (K) and the contour denotes meridional wind (solid for positive, dashed for negative, and zero omitted; contour interval is 10 m s^{−1}). (c) The probability of pressure surface above the ground at 40°N,

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

An example of isobaric and isentropic mass-weighted means at 40°N at 0000 UTC 28 Dec 2010. (a) MSLP (hPa) is shown in shading, and the dashed line denotes 40°N. (b) The longitude–pressure cross section at 40°N, where the shading denotes potential temperature (K) and the contour denotes meridional wind (solid for positive, dashed for negative, and zero omitted; contour interval is 10 m s^{−1}). (c) The probability of pressure surface above the ground at 40°N,

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

An example of isobaric and isentropic mass-weighted means at 40°N at 0000 UTC 28 Dec 2010. (a) MSLP (hPa) is shown in shading, and the dashed line denotes 40°N. (b) The longitude–pressure cross section at 40°N, where the shading denotes potential temperature (K) and the contour denotes meridional wind (solid for positive, dashed for negative, and zero omitted; contour interval is 10 m s^{−1}). (c) The probability of pressure surface above the ground at 40°N,

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

## 3. Static stability and meridional wind

Here we discuss the isentropic static stability and meridional wind in comparison with their corresponding isobaric means, especially in the context of the near-surface stability and circulation. Figure 3 shows the isentropic mass-weighted mean potential temperature (dash–dotted) and the isentropic mean minus isobaric mean temperature (i.e.,

Time and isentropic mass-weighted mean potential temperature (dash–dotted, contour interval is 20 K) and the difference between isentropic and isobaric means (i.e.,

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Time and isentropic mass-weighted mean potential temperature (dash–dotted, contour interval is 20 K) and the difference between isentropic and isobaric means (i.e.,

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Time and isentropic mass-weighted mean potential temperature (dash–dotted, contour interval is 20 K) and the difference between isentropic and isobaric means (i.e.,

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

*δφ*and

*δp*denote the meridional and vertical parcel displacements, respectively. Figure 4 displays an idealized temperature profile with a wavenumber-1 perturbation, which may be attributed to either meridional and/or vertical parcel perturbation. The isentropic mass density,

(left) Idealized profile of potential temperature (K) with uniform stratification and a wavenumber-1 perturbation: *θ*(*λ*, *p*) = 280 + 60(1 − *p*/10^{5}) + 10 sin(*λ*). The zonalization of the 280-K isentrope is denoted by shading as the area integral above the 280-K isentrope minus the integral above the corresponding *p _{e}* ≈ 950 hPa, as exemplified in Fig. 1c. (right) The corresponding (solid) isobaric mean and (dashed) isentropic mean of the idealized temperature profile.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

(left) Idealized profile of potential temperature (K) with uniform stratification and a wavenumber-1 perturbation: *θ*(*λ*, *p*) = 280 + 60(1 − *p*/10^{5}) + 10 sin(*λ*). The zonalization of the 280-K isentrope is denoted by shading as the area integral above the 280-K isentrope minus the integral above the corresponding *p _{e}* ≈ 950 hPa, as exemplified in Fig. 1c. (right) The corresponding (solid) isobaric mean and (dashed) isentropic mean of the idealized temperature profile.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

(left) Idealized profile of potential temperature (K) with uniform stratification and a wavenumber-1 perturbation: *θ*(*λ*, *p*) = 280 + 60(1 − *p*/10^{5}) + 10 sin(*λ*). The zonalization of the 280-K isentrope is denoted by shading as the area integral above the 280-K isentrope minus the integral above the corresponding *p _{e}* ≈ 950 hPa, as exemplified in Fig. 1c. (right) The corresponding (solid) isobaric mean and (dashed) isentropic mean of the idealized temperature profile.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

The increased stratification in the isentropic mean near the boundaries results from an entropy-conserved zonalization, which retains the extreme value of entropy, whereas the extreme value is smoothed in the zonal mean. For example, while the 280-K value in the zonal-mean temperature corresponds to the pressure level 1000 hPa, the reversible zonalization of the 280-K isentrope corresponds to the equivalent pressure *p _{e}* ≈ 950 hPa. The reversible zonalization of the 280-K isentrope is demonstrated by the shading in Fig. 4 as the area integral above the 280-K isentrope minus the integral above the corresponding equivalent pressure 950 hPa. Interestingly, Eq. (8) becomes

*p*=

*p*in Fig. 1b, corresponding to the mass-weighted integral of light shading minus the integral of dark shading in Fig. 1c:

_{e}*X*) = 0. Therefore, Δ

*X*) is an eddy quantity related to the thermal perturbation.

*υ*, Δ

*υ*) denotes an eddy quantity analogous to the wave activity defined in Nakamura and Zhu (2010). As exemplified in Fig. 1c, high potential temperature in the light shading is carried poleward and low potential temperature in the dark shading is transported equatorward, and therefore, Δ

*υ*) describes the poleward transport of heat. In the small-amplitude limit, it is related to eddy heat flux [derivations are shown in Eq. (A4)] by

*p*/∂Θ gradually approaches zero, and therefore Δ

_{e}*υ*) is closed at the surface. Also as noted in appendix A, the small-amplitude form is not exact when the isentropic layer intersects the ground under the adiabatic condition. Diabatic heating or diffusion is expected near the ground in this case, as discussed in context of the diabatic surface layer in Plumb and Ferrari (2005).

*m*and using Eqs. (4), (6), and (7), one can obtain a simple relationship between the isentropic and isobaric mean meridional winds:

*σ*= 1 and

*υ*) is a finite-amplitude Stokes correction between the isobaric and isentropic mean winds.

We compare Δ*υ*) with its small-amplitude approximations in Fig. 5, using the isobaric and isentropic static stability, respectively, as the two differ near the ground. Similar to the seasonal cycle of eddy heat fluxes, Δ*υ*) is greatest in the winter hemispheres and smallest in the NH summer. Above the near-surface layer, as can be also inferred from the diagnostics of Tanaka et al. (2004), Δ*υ*) and its small-amplitude approximations agree very well, although using the isobaric static stability results in some subtle differences for the 6-hourly data we analyzed. As one moves into the near-surface layers, the formulations start to diverge. While *σ* leads to zero heat flux when the pressure surface is completely below the ground, the nearly neutral Eulerian static stability results in unrealistically large values of

A comparison of Stokes correction Δ*υ*) with its small-amplitude approximations in (left) DJF and (right) JJA. (a),(b) Isobaric means, *υ*). Contour intervals are 10^{4} m s^{−1} Pa. Light shading indicates the pressure surface where the percentage of mass above the ground is between 10% and 90%.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

A comparison of Stokes correction Δ*υ*) with its small-amplitude approximations in (left) DJF and (right) JJA. (a),(b) Isobaric means, *υ*). Contour intervals are 10^{4} m s^{−1} Pa. Light shading indicates the pressure surface where the percentage of mass above the ground is between 10% and 90%.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

A comparison of Stokes correction Δ*υ*) with its small-amplitude approximations in (left) DJF and (right) JJA. (a),(b) Isobaric means, *υ*). Contour intervals are 10^{4} m s^{−1} Pa. Light shading indicates the pressure surface where the percentage of mass above the ground is between 10% and 90%.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Also, Δ*υ*) can account for the meridional mass flux in a nonmonotonic vertical coordinate, in which it is more than an eddy term. As in Pauluis et al. (2008) and Pauluis et al. (2010), the mass streamfunction can be obtained by using the mass above moist isentropes as the vertical coordinate, and the magnitude of moist isentropic streamfunction is found to differ from dry isentropic streamfunction (not shown). Although our definitions in Eqs. (2) and (5) do not require the monotonicity of the vertical coordinate, the mean circulation in a nonmonotonic coordinate is more complex to interpret, as the change from the new coordinate *p _{e}* back to the geometric coordinate

*p*is not a one-to-one correspondence.

For example, Fig. 6 depicts idealized zonally symmetric profiles of meridional wind and equivalent potential temperature, with a moist entropy minimum in the lower troposphere. The isobaric means are denoted by solid lines and isentropic means by dashed lines. While there is no one-to-one correspondence from Θ to *p* because of the nonmonotonic profile, the correspondence from Θ to *p _{e}* can be established, as

*p*corresponds to the total mass of the air parcel with the potential temperature larger than Θ in Eq. (2). In this idealized case, the isobaric mean static stability vanishes at the local minimum

_{e}*p*is nonzero by construction, as shown in the dashed line of Fig. 6. While the circulation is assumed zonally symmetric, the isentropic mean redistributes the meridional wind vertically, as a result of adjusting the nonmonotonic entropy profile to a monotonic basic state. In the small-amplitude form of Eq. (10), the eddy term would vanish. Therefore, the difference of the isentropic and isobaric means in this case does not represent an eddy term, but the zonally symmetric convective adjustment of the atmosphere. This is analogous to the discussion of wave activity associated with nonmonotonic zonally symmetric PV gradient in Fig. 3 of Nakamura and Zhu (2010). Despite the added complexity in interpretation, this demonstrates an advantage of the hybrid isobaric–isentropic formulation in understanding the mean circulation in nonmonotonic vertical coordinates.

_{e}Idealized zonally symmetric profiles of (left) meridional wind and (right) equivalent potential temperature as a function of pressure (hPa). The isobaric mean meridional wind is specified as *υ*(*p*) = −2cos(*πp*/10^{5}) and equivalent potential temperature as *θ _{e}*(

*p*) = 280 + 60(1 −

*p*/10

^{5}) − 30 sin(

*πp*/10

^{5}), denoted by solid lines. The corresponding isentropic means are denoted by dashed lines.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Idealized zonally symmetric profiles of (left) meridional wind and (right) equivalent potential temperature as a function of pressure (hPa). The isobaric mean meridional wind is specified as *υ*(*p*) = −2cos(*πp*/10^{5}) and equivalent potential temperature as *θ _{e}*(

*p*) = 280 + 60(1 −

*p*/10

^{5}) − 30 sin(

*πp*/10

^{5}), denoted by solid lines. The corresponding isentropic means are denoted by dashed lines.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Idealized zonally symmetric profiles of (left) meridional wind and (right) equivalent potential temperature as a function of pressure (hPa). The isobaric mean meridional wind is specified as *υ*(*p*) = −2cos(*πp*/10^{5}) and equivalent potential temperature as *θ _{e}*(

*p*) = 280 + 60(1 −

*p*/10

^{5}) − 30 sin(

*πp*/10

^{5}), denoted by solid lines. The corresponding isentropic means are denoted by dashed lines.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

## 4. Mean circulation and thermodynamic balance

*θ*:

*u*,

*υ*, and

*ω*are zonal, meridional, and pressure velocities;

*θ*due to diabatic heating or turbulent diffusion; and

**∇**

_{φ}is a gradient operator in the longitude–pressure plane defined as

*σ*, the continuity equation in the pressure coordinate can be written as

^{1}

*σ*≡ 1), the continuity equation and streamfunction reduce to their familiar isobaric forms (e.g., Holton 2004).

*θ*= Θ in a longitude–pressure plane is illustrated in Fig. 7. The mass fluxes across an isentropic surface at a latitude are equal to the meridional fluxes above the same isentrope by mass continuity. Following Nakamura (1995) (the derivations are offered in appendix B), the continuity equation in the coordinates (

*φ*, Θ,

*t*) can be written as

The budget of mass above *θ* = Θ at a latitude. Shaded arrows denote the mass fluxes across the Θ surface in the longitude–pressure plane and white arrows denote the corresponding mass fluxes across the latitude by mass continuity.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

The budget of mass above *θ* = Θ at a latitude. Shaded arrows denote the mass fluxes across the Θ surface in the longitude–pressure plane and white arrows denote the corresponding mass fluxes across the latitude by mass continuity.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

The budget of mass above *θ* = Θ at a latitude. Shaded arrows denote the mass fluxes across the Θ surface in the longitude–pressure plane and white arrows denote the corresponding mass fluxes across the latitude by mass continuity.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

*φ*, Θ,

*t*) to the equivalent pressure (

*φ*,

*p*,

_{e}*t*),

*m*/∂Θ)

_{φ,t}yields the thermodynamic energy equation in the coordinates (

*φ*,

*p*,

_{e}*t*):

*φ*,

*p*,

_{e}*t*) are defined as

Equation (18) is an exact balance between the mean circulations in equivalent-pressure coordinates *D*/*Dt* in the (*φ*, *p _{e}*) plane. Since there is no eddy term in Eq. (18), the mean circulation is essentially equivalent to the isentropic circulation in Eq. (2.24) of Tanaka et al. (2004) in the interior of the atmosphere, which is derived by a change in vertical coordinate from the isentropic circulation of Andrews (1983).

*σ*in the definition of

*p*. The nondivergent mass streamfunction can be obtained by the vertical integral of the meridional mass flux from the top of the atmosphere to the Θ contour:

_{e}^{Θ}is a finite-amplitude extension of the conventional TEM residual streamfunction (e.g., Andrews and McIntyre 1976; Edmon et al. 1980), and

*σ*= 1.

The mean mass streamfunction in equivalent pressure is displayed in Fig. 8 for DJF and JJA. As the TEM residual circulation (Edmon et al. 1980) or the isentropic mean circulation (e.g., Held and Schneider 1999; Tanaka et al. 2004), there exists a single-cell circulation in each hemisphere. The tropical circulation describes a solstitial Hadley cell circulation, characterized by an intense updraft within 10° latitude of the summer hemisphere, crossing the equator in the upper troposphere, and then descending in the subtropics of the winter hemisphere. The intensity^{2} of tropical circulation is stronger in JJA than DJF, possibly owing to larger diabatic heating associated with the monsoons in the NH. The extratropical circulation rises approximately along the isentropic surface, for example, moving along the 300-K isentropic surface from the subtropical lower troposphere to extratropical upper troposphere. This suggests that the entropy source/sink in Eq. (18) is secondary in the extratropical middle-to-upper troposphere. The intensity of extratropical circulation peaks at a similar magnitude in both winter hemispheres, and the intensity is weakest in the NH summer because of weak baroclinic eddies (cf. the heat fluxes in Fig. 5). The high-latitude subsidence moves equatorward near the surface where the isentropes intersect the ground, completing the mean circulation within a finite-depth near-surface layer. These circulation patterns agree well with those in the isentropic coordinates.

Isentropic mass-weighted mean potential temperature (dash–dotted; contour interval is 20 K) and mass streamfunction (solid for positive and dashed for negative; contour interval is 20 × 10^{9} kg s^{−1}) in (left) DJF and (right) JJA. Light shading indicates the pressure surface where the percentage of mass above the ground is between 10% and 90%.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Isentropic mass-weighted mean potential temperature (dash–dotted; contour interval is 20 K) and mass streamfunction (solid for positive and dashed for negative; contour interval is 20 × 10^{9} kg s^{−1}) in (left) DJF and (right) JJA. Light shading indicates the pressure surface where the percentage of mass above the ground is between 10% and 90%.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Isentropic mass-weighted mean potential temperature (dash–dotted; contour interval is 20 K) and mass streamfunction (solid for positive and dashed for negative; contour interval is 20 × 10^{9} kg s^{−1}) in (left) DJF and (right) JJA. Light shading indicates the pressure surface where the percentage of mass above the ground is between 10% and 90%.

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

As in Tanaka et al. (2004), the key distinction of our hybrid isobaric–isentropic diagnostics from the isentropic diagnostics lies in that the time or latitudinal integral/derivative is evaluated at the surface of constant mass rather than constant entropy, although the zonal average is both calculated on constant potential temperature surface. When using mass as the vertical coordinate, the majority of air mass lies above the mean surface pressure; it is then a more convenient coordinate to visualize the mean circulation than the isentropic coordinate.

## 5. Angular momentum balance and Eliassen–Palm flux

*N*denotes the surface friction and other momentum source/sink. The absolute vorticity

*υ*), an eddy quantity related to eddy heat fluxes, is calculated from the difference between the isentropic and isobaric means.

*σ*= 1 and

The angular momentum balance is presented for the reanalysis in DJF and JJA (Fig. 9). The EP flux divergence is mostly confined in the near-surface layers, and the corresponding convergence lies within the isentropic layers aloft. As the EP flux convergence in the interior of the atmosphere reduces to the TEM formulation in the small-amplitude limit, it may be regarded as a finite-amplitude extension of the PV flux. The EP flux convergence aloft aligns approximately along the isentropic surfaces, which can be thought of as the result of the isentropic mixing of PV, analogous to the isentropic mixing of passive tracers found in geometric coordinates (e.g., Plumb and Mahlman 1987). As the EP flux divergence and convergence largely cancel each other in the mass-weighted vertical average, their magnitudes must be inversely proportional to their masses. Indeed, the ratio of the divergence versus convergence is about a factor of 4–5 in the SH and 2–3 in the NH, and this is consistent the ratio of masses between the free troposphere and near-surface layers.

The angular momentum balance in (left) DJF and (right) JJA. (a),(b) Isentropic mass-weighted mean potential temperature (dash–dotted; contour interval is 20 K) and EP flux divergence (solid for positive and dashed for negative; contour interval is 4 m s^{−1} day^{−1} and, for values larger than 30 m s^{−1} day^{−1}, is 8 m s^{−1} day^{−1}). (c),(d) As in Fig. 8, but for potential temperature and mass streamfunction diagnosed from the EP flux divergence. (e),(f) Zonal-mean zonal wind (dash–dotted; contour interval is 5 m s^{−1}) and the residual of angular momentum balance (solid for positive and dashed for negative; contour interval is 4 m s^{−1} day^{−1}). Both the EP flux divergence and residual are multiplied by

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

The angular momentum balance in (left) DJF and (right) JJA. (a),(b) Isentropic mass-weighted mean potential temperature (dash–dotted; contour interval is 20 K) and EP flux divergence (solid for positive and dashed for negative; contour interval is 4 m s^{−1} day^{−1} and, for values larger than 30 m s^{−1} day^{−1}, is 8 m s^{−1} day^{−1}). (c),(d) As in Fig. 8, but for potential temperature and mass streamfunction diagnosed from the EP flux divergence. (e),(f) Zonal-mean zonal wind (dash–dotted; contour interval is 5 m s^{−1}) and the residual of angular momentum balance (solid for positive and dashed for negative; contour interval is 4 m s^{−1} day^{−1}). Both the EP flux divergence and residual are multiplied by

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

The angular momentum balance in (left) DJF and (right) JJA. (a),(b) Isentropic mass-weighted mean potential temperature (dash–dotted; contour interval is 20 K) and EP flux divergence (solid for positive and dashed for negative; contour interval is 4 m s^{−1} day^{−1} and, for values larger than 30 m s^{−1} day^{−1}, is 8 m s^{−1} day^{−1}). (c),(d) As in Fig. 8, but for potential temperature and mass streamfunction diagnosed from the EP flux divergence. (e),(f) Zonal-mean zonal wind (dash–dotted; contour interval is 5 m s^{−1}) and the residual of angular momentum balance (solid for positive and dashed for negative; contour interval is 4 m s^{−1} day^{−1}). Both the EP flux divergence and residual are multiplied by

Citation: Journal of the Atmospheric Sciences 70, 7; 10.1175/JAS-D-12-0239.1

Figures 9c and 9d show the mass streamfunction diagnosed from the EP flux divergence from the angular momentum balance above. It agrees well with the mass streamfunction calculated directly except in the deep tropics and near-surface layer (cf. Figs. 8 and 9). In the deep tropics, the downward control diagnostic breaks down, since

## 6. Conclusions and discussion

In this paper, we present the mean meridional circulation of the atmosphere using the mass above the isentrope of interest as the vertical coordinate. In this vertical coordinate, the mass-weighted mean circulation is exactly balanced by entropy sources and sinks with no eddy flux contribution as in the isentropic coordinate (e.g., Andrews 1983; Tung 1986; Iwasaki 1998), and the coordinate can be readily generalized to the mass above moist isentropes or other quasi-conservative tracers, as in the tracer-based coordinate in the modified Lagrangian-mean diagnostic of Nakamura (1995). We also illustrate the applicability of this framework to an idealized nonmonotonic moist entropy profile, and it would be interesting to apply this diagnostic to the mean circulation in the moist isentropic coordinates of Pauluis et al. (2008, 2010). The new framework here is not restricted by the Gaussian statistics assumed by Pauluis et al. (2011) in dealing with nonmonotonic vertical coordinates.

It is shown in the NCEP–DOE Reanalysis 2 that the new formulation resolves the deficiency of the conventional TEM formulation for the near-surface flow as well as converges to the conventional TEM in the free troposphere. In the small-amplitude limit, the hybrid isobaric–isentropic formulation reduces to the TEM formulation. Therefore, the key improvement near the surface can be partially attributed to the isentropic static stability [the isentropic mass density

The corresponding EP flux divergence for the zonal-mean angular momentum is formulated in a hybrid isobaric–isentropic form, extending the conventional TEM formulation (Andrews and McIntyre 1976; Edmon et al. 1980) to finite-amplitude nongeostrophic eddies on the sphere. This has a straightforward connection to the TEM counterpart in comparison with the isentropic formulation (Andrews 1983; Iwasaki 1998) or the geometric counterpart (Plumb and Ferrari 2005). Following the downward control diagnostic of Haynes et al. (1991), the mean mass streamfunction can be approximately obtained from the EP flux divergence except for the deep tropics or the near-surface flow, highlighting the dominant control of potential vorticity mixing for the subtropics-to-pole mean circulations. It is then expected to provide a valuable diagnostic framework not only for global circulation theory but also for atmospheric transport in the troposphere.

## Acknowledgments

I thank Alan Plumb and Noboru Nakamura for valuable discussions throughout this work. I am also grateful for Olivier Pauluis and anonymous reviewers whose comments led to substantial improvements of the manuscript. The author is supported by the National Science Foundation (NSF) Climate and Large-Scale Dynamics program under Grants AGS-1042787 and AGS-1064079.

## APPENDIX A

### Small-Amplitude Limit of Δ (*υ*)

*p*by a small displacement

_{e}*δp*(

_{e}*λ*) =

*p*(

*λ*, Θ) −

*p*, where

_{e}*p*(

*λ*, Θ) is the pressure at Θ, as illustrated by Fig. 1. The isentropic line can intersect the ground or not. For a conservative perturbation, the temperature anomaly can be approximated as

*p*are equal. Using Eq. (9), we have

_{e}*υ*) in Eq. (9) can be written as

*υ*) is proportional to the eddy heat flux divided by the isentropic static stability parameter in the small-amplitude limit. In this formula, the static stability is evaluated normal to an isentropic surface rather than an isobaric surface, but the eddy flux is evaluated at constant pressure. The difference between the isentropic and isobaric static stability is discussed in section 3.

It should be noted that the approximation in Eq. (A1) is not exact when the isentropic surface intersects the ground under the adiabatic condition. For example, in the left panel of Fig. 4, *δp _{e}* is constant for most of the light-shaded area, and therefore Eq. (A1) is insufficient to describe the temperature anomaly. Diabatic heating or diffusion is expected in this case, as discussed in context of the diabatic surface layer in Plumb and Ferrari (2005).

## APPENDIX B

### Continuity Equation [Eq. (16)]

*Q*,

*θ*,

*t*), where

*Q*is a tracer for the meridional coordinate and

*φ*, Θ,

*t*), where the potential temperature Θ is for the vertical coordinate and

*υ*is the meridional velocity. The rate of change in mass above the Θ isentrope at a fixed latitude and time is

*φ*, Θ,

*t*) coordinates

## APPENDIX C

### Mean Circulation in Equivalent Pressure [Eq. (19)] and Eliassen–Palm Flux [Eq. (25)]

## REFERENCES

Andrews, D. G., 1983: A finite-amplitude Eliassen–Palm theorem in isentropic coordinates.

,*J. Atmos. Sci.***40**, 1877–1883.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.Andrews, D. G., and M. E. McIntyre, 1978: Generalized Eliassen–Palm and Charney–Drazin theorems for waves on axismmetric mean flows in compressible atmospheres.

,*J. Atmos. Sci.***35**, 175–185.Andrews, D. G., J. R. Holton, and C. B. Leovy, 1987:

*Middle Atmosphere Dynamics.*Academic Press, 489 pp.Butchart, N., and E. E. Remsberg, 1986: The area of the stratospheric polar vortex as a diagnostic for tracer transport on an isentropic surface.

,*J. Atmos. Sci.***43**, 1319–1339.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.Czaja, A., and J. Marshall, 2006: The partitioning of poleward heat transport between the atmosphere and ocean.

,*J. Atmos. Sci.***63**, 1498–1511.Edmon, H. J., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen–Palm cross sections for the troposphere.

,*J. Atmos. Sci.***37**, 2600–2616.Haynes, P. H., C. J. Marks, M. E. McIntyre, T. G. Shepherd, and K. P. Shine, 1991: On the “downward control” of extratropical diabatic circulations by eddy-induced mean zonal forces.

,*J. Atmos. Sci.***48**, 651–678.Held, I. M., and T. Schneider, 1999: The surface branch of the zonally averaged mass transport circulation in the troposphere.

,*J. Atmos. Sci.***56**, 1688–1697.Holton, J. R., 2004:

*An Introduction to Dynamic Meteorology.*Academic Press, 535 pp.Iwasaki, T., 1989: A diagnostic formulation for wave-mean flow interactions and Lagrangian-mean circulation with a hybrid vertical coordinate of pressure and isentropes.

,*J. Meteor. Soc. Japan***67**, 293–312.Iwasaki, T., 1998: A set of zonal mean equations in a pressure–isentrope hybrid vertical coordinate.

,*J. Atmos. Sci.***55**, 3000–3002.Juckes, M., 2001: A generalization of the transformed Eulerian-mean meridional circulation.

,*Quart. J. Roy. Meteor. Soc.***127**, 147–160.Juckes, M., I. N. James, and M. Blackburn, 1994: The influence of Antarctica on the momentum budget of the southern extratropics.

,*Quart. J. Roy. Meteor. Soc.***120**, 1017–1044.Kanamitsu, M., W. Ebisuzaki, J. Woollen, S.-K. Yang, J. J. Hnilo, M. Fiorino, and G. L. Potter, 2002: NCEP–DOE AMIP-II Reanalysis (R-2).

,*Bull. Amer. Meteor. Soc.***83**, 1631–1643.Koh, T.-Y., and R. A. Plumb, 2004: Isentropic zonal average formalism and the near-surface circulation.

,*Quart. J. Roy. Meteor. Soc.***130**, 1631–1653.Kuo, A., R. A. Plumb, and J. Marshall, 2005: Transformed Eulerian-mean theory. Part II: Potential vorticity homogenization and the equilibrium of a wind- and buoyancy-driven zonal flow.

,*J. Phys. Oceanogr.***35**, 175–187.McIntosh, P. C., and T. J. McDougall, 1996: Isopycnal averaging and the residual mean circulation.

,*J. Phys. Oceanogr.***26**, 1655–1660.Nakamura, N., 1995: Modified Lagrangian-mean diagnostics of the stratospheric polar vortices. Part I. Formulation and analysis of GFDL SKYHI GCM.

,*J. Atmos. Sci.***52**, 2096–2108.Nakamura, N., and D. Zhu, 2010: Finite-amplitude wave activity and diffusive flux of potential vorticity in eddy-mean flow interaction.

,*J. Atmos. Sci.***67**, 2701–2716.Pauluis, O., A. Czaja, and R. Korty, 2008: The global atmospheric circulation on moist isentropes.

,*Science***321**, 1075–1078.Pauluis, O., A. Czaja, and R. Korty, 2010: The global atmospheric circulation in moist isentropic coordinates.

,*J. Climate***23**, 3077–3093.Pauluis, O., T. Shaw, and F. Laliberté, 2011: A statistical generalization of the transformed Eulerian-mean circulation for an arbitrary vertical coordinate system.

,*J. Atmos. Sci.***68**, 1766–1783.Peixoto, J. P., and A. H. Oort, 1992:

*Physics of Climate.*AIP Press, 520 pp.Plumb, R. A., and J. D. Mahlman, 1987: The zonally averaged transport characteristics of the GFDL general circulation/transport model.

,*J. Atmos. Sci.***44**, 298–327.Plumb, R. A., and R. Ferrari, 2005: Transformed Eulerian-mean theory. Part I: Nonquasigeostrophic theory for eddies on a zonal-mean flow.

,*J. Phys. Oceanogr.***35**, 165–174.Schneider, T., 2005: Zonal momentum balance, potential vorticity dynamics, and mass fluxes on near-surface isentropes.

,*J. Atmos. Sci.***62**, 1884.Simmons, A. J., and D. M. Burridge, 1981: An energy and angular-momentum conserving vertical finite-difference scheme and hybrid vertical coordinates.

,*Mon. Wea. Rev.***109**, 758–766.Stachnik, J. P., and C. Schumacher, 2011: A comparison of the Hadley circulation in modern reanalyses.

,*J. Geophys. Res.***116**, D22102, doi:10.1029/2011JD016677.Tanaka, D., T. Iwasaki, S. Uno, M. Ujiie, and K. Miyazaki, 2004: Eliassen–Palm flux diagnosis based on isentropic representation.

,*J. Atmos. Sci.***61**, 2370–2383.Townsend, R. D., and D. R. Johnson, 1985: A diagnostic study of the isentropic zonally averaged mass circulation during the First GARP Global Experiment.

,*J. Atmos. Sci.***42**, 1565–1579.Tung, K. K., 1986: Nongeostrophic theory of zonally averaged circulation. Part I: Formulation.

,*J. Atmos. Sci.***43**, 2600–2618.

^{1}

More generally, the continuity equation in the *η* (a monotonic function of pressure *p*) vertical coordinate on the sphere can be written as

^{2}

One should note a difference in magnitude between this study and Fig. 1 of Tanaka et al. (2004). This is due to the difference between the NCEP–NCAR and NCEP–DOE reanalyses (cf. Stachnik and Schumacher 2011).