Publications

The equatorial jets dominating the dynamics of the Jovian planets exhibit two distinct types of zonal flows: strongly eastward in the gas giants (superrotation) and strongly westward in the ice giants (subrotation). Existing theories propose different mechanisms for these patterns, but no single mechanism has successfully explained both. However, the planetary parameters of the four Solar System giant planets suggest that a fundamentally different mechanism is unlikely. In this study, we show that convection-driven columnar structures can account for both eastward and westward equatorial jets, framing the phenomenon as a bifurcation. Consequently, both superrotation and subrotation emerge as stable branches of the same mechanistic solution. Our analysis of these solutions uncovers similarities in the properties of equatorial waves and the leading-order momentum balance. This study suggests that the fundamental dynamics governing equatorial jet formation may be more broadly applicable across the Jovian planets than previously believed, offering a unified explanation for their two distinct zonal wind patterns. A bifurcation in equatorial jet formation can explain the differences between the gas and ice giants.

This paper presents an analysis of Juno's first radio occultation experiments. Relying on two-way radio links in the X- and Ka-bands, we processed data from NASA's Deep Space Network antennas through a ray-tracing inversion algorithm. By effectively isolating dispersive effects, we obtained measurements of the neutral atmosphere's characteristics. This enabled the derivation of pressure and temperature profiles from the recorded frequencies. These results complement prior data from Voyager occultations and CIRS observations, providing valuable contributions to our understanding of Jupiter's atmospheric dynamics.

Jupiter, the fastest-rotating planet in the Solar System, exhibits a pronounced equatorial bulge, with its equatorial radius exceeding the polar radius by approximately 7%. This oblate shape reflects the combined effects of rapid rotation, complex internal structure, and atmospheric winds. Existing estimates of Jupiter’s shape, with uncertainties of about 4 km, are based on a single analysis of Voyager and Pioneer radio occultations from nearly five decades ago and do not account for the influence of Jupiter’s strong winds. The Juno spacecraft has recently returned numerous high-precision radio-occultation measurements, enabling a more accurate determination. Incorporating the effects of zonal winds, we derive Jupiter’s shape with an order-of-magnitude reduction in uncertainty. At the 1-bar pressure level, we find a polar radius of 66,842 p/m 0.4 km, an equatorial radius of 71,488 p/m 0.4 km, and a mean radius of 69,886 p/m 0.4 km, which are 12 km, 4 km, and 8 km smaller than previous estimates, respectively. The results indicate that winds above the visible cloud tops are largely barotropic, showing minimal vertical variation. The updated shape has important implications for interior structure models, supporting a metal-enriched and cooler atmosphere, thereby helping reconcile discrepancies between models, Galileo probe measurements, and Voyager-derived temperatures. The refined radius profile also improves spatial referencing for pressure-dependent measurements, offering a more precise context for interpreting Jupiter’s atmospheric dynamics.

The eddy-driven jet meandering has been hypothesized to increase due to climate change. This meandering frequently induces slow-moving patterns of surface low and high-pressure anomalies, potentially causing extreme weather events such as droughts, flooding, heat waves and cold spells. However, the quantitative link between jet meandering and storm development is still lacking, as well as a conclusive mechanism for the effect of climate change on jet meandering. In this study, we separate the physical components in the atmospheric complex system using an idealized moist global circulation model (GCM), performing a series of global warming simulations. We outline the connection between the decreasing equator-to-pole temperature gradient due to Arctic amplification and the meandering of the jet. We derive a theoretical relation between the 500 hPa geopotential height meridional gradient and jet meandering, indicating that the leading reason for the increase in jet meandering is the flattening of the mid-atmospheric meridional temperature gradient, constraining the variations of the streamfunction anomaly with latitude. As it decreases, the eddy-driven jet slows and its meridional layout widens. By Lagrangian tracking of cyclones and anticyclones, we link the jet meandering to the development of cyclones and anticyclones, showing that intense storms are more likely to be associated with jet meanders during their development stage than weak storms. Furthermore, jet-stream meanders with greater meridional extent will be more likely associated with intense storms. In simulations with increased mean surface temperature, our results demonstrate an increase in both jet meandering and storm intensity.

Context. Jupiter’s polar upper troposphere and stratosphere host a persistent cold vortex poleward of 65°N, but its detailed structure and dynamics have remained difficult to resolve.

Aims. The goal is to characterize the thermal structure and dynamics of the polar vortex using new and complementary remote sensing techniques.

Methods. We use the combination of high-resolution vertical profiles derived from Juno’s recent radio occultation measurements with mid-infrared imaging from the VLT/VISIR instrument. The former provided direct retrievals of temperature and density near and within the vortex, while VISIR imaging reveals spatial thermal contrasts across the region.

Results. Our analysis confirms the presence of a steep meridional temperature jump at 65°N, of about 7 K at 100 mbar, consistent with strong vertical wind shear and a prograde polar stratospheric jet reaching up to 80 m s−1 at the 10 mbar level. We find the atmosphere to be thermally stable above 0.55 bar reaching a Brunt–Väisälä frequency of 0.025 in the mid-stratosphere. Thermal contrasts observed in the infrared data align with the vertical structures inferred from radio occultations, validating the presence and extent of the cold vortex.

Conclusions. These findings offer a quantitative analysis of the thermal structure and the dynamical behavior of Jupiter’s polar atmosphere and demonstrate the diagnostic power of combining radio occultation and thermal infrared techniques in planetary atmospheric studies.

Europa, Jupiters second Galilean moon, is believed to host a subsurface ocean in contact with a rocky mantle, where hydrothermal activity may drive the synthesis of organic molecules. Among these possible organic molecules, abiotic synthesis of aromatic amino acids is unlikely, so their detection on planetary surfaces such as Europa suggests that they could be considered a potential biosignature. Fluorescence from aromatic amino acids, with characteristic emissions in the 200-400 nm wavelength range, can be induced by a laser and may be detectable where ocean material has been relatively recently emplaced on Europas surface, as indicated by geologically young terrain and surface features. However, surface bombardment by charged particles from the jovian magnetosphere and solar ultraviolet (UV) radiation degrades organic molecules and limits their longevity. We model radiolysis and photolysis of aromatic amino acids embedded in ice. Our model shows dependencies on hemispheric and latitudinal patterns of charged particle bombardment and ice phase. We demonstrate that such molecules contained within freshly deposited ice in high-latitude regions on the surface of Europa are detectable using laser-induced UV fluorescence, even from an orbiting spacecraft.

