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

At wavelengths between 0.15 and 0.19 $\mu$m, the far-ultraviolet spectrum of Jupiter is dominated by the scattered solar spectrum, attenuated by molecular absorptions primarily by acetylene and ethane, and to a lesser extent ammonia and phosphine. We describe the development of our radiative transfer code that enables the retrieval of abundances of these molecular species from ultraviolet reflectance spectra. As a proof-of-concept we present an analysis of Cassini Ultraviolet Imaging Spectrograph (UVIS) observations of the disk of Jupiter during the 2000/2001 flyby. The ultraviolet-retrieved acetylene abundances in the upper stratosphere are lower than those predicted by models based solely on infrared thermal emission from the mid-stratosphere observed by the Composite Infrared Spectrometer (CIRS), requiring an adjustment to the vertical profiles above 1 mbar. We produce a vertical acetylene abundance profile that is compatible with both CIRS and UVIS, with reduced abundances at pressures $<$1 mbar: the 0.1 mbar abundances are $1.21 \pm 0.07$ ppm for acetylene and $20.8 \pm 5.1$ ppm for ethane. Finally, we perform a sensitivity study for the JUICE UVS instrument, which has extended wavelength coverage out to 0.21 $\mu$m, enabling the retrieval of ammonia and phosphine abundances, in addition to acetylene and ethane.

We present multi-wavelength measurements of the thermal, chemical, and cloud contrasts associated with the visibly dark formations (also known as 5-$\mu$m hot spots) and intervening bright plumes on the boundary between Jupiter's Equatorial Zone (EZ) and North Equatorial Belt (NEB). Observations made by the TEXES 5-20 $\mu$m spectrometer at the Gemini North Telescope in March 2017 reveal the upper-tropospheric properties of 12 hot spots, which are directly compared to measurements by Juno using the Microwave Radiometer (MWR), JIRAM at 5 $\mu$m, and JunoCam visible images. MWR and thermal-infrared spectroscopic results are consistent near 0.7 bar. Mid-infrared-derived aerosol opacity is consistent with that inferred from visible-albedo and 5-$\mu$m opacity maps. Aerosol contrasts, the defining characteristics of the cloudy plumes and aerosol-depleted hot spots, are not a good proxy for microwave brightness. The hot spots are neither uniformly warmer nor ammonia-depleted compared to their surroundings at $p<1$ bar. At 0.7 bar, the microwave brightness at the edges of hot spots is comparable to other features within the NEB. Conversely, hot spots are brighter at 1.5 bar, signifying either warm temperatures and/or depleted NH$_3$ at depth. Temperatures and ammonia are spatially variable within the hot spots, so the precise location of the observations matters to their interpretation. Reflective plumes sometimes have enhanced NH$_3$, cold temperatures, and elevated aerosol opacity, but each plume appears different. Neither plumes nor hot spots had microwave signatures in channels sensing $p>10$ bars, suggesting that the hot-spot/plume wave is a relatively shallow feature.

Jupiter's banded structure undergoes strong temporal variations, changing the visible and infrared appearance of the belts and zones in a complex and turbulent way due to physical processes that are not yet understood. In this study we use ground-based 5-$\mu$m infrared data captured between 1984 and 2018 by 8 different instruments mounted on the Infrared Telescope Facility in Hawai'i and on the Very Large Telescope in Chile to analyze and characterize the long-term variability of Jupiter's cloud-forming region at the 1-4 bar pressure level. The data show a large temporal variability mainly at the equatorial and tropical latitudes, with a smaller temporal variability at mid-latitudes. We also compare the 5-$\mu$m-bright and -dark regions with the locations of the visible zones and belts and we find that these regions are not always co-located, specially in the southern hemisphere. We also present Lomb-Scargle and Wavelet Transform analyzes in order to look for possible periodicities of the brightness changes that could help us understand their origin and predict future events. We see that some of these variations occur periodically in time intervals of 4-8 years. The reasons of these time intervals are not understood and we explore potential connections to both convective processes in the deeper weather layer and dynamical processes in the upper troposphere and stratosphere. Finally we perform a Principal Component analysis to reveal a clear anticorrelation on the 5-$\mu$m brightness changes between the North Equatorial Belt and the South Equatorial Belt, suggesting a possible connection between the changes in these belts.

