1. Introduction
The global atmospheric circulation can transport air masses and constituents by either mean advection or eddy diffusion. Using a Green’s function diagnostic, Bowman and collaborators (Bowman and Carrie 2002; Bowman and Erukhimova 2004; Erukhimova and Bowman 2006) concluded that tropospheric transport can be described qualitatively by the slow mean diabatic circulation and rapid isentropic mixing. In the observations, the mean diabatic circulation exhibits an equator-to-pole cell, with poleward circulations in the stratosphere and in the upper troposphere and an equatorward return flow near the surface (e.g., Tanaka et al. 2004).
The isentropic mixing in tropospheric transport is not well understood. This is largely due to the complex boundary condition that nearly half of the isentropes in the troposphere frequently intersect the ground. This near-surface branch of isentropic flow is important for the transport of air pollutants from midlatitude industrialized regions to Arctic upper troposphere. Furthermore, even for the isentropic flow that does not intersect the ground, the recent study by Birner et al. (2013) emphasized the upgradient potential vorticity (PV) fluxes near the subtropical jet, which poses a challenge to the diffusive closure theory of eddy PV fluxes (e.g., Held and Larichev 1996; Lapeyre and Held 2003). The focus of this study is to develop a quantitative framework to assess the effective diffusivity of isentropic mixing in tropospheric transport, especially for the near-surface isentropic flow that frequently intersects the ground.
Much of our understanding of isentropic mixing has been developed by using a quasi-conservative tracer as the meridional coordinate, which separates irreversible mixing from reversible undulation of tracer contours (e.g., Butchart and Remsberg 1986; Nakamura 1995). More specifically, Nakamura (1996) showed that a passive tracer advected by a 2D nondivergent flow with pure horizontal diffusion can be formulated in a diffusive form and that the “effective diffusivity” can be defined by the mass flux across the tracer contours. Using a passive tracer driven by the nondivergent isentropic winds in the stratosphere and in the upper troposphere, the effective diffusivity can characterize enhanced mixing within the surf zone and mixing barriers at the subtropics and at the polar vortex edge (Haynes and Shuckburgh 2000a,b; Allen and Nakamura 2001). This diagnostic can be further generalized to the 3D flow, leading to an exact advection–diffusion equation with mass and entropy as the meridional and vertical coordinates, respectively (Nakamura 1998). More recently, Leibensperger and Plumb (2014) extended Nakamura’s formalism of effective diffusivity from pure horizontal diffusion to both isentropic and diabatic diffusions. This extension takes into account the vertical parcel displacement for the small-scale diffusion, which is deemed crucial, since the stirred tracer filaments become thinner and shallower in typical baroclinic flows (Haynes and Anglade 1997). Despite these successes in quantifying isentropic mixing, most studies have assumed that the latitudinal change of isentropic static stability is small, which cannot be applied to the near-surface isentropic flow, where the isentropic air mass converges to zero at the ground.
To circumvent the difficulty in dealing with the near-surface isentropic flow, quasigeostrophic (QG) scaling, in which the isentropic slope is small, is often used. Indeed, studies of the equilibration of baroclinic eddies in the two-layer QG models (e.g., Lee and Held 1993) reveal that the upper troposphere is characterized by wavelike disturbances with enhanced mixing on the jet’s flanks, analogous to the isentropic mixing at the polar night jet. The lower troposphere, in contrast, is dominated by chaotic mixing throughout the baroclinic zone, and the eddy diffusivity can be parameterized by the mean state (e.g., Held and Larichev 1996; Lapeyre and Held 2003). Using the concept of effective diffusivity, Greenslade and Haynes (2008) quantified the vertical transition from the wavelike upper-level disturbance to lower-level chaotic mixing; they emphasized a sharp vertical transition in effective diffusivity and that the transitional height varies with external parameters such as the beta effect, thermal relaxation rate, and surface friction. Furthermore, Greenslade and Haynes (2008) argued that the vertical transition of diffusivity in a QG model, depicted in Fig. 1a, is analogous to that in a primitive equation model, as illustrated in Fig. 1b. These characteristics are generally consistent with the structure of eddy diffusivity inferred from a gradient–flux relationship of tracer transport (e.g., Plumb and Mahlman 1987) or PV mixing (e.g., Jansen and Ferrari 2013).
Schematic diagrams of eddy mixing in (a) quasigeostrophic flow and (b) adiabatic flow. Solid lines denote isobars in (a) and isentropes in (b). Dashed lines and light shading indicate a westerly jet. Arrows denote the directions and strengths of eddy mixing.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1


We first extend previous formalism of effective diffusivity (Haynes and Shuckburgh 2000a,b; Allen and Nakamura 2001) to deal with the intersections of isentropes with the ground. A new definition of effective diffusivity from the small-scale diffusion experienced by a tracer is introduced to account for both horizontal and vertical parcel displacements. Next, passive tracers are used to diagnose the effective diffusivities induced by the QG flow and adiabatic flow simulated in a dry primitive equation model on the sphere. More specifically, the QG diffusivity is obtained by advecting a tracer with the geostrophic flow at each model level. The adiabatic flow is approximated by eliminating the transformed Eulerian-mean (TEM) residual circulation from the advecting velocities, and therefore the residual tracer transport is expected to be parallel to the isentropic surfaces (e.g., Andrews and McIntyre 1978; Plumb and Ferrari 2005). It is found that while the directions of eddy mixing and the resulting tracer isopleth slopes are different for the QG and adiabatic flows, the patterns of the QG and isentropic effective diffusivities agree quite well. This directly supports the analogy between the QG and isentropic diffusivities. Furthermore, in spite of the upgradient QG PV fluxes near the jet core, the time-mean diffusivities are positive near the jet, implying downgradient eddy tracer fluxes.
The paper is organized as follows. The theory of effective diffusivity is described in section 2. The primitive equation model and idealized numerical experiments are presented in section 3. The diffusivity diagnostics are applied to the tracer transport in the idealized simulations in section 4. Concluding remarks are provided in section 5. Some detailed derivations are provided in the appendix.
2. Theory of effective diffusivity
a. Modified Lagrangian-mean (MLM) diagnostics




