Turbulence in the lower tachocline will be highly anisotropic due to the strong stable stratification and the large rotational influence. Both effects tend to make the dynamics quasi2D, but in very different ways. The rotational influence will induce vertical coherence, organizing the flow into vortex columns aligned with the rotation vector (e.g., Bartello et al., 1994; Cambon et al., 1997). This is another manifestation of the TaylorProudman theorem which was also discussed in Section 4.3.2. Meanwhile, stable stratification inhibits vertical flows and tends to decouple horizontal layers, favoring pancakelike vortices with large vertical shear (e.g., Métais and Herring, 1989; Riley and Lelong, 2000; GodoyDiana et al., 2004). The relative influence of these two competing effects can be gauged by the Rossby deformation radius, , defined by
where is the BruntVäisälä frequency. The Rossby number and Froude number are defined as and where is a characteristic velocity scale. For motions on scales , stratification breaks the vertical coherence induced by rotation (Dritschel et al., 1999). In the lower tachocline so stratification dominates but approaches zero at the base of the convection zone.Twodimensional turbulence has been studied extensively both theoretically and numerically. It is now well known that nonlinear interactions involving triads of wavevectors in 2D turbulence conserve enstrophy (vorticity squared) as well as energy, and that this gives rise to an inverse cascade of energy from small to large scales (e.g., Lesieur, 1997; Pope, 2000)^{20}. This is in stark contrast to 3D turbulence which exhibits a forward cascade of energy from large to small scales where dissipation occurs. The inverse cascade is manifested as small vortices interact and coalesce into larger vortices. The inverse cascade in 2D turbulence will proceed to the largest scales unless some mechanism suppresses it, such as surface drag in the oceans and atmosphere. Another mechanism for halting the inverse cascade which is more relevant for solar applications occurs in geometries which admit Rossby waves such as rotating spherical shells or planes. If the rotation is rapid enough, patches of vorticity can propagate as Rossby wave packets and disperse before they coalesce. Since the phase speed of a Rossby wave increases with the wavelength (see Appendix A.6), this occurs only for wavenumbers below a critical value , often referred to as the Rhines wavenumber after Rhines (1975). At scales above , the flow has a Rossbywave character and at scales below , it has the character of 2D turbulence.
The most notable thing about the arrest of the inverse cascade by Rossby wave dispersion is that it is anisotropic (Rhines, 1975; Vallis and Maltrud, 1993). Low latitudinal wavenumbers are suppressed, but the cascade can proceed to low longitudinal wavenumbers^{21}. This tends to produce banded zonal flows as observed in the jovian planets (Yoden and Yamada, 1993; Nozawa and Yoden, 1997; Huang and Robinson, 1998; Danilov and Gurarie, 2004). Similar processes also occur in shallowwater and twolayer systems, in both freely decaying and forced configurations (Panetta, 1993; Rhines, 1994; Cho and Polvani, 1996a,b; Peltier and Stuhne, 2002; Kitamura and Matsuda, 2004). The number of bands, or jets, is roughly given by . Taking yields in the solar tachocline, which implies as many as 15 jets. Does a quasi2D inverse cascade occur in real 3D flows? It does in the socalled quasigeostrophic limit first studied by Charney (1971). He showed that in the limit of strong stratification and rapid rotation (), nonlinear interactions conserve potential enstrophy (potential vorticity squared; see Appendix A.6) as well as energy, again giving rise to an inverse cascade of energy (see also Salmon, 1978; Vallis, 2005). This has been demonstrated in 3D simulations by Métais et al. (1996). However, the quasigeostrophic limit does not strictly apply to globalscale motions in spherical shells. It is plausible that similar dynamics occur in spherical systems but this has not yet been rigorously demonstrated.
Thus far in our discussion we have neglected magnetic fields, which can have a profound influence on selforganization processes in turbulent fluids. MHD turbulence does not conserve enstrophy even in the 2D limit, so there is nothing to inhibit a forward cascade of kinetic energy (e.g., Biskamp, 1993; Kim and Dubrulle, 2002). An inverse cascade does occur, but it involves a different ideal invariant, namely magnetic helicity (or in 2D, the magnetic potential). Thus, the physical mechanisms described above which can create banded zonal flows probably do not operate globally in the tachocline, although related dynamics likely occur in relatively fieldfree regions. Selforganization in MHD turbulence generally proceeds by creating largescale magnetic structures which can then feed back on mean flows through the Maxwell stress.
Another important factor in a tachocline context is the presence of rotational shear imposed by largescale stresses from the overlying convective envelope. If the turbulence is itself driven by instabilities of this rotational shear, one may expect it to have a diffusive influence, extracting energy from the shear flow by reducing its amplitude. One may also expect a diffusive behavior if the turbulence is smallscale, isotropic, and homogeneous across horizontal surfaces. In other words, if there is a scale separation with local turbulent mixing. Alternatively, if the flow is dominated by waves, one might expect nonlocal transport which is in general nondiffusive (e.g., McIntyre, 1998, 2003).
The influence of an imposed differential rotation on 2D turbulence in a plane was studied by Shepherd (1987). He found that the shearing of vortices by the differential rotation substantially altered the nonlinear transfer rates among spectral modes. In forceddissipative simulations, smallscale turbulent motions tended to extract energy from the mean shear but the shearinduced Reynolds stress from the largerscale wave field () tended to amplify the mean flow. The net transfer between the mean flow and the fluctuations about it depended sensitively on the parameters of the problem. Shepherd concluded that this complex interaction could not be modeled with a simple linear parameterization, diffusive or otherwise. More recent simulations by Williams (2003) in 2D spherical shells have also shown that the interaction between Rossby wave turbulence and horizontal shear flows can act either to suppress or enhance the shear, depending on the particular details of the problem.
Research into the interaction between a shear flow and 3D, stablystratified turbulence has focused mainly on the case of nonrotating Cartesian domains with vertical shear. Here an important parameter is the Richardson number where is the mean shear (cf. Section 8.2). At small (sheardominated), the turbulent transport of momentum and buoyancy tends to be downgradient (diffusive) but at large (buoyancydominated), turbulent transport is generally oscillatory and can be persistently countergradient (Holt et al., 1992; Galmiche et al., 2002; Jacobitz, 2002). These studies are based on numerical simulations of freelyevolving (decaying) turbulence with homogeneous and isotropic initial conditions and an imposed shear. An effective timedependent viscosity and diffusivity can be defined based on the instantaneous turbulent fluxes and the mean gradients as shown in Figure 26. Countergradient transport is manifested as a negative turbulent viscosity after about 2.5 shear timescales in the stronglystratified simulation (Figure 26, panel b). Although oscillatory, countergradient transport is a robust result of these numerical experiments, it may be a consequence of how they are set up; turbulent fluctuations are sheared by a mean flow which is switched on at some arbitrary time. An analysis in terms of rapid distortion theory by Hanazaki and Hunt (2004) suggests that the countergradient fluxes become very weak as and may be absent altogether in statistically steady flows.

