Ice Giant Circulation Patterns: Implications for Atmospheric Probes

TL;DR

This review synthesizes multi-spectral observations of Ice Giant atmospheres, revealing complex, layered circulation patterns with distinct latitude domains, informing future probe mission site selection and atmospheric modeling.

Contribution

It provides a comprehensive observational framework for Ice Giant circulation, integrating diverse measurements to identify layered cells and potential vertical separation mechanisms.

Findings

Identifies three latitude domains with distinct circulation characteristics.

Suggests a tiered structure of stacked circulation cells in the atmosphere.

Proposes mechanisms for vertical separation of circulation layers.

Abstract

Atmospheric circulation patterns derived from multi-spectral remote sensing can serve as a guide for choosing a suitable entry site for a future in situ probe mission. Since the Voyager-2 flybys in the 1980s, three decades of observations from ground- and space-based observatories have generated a picture of Ice Giant circulation that is complex, perplexing, and altogether unlike that seen on the Gas Giants. This review seeks to reconcile the various competing circulation patterns from an observational perspective, accounting for spatially-resolved measurements of: zonal albedo contrasts and banded appearances; cloud-tracked zonal winds; temperature and para-H$_2$ measurements above the condensate clouds; and equator-to-pole contrasts in condensable volatiles (methane and hydrogen sulphide) in the deeper troposphere. These observations identify three distinct latitude domains: an equatorial domain of deep upwelling and upper-tropospheric subsidence, potentially bounded by peaks in the retrograde zonal jet and analogous to Jovian cyclonic belts; a mid-latitude transitional domain of upper-tropospheric upwelling, vigorous cloud activity, analogous to Jovian anticyclonic zones; and a polar domain of strong subsidence, volatile depletion, and small-scale (and potentially seasonally-variable) convective activity. Taken together, the multi-wavelength observations suggest a tiered structure of stacked circulation cells (at least two in the troposphere and one in the stratosphere), potentially separated in the vertical by (i) strong molecular weight gradients associated with cloud condensation, and by (ii) transitions from a thermally-direct circulation regime at depth to a wave-driven circulation regime at high altitude. The inferred circulation can be tested in the coming decade by 3D simulations and by observations from future world-class facilities. [Abridged]

Full text

††journal: Space Science Reviews

1. 1 Introduction

2. 2 Upper Tropospheric Circulation: Winds, Temperatures, and Reflectivity

3. 2.1 Aerosols and Albedo

4. 2.1.1 Ice Giant Banding

5. 2.1.2 Temporal Variability in Reflectivity

6. 2.2 Zonal Winds

7. 2.3 Atmospheric Temperatures and Ortho/Para-H2

8. 3 Mid-Tropospheric Circulation: Equator-to-Pole Contrasts

9. 3.1 Atmospheric Composition

10. 3.2 Methane

11. 3.3 H2S

12. 4 Stacked Circulation Cells

13. 4.1 Tropospheric Circulation

14. 4.2 Connection to the Interior

15. 5 Stratospheric Circulation: Chemical Tracers

16. 6 Unresolved Equatorial Winds

17. 7 Conclusion: Where to target a probe?

1 Introduction

Although three decades have passed since the Voyager 2 spacecraft encountered Uranus and Neptune, our understanding of Ice Giant meteorology and atmospheric circulation remains in its infancy, largely due to the tremendous challenge of observing these distant worlds. The scientific case for an in situ entry probe for one or both of these worlds (Mousis et al., 2018) rests on its ability to uniquely determine atmospheric composition and structure during its descent, providing access to chemical species and altitude domains that are inaccessible to remote sensing. However, our only previous experience with Giant Planet entry probes (the descent of the Galileo probe into Jupiter in 1995) demonstrated that the interpretation of measurements of condensable species required a good understanding of the local meteorology, which is in turn determined by larger-scale atmospheric circulation patterns (Orton et al., 1998). The experience with Galileo argues in favour of multi-probe missions, but it is clear that any future atmospheric probes must be carefully targeted to maximise the scientific return - such as the deliberate targeting of strong regions of upwelling, or the desire to sample more ‘representative‘ regions of an Ice Giant. This article reviews our current understanding of Ice Giant circulation patterns as determined from remote sensing, revealing that different altitude domains (the stratosphere, upper troposphere, and deeper troposphere below the clouds) may exhibit different patterns, leading to a stacked tier of different - but connected - circulation cells.

How do we explore Ice Giant circulation patterns? Remote sensing is required to diagnose these circulation regimes, to serve as a guide for the scientifically optimal locations for probe entry. These observations can be subdivided into three categories: (1) observations of reflected sunlight in the ultraviolet, visible, and near-infrared, probing in and out of strong methane absorption bands to sense the aerosol distribution (condensate clouds and photochemical hazes) as a function of altitude; (2) observations of thermal emission in the mid-infrared, far-infrared, sub-millimetre and radio, sensing atmospheric temperatures and gaseous composition from the stratosphere to the deep troposphere; and (3) observations of thermospheric emission from H2 and H${}_{3}^{+}$ to determine the circulation patterns in the upper atmosphere. This review will focus on the first two categories. Diagnosing atmospheric circulation requires high spatial resolution, which for reflected light and thermal infrared observations either demands the use of 8-to-10-m diameter observatories on the ground (e.g., Gemini, Keck, Subaru and the Very Large Telescope, VLT), or the stable conditions of the Hubble Space Telescope (HST). Furthermore, amateur observers have recently begun tracking prominent atmospheric storms and bands on both Uranus and Neptune (Hueso et al., 2017). Although these are not yet of sufficient quality to aid in the exploration of Ice Giant circulation patterns, we might reasonably expect significant improvements in the decade before any Ice Giant mission.

As spatial resolution decreases with increasing wavelength, the majority of the literature deals with imaging and spectroscopy at visible and near-infrared wavelengths, providing insights into the distribution of clouds and hazes as a function of time, and revealing the banded structure of the Ice Giant clouds (see the recent review by Sanchez-Lavega et al., 2019, and references therein). Thermal-infrared observations from Voyager 2 (Conrath et al., 1998) and ground-based facilities (Orton et al., 2007; Fletcher et al., 2014; Orton et al., 2015; Roman et al., 2020) have revealed structures on large scales, but not at the same resolution as the cloud banding. Finally, centimetre and millimetre-wavelength arrays like the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter Array (ALMA) (de Pater and Gulkis, 1988; Hofstadter and Butler, 2003; de Pater et al., 2014) provide maps of gaseous contrasts at greater depths, beneath the levels sampled by images in reflected sunlight. In particular in the past few years, after a sensitivity upgrade to the VLA and commissioning of most ALMA antennas, such maps rival the reflected sunlight images in terms of spatial resolution (Tollefson et al., 2019; de Pater et al., 2018). This review will seek to piece together these different strands of observational evidence to generate a picture of Ice Giant circulations that can be tested by future observing campaigns.

What drives atmospheric circulation? Planetary atmospheres respond to differences in energy inputs as a function of altitude, location, and time, resulting in a circulation that is a delicate balance between solar inputs from above (with the axial tilt generating seasonally-dependent hemispheric contrasts) and spatially-variable heating from internal sources (e.g., residual energy from planetary formation, rain-out, or other ongoing gravitational settling driving convective motions, atmospheric instabilities, etc.). This circulation is largely axisymmetric as a result of the planetary rotation, and the latitudinal extent of the circulation cells depends on the rotation period (e.g., Held and Hou, 1980). Our terrestrial troposphere features a thermally-direct Hadley circulation cell in the tropics (i.e., air rises near the equator where it is warm, and sinks where it is cold), which is prevented from extending all the way to the poles by the angular momentum of the rotating Earth. In addition, a weaker thermally-indirect eddy-driven Ferrel circulation cell exists in the extra-tropics, which exhibits rising air at its polar-boundary, equatorward transport at altitude, and subsidence at the edge of the thermally-direct Hadley cell (e.g., Showman et al., 2013). For the latter, the generation of eddies on small scales provides a ‘stirring’ mechanism to generate larger-scale Rossby waves (Vallis, 2006), which are able to propagate energy latitudinally away from their generation region, which leads to the flux of momentum back into their generation region to accelerate the extra-tropical jets. Note that the overall circulation at Earth’s mid-latitudes is still thermally direct (e.g., heat is transported polewards), because the heat transport by the Ferrel circulation is overwhelmed by the energy transport associated with the eddies. A similar “Ferrel-like” process may be responsible for the formation of the multiple belts and zones (and their associated zonal jets) on the Gas Giants, Jupiter and Saturn (see recent reviews by Sanchez-Lavega et al., 2019; Showman et al., 2018; Fletcher et al., 2020). The Gas Giants rotate much faster than the Earth ($\sim 10$-hour rotation periods), resulting in finer-scale banding than the terrestrial case. These belts and zones are associated with overturning circulations (upwelling and subsidence) at the scales of the individual bands. Eddy momentum fluxes into the prograde zonal jets have been observed on both Jupiter (e.g., Salyk et al., 2006) and Saturn (Del Genio and Barbara, 2012), and are hypothesised to produce super-rotating equatorial jets on hot Jupiters (Showman and Polvani, 2011). The different ‘flavours’ of the bands (cyclonic or anticyclonic) also exhibit characteristic meteorology on smaller scales.

As Uranus and Neptune have similar 16-to-17-hour rotation periods that shape their global dynamics, one might expect their atmospheric circulation regimes to be intermediate between the Gas Giant and terrestrial cases. However, our knowledge of the formation of the Ice Giant bands, and the relation with the observed wind field and storms, is nowhere near the maturity of our understanding of the Gas Giants from the Voyager, Galileo, Juno, and Cassini missions. Given the appearance of the planets in reflected sunlight, it is not immediately clear even how to define an Ice Giant belt or zone. As on Saturn, where the relationship between the temperature/wind fields and the visible banding is unclear (Showman et al., 2018), we will define Ice Giant bands as cyclonic belts (where the circulation between the peripheral jets is in the same sense as the planetary rotation) and anticyclonic zones (where the circulation opposes the planetary rotation), based on their measured temperature contrasts and zonal winds. Axisymmetric albedo contrasts within these cyclonic and anticyclonic regions will be discussed below. Furthermore, we might expect stark differences between the circulation patterns on Uranus and Neptune. Uranus’ atmosphere is unique in the Solar System, receiving negligible heat flux from the deep interior and experiencing extreme seasonal forcing due to the $98^{\circ}$ obliquity. Neptune, with its powerful self-luminosity ($2.61\pm 0.28\times$ the solar input, Pearl and Conrath, 1991) and strong meteorological activity (dark and drifting ovals, rapidly-evolving bright clouds), provides an important counter-example of a convectively-active Ice Giant weather layer. Taken together, these two worlds may be end-members of a whole category of astrophysical object (the Neptunes and sub-Neptunes, Fulton and Petigura, 2018), and provide an extreme test of our understanding of atmospheric circulation.