The Jupiter Icy Moons Explorer (JUICE) is a European Space Agency mission to explore Jupiter and its three icy Galilean moons: Europa, Ganymede, and Callisto. Numerous JUICE investigations concern the magnetised space environments containing low-density populations of charged particles that surround each of these bodies. In the case of both Jupiter and Ganymede, the magnetic field generated internally produces a surrounding volume of space known as a magnetosphere. All these regions are natural laboratories where we can test and further our understanding of how such systems work, and improved knowledge of the environments around the moons of interest is important for probing sub-surface oceans that may be habitable. Here we review the magnetosphere and plasma science that will be enabled by JUICE from arrival at Jupiter in July 2031. We focus on the specific topics where the mission will push forward the boundaries of our understanding through a combination of the spacecraft trajectory through the system and the measurements that will be made by its suite of scientific instruments. Advances during the initial orbits around Jupiter will include construction of a comprehensive picture of the poorly understood region of Jupiters magnetosphere where rigid plasma rotation with the planet breaks down, and new perspectives on how Jupiters magnetosphere interacts with both Europa and Callisto. The later orbits around Ganymede will dramatically improve knowledge of this moons smaller magnetosphere embedded within the larger magnetosphere of Jupiter. We conclude by outlining the high-level operational strategy that will support this broad science return.