Since its discovery in the aurorae of Jupiter ~30 years ago, the H$_{3}^{+}$ ion has served as an invaluable probe of giant planet upper atmospheres. However, the vast majority of monitoring of planetary H$_{3}^{+}$ radiation has followed from observations that rely on deriving parameters from column-integrated paths through the emitting layer. Here, we investigate the effects of density and temperature gradients along such paths on the measured H$_{3}^{+}$ spectrum and its resulting interpretation. In a non-isothermal atmosphere, H$_{3}^{+}$ column densities retrieved from such observations are found to represent a lower limit, reduced by 20% or more from the true atmospheric value. Global simulations of Uranus' ionosphere reveal that measured H$_{3}^{+}$ temperature variations are often attributable to well-understood solar zenith angle effects rather than indications of real atmospheric variability. Finally, based on these insights, a preliminary method of deriving vertical temperature structure is demonstrated at Jupiter using model reproductions of electron density and H$_{3}^{+}$ measurements. The sheer diversity and uncertainty of conditions in planetary atmospheres prohibits this work from providing blanket quantitative correction factors; nonetheless, we illustrate a few simple ways in which the already formidable utility of H$_{3}^{+}$ observations in understanding planetary atmospheres can be enhanced.

We characterise the origin and evolution of a mesoscale wave pattern in Jupiter's North Equatorial Belt (NEB), detected for the first time at 5 $\mu$m using a 2016-17 campaign of `lucky imaging' from the VISIR instrument on the Very Large Telescope and the NIRI instrument on the Gemini observatory, coupled with M-band imaging from Juno's JIRAM instrument during the first seven Juno orbits. The wave is compact, with a $1.1-1.4^\circ$ longitude wavelength (wavelength 1,300-1,600 km, wavenumber 260-330) that is stable over time, with wave crests aligned largely north-south between $14$ and $17^\circ$N (planetographic). The waves were initially identified in small ($10^\circ$ longitude) packets immediately west of cyclones in the NEB at $16^\circ$N, but extended to span wider longitude ranges over time. The waves exhibit a 7-10 K brightness temperature amplitude on top of a $\sim210$-K background at 5 $\mu$m. The thermal structure of the NEB allows for both inertio-gravity waves and gravity waves. Despite detection at 5 $\mu$m, this does not necessarily imply a deep location for the waves, and an upper tropospheric aerosol layer near 400-800 mbar could feature a gravity wave pattern modulating the visible-light reflectivity and attenuating the 5-$\mu$m radiance originating from deeper levels. Strong rifting activity appears to obliterate the pattern, which can change on timescales of weeks. The NEB underwent a new expansion and contraction episode in 2016-17 with associated cyclone-anticyclone formation, which could explain why the mesoscale wave pattern was more vivid in 2017 than ever before.

Acetylene (C$_2$H$_2$) and ethane (C$_2$H$_6$) are both produced in the stratosphere of Jupiter via photolysis of methane (CH$_4$). Despite this common source, the latitudinal distribution of the two species is radically different, with acetylene decreasing in abundance towards the pole, and ethane increasing towards the pole. We present six years of NASA IRTF TEXES mid-infrared observations of the zonally-averaged emission of methane, acetylene and ethane. We confirm that the latitudinal distributions of ethane and acetylene are decoupled, and that this is a persistent feature over multiple years. The acetylene distribution falls off towards the pole, peaking at $\sim$30$^{\circ}$N with a volume mixing ratio (VMR) of $\sim$0.8 parts per million (ppm) at 1 mbar and still falling off at $\pm70^\circ$ with a VMR of $\sim$0.3 ppm. The acetylene distributions are asymmetric on average, but as we move from 2013 to 2017, the zonally-averaged abundance becomes more symmetric about the equator. We suggest that both the short term changes in acetylene and its latitudinal asymmetry is driven by changes to the vertical stratospheric mixing, potentially related to propagating wave phenomena. Unlike acetylene, ethane has a symmetric distribution about the equator that increases toward the pole, with a peak mole fraction of $\sim$18 ppm at about $\pm50^{\circ}$ latitude, with a minimum at the equator of $\sim$10 ppm at 1 mbar. [...] The equator-to-pole distributions of acetylene and ethane are consistent with acetylene having a shorter lifetime than ethane that is not sensitive to longer advective timescales, but is augmented by short-term dynamics, such as vertical mixing. Conversely, the long lifetime of ethane allows it to be transported to higher latitudes faster than it can be chemically depleted.