What are Ice Giants made of? Before exploring Ice Giant circulation, we first briefly review their composition and cloud structure. While the Gas Giants Jupiter and Saturn have a composition that is, broadly speaking, similar to that of our Sun, the Ice Giants Uranus and Neptune are much more enriched in “heavy” elements (more massive than He). In all Giant Planet atmospheres, we expect the elements O, N, S, and C to be in the form of H2O, NH3, H2S, and CH4. At the higher altitudes, where the temperature is lower, these gases will condense when their partial pressure exceeds their saturated vapour pressure curve. Moreover, some NH3 and H2S will be dissolved in a deep water (solution) cloud, and at higher altitudes we expect NH3 and H2S to combine to form a solid NH4SH cloud layer. For near-solar composition atmospheres (or atmospheres in which heavy elements are enhanced uniformly over solar values, as on Jupiter) we thus expect a deep water or solution cloud to form, topped off by water ice, an NH4SH cloud, and above that an NH3-ice cloud. If it gets cold enough, like on the Ice Giants, the topmost layer will be CH4-ice (Lewis, 1969; Weidenschilling and Lewis, 1973; Atreya and Romani, 1985). Microwave observations suggested that H2S on Uranus and Neptune was much more enriched than NH3, when compared to solar composition (Gulkis et al., 1978; de Pater et al., 1991), an observation that was confirmed both by measuring the microwave opacity of H2S in the lab (DeBoer and Steffes, 1996) and by detecting H2S directly in the near-infrared (Irwin et al., 2018, 2019b). Throughout this review we will therefore assume that H2S gas (and not NH3) is present above the NH4SH cloud layer, and that it will form an H2S-ice cloud near the 3-6 bar level.

The structure of this review is as follows. Section 2 explores circulation hypotheses related to reflected sunlight observations of planetary banding, cloud-tracking of zonal winds, and temperature measurements in the thermal-infrared, which all broadly sound the upper troposphere above the condensate clouds (methane and H2S ice). Section 3 then explores observations of equator-to-pole contrasts in the primary gaseous volatiles (methane, ammonia, and hydrogen sulphide) inferred from both near-infrared and microwave spectra, and Section 4 shows how the upper-tropospheric and mid-tropospheric circulations might be reconciled. Section 5 looks at the large-scale overturning circulation in the stratosphere, and how it might be coupled to the circulation patterns and winds in the troposphere. Finally, Section 7 attempts to draw these hypotheses together to understand Ice Giant circulation as a function of depth, concluding with an assessment of where a planetary entry probe might be best targeted.

2 Upper Tropospheric Circulation: Winds, Temperatures, and Reflectivity

Figure 1 presents a comparison of the latitudinal distributions of albedo (Section 2.1), zonal winds (Section 2.2), temperatures and para-H2 (Section 2.3), methane (Section 3.2), and microwave absorbers like H2S and NH3 (Section 3.3), which will be referred to repeatedly through the following sections. We note that temporal variations in albedo have been observed on both worlds, related to weather phenomena, seasonal changes, and solar-driven variability. We mention them briefly in the sections that follow, but given that we have only been able to resolve atmospheric features for about 1.5 Uranian seasons (each season lasting 21 years), and have yet to see in detail even a complete 41-year Neptunian season, our understanding of these variations is necessarily limited.

2.1 Aerosols and Albedo

2.1.1 Ice Giant Banding

Jupiter provides the archetype for understanding planetary banding, with a clear distinction between reflective zones and darker belts, and a close correlation between albedo and the latitudinal variations of temperatures and winds. However, we should not let our historical familiarity with Jupiter’s bands bias our understanding of atmospheric circulation on the Ice Giants. The atmospheres of Uranus (Smith et al., 1986; Lindal et al., 1987; Tyler et al., 1986) and Neptune (Smith et al., 1989; Conrath et al., 1989; Broadfoot et al., 1989; Tyler et al., 1989; Lindal, 1992) were first explored by the Voyager-2 spacecraft in 1986 and 1989, respectively, and revealed worlds that were dramatically different from Jupiter. Planetary banding on Neptune appeared to be much more subtle than Jupiter’s, and the bland appearance of Uranus suggested a wholly different type of circulation regime to the other giants (e.g., see the comprehensive review by Sanchez-Lavega et al., 2019). However, modern image processing techniques, and improvements in Earth-based observations, are beginning to shift this view.

Voyager 2 encountered Uranus near its southern summer solstice, with the northern hemisphere hidden in winter darkness. Visible-light maps of the southern hemisphere reprocessed by Karkoschka (2015) showed that Uranus was not so bland after all: although somewhat subjective, these albedo maps showed narrow, reflective bands near $25-30^{\circ}$S, $52^{\circ}$S, $68^{\circ}$S and $77^{\circ}$S, with finer-scale dark banding observed equatorward of $45^{\circ}$S and between $68-77^{\circ}$S. The darkest and blandest band was observed between $78-83^{\circ}$S, and can be seen in Fig. 1A. Uranus’ zonal banding was found to be probably due to variations in aerosol optical depth, although one of the bands could have been caused by variations in aerosol absorption, suggestive of different materials in some bands (Karkoschka, 2015).

As Uranus passed through northern spring equinox in 2007, the northern hemisphere came into view and could be captured via advanced ground-based imaging techniques, providing the unprecedented views of fine-scale banding shown in Fig. 1A. Bright storm clouds reaching high above the surrounding clouds (presumably of methane ice) became increasingly visible in the years following equinox, and have been tracked with high-resolution ground-based imaging (de Pater et al., 2015; Irwin et al., 2017; Sromovsky et al., 2015). Bright features can be readily seen in methane bands and at red wavelengths (Sromovsky et al., 1995). Long exposures in the near-infrared (H band, 1.4-1.8 $\mu$m) using instruments on Keck and Gemini in 2012-14 were ‘derotated‘ using a knowledge of the zonal wind field in order to enhance the visibility of features in Fig. 1A (Fry et al., 2012; Sromovsky et al., 2015). Images at these wavelengths reveal Uranus’ banded structure, albeit interrupted by bright clouds, vortices, and ephemeral storms. Zonal medians of the near-infrared reflectivity maps (Sromovsky et al., 2015) indicate that these zonal contrasts are not static, but change from observation to observation, potentially due to obscuration of the banded structure by discrete features. Some persist - brighter bands in Fig. 1A could be seen between $40-50^{\circ}$S, $10-20^{\circ}$S, $0-8^{\circ}$N, $10-12^{\circ}$N, $18-31^{\circ}$N, $38-42^{\circ}$N, and $48-52^{\circ}$N, along with fainter zonal albedo contrasts on a $\sim 5^{\circ}$-latitude scale. Sromovsky et al. (2015) also found an equatorial wave feature with diffuse bright features every $30-40^{\circ}$ longitude. In the north polar region, brightness minima occurred near $53-54^{\circ}$N, $60-61^{\circ}$N, $70-71^{\circ}$N, and $76-80^{\circ}$N, and Sromovsky et al. (2015) showed that the springtime polar region beyond $55^{\circ}$N was characterised by small, bright cloud features with 600-800 km diameters, potentially of convective origin. This also corresponded to a region of solid body rotation, as explained in Section 2.2. Taken together, the contrast-enhanced views of Uranus’ southern (Karkoschka, 2015) and northern (Sromovsky et al., 2015) hemispheres revealed banded albedo patterns on a much finer scale than the temperature and wind fields, which is totally unlike Jupiter and Saturn.

Neptune’s banded pattern was evident in Voyager-2 flyby images (Smith et al., 1989; Limaye and Sromovsky, 1991; Sromovsky et al., 2001), although the overall morphology changed considerably between 1989 and 2001. In particular, the bright bands with fine-scale zonal structures at southern and northern mid-latitudes became very prominent only after the Voyager flyby (Sromovsky et al., 2002; Karkoschka, 2011; Martin et al., 2012). Given Neptune’s long orbital period, we have only ever been able to observe the southern hemisphere (summer solstice was in 2005). Neptune in Fig. 1A exhibits bright and variable cloud activity at mid-latitudes, with a number of narrow, bright bands near $25-50^{\circ}$S and $25-45^{\circ}$N, but bands of lower reflectivity equatorward of $\pm 25^{\circ}$ and poleward of $50^{\circ}$S (de Pater et al., 2014). Neptune’s equator is relatively quiescent, lacking the bright cloud activity compared to the stormy mid-latitudes, where cloud activity seems to peak in latitude bands centred near $25^{\circ}$S and $30^{\circ}$N. Large anticyclones on Uranus and Neptune can form and dissipate on the timescales of years (Hammel et al., 1995). A small dark spot was observed on Uranus in 2006 at $28^{\circ}$N (Hammel et al., 2009), and since that time several more have been observed (e.g., Sromovsky et al., 2015). Before disappearing by 1994, the Great Dark Spot dominated the Voyager-2 flyby of Neptune (Smith et al., 1989; Baines and Hammel, 1994; Sromovsky et al., 1993; LeBeau and Dowling, 1998), and drifted equatorward before disappearing. Dark ovals on both Ice Giants are usually accompanied by bright companion clouds due to air being forced upwards over the underlying vortex (orographic clouds, Stratman et al., 2001). Since the Voyager flyby, several new dark spots have been observed on Neptune by Hubble, allowing a comprehensive exploration of their drift rates and lifetimes (Sromovsky et al., 2001; Wong et al., 2018; Simon et al., 2019). Since the dark spots on Uranus are much smaller than those on Neptune, they are harder to observe, and may be more rare on this planet. The latitudinal drifting of these dark anticyclones sets the Ice Giants apart from the Gas Giants, where vortices remain in their latitude bands due to the strong shears associated with the zonal winds. This again hints at atmospheric circulation patterns that differ between the Gas and Ice Giants.

2.1.2 Temporal Variability in Reflectivity

