Convective Models of Jupiter’s Zonal Jets with Realistic and Hyper-Energetic Excitation Source

Numerical simulations of Jupiter’s zonal jets are presented, which are generated with realistic and hyper energetic source. The models are three dimensional and nonlinear, applied to a gas that is convective, stratified and compressible. Two solutions are presented, one for a shallow 0.6% envelope, the other one 5% deep. For the shallow model (SM), Jupiter’s small energy flux was applied with low kinematic viscosity. For the deep model (DM), the energy source and viscosity had to be much larger to obtain a solution with manageable computer time. Alternating zonal winds are generated of order 100 m/s, and the models reproduce the observed width of the prograde equatorial jet and adjacent retrograde jets at 20˚ latitude. But the height variations of the zonal winds differ markedly. In SM the velocities vary radially with altitude, but in DM Taylor columns are formed. The dynamical properties of these divergent model results are discussed in light of the computed meridional wind velocities. With large planetary rotation rate Ω, the zonal winds are close to geostrophic, and a quantitative measure of that property is the meridional Rossby number, Ro m . In the meridional momentum balance, the ratio between inertial and Coriolis forces produces Ro m = V 2 /ΩLU, U zonal, V meridional winds, L horizontal length scale. Our analysis shows that the meridional winds vary with the viscosity U/ΩL, of Order 1, which is commonly used as a dimensionless measure of the zonal wind velocities.

In this paper, we present simulations of Jupiter's alternating wind bands generated by convection. The numerical models are nonlinear and three dimensional, applied to a gas that is stratified and compressible. Two solutions are discussed, one for a shallow envelope 0.6% of planetary radius, the other one 5% deep more commensurate with reality. For the shallow model (SM), the small planetary energy flux from the interior was applied together with the corresponding low kinematic viscosity. For the deep model (DM), the applied energy and viscosity had to be much larger to achieve sufficient fast thermal relaxation with manageable computer time, conceptually similar to the deep convective models that have appeared in print.
The mean zonal winds generated with SM and DM differ markedly, varying radially with altitude and aligned along Taylor columns, respectively, and they are discussed in light of the computed meridional winds.

Numerical Models
The numerical models discussed are based on a series of earlier studies (Chan and Sofia [18] [19]; Chan et al. [20]; Chan [21] [22]; Chan and Mayr [15]). Jupiter's atmosphere is treated as an ideal gas, and a uniform energy flux is applied at the bottom. The energy is carried by convection over 95% of the layer, and it is emitted by radiative diffusion from a thin stable layer at the top. The vertical component of the Coriolis force is accounted for, which figures prominently in the momentum budget. A uniform kinematic viscosity is adopted to dissipate the kinetic energy generated by convection, with the Prandtl number set to 1/3.

Shallow Model (SM)
The numerical code employs the transformed spectral procedure with associated vector spherical harmonics and solves the time-dependent 3D nonlinear Navier Stokes equations (Chan et al. [20]; Chan and Mayr [15]). The solution procedure is carried out in two stages.   Figure   1(b) is taken from Chan and Mayr [15] and shows the winds at 4 different altitude levels, identified in the underlying scale of the fractional radius. In qualitative agreement with observations, alternating wind bands are generated. A dominant prograde equatorial jet is produced with velocities close to 70 m/s, which is within a factor of 2 of the observed values. Adjacent to the equatorial jet, strong retrograde winds are generated at 20˚ latitudes, in agreement with the Voyager observations (Smith et al. [2]). But the alternating jets at higher latitudes are too wide and too few in number.
The Taylor-Proudman theorem predicts that for a fluid that is geostrophic and incompressible, the zonal winds become aligned along Taylor columns (TC), and this is the prevailing picture of convective models of the Jupiter atmosphere (e.g., Christensen [12]; Aurnou and Olson [13]; Heimpel and Aurnou [14]). In typical TC models, the tangent cylinder to the inner boundaries determines the width of the equatorial jet, which is defined by the location of the maximum retrograde zonal winds. As illustrated in Figure 1(a), the angular span for the intersections of the tangent cylinder extends to 6˚ latitudes. In contrast, the equatorial jet in Figure 1(b) has a much larger 20˚ span, which demonstrates that SM does not produce a TC pattern at low latitudes. Figure 1 shows that the zonal winds at different levels vary with altitude. This picture is brought into focus in Figure 1 [15] and shows the winds at 4 different altitude levels, identified in the underlying scale of the fractional radius. As illustrated in (a), the angular span for the intersections of the tangent cylinder extends to 6˚ latitudes. In contrast, the equatorial jet in (b) has a much larger 20˚ span, which demonstrates that SM does not produce a Taylor column pattern parallel to the rotation axis. This picture is brought into focus in (c), where the zonal wind pattern of the southern hemisphere is displayed on an expanded altitude scale.
The zonal winds of SM are not aligned along the rotation axis but vary radially with altitude, and the question is whether this property will survive in deeper models with realistic planetary parameters and sufficient low viscosity. Short of the results from such a computationally demanding study, it is instructive to examine a simulation from a numerical model, stratified and compressible, which is much deeper but employs a much larger energy source and viscosity.