Chang (2024, https://doi.org/10.1029/2024gl110011) challenged the methodology proposed recently by Okajima et al. for evaluating cyclonic and anticyclonic contributions to Eulerian eddy statistics and atmospheric energetics based on the local flow curvature. He argued that using the local wind curvature to separate energetic contributions from cyclones and anticyclones is not physically meaningful. Here we argue that his claims are based on an unrealistic assumption of monopolar relative vorticity in an entire storm-track domain and a meridionally uniform zonal background flow atypical to midlatitudes. We also demonstrate that the error in attributing eddy statistics to cyclones and anticyclones is significantly smaller than his estimation. Rather, we further demonstrate that the curvature-based methodology effectively eliminates the shear influence to identify cyclonic and anticyclonic regions, which is dismissed in his argument. We conclude that the curvature-based methodology is beneficial in evaluating distinct cyclonic and anticyclonic contributions to atmospheric energetics in realistic conditions.

Jupiter's equatorial eastward zonal flows reach wind velocities of ∼100 m s⁻¹, while on Saturn they are three times as strong and extend about twice as wide in latitude, despite the two planets being overall dynamically similar. Recent gravity measurements obtained by the Juno and Cassini spacecraft uncovered that the depth of zonal flows on Saturn is about three times greater than on Jupiter. Here we show, using 3D deep convection simulations, that the atmospheric depth is the determining factor controlling both the strength and latitudinal extent of the equatorial zonal flows, consistent with the measurements for both planets. We show that the atmospheric depth is proportional to the convectively driven eddy momentum flux, which controls the strength of the zonal flows. These results provide a mechanistic explanation for the observed differences in the equatorial regions of Jupiter and Saturn, and offer new understandings about the dynamics of gas giants beyond the Solar System.

Climate models generally predict a poleward shift of the midlatitude circulation in response to climate change induced by increased greenhouse gas concentration, but the intermodel spread of the eddy-driven jet shift is large and poorly understood. Recent studies point to the significance of midlatitude midtropospheric diabatic heating for the intermodel spread in the jet latitude. To examine the role of diabatic heating in the jet response to climate change, a series of simulations are performed using an idealized aquaplanet model. It is found that both increased CO₂ concentration and increased saturation vapor pressure induce a similar warming response, leading to a poleward and upward shift of the midlatitude circulation. An exception to this poleward shift is found for a certain range of temperatures, where the eddy-driven jet shifts equatorward, while the latitude of the eddy heat flux remains essentially unchanged. This equatorward jet shift is explained by the connection between the zonal-mean momentum and heat budgets: increased diabatic heating in the midlatitude midtroposphere balances the cooling by the Ferrel cell ascending branch, enabling an equatorward shift of the Ferrel cell streamfunction and eddy-driven jet, while the latitude of the eddy heat flux remains unchanged. The equatorward jet shift and the strengthening of the midlatitude diabatic heating are found to be sensitive to the model resolution. The implications of these results for a potential reduction in the jet shift uncertainty through the improvement of convective parameterizations are discussed.

Letter

ESAs Jupiter Icy Moons Explorer (JUICE) will provide a detailed investigation of the Jovian system in the 2030s, combining a suite of state-of-the-art instruments with an orbital tour tailored to maximise observing opportunities. We review the Jupiter science enabled by the JUICE mission, building on the legacy of discoveries from the Galileo, Cassini, and Juno missions, alongside ground- and space-based observatories. We focus on remote sensing of the climate, meteorology, and chemistry of the atmosphere and auroras from the cloud-forming weather layer, through the upper troposphere, into the stratosphere and ionosphere. The Jupiter orbital tour provides a wealth of opportunities for atmospheric and auroral science: global perspectives with its near-equatorial and inclined phases, sampling all phase angles from dayside to nightside, and investigating phenomena evolving on timescales from minutes to months. The remote sensing payload spans far-UV spectroscopy (50-210 nm), visible imaging (340-1080 nm), visible/near-infrared spectroscopy (0.49-5.56 μm), and sub-millimetre sounding (near 530-625 GHz and 1067-1275 GHz). This is coupled to radio, stellar, and solar occultation opportunities to explore the atmosphere at high vertical resolution; and radio and plasma wave measurements of electric discharges in the Jovian atmosphere and auroras. Cross-disciplinary scientific investigations enable JUICE to explore coupling processes in giant planet atmospheres, to show how the atmosphere is connected to (i) the deep circulation and composition of the hydrogen-dominated interior; and (ii) to the currents and charged particle environments of the external magnetosphere. JUICE will provide a comprehensive characterisation of the atmosphere and auroras of this archetypal giant planet.

The shape of the two gas giants, Jupiter and Saturn, is determined primarily by their rotation rate, and interior density distribution. It is also affected by their zonal winds, causing an anomaly of O(10 km) at low latitudes. However, uncertainties in the observed cloud-level wind and the polar radius, translate to an uncertainty in the shape with the same order of magnitude. The Juno (Jupiter) and Cassini (Saturn) missions gave unprecedented accurate gravity measurements, constraining better the uncertainty in the wind structure. Using an accurate shape calculation, and a joint optimization, given both gravity and radio-occultation measurements, we calculate the possible range of dynamical height for both planets. We find that for Saturn there is an excellent match to the radio-occultation measurements, while at Jupiter such a match is not achieved. This may point to deviations from a barotropic flow above the cloud level, which might be tested with the forthcoming radio-occultation measurements by Juno.

Clouds are one of the most influential components of Earths climate system. Specifically, the midlatitude clouds play a vital role in shaping Earths albedo. This study investigates the connection between baroclinic activity, which dominates the midlatitude climate, and cloud-albedo and how it relates to Earths existing hemispheric albedo symmetry. We show that baroclinic activity and cloud-albedo are highly correlated. By using Lagrangian tracking of cyclones and anticyclones and analyzing their individual cloud properties at different vertical levels, we explain why their cloud-albedo increases monotonically with intensity. We find that while for anticyclones, the relation between strength and cloudiness is mostly linear, for cyclones, in which clouds are more prevalent, the relation saturates with strength. Using the cloud-albedo strength relationships and the climatology of baroclinic activity, we demonstrate that the observed hemispheric difference in cloud-albedo is well explained by the difference in the population of cyclones and anticyclones, which counter-balances the difference in clear-sky albedo. Finally, we discuss the robustness of the hemispheric albedo symmetry in the future climate. Seemingly, the symmetry should break, as the northern hemispheres storm track response differs from that of the southern hemisphere due to Arctic amplification. However, we show that the saturation of the cloud response to storm intensity implies that the increase in the skewness of the southern hemisphere storm distribution toward strong storms will decrease future cloud-albedo in the southern hemisphere. This complex response explains how albedo symmetry might persist even with the predicted asymmetric hemispheric change in baroclinicity under climate change.

This chapter reviews the way the six key questions about planetary systems, from their origins to the way they work and their habitability, identified in Chapter 1 (Blanc et al., 2021), can be addressed by means of solar system exploration, and how one can find partial answers to these six questions by flying to the different provinces to the solar system: terrestrial planets, giant planets, small bodies, and up to its interface with the local interstellar medium. It derives from this analysis a synthetic description of the most important space observations to be performed at the different solar system objects by future planetary exploration missions. These \u201cobservation requirements\u201d illustrate the diversity of measurement techniques to be used as well as the diversity of destinations where these observations must be made. They constitute the base for the identification of the future planetary missions we need to fly by 2061, which are described in Chapter 4.

Context. The Juno mission has provided measurements of Jupiter's gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. Aims. In this study, we investigated the differences between the state-of-the-art equations of state and their impact on the properties of interior models. Accounting for uncertainty on the hydrogen and helium equation of state, we assessed the span of the interior features of Jupiter. Methods. We carried out an extensive exploration of the parameter space and studied a wide range of interior models using Markov chain Monte Carlo simulations. To consider the uncertainty on the equation of state, we allowed for modifications of the equation of state in our calculations. Results. Our models harbour a dilute core and indicate that Jupiter's internal entropy is higher than what is usually assumed from the Galileo probe measurements. We obtain solutions with extended dilute cores, but contrary to other recent interior models of Jupiter, we also obtain models with small dilute cores. The dilute cores in such solutions extend to ∼20% of Jupiter's mass, leading to better agreement with formation-evolution models. Conclusions. We conclude that the equations of state used in Jupiter models have a crucial effect on the inferred structure and composition. Further explorations of the behaviour of hydrogen-helium mixtures at the pressure and temperature conditions in Jupiter will help to constrain the interior of the planet, and therefore its origin.

The large-scale Hadley circulation is a key element in the global heat and moisture transport. It is traditionally defined as the zonally averaged meridional circulation in the tropics, but was shown to have a strong longitudinal dependence, as seen in a decomposition of the three-dimensional atmospheric flow into spatially dependent meridional and zonal circulations. Recent studies provided a useful analysis of the regional strengthening/weakening of the decomposed circulation but not its patterns. Here, we study the interannual variability of the longitudinally dependent meridional circulation (LMC), with a focus on its spatial patterns. We use hierarchical clustering to objectively determine the four main modes of the LMC interannual variability, and apply a Lagrangian air parcel tracking method to reveal the full circulation patterns. While El Niño and La Niña are found, as in previous studies, to play a role in setting these patterns, we find the patterns are not uniquely characterized by standard El Niño-Southern Oscillation (ENSO) indices (Nino3.4 or Southern Oscillation Index). Instead, ENSO flavors (i.e., East Pacific vs. Central Pacific) have different effects on the LMC. The most prominent interannual variability of the LMC is an east-west shift. Latitudinal shifts, as well as contraction/expansion in both latitude and longitude are also identified. Multiple linear regression analysis shows that while a large fraction of the LMC variance is explained by Sea Surface Temperature, the Madden-Julian Oscillation makes a nonnegligible independent contribution. The clustering patterns are also used to study the remote precipitation and surface air temperature teleconnections.

Of all the myriad environments in our Solar System, the least explored are the distant Ice Giants Uranus and Neptune, and their diverse satellite and ring systems. These intermediate-sized worlds are the last remaining class of Solar System planet to be characterised by a dedicated robotic mission, and may shape the paradigm for the most common outcome of planetary formation throughout our galaxy. In response to the 2019 European Space Agency call for scientific themes in the 2030s and 2040s (known as Voyage 2050), we advocated that an international partnership mission to explore an Ice Giant should be a cornerstone of ESAs science planning in the coming decade, targeting launch opportunities in the early 2030s. This article summarises the inter-disciplinary science opportunities presented in that White Paper [1], and briefly describes developments since 2019.

Clouds are primary modulators of Earth's energy balance. It is thus important to understand the links connecting variabilities in cloudiness to variabilities in other state variables of the climate system, and also describe how these links would change in a changing climate. A conceptual model of global cloudiness can help elucidate these points. In this work we derive simple representations of cloudiness, that can be useful in creating a theory of global cloudiness. These representations illustrate how both spatial and temporal variability of cloudiness can be expressed in terms of basic state variables. Specifically, cloud albedo is captured by a nonlinear combination of pressure velocity and a measure of the low-level stability, and cloud longwave effect is captured by surface temperature, pressure velocity, and standard deviation of pressure velocity. We conclude with a short discussion on the usefulness of this work in the context of global warming response studies.

Context. While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-element distribution in Jupiter's envelope. Aims. In this paper, we assemble the most comprehensive and diverse collection of Jupiter interior models to date and use it to study the distribution of heavy elements in the planet's envelope. Methods. We apply a Bayesian statistical approach to our interior model calculations, reproducing the Juno gravitational and atmospheric measurements and constraints from the deep zonal flows. Results. Our results show that the gravity constraints lead to a deep entropy of Jupiter corresponding to a 1 bar temperature that is 515 K higher than traditionally assumed. We also find that uncertainties in the equation of state are crucial when determining the amount of heavy elements in Jupiter's interior. Our models put an upper limit to the inner compact core of Jupiter of 7 MEarth, independently of the structure model (with or without a dilute core) and the equation of state considered. Furthermore, we robustly demonstrate that Jupiter's envelope is inhomogeneous, with a heavy-element enrichment in the interior relative to the outer envelope. This implies that heavy-element enrichment continued through the gas accretion phase, with important implications for the formation of giant planets in our Solar System and beyond.

Storm-track activity over the North Pacific (NP) climatologically exhibits a clear minimum in midwinter, when the westerly jet speed sharply maximizes. This counterintuitive phenomenon, referred to as the “midwinter minimum (MWM),” has been investigated from various perspectives, but the mechanisms are still to be unrevealed. Toward better understanding of this phenomenon, the present study delineates the detailed seasonal evolution of climatological-mean Eulerian statistics and energetics of migratory eddies along the NP storm track over 60 years. As a comprehensive investigation of the mechanisms for the MWM, this study has revealed that the net eddy conversion/generation rate normalized by the eddy total energy, which is independent of eddy amplitude, is indeed reduced in midwinter. The reduction from early winter occurs mainly due to the decreased effectiveness of the baroclinic energy conversion through seasonally weakened temperature fluctuations and the resultant poleward eddy heat flux. The reduced net normalized conversion/generation rate in midwinter is also found to arise in part from the seasonally enhanced kinetic energy conversion from eddies into the strongly diffluent Pacific jet around its exit. The seasonality of the net energy influx also contributes especially to the spring recovery of the net normalized conversion/generation rate. The midwinter reduction in the normalized rates of both the net energy conversion/generation and baroclinic energy conversion was more pronounced in the period before the late 1980s, during which the MWM of the storm-track activity was climatologically more prominent.

Jupiters atmosphere has a system of zones and belts punctuated by small and large vortices, the largest being the Great Red Spot. How these features change with depth is unknown, with theories of their structure ranging from shallow meteorological features to surface expressions of deep-seated convection. We present observations of atmospheric vortices using the Juno spacecrafts Microwave Radiometer. We found vortex roots that extend deeper than the altitude at which water is expected to condense, and we identified density inversion layers. Our results constrain the three-dimensional structure of Jupiters vortices and their extension below the clouds.

Zonal jets are common in planetary atmospheres. Their character, structure, and seasonal variability depend on the planetary parameters. During solstice on Earth and Mars, there is a strong westerly jet in the winter hemisphere and weak, low-level westerlies in the ascending regions of the Hadley cell in the summer hemisphere. This summer jet has been less explored in a broad planetary context, both due to the dominance of the winter jet and since the balances controlling it are more complex, and understanding them requires exploring a broader parameter regime. To better understand the jet characteristics on terrestrial planets and the transition between winter- and summer-dominated jet regimes, we explore the jet's dependence on rotation rate and obliquity. Across a significant portion of the parameter space, the dominant jet is in the winter hemisphere, and the summer jet is weaker and restricted to the boundary layer. However, we show that for slow rotation rates and high obliquities, the strongest jet is in the summer rather than the winter hemisphere. Analysis of the summer jet's momentumbalance reveals that the balance is not simply cyclostrophic and that both boundary layer drag and vertical advection are essential. At high obliquities and slow rotation rates, the cross-equatorial winter cell is wide and strong. The returning poleward flow in the summer hemisphere is balanced by low-level westerlies through an Ekman balance and momentum is advected upward close to the ascending branch, resulting in a midtroposphere summer jet.

Juno microwave radiometer (MWR) observations of Jupiter's midlatitudes reveal a strong correlation between brightness temperature contrasts and zonal winds, confirming that the banded structure extends throughout the troposphere. However, the microwave brightness gradient is observed to change sign with depth: the belts are microwave-bright in the p10 bar range. The transition level (which we call the "jovicline") is evident in the MWR 11.5 cm channel, which samples the 5-14 bar range when using the limb-darkening at all emission angles. The transition is located between 4 and 10 bars, and implies that belts change with depth from being NH3-depleted to NH3-enriched, or from physically warm to physically cool, or more likely a combination of both. The change in character occurs near the statically stable layer associated with water condensation. The implications of the transition are discussed in terms of ammonia redistribution via meridional circulation cells with opposing flows above and below the water condensation layer, and in terms of the "mushball" precipitation model, which predicts steeper vertical ammonia gradients in the belts versus the zones. We show via the moist thermal wind equation that both the temperature and ammonia interpretations can lead to vertical shear on the zonal winds, but the shear is similar to 50x weaker if only NH3 gradients are considered. Conversely, if MWR observations are associated with kinetic temperature gradients then it would produce zonal winds that increase in strength down to the "jovicline", consistent with Galileo probe measurements; then decay slowly at higher pressures.

During the past few years, both the Cassini mission at Saturn and the Juno mission at Jupiter provided measurements with unprecedented accuracy of the gravity and magnetic fields of the two gas giants. Using the gravity measurements, it was found that the strong zonal flows observed at the cloud level of the gas giants are likely to extend thousands of kilometres deep into the planetary interior. However, the gravity measurements alone, which are by definition an integrative measure of mass, cannot constrain with high certainty the exact vertical structure of the flow. Taking into account the recent Cassini magnetic field measurements of Saturn, and past secular variations of Jupiter's magnetic field, we obtain an additional physical constraint on the vertical decay profile of the observed zonal flows on these planets. Our combined gravity-magnetic analysis reveals that the cloud-level winds on Saturn (Jupiter) extend with very little decay, i.e. barotropically, down to a depth of around 7000 km (2000 km) and then decay rapidly in the semiconducting region, so that within the next 1000 km (600 km) their value reduces to about 1 per cent of that at the cloud level. These results indicate that there is no significant mechanism acting to decay the flow in the outer neutral region, and that the interaction with the magnetic field in the semiconducting region might play a central role in the decay of the flows.

Uranus and Neptune are the archetypes of \u201cice giants\u201d, a class of planets that may be among the most common in the Galaxy. They are the last unexplored planets of the Solar System, yet they hold the keys to understand the atmospheric dynamics and structure of planets with hydrogen atmospheres inside and outside the solar system.

Jupiters Great Red Spot (GRS) is the largest atmospheric vortex in the Solar System and has been observed for at least two centuries. It has been unclear how deep the vortex extends beneath its visible cloud tops. We examined the gravity signature of the GRS using data from 12 encounters of the Juno spacecraft with the planet, including two direct overflights of the vortex. Localized density anomalies due to the presence of the GRS caused a shift in the spacecraft line-of-sight velocity. Using two different approaches to infer the GRS depth, which yielded consistent results, we conclude that the GRS is contained within the upper 500 kilometers of Jupiters atmosphere.

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. Remote sensing observations, from both visiting spacecraft and Earth-based astronomical facilities, have revealed the significant variation in environmental conditions from one band to the next. 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. Numerical simulations successfully reproduce the former, whereas there is a wealth of observational evidence in support of the latter. 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," with a natural transition from zonal jet pumping to dissipation as we move from the convectively-unstable mid-troposphere into the stably-stratified upper troposphere.

The Ferrel cell consists of the zonal mean vertical and meridional winds in the midlatitudes. The continuity of the Ferrel circulation and the zonal mean momentum and heat budgets imply a collocation of the eddy-driven jet and poleward eddy heat flux maxima, under certain assumptions, including the negligibility of diabatic heating. The latter assumption is questioned, since midlatitude storms are associated with latent heating in the midtroposphere. In this study, the heat budget of the Ferrel cell in both hemispheres is examined, using the JRA55 reanalysis data set. The diabatic heating rate is significant close to the center of the Ferrel cell during winter and at the ascending branch during summer in both hemispheres. The interannual variability shows a positive correlation between the diabatic heating rate in the midlatitude midtroposphere and the latitudinal separation between the eddy heat flux and the eddy-driven jet maxima during winter in both hemispheres.

The Hadley circulation is a key element of the climate system. It is traditionally defined as the zonally averaged meridional circulation in the tropics, therefore treated as a zonally symmetric phenomenon. However, differences in temperature between land and sea cause zonal asymmetries on Earth, dramatically affecting the circulation. This longitudinal dependence of the meridional circulation evokes questions about where and when the actual large-scale tropical circulation occurs. Here, we look into the connection between the longitudinally dependent meridional circulation, and the actual large-scale transport of air in the tropics using a coupled Eulerian and Lagrangian approach. Decomposing the velocity field into rotational and divergent components, we identify how each component affects the actual circulation. We propose an alternative definition for the circulation, which describes the actual path of air parcels in the tropics, as a tropical atmospheric conveyor belt.

The insolation a planet receives from its parent star is the main driver of the climate and depends on the planet's orbital configuration. Planets with nonzero obliquity and eccentricity experience variations in seasonal insolation. As a result, the climate exhibits 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 nonzero 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 nonzero 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 the maximum temperature latitude. Finally, we discuss the importance of different planetary parameters that affect the climate response to orbital configuration variations.

Superrotation is a dynamical regime where the atmosphere circulates around the planet in the direction of planetary rotation with excess angular momentum in the equatorial region. Superrotation is known to exist in the atmospheres of Venus, Titan, Jupiter, and Saturn in the solar system. Some of the exoplanets also exhibit superrotation. Our understanding of superrotation in a framework of circulation regimes of the atmospheres of terrestrial planets is in progress thanks to the development of numerical models; a global instability involving planetary-scale waves seems to play a key role, and the dynamical state depends on the Rossby number, a measure of the relative importance of the inertial and Coriolis forces, and the thermal inertia of the atmosphere. Recent general circulation models of Venus's and Titan's atmospheres demonstrated the importance of horizontal waves in the angular momentum transport in these atmospheres and also an additional contribution of thermal tides in Venus's atmosphere. The atmospheres of Jupiter and Saturn also exhibit strong superrotation. Recent gravity data suggests that these superrotational flows extend deep into the planet, yet currently no single mechanism has been identified as driving this superrotation. Moreover, atmospheric circulation models of tidally locked, strongly irradiated exoplanets have long predicted the existence of equatorial superrotation in their atmospheres, which has been attributed to the result of the strong day-night thermal forcing. As predicted, recent Doppler observations and infrared phase curves of hot Jupiters appear to confirm the presence of superrotation on these objects.

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.

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 optimization for the flow field using thermal wind balance, we calculate the 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.Plain Language Summary Observations show strong east-west flows at the cloud level of Saturn. These winds are strongest at the equatorial regions, reaching up to 400 m/s, about 4 times stronger than tornado strength winds on Earth. Yet until now we had no knowledge on how deep these winds penetrate into the interior of the gas giant. The gravity experiment executed during the Grand Finale stage (May-August 2017) of the NASA Cassini mission helps answering this question. It is well established that any large-scale motion of the fluid would have a signature in the density distribution and therefore in the planet gravity field. If we can estimate the internal structure and shape of the planet, we might be able to decipher the depth of the winds from its signal in the gravity measurements. Moreover, the rigid-body and flow contribution to gravity field are entangled together, therefore it is necessary to use a wide range of rigid-body models in order to define the wind-induced gravity signal. In this work we propose a solution to the problem. We find that the gravity measurements can be explained with a flow pattern, similar to that observed at the cloud level, penetrating to depths of more than 8,000 km into the planet interior. This has important implications not only for the atmospheric dynamics but also for the interior density structure, composition, magnetic field, and core mass.

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° and a 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 seasonally dependent Hadley circulation.

The depth to which Jupiter's observed east-west jet streams extend has been a long-standing question(1,2). Resolving this puzzle has been a primary goal for the Juno spacecraft(3,4), which has been in orbit around the gas giant since July 2016. Juno's gravitational measurements have revealed that Jupiter's gravitational field is north-south asymmetric(5), which is a signature of the planet's atmospheric and interior flows(6). Here we report that the measured odd gravitational harmonics J(3), J(5), J(7) and J(9) indicate that the observed jet streams, as they appear at the cloud level, extend down to depths of thousands of kilometres beneath the cloud level, probably to the region of magnetic dissipation at a depth of about 3,000 kilometres(7,8). By inverting the measured gravity values into a wind field(9), we calculate the most likely vertical profile of the deep atmospheric and interior flow, and the latitudinal dependence of its depth. Furthermore, the even gravity harmonics J(8) and J(10) resulting from this flow profile also match the measurements, when taking into account the contribution of the interior structure(10). These results indicate that the mass of the dynamical atmosphere is about one per cent of Jupiter's total mass.

The gravity harmonics of a fluid, rotating planet can be decomposed into static components arising from solid-body rotation and dynamic components arising from flows. In the absence of internal dynamics, the gravity field is axially and hemispherically symmetric and is dominated by even zonal gravity harmonics J 2n that are approximately proportional to q n, where q is the ratio between centrifugal acceleration and gravity at the planet's equator. Any asymmetry in the gravity field is attributed to differential rotation and deep atmospheric flows. The odd harmonics, J 3, J 5, J 7, J 9 and higher, are a measure of the depth of the winds in the different zones of the atmosphere. Here we report measurements of Jupiter's gravity harmonics (both even and odd) through precise Doppler tracking of the Juno spacecraft in its polar orbit around Jupiter. We find a north-south asymmetry, which is a signature of atmospheric and interior flows. Analysis of the harmonics, described in two accompanying papers, provides the vertical profile of the winds and precise constraints for the depth of Jupiter's dynamical atmosphere.

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 cases 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 Hadley cell ascending branch on Mars and Titan.Plain Language Summary The Hadley circulation is a thermally driven circulation, meaning that air raises at warm latitudes and descends at colder ones. As the solar forcing is seasonal this cell has a seasonal cycle as well, with typically the winter cell being stronger. Previous studies showed that under planetary conditions where the maximum temperature at solstice is at the summer pole, the ascending branch does not necessarily follow the maximum surface temperature; however, slowing down the rotation rate allows the ascending branch to follow the warmest latitude. In this study, we aim to explain this rotation rate dependence and the Hadley circulation on planets that exhibit strong seasonality like Titan and Mars, by using both theoretical arguments based on angular momentum conservation and idealized 3D model simulations. We find that the rotation rate is the main factor controlling the ascending branch of the circulation.

The atmosphere exhibits two distinct types of jets: the thermally driven subtropical jet and the more poleward eddy-driven jet. Depending on location and season, these jets are often merged or separated, and their position, structure, and intensity strongly influence the eddy fields. Here, the authors study the sensitivity of eddies to changes in the jets' amplitudes and positions in an idealized general circulation model. A modified Newtonian relaxation scheme that has a very short relaxation time for the mean state and a long relaxation time for eddies is used. This scheme makes it possible to obtain any zonally symmetric temperature distribution and is used to systematically modify the jets' amplitudes and locations. It is found that eddies are more sensitive to changes in the amplitude of the eddy-driven jet than to changes in the amplitude of the subtropical jet. Furthermore, when the eddy-driven jet is shifted poleward, eddies tend to intensify. These results are tested for robustness in two different reference simulations: one resembling a situation where the subtropical and eddy-driven jets are nearly merged and one when they are separated.

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-2018. 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 baroclinic wind structures to be analysed, but have made simplifying assumptions that neglected several physical effects. These 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 of these neglected terms, but lacked an overall approach that is able to include all effects 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 analysed, and the results are compared in the barotropic limit with 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 yet can affect the gravitational moments to some degree. However, none of these terms is dominant so any approximation attempting to improve over the simplified thermal wind approach needs to include all other terms.

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 many recently discovered terrestrial exoplanets are expected to hold a wide range of atmospheric masses. Here the dynamic-thermodynamic effects of atmospheric mass on atmospheric circulation are studied using an idealized global circulation model by systematically varying the atmospheric surface pressure. On an Earth analog planet, an increase in atmospheric mass weakens the Hadley circulation and decreases its latitudinal extent. These changes are found to be related to the reduction of the convective fluxes and net radiative cooling (due to the higher atmospheric heat capacity), which, respectively, cool the upper troposphere at mid-low latitudes and warm the troposphere at high latitudes. These together decrease the meridional temperature gradient, tropopause height and static stability. The reduction of these parameters, which play a key role in affecting the flow properties of the tropical circulation, weakens and contracts the Hadley circulation. The reduction of the meridional temperature gradient also decreases the extraction of mean potential energy to the eddy fields and the mean kinetic energy, which weakens the extratropical circulation. The decrease of the eddy kinetic energy decreases the Rhines wavelength, which is found to follow the meridional jet scale. The contraction of the jet scale in the extratropics results in multiple jets and meridional circulation cells as the atmospheric mass increases.

Earth's midlatitudes are dominated by regions of large atmospheric weather variability-often referred to as storm tracks-which influence the distribution of temperature, precipitation and wind in the extratropics. Comprehensive climate models forced by increased greenhouse gas emissions suggest that under global warming the storm tracks shift poleward. While the poleward shift is a robust response across most models, there is currently no consensus on what the underlying dynamical mechanism is. Here we present a new perspective on the poleward shift, which is based on a Lagrangian view of the storm tracks. We show that in addition to a poleward shift in the genesis latitude of the storms, associated with the shift in baroclinicity, the latitudinal displacement of cyclonic storms increases under global warming. This is achieved by applying a storm-tracking algorithm to an ensemble of CMIP5 models. The increased latitudinal propagation in a warmer climate is shown to be a result of stronger upper-level winds and increased atmospheric water vapour. These changes in the propagation characteristics of the storms can have a significant impact on midlatitude climate.

The close-by orbits of the ongoing Juno mission allow measuring with unprecedented accuracy Jupiter's low-degree even gravity moments J₂, J₄, J₆, and J₈. These can be used to better determine Jupiter's internal density profile and constrain its core mass. Yet the largest unknown on these gravity moments comes from the effect of differential rotation, which gives a degree of freedom unaccounted for by internal structure models. Here considering a wide range of possible internal flow structures and dynamical considerations, we provide upper bounds to the effect of dynamics (differential rotation) on the low-degree gravity moments. In light of the recent Juno gravity measurements and their small uncertainties, this allows differentiating between the various models suggested for Jupiter's internal structure.

Comprehensive models of climate change projections have shown that the latitudinal band of extratropical storms will likely shift poleward under global warming. Here we study this poleward shift from a Lagrangian storm perspective, through simulations with an idealized general circulation model. By employing a feature tracking technique to identify the storms, we demonstrate that the poleward motion of individual cyclones increases with increasing global mean temperature. A potential vorticity tendency analysis of the cyclone composites highlights two leading mechanisms responsible for enhanced poleward motion: nonlinear horizontal advection and diabatic heating associated with latent heat release. Our results imply that for a 4 K rise in the global mean surface temperature, the mean poleward displacement of cyclones increases by about 0.85 degrees of latitude, and this occurs in addition to a poleward shift of about 0.6 degrees in their mean genesis latitude. Changes in cyclone tracks may have a significant impact on midlatitude climate, especially in localized storm tracks such as the Atlantic and Pacific storm tracks, which may exhibit a more poleward deflected shape.

The Atlantic and Pacific storm tracks in the Northern Hemisphere are characterized by a downstream poleward deflection, which has important consequences for the distribution of heat, wind, and precipitation in the midlatitudes. In this study, the spatial structure of the storm tracks is examined by tracking transient cyclones in an idealized GCM with a localized ocean heat flux. The localized atmospheric response is decomposed in terms of a time- and zonal-mean background flow, a stationary wave, and a transient eddy field. The Lagrangian tracks are used to construct cyclone composites and perform a spatially varying PV budget. Three distinct mechanisms that contribute to the poleward tilt emerge: transient nonlinear advection, latent heat release, and stationary advection. The downstream evolution of the PV composites shows the different role played by the stationary wave in each region. In the region where the tilt is maximized, all three mechanisms contribute to the poleward propagation of the low-level PV anomaly associated with the cyclone. Upstream of that region, the stationary wave is opposing the former two, and the poleward tendency is therefore reduced. Finally, through repeated experiments with enhanced strength of the heating source, it is shown that the poleward deflection of the storms enhances when the amplitude of the stationary wave increases.

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 wind is indeed superficial, yet there exists a completely decoupled deep flow. 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 cloud-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, and discuss 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. (C) 2017 Elsevier Inc. All rights reserved.

The effect of eddy-eddy interactions on zonal and meridional macroturbulent scales is investigated over a wide range of eddy scales, using high-resolution idealized GCM simulations with and without eddy-eddy interactions. The wide range of eddy scales is achieved through systematic variation of the planetary rotation rate and thus multiple-jet planets. It is found that not only are eddy-eddy interactions not essential for the formation of jets, but the existence of eddy-eddy interactions decreases the number of eddy-driven jets in the atmosphere. The eddy-eddy interactions have little effect on the jet scale, which in both types of simulations coincides with the Rhines scale through all latitudes. The decrease in the number of jets in the presence of eddy-eddy interactions occurs because of the narrowing of the latitudinal region where zonal jets appear. This narrowing occurs because eddy-eddy interactions are mostly important at latitudes poleward of where the Rhines scale is equal to the Rossby deformation radius. Thus, once eddy-eddy interactions are removed, the conversion from baroclinic to barotropic eddy kinetic energy increases, and eddy-mean flow interactions intrude into these latitudes and maintain additional jets there. The eddy-eddy interactions are found to increase the energy-containing zonal scale so it coincides with the jets' scale and thus make the flow more isotropic. While the conversion scale coincides with the most unstable scale, the Rossby deformation radius does not provide a good indication to these scales in both types of simulations.

The poleward propagation of midlatitude storms is studied using a potential vorticity (PV) tendency analysis of cyclone-tracking composites, in an idealized zonally symmetric moist GCM. A detailed PV budget reveals the important role of the upper-level PV and diabatic heating associated with latent heat release. During the growth stage, the classic picture of baroclinic instability emerges, with an upper-level PV to the west of a low-level PV associated with the cyclone. This configuration not only promotes intensification, but also a poleward tendency that results from the nonlinear advection of the low-level anomaly by the upper-level PV. The separate contributions of the upper- and lower-level PV as well as the surface temperature anomaly are analyzed using a piecewise PV inversion, which shows the importance of the upper-level PV anomaly in advecting the cyclone poleward. The PV analysis also emphasizes the crucial role played by latent heat release in the poleward motion of the cyclone. The latent heat release tends to maximize on the northeastern side of cyclones, where the warm and moist air ascends. A positive PV tendency results at lower levels, propagating the anomaly eastward and poleward. It is also shown here that stronger cyclones have stronger latent heat release and poleward advection, hence, larger poleward propagation. Time development of the cyclone composites shows that the poleward propagation increases during the growth stage of the cyclone, as both processes intensify. However, during the decay stage, the vertical alignment of the upper and lower PV anomalies implies that these processes no longer contribute to a poleward tendency.

The relation between the mean meridional temperature gradient and eddy fluxes has been addressed by several eddy flux closure theories. However, these theories give little information on the dependence of eddy fluxes on the vertical structure of the temperature gradient. The response of eddies to changes in the vertical structure of the temperature gradient is especially interesting since global circulation models suggest that as a result of greenhouse warming, the lower-tropospheric temperature gradient will decrease whereas the upper-tropospheric temperature gradient will increase. The effects of the vertical structure of baroclinicity on atmospheric circulation, particularly on the eddy activity, are investigated. An idealized global circulation model with a modified Newtonian relaxation scheme is used. The scheme allows the authors to obtain a heating profile that produces a predetermined mean temperature profile and to study the response of eddy activity to changes in the vertical structure of baroclinicity. The results indicate that eddy activity is more sensitive to temperature gradient changes in the upper troposphere. It is suggested that the larger eddy sensitivity to the upper-tropospheric temperature gradient is a consequence of large baroclinicity concentrated in upper levels. This result is consistent with a 1D Eady-like model with nonuniform shear showing more sensitivity to shear changes in regions of larger baroclinicity. In some cases, an increased temperature gradient at lower-tropospheric levels might decrease the eddy kinetic energy, and it is demonstrated that this might be related to the midwinter minimum in eddy kinetic energy observed above the northern Pacific.

Geostrophic turbulence theory predicted already a few decades ago an inverse energy cascade in the barotropic mode, yet there has been limited evidence for it in the ocean. In this study, the latitudinal behavior of the oceanic barotropic energy balance and macroturbulent scales is studied using the ECCO2 (Estimating the Circulation and Climate of the Ocean) state estimate, which synthesizes satellite data and in situ measurements with a high-resolution general circulation model containing realistic bathymetry and wind forcing. It is found that inverse energy cascade occurs at high latitudes, as eddy-eddy interactions spread the conversion of eddy kinetic energy from the baroclinic to the barotropic mode, both upscale and downscale. At these latitudes, the conversion scale of baroclinic eddy kinetic energy and the energy-containing scale follow the most unstable and Rhines scales, respectively. Even though an inverse energy cascade occurs at high latitudes, the energy spectrum follows a steeper slope than the -5/3 slope. Different than classic arguments, the Rossby deformation radius does not follow the baroclinic conversion and most unstable scales.

Slantwise convection should be ubiquitous in the atmospheres of rapidly rotating fluid planets. We argue that convectively adjusted lapse rates should be interpreted along constant angular momentum surfaces instead of lines parallel to the local gravity vector. Using Cassini wind observations of Jupiter and different lapse rates to construct toy atmospheres, we explore parcel paths in symmetrically stable and unstable weather layers by the numerically modeled insertion of negatively buoyant bubbles. Low-Richardson number atmospheres are very susceptible to transient symmetric instability upon local diabatic forcing, even outside of the tropics. We explore parcel paths in symmetrically stable and unstable weather layer environments, the latter by adding thermal bubbles to the weather layer. Parcels that cool in Jupiter's belt regions have particularly horizontal paths, with implications for jetward angular momentum fluxes. These considerations may be relevant to the interpretation of Juno's ongoing observations of Jupiter's weather layer.

Poleward migration of eddy-driven jets is found to occur in the extratropics when the subtropical and eddy-driven jets are clearly separated, as achieved by simulations at high-rotation rates. The poleward migration of these eddy-driven baroclinic jets over time is consistent with variation of eddy momentum flux convergence and baroclinicity across the width of the jet. We demonstrate this using a high-resolution idealized GCM where we systematically examine the eddy-driven jets over a wide range of rotation rates (up to 16 times the rotation rate of Earth). At the flanks of the jets, the poleward migration is caused by a poleward bias in baroclinicity across the width of the jet, estimated through measures such as Eady growth rate and supercriticality. The poleward biased baroclinicity is due to the meridional variation of the Coriolis parameter, which causes a poleward bias of the eddy momentum flux convergence. At the core of the jets, the poleward biased eddy momentum flux convergence relative to the mean jet deflects over time the baroclinicity and the jets poleward. As the rotation rate is increased, and more (narrower) jets emerge the migration rate becomes smaller due to less eddy momentum flux convergence over the narrower baroclinic zones. We find a linear relation between the migration rate of the jets and the net eddy momentum flux convergence across the jets. This poleward migration might be related to the slow poleward propagation of temporal anomalies of zonal winds observed in the upper troposphere.

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 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 planetsbroadly definedand 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 \u201cextratropics\u201d) 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 \u201ctropics\u201d). 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.

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 three-dimensional, 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 s^-1, and vertical velocities that advect air over a scale height in ∼10⁵-10⁶ s. 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.

