The angular momentum per unit mass is defined as
The first term represents the advection of angular momentum by the mean meridional circulation, having the form . The uniform rotation component of this, , represents the Coriolis force which redirects meridional flows into zonal flows. Within the anelastic approximation, the divergence of may also be expressed as
The net angular momentum transport through any closed surface of constant must vanish due to the divergenceless nature of the mass flux. For similar reasons, the component of due to the uniform rotation, , cannot transport angular momentum across cylindrical surfaces aligned with the rotation axis. This result also applies to the more general case of a cylindrical rotation profile . Any net transport of angular momentum toward or away from the rotation axis by meridional circulation must come from the advection of the non-cylindrical component of the differential rotation (see also Section 4.3.2).
It may also be noted that angular momentum transport by meridional circulation alone cannot produce localized minima or maxima in . This follows from Equation (6), since vanishes at local extrema. Isolated features in the differential rotation profile such as jets must be produced by other means.
The main driver in maintaining the solar rotation profile is thought to be the Reynolds stress, . This term represents the redistribution of angular momentum by non-axisymmetric motions, particularly convection. Rotation, stratification, magnetic fields, and the spherical shell geometry all introduce anisotropies into the flow which give rise to systematic correlations between the fluctuating velocity components. Horizontal velocity correlations produce latitudinal angular momentum transport whereas correlations produce radial transport. Elucidating the nature of these correlations ranks among the greatest challenges in solar interior dynamics.
In the solar envelope, the Reynolds stress is dominated by turbulent convection, but other motions may also contribute in the tachocline and radiative interior. Convective overshoot excites a spectrum of internal wave modes, most notably gravity waves, which propagate throughout the radiative interior (see Section 8.4). In the absence of dissipation, linear waves cannot redistribute angular momentum. However, dissipation by thermal diffusion or wave breaking can induce a net angular momentum transport via the Reynolds stress which is generally long-range (non-local) and therefore difficult to model. A reliable model of wave transport requires a realistic depiction of wave generation, propagation, and dissipation, which is a formidable task due to the wide range of spatial scales involved. Other potential sources of Reynolds and Maxwell stresses include shear instabilities (see Section 8.2).
Magnetism can alter the rotation profile either by altering the Reynolds stress or by redistributing angular momentum directly via the Lorentz force. The angular momentum flux by the Lorentz force is here decomposed into contributions from fluctuating (non-axisymmetric) fields, , and mean (axisymmetric) fields, . The fluctuating component is known as the Maxwell stress and involves the nonlinear correlations and . Like the Reynolds stress, these may arise from turbulent convection, waves, or instabilities, and understanding their nature is every bit as challenging. The mean-field contribution is more straightforward and can be expressed as
The viscous contribution, , is negligible in the Sun but can be significant in numerical and theoretical models (see Section 6.3). This term opposes angular velocity gradients, , driving the system toward uniform rotation.
The primary angular momentum balance in the Sun is thought to be between the Reynolds stress and meridional circulation, with a lesser role played by the Lorentz force. Thus, if the differential rotation is in a statistically steady state, we expect the following to hold, at least in an approximate and time-averaged sense:
Although the angular velocity in the solar envelope, , varies by from equator to pole and exhibits nearly radial contours at mid-latitudes (Figure 6, panel a), the corresponding specific angular momentum, , is approximately cylindrical (Figure 6, panel b). The hypothetical meridional circulation pattern shown in panel c of Figure 6 would redistribute this angular momentum as shown in panel d of Figure 6. Thus, if the balance expressed in Equation (8) holds, the Reynolds stress must act to accelerate the lower convection zone and equatorial regions and to decelerate the upper convection zone in order to offset the advection of angular momentum by the meridional circulation. Any self-consistent mean-field model which exhibits a solar-like differential rotation profile as shown in panel a of Figure 6 and a single-celled meridional circulation pattern as shown in panel c of Figure 6 must include a Reynolds stress parameterization which redistributes angular momentum as shown in panel d of Figure 6 (unless the Lorentz force plays a significant role).
The results shown in Figure 6 are easily generalized to more complicated circulation patterns. If the angular momentum transport by Reynolds stress is to maintain a balance, it must converge wherever the circulation is away from the rotation axis and diverge wherever it is toward the rotation axis. This is best demonstrated by expressing the meridional circulation flux divergence as in Equation (6) and by noting that is directed away from the rotation axis. Another perspective can be gained by turning the problem around. For a given model of the Reynolds stress, helioseismic rotation profiles can be used to deduce the meridional circulation needed to maintain an equilibrium. This has been done by Durney (2000a).
In the previous section we discussed the mechanisms which can redistribute angular momentum in the solar interior, giving rise to differential rotation. There is more we can say about the angular momentum balance which may eventually be achieved if we consider the limit of rapid rotation3 such that , where the Rossby number is defined as
If the rotation profile satisfies Equation (11) it is said to be in thermal wind balance, in analogy with the thermal wind of geophysical fluid dynamics (Pedlosky, 1987). More specifically, the thermal wind component of the differential rotation is the component which is non-cylindrical and which satisfies Equation (11).
In a thermal wind, departures from cylindrical symmetry are maintained by latitudinal entropy gradients. This is consistent with the angular momentum Equation (5) because if the Taylor-Proudman balance is satisfied perfectly, then both the Reynolds stress and the meridional circulation are negligible (as are the Lorentz and viscous forces), so Equation (5) becomes degenerate. However, meridional circulations are the means by which the thermal wind balance is established and maintained in a rapidly-rotating fluid shell. An imbalance in Equation (11) will drive circulations which will redistribute angular momentum until balance is achieved.
In the solar envelope, latitudinal entropy gradients may be established by the influence of rotation on the efficiency of the convection. For example, if convection is more efficient in the polar regions where the rotation vector is nearly vertical, then these regions will be relatively warm. In radiative equilibrium (the net energy flux into the convection zone equals the net flux out through the surface), such efficiency variations must be balanced by latitudinal energy transport as reflected by Equation (2). Thus, the role played by anisotropic energy transport in maintaining the solar differential rotation may potentially be as important as that played by the Reynolds stress, and may be just as enigmatic.
© Max Planck Society and the author(s)