Observations of Uranus now span 1.5 Uranian seasons (Uranus’ year is 84 Earth years long), which has revealed seasonal changes in its albedo. High southern latitudes (poleward of $45^{\circ}$S) were found to be more reflective than mid-latitudes during the Voyager encounter (potentially due to the increased optical depth of the methane cloud near 1.2-1.3 bar, Rages et al., 1991), with an absence of any small-scale convective structures (Sromovsky et al., 2012) suggesting suppressed convection at the summer pole under a ‘south polar cap’, or convective features being hidden by this cap of high-altitude aerosols. This albedo pattern changed over the subsequent years: Hubble 700-1000 nm images between 1994-2002 (Rages et al., 2004) revealed a darkening of the south pole, the formation of a bright ring near $70^{\circ}$S, and a south polar collar at $45^{\circ}$S. As the northern hemisphere came into view after 2007, the south polar collar diminished in brightness (Irwin et al., 2010; Roman et al., 2018), and a north polar collar became steadily brighter at $45^{\circ}$N (Irwin et al., 2012), with a bright ‘north polar cap’ observed after 2014 (de Pater et al., 2015; Sromovsky et al., 2015; Toledo et al., 2018), as shown in the centre of Fig. 2. Unlike the southern polar region, the north polar region exhibited numerous small-scale features in 2012-14, suggestive of convective clouds and an asymmetry between the southern summer and northern spring poles (Sromovsky et al., 2012). These cloud features are still visible at the time of writing (2019).

These changes in Uranus’ albedo may be due to both a high-latitude depletion of methane (Section 3.2), coupled with seasonally-changing aerosols. Stratospheric aerosols were revealed by Voyager 2 (Pollack et al., 1987; Rages et al., 1991), most likely related to photochemistry and the resulting condensation of hydrocarbon ices. But the rapid decline of the south polar cap and the development of a cap in the north is much faster than would be expected from radiative timescales (Conrath et al., 1990), photochemical timescales (Moses et al., 2018), and aerosol microphysical timescales (Toledo et al., 2019). Furthermore, stratospheric hazes formed from the condensation of photochemically-produced hydrocarbons would need to be transported meridionally from upwelling regions to subsiding regions in the upper troposphere to influence the reflectivity of the $p\sim 1$-bar cloud layer (there is insufficient UV-penetration to the 1-2 bar region to drive photochemical haze production in situ). Accumulation, sedimentation and coagulation are relatively slow processes that would produce a substantial seasonal lag, again inconsistent with the rapid changes observed on Uranus (Toledo et al., 2019). This all hints at meridional circulation patterns in the stratosphere and upper troposphere, to which we will return in the coming sections. Finally, we note that Uranus is not alone in displaying albedo variations - long-term monitoring of Neptune’s disc-averaged visible magnitude shows a consistent brightening over many decades (Lockwood and Jerzykiewicz, 2006; Hammel et al., 2007; Lockwood, 2019) and a modulation by solar ultraviolet and galactic cosmic rays that drive aerosol nucleation via the production of ions (Moses et al., 1992; Aplin and Harrison, 2016). Furthermore, Karkoschka (2011) used fourteen years of HST observations to identify intriguing periodicities of $\sim 5$ years in Neptune’s discrete cloud activity, suggesting that the Ice Giants might exhibit atmospheric cycles on timescales much shorter than a season.

2.2 Zonal Winds

The measurement of Ice Giant winds via cloud tracking probably accounts for the largest body of literature for these two worlds. Winds have been measured from Voyager 2 (Smith et al., 1986, 1989; Hammel et al., 1989; Limaye and Sromovsky, 1991; Sromovsky et al., 1993; Karkoschka, 2015), the Hubble Space Telescope (Sromovsky et al., 1995; Hammel and Lockwood, 1997; Sromovsky et al., 2001; Karkoschka, 1998; Hammel et al., 2001), and ground-based facilities like Keck and Gemini (Fry and Sromovsky, 2004; Sromovsky and Fry, 2005; Hammel et al., 2005; Sromovsky et al., 2009; Martin et al., 2012; Sromovsky et al., 2015; Tollefson et al., 2018). In contrast to Uranus’ stable zonal wind patterns (Sromovsky et al., 2015), there is a good deal of dispersion in the motions of the cloud features on Neptune: it is not always clear that the features are genuinely tracking the underlying winds, and the different features are also likely representative of different altitudes (Fitzpatrick et al., 2014; Tollefson et al., 2018). Nevertheless, these studies have shown that both worlds feature retrograde jets at their equators and a single prograde jet at high latitudes in each hemisphere (Fig. 1C). Neptune’s prograde jet peaks between $70-80^{\circ}$S (a symmetric jet is expected, but not yet observed, in the north), whereas Uranus’ jets peak at $50-60^{\circ}$. These wind patterns are very different from the multi-jet circulation patterns on Jupiter and Saturn (see the comprehensive review by Sanchez-Lavega et al., 2019), leading to the suggestion that some latitudinal gradients of the zonal winds might not have been fully resolved by the measurements to date (see Section 6). Nevertheless, the winds on the two Ice Giants are surprisingly similar considering their different spin axes and internal heat fluxes.

Most of the tracked cloud features on Uranus are near the 1.2-bar methane condensation level or in the deeper 2-3 bar main clouds of H2S ice (Irwin et al., 2018)111Note that infrared observations are most sensitive to the cloud tops, whereas the base pressure could be at higher pressures., but some can reach the 250-600 mbar region in the upper troposphere (Sromovsky et al., 2007, 2012; Roman et al., 2018). A long-lived feature known as the ‘Berg‘ had bright features near 550-750 mbar, but with the main parts of the structure near 1.7-3.5 bars (de Pater et al., 2011). On Neptune, clouds are seen in the stratosphere at 20-60 mbar, at altitudes much higher than those on Uranus (Gibbard et al., 2003; Irwin et al., 2016), making them visible in the K’ filter (sensing strong methane absorption, and therefore high altitudes, between 2.0-2.4 $\mu$m), tracking features in the 20-300 mbar range (Sromovsky et al., 2001), well above the methane condensation level near 1.4 bar. At these altitudes, windshear could add several tens of m/s dispersion to the results in Fig. 1C. Adding to this complexity, Neptune’s clouds evolve over very short timescales (Martin et al., 2012; Simon et al., 2016; Stauffer et al., 2016), displaying both episodic and continuous cloud activity (Baines and Hammel, 1994; Sromovsky et al., 1995; Molter et al., 2019) that can even be observed from the ground by amateur astronomers (Hueso et al., 2017). Neptune’s rapidly evolving convective clouds seem to prevail at mid-latitudes, and clouds in the main storm band at $20-40^{\circ}$S may have become increasingly vigorous since the Voyager flyby (Hammel and Lockwood, 2007). Similar rapid evolution is observed in Uranus’ small-scale clouds (Irwin et al., 2017), making it challenging to track them for long periods, with the highest-intensity storms occurring near $30-40^{\circ}$N (de Pater et al., 2015). The canonical zonal wind profiles in Fig. 1C are therefore subject to considerable uncertainty, but the jets do demarcate a polar domain, mid-latitude domain, and equatorial domain on both worlds.

A crucial open question is how the strength of the winds varies with altitude. Section 2.3 shows how the zonal jets should decay with altitude through the upper troposphere ($p<1$ bar). Analyses of Uranus’ and Neptune’s gravity fields out to the fourth-order harmonic (Hubbard et al., 1991; Kaspi et al., 2013) suggest that the zonal wind patterns are restricted to the outermost 1000 km of the planet’s radii, so the winds must ultimately be decaying with depth too, although the decay function is poorly constrained (Kaspi et al., 2013). Exploring the interface between these two domains is vitally important for an understanding of the processes driving and maintaining the zonal winds. Tollefson et al. (2018) detected vertical wind shear at Neptune’s equator by tracking bright cloud features in the H- (1.4-1.8 $\mu$m) and K’ (2.0-2.4 $\mu$m) bands with Keck in 2013-14. The higher-altitude features at K’ (sensing features at $\sim 10$ mbar) showed stronger retrograde velocities than deeper features seen in the H-band ($p>1$ bar), suggesting that the winds increase in strength with height, opposite to that inferred from the thermal field (Section 2.3). The deep features have a less negative retrograde velocity (by $\sim 100$ m/s) than the higher-altitude features, suggesting a significant windshear, at least at Neptune’s low latitudes (the equatorial region, between $\pm 25^{\circ}$ latitude). A potential reconciliation of the windshear and temperature fields will be presented in Section 2.3.

A final complication is that the zonal wind field is potentially variable with the seasons. On Uranus, the only Ice Giant where winds have been observable for more than a single season (and even here it is only  35 years or 1.5 seasons), an asymmetry between northern spring in 2012-14 (where a broad region of solid-body rotation from 62-83∘N was identified, Sromovsky et al., 2015) and southern summer in 1986 (with a much smaller region of solid-body rotation and a large gradient in drift rates, Karkoschka, 2015) might potentially reverse as Uranus approaches northern summer solstice in 2028. The zonal motions, as well as the contrasting distributions of small-scale convective clouds, might be intricately linked with the seasonal insolation and thus, with any large-scale hemispheric circulation.

2.3 Atmospheric Temperatures and Ortho/Para-H2

Cloud-tracked zonal winds are typically relevant to a narrow altitude range, with the precise cloud-top and base pressures subject to large uncertainties. Thermal-infrared remote sensing, on the other hand, can provide vertical profiles of atmospheric temperature as a function of position, extending the two-dimensional windfields into three dimensions. The upper tropospheric temperatures derived from Voyager 2 in Fig. 1B show cool mid-latitudes in the 80-800 mbar range, contrasted with the warmer equator and poles, as revealed by both Voyager 2 (Flasar et al., 1987; Conrath et al., 1991, 1998) and subsequent ground-based mid-infrared imaging in Fig. 2 (Hammel et al., 2006; Orton et al., 2007; Fletcher et al., 2014; Orton et al., 2015). In the years after the Voyager-2 encounter, a warm summertime vortex developed over Neptune’s southern pole in the upper troposphere and stratosphere (Orton et al., 2007; Fletcher et al., 2014). Note that, like the zonal winds, the temperature bands also demarcate polar, mid-latitude, and equatorial domains, but do not display the fine-scale banding observed in the albedo in Fig. 1A. In the absence of any appreciable latitudinal differences in heating, horizontal temperature contrasts can be used as an indirect measure of vertical motions. For example, if we assume latitudinal uniformity in radiative heating (in gases and aerosol layers), exothermic chemical reactions, latent heat of volatile condensation and/or para-hydrogen equilibriation, then cool temperatures suggest rising motion with adiabatic cooling at mid-latitudes, accompanied by subsidence and adiabatic warming at the equator and poles (Flasar et al., 1987; Conrath et al., 1991; Bézard et al., 1991; Conrath et al., 1998).

