The geology of Pluto and Charon through the eyes of New Horizons

Chapter 7: Assessing Habitability Beyond Earth

All known life on Earth inhabits environments that maintain conditions between certain extremes of temperature, chemical composition, energy availability, and so on (Chapter 6). Life may have emerged in similar environments elsewhere in the Solar System and beyond. The ongoing search for life elsewhere mainly focuses on those environments most likely to support life, now or in the past-that is, potentially habitable environments. Discussion of habitability is necessarily based on what we know about life on Earth, as it is our only example. This chapter gives an overview of the known and presumed requirements for life on Earth and discusses how these requirements can be used to assess the potential habitability of planetary bodies across the Solar System and beyond. We first consider the chemical requirements of life and potential feedback effects that the presence of life can have on habitable conditions, and then the planetary, stellar, and temporal requirements for habitability. We then review the state of knowledge on the potential habitability of bodies across the Solar System and exoplanets, with a particular focus on Mars, Venus, Europa, and Enceladus. While reviewing the case for the potential habitability of each body, we summarize the most prominent and impactful studies that have informed the perspective on where habitable environments are likely to be found.

Chapter 1: The Astrobiology Primer 3.0

The Astrobiology Primer 3.0 (ABP3.0) is a concise introduction to the field of astrobiology for students and others who are new to the field of astrobiology. It provides an entry into the broader materials in this supplementary issue of Astrobiology and an overview of the investigations and driving hypotheses that make up this interdisciplinary field. The content of this chapter was adapted from the other 10 articles in this supplementary issue and thus represents the contribution of all the authors who worked on these introductory articles. The content of this chapter is not exhaustive and represents the topics that the authors found to be the most important and compelling in a dynamic and changing field.

Pluto’s Surface Mapping Using Unsupervised Learning from Near-infrared Observations of LEISA/Ralph

We map the surface of Pluto using an unsupervised machine-learning technique using the near-infrared observations of the LEISA/Ralph instrument on board NASA’s New Horizons spacecraft. The principal-component-reduced Gaussian mixture model was implemented to investigate the geographic distribution of the surface units across the dwarf planet. We also present the likelihood of each surface unit at the image pixel level. Average I/F spectra of each unit were analyzed—in terms of the position and strengths of absorption bands of abundant volatiles such as N2, CH4, and CO and nonvolatile H2O—to connect the unit to surface composition, geology, and geographic location. The distribution of surface units shows a latitudinal pattern with distinct surface compositions of volatiles—consistent with the existing literature. However, previous mapping efforts were based primarily on compositional analysis using spectral indices (indicators) or implementation of complex radiative transfer models, which need (prior) expert knowledge, label data, or optical constants of representative end-members. We prove that an application of unsupervised learning in this instance renders a satisfactory result in mapping the spatial distribution of ice compositions without any prior information or label data. Thus, such an application is specifically advantageous for a planetary surface mapping when label data are poorly constrained or completely unknown, because an understanding of surface material distribution is vital for volatile transport modeling at the planetary scale. We emphasize that the unsupervised learning used in this study has wide applicability and can be expanded to other planetary bodies of the solar system for mapping surface material distribution.

Endogenically sourced volatiles on Charon and other Kuiper belt objects

Kuiper belt objects (KBOs) have diverse surface compositions, and the New Horizons mission to the Pluto-Charon system allows us to test hypotheses on the origin and evolution of these KBO surfaces. Previous work proposed that Charon’s organic-rich north pole formed from radiolytically processed volatiles sourced from Pluto’s escaping atmosphere. Here, we show an endogenic source of volatiles from Charon’s interior is plausible. We calculate that cryovolcanic resurfacing released 1.29 × 10¹⁵–3.47 × 10¹⁵ kg of methane to Charon’s surface from its interior. We modeled volatile transport and found the vast majority of this volcanically released methane migrates to Charon’s poles, with deposition rates sufficient to be processed into the observed organic compounds. Irradiated methane products appear on similarly sized KBOs that do not orbit a Pluto-sized object to draw an escaping atmosphere from, so interior-sourced volatiles could be a common and important process across the Kuiper belt.