Next, we formulate an MLM tracer equation as in Nakamura (1998). As the isentropic MLM formalism of Nakamura is based on mass continuity, which does not depend on the vertical coordinate used, the MLM formalism of the η coordinates can be obtained by changing the symbols from potential temperature θ and
Schematics of a quasi-conservative tracer q associated with an idealized wave-breaking event in the polar stereographic projection. The tracer is advected by pure horizontal motions on a constant η surface (i.e.,
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1




b. Definition of effective diffusivity with tracer gradients



c. Definition of effective eddy diffusivity with small-scale diffusion





This effective eddy diffusivity is explained schematically in Fig. 2d. The gray arrows depict the diffusive mass flux across the tracer contour, (Keff/a)(∂Q/∂ϕe), and the white arrows depict the zonal mean diffusive flux across the latitude,


For adiabatic flow, Nakamura (1998) also defined effective diffusivity using the small-scale diffusion experienced by a tracer. However, our definition separates the diffusivity acting on the mean flow from the diffusivity on the eddy, which is important for a tracer dominated by the small-scale vertical diffusion. That is, the effective isentropic diffusivity should be separate from the effective diabatic diffusivity. More importantly for the small-scale vertical (diabatic) diffusivity, we have assumed that (i) the antidiagonal terms in the effective diffusivity tensor, as formulated in Leibensperger and Plumb (2014), can be ignored, and (ii) the effective diabatic diffusivity is not augmented from the small-scale diabatic diffusivity even in the presence of stirring [formulated as




This can be thought of as the small-amplitude approximation of the effective diffusivities for the horizontal diffusion in Eq. (17) or for the vertical diffusion examined in Leibensperger and Plumb (2014). In this limit, the effective diffusivity appears to be nonnegative. However, since the variance of
d. Eddy diffusive closure

This can be used to define the mean advection of wave activity, shown schematically in Fig. 2c. The mean advection across ϕe, in the flux form of
Hereby we have extended the QG formalism [Eq. (2)] to the effective eddy diffusivity with pure horizontal motions in a general vertical coordinate η. If q denotes the Ertel PV, this also completes the isentropic form of the nonacceleration theorem in Nakamura and Solomon (2011) by the concept of effective diffusivity. Compared with the Eulerian form (e.g., Jansen and Ferrari 2013), the advective flux in the mass equivalent latitude is expressed in a more straight-forward form with the mean meridional mass flux. In the time average, we see that the isentropic eddy flux is maintained by the diffusive flux along the Lagrangian gradient plus an additional advective flux of wave activity. By corollary, the eddy flux is downgradient provided that the advective flux of wave activity is small.
3. Dry primitive equation model
The model is run at the horizontal resolution T42 and 20 equally spaced sigma vertical levels, and the results are presented in the figures with the sigma level multiplied by 1000 hPa as the vertical coordinate. The model is integrated for 2000 days to reach a statistical equilibrium of baroclinic waves and the next 300 days of the simulation is used to advect passive tracers for diffusivity diagnoses. Passive spectral tracers, advected by the velocities created by the model, are diffused by the same ∇8 hyperdiffusion as the dynamical fields such as PV or potential temperature. This helps to determine whether the upgradient eddy PV fluxes, if simulated in the model, are associated with negative effective diffusivity or not, since the effective diffusivity of the regular diffusion is positive definite by design. Furthermore, the evolution of the spectral tracer is checked with the grid tracer advected by the same flow with a semi-Lagrangian scheme for horizontal advection (Lin et al. 1994) and a finite-volume parabolic scheme for vertical advection. Since the small-scale diffusion experienced by the grid tracer is implicit, we can only calculate the equivalent length of the tracer contours. Therefore, most of our results are shown for the spectral tracer subject to the hyperdiffusion (except for the grid tracer results in Figs. 12 and 13). A summary of the idealized experiments of tracer transport is documented in Table 1 and described in detail as follows.
Next, the geostrophic flow is obtained by taking the nondivergent component of the pressure-weighted horizontal velocity (i.e.,
The effective diffusivities are calculated from the tracers advected by the QG flow and adiabatic flow, and additional sensitivity experiments are performed with respect to initial tracer distributions and advecting flows as in Haynes and Shuckburgh (2000a). Three initial distributions [i.e., sin(ϕ), sin3(ϕ), and tanh(2ϕ)] are compared in experiments 3a–3c and experiments 4a–4c. As discussed in Haynes and Shuckburgh (2000a), the initial tracer distribution will adjust to align with the dominant quasi-zonal flow quickly, and the effective diffusivity thereafter is expected to be independent of initial conditions. Furthermore, the strengths of the jets are altered by changing the parameter of hemispherically asymmetric heating from ϵ = 0 to ϵ = 40 K with an interval of 10 K (experiments 5a–5d and 6a–6d). As ϵ gets larger, the jet will be stronger in the winter hemisphere and weaker in the summer hemisphere.
Following Eq. (19), the effective eddy diffusivity is diagnosed with the small-scale diffusion experienced by the spectral tracers. The area integral with respect to tracer contours is performed with a box counting method as Nakamura and Solomon (2011). For the geostrophic flow, the surface pressure is used as the mass density in the sigma coordinates. As the diffusion is applied along the sigma surface, the diffusivity is also confirmed with the diffusivity definition in Eq. (17) using the tracer gradients. Next, for the isentropic diffusivity, all the variables are first interpolated onto 71 isentropic levels ranging from 240 to 380 K with an increment of 2 K. The mass equivalent latitude facilitates a simple treatment of the surface boundary condition for the isentropic surfaces intersecting the ground (Nakamura and Solomon 2011). The effective diffusivity is diagnosed with the interpolated small-scale diffusion, which includes both isentropic and diabatic diffusions. Finally, as the small-scale diffusion for the grid tracer is implicit, the equivalent length is calculated with the tracer contours using Eq. (23).
4. Tracer transport diagnosis
a. Zonal mean tracer isopleths
We first examine the effects of advecting flows on the zonal mean tracer isopleths. Figure 3 shows the climatological-mean zonal wind, potential temperature, TEM residual circulation, and QG PV flux for the experiments with ϵ = 20 K in the radiative equilibrium temperature. In this configuration, the westerly jet is stronger and slightly more equatorward in the NH (~43 m s−1 and 43°N) than in the SH (~24 m s−1 and 46°S). The 300- or 310-K isentropic surface extends approximately from the subtropics at the surface to the high latitudes near the tropopause. The TEM residual circulation rises near 15°S in the summer hemisphere. The northern branch of the residual circulation crosses the equator into the winter hemisphere, descends slightly in the subtropics, and finally moves down to the surface in the high latitudes. The summer cell is similar in structure to the winter cell, but the strength is about half of its winter cell. Furthermore, the QG PV fluxes are predominately positive in the lower troposphere and negative in the upper troposphere, as expected from eddy mixing of the background PV gradient. Interestingly, there are noticeable upgradient PV fluxes near the jet core of the winter hemisphere, as found in Birner et al. (2013). All in all, the idealized simulation captures the basic structures of the observed mean circulations in the solstitial season (cf. Haynes and Shuckburgh 2000b; Tanaka et al. 2004; Chen 2013; Birner et al. 2013), while in the observations the wintertime subtropical jet is more equatorward and the subtropical descent of the residual circulation is more distinct from the midlatitude maximum.
The climatological-mean (a) zonal wind, (b) potential temperature, (c) TEM residual circulation, and (d) quasigeostrophic potential vorticity flux. The contour intervals are 5 m s−1, 10 K, 20 × 109 kg s−1, and 2 m s−1 day−1, respectively. The parameter for radiative equilibrium temperature is ϵ = 20 K. The dashed–dotted lines in (b) indicate the tropopause defined by the WMO lapse-rate definition and the 95% quantile of the surface temperature distribution. Pressure is the vertical coordinate and hereafter corresponds to the sigma level (p/ps) multiplied by 1000 hPa.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
The total velocities are used to advect a passive spectral tracer initialized with a sinϕ distribution, and the zonal mean isopleths on selected days are shown in the left panels of Fig. 4. From days 3 to 30, we see that the initial vertically aligned isopleths (dashed and solid lines) in the extratropics are tilted toward local isentropic surfaces (dashed–dotted lines), and the isopleths in the tropics move upward at 20°S–0° because of the tropical rising branch of the residual circulation (Fig. 3). This results in a local minimum of tracer mixing ratio near 10°N in the lower troposphere. The tracer isopleth slopes remain roughly unchanged from days 30 to 90, although the meridional gradient of tracer mixing ratio on day 90 is much reduced from that on day 30. This suggests that the tracer isopleth slopes may be explained by a balance between the mean advection and eddy mixing as the long-lived tracers in the stratosphere (e.g., Holton 1986; Mahlman et al. 1986). The time scales are consistent with Haynes and Shuckburgh (2000a), who found that the effective diffusivity of the stratosphere is independent of the initial conditions of tracers after a 1-month spinup, and that the effective diffusivity of lower-stratospheric eddy mixing allows the tracer gradients to survive a seasonal cycle. Note that there is no small-scale convection and associated vertical diffusion in this model, and their effects are beyond the scope of this study.
The zonal mean spectral tracer mixing ratio on (top) day 3, (middle) day 30, and (bottom) day 90, advected by (left) the total flow (expt 1) and (right) the TEM residual circulation (expt 2). The interval of (solid and dashed) contours is 0.2. The dashed–dotted lines in the left panels denote the mean potential temperature, as in Fig. 3. See Table 1 for descriptions of experiments.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
The separate effects of mean advection and mixing are illustrated by advecting the tracer with different components of the flow. When the tracer is advected by the TEM residual circulation shown in the right panels of Fig. 4, the zonal mean tracer isopleths from days 3 to 30 are similar to those advected by the total flow. However, the tracer isopleths on day 90 advected by the TEM circulation differ considerably from those advected by the total flow. If the advecting circulation is steady with no diffusion, the tracer mixing ratio would be homogenized along the streamlines in the steady state, and the equilibrium tracer isopleths would resemble the streamlines. Although the TEM circulation is not exactly steady, the tracer isopleths on day 90 approximately follow the streamlines of the time averaged TEM circulation (cf. Figs. 4 and 3). Interestingly, while the initial tracer mixing ratio is positive everywhere in the NH, negative mixing ratios are found in the NH midlatitudes on day 90, which can be traced back to the upper-level air masses crossing the equator on day 30. It should be noted that the same ∇8 numerical diffusion is applied to this experiment, but little effective diffusion occurs since the tracer distribution remains zonally symmetric in the absence of eddy stirring. In comparison with the transport by the total flow, we see that the mean mass circulation is important for the tropical ascent and interhemispheric exchange of tracers, but the extratropical tracer isopleths depend critically on the eddy mixing.
Turning our attention to eddy mixing, the evolution of the zonal mean tracer isopleths are shown in Fig. 5 for the QG flow in the left panels and for the adiabatic flow in the right panels. In both experiments, the tracer isopleth slopes remain approximately unchanged after day 30. For the QG flow, the primary direction of parcel displacements takes place along the model sigma levels, and for the adiabatic flow, the parcel displacements along isentropic surfaces. As such, a weak meridional tracer gradient along the mixing direction [i.e.,
The zonal mean spectral tracer mixing ratio, as in Fig. 4, but advected by (left) the QG flow (expt 3a) and (right) the adiabatic flow (expt 4a). Contour interval is 0.2. The region of weak meridional tracer gradient is shaded for
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
Figure 6 displays snapshots of the horizontal distribution of tracers on day 90 for the QG flow in the left panels and for the adiabatic flow in the right panels. The gray shading highlights the area bounded by the tracer contour with the mass equivalent latitude of 40°N. Consistent with the meridional tracer gradient, the tracer distribution in the upper troposphere (at 225 hPa or 330 K) is wavelike near the jet, and the tracer distribution near the surface (at 775 hPa or 290 K), in contrast, is turbulent with frequent overturning tracer contours. Given the complexity in tracer geometry, a quantitative measure of the flow is needed to understand the effect on tracer transport. It is also important to note that the actual latitude of a tracer contour can be far away from its mass equivalent latitude owing to large displacements of eddies, as exemplified by the mass equivalent latitude of 40°N. The mass equivalent latitude also helps to simplify the treatment of the intersections of isentropes with the ground.
Snapshot of the horizontal distribution of tracers on day 90. (left) The tracers advected by the QG flow at (a) 225 and (c) 775 hPa (expt 3a). (right) The tracers advected by the adiabatic flow at (b) 330 and (d) 290 K (expt 4a). Contour intervals are 0.1 in (a),(b) and 0.05 in (c),(d). The gray shading indicates the area bounded by the tracer contour that corresponds to the mass equivalent latitude of 40°N. The black shading in (d) indicates the underground.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
In summary, while the tropical ascent of isopleths and interhemispheric exchange of tracers in the total flow are captured only by the simulation advected by the TEM residual circulation, the extratropical tracer isopleth slopes are reproduced only in the simulation advected by the adiabatic flow. The simulation with the QG flow cannot yield either the tropical or extratropical tracer isopleth slopes, but the regions of strong and weak mixing are consistent with those advected by the adiabatic flow. These mixing characteristics are investigated in detail in the following section.
b. Effective eddy diffusivities
Here, we compare the effective eddy diffusivity [Eq. (19)] calculated from the spectral tracers advected by the QG and adiabatic flows. As the primary direction of mixing is different, the effective diffusivity is separately evaluated along the sigma level for the QG flow and along the isentropic surface for the adiabatic flow. Figure 7 shows the temporal evolution of the QG and isentropic diffusivities averaged during selected days. In the upper troposphere, the QG diffusivity at 225 hPa in the left panels roughly corresponds to the isentropic diffusivity at 330 K in the right panels (cf. Fig. 3), and the two diffusivities agree, at least qualitatively. When averaged over days 100–300, the diffusivity in each hemisphere displays a midlatitude minimum at the jet core (~45°N or 50°S) with two maxima at the jet’s flanks, consistent with the observed pattern of upper-tropospheric diffusivity (Haynes and Shuckburgh 2000b). The diffusivity minimum near 45°N does not emerge during days 1–5, likely influenced by the tracer initial conditions. After day 30, the midlatitude minimum is persistent and the meridional structure of diffusivity remains qualitatively similar, although the values in selected periods differ largely from the 200-day mean. Furthermore, the weak midlatitude mixing and strong high-latitude mixing are consistent with the strong midlatitude tracer gradients and weak high-latitude gradients in the upper level on day 90 in Figs. 5 and 6. Note that effective eddy diffusivity is presented at the mass equivalent latitude rather than the actual latitude (see Fig. 6 for their difference).
Effective eddy diffusivity (m2 s−1) averaged for days 1–5, 28–32, 88–92, and 100–300. (left) Values are calculated from the tracers advected by the QG flow at (top) 225 and (bottom) 775 hPa (expt 3a). (right) Values are calculated from the tracers advected by the adiabatic flow at (top) 330 K and (bottom) 95% quantile of surface potential temperature distribution (expt 4a). The effective diffusivity is calculated from spectral tracers with Eq. (19).
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
In the lower troposphere, the QG diffusivity at 775 hPa in the left panels is compared with the isentropic diffusivity at the top of the surface layer (i.e., 95% quantile of surface potential temperature distribution) in the right panels (cf. Fig. 3). The two diffusivities are qualitatively similar. When averaged over days 100–300, both the QG and isentropic diffusivities are large throughout the baroclinic zone with peaks near 45°S and 40°N, consistent with the weak tracer gradients in the lower level on day 90 in Figs. 5 and 6. The isentropic diffusivity seems to suggest a second maximum in the high latitudes (~70°S/N), and the two diffusivities averaged over selected periods disagree quantitatively especially in the middle and high latitudes. These discrepancies may be attributed to the mixing associated with the vertical parcel displacements along near-surface isentropes (where the isentropic slopes are large) that are ignored in the QG flow.
The vertical structures of the QG and isentropic diffusivities are compared in Fig. 8, and the WMO tropopause and the top of surface layer are denoted by the green lines. In the NH, the QG diffusivity displays a maximum near 700 hPa and 45°N, corresponding to the isentropic diffusivity maximum near 290 K. The diffusivity decreases with altitude from above the top of the surface layer, reaching a minimum at the jet core or at the tropopause. The decrease with altitude is slower at the jet’s flanks. The diffusivity structure in the SH is similar to that in the NH, but the amplitude is weaker. The characteristics of the QG and isentropic diffusivities in this primitive equation model are consistent with those in the QG simulations (Greenslade and Haynes 2008). The diffusivity in the deep tropics is very weak in this dry model owing to the lack of small-scale convection or large-scale stirring, since the equatorward-propagating Rossby waves would break down and get dissipated near their subtropical critical latitudes (e.g., Randel and Held 1991).
Effective eddy diffusivity (shading; m2 s−1) and zonal wind (black contours; m s−1) averaged for days 100–300. Values are calculated from the tracers advected by (left) the QG flow (expt 3a) and (right) the adiabatic flow (expt 4a). The effective diffusivity is calculated from spectral tracers by Eq. (19). Contour intervals for the zonal winds are 5 m s−1. The green lines indicate the WMO tropopause and the 95% quantile of surface potential temperature distribution θs. The diffusivity for the isentropes less than the 10% quantile of θs is not shown.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
Some quantitative differences are also noticeable between the two diffusivities. The low values of effective diffusivity near the midlatitude jet in the QG case are closely tied to the zonal winds (i.e., the diffusivity minimum in blue corresponds to the zonal winds greater than 30 m s−1), whereas in the adiabatic case, the values of effective diffusivity near the midlatitude jet are high below about 310 K, even when the zonal winds are large.
c. Sensitivities to initial tracer distributions
Shuckburgh and Haynes (2003) showed generally for time-periodic flows, the effective diffusivity of a nondivergent advecting flow is independent of initial conditions of tracers except for very weak initial tracer gradients. The effective diffusivities diagnosed with three different initial tracer distributions are compared in Fig. 9. Indeed, different initial distributions yield nearly identical effective diffusivities for the QG flow and for the adiabatic flow in the upper troposphere. The near-surface isentropic diffusivities show some sensitivities to initial conditions, which may be attributed to strong thermal damping near the surface in the model. As pointed out by Plumb and Ferrari (2005), under strong diabatic heating in the surface layer, the residual eddy tracer transport from removing the TEM circulation may not be exactly along the isentropic surface. In spite of this limitation, the structures of the near-surface isentropic diffusivity with different initial conditions agree qualitatively.
Effective diffusivity (m2 s−1) averaged for days 100–300, as in Fig. 7, but for three initial tracer distributions: sin(ϕ), sin3(ϕ), and tanh(2ϕ). (left) Values are calculated from spectral tracers advected by the QG flow at (top) 225 and (bottom) 775 hPa (expts 3a–3c), and (right) values are calculated from the tracers advected by the adiabatic flow at (top) 330 K and (bottom) 95% quantile of surface potential temperature (expts 4a–4c).
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
d. Sensitivities to jet strength
The impact of the jet strength on the effective diffusivity is also investigated by altering the parameter for hemispherically asymmetric heating (Fig. 10). For ϵ = 0 K, the time-mean zonal wind is symmetric about the equator with the jet maximum around 33 m s−1. For ϵ = 40 K, the zonal wind peaks around 51 m s−1 in the NH and around 18 m s−1 in the SH. In spite of the large difference in the zonal jet speed between the two simulations, we see a robust midlatitude diffusivity minimum at the jet core for both the QG flow and adiabatic flow, surrounded by larger values in diffusivity at the jet’s two flanks. Interestingly, for the case of ϵ = 40 K, the diffusivity maximum on the equatorward flank of the wintertime subtropical jet encroaches into the regions of tropical easterlies, indicating a vigorous nonlinear critical layer. In the lower troposphere (below 500 hPa or 310 K), the diffusivity is characterized by a midlatitude maximum below the jet and the diffusivity spreads throughout the baroclinic zone. When ϵ is not zero, the effective diffusivity is larger in the winter hemisphere than the summer hemisphere, consistent with the observed seasonal difference in eddy diffusivity inferred from the energy transport (e.g., Stone and Miller 1980).
Effective eddy diffusivity (shading; m2 s−1) and zonal winds (black contours; m s−1), as in Fig. 8, but for (top) ϵ = 0 and (bottom) ϵ = 40 K. (left) Values are calculated from spectral tracers advected by (left) the QG flow (expts 5a and 5d), and (right) the adiabatic flow (expts 6a and 6d). Contour intervals for the zonal winds are 5 m s−1. The green lines indicate the WMO tropopause and the 95% quantile of θs.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
The diffusivity maxima at the jet’s flanks and the minimum at the jet core are compared with the zonal jet speed in Fig. 11 for ϵ varied from 0 to 40 K with an interval of 10 K. The contrast between the midlatitude minimum and subtropical and/or subpolar values is larger for the jet with a larger zonal speed, as expected from larger baroclinicity and more energetic baroclinic eddies that break down on the jet’s flanks. Also, a stronger zonal jet can suppress the mixing at the jet core by acting as a more effective mixing barrier [see also Haynes and Shuckburgh (2000b)].
Effective eddy diffusivities at the jet’s flanks and at the jet core vs the zonal jet speed maximum at (left) 225 hPa and (right) 330 K. The results are for the experiments in which ε is 0, 10, 20, 30, and 40 K. Both hemispheres are included, with larger zonal jet speeds for the winter hemisphere. Plus signs and triangles indicate the effective diffusivity maxima on the jet’s equatorward and poleward flanks, respectively, and circles indicate the diffusivity minimum at the jet core [see the top panels of Fig. 9 for examples]. Values are calculated from spectral tracers advected by (left) the QG flow (expts 3a and 5a–5d) and from tracers advected by (right) the adiabatic flow (expts 4a and 6a–6d).
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
While our results agree qualitatively with the observed seasonal cycle of effective diffusivity in Haynes and Shuckburgh (2000b), some quantitative differences in the region of weak mixing can be attributed to the differences in the diagnostic method. Haynes and Shuckburgh (2000b) use the regular diffusion with κ = 3.24 × 105 m2 s−1, and this sets a minimum effective diffusivity in the absence of eddy stirring. As shown in Figs. 8 and 10, the effective eddy diffusivity in our diagnostic can be much smaller than this minimum, and it would be zero without stirring [see Eq. (22)]. Also, when the hyperdiffusion is used rather than the regular diffusion, the effective diffusivity is not restricted by such a diffusivity minimum [see also Nakamura and Zhu (2010)]. Despite these differences, the time-mean effective eddy diffusivity is positive, indicating downgradient eddy tracer fluxes.
e. Sensitivities to numerical schemes of small-scale diffusion
The spectral diffusion is often not desirable for tracer transport because of the nonlocal Gibbs effect. For all the tracer transport experiments, we have repeated the transport calculation with a grid tracer advected by a semi-Lagrangian scheme for horizontal advection, and the grid tracer generally agrees with the spectral tracer advected by the same circulation. To show this, we calculate eddy mixing by the equivalent length [Eq. (23)] for the spectral and grid tracers advected by the same adiabatic flow for ϵ = 20 K. Despite different numerical schemes, Fig. 12 shows a very good agreement in eddy mixing in terms of the tracer equivalent length, with consistent eddy diffusivity maxima and minima between the two schemes. A careful analysis of the diffusivity magnitude for the 330-K isentrope and the near-surface flow indicates that the equivalent length of the spectral tracer is larger than that of the grid tracer (Fig. 13). As the equivalent length quantifies the complexity in tracer geometry, this is consistent with the snapshot of tracer contours that the spectral tracer has slightly more finescale chaotic structures than the grid tracer (not shown). This suggests that the ∇8 hyperdiffusion is slightly less diffusive than the semi-Lagrangian scheme.
Eddy equivalent length ratio (shading) and zonal wind (black contours) averaged for days 100–300, as in Fig. 8. The eddy equivalent length ratio is calculated by Eq. (23) from (left) spectral tracer (expt 4a) and (right) grid tracer (expt 7) advected by the same adiabatic flow. Contour intervals for the zonal winds are 5 m s−1.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
Eddy equivalent length ratio averaged for days 100–300, as in Fig. 12, but for (top) 330 K and (bottom) 95% quantile of surface potential temperature. Solid (dashed) lines denote the spectral (grid) tracer.
Citation: Journal of the Atmospheric Sciences 71, 9; 10.1175/JAS-D-13-0333.1
Using a global tracer variance budget described in Allen and Nakamura (2001), we can parameterize the decay of the global-mean tracer variance by a regular diffusion, and the diffusion coefficient is estimated as κ ≈ 1 × 105 m2 s−1 for both spectral and grid tracers. This leads to an approximate effective isentropic diffusivity as
Shuckburgh and Haynes (2003) and Marshall et al. (2006) found that the effective diffusivity of tracer transport becomes independent of the small-scale diffusion in the circulation regime where the Peclet number (approximately the square of the ratio of eddy length scale and Batchelor scale) exceeds 50. In this circulation regime, the Peclet number scales with the equivalent length ratio of tracer contours, which should also be greater than 50. Figure 13 suggests the current model resolution is not sufficient for this circulation regime. Also, Eq. (22) indicates a weak sensitivity of eddy diffusivity to the diffusion coefficient even when the effective diffusivity becomes independent of the coefficient. Therefore, the exact magnitude of effective eddy diffusivity diagnosed in our study is likely sensitive to the resolution of the model examined, although the pattern is expected to remain unchanged.
5. Conclusions and discussion
A theoretical framework for the effective isentropic diffusivity of tropospheric transport is presented. Compared with previous isentropic analysis (Haynes and Shuckburgh 2000a,b; Allen and Nakamura 2001), we develop a diffusivity diagnostic that allows the isentropes to intersect with the ground frequently. Using the small-scale diffusion experienced by the tracer, the effective diffusivity can account for both the horizontal and vertical parcel displacements in the near-surface isentropic flow. Furthermore, this effective diffusivity is directly linked to a diffusive closure of eddy fluxes through a finite-amplitude wave activity equation.
The role of atmospheric circulation in the tracer transport is investigated in the idealized simulations of a primitive equation model, and the effects of the advecting flow are separated with different components of the flow [i.e., total, TEM, quasigeostrophic (QG), and adiabatic flows]. While the tropical ascent of tracer isopleths and interhemispheric exchange of tracers in the total flow are captured by the simulation advected by the TEM residual circulation, the extratropical tracer isopleth slopes are reproduced only in the simulation advected by the adiabatic flow. The simulation with the QG flow cannot yield either the tropical or extratropical tracer isopleth slopes, but the regions of strong and weak mixing agree with those advected by the adiabatic flow. These individual roles of diabatic circulation and eddy mixing in tropospheric transport are consistent with our understandings of their roles in stratospheric transport (e.g., Holton 1986; Mahlman et al. 1986).
The effective diffusivity is evaluated for the QG flow along the model level and for the adiabatic flow along the isentropic surface, respectively. Although the directions of eddy mixing and the resulting tracer isopleth slopes are different for the QG flow and adiabatic flow, the patterns of effective diffusivities agree well: the upper troposphere is characterized by a diffusivity minimum at the jet’s center with enhanced mixing at the jet’s flanks, and the lower troposphere is dominated by stronger mixing throughout the baroclinic zone. These diffusivity patterns are relatively independent of the initial tracer gradients and are qualitatively similar with altered jet speeds. Furthermore, the isentropic diffusivity agrees broadly with the equivalent lengths of tracer contours by using either a spectral diffusion scheme or a semi-Lagrangian advection scheme. These results provide additional evidence to tropospheric mixing characteristics (Fig. 1) in previous studies with the single isentropic layer model (e.g., Haynes and Shuckburgh 2000b; Allen and Nakamura 2001) or the QG model (e.g., Greenslade and Haynes 2008).
These results above also indicate that our understandings of eddy mixing in the troposphere from the QG theory can be applied, at least qualitatively, to the isentropic flow that intersects the ground. However, this does not mean that the vertical parcel displacement is unimportant. We have mainly considered quasi-horizontal diffusion in our numerical calculation, but in reality the vertical diffusion associated with the vertical parcel displacement can be important (Haynes and Anglade 1997). Leibensperger and Plumb (2014) considered both horizontal and vertical diffusions and found that the resulting effective diffusivities for horizontal and vertical diffusions are qualitatively similar.
It should be noted that the idealized dry model examined has not included the effects of small-scale convection and associated vertical diffusion on tracer transport. It is expected that convection can generate stirring (e.g., Rossby waves) in addition to the baroclinic instability in the dry model. However, the resemblance of our results to the seasonal cycle of effective diffusivity in observations (e.g., Haynes and Shuckburgh 2000b) indicates that our results can hold qualitatively even in the presence of moisture. This is consistent with the results that global-scale tracer transport in the model simulations with and without convective transport are qualitatively similar (Erukhimova and Bowman 2006; Hess 2005).
While previous works on effective diffusivity focus on long-lived tracers, the wave activity budget in Eq. (29) can be readily extended to include nonnegligible sources or sinks of tracers other than the small-scale diffusion. This can be particularly useful to understand the upgradient PV fluxes near the subtropical jet. While our idealized model can simulate the upgradient eddy PV flux near the jet (Fig. 3d), the time-mean effective diffusivity of the adiabatic flow is positive (Fig. 8), implying that the adiabatic eddy tracer fluxes are downgradient. This indicates that diabatic heating may be important for these upgradient fluxes, although a complete wave activity budget of the isentropic PV is needed to resolve this issue. Similarly, one may analyze the wave activity budget of air pollutants to understand the effect of diabatic heating, isentropic mixing, emission, and deposition on their global-scale transport.
The gradient–flux relationships for quasi-conservative tracers have been derived for small-amplitude waves (Plumb 1979; Matsuno 1980; Holton 1981) and for finite-amplitude eddies in the QG flow (Nakamura and Zhu 2010). Our theoretical framework has extended these works to pure horizontal flow in a general vertical coordinate. An improved understanding of the eddy diffusivity is valuable for many aspects of the global atmospheric circulation such as atmospheric energy transport (e.g., Held and Larichev 1996; Lapeyre and Held 2003), the height of extratropical tropopause (e.g., Haynes and Shuckburgh 2000a; Jansen and Ferrari 2013), or the latitude of surface westerly winds (e.g., Chen et al. 2013).
Acknowledgments
We thank Eric Leibensperger, Noboru Nakamura, Peter Hess, and Lantao Sun for valuable discussion. We are also grateful for three anonymous reviewers for constructive comments. GC is supported by the National Science Foundation (NSF) Climate and Large-Scale Dynamics program under Grants AGS-1042787 and AGS-1064079.
APPENDIX
Small-Amplitude Limit of Effective Eddy Diffusivity [Eq. (26)]





Substituting Eqs. (A6) and (A7) into Eq. (A5) yields Eq. (26).
REFERENCES
Allen, D. R., and N. Nakamura, 2001: A seasonal climatology of effective diffusivity in the stratosphere. J. Geophys. Res., 106, 7917–7935, doi:10.1029/2000JD900717.
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. International Geophysics Series, Vol. 40, Academic Press, 489 pp.
Birner, T., D. W. J. Thompson, and T. G. Shepherd, 2013: Up-gradient eddy fluxes of potential vorticity near the subtropical jet. Geophys. Res. Lett., 40, 5988–5993, doi:10.1002/2013GL057728.
Bowman, K. P., and G. D. Carrie, 2002: The mean-meridional transport circulation of the troposphere in an idealized GCM. J. Atmos. Sci., 59, 1502–1514, doi:10.1175/1520-0469(2002)059<1502:TMMTCO>2.0.CO;2.
Bowman, K. P., and T. Erukhimova, 2004: Comparison of global-scale Lagrangian transport properties of the NCEP reanalysis and CCM3. J. Climate, 17, 1135–1146, doi:10.1175/1520-0442(2004)017<1135:COGLTP>2.0.CO;2.
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, doi:10.1175/1520-0469(1986)043<1319:TAOTSP>2.0.CO;2.
Chen, G., 2013: The mean meridional circulation of the atmosphere using the mass above isentropes as the vertical coordinate. J. Atmos. Sci., 70, 2197–2213, doi:10.1175/JAS-D-12-0239.1.
Chen, G., J. Lu, and L. Sun, 2013: Delineating the eddy–zonal flow interaction in the atmospheric circulation response to climate forcing: Uniform SST warming in an idealized aquaplanet model. J. Atmos. Sci., 70, 2214–2233, doi:10.1175/JAS-D-12-0248.1.
Erukhimova, T., and K. P. Bowman, 2006: Role of convection in global-scale transport in the troposphere. J. Geophys. Res., 111, D03105, doi:10.1029/2005JD006006.
Greenslade, M. D., and P. H. Haynes, 2008: Vertical transition in transport and mixing in baroclinic flows. J. Atmos. Sci., 65, 1137–1157, doi:10.1175/2007JAS2236.1.
Haynes, P. H., and J. Anglade, 1997: The vertical-scale cascade in atmospheric tracers due to large-scale differential advection. J. Atmos. Sci., 54, 1121–1136, doi:10.1175/1520-0469(1997)054<1121:TVSCIA>2.0.CO;2.
Haynes, P. H., and E. F. Shuckburgh, 2000a: Effective diffusivity as a diagnostic of atmospheric transport: 1. Stratosphere. J. Geophys. Res., 105, 22 777–22 794, doi:10.1029/2000JD900093.
Haynes, P. H., and E. F. Shuckburgh, 2000b: Effective diffusivity as a diagnostic of atmospheric transport: 2. Troposphere and lower stratosphere. J. Geophys. Res., 105, 22 795–22 810, doi:10.1029/2000JD900092.
Held, I. M., and M. J. Suarez, 1994: A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models. Bull. Amer. Meteor. Soc., 75, 1825–1830, doi:10.1175/1520-0477(1994)075<1825:APFTIO>2.0.CO;2.
Held, I. M., and V. D. Larichev, 1996: A scaling theory for horizontally homogeneous, baroclinically unstable flow on a beta plane. J. Atmos. Sci., 53, 946–952, doi:10.1175/1520-0469(1996)053<0946:ASTFHH>2.0.CO;2.
Hess, P. G., 2005: A comparison of two paradigms: The relative global roles of moist convective versus nonconvective transport. J. Geophys. Res., 110, D20302, doi:10.1029/2004JD005456.
Holton, J. R., 1981: An advective model for two-dimensional transport of stratospheric trace species. J. Geophys. Res., 86, 11 989–11 994, doi:10.1029/JC086iC12p11989.
Holton, J. R., 1986: A dynamically based transport parameterization for one-dimensional photochemical models of the stratosphere. J. Geophys. Res., 91, 2681–2686, doi:10.1029/JD091iD02p02681.
Jansen, M., and R. Ferrari, 2013: The vertical structure of the eddy diffusivity and the equilibration of the extratropical atmosphere. J. Atmos. Sci., 70, 1456–1469, doi:10.1175/JAS-D-12-086.1.
Lapeyre, G., and I. M. Held, 2003: Diffusivity, kinetic energy dissipation, and closure theories for the poleward eddy heat flux. J. Atmos. Sci., 60, 2907–2916, doi:10.1175/1520-0469(2003)060<2907:DKEDAC>2.0.CO;2.
Lee, S., and I. M. Held, 1993: Baroclinic wave packets in models and observations. J. Atmos. Sci., 50, 1413–1428, doi:10.1175/1520-0469(1993)050<1413:BWPIMA>2.0.CO;2.
Leibensperger, E. M., and R. A. Plumb, 2014: Effective diffusivity in baroclinic flow. J. Atmos. Sci., 71, 972–984, doi:10.1175/JAS-D-13-0217.1.
Lin, S.-J., W. C. Chao, Y. C. Sud, and G. K. Walker, 1994: A class of the van Leer-type transport schemes and its application to the moisture transport in a general circulation model. Mon. Wea. Rev., 122, 1575–1593, doi:10.1175/1520-0493(1994)122<1575:ACOTVL>2.0.CO;2.
Mahlman, J. D., H. Levy II, and W. J. Moxim, 1986: Three-dimensional simulations of stratospheric N2O: Predictions for other trace constituents. J. Geophys. Res., 91, 2687–2707, doi:10.1029/JD091iD02p02687.
Marshall, J., E. F. Shuckburgh, H. Jones, and C. Hill, 2006: Estimates and implications of surface eddy diffusivity in the Southern Ocean derived from tracer transport. J. Phys. Oceanogr., 36, 1806–1821, doi:10.1175/JPO2949.1.
Matsuno, T., 1980: Lagrangian motion of air parcels in the stratosphere in the presence of planetary waves. Pure Appl. Geophys. (PAGEOPH), 118, 189–216, doi:10.1007/BF01586451.
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, doi:10.1175/1520-0469(1995)052<2096:MLMDOT>2.0.CO;2.
Nakamura, N., 1996: Two-dimensional mixing, edge formation, and permeability diagnosed in an area coordinate. J. Atmos. Sci., 53, 1524–1537, doi:10.1175/1520-0469(1996)053<1524:TDMEFA>2.0.CO;2.
Nakamura, N., 1998: Leaky containment vessels of air: A Lagrangian-mean approach to the stratospheric tracer transport. Dynamics of Atmospheric Flows: Atmospheric Transport and Diffusion Processes, M. P. Singh and S. Raman, Eds., WIT Press, 193–246.
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, doi:10.1175/2010JAS3432.1.
Nakamura, N., and A. Solomon, 2011: Finite-amplitude wave activity and mean flow adjustments in the atmospheric general circulation. Part II: Analysis in the isentropic coordinate. J. Atmos. Sci., 68, 2783–2799, doi:10.1175/2011JAS3685.1.
Plumb, R. A., 1979: Eddy fluxes of conserved quantities by small-amplitude waves. J. Atmos. Sci., 36, 1699–1704, doi:10.1175/1520-0469(1979)036<1699:EFOCQB>2.0.CO;2.
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, doi:10.1175/1520-0469(1987)044<0298:TZATCO>2.0.CO;2.
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, doi:10.1175/JPO-2669.1.
Randel, W. J., and I. M. Held, 1991: Phase speed spectra of transient eddy fluxes and critical layer absorption. J. Atmos. Sci., 48, 688–697, doi:10.1175/1520-0469(1991)048<0688:PSSOTE>2.0.CO;2.
Shuckburgh, E. F., and P. H. Haynes, 2003: Diagnosing transport and mixing using a tracer-based coordinate system. Phys. Fluids, 15, 3342, doi:10.1063/1.1610471.
Stone, P. H., and D. Miller, 1980: Empirical relations between seasonal changes in meridional temperature gradients and meridional fluxes of heat. J. Atmos. Sci., 37, 1708–1721, doi:10.1175/1520-0469(1980)037<1708:ERBSCI>2.0.CO;2.
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, doi:10.1175/1520-0469(2004)061<2370:EFDBOI>2.0.CO;2.
Tung, K. K., 1986: Nongeostrophic theory of zonally averaged circulation. Part I: Formulation. J. Atmos. Sci., 43, 2600–2618, doi:10.1175/1520-0469(1986)043<2600:NTOZAC>2.0.CO;2.
We have performed direct calculations with the pressure coordinates for the QG flow. The results in pressure coordinates are very similar to those in sigma coordinates, but the pressure coordinates have to deal with missing values near the surface when the ground does not coincide with an isobaric surface.