Tropospheric temperatures were derived from the 17-50 $\mu$m region of the Ice Giant spectra, which is dominated by a collision-induced continuum of hydrogen and helium. The shape of this continuum is also governed by the ratio of the two spin isomers of H2 (Massie and Hunten, 1982; Conrath et al., 1998): the S(1) absorption near 17 $\mu$m is formed from transitions within ortho-H2 (the odd spin state of H2 with parallel spins), whereas S(0) near 28 $\mu$m is formed from transitions within para-H2 (the even spin state of H2 with anti-parallel spins). Populating or depopulating the S(0) states therefore affects the shape and gradient of the far-IR continuum, with cooler regions near the tropopause having a higher fraction of para-H2. This allows us to identify regions of sub-equilibrium conditions, where the para-H2 fraction is lower than would be expected from thermal equilibrium due to upwelling of low-para-H2 air from the deeper troposphere. Conversely, super-equilibrium conditions are consistent with the subsidence of high-para-H2 air from the tropopause region. Fig. 1D shows the disequilibrium of para-H2, derived from Voyager/IRIS measurements using modern opacities for the collision-induced absorption (Fletcher et al., 2014; Orton et al., 2015; Fletcher et al., 2018), and confirming the presence of mid-latitude upwelling and equatorial/polar subsidence (Flasar et al., 1987; Conrath et al., 1998). However, we caution that the IRIS data were only sensitive to the 100-800 mbar range (i.e., at depths below the tropopause, Conrath et al., 1998), with speculative results at higher altitudes in Fig. 1D resulting from smooth relaxation of retrieved profiles to their priors. At present, para-H2 does not provide any constraint on lower stratospheric circulation in the text that follows.

The Ice Giant circulations inferred from the temperature, para-H2 (below the tropopause), and winds are shown in schematic form in Fig. 3, representing the 0.1-1.0 bar range, approximately. Geostrophy implies that temperatures and winds are in balance with one another via the thermal windshear equation (Holton, 2004), and Flasar et al. (1987) proposed a model where the acceleration of the zonal flow due to conservation of angular momentum was balanced by frictional damping (i.e., vertical shear), potentially due to the breaking of vertically-propagating waves in the upper troposphere. The upper tropospheric temperature gradients in the 80-800 mbar range imply maximum positive windshears near $\pm(15-30)^{\circ}$ latitude (i.e., on the flanks of the equatorial retrograde jet) and maximum negative windshear near $\pm(60-75)^{\circ}$ (i.e., near the prograde jets at high latitudes). This is shown by coloured bars of decreasing size and contrast in Fig. 3. The windshear is minimal (i.e., close to barotropic, with no wind variability with height) in the $\pm 30-50^{\circ}$ latitude range associated with the mid-latitude temperature minima. It is in this mid-latitude domain, with the coolest temperatures and sub-equilibrium para-H2, that the most frequent and notable storm activity occurs, and where the seasonal polar collars emerge on Uranus. The upper tropospheric circulation in Fig. 3 might also suggest a meridional transport of aerosols from regions of mid-latitude upwelling to a region of equatorial subsidence. This is consistent with a reflective band seen near the equator in Hubble and VLT imaging (Toledo et al., 2018), suggesting an accumulation of haze at the equator.

It is interesting to consider whether the symmetric circulation depicted in Fig. 3 could be seasonally variable. For example, averaged over a Uranian year, Uranus receives more heat at the poles than at the equator, so that energy must be transported equatorward. Following earlier work by Friedson and Ingersoll (1987) and Flasar et al. (1987), Conrath et al. (1990) used radiative-dynamical balance equations to explore the relationship between temperatures and winds on Uranus and Neptune as a function of time. Their model features atmospheric motions forced by radiative heating in the stratospheres (by weak methane absorption bands), or by the mechanical forcing from zonal winds imposed at the lower boundaries, without any momentum convergence to accelerate the flow, and with a simplified frictional drag to decay the winds with altitude. For Uranus, the radiative time constant ($\sim 130$ years) was longer than the orbital period, such that the atmospheric temperatures remained close to the annual-average radiative equilibrium values, despite the large amplitude of the radiative forcing, consistent with the lack of observed temperature variations over the full season between solstice (1986), equinox (2007), and the present day (Orton et al., 2015; Roman et al., 2020). Polar regions were only slightly warmer than the equator, depending only on the integrated solar irradiance at the top of the atmosphere, which is some $\sim 30$% larger at the poles. This result would be changed if a latitude-dependent absorbing aerosol were present (see Section 2.1), but the small magnitude of the temperature differences due to Uranus’ small haze opacity should result in a weak annual-mean meridional circulation, with a low-latitude cell of rising motion between $10-30^{\circ}$ in both hemispheres and subsidence at the equator. Li et al. (2018) calculated updated estimates of the radiative heating/cooling rates, using modern estimates of temperature and hydrocarbon profiles. Uranus exhibits a longer radiative time constant than at Neptune because it has the coldest atmosphere and the lowest methane abundance (and resulting photochemical products). As a result, cooling of Uranus is primarily due to the collision-induced opacity of H2. The calculated cooling rates are larger than the heating rates on both worlds, creating the stratospheric energy crisis (i.e., the stratospheres are much warmer than would be calculated from a pure radiative model). Li et al. (2018) noted that their radiative time constants were much shorter than the older results of Conrath et al. (1990), indicating that seasonal effects might be more important than previously ascertained - however, this is inconsistent with the absence of seasonal change on Uranus noted by Orton et al. (2015) and Roman et al. (2020), perhaps suggesting much longer radiative time constants or efficient transport at these atmospheric depths.

We now return to the conundrum in Section 2.2, where Tollefson et al. (2018) had identified a strengthening of the retrograde flow with increasing altitude, counter to the sense of the windshear in Fig. 3. Using a modified thermal wind equation, Tollefson et al. (2018) showed that the required windshear could be produced either by (i) having a cool equator compared to mid-latitudes at $p>1$ bar (i.e., in the opposite sense to that in Fig. 1B and Fig. 2); (ii) by having methane strongly enriched at the equator compared to latitudes $>\pm 40^{\circ}$; or (iii) some combination of the two. They showed that this could be made consistent with the Voyager/IRIS measurements, provided this ‘reversed pattern’ of a cool equator and warm mid-latitudes were restricted to pressures greater than 1 bar. The implied circulation is shown schematically in Fig. 4, showing the retrograde jet weakening into the deeper atmosphere, and a cool equator (upwelling) at $p>1$ bar sitting beneath the warm subsiding region at $p<1$ bar. This deeper circulation is more reminiscent of the cool, upwelling equators of Jupiter and Saturn. Stacked tiers of opposing circulation regimes have been proposed for the gas giants too (see review by Fletcher et al., 2020), to resolve the discrepancies between the meridional motions at the cloud tops associated with eddy momentum convergence on the prograde jets, and meridional motions observed in the upper troposphere (Ingersoll et al., 2000; Showman and de Pater, 2005; Fletcher et al., 2011). We shall return to this interpretation in Section 3, but first we explore further observational constraints on atmospheric circulation in the troposphere and stratosphere.

3 Mid-Tropospheric Circulation: Equator-to-Pole Contrasts

3.1 Atmospheric Composition

The visible and near-infrared reflectivity explored in Section 2.1 depends on the vertical distribution of aerosols in the tropospheres of Uranus and Neptune, which in turn depend on the distributions of condensable volatile species. Based on cosmochemical abundances, and observations of Jupiter and Saturn, the most abundant condensable species in Ice Giant tropospheres are expected to be CH4, H2O, NH3, and H2S. Of these, only CH4 has been directly measured (Karkoschka and Tomasko, 2011, 2009), and is found to be enriched relative to H2 by factors of 10 to 100 times over solar abundance ratios. The bulk mixing ratios of most other condensable species are thought to be similarly enriched, while H2O is expected to be even higher, accounting for  60% of the planet’s total mass (e.g., Guillot, 1995). One must remember, however, that the bulk abundances in the planet may not be reflected in abundances within the troposphere. For example, radio observations found, surprisingly, that NH3 is strongly depleted in at least the upper $\sim 30$ bars of Uranus’ atmosphere, and perhaps deeper (Gulkis et al., 1978). Whether this reflects a misunderstanding of the overall composition of planets (e.g., the lack of nitrogen in Uranus and Neptune could be caused by the inefficient trapping of N2 by clathration in water ice, Hersant et al., 2004), or is due to chemical trapping of NH3 in the interior, is not yet clear.

Based on these chemical abundances, the topmost clouds are expected to be comprised of a thin cloud of CH4 ice with a base near 1.3 bar, where latent heat released by condensation was observed to modify the temperature lapse rate observed in Voyager radio occultations (Lindal et al., 1987; Lindal, 1992). Below, at the 3-6 bar level, lie clouds of H2S ice (de Pater et al., 1991), whereas H2S vapour has been detected above these cloud decks (Irwin et al., 2018, 2019b). Clouds of NH4SH may exist in the 30-to-40-bar region (Weidenschilling and Lewis, 1973; Atreya and Wong, 2005), and are expected to lock away any remaining tropospheric NH3 that is not trapped at much deeper levels (de Pater et al., 1991; Hofstadter and Butler, 2003). Dry or wet adiabatic extrapolation of the temperature profiles from 2 bars (where the temperature was determined by radio occultations, Lindal et al., 1987; Lindal, 1992) to deeper pressure levels are used to estimate where the water cloud base forms. If the atmospheric mixing ratio of H2O were solar, it would occur near 50 bars. Both planetary interior models (Guillot, 1995) and recent measurements of CO (Luszcz-Cook et al., 2013; Cavalié et al., 2017) suggest that the H2O abundance is much higher, potentially several hundred times solar, putting the water cloud base at pressures exceeding 200 bars. This condensation region may be stable against convection, and could separate the troposphere above it from the much deeper atmosphere and the suspected deep, superionic, watery ocean with its own interior circulation patterns at great depths (see Section 4.2). If the condensables are sufficiently abundant, then density stratification (i.e., the mean molecular weight gradient) could have a stabilising influence on convective motions (Gierasch and Conrath, 1987; Guillot, 1995), producing distinct vertical boundaries to heat transport that act like an insulating layer and trap energy. The spatial distribution of volatiles, particularly CH4 (Section 3.2) and H2S (Section 3.3), are therefore crucial to our assessment of atmospheric circulation.

We briefly note that measurements of disequilibrium tropospheric species like CO and PH3 could also provide an indirect way of measuring atmospheric circulation and constraining the deep oxygen abundances (Visscher and Fegley, 2005; Cavalié et al., 2017). PH3 is used as a tracer of vertical mixing on Jupiter and Saturn, but only upper limits are available from millimetre observations on Uranus and Neptune (Moreno et al., 2009; Teanby et al., 2019). Maybe the large water abundance is responsible for converting all the PH3 into other compounds in the 1000-bar region (Fegley and Prinn, 1985; Visscher and Fegley, 2005). Tropospheric CO has been detected on both planets, but Neptune has an extremely high abundance (Marten et al., 1993, 2005; Lellouch et al., 2005, 2010; Fletcher et al., 2010), implying that its deep atmosphere is probably fully convective and well-mixed. Luszcz-Cook et al. (2013) used the far wings of CO lines to show that a brightness temperature increase from Neptune’s southern mid-latitudes to the south pole could be explained by a drop in CO opacity. However, for both PH3 and CO, we have only a very limited sense of their spatial distributions with which to derive atmospheric circulation - they will not be considered again in this review. Instead, we turn to what can be learned from the spatial distributions of CH4 and H2S.