We show cryovolcanic eruptions released sufficient methane to source volatile products on Charon. Irradiated methane products are found on other Kuiper belt objects, so endogenically sourced volatiles could be important across the Kuiper belt.

Ab initio simulations of α- and β-ammonium carbamate (NH4·NH2CO2), and the thermal expansivity of deuterated α-ammonium carbamate from 4.2 to 180 K by neutron powder diffraction

Experimental and computational studies of ammonium carbamate have been carried out, with the objective of studying the elastic anisotropy of the framework manifested in (i) the thermal expansion and (ii) the compressibility; furthermore, the relative thermodynamic stability of the two known polymorphs has been evaluated computationally. Using high-resolution neutron powder diffraction data, the crystal structure of α-ammonium carbamate (ND₄·ND₂CO₂) has been refined [space group Pbca, Z = 8, with a = 17.05189 (15), b = 6.43531 (7), c = 6.68093 (7) Å and V = 733.126 (9) ų at 4.2 K] and the thermal expansivity of α-ammonium carbamate has been measured over the temperature range 4.2-180 K. The expansivity shows a high degree of anisotropy, with the b axis most expandable. The ab initio computational studies were carried out on the α- and β-polymorphs of ammonium carbamate using density functional theory. Fitting equations of state to the P(V) points of the simulations (run athermally) gave the following values: V₀ = 744 (2) ų and bulk modulus K₀ = 16.5 (4) GPa for the α-polymorph, and V₀ = 713.6 (5) ų and K₀ = 24.4 (4) GPa for the β-polymorph. The simulations show good agreement with the thermoelastic behaviour of α-ammonium carbamate. Both phases show a high-degree of anisotropy; in particular, α-ammonium carbamate shows unusual compressive behaviour, being determined to have negative linear compressibility (NLC) along its a axis above 5 GPa. The thermodynamically stable phase at ambient pressure is the α-polymorph, with a calculated enthalpy difference with respect to the β-polymorph of 0.399 kJ mol^-1; a transition to the β-polymorph could occur at ∼0.4 GPa.

Large-scale cryovolcanic resurfacing on Pluto

The New Horizons spacecraft returned images and compositional data showing that terrains on Pluto span a variety of ages, ranging from relatively ancient, heavily cratered areas to very young surfaces with few-to-no impact craters. One of the regions with very few impact craters is dominated by enormous rises with hummocky flanks. Similar features do not exist anywhere else in the imaged solar system. Here we analyze the geomorphology and composition of the features and conclude this region was resurfaced by cryovolcanic processes, of a type and scale so far unique to Pluto. Creation of this terrain requires multiple eruption sites and a large volume of material (>10⁴ km³) to form what we propose are multiple, several-km-high domes, some of which merge to form more complex planforms. The existence of these massive features suggests Pluto’s interior structure and evolution allows for either enhanced retention of heat or more heat overall than was anticipated before New Horizons, which permitted mobilization of water-ice-rich materials late in Pluto’s history.

Giant icy volcanos (cryovolcanos) on Pluto are unique in the imaged solar system and provide evidence for unexpected, active geology late in Pluto’s history.

A Geoscientific Review on CO and CO2 Ices in the Outer Solar System