Emission from the molecular ion H$_3^+$ is a powerful diagnostic of the upper atmosphere of Jupiter, Saturn, and Uranus, but it remains undetected at Neptune. In search of this emission, we present near-infrared spectral observations of Neptune between 3.93 and 4.00 $\mu$m taken with the newly commissioned iSHELL instrument on the NASA Infrared Telescope Facility in Hawaii, obtained 17-20 August 2017. We spent 15.4 h integrating across the disk of the planet, yet were unable to unambiguously identify any H$_3^+$ line emissions. Assuming a temperature of 550 K, we derive an upper limit on the column integrated density of $1.0^{+1.2}_{-0.8}\times10^{13}$ m$^{-2}$, which is an improvement of 30\% on the best previous observational constraint. This result means that models are over-estimating the density by at least a factor of 5, highlighting the need for renewed modelling efforts. A potential solution is strong vertical mixing of polyatomic neutral species from Neptune's upper stratosphere to the thermosphere, reacting with H$_3^+$, thus greatly reducing the column integrated H$_3^+$ densities. This upper limit also provide constraints on future attempts at detecting H$_3^+$ using the James Webb Space Telescope.

The dark colors of Jupiter's North Equatorial Belt (NEB, $7-17^\circ$N) appeared to expand northward into the neighboring zone in 2015, consistent with a 3-5 year cycle of activity in the NEB. Inversions of thermal-IR imaging from the Very Large Telescope revealed a moderate warming and reduction of aerosol opacity at the cloud tops at $17-20^\circ$N, suggesting subsidence and drying in the expanded sector. Two new thermal waves were identified during this period: (i) an upper tropospheric thermal wave (wavenumber 16-17, amplitude 2.5 K at 170 mbar) in the mid-NEB that was anti-correlated with haze reflectivity; and (ii) a stratospheric wave (wavenumber 13-14, amplitude 7.3 K at 5 mbar) at $20-30^\circ$N. Both were quasi-stationary, confined to regions of eastward zonal flow, and are morphologically similar to waves observed during previous expansion events.

The aurorae of Uranus were recently detected in the far ultraviolet with the Hubble Space Telescope (HST) providing a new, so far unique, means to remotely study the asymmetric Uranian magnetosphere from Earth. We analyze here two new HST Uranus campaigns executed in Sept. 2012 and Nov. 2014 with different temporal coverage and under variable solar wind conditions numerically predicted by three different MHD codes. Overall, HST images taken with the Space Telescope Imaging Spectrograph reveal auroral emissions in three pairs of successive images (one pair acquired in 2012 and two in 2014), hence six additional auroral detections in total, including the most intense Uranian aurorae ever seen with HST. The detected emissions occur close the expected arrival of interplanetary shocks. They appear as extended spots at southern latitudes, rotating with the planet. They radiate 5-24 kR and 1.3-8.8 GW of ultraviolet emission from H$_2$, last for tens of minutes and vary on timescales down to a few seconds. Fitting the 2014 observations with model auroral ovals constrains the longitude of the southern (northern) magnetic pole to 104+\-26{\deg} (284-\-26{\deg}) in the Uranian Longitude System. We suggest that the Uranian near-equinoctial aurorae are pulsed cusp emissions possibly triggered by large-scale magnetospheric compressions.