3.2 Methane

Fig. 1E shows the similarities in the spatial distributions of tropospheric methane on both Ice Giants, where there is a significant decline from the equatorial domain to the polar domain. Rather than being a gradual decline, there are hints of a step-like structure, with a transition in the $\pm 30-40^{\circ}$ latitude region on both worlds. Derivations of the methane abundance are usually entangled with assumptions about the scattering properties and spatial distributions of aerosols, but at some near-infrared wavelengths the collision-induced absorption of H2 (near 825 and 1080 nm) can be used to separate these variables, allowing a direct measurement of the methane abundance. Narrow-band observations at these H2-dominated wavelengths, in conjunction with data at 1.29 $\mu$m probing similar depths, have been exploited using data from the Hubble Space Telescope and ground-based near-IR spectroscopy (Sromovsky et al., 2014, 2019).

Uranus’ south polar depletion of methane was first detected via Hubble 300-1000 nm (STIS) spectroscopy in 2002 (Karkoschka and Tomasko, 2009), as the southern hemisphere was approaching autumn conditions. The suggested polar subsidence was thought to be responsible for the notable absence of any small-scale convective structures near the south pole (Sromovsky et al., 2012). However, further Hubble observations in 2012 (Sromovsky et al., 2014) and 2015 (Sromovsky et al., 2019), combined with ground-based data in the near-IR, revealed that the north polar region, now emerging into spring, was similarly depleted in methane and yet did display small-scale convective clouds. The symmetric polar depletions might therefore be separated from the seasonal differences in the appearance of small-scale clouds in the polar domain. These STIS observations, depicted in Fig. 1E, also show that the methane drop from $\sim 4$% near the equator to $\sim 2$% at high latitudes is relatively stable over time. The ongoing brightening of the north polar region during spring (Section 2.1, and shown in the centre of Fig. 2) must therefore be related to changes in aerosol scattering properties, rather than changes in the methane depletion. Uranus’ polar depletion was also observed from ground-based facilities in the near-infrared, including the IRTF (Tice et al., 2013), Keck (de Kleer et al., 2015), Palomar (Roman et al., 2018), the VLT and Gemini in the H (1.4-1.8 $\mu$m) band (Toledo et al., 2018; Irwin et al., 2019b).

Models reproducing the spectroscopy suggest that Uranus’ methane variability is restricted to the upper troposphere (i.e., relatively shallow), with a latitudinally variable decrease in abundance from $\sim 3-5$ bar to the condensation altitude near 1.1 bar (Sromovsky et al., 2014, 2019). The latitudinally-uniform deep abundance for $p>5$ bar depends on the assumed aerosol scattering properties, but ranges from $\sim 3.5$% to $\sim 2.7$% for Uranus. Although the equator-to-pole structure is statistically significant, the methane abundances derived in these studies also show smaller-scale latitudinal structure at a lower significance level. For example, the STIS observations in 2012 and 2015 suggest a local minimum in CH4 abundance right at the equator, flanked by (weak) local maxima at $\pm 5-20^{\circ}$ latitude (Figs. 24-25 of Sromovsky et al., 2019), then a decrease by a factor of 2 to $45^{\circ}$N and a factor of 3 by $60^{\circ}$N, but with abundances that are ultimately sensitive to the aerosol properties. A principal component analysis of Gemini H-band spectra by Irwin et al. (2019b) reproduced the general equator-to-pole contrast, and also showed tentative hints that CH4 is not maximal right at the equator (their Fig. 13). The suggested circulation pattern consists of low-latitude upwelling of moist, CH4-rich air; poleward motion as the air dries by condensation, precipitation and sedimentation; and high-latitude subsidence of CH4-depleted air (and modified by more complex equatorial circulation). As will be described below, this requires alterations to the $p>1$-bar schematic in Fig. 4.

Neptune’s methane distribution, at least for the southern hemisphere visible from Earth, is shown in Fig. 1E, as derived from Hubble STIS 300-1000 nm spectroscopy (Karkoschka and Tomasko, 2011). These results are consistent with a deep mixing ratio of $\sim 4$% at $p>3$ bar, but a decrease by a factor of three from the equator to mid-latitudes in the $p<3$ bar region. Fig. 14 of Karkoschka and Tomasko (2011) shows a transition between elevated CH4 equatorward of $\pm 20^{\circ}$ latitude to depleted CH4 poleward of $45^{\circ}$S, again depicting three distinct regions: a well-mixed equatorial region, a region of decline at mid-latitudes, then a well-mixed polar region with depressed abundances. As on Uranus, smaller-scale and low-significance methane variations are indicated in Fig. 1E, with suggestions of depletion at $25^{\circ}$S and $45-55^{\circ}$S and enhancements near $35^{\circ}$S. Neptune’s $45-55^{\circ}$S region is one of the most active regions of cloud formation (Karkoschka, 2011), but subsidence would be inferred from the methane depletion results there (Karkoschka and Tomasko, 2011). Using millimetre observations sensing the 1-10 bar range, Tollefson et al. (2019) find elevated methane between $12-32^{\circ}$S and $2-20^{\circ}$N, but depressed abundances poleward of $66^{\circ}$S and near the equator at $12^{\circ}$S-$2^{\circ}$N - i.e., a local equatorial minimum flanked by weak off-equatorial maxima, similar to that seen on Uranus. However, Tollefson et al. (2019) stress that their results are not entirely consistent with those of Karkoschka and Tomasko (2011), potentially due to the different modelling approaches in the near-infrared and microwave. A recent principal-component analysis (PCA) of VLT MUSE visible-light spectroscopy of Neptune confirms the equator-to-pole gradient in Neptune’s methane observed by Karkoschka and Tomasko (2011), alongside a slight lowering of methane at the equator consistent with the millimetre observations (Irwin et al., 2019a). It is clear that methane is tracking both a global circulation pattern, alongside more local meteorology on a finer zonal scale.

The connection between the equator-to-pole gradient in methane gas, and the observed distributions of clouds potentially associated with methane ice, is perplexing. Strong subsidence in the polar domain would tend to inhibit convection and CH4 cloud formation (Sromovsky et al., 2014), and yet we see discrete clouds at the poles, and a seasonally-changing opacity of aerosols in Uranus’ 1-2 bar region (Section 2.1). Furthermore, the CH4 distributions do not seem to be tracing the upper tropospheric ($p<1$ bar) mid-latitude upwelling responsible for the cool temperatures, sub-equilibrium para-H2, and polar collars of aerosols in Uranus’ $40-50^{\circ}$ domain. Strong equatorial upwelling might promote cloud and haze formation around the equator, coupled with cool temperatures from adiabatic expansion - neither of which are observed. The hints of small-scale methane variations near the equator might suggest that the upper tropospheric equatorial subsidence (and warming $p<1$ bar) is modulating the $p>1$-bar methane distribution, but the two circulation regimes do seem somewhat separate. The near-IR spectroscopy also suggests that the methane variability is restricted to $p<3-5$ bar. To further investigate this deeper circulation, we now look at the motions inferred from millimetre- and centimetre-wave observations.

3.3 H2S

Millimetre- and centimetre-wave observations (i.e., microwave) of the Ice Giants (de Pater et al., 1991) sense thermal emission modulated by the pressure-broadened wings of NH3 and H2S (along with other contributions from CH4, CO, H2O, potentially PH3, and the hydrogen-helium continuum). Uniquely distinguishing between these contributions is a challenge, and spatial variations of temperatures in the deep troposphere (including the lapse-rate changes associated with cloud condensation) could also contribute to the observed emissions in Fig. 2. Nevertheless, spectral fitting appears consistent with H2S as the dominant absorber in the upper $\sim 20$ bars of the troposphere, with NH3 removed via the formation of the NH4SH cloud just below this level and (probably) dissolution in the aqueous layers at higher pressures. This was recently supported by the detection of H2S absorption features in the H-band near 1.57-1.58 $\mu$m on both Uranus (Irwin et al., 2018) and Neptune (Irwin et al., 2019b) using Gemini-North/NIFS - this suggests that the atmospheric S/N ratio is greater than unity, implying that the main cloud deck visible on both worlds comprise H2S ices. As millimetre and centimetre wavelengths probe beneath the clouds down to 50-100 bars, recent advances in spatially-resolved microwave observations (Fig. 2) have provided an invaluable tool for tracing deep atmospheric circulation.

VLA centimetre observations of Uranus between 1982 and 1994 (de Pater and Gulkis, 1988; de Pater, 1991; Hofstadter and Butler, 2003), probing down to 50 bar, showed an equator-to-pole gradient, with the south pole considerably brighter than the equator, and a boundary somewhere near $45^{\circ}$S. VLA observations since 2003 (Fig. 2) also indicated that the north polar region was similarly bright (Hofstadter et al., 2004; de Pater et al., 2018). ALMA observations in Fig. 2 taken in 2015-2018 at 1-3 mm, probing the 1-10 bar range, also show a bright north polar region (Molter et al., 2020). The latitudinal distribution of H2S, as derived from a combined analysis of ALMA and VLA data, is shown in Fig. 1F (Molter et al., 2020). The low opacity over Uranus’ poles is explained as a complete absence of H2S down to $p\sim 35$ bar, although a small amount of NH3 gas (few times $10^{-7}$) appear to be required (Molter et al., 2020). The NH3 abundance is essentially zero at mid-latitudes and the equator. The images in Fig. 2 further show relatively bright bands at latitudes of $\sim 20^{\circ}$N and $\sim 20^{\circ}$S at 1-3 mm and 1-2 cm, while a third band is visible at the ALMA millimetre wavelengths at the equator. These observations are indicative of a lower opacity at these latitudes at 0.5-5 bar, likely due to variations in the H2S relative humidity and/or CH4 abundance.

These observations show a morphological similarity to the distribution of CH4 in the 1-5 bar range (Sromovsky et al., 2014), suggesting that the polar subsidence (volatile depletion) and equatorial upwelling (volatile enrichment) might extend over great depths. However, the latitudinal distribution may be subtly different at the altitudes of the H2S ice cloud (2-4 bars): a PCA analysis of Gemini H-band observations by Irwin et al. (2019b) suggested that Uranus’ H2S column abundance above the 2-to-4 bar cloud is largest at mid-latitudes, and displays local minima at the equator and north pole (the south pole was not in view). At first glance, this might be more consistent with modulation by the upper tropospheric circulation pattern shown in Fig. 3, with time-variable mid-latitude upwelling being responsible for Uranus’ polar collars.