Ground-based telescopes and space exploration have provided outstanding observations of the complexity of icy planetary surfaces. This work presents our review of the varying nature of carbon dioxide (CO2) and carbon monoxide (CO) ices from the cold traps on the Moon to Pluto in the Kuiper Belt. This review is organized into five parts. First, we review the mineral physics (e.g., rheology) relevant to these environments. Next, we review the radiation-induced chemical processes and the current interpretation of spectral signatures. The third section discusses the nature and distribution of CO2 in the giant planetary systems of Jupiter and Saturn, which are much better understood than the satellites of Uranus and Neptune, discussed in the subsequent section. The final sections focus on Pluto in comparison to Triton, having mainly CO, and a brief overview of cometary materials. We find that CO2 ices exist on many of these icy bodies by way of magnetospheric influence, while intermixing into solid ices with CH4 (methane) and N2 (nitrogen) out to Triton and Pluto. Such radiative mechanisms or intermixing can provide a wide diversity of icy surfaces, though we conclude where further experimental research of these ices is still needed.

Excited state photochemically driven surface formation of benzene from acetylene ices on Pluto and in the outer solar system

NASA's New Horizons mission unveiled a diverse landscape of Pluto's surface with massive regions being neutral in color, while others like Cthulhu Macula range from golden-yellow to reddish comprising up to half of Pluto's carbon budget. Here, we demonstrate in laboratory experiments merged with electronic structure calculations that the photolysis of solid acetylene - the most abundant precipitate on Pluto's surface - by low energy ultraviolet photons efficiently synthesizes benzene and polycyclic aromatic hydrocarbons via excited state photochemistry thus providing critical molecular building blocks for the colored surface material. Since low energy photons deliver doses to Pluto's surface exceeding those from cosmic rays by six orders of magnitude, these processes may significantly contribute to the coloration of Pluto's surface and of hydrocarbon-covered surfaces of Solar System bodies such as Triton in general. This discovery critically enhances our perception of the distribution of aromatic molecules and carbon throughout our Solar System.

Sublimation-driven convection in Sputnik Planitia on Pluto

Sputnik Planitia is a nitrogen-ice-filled basin on Pluto¹. Its polygonal surface patterns² have been previously explained as a result of solid-state convection with either an imposed heat flow³ or a temperature difference within the 10-km-thick ice layer⁴. Neither explanation is satisfactory, because they do not exhibit surface topography with the observed pattern: flat polygons delimited by narrow troughs⁵. Internal heating produces the observed patterns⁶, but the heating source in such a setup remains enigmatic. Here we report the results of modelling the effects of sublimation at the surface. We find that sublimation-driven convection readily produces the observed polygonal structures if we assume a smaller heat flux (~0.3 mW m^-2) at the base of the ice layer than the commonly accepted value of 2-3 mW m^-2 (ref. ⁷). Sustaining this regime with the latter value is also possible, but would require a stronger viscosity contrast (~3,000) than the nominal value (~100) considered in this study.

Topography and geology of Uranian mid-sized icy satellites in comparison with Saturnian and Plutonian satellites

Newly processed global imaging and topographic mapping of Uranus's five major satellites reveal differences and similarities to mid-sized satellites at Saturn and Pluto. Three modes of internal heat redistribution are recognized. The broad similarity of Miranda's three oval resurfacing zones to those mapped on Enceladus and (subtly) on Dione are likely due to antipodal diapiric upwelling. Conversely, break-up and foundering of crustal blocks accompanied by extensive (cryo)volcanism is the dominant mode on both Charon and Ariel. Titania's fault network finds parallels on Rhea, Dione, Tethys and possibly Oberon. Differences in the geologic style of resurfacing in the satellite systems (e.g. plains on Charon, Dione, Tethys and perhaps Titania versus ridges on Miranda and Ariel) may be driven by differences in ice composition. Surface processes such as volatile transport may also be indicated by bright and dark materials on Oberon, Umbriel and Charon. The more complete and higher quality observations of the Saturnian and Plutonian mid-sized icy satellites by Cassini and New Horizons reveal a wealth of features and phenomena that cannot be perceived in the more limited Voyager coverage of the Uranian satellites, harbingers of many discoveries awaiting us on a return to Uranus. This article is part of a discussion meeting issue 'Future exploration of ice giant systems'.

