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.

Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of satellites that might harbour sub-surface oceans). How can monitoring of auroral emission help inform future remote observations of emission from exoplanets? Our Solar System provides the only laboratory in which we can perform in-situ experiments to understand exoplanet formation, dynamos, systems and magnetospheres.

We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterising 133P in detail, solving the puzzle of the MBC's activity, and making the first in situ measurements of water in the asteroid belt. In many ways a successor to ESA's highly successful Rosetta mission, Castalia will allow direct comparison between very different classes of comet, including measuring critical isotope ratios, plasma and dust properties. It will also feature the first radar system to visit a minor body, mapping the ice in the interior. Castalia was proposed, in slightly different versions, to the ESA M4 and M5 calls within the Cosmic Vision programme. We describe the science motivation for the mission, the measurements required to achieve the scientific goals, and the proposed instrument payload and spacecraft to achieve these.

A leading explanation for the origin of Galactic cosmic rays is acceleration at high-Mach number shock waves in the collisionless plasma surrounding young supernova remnants. Evidence for this is provided by multi-wavelength non-thermal emission thought to be associated with ultrarelativistic electrons at these shocks. However, the dependence of the electron acceleration process on the orientation of the upstream magnetic field with respect to the local normal to the shock front (quasi-parallel/quasi-perpendicular) is debated. Cassini spacecraft observations at Saturn's bow shock has revealed examples of electron acceleration under quasi-perpendicular conditions, and the first in situ evidence of electron acceleration at a quasi-parallel shock. Here we use Cassini data to make the first comparison between energy spectra of locally accelerated electrons under these differing upstream magnetic field regimes. We present data taken during a quasi-perpendicular shock crossing on 2008 March 8 and during a quasi-parallel shock crossing on 2007 February 3, highlighting that both were associated with electron acceleration to at least MeV energies. The magnetic signature of the quasi-perpendicular crossing has a relatively sharp upstream-downstream transition, and energetic electrons were detected close to the transition and immediately downstream. The magnetic transition at the quasi-parallel crossing is less clear, energetic electrons were encountered upstream and downstream, and the electron energy spectrum is harder above ~100 keV. We discuss whether the acceleration is consistent with diffusive shock acceleration theory in each case, and suggest that the quasi-parallel spectral break is due to an energy-dependent interaction between the electrons and short, large-amplitude magnetic structures.

Low frequency turbulence in Saturn's magnetosheath is investigated using in-situ measurements of the Cassini spacecraft. Focus is put on the magnetic energy spectra computed in the frequency range $\sim[10^{-4}, 1]$Hz. A set of 42 time intervals in the magnetosheath were analyzed and three main results that contrast with known features of solar wind turbulence are reported: 1) The magnetic energy spectra showed a $\sim f^{-1}$ scaling at MHD scales followed by an $\sim f^{-2.6}$ scaling at the sub-ion scales without forming the so-called inertial range; 2) The magnetic compressibility and the cross-correlation between the parallel component of the magnetic field and density fluctuations $ C(\delta n,\delta B_{||}) $ indicates the dominance of the compressible magnetosonic slow-like modes at MHD scales rather than the Alfv\'en mode; 3) Higher order statistics revealed a monofractal (resp. multifractal) behaviour of the turbulent flow downstream of a quasi-perpendicular (resp. quasi-parallel) shock at the sub-ion scales. Implications of these results on theoretical modeling of space plasma turbulence are discussed.