VLA 1.3-6 cm maps of Neptune (lower part of Fig. 2) probe the 10-50 bar range (de Pater et al., 2014) and showed strong depletion of volatiles at the south pole (Butler et al., 2012; Luszcz-Cook et al., 2013; de Pater et al., 2014) in a region coinciding with a warm summertime polar vortex in the upper troposphere and stratosphere at $p<1$ bar (Hammel et al., 2007; Orton et al., 2007; Fletcher et al., 2014). The warm south polar region extends to approximately $65^{\circ}$S with a low abundance of H2S down to $\sim 40$ bar - this is a smaller region of volatile depletion than at Uranus’ poles, but consistent with Neptune’s prograde jet being at a higher latitude than Uranus’ prograde jet. The centimetre maps from de Pater et al. (2014) do not show an equatorial brightening, but 0.9-3.0 cm observations from the upgraded VLA shown in Fig. 2 (Butler et al., 2012; de Pater et al., 2018) support a picture of equatorial subsidence to the 5-10 bar level, modulating the general equator-to-pole gradient in H2S. Recent ALMA millimetre observations of Neptune (Tollefson et al., 2019), probing the 1-10 bar range, have identified distinct bands of warmer brightness temperatures, explained by H2S abundance variations at $p<10$ bar, with a severe depletion at the south pole. The circulation inferred from these millimetre observations suggests air rising at mid-latitudes ($12-32^{\circ}$S) and north of the equator ($2-10^{\circ}$N), and sinking in the $2^{\circ}$N-$12^{\circ}$S region and poleward of $66^{\circ}$S. As on Uranus, this suggests equator-to-pole transport in the mid-troposphere, coupled with a complicated pattern of equatorial subsidence and near-equatorial upwelling being modulated by the upper tropospheric circulation in Fig. 3. However, the recent detection of H2S in the near-IR (Irwin et al., 2019b) has complicated this picture, showing an enhanced H2S relative humidity in the south polar region above the clouds, potentially associated with the warmer temperatures of the polar vortex, with local microphysical effects at the cloud tops, or with an absence of aerosols permitting longer path lengths through the polar atmosphere. A reconciliation of this cloud-top polar enhancement with the stark polar depletion observed in the microwave has yet to be performed.

4 Stacked Circulation Cells

4.1 Tropospheric Circulation

Fig. 5 updates our previous schematics of the upper-tropospheric circulation (Figs. 3-4) to incorporate the findings from the CH4 and H2S distributions. This follows similar attempts to construct two-dimensional tropospheric circulations by previous authors (Hofstadter and Butler, 2003; de Pater et al., 2014; Sromovsky et al., 2014; Tollefson et al., 2019). The equatorial upwelling and $p\sim 1$ bar cooling required to produce the windshear measured by Tollefson et al. (2018) is included222We note that the situation is more complex than this, as a cold equator and warm mid-latitudes are only needed if methane (and hence atmospheric density) are meridionally uniform. However, the results of Tollefson et al. (2018) are somewhat degenerate, as the inclusion of an equator-to-pole depletion in methane by a factor of four allows for the warm equator and cool mid-latitudes that are actually observed, consistent with Neptune’s “upper cell” of air rising at mid-latitudes and sinking over the equator. This study makes it clear that both the temperature gradients and density gradients should be accounted for when trying to interpret vertical wind shear., but new horizontal arrows depict the general equator-to-pole transport of both methane and H2S, and the strong polar depletions are depicted by downward arrows to the bottom of our domain (50-100 bars). While not indicated in the figure, it should be kept in mind that the polar region of depletion is much larger on Uranus, extending down to $\pm 45^{\circ}$ in both hemispheres, than it is on Neptune, where it extends to $\sim 65^{\circ}$S (the only hemisphere currently observed). The $p>1$-bar equatorial upwelling meets the $p<1$-bar equatorial subsidence somewhere near the methane-ice cloud tops, where a meridional divergence towards the poles would transport material to higher, off-equatorial latitudes. This could be responsible for the local minima in CH4 and H2S at the equator (Sromovsky et al., 2019; Irwin et al., 2019b; Tollefson et al., 2019), flanked by equatorial maxima where the upper tropospheric subsidence is weaker.

At mid-latitudes we retain the upwelling inferred from the temperature and para-H2 distributions for $p<1$ bar, which may also be responsible for mid-latitude peaks in Neptune’s H2S abundance333We caution the reader that these observations measure column abundances above the clouds, so if the cloud tops are lower then the column abundance will increase without any need to increase the H2S abundance. at these altitudes (Irwin et al., 2019b) and the formation of the polar collars on Uranus. The equator-to-pole gradients in methane and H2S mixing ratios could have a significant effect on the zonal windshear through density changes (Tollefson et al., 2018), which is quite unlike anything found on the Gas Giants, and the decrease of density with latitude should supplement the observed thermal wind shear, opposing the observed prograde jets as shown in Fig. 5. Such shears should persist to great depths, following the contrasts in volatile absorbers shown in Fig. 1F.

The mid-latitude subsidence at $p>1$ bar is more speculative, but could be responsible for the mid-latitude transitions (from strongly enriched to strongly depleted) observed in both the CH4 and H2S distributions in Fig. 1E-F. Furthermore, the $p>1$ bar circulation cell across the prograde jets in Fig. 5 are highly speculative, suggesting equatorward flow across the jet near 1 bar, balanced by poleward flow at some uncertain depth. This structure has been hypothesised by analogy to Jupiter’s prograde jets, which are forced by eddy-momentum fluxes converging and transporting energy into the jet streams (Ingersoll et al., 2004; Salyk et al., 2006). This convergence is balanced by a meridional circulation from the cyclonic belts to anticyclonic zones at the cloud-tops, and from the zones to the belts at depth (see the review by Fletcher et al., 2020). In this picture, the Ice Giant poles are analogous to jovian belts, the mid-latitudes to jovian zones, but we stress that there is currently no observational evidence for eddies forcing the prograde jets on the Ice Giants. This scenario is also qualitatively consistent with zones (cold and volatile-enriched) and belts (warm and volatile-depleted) in Jupiter’s upper troposphere at $p<1$ bar. Nevertheless, there are several problems with this conceptual picture. Firstly, the equatorward flow near $p\sim 1$ bar is in the opposite direction to the strong CH4 and H2S gradients. Secondly, it is unclear how any polar upwelling at depth (associated with the cross-jet flow) is mixed with the strong volatile depletions shown in the microwave observations. And thirdly, it remains to be seen how the circulation patterns modulate the fine-scale albedo bands described in Section 2.1. It is possible to envision a scenario where all of these processes are at work, but that one dominates the equator-to-pole gradients in composition, and the others contribute to local meteorological and cloud phenomena. Furthermore, this circulation pattern may ultimately be unnecessary if the equator-to-pole contrasts are shown to be driven by local microphysical and chemical processes, rather than large-scale transport.

By treating the polar domains as analogous to jovian belts, we might gain insights into a further conundrum: large-scale polar subsidence should be inhibiting cloud formation and small-scale convection, but we see evidence for both of these on Uranus. Jupiter’s cyclonic belts exhibit frequent small-scale moist convective activity over small areas, despite being locations of net subsidence in the upper troposphere. Indeed, maps of the distribution of lightning show a higher occurrence in cyclonic domains, and cyclonic features in particular (Little et al., 1999; Gierasch et al., 2000). Although lightning has been detected on both Uranus (Zarka and Pedersen, 1986) and Neptune (Gurnett et al., 1990), requiring the presence of polarisable, mixed-phase materials like water, the spatial distribution of Ice Giant lighting is unknown. Jupiter’s belts are strongly depleted in ammonia (Achterberg et al., 2006; Fletcher et al., 2016; de Pater et al., 2016; Li et al., 2017), in the same way as Ice Giant poles are strongly depleted in H2S.

But the exact mechanisms for moist convection in these domains of net subsidence and depleted volatiles remains unclear without some supply of condensable volatiles (and their latent heat) to the cloud base, possibly via the deep, poleward transport depicted in Fig. 5. The need for a deeper circulation cell, transporting H2S to high latitudes at the cloud level, was also noted by Sromovsky et al. (2014) as part of their multi-tiered structure for Uranus, and may be responsible for the enhanced H2S relative humidity at the cloud tops of Neptune observed by Irwin et al. (2019b). Maybe the strength of the polar convection changes with season, due to a destabilising effect as the troposphere warms in the spring, resulting in more convective transport of volatile-laden air at $p<1$ bar to generate the bright polar caps described in Section 2.1.

However, these multi-tiered structures in Fig. 5 are seemingly at odds with the single deep cells proposed for Neptune by de Pater et al. (2014) and Tollefson et al. (2019), and for Uranus by Molter et al. (2020), extending from the stratosphere down into the deep atmosphere. Maybe strong molecular weight gradients at the cloud condensation altitudes (Guillot, 1995) could serve to keep circulation cells separated in layers that have yet to be properly resolved by any of the remote-sensing investigations presented so far. Nevertheless, the multi-tiered stack of circulation cells provides a hypothesis to be tested by future measurements and numerical simulations.

4.2 Connection to the Interior

It is reasonable to ask how far down the atmospheric circulation in Fig. 5 extends - we currently assume that the bottom of this figure is at the 50-100 bar level, below the NH4SH clouds but still above the H2O cloud base. We have seen that these condensation zones for CH4 (and, to a lesser extent, H2S) may be accompanied by significant molecular weight gradients and alterations to the temperature lapse rate (Lindal et al., 1987; Guillot, 1995), which may serve to separate different layers or strata of circulation and convection (Gierasch and Conrath, 1987) and act as insulators for the escape of internal heat (e.g., Fortney et al., 2011). In particular, a stable layer at the depth of the water clouds could effectively separate the water-rich interior from the dry exterior of an Ice Giant. Episodic erosion of such a stable layer might be responsible for intermittent outbursts of storm activity or release of internal energy (Gierasch and Conrath, 1987; Li and Ingersoll, 2015; Friedson and Gonzales, 2017). And the potential stable layers within Uranus might be preventing it from following the expected thermal evolution, partially explaining its apparent absence of internal heat (Pearl and Conrath, 1991).