Equatorial mountains on Pluto are covered by methane frosts resulting from a unique atmospheric process

Pluto is covered by numerous deposits of methane, either diluted in nitrogen or as methane-rich ice. Within the dark equatorial region of Cthulhu, bright frost containing methane is observed coating crater rims and walls as well as mountain tops, providing spectacular resemblance to terrestrial snow-capped mountain chains. However, the origin of these deposits remained enigmatic. Here we report that they are composed of methane-rich ice. We use high-resolution numerical simulations of Pluto’s climate to show that the processes forming them are likely to be completely different to those forming high-altitude snowpack on Earth. The methane deposits may not result from adiabatic cooling in upwardly moving air like on our planet, but from a circulation-induced enrichment of gaseous methane a few kilometres above Pluto’s plains that favours methane condensation at mountain summits. This process could have shaped other methane reservoirs on Pluto and help explain the appearance of the bladed terrain of Tartarus Dorsa.

Pluto is covered by numerous deposits of methane. Here, the authors show that the formation of methane frost on mountain tops and crater rims in Pluto’s equatorial regions completely differ from those forming snow-capped mountains on Earth.

The crystallography of Pluto

Maynard-Casely and co-workers [IUCrJ (2020). 7, 844-851.] investigate two of Pluto's most abundant minerals with neutron diffraction. The new results will be key to understanding the geology of our distant neighbour and represent a significant advance in the emerging field of small-molecule geology.

Re-examining the crystal structure behaviour of nitrogen and methane

In the light of NASA's New Horizons mission, the solid-phase behaviour of methane and nitrogen has been re-examined and the thermal expansion coefficients of both materials have been determined over their whole solid temperature range for the first time. Neutron diffraction results indicate that the symmetric Pa 3 space group is the best description for the α-nitrogen structure, rather than the long-accepted P213. Furthermore, it is also observed that β-nitrogen and methane phase I show changes in texture on warming, indicating grain growth.

Prebiotic Chemistry of Pluto

We present the case for the presence of complex organic molecules, such as amino acids and nucleobases, formed by abiotic processes on the surface and in near-subsurface regions of Pluto. Pluto's surface is tinted with a range of non-ice substances with colors ranging from light yellow to red to dark brown; the colors match those of laboratory organic residues called tholins. Tholins are broadly characterized as complex, macromolecular organic solids consisting of a network of aromatic structures connected by aliphatic bridging units (e.g., Imanaka et al., 2004; Materese et al., 2014, 2015). The synthesis of tholins in planetary atmospheres and in surface ices has been explored in numerous laboratory experiments, and both gas- and solid-phase varieties are found on Pluto. A third variety of tholins, exposed at a site of tectonic surface fracturing called Virgil Fossae, appears to have come from a reservoir in the subsurface. Eruptions of tholin-laden liquid H₂O from a subsurface aqueous repository appear to have covered portions of Virgil Fossae and its surroundings with a uniquely colored deposit (D.P. Cruikshank, personal communication) that is geographically correlated with an exposure of H₂O ice that includes spectroscopically detected NH₃ (C.M. Dalle Ore, personal communication). The subsurface organic material could have been derived from presolar or solar nebula processes, or might have formed in situ. Photolysis and radiolysis of a mixture of ices relevant to Pluto's surface composition (N₂, CH₄, CO) have produced strongly colored, complex organics with a significant aromatic content having a high degree of nitrogen substitution similar to the aromatic heterocycles pyrimidine and purine (Materese et al., 2014, 2015; Cruikshank et al., 2016). Experiments with pyrimidines and purines frozen in H₂O-NH₃ ice resulted in the formation of numerous nucleobases, including the biologically relevant guanine, cytosine, adenine, uracil, and thymine (Materese et al., 2017). The red material associated with the H₂O ice may contain nucleobases resulting from energetic processing on Pluto's surface or in the interior. Some other Kuiper Belt objects also exhibit red colors similar to those found on Pluto and may therefore carry similar inventories of complex organic materials. The widespread and ubiquitous nature of similarly complex organic materials observed in a variety of astronomical settings drives the need for additional laboratory and modeling efforts to explain the origin and evolution of organic molecules. Pluto observations reveal complex organics on a small body that remains close to its place of origin in the outermost regions of the Solar System.