The leading explanation for the origin of galactic cosmic rays is particle acceleration at the shocks surrounding young supernova remnants (SNRs), although crucial aspects of the acceleration process are unclear. The similar collisionless plasma shocks frequently encountered by spacecraft in the solar wind are generally far weaker (lower Mach number) than these SNR shocks. However, the Cassini spacecraft has shown that the shock standing in the solar wind sunward of Saturn (Saturn's bow shock) can occasionally reach this high-Mach number astrophysical regime. In this regime Cassini has provided the first in situ evidence for electron acceleration under quasi-parallel upstream magnetic conditions. Here we present the full picture of suprathermal electrons at Saturn's bow shock revealed by Cassini. The downstream thermal electron distribution is resolved in all data taken by the low-energy electron detector (CAPS-ELS, <28 keV) during shock crossings, but the higher energy channels were at (or close to) background. The high-energy electron detector (MIMI-LEMMS, >18 keV) measured a suprathermal electron signature at 31 of 508 crossings, where typically only the lowest energy channels (<100 keV) were above background. We show that these results are consistent with theory in which the "injection" of thermal electrons into an acceleration process involves interaction with whistler waves at the shock front, and becomes possible for all upstream magnetic field orientations at high Mach numbers like those of the strong shocks around young SNRs. A future dedicated study will analyze the rare crossings with evidence for relativistic electrons (up to ~1 MeV).

Shock waves exist throughout the universe and are fundamental to understanding the nature of collisionless plasmas. Reformation is a process, driven by microphysics, which typically occurs at high Mach number supercritical shocks. While ongoing studies have investigated this process extensively both theoretically and via simulations, their observations remain few and far between. In this letter we present a study of very high Mach number shocks in a parameter space that has been poorly explored and we identify reformation using in situ magnetic field observations from the Cassini spacecraft at 10 AU. This has given us an insight into quasi-perpendicular shocks across two orders of magnitude in Alfven Mach number (MA) which could potentially bridge the gap between modest terrestrial shocks and more exotic astrophysical shocks. For the first time, we show evidence for cyclic reformation controlled by specular ion reflection occurring at the predicted timescale of ~0.3 {\tau}c, where {\tau}c is the ion gyroperiod. In addition, we experimentally reveal the relationship between reformation and MA and focus on the magnetic structure of such shocks to further show that for the same MA, a reforming shock exhibits stronger magnetic field amplification than a shock that is not reforming.

A planet's magnetosheath extends from downstream of its bow shock up to the magnetopause where the solar wind flow is deflected around the magnetosphere and the solar wind embedded magnetic field lines are draped. This makes the region an important site for plasma turbulence, instabilities, reconnection and plasma depletion layers. A relatively high Alfv\'en Mach number solar wind and a polar-flattened magnetosphere make the magnetosheath of Saturn both physically and geometrically distinct from the Earth's. The polar flattening is predicted to affect the magnetosheath magnetic field structure and thus the solar wind-magnetosphere interaction. Here we investigate the magnetic field in the magnetosheath with the expectation that polar flattening is manifested in the overall draping pattern. We compare an accumulation of Cassini data between 2004 and 2010 with global magnetohydrodynamic (MHD) simulations and an analytical model representative of a draped field between axisymmetric boundaries. The draping patterns measured are well captured and in broad agreement for given upstream conditions with those of the MHD simulations (which include polar flattening). The deviations from the analytical model, based on no polar flattening, suggest that non-axisymmetry is invariably a key feature of the magnetosphere's global structure. Our results show a comprehensive overview of the configuration of the magnetic field in a non-axisymmetric magnetosheath as revealed by Cassini. We anticipate our assessment to provide an insight to this barely studied interface between a high Alfv\'enic bow shock and a dynamic magnetosphere.

Plasma rotation plays a large role in determining the size and shape of Saturn's disc-like magnetosphere. A magnetosphere more confined to the equator in the polar regions is expected as a result of the interaction between this type of obstacle and the solar wind. In addition, at times away from equinox, a north-south asymmetry is expected where the magnetopause will be further confined in one hemisphere but less confined in the opposite hemisphere. Examining the extent of this confinement has been limited by a lack of high-latitude spacecraft observations. Here, for the first time, direct evidence for polar confinement of Saturn's magnetopause has been observed using in-situ data obtained by the Cassini spacecraft during a series of high-inclination orbits between 2007 and 2009. Following techniques established by previous authors, we assume an equilibrium between the solar wind dynamic pressure (which Cassini is generally unable to measure directly), and the magnetic plus plasma pressure inside the magnetosphere. This assumption thus allows us to estimate the upstream solar wind dynamic pressure ($D_P$) for a series of magnetopause crossings, and hence to determine the expected location and global shape of the magnetopause as a function of $D_P$. A clear divergence from the familiar axisymmetric models of the magnetosphere is observed which may be characterised by an 'apparent flattening parameter' of 0.81+0.03/-0.06 (representing a simple dilation of the nominal axisymmetric boundary along the $Z_{KSM}$ axis such that the extent is reduced by 19% in this direction). This figure is insensitive to variations in $D_P$.

Collisionless shocks are loosely defined as shocks where the transition between pre-and post-shock states happens on a length scale much shorter than the collisional mean free path. In the absence of collision to enforce thermal equilibrium post-shock, electrons and ions need not have the same temperatures. While the acceleration of electrons for injection into shock acceleration processes to produce cosmic rays has received considerable attention, the related problem of the shock heating of quasi-thermal electrons has been relatively neglected. In this paper we review that state of our knowledge of electron heating in astrophysical shocks, mainly associated with supernova remnants (SNRs), shocks in the solar wind associated with the terrestrial and Saturnian bowshocks, and galaxy cluster shocks. The solar wind and SNR samples indicate that the ratio of electron temperature to ion temperature declines with increasing shock speed or Alfvenic Mach number. We discuss the extent to which such behavior can be understood via cosmic ray-generated waves in a shock precursor, which then subsequently damp by heating electrons. Finally, we speculate that a similar mechanism may be at work for both solar wind and SNR shocks.

Electrons can be accelerated to ultrarelativistic energies at strong (high-Mach number) collisionless shock waves that form when stellar debris rapidly expands after a supernova. Collisionless shock waves also form in the flow of particles from the Sun (the solar wind), and extensive spacecraft observations have established that electron acceleration at these shocks is effectively absent whenever the upstream magnetic field is roughly parallel to the shock surface normal (quasi-parallel conditions). However, it is unclear whether this magnetic dependence of electron acceleration also applies to the far stronger shocks around young supernova remnants, where local magnetic conditions are poorly understood. Here we present Cassini spacecraft observations of an unusually strong solar system shock wave (Saturn's bow shock) where significant local electron acceleration has been confirmed under quasi-parallel magnetic conditions for the first time, contradicting the established magnetic dependence of electron acceleration at solar system shocks. Furthermore, the acceleration led to electrons at relativistic energies (~MeV), comparable to the highest energies ever attributed to shock-acceleration in the solar wind. These observations demonstrate that at high-Mach numbers, like those of young supernova remnant shocks, quasi-parallel shocks become considerably more effective electron accelerators.