We use high-resolution ground-based observations from the VLT CRIRES instrument in December 2012 to identify sixteen previously undetected H3+ emission lines from Jupiter's ionosphere. These emission lines are located in Jupiter's 5-micron window (4.5-5.2 {\mu}m), an optically-thin region of the planet's spectrum where the radiation mostly originates from the deep troposphere. The H3+ emission lines are so strong that they are visible even against this bright background. We measure the Doppler broadening of the H3+ emission lines in order to evaluate the kinetic temperature of the molecules, and we obtain a value of 1390$\pm$160 K. We also measure the relative intensities of lines in the {\nu}2 fundamental in order to calculate the rotational temperature, obtaining a value of 960$\pm$40 K. Finally, we use the detection of an emission line from the 2{\nu}2(2)-{\nu}2 overtone to measure a vibrational temperature of 925$\pm$25 K. We use these three independent temperature estimates to discuss the thermodynamic equilibrium of Jupiter's ionosphere.

We present an analysis of recent high spatial and spectral resolution ground-based infrared observations of H3+ spectra obtained with the 10-metre Keck II telescope in April 2011. We observed H3+ emission from Saturn's northern and southern auroral regions, simultaneously, over the course of more than two hours, obtaining spectral images along the central meridian as Saturn rotates. Previous ground-based work has derived only an average temperature of an individual polar region, summing an entire night of observations. Here we analyse 20 H3+ spectra, 10 for each hemisphere, providing H3+ temperature, column density and total emission in both the northern and southern polar regions simultaneously, improving on past results in temporal cadence and simultaneity. We find that: 1) the average thermospheric temperatures are 527+/-18 K in northern Spring and 583+/-13 K in southern Autumn, respectively; 2) this asymmetry in temperature is likely to be the result of an inversely proportional relationship between the total thermospheric heating rate (Joule heating and ion drag) and magnetic field strength - i.e. the larger northern field strength leads to reduced total heating rate and a reduced temperature, irrespective of season, and 3) this implies that thermospheric heating and temperatures are relatively insensitive to seasonal effects.

Saturn's ionosphere is produced when the otherwise neutral atmosphere is exposed to a flow of energetic charged particles or solar radiation. At low latitudes the latter should result in a weak planet-wide glow in infrared (IR), corresponding to the planet's uniform illumination by the Sun. The observed low-latitude ionospheric electron density is lower and the temperature higher than predicted by models. A planet-ring magnetic connection has been previously suggested in which an influx of water from the rings could explain the lower than expected electron densities in Saturn's atmosphere. Here we report the detection of a pattern of features, extending across a broad latitude band from ~25 to 60 degrees, that is superposed on the lower latitude background glow, with peaks in emission that map along the planet's magnetic field lines to gaps in Saturn's rings. This pattern implies the transfer of charged water products from the ring-plane to the ionosphere, revealing the influx on a global scale, flooding between 30 to 43% of the planet's upper-atmospheric surface. This ring `rain' plays a fundamental role in modulating ionospheric emissions and suppressing electron densities.

From the 27th to the 28th January 2009, the Cassini spacecraft remotely acquired combined observations of Saturns southern aurorae at radio, ultraviolet and infrared wavelengths, while monitoring ion injections in the middle magnetosphere from energetic neutral atoms. Simultaneous measurements included the sampling of a full planetary rotation, a relevant timescale to investigate auroral emissions driven by processes internal to the magnetosphere. In addition, this interval coincidently matched a powerful substorm-like event in the magnetotail, which induced an overall dawnside intensification of the magnetospheric and auroral activity. We comparatively analyze this unique set of measurements to reach a comprehensive view of kronian auroral processes over the investigated timescale. We identify three source regions in atmospheric aurorae, including a main oval associated with the bulk of Saturn Kilometric Radiation (SKR), together with polar and equatorward emissions. These observations reveal the co-existence of corotational and sub-corototational dynamics of emissions associated with the main auroral oval. Precisely, we show that the atmospheric main oval hosts short-lived sub-corotating isolated features together with a bright, longitudinally extended, corotating region locked at the southern SKR phase. We assign the susbtorm-like event to a regular, internally driven, nightside ion injection possibly triggerred by a plasmoid ejection. We also investigate the total auroral energy budget, from the power input to the atmosphere, characterized by precipitating electrons up to 20 keV, to its dissipation through the various radiating processes. Finally, through simulations, we confirm the search-light nature of the SKR rotational modulation and we show that SKR arcs relate to isolated auroral spots. The resulting findings are discussed in the frame of pending questions.