Ferrovolcanism: Iron Volcanism on Metallic Asteroids

Metallic asteroids, the exposed cores of disrupted planetesimals, are expected to have been exposed while still molten. Some would have cooled from the outside in, crystallizing a surface crust which would then grow inward. Because the growing crust is expected to be more dense than the underlying melt, this melt will tend to migrate toward the surface whenever it is able. Compressional stresses produced in the crust while it cools will be relieved locally by thrust faulting, which will also provide potential conduits for melt to reach the surface. We predict iron volcanism to have occurred on metallic asteroids as they cooled and discuss the implications of this process for both the evolution and the modern appearance of these bodies.

Detection of ammonia on Pluto's surface in a region of geologically recent tectonism

We report the detection of ammonia (NH3) on Pluto's surface in spectral images obtained with the New Horizons spacecraft that show absorption bands at 1.65 and 2.2 μm. The ammonia signature is spatially coincident with a region of past extensional tectonic activity (Virgil Fossae) where the presence of H2O ice is prominent. Ammonia in liquid water profoundly depresses the freezing point of the mixture. Ammoniated ices are believed to be geologically short lived when irradiated with ultraviolet photons or charged particles. Thus, the presence of NH3 on a planetary surface is indicative of a relatively recent deposition or possibly through exposure by some geological process. In the present case, the areal distribution is more suggestive of cryovolcanic emplacement, however, adding to the evidence for ongoing geological activity on Pluto and the possible presence of liquid water at depth today.

Impact craters on Pluto and Charon indicate a deficit of small Kuiper belt objects

The flyby of Pluto and Charon by the New Horizons spacecraft provided high-resolution images of cratered surfaces embedded in the Kuiper belt, an extensive region of bodies orbiting beyond Neptune. Impact craters on Pluto and Charon were formed by collisions with other Kuiper belt objects (KBOs) with diameters from ~40 kilometers to ~300 meters, smaller than most KBOs observed directly by telescopes. We find a relative paucity of small craters ≲13 kilometers in diameter, which cannot be explained solely by geological resurfacing. This implies a deficit of small KBOs (≲1 to 2 kilometers in diameter). Some surfaces on Pluto and Charon are likely ≳4 billion years old, thus their crater records provide information on the size-frequency distribution of KBOs in the early Solar System.

Dunes on Pluto

The surface of Pluto is more geologically diverse and dynamic than had been expected, but the role of its tenuous atmosphere in shaping the landscape remains unclear. We describe observations from the New Horizons spacecraft of regularly spaced, linear ridges whose morphology, distribution, and orientation are consistent with being transverse dunes. These are located close to mountainous regions and are orthogonal to nearby wind streaks. We demonstrate that the wavelength of the dunes (~0.4 to 1 kilometer) is best explained by the deposition of sand-sized (~200 to ~300 micrometer) particles of methane ice in moderate winds (<10 meters per second). The undisturbed morphology of the dunes, and relationships with the underlying convective glacial ice, imply that the dunes have formed in the very recent geological past.

Vacuum ultraviolet photoabsorption of prime ice analogues of Pluto and Charon

Here we present the first Vacuum UltraViolet (VUV) photoabsorption spectra of ice analogues of Pluto and Charon ice mixtures. For Pluto the ice analogue is an icy mixture containing nitrogen (N₂), carbon monoxide (CO), methane (CH₄) and water (H₂O) prepared with a 100:1:1:3 ratio, respectively. Photoabsorption of icy mixtures with and without H₂O were recorded and no significant changes in the spectra due to presence of H₂O were observed. For Charon a VUV photoabsorption spectra of an ice analogue containing ammonia (NH₃) and H₂O prepared with a 1:1 ratio was recorded, a spectrum of ammonium hydroxide (NH₄OH) was also recorded. These spectra may help to interpret the P-Alice data from New Horizons.

