The asymmetric gravity field measured by the Juno spacecraft allowed estimation of the depth of Jupiter's zonal jets, showing that the winds extend approximately $3000$ km beneath the cloud-level. This estimate was based on an analysis using a combination of all measured odd gravity harmonics $J_{3}$, $J_{5}$, $J_{7}$, and $J_{9}$, but the wind profile dependence on each of them separately has not been investigated. Furthermore, these calculations assumed the meridional profile of the cloud-level wind extends to depth. However, it is possible that the interior jet profile varies from that of the cloud-level as hinted by the Juno microwave measurement that find a smoother nadir brightness temperature profile at depth compared to the cloud-level. Here we analyze in detail the possible meridional and vertical structure of Jupiter's deep jet-streams. We find that each odd gravity harmonic constrains the flow at a different depth, with $J_{3}$ being the most dominant at depths below $3000$ km, $J_{5}$ being the most restrictive overall, and $J_{9}$ not constraining the flow at all if the other odd harmonics are within the measurement range. Interior flow profiles constructed from perturbations to the cloud-level winds allow a more extensive range of vertical wind profiles, yet when the profiles differ substantially from the cloud-level, the ability to match the gravity data reduces significantly. Overall, we find that while interior wind profiles that do not resemble the cloud-level are possible, they are statistically unlikely. However, slightly smoother profiles, which resemble the Juno's microwave radiometer temperature profile at depth, are still compatible with the gravity measurements.

Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

The insolation a planet receives from its parent star is the main engine of the climate and depends on the planet's orbital configuration. Planets with non-zero obliquity and eccentricity will experience seasonal insolation variations. As a result, the climate will have a seasonal cycle, with its strength depending on the orbital configuration, and atmospheric characteristics. In this study, using an idealized general circulation model, we examine the climate response to changes in eccentricity for both zero and non-zero obliquity planets. In the zero obliquity case, a comparison between the seasonal response to changes in eccentricity and perpetual changes in the solar constant shows that the seasonal response strongly depends on the orbital period and radiative timescale. More specifically, using a simple energy balance model, we show the importance of the latitudinal structure of the radiative timescale in the climate response. We also show that the response strongly depends on the atmospheric moisture content. The combination of an eccentric orbit with non-zero obliquity is complex, as the insolation also depends on the perihelion position. Although the detailed response of the climate to variations in eccentricity, obliquity, and perihelion is involved, the circulation is constrained mainly by the thermal Rossby number and latitude of maximum temperature. Finally, we discuss the importance of different planetary parameters that affect the climate response to orbital configuration variations.