Deep Model (DM)
Considering Ohmic dissipation associated with Jupiter's magnetic field and measured conductivity, Liu et al. [23] estimated that the zonal winds cannot penetrate below 0.96 radius. This is on the order of depth employed in present convective models (e.g., Christensen [12]; Aurnou and Olson [13]; Heimpel and Aurnou [14]), referred to as CAOH. CAOH simulate the prograde equatorial jet and alternating wind bands with large viscosity, applying much larger energy flux than Jupiter's.
Following CAOH, a deep model (DM) was constructed with relative depth Δr = 5%, which was presented by Chan and Mayr [16]. For this model, a grid point formulation was employed to integrate the time-dependent 3D nonlinear H. G. Mayr, K. L. Chan International Journal of Astronomy and Astrophysics Navier Stokes equations across the globe (Chan [22]). The model was run with an energy flux, F = 7.8 × 10 4 W/m 2 , a factor of 1.4 × 10 4 larger than that of Jupiter. For the kinematic viscosity the value ν = 2.3 × 10 5 m 2 /s was chosen, which produces the Ekman number E = ν/ΩΔr 2 = 1.0 × 10 −4 . With this large viscosity, the model ran 2 months to produce the numerical results.
Analogous to Figure 1, we present in Figure 2 Figure 1(b) is taken from Chan and Mayr [16] and shows the winds at 4 different altitude levels, identified in the underlying scale of the fractional radius. (a) shows that the angular span of the tangent cylinder extends to 20˚ latitudes, in agreement with the width of the equatorial jet in (b). The zonal winds form a Taylor column pattern, brought into focus in (c). International Journal of Astronomy and Astrophysics