Vital Signs: Seismology of Icy Ocean Worlds

Ice-covered ocean worlds possess diverse energy sources and associated mechanisms that are capable of driving significant seismic activity, but to date no measurements of their seismic activity have been obtained. Such investigations could reveal the transport properties and radial structures, with possibilities for locating and characterizing trapped liquids that may host life and yielding critical constraints on redox fluxes and thus on habitability. Modeling efforts have examined seismic sources from tectonic fracturing and impacts. Here, we describe other possible seismic sources, their associations with science questions constraining habitability, and the feasibility of implementing such investigations. We argue, by analogy with the Moon, that detectable seismic activity should occur frequently on tidally flexed ocean worlds. Their ices fracture more easily than rocks and dissipate more tidal energy than the <1 GW of the Moon and Mars. Icy ocean worlds also should create less thermal noise due to their greater distance and consequently smaller diurnal temperature variations. They also lack substantial atmospheres (except in the case of Titan) that would create additional noise. Thus, seismic experiments could be less complex and less susceptible to noise than prior or planned planetary seismology investigations of the Moon or Mars. Key Words: Seismology-Redox-Ocean worlds-Europa-Ice-Hydrothermal. Astrobiology 18, 37-53.

Global drainage patterns and the origins of topographic relief on Earth, Mars, and Titan

Rivers have eroded the topography of Mars, Titan, and Earth, creating diverse landscapes. However, the dominant processes that generated topography on Titan (and to some extent on early Mars) are not well known. We analyzed drainage patterns on all three bodies and found that large drainages, which record interactions between deformation and erosional modification, conform much better to long-wavelength topography on Titan and Mars than on Earth. We use a numerical landscape evolution model to demonstrate that short-wavelength deformation causes drainage directions to diverge from long-wavelength topography, as observed on Earth. We attribute the observed differences to ancient long-wavelength topography on Mars, recent or ongoing generation of long-wavelength relief on Titan, and the creation of short-wavelength relief by plate tectonics on Earth.

Flow and evolution of ice-sucrose crystal mushes

We study the rheology of suspensions of ice crystals at moderate to high volume fractions in a sucrose solution in which they are partially soluble, a model system for a wide class of crystal mushes or slurries. Under step changes in shear rate, the viscosity changes to a relaxed value over several minutes, in a manner well fitted by a single exponential. The behavior of the relaxed viscosity is power-law shear thinning with shear rate, with an exponent of -1.76±0.25, so that shear stress falls with increasing shear rate. On longer time scales, the crystals ripen (leading to a falling viscosity) so that the mean radius increases with time to the power 0.14±0.07. We speculate that this unusually small exponent is due to the interaction of classical ripening dynamics with abrasion or breakup under flow. We compare the rheological behavior to mechanistic models based on flow-induced aggregation and breakup of crystal clusters, finding that the exponents can be predicted from liquid phase sintering and breakup by brittle fracture.

Charon tectonics

New Horizons images of Pluto's companion Charon show a variety of terrains that display extensional tectonic features, with relief surprising for this relatively small world. These features suggest a global extensional areal strain of order 1% early in Charon's history. Such extension is consistent with the presence of an ancient global ocean, now frozen.

Penitentes as the origin of the bladed terrain of Tartarus Dorsa on Pluto