The atmospheres of the four giant planets of our Solar System share a common and well-observed characteristic: they each display patterns of planetary banding, with regions of different temperatures, composition, aerosol properties and dynamics separated by strong meridional and vertical gradients in the zonal (i.e., east-west) winds. On Jupiter, the reflective white bands of low temperatures, elevated aerosol opacities, and enhancements of quasi-conserved chemical tracers are referred to as 'zones.' Conversely, the darker bands of warmer temperatures, depleted aerosols, and reductions of chemical tracers are known as `belts.' On Saturn, we define cyclonic belts and anticyclonic zones via their temperature and wind characteristics, although their relation to Saturn's albedo is not as clear as on Jupiter. On distant Uranus and Neptune, the exact relationships between the banded albedo contrasts and the environmental properties is a topic of active study. This review is an attempt to reconcile the observed properties of belts and zones with (i) the meridional overturning inferred from the convergence of eddy angular momentum into the eastward zonal jets at the cloud level on Jupiter and Saturn and the prevalence of moist convective activity in belts; and (ii) the opposing meridional motions inferred from the upper tropospheric temperature structure, which implies decay and dissipation of the zonal jets with altitude above the clouds. These two scenarios suggest meridional circulations in opposing directions, the former suggesting upwelling in belts, the latter suggesting upwelling in zones. This presents an unresolved paradox for our current understanding of the banded structure of giant planet atmospheres, that could be addressed via a multi-tiered vertical structure of 'stacked circulation cells.' [Abridged]

The nature and structure of the observed east-west flows on Jupiter and Saturn has been one of the longest-lasting mysteries in planetary science. This mystery has been recently unraveled due to the accurate gravity measurements provided by the Juno mission to Jupiter and the Grand Finale of the Cassini mission to Saturn. These two experiments, which coincidentally happened around the same time, allowed determination of the vertical and meridional profiles of the zonal flows on both planets. This paper reviews the topic of zonal jets on the gas giants in light of the new data from these two experiments. The gravity measurements not only allow the depth of the jets to be constrained, yielding the inference that the jets extend roughly 3000 and 9000 km below the observed clouds on Jupiter and Saturn, respectively, but also provide insights into the mechanisms controlling these zonal flows. Specifically, for both planets this depth corresponds to the depth where electrical conductivity is within an order of magnitude of 1 S/m, implying that the magnetic field likely plays a key role in damping the zonal flows.

Jupiter's internal flow structure is still not fully known, but can be now better constrained due to Juno's high-precision measurements. The recently published gravity and magnetic field measurements have led to new information regarding the planet and its internal flows, and future magnetic measurements will allow taking another step in resolving this puzzle. In this study, we propose a new method to better constrain Jupiter's internal flow field using the Juno gravity measurements combined with the expected measurements of magnetic secular variation. Based on a combination of hydrodynamical and magnetic field considerations we show that an optimized vertical profile of the zonal flows that fits both measurements can be obtained. Incorporating the magnetic field effects on the flow better constraints the flow decay profile. This will allow getting closer to answering the long-lived question regarding the depth and nature of the flows on Jupiter.

Thousands of exoplanets have been detected to date, and with future planned missions this tally will increase. Understanding the climate dependence on the planetary parameters is vital for the study of terrestrial exoplanet habitability. Using an idealized general circulation model with a seasonal cycle, we study the seasonal response of the surface temperature and Hadley circulation to changes in the orbital, rotational and radiative timescales. Analyzing the climate's seasonal response to variations in these timescales, we find a regime transition between planets controlled by the annual mean insolation to planets controlled by the seasonal variability depending on the relation between the length of the orbital period, obliquity and radiative timescale. Consequently, planets with obliquity greater than $54^{\circ}$ and short orbital period will have a minimum surface temperature at the equator. We also show that in specific configurations, mainly high atmospheric mass and short orbital periods, high obliquity planets can still have an equable climate. Based on the model results, we suggest an empirical power law for the ascending and descending branches of the Hadley circulation and its strength. These power laws show that the Hadley circulation becomes wider and stronger by increasing the obliquity and orbital period or by decreasing the atmospheric mass and rotation rate. Consistent with previous studies, we show that the rotation rate plays an essential role in dictating the width of the Hadley circulation.

Hadley cells dominate the meridional circulation of terrestrial atmospheres. The Solar System terrestrial atmospheres, Venus, Earth, Mars and Titan, exhibit a large variety in the strength, width and seasonality of their Hadley circulation. Despite the Hadley cell being thermally driven, in all planets, the ascending branch does not coincide with the warmest latitude, even in cases with very long seasonality (e.g., Titan) or very small thermal inertia (e.g., Mars). In order to understand the characteristics of the Hadley circulation in case of extreme planetary characteristics, we show both theoretically, using axisymmetric theory, and numerically, using a set of idealized GCM simulations, that the thermal Rossby number dictates the character of the circulation. Given the possible variation of thermal Rossby number parameters, the rotation rate is found to be the most critical factor controlling the circulation characteristics. The results also explain the location of the ascending branch on Mars and Titan.

One of Jupiter's most prominent atmospheric features, the Great Red Spot (GRS), has been observed for more than two centuries, yet little is known about its structure and dynamics below its observed cloud-level. While its anticyclonic vortex appearance suggests it might be a shallow weather-layer feature, the very long time span for which it was observed implies it is likely deeply rooted, otherwise it would have been sheared apart by Jupiter's turbulent atmosphere. Determining the GRS depth will shed light not only on the processes governing the GRS, but on the dynamics of Jupiter's atmosphere as a whole. The Juno mission single flyby over the GRS (PJ7) discovered using microwave radiometer measurements that the GRS is at least a couple hundred kilometers deep (Li et al. 2017). The next flybys over the GRS (PJ18 and PJ21), will allow high-precision gravity measurements that can be used to estimate how deep the GRS winds penetrate below the cloud-level. Here we propose a novel method to determine the depth of the GRS based on the new gravity measurements and a Slepian function approach that enables an effective representation of the wind-induced spatially-confined gravity signal, and an efficient determination of the GRS depth given the limited measurements. We show that with this method the gravity signal of the GRS should be detectable for wind depths deeper than 300 kilometers, with reasonable uncertainties that depend on depth (e.g., $\pm$100km for a GRS depth of 1000km).

How deep do Saturn's zonal winds penetrate below the cloud-level has been a decades-long question, with important implications not only for the atmospheric dynamics, but also for the interior density structure, composition, magnetic field and core mass. The Cassini Grand Finale gravity experiment enables answering this question for the first time, with the premise that the planet's gravity harmonics are affected not only by the rigid body density structure but also by its flow field. Using a wide range of rigid body interior models and an adjoint based thermal wind balance, we calculate the optimal flow structure below the cloud-level and its depth. We find that with a wind profile, largely consistent with the observed winds, when extended to a depth of around 8,800 km, all the gravity harmonics measured by Cassini are explained. This solution is in agreement with considerations of angular momentum conservation, and is consistent with magnetohydrodynamics constraints.

Deciphering the flow below the cloud-level of Jupiter remains a critical milestone in understanding Jupiter's internal structure and dynamics. The expected high-precision Juno measurements of both the gravity field and the magnetic field might help to reach this goal. Here we propose a method that combines both fields to constrain the depth dependent flow field inside Jupiter. This method is based on a mean-field electrodynamic balance that relates the flow field to the anomalous magnetic field, and geostrophic balance that relates the flow field to the anomalous gravity field. We find that the flow field has two distinct regions of influence - an upper region in which the flow affects mostly the gravity field, and a lower region in which the flow affects mostly the magnetic field. An optimization procedure allows to reach a unified flow structure that is consistent with both the gravity and the magnetic fields.

The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction of the planet's radius, is helpful in reconciling the calculated $J_n$ with Juno's observations. Although model predictions are strongly affected by the chosen equation of state, the prediction of an enrichment of $Z$ in the deep, metallic envelope over that in the shallow, molecular envelope holds. We estimate Jupiter's core to contain an 7--25 Earth mass of heavy elements. We discuss the current difficulties in reconciling measured $J_n$ with the equations of state, and with theory for formation and evolution of the planet.

The nature of the flow below the cloud level on Jupiter and Saturn is still unknown. Relating the flow on these planets to perturbations in their density field is key to the analysis of the gravity measurements expected from both the Juno (Jupiter) and Cassini (Saturn) spacecrafts during 2016-17. Both missions will provide latitude-dependent gravity fields, which in principle could be inverted to calculate the vertical structure of the observed cloud-level zonal flow on these planets. Theories to date connecting the gravity field and the flow structure have been limited to potential theories under a barotropic assumption, or estimates based on thermal wind balance that allow analyzing baroclinic wind structures, but have made simplifying assumptions. Those include the effects of the deviations from spherical symmetry, the centrifugal force due to density perturbations, and self-gravitational effects of the density perturbations. Recent studies attempted to include some effects but not in a self-consistent manner. The present study introduces such a self-consistent perturbation approach to the thermal wind balance that incorporates all physical effects, and applies it to several example wind structures, both barotropic and baroclinic. The contribution of each term is analyzed, and the results are compared in the barotropic limit to those of potential theory. It is found that the dominant balance involves the original simplified thermal wind approach. This balance produces a good order-of-magnitude estimate of the gravitational moments, and is able, therefore, to address the order one question of how deep the flows are given measurements of gravitational moments. The additional terms are significantly smaller and none of these terms is dominant, so any approximation attempting to improve over the simplified thermal wind approach needs to include all other terms.

Observations of the flow on Jupiter exists essentially only for the cloud-level, which is dominated by strong east-west jet-streams. These have been suggested to result from dynamics in a superficial thin weather-layer, or alternatively be a manifestation of deep interior cylindrical flows. However, it is possible that the observed winds are indeed superficial, yet there exists deep flow that is completely decoupled from it. To date, all models linking the wind, via the induced density anomalies, to the gravity field, to be measured by Juno, consider only flow that is a projection of the observed could-level wind. Here we explore the possibility of complex wind dynamics that include both the shallow weather-layer wind, and a deep flow that is decoupled from the flow above it. The upper flow is based on the observed cloud-level flow and is set to decay with depth. The deep flow is constructed to produce cylindrical structures with variable width and magnitude, thus allowing for a wide range of possible scenarios for the unknown deep flow. The combined flow is then related to the density anomalies and gravitational moments via a dynamical model. An adjoint inverse model is used for optimizing the parameters controlling the setup of the deep and surface-bound flows, so that these flows can be reconstructed given a gravity field. We show that the model can be used for examination of various scenarios, including cases in which the deep flow is dominating over the surface wind. We discuss extensively the uncertainties associated with the model solution. The flexibility of the adjoint method allows for a wide range of dynamical setups, so that when new observations and physical understanding will arise, these constraints could be easily implemented and used to better decipher Jupiter flow dynamics.

During 2016-17 the Juno and Cassini spacecraft will both perform close eccentric orbits of Jupiter and Saturn, respectively, obtaining high-precision gravity measurements for these planets. This data will be used to estimate the depth of the observed surface flows on these planets. All models to date, relating the winds to the gravity field, have been in the forward direction, thus allowing only calculation of the gravity field from given wind models. However, there is a need to do the inverse problem since the new observations will be of the gravity field. Here, an inverse dynamical model, is developed to relate the expected measurable gravity field, to perturbations of the density and wind fields, and therefore to the observed cloud-level winds. In order to invert the gravity field into the 3D circulation, an adjoint model is constructed for the dynamical model, thus allowing backward integration. This tool is used for examination of various scenarios, simulating cases in which the depth of the wind depends on latitude. We show that it is possible to use the gravity measurements to derive the depth of the winds, both on Jupiter and Saturn, taking into account also measurement errors. Calculating the solution uncertainties, we show that the wind depth can be determined more precisely in the low-to-midlatitudes. In addition, the gravitational moments are found to be particularly sensitive to flows at the equatorial intermediate depths, therefore we expect that if deep winds exist on these planets they will have a measurable signature by Juno and Cassini.

The recent discoveries of terrestrial exoplanets and super-Earths extending over a broad range of orbital and physical parameters suggest that these planets will span a wide range of climatic regimes. Characterization of the atmospheres of warm super-Earths has already begun and will be extended to smaller and more distant planets over the coming decade. The habitability of these worlds may be strongly affected by their three-dimensional atmospheric circulation regimes, since the global climate feedbacks that control the inner and outer edges of the habitable zone including transitions to Snowball-like states and runaway-greenhouse feedbacks depend on the equator-to-pole temperature differences, patterns of relative humidity, and other aspects of the dynamics. Here, using an idealized moist atmospheric general circulation model including a hydrological cycle, we study the dynamical principles governing the atmospheric dynamics on such planets. We show how the planetary rotation rate, stellar flux, atmospheric mass, surface gravity, optical thickness, and planetary radius affect the atmospheric circulation and temperature distribution on such planets. Our simulations demonstrate that equator-to-pole temperature differences, meridional heat transport rates, structure and strength of the winds, and the hydrological cycle vary strongly with these parameters, implying that the sensitivity of the planet to global climate feedbacks will depend significantly on the atmospheric circulation. We elucidate the possible climatic regimes and diagnose the mechanisms controlling the formation of atmospheric jet streams, Hadley and Ferrel cells, and latitudinal temperature differences. Finally, we discuss the implications for understanding how the atmospheric circulation influences the global climate.

The alignment of Saturn's magnetic pole with its rotation axis precludes the use of magnetic field measurements to determine its rotation period. The period was previously determined from radio measurements by the Voyager spacecraft to be 10h 39m 22.4s. When the Cassini spacecraft measured a period of 10h 47m 6s, which was additionally found to change between sequential measurements, it became clear that the radio period could not be used to determine the bulk planetary rotation period. Estimates based upon Saturn's measured wind fields have increased the uncertainty even more, giving numbers smaller than the Voyager rotation period, and at present Saturn's rotation period is thought to be between 10h 32m and 10h 47m, which is unsatisfactory for such a fundamental property. Here we report a period of 10h 32m 45s +- 46s, based upon an optimization approach using Saturn's measured gravitational field and limits on the observed shape and possible internal density profiles. Moreover, even when solely using the constraints from its gravitational field, the rotation period can be inferred with a precision of several minutes. To validate our method, we applied the same procedure to Jupiter and correctly recovered its well-known rotation period.

A variety of observations provide evidence for vigorous motion in the atmospheres of brown dwarfs and directly imaged giant planets. Motivated by these observations, we examine the dynamical regime of the circulation in the atmospheres and interiors of these objects. Brown dwarfs rotate rapidly, and for plausible wind speeds, the flow at large scales will be rotationally dominated. We present 3D, global, numerical simulations of convection in the interior, which demonstrate that, at large scales, the convection aligns in the direction parallel to the rotation axis. Convection occurs more efficiently at high latitudes than low latitudes, leading to systematic equator-to-pole temperature differences that may reach ~1 K near the top of the convection zone. The interaction of convection with the overlying, stably stratified atmosphere will generate a wealth of atmospheric waves, and we argue that, as in the stratospheres of planets in the solar system, the interaction of these waves with the mean flow will cause a significant atmospheric circulation at regional to global scales. At large scales, this should consist of stratified turbulence (possibly organizing into coherent structures such as vortices and jets) and an accompanying overturning circulation. We present an approximate analytic theory of this circulation, which predicts characteristic horizontal temperature variations of several to ~50 K, horizontal wind speeds of ~10-300 m/sec, and vertical velocities that advect air over a scale height in ~10^5-10^6 sec. This vertical mixing may help to explain the chemical disequilibrium observed on some brown dwarfs. Moreover, the implied large-scale organization of temperature perturbations and vertical velocities suggests that, near the L/T transition, patchy clouds can form near the photosphere, helping to explain recent observations of brown-dwarf variability in the near-IR.

The investigation of planets around other stars began with the study of gas giants, but is now extending to the discovery and characterization of super-Earths and terrestrial planets. Motivated by this observational tide, we survey the basic dynamical principles governing the atmospheric circulation of terrestrial exoplanets, and discuss the interaction of their circulation with the hydrological cycle and global-scale climate feedbacks. Terrestrial exoplanets occupy a wide range of physical and dynamical conditions, only a small fraction of which have yet been explored in detail. Our approach is to lay out the fundamental dynamical principles governing the atmospheric circulation on terrestrial planets--broadly defined--and show how they can provide a foundation for understanding the atmospheric behavior of these worlds. We first survey basic atmospheric dynamics, including the role of geostrophy, baroclinic instabilities, and jets in the strongly rotating regime (the "extratropics") and the role of the Hadley circulation, wave adjustment of the thermal structure, and the tendency toward equatorial superrotation in the slowly rotating regime (the "tropics"). We then survey key elements of the hydrological cycle, including the factors that control precipitation, humidity, and cloudiness. Next, we summarize key mechanisms by which the circulation affects the global-mean climate, and hence planetary habitability. In particular, we discuss the runaway greenhouse, transitions to snowball states, atmospheric collapse, and the links between atmospheric circulation and CO2 weathering rates. We finish by summarizing the key questions and challenges for this emerging field in the future.

Three-dimensional studies of convection in deep spherical shells have been used to test the hypothesis that the strong jet streams on Jupiter, Saturn, Uranus, and Neptune result from convection throughout the molecular envelopes. Due to computational limitations, these simulations must adopt viscosities and heat fluxes many orders of magnitude larger than the planetary values. Several numerical investigations have identified trends for how the mean jet speed varies with heat flux and viscosity, but no previous theories have been advanced to explain these trends. Here, we show using simple arguments that if convective release of potential energy pumps the jets and viscosity damps them, the mean jet speeds split into two regimes. When the convection is weakly nonlinear, the equilibrated jet speeds should scale approximately with F/nu, where F is the convective heat flux and nu is the viscosity. When the convection is strongly nonlinear, the jet speeds are faster and should scale approximately as (F/nu)^{1/2}. We demonstrate how this regime shift can naturally result from a shift in the behavior of the jet-pumping efficiency with heat flux and viscosity. Moreover, the simulations hint at a third regime where, at sufficiently small viscosities, the jet speed becomes independent of the viscosity. We show based on mixing-length estimates that if such a regime exists, mean jet speeds should scale as heat flux to the 1/4 power. Our scalings provide a good match to the mean jet speeds obtained in previous Boussinesq and anelastic, three-dimensional simulations of convection within giant planets over a broad range of parameters. When extrapolated to the real heat fluxes, these scalings suggest that the mass-weighted jet speeds in the molecular envelopes of the giant planets are much weaker--by an order of magnitude or more--than the speeds measured at cloud level.