[MGI40-P02] Surface zonal flows induced by anelastic thermal convection in a rapidly rotating spherical shell
Keywords:Jovian planets, Anelastic system
Banded structures and alternating zonal jets observed in the surface atmospheres of Jupiter and Saturn have attracted many researchers in planetary atmospheric sciences, however, satisfactory physical explanations and understandings are not yet obtained. In this study, we perform massive parallel numerical experiments treating both small scale convection and planetary scale flows simultaneously, solve fine structures of turbulent motions which have not yet been resolved by the previous numerical models so far, and try to illustrate dynamical origin of global scale structures of surface flows of Jovian planets.
For this purpose, we developed and parallelized an anelastic model of thermal convection in a rotating spherical shell considering basic radial density variation. The spectral transformation library used in this model was improved to introduce MPI parallelization not only in the latitude direction but also in the radial direction. As a result, we succeeded in increasing the number of parallel processes which had been limited by the number of latitudinal grid points, and more massive parallel numerical experiments became possible.
Followup and extended numerical experiments of Case 2 in Gastine al. (2014) were performed with the developed anelastic convection model. When the time integration is as short as viscous diffusion time, the latitudinal profile of induced surface zonal flow is consistent with that of Gastine et al. (2014). However, when time integration is further extended, prograde surface zonal flows in high latitudes start to develop, and become prominent. This result suggests that surface zonal flows obtained by the previous studies may not be statistically steady states.
We started to perform highly resolved numerical simulation of anelastic rotating spherical convection with full spherical domain. The numerical setup follows the case5 of Gastine eta al. (2014), which is one of the highest resolution and most rapid rotation cases performed so far, however, Gastine et al. do not treat full spherical domain but assumed eight-fold symmetry in the longitudinal direction. In the present stage, the broad and strong equatorial prograde jet emerges. The weak but alternating zonal jets are also generated in mid and high latitudes. However, kinetic energy does not steady but still increases, meaning that the system does not reach statistically steady state. It is necessary to continue the time integration and to observe the transition and characteristics of banded structures.
For this purpose, we developed and parallelized an anelastic model of thermal convection in a rotating spherical shell considering basic radial density variation. The spectral transformation library used in this model was improved to introduce MPI parallelization not only in the latitude direction but also in the radial direction. As a result, we succeeded in increasing the number of parallel processes which had been limited by the number of latitudinal grid points, and more massive parallel numerical experiments became possible.
Followup and extended numerical experiments of Case 2 in Gastine al. (2014) were performed with the developed anelastic convection model. When the time integration is as short as viscous diffusion time, the latitudinal profile of induced surface zonal flow is consistent with that of Gastine et al. (2014). However, when time integration is further extended, prograde surface zonal flows in high latitudes start to develop, and become prominent. This result suggests that surface zonal flows obtained by the previous studies may not be statistically steady states.
We started to perform highly resolved numerical simulation of anelastic rotating spherical convection with full spherical domain. The numerical setup follows the case5 of Gastine eta al. (2014), which is one of the highest resolution and most rapid rotation cases performed so far, however, Gastine et al. do not treat full spherical domain but assumed eight-fold symmetry in the longitudinal direction. In the present stage, the broad and strong equatorial prograde jet emerges. The weak but alternating zonal jets are also generated in mid and high latitudes. However, kinetic energy does not steady but still increases, meaning that the system does not reach statistically steady state. It is necessary to continue the time integration and to observe the transition and characteristics of banded structures.