The deep water-ice clouds may be at the top of a massive aqueous water cloud extending to thousands of bars. These values are highly uncertain, given the lack of knowledge of the temperature lapse rate, particularly in the zone of water condensation where there might be significant departures from adiabatic behaviour (Guillot, 1995; Leconte et al., 2017). The bulk water abundance, inferred from the planet’s bulk density and using thermochemical modelling with disequilibrium trace species, may be enhanced by hundreds of times the solar abundance (Lodders and Fegley, 1994; Luszcz-Cook et al., 2013; Cavalié et al., 2017). Molecular dynamics calculations suggest that water may form a superionic icy ocean at great depths (Goldman, 2005; Millot et al., 2019), separated from the molecular envelope of H2 and removing (via dissolution) ammonia from the upper levels of the troposphere (de Pater et al., 1991). This fluid ionic “watery” ocean is likely the region that generates each planet’s internal dynamo (Ness et al., 1986, 1989). Intriguingly, dynamo simulations predict large circulation cells in the deep interior with upwelling (and a peak in the radial heat flux) near the equator and, depending on the thickness of the convecting layer, polar meridional cells may also be present (Soderlund et al., 2013). The connection between interior circulation and atmospheric circulation would be an intriguing subject for future work.

5 Stratospheric Circulation: Chemical Tracers

Having explored the potentially multi-tiered circulation cells in the upper troposphere and mid-troposphere, we now turn our attention to atmospheric motions above the tropopause. On the Gas Giants, the zonal winds and banded structures persist high into the stratosphere, and are superimposed onto larger interhemispheric circulations, both non-seasonal on Jupiter (Nixon et al., 2007; Zhang et al., 2013; Melin et al., 2018) and seasonal on Saturn (see review by Fletcher et al., 2018). The temperature and composition of the Gas Giant stratospheres is also modulated by the presence of polar vortices and wave-driven equatorial oscillations. In contrast, our understanding of spatial contrasts in the Ice Giant stratospheres is in its infancy due to their great distance from ground- and space-based telescopes, and the lack of a long-lived orbiter observing them from close range.

Stratospheric composition is driven by the UV photolysis of methane, resulting in a complicated mix of hydrocarbons (Atreya and Ponthieu, 1983; Moses et al., 2005; Dobrijevic et al., 2010; Moses et al., 2018) that can be investigated via mid-infrared remote sensing (Orton et al., 1987; Feuchtgruber et al., 1997; Encrenaz et al., 1998; Fletcher et al., 2010; Greathouse et al., 2011; Orton et al., 2014; Fletcher et al., 2014; Roman et al., 2020) and UV occultations (Herbert et al., 1987; Bishop et al., 1990). In addition, oxygen species like CO, CO2 and H2O are also present in the upper stratosphere, originating from cometary impacts, satellite debris, and ablation of interplanetary dust and ring particles. These compounds potentially play an important role in the photochemical pathways (Feuchtgruber et al., 1997; Lellouch et al., 2005, 2010; Luszcz-Cook et al., 2013; Cavalié et al., 2014; Orton et al., 2014; Poppe, 2016; Moses and Poppe, 2017). The abundances are sensitive to the strength of atmospheric mixing, so Uranus’ sluggish motions produces a lower methane homopause than on Neptune (Herbert et al., 1987; Bishop et al., 1990), ensuring that photochemistry on Uranus occurs in a different physical regime (higher pressures) than on any other giant planet, suppressing photochemical networks (Atreya et al., 1991). The effects of this are apparent from our observations: Uranus’ hydrocarbons are confined to altitudes below the 0.1-mbar level, and the ratio of ethane to acetylene is very different on Uranus compared to all the other giants (Orton et al., 2014). Furthermore, photolysis of CO and CO2 can lead to secondary peaks of hydrocarbon production at higher altitudes (Moses et al., 2018). The spatial distribution of these stratospheric species control the local radiative balance (ethane and acetylene are excellent coolants, but their efficiency leads to a stratospheric energy crisis, e.g., Li et al., 2018) and the condensation of thin stratospheric haze layers in the 0.1-30 mbar range (Rages et al., 1991; Romani et al., 1993; Moses and Poppe, 2017; Toledo et al., 2018). Aerosol layers of water, benzene, CO2, acetylene, ethane and propane are just some of the various condensed layers that might be expected at these low temperatures. As they sediment downwards into the troposphere, they can also modify the optical properties (or serve as nucleation sites for) the tropospheric aerosols (see Section 2.1). Thus any redistribution of the hydrocarbons via atmospheric circulation would have important implications for the energetics and hazes of the stratosphere.

The terrestrial stratosphere exhibits a wave-driven ‘Brewer-Dobson’ circulation (BDC, Andrews et al., 1987), transporting air (and ozone) from the equator to the pole. The circulation, with air rising at low latitudes and descending at mid- to high-latitudes, is driven primarily by Rossby (planetary) waves from the troposphere. At even higher altitudes, an upward motion in the summer hemisphere and downward motion in the winter hemisphere is known as the solsticial mesospheric circulation, and is primarily driven by gravity waves. Wave propagation and breaking on the Ice Giants may be a key mechanism for energy transport to partially resolve the stratospheric and thermospheric energy crises, where solar heating alone is insufficient to explain the high temperatures (Herbert et al., 1987; Stevens et al., 1993; Li et al., 2018; Melin et al., 2019). So wave-driven circulations might be at work in giant planet stratospheres, in addition to thermally-driven Hadley-like circulations such as those modelled by Conrath et al. (1990), who predicted mid-latitude rising and equatorial subsidence on Uranus. For example, Friedson and Moses (2012) provided models for a seasonally reversing circulation in Saturn’s stratosphere with rising motion in the summer hemisphere and sinking motion in the winter hemisphere, and Guerlet et al. (2009) interpreted a local maximum in Saturn’s hydrocarbons at $25^{\circ}$N as evidence for the descending branch of a Hadley-like circulation reaching into the stratosphere.

Assessing these circulations requires spatially-resolved observations of stratospheric temperatures, composition, and hazes, but these are extremely challenging. The stratospheric distribution of methane is particularly uncertain - it could be uniformly mixed; it could be reaching the stratosphere via convective overshooting from storm systems at mid-latitudes; it could be enhanced at the equator due to the tropospheric circulation in Section 3.2; or it could be leaking through the warm polar vortices (where the cold-trap is less efficient, Orton et al., 2007). Distinguishing between these possibilities remains a considerable challenge. Ground-based observations in the thermal-infrared have revealed relatively uniform stratospheric temperatures and composition on Neptune (Bézard et al., 1991; Greathouse et al., 2011; Fletcher et al., 2014), but with suggestions of a mid-latitude minimum and a rise in emission over the warm summertime polar vortex (Hammel et al., 2006; Orton et al., 2007, 2012; de Pater et al., 2014; Fletcher et al., 2014) suggestive of polar subsidence within $\sim 30^{\circ}$ of the south pole.

The colder temperatures of Uranus have made similar measurements of the Uranian stratosphere more challenging, but imaging data near 13 $\mu$m (sensitive to stratospheric acetylene) have revealed a contrasting pattern of emission with a distinct equatorial minimum in emission that rises sharply poleward of $\sim 25^{\circ}$ in both hemispheres (Orton et al., 2018; Roman et al., 2020). The observed gradient is inconsistent with the smooth trend predicted by radiative and photochemical modelling (Moses et al., 2018), and is the complete opposite of the warm equator observed in the upper troposphere (Section 2.3). If the emission contrast is due a sharp latitudinal gradient in stratospheric temperatures, this implies an additional equator-to-pole circulation above the tropopause, with low-latitude upwelling and high-latitude subsidence, as depicted in Fig. 6. The strong mid-latitude temperature gradient would also have implications for the thermal wind balance (see Section 6). However, if Uranus’ stratospheric warming poleward of $\sim 25^{\circ}$ is caused by subsidence, it implies an extremely broad area of polar downwelling compared to that observed on Neptune. Furthermore, Orton et al. (2018) reported an absence of major meridional structure in the stratospheric temperatures unambiguously probed using the H2 quadrupole lines, albeit sensing slightly deeper than the acetylene emission. We also note that the stratospheric circulation in Fig. 6 might appear, at first glance, to be inconsistent with the distribution of para-H2 disequilibrium in Fig. 1D. However, Voyager IRIS spectroscopy only constrained tropospheric para-H2, with lower-stratospheric distributions being the result of a speculative relaxation to uniform priors.

An alternative explanation was proposed by Roman et al. (2020) and is depicted as a schematic in Fig. 7. They showed that the observed stratospheric emission may, in part, be due to a mid-latitude enhancement of acetylene rather than temperatures. Noting that the structure of zonally averaged stratospheric emission over most latitudes appeared to be the inverse of what is observed in the upper troposphere (i.e., the stratosphere was brightest where the troposphere was coldest), Roman et al. (2020) argued that the resulting stratospheric structure was likely an extension of the upper-tropospheric circulation. In this scenario, upwelling at mid-latitudes transports methane from the troposphere, through the cold trap, and into the stratosphere for subsequent photolysis to acetylene. A similar explanation had been previously proposed to explain inferences of enhanced stratospheric hydrocarbons found in UVS data by Yelle et al. (1989) and McMillan (1992). The weaker emission at the equator compared to the pole (which tentatively appears bright in both the stratosphere and troposphere) would imply that either very little acetylene survives the journey to the equator or that there is an asymmetry in the meridional transport that leads to relatively greater acetylene abundances and/or temperatures at the poles. Data are still sparse and details of this model remain to be worked out, but more stringent constraints on the stratospheric circulation of both Ice Giants are expected from the James Webb Space Telescope (Moses et al., 2018).

6 Unresolved Equatorial Winds

Before discussing the different options for probe entry locations, we make one final speculative modification to the schematics presented in Figs. 3, 5, 6 and 7. On Jupiter and Saturn the relationship between the atmospheric temperatures and winds, via the thermal wind equation, is well established: prograde zonal jets occur poleward of cool anticyclonic zones and equatorward of warm, cyclonic belts; retrograde jets occur poleward of warm, cyclonic belts and equatorward of cool anticyclonic zones. In all cases, the maximum temperature gradient (and hence windshear) is co-located with the peak of a zonal jet, either eastward or westwards. At high latitudes on Uranus and Neptune, this relationship also seems to be true, with the prograde jet separating the cyclonic polar domain from the anticyclonic mid-latitude domain. But at low latitudes, both Voyager/IRIS (Fig. 1B) and ground-based observations of the tropospheres (Flasar et al., 1987; Conrath et al., 1998; Fletcher et al., 2014; Orton et al., 2015), suggest a strong temperature gradient in the $15-30^{\circ}$ latitude range, spatially coincident with the strong gradient in hydrocarbons or temperatures in the stratosphere (Orton et al., 2018; Roman et al., 2020). This is partially supported by Flasar et al. (1987), who computed the vertical windshear on Uranus to show maximal shear in the $\pm 15^{\circ}$ latitude region.