Discussion
For model simulations of a convective Jovian atmosphere that is both stratified and compressible, it is remarkable that the resulting alternating wind bands are formed with such different altitude patterns. In the shallow model the zonal velocities vary radially, but in the deeper model the variations are aligned along the rotation axis to form Taylor columns. Apart from the different vertical domains, the applied energy source and related viscosity must come into play.
The Taylor-Proudman theorem applies if the zonal winds are in geostrophic balance. For the shallow model (SH), geostrophy was explicitly demonstrated by comparing the meridional pressure gradient with the Coriolis force (Chan and Mayr [15]). A quantitative measure of geostrophic balance is the meridional Rossby number, Ro m , which is defined as the ratio between inertial and Coriolis forces. Among the nonlinear inertial accelerations in the meridional momentum equation that describes the mean zonal wind, U, the term V∂V/∂θ is the largest, θ latitude and V mean horizontal meridional wind. Compared with the Coriolis force term ΩU, the meridional Rossby number then can be estimated, Ro m = V 2 /LΩU, with the characteristic horizontal length scale, L = λ/2π (λ, horizontal wavelength) that is related to the planetary radius, r. Figure 3 shows the meridional winds, on the left for the shallow model (SM), on the right for the deep model (DM). At any given latitude, the converging and diverging velocities have opposite directions in the northern and southern hemispheres. Apart from that, the different wind patterns mirror those of the zonal velocities ( Figure 1, Figure 2): a dominant single cell circulation for SM with dominant equatorial jet, and a multi-cellular circulation for the alternating wind bands of DM. Obeying flow continuity, the meridional velocities are much larger  [16]). International Journal of Astronomy and Astrophysics at the top of the domain (black versus green/red), where the ambient densities are much smaller (identical in both models). Apart from the contrasting wind patterns, the meridional winds feature large differences in magnitude. For SM at the top of the domain, the maximum wind velocities are very small less than 0.10 m/s, in contrast to DM with velocities close to 30 m/s. Given the wind velocities, the corresponding model parameters are listed in Table 1 for L = r. The applied viscosities, ν, are sufficiently small in both models to assure that rotation dominates and the Ekman numbers are very small, E << 1. This holds also for the deep model with much larger ν. The meridional Rossby numbers are very small as well, Ro m << 1, demonstrating that both model results are approximately in geostrophic balance. But Ro m is much larger for DM with Taylor columns. Our model results reveal an intriguing relationship between the viscosity and meridional winds. As shown in Table 1, last column, the dimensionless numbers for the two models, V 2 /νΩ differ only by a factor of two, which is remarkable considering that the input parameters for the energy source and viscosity differ by orders of magnitude. The chosen viscosity apparently determines the magnitude of the meridional wind, V varying with (νΩ) 1/2 .
Another intriguing property of the numerical results is the invariance of the zonal velocities in relation to the energy source. As shown in Table 1 Mayr et al. [24] addressed this problem with a simplified 2D scale analysis of the zonal momentum budget, where the Coriolis force, ΩV, is balanced by the viscous stress of the zonal winds, UK/L 2 , with K the eddy viscosity and L the horizontal scale of the circulation. Applying mixing length theory, K = VL, one obtains U = ΩL, which produces for the Jovian circulation zonal winds of order 100 m/s. Essentially, the solution is provided by the nonlinear zonal momentum budget. For the zonal-mean circulation, the inertial force, V∂U/∂λ, dominates, and the balance with the Coriolis force yields, VU/L = ΩV, to produce U = ΩL, independent of the energy source. The same relationship produces the zonal Rossby number, Ro z = U/ΩL, of order 1, which is commonly used as dimensionless measure of the zonal wind velocities. In geostrophic balance, the zonal winds are produced by latitudinal temperature/pressure variations. But the source that produces the temperature variations is also generating the meridional circulation that redistributes or dissipates the kinetic energy. The zonal velocities thus tend to be independent of the energy source.

Summary
Numerical simulations of Jupiter's zonal jets are discussed, which are generated with models that are three dimensional and fully nonlinear, applied to a gas that is convective, stratified and compressible. Solutions are presented for shallow and deep atmospheric envelopes, generated with realistic and hyper-energetic source. In the shallow model (SM) with realistic energy source, the zonal winds vary radially with altitude, in contrast to the energetic deep model (DM) where the winds are aligned along the rotation axis to form Taylor columns (TC). In agreement with observations, both models produce prograde equatorial jets of order 100 m/s. Both models also reproduce the observed width of the equatorial jet with adjacent retrograde jets at 20˚ latitude-a natural outcome for SM, but determined by the chosen 5% depth of DM with TC.
The dynamical properties of these divergent model results are discussed in light of the meridional winds, which are small in magnitude compared with the zonal winds. But unlike the rotational zonal winds, the meridional winds have divergence, and thus are involved with energy and momentum transport, which is of central importance for understanding the zonal mean circulation.
The Rossby number, Ro m , for the meridional momentum balance is the quantitative measure of geostrophy, and it is a quadratic function of the meridional winds. For DM with large viscosity and TC, the meridional winds and Ro m are orders of magnitude larger compared to SM. DM is much less geostrophic. Ranking geostrophy cannot explain the difference between SM and DM. For Taylor columns to form, Taylor-Proudman requires that the gas is also incompressible, in addition to geostrophic. But both models treat the atmosphere as compressible. In models like DM with large viscosity, the enhanced energy transport by the meridional winds has the capacity to reduce the vertical variations in the latitudinal temperature distribution to produce a barotropic environment that favors the formation of TC. The question is whether deeper models, with realistic energy flux and low viscosity, will produce zonal winds that vary radially with altitude like SM.
The numerical results presented highlight an important property of planetary atmospheres, the invariance of the zonal winds in relation to the energy source. With 10 4 times larger source, the velocities of DM increase only by a factor of 3. And the Jupiter models in the literature with energies far exceeding the planetary value all feature zonal winds comparable to those observed. Following up on an earlier paper (Mayr et al. [24]), the solution of this problem is provided by the nonlinear zonal momentum budget with meridional winds. It produces zonal winds varying with the planetary rotation rate and horizontal scale of the circu-International Journal of Astronomy and Astrophysics lation, independent of the energy source.