Penitentes are snow and ice features formed by erosion that, on Earth, are characterized by bowl-shaped depressions several tens of centimetres across, whose edges grade into spires up to several metres tall. Penitentes have been suggested as an explanation for anomalous radar data on Europa, but until now no penitentes have been identified conclusively on planetary bodies other than Earth. Regular ridges with spacings of 3,000 to 5,000 metres and depths of about 500 metres with morphologies that resemble penitentes have been observed by the New Horizons spacecraft in the Tartarus Dorsa region of Pluto (220°-250° E, 0°-20° N). Here we report simulations, based upon a recent model representing conditions on Pluto, in which deepening penitentes reproduce both the tri-modal (north-south, east-west and northeast-southwest) orientation and the spacing of the ridges of this bladed terrain. At present, these penitentes deepen by approximately one centimetre per orbital cycle and grow only during periods of relatively high atmospheric pressure, suggesting a formation timescale of several tens of millions of years, consistent with crater ages. This timescale implies that the penitentes formed from initial topographic variations of no more than a few tens of metres, consistent with Pluto's youngest terrains.

Reorientation and faulting of Pluto due to volatile loading within Sputnik Planitia

Pluto is an astoundingly diverse, geologically dynamic world. The dominant feature is Sputnik Planitia-a tear-drop-shaped topographic depression approximately 1,000 kilometres in diameter possibly representing an ancient impact basin. The interior of Sputnik Planitia is characterized by a smooth, craterless plain three to four kilometres beneath the surrounding rugged uplands, and represents the surface of a massive unit of actively convecting volatile ices (N₂, CH₄ and CO) several kilometres thick. This large feature is very near the Pluto-Charon tidal axis. Here we report that the location of Sputnik Planitia is the natural consequence of the sequestration of volatile ices within the basin and the resulting reorientation (true polar wander) of Pluto. Loading of volatile ices within a basin the size of Sputnik Planitia can substantially alter Pluto's inertia tensor, resulting in a reorientation of the dwarf planet of around 60 degrees with respect to the rotational and tidal axes. The combination of this reorientation, loading and global expansion due to the freezing of a possible subsurface ocean generates stresses within the planet's lithosphere, resulting in a global network of extensional faults that closely replicate the observed fault networks on Pluto. Sputnik Planitia probably formed northwest of its present location, and was loaded with volatiles over million-year timescales as a result of volatile transport cycles on Pluto. Pluto's past, present and future orientation is controlled by feedbacks between volatile sublimation and condensation, changing insolation conditions and Pluto's interior structure.

Observed glacier and volatile distribution on Pluto from atmosphere-topography processes

Pluto has a variety of surface frosts and landforms as well as a complex atmosphere. There is ongoing geological activity related to the massive Sputnik Planitia glacier, mostly made of nitrogen (N₂) ice mixed with solid carbon monoxide and methane, covering the 4-kilometre-deep, 1,000-kilometre-wide basin of Sputnik Planitia near the anti-Charon point. The glacier has been suggested to arise from a source region connected to the deep interior, or from a sink collecting the volatiles released planetwide. Thin deposits of N₂ frost, however, were also detected at mid-northern latitudes and methane ice was observed to cover most of Pluto except for the darker, frost-free equatorial regions. Here we report numerical simulations of the evolution of N₂, methane and carbon monoxide on Pluto over thousands of years. The model predicts N₂ ice accumulation in the deepest low-latitude basin and the threefold increase in atmospheric pressure that has been observed to occur since 1988. This points to atmospheric-topographic processes as the origin of Sputnik Planitia's N₂ glacier. The same simulations also reproduce the observed quantities of volatiles in the atmosphere and show frosts of methane, and sometimes N₂, that seasonally cover the mid- and high latitudes, explaining the bright northern polar cap reported in the 1990s and the observed ice distribution in 2015. The model also predicts that most of these seasonal frosts should disappear in the next decade.

The rapid formation of Sputnik Planitia early in Pluto's history