Could this be associated with peaks in the retrograde flow that have yet to be resolved in the zonal winds shown in Fig. 1C? Could such an unresolved jet be separating the mid-latitude anticyclonic domain (upwelling) from the equatorial cyclonic domain (subsidence)? We depict this in the schematic in Fig. 8. The quality of the available zonal wind data to date (e.g., Tollefson et al., 2018) do not preclude this idea, and the albedo patterns do show a good deal of zonal banding at low latitudes that could be connected with unresolved zonal jets. The jets need not be symmetric about the equator, given that the methane and H2S distributions in Section 5 show some asymmetries in terms of upwelling and subsiding regions. The existence of such off-equatorial zonal wind maxima would be analogous to those on Jupiter (where prograde winds peak at the edges of the Equatorial Zone at $\pm 7^{\circ}$ latitude) and Saturn (e.g., García-Melendo et al., 2011). Future observations of cloud-tracked winds will be needed to identify whether the equatorial retrograde jet shows a splitting in this manner.

7 Conclusion: Where to target a probe?

Our insights into Ice Giant circulation patterns have progressed slowly over the three decades since the Voyager observations, but enhanced techniques at visible wavelengths (including the Hubble Space Telescope) have tracked clouds and provided spatial distributions of methane; improvements in thermal-infrared observations allow us to characterise upper tropospheric and stratospheric temperatures; and the development of millimetre and centimetre capabilities now allow us to study the distribution of important volatiles (H2S and NH3) in the atmosphere well below the top-most clouds. Taken together, the distributions of these species present a puzzling, and sometimes conflicting, picture of meridional circulations on Uranus and Neptune. Indeed, many of our inferences may be biased by the more comprehensive studies of Jupiter and Saturn. An orbital mission in the coming decades, capable of multi-spectral remote sensing, will be essential to further reveal these circulation patterns. In addition, an in situ atmospheric probe falling to the 10-bar level (or deeper) would provide essential ‘ground-truth’ for the remote measurements, providing vertical profiles of temperature, density, gaseous abundances, para-hydrogen, aerosol properties, and Doppler-tracked horizontal winds (Mousis et al., 2018). As we described in Section 1, knowledge of the local circulation and meteorology of a probe entry location will be essential for the interpretation of the compositional and atmospheric structure measurements.

The selection of a probe entry location will depend upon a multitude of factors, primarily orbital mechanics and the requirement to have direct communications between a probe and its carrier spacecraft (e.g., Mousis et al., 2018). Nevertheless, the circulation patterns inferred in the previous sections can tell us about the expected conditions in different atmospheric domains (Fig. 7-8). Note that in this section, we are discussing how probe data can advance our knowledge of condensable species and atmospheric circulation. Measurements of noble gases and their isotopes, which are critical to testing planetary formation models, can be made at any latitude without regard to atmospheric patterns and weather due to the non-reactive nature of those species.

• •

Equatorial Domain from $0^{\circ}$ to $15-30^{\circ}$N/S: A deep atmosphere enriched in methane and H2S (and possibly other volatiles such as NH3 and H2O) by rising motion from the hundred-bar level or deeper. This rising motion meets a region of atmospheric subsidence in the upper troposphere, characterised by warm temperatures, super-equilibrium para-H2 fractions, and small-scale modulations of the CH4, H2S, and aerosol albedo. A region of vertical convergence must exist somewhere in the CH4- and H2S cloud forming regions (1-4 bar), with air moving meridionally away from the equator towards mid-latitudes. The meridional temperature gradient for $p>1$ bar provides a windshear that makes retrograde winds stronger with altitude, whereas the temperature gradient for $p<1$ bar makes the retrograde winds decay with height in the upper troposphere. Above the tropopause we may transition into a region dominated by wave-driven circulation, where recent work finds suppressed stratospheric emission associated with either low-latitude upwelling and adiabatic cooling, or with an absence of stratospheric hydrocarbons (Roman et al., 2020). This multi-tier structure, with vertical convergence in the cloud-forming region leading to meridional divergence and poleward flow, might be analogous to cyclonic jovian belts: regions of net subsidence in the upper troposphere, but net upwelling (and associated lightning from moist convection) in the deeper troposphere. The equatorial domain might be bounded by off-equatorial maxima in the retrograde zonal jet, which are currently unresolved in the available cloud tracking measurements. A probe descending right at the equator might encounter local minima in volatile species in the upper troposphere, but will ultimately reach the deeper cell where the abundances of CH4 and H2S are at their maximum. Local meteorological features, such as Uranus’ equatorial wave phenomena (Sromovsky et al., 2015), should be avoided to mitigate the risks of encountering strong downdrafts like those encountered by the Galileo probe (Orton et al., 1998).

• •

Mid-Latitude Domain from $15-30^{\circ}$N/S to $60-75^{\circ}$N/S: A transitional domain where the abundance of CH4 and H2S declines dramatically with increasing latitude due to increasingly strong deep-atmosphere subsidence. This domain is situated between the low-latitude retrograde jet and the high-latitude prograde jet, both of which decay with increasing depth below the clouds, and decay with increasing altitude above the clouds. Cold temperatures, sub-equilibrium para-H2 fractions, sporadic convective storm activity, and ‘polar collars‘ of aerosols all imply upwelling motions in the upper troposphere. This suggests vertical divergence somewhere near the cloud level (between upper tropospheric upwelling, and deep tropospheric subsidence), and thus meridional convergence of the air from the equatorial and polar domains. In this regard, the Ice Giant mid-latitudes are analogous to anticyclonic jovian zones, which are regions of net rising in the upper troposphere. Jovian zones also exhibit equatorward flow at cloud level across their bounding prograde jets (due to the meridional circulation balancing eddy momentum flux converging on the prograde jets), although such motions have not been observed on the Ice Giants. Probes entering these mid-latitude domains might encounter localised storm activity in the upper-tropospheric upwelling, enhanced optical depths of clouds and hazes, but negligible vertical shear on the zonal winds.

• •

Polar domain poleward of $60-75^{\circ}$N/S: The polar domains, bounded by the high-latitude prograde jets, are most depleted in volatiles CH4, H2S, NH3, and perhaps H2O due to strong atmospheric subsidence. Thus a probe entering this location might only be capable of returning upper limits on key elemental abundances and isotopic ratios (we note that noble gas measurements, however, can be made in this region). However, a polar probe would also sample a unique region where small-scale convective activity persists in a region of net subsidence, potentially leading to enhanced humidity of H2S immediately above the clouds. This small-scale convective activity is reminiscent of the plumes and lightning observed in Jupiter’s cyclonic belts, and maybe the puffy clouds observed at high latitudes on Saturn. It would also sample warm polar vortices evident in the troposphere and stratosphere (although more spatially confined to the pole on Neptune than on Uranus), potentially accessing unique chemical domains not found elsewhere on the planet.

Given the primary goal of an entry probe would be the measurement of bulk chemical enrichment and isotopic ratios, targeting the peak abundances at low latitudes would be a sensible first step. However, avoidance of low-latitude meteorological features is key to sampling a ‘representative’ region of an Ice Giant, so an off-equatorial entry site might be optimal, before encountering the storm bands and strong upwelling of the mid-latitudes. A secondary probe, if available, could target the polar domain to provide an extreme counterpoint for the low-latitude measurement. Such a comparison would provide the much-needed ground truth for remote-sensing investigations, able to sample the vertical distributions of temperatures, clouds, aerosols and gaseous species with a far greater vertical resolution than could ever be achieved from orbit.

Spatial variations in temperatures, clouds, and composition are extremely challenging to monitor from 20 or 30 AU away, even with future facilities like the 30-40-m class observatories or the James Webb Space Telescope. Multi-spectral remote sensing of the probe entry site will be invaluable to diagnose the probe results, including (i) tracking of winds and cloud features; (ii) UV-visible-near-infrared spectroscopy to assess aerosol distributions with height; (iii) thermal emission from the infrared to the sub-millimetre to determine temperatures and upwelling/subsiding motions; and (iv) microwave observations to connect the cloud-top meteorology to circulation patterns at great depth. Carefully planned, but ultimately serendipitous, ground-based observations were able to provide some of this information for the Galileo probe (Orton et al., 1998), but any probe mission should carefully consider having these capabilities on an associated orbiter.

In seeking to consider the observational evidence for atmospheric circulation on the Ice Giants, we have hypothesised a multi-tiered structure of stacked circulation cells, with motions that are potentially in opposition to one another. Can reality really be this complex? The case for stacked cells has been postulated on Jupiter and Saturn for some time (Ingersoll et al., 2000; Showman and de Pater, 2005; Fletcher et al., 2011), and for Uranus (Sromovsky et al., 2014) and Neptune (Tollefson et al., 2018), and was reviewed by Fletcher et al. (2020). The terrestrial atmosphere exhibits a transition from the troposphere (thermally-direct Hadley cells and mid-latitude heat transport by eddies) to the stratosphere (wave-driven thermally indirect circulation), with an associated reversal in the sense of the annual mean temperature gradient that might be considered as a multi-tiered circulation structure. The enrichments in volatiles in the Ice Giant atmospheres may generate substantial density stratifications where they condense, leading to circulation regimes that are only weakly connected to one another and a natural transition point between the stacked cells. However, there is an absence of numerical simulations of Ice Giant atmospheres (which must crucially include condensable species as active agents in the circulation) against which such inferences can be tested. We hope that any numerical simulations developed in the coming decade take all the observations presented in Fig. 1 into account as a test of their credibility, as well as future proposed observations from next-generation ground- and space-based observatories. This may well lead to a comprehensive rejection of the qualitative and complex multi-tiered circulations in Fig. 8, but would be invaluable in expanding and maturing our understanding of Ice Giant circulation patterns, and to guide our targeting of humankind’s first atmospheric probe for these distant and enigmatic worlds.

Acknowledgements.

Fletcher was supported by a Royal Society Research Fellowship at the University of Leicester. Fletcher and Roman also received funding from a European Research Council Consolidator Grant (under the European Union’s Horizon 2020 research and innovation programme, grant agreement No 723890). A portion of this work was performed by Orton at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA. de Pater was in part supported by the National Science Foundation, NSF grant AST-1615004 and NASA Grant NNX16AK14G through the Solar System Observations (SSO) program to the University of California, Berkeley. Irwin and Toledo acknowledge support from the UK Science and Technologies Facilities Council (STFC). We are grateful for support from the In Situ Exploration of the Ice Giants: Science and Technology workshop held in February 2019 at the Laboratoire d’Astrophysique de Marseille for inspiring this review.