This study demonstrates that water vapor transport and precipitation are largely modulated by the intensity of the subtropical jet, transient eddies, and the location of wave breaking events during the different phases of ENSO. Clear differences are found in the potential vorticity (PV), meteorological fields, and trajectory pathways between the two different phases. Rossby wave breaking events have cyclonic and anticyclonic regimes, with associated differences in the frequency of occurrence and the dynamic response. During La Ni~na, there is a relatively weak subtropical jet allowing PV to intrude into lower latitudes over the western United States. This induces a large amount of moisture transport inland ahead of the PV intrusions, as well as northward transport to the west of a surface anticyclone. During El Ni~no, the subtropical jet is relatively strong and is associated with an enhanced cyclonic wave breaking. This is accompanied by a time-mean surface cyclone, which brings zonal moisture transport to the western United States. In both (El Ni~no and La Nĩna) phases, there is a high correlation (0.3-0.7) between upper-level PV at 250 hPa and precipitation over the west coast of the United States with a time lag of 0-1 days. Vertically integrated water vapor fluxes during El Nĩno areup to 70kgm^-1s^-1larger than those during La Nĩna along the west coast of the United States. The zonal and meridional moist static energy flux resembles wave vapor transport patterns, suggesting that they are closely controlled by the large-scale flows and location of wave breaking events during the different phase of ENSO.

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 be performed at parameter settings far from jovian values and generally adopt heat fluxes 5-10 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 in these models, 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 Fly. where F is the convective heat flux and v is the viscosity. When the convection is strongly nonlinear, the jet speeds are faster and should scale approximately as (F/v)(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, both Boussinesq and anelastic simulations hint at the existence of a third regime where, at sufficiently high heat fluxes or 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. (C) 2010 Elsevier Inc. All rights reserved.

In winter, northeastern North America and northeastern Asia are both colder than other regions at similar latitudes. This has been attributed to the effects of stationary weather systems set by elevated terrain (orography)(1), and to a lack of maritime influences from the prevailing westerly winds(2). However, the differences in extent and orography between the two continents suggest that further mechanisms are involved. Here we show that this anomalous winter cold can result in part from westward radiation of large-scale atmospheric waves-nearly stationary Rossby waves-generated by heating of the atmosphere over warm ocean waters. We demonstrate this mechanism using simulations with an idealized general circulation model(3-5), with which we show that the extent of the cold region is controlled by properties of Rossby waves, such as their group velocity and its dependence on the planetary rotation rate. Our results show that warm ocean waters contribute to the contrast in mid-latitude winter temperatures between eastern and western continental boundaries not only by warming western boundaries, but also by cooling eastern boundaries.

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.

In this paper it is proposed that baroclinic instability of even a weak shear may play an important role in the generation and stability of the strong zonal jets observed in the atmospheres of the giant planets. The atmosphere is modeled as a two-layer structure, where the upper layer is a standard quasigeostrophic layer on a beta plane and the lower layer is parameterized to represent a deep interior convective columnar structure using a negative beta plane as in Ingersoll and Pollard. Linear stability theory predicts that the high wave-number perturbations will be the dominant unstable modes for a small vertical wind shear like that inferred from observations. Here a nonlinear analytical model is developed that is truncated to one growing mode that exhibits a multiple jet meridional structure, driven by the nonlinear interaction between the eddies. In the weakly supercritical limit, this model agrees with previous weakly nonlinear theory, but it can be explored beyond this limit allowing the multiple jet-induced zonal flow to be stronger than the eddy field. Calculations with a fully nonlinear pseudospectral model produce stable meridional multijet structures when beginning from a random potential vorticity perturbation field. The instability removes energy-from the background weak baroclinic shear and generates turbulent eddies that undergo an inverse energy cascade and form multijet zonal winds. The jets are the dominant feature in the instantaneous upper-layer flow. with the eddies being relatively weak. The jets scale with the Rhines length, but are strong enough to violate the barotropic stability criterion. It is shown that the basic physical mechanism for the generation and stability of the jets in the full numerical model is similar to that of the truncated model.