Pluto's Sputnik Planitia is a bright, roughly circular feature that resembles a polar ice cap. It is approximately 1,000 kilometres across and is centred on a latitude of 25 degrees north and a longitude of 175 degrees, almost directly opposite the side of Pluto that always faces Charon as a result of tidal locking. One explanation for its location includes the formation of a basin in a giant impact, with subsequent upwelling of a dense interior ocean. Once the basin was established, ice would naturally have accumulated there. Then, provided that the basin was a positive gravity anomaly (with or without the ocean), true polar wander could have moved the feature towards the Pluto-Charon tidal axis, on the far side of Pluto from Charon. Here we report modelling that shows that ice quickly accumulates on Pluto near latitudes of 30 degrees north and south, even in the absence of a basin, because, averaged over its orbital period, those are Pluto's coldest regions. Within a million years of Charon's formation, ice deposits on Pluto concentrate into a single cap centred near a latitude of 30 degrees, owing to the runaway albedo effect. This accumulation of ice causes a positive gravity signature that locks, as Pluto's rotation slows, to a longitude directly opposite Charon. Once locked, Charon raises a permanent tidal bulge on Pluto, which greatly enhances the gravity signature of the ice cap. Meanwhile, the weight of the ice in Sputnik Planitia causes the crust under it to slump, creating its own basin (as has happened on Earth in Greenland). Even if the feature is now a modest negative gravity anomaly, it remains locked in place because of the permanent tidal bulge raised by Charon. Any movement of the feature away from 30 degrees latitude is countered by the preferential recondensation of ices near the coldest extremities of the cap. Therefore, our modelling suggests that Sputnik Planitia formed shortly after Charon did and has been stable, albeit gradually losing volume, over the age of the Solar System.

The formation of Charon's red poles from seasonally cold-trapped volatiles

A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to tholin-like organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons in the production of this material. The escape of Pluto's atmosphere provides a potential feedstock for a complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon's winter pole was proposed as an explanation for the dark coloration on the basis of an image of Charon's northern hemisphere, but not modelled quantitatively. Here we report images of the southern hemisphere illuminated by Pluto-shine and also images taken during the approach phase that show the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon and the supply and temporary cold-trapping of material escaping from Pluto, as well as the photolytic processing of this material into more complex and less volatile molecules while cold-trapped. The model results are consistent with the proposed mechanism for producing the observed colour pattern on Charon.

Convection in a volatile nitrogen-ice-rich layer drives Pluto's geological vigour

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's vigorous geological activity. Composed of molecular nitrogen, methane, and carbon monoxide ices, but dominated by nitrogen ice, this layer is organized into cells or polygons, typically about 10 to 40 kilometres across, that resemble the surface manifestation of solid-state convection. Here we report, on the basis of available rheological measurements, that solid layers of nitrogen ice with a thickness in excess of about one kilometre should undergo convection for estimated present-day heat-flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-kilometre-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of nitrogen-ice viscosity implies that the ice layer convects in the so-called sluggish lid regime, a unique convective mode not previously definitively observed in the Solar System. Average surface horizontal velocities of a few centimetres a year imply surface transport or renewal times of about 500,000 years, well under the ten-million-year upper-limit crater retention age for Sputnik Planum. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help to explain the high albedos shown by some of these bodies.

The small satellites of Pluto as observed by New Horizons

The New Horizons mission has provided resolved measurements of Pluto's moons Styx, Nix, Kerberos, and Hydra. All four are small, with equivalent spherical diameters of ~40 kilometers for Nix and Hydra and ~10 kilometers for Styx and Kerberos. They are also highly elongated, with maximum to minimum axis ratios of ~2. All four moons have high albedos (~50 to 90%) suggestive of a water-ice surface composition. Crater densities on Nix and Hydra imply surface ages of at least 4 billion years. The small moons rotate much faster than synchronous, with rotational poles clustered nearly orthogonal to the common pole directions of Pluto and Charon. These results reinforce the hypothesis that the small moons formed in the aftermath of a collision that produced the Pluto-Charon binary.