Countergradient transport in stablystratified flows is often associated with the presence of waves (although this is not the only mechanism, (e.g., Holt et al., 1992; Galmiche and Hunt, 2002)). Waves carry pseudomomentum which is conserved until they dissipate, giving rise to longrange transport as described in Section 8.4.
Magnetic fields generally tend to induce downgradient momentum transport in turbulent shear flows by suppressing upscale kinetic energy transfer (cf. inverse cascades) and by imposing rigidity via magnetic tension. However, the transport efficiency can be reduced due to the partial offset of Reynolds and Maxwell stresses, which often have opposite senses (e.g., Kim et al., 2001). Magnetic fields can also suppress turbulent magnetic diffusion (Cattaneo and Vainshtein, 1991; Yousef et al., 2003). Still, magnetism can also have nondiffusive effects. For example, the balance between the Lorentz and Coriolis forces in toroidal field bands can induce zonal jets (see Section 8.2).
In summary, turbulent transport and selforganization in the tachocline is complex and not well understood. A variety of processes can act to establish or to suppress mean flows. Which of these prevail will depend on the subtleties of how the tachocline couples to the convection zone and radiative interior, a topic which will likely occupy researchers for many years to come.
http://www.livingreviews.org/lrsp20051 
© Max Planck Society and the author(s)
Problems/comments to 