Isotopic evidence of long-lived volcanism on Io.

Jupiter's moon Io hosts extensive volcanism, driven by tidal heating. The isotopic composition of Io's inventory of volatile chemical elements, including sulfur and chlorine, reflects its outgassing and mass-loss history and thus records information about its evolution. We used submillimeter observations of Io's atmosphere to measure sulfur isotopes in gaseous sulfur dioxide and sulfur monoxide, and chlorine isotopes in gaseous sodium chloride and potassium chloride. We find 34S/32S = 0.0595 ± 0.0038 (equivalent to δ34S = +347 ± 86‰), which is highly enriched compared to average Solar System values and indicates that Io has lost 94 to 99% of its available sulfur. Our measurement of 37Cl/35Cl = 0.403 ± 0.028 (δ37Cl = +263 ± 88‰) shows that chlorine is similarly enriched. These results indicate that Io has been volcanically active for most (or all) of its history, with potentially higher outgassing and mass-loss rates at earlier times.

Investigating Europa's Habitability with the Europa Clipper.

The habitability of Europa is a property within a system, which is driven by a multitude of physical and chemical processes and is defined by many interdependent parameters, so that its full characterization requires collaborative investigation. To explore Europa as an integrated system to yield a complete picture of its habitability, the Europa Clipper mission has three primary science objectives: (1) characterize the ice shell and ocean including their heterogeneity, properties, and the nature of surface-ice-ocean exchange; (2) characterize Europa's composition including any non-ice materials on the surface and in the atmosphere, and any carbon-containing compounds; and (3) characterize Europa's geology including surface features and localities of high science interest. The mission will also address several cross-cutting science topics including the search for any current or recent activity in the form of thermal anomalies and plumes, performing geodetic and radiation measurements, and assessing high-resolution, co-located observations at select sites to provide reconnaissance for a potential future landed mission. Synthesizing the mission's science measurements, as well as incorporating remote observations by Earth-based observatories, the James Webb Space Telescope, and other space-based resources, to constrain Europa's habitability, is a complex task and is guided by the mission's Habitability Assessment Board (HAB).

The internal structure of Eris inferred from its spin and orbit evolution.

The large Kuiper Belt object Eris is tidally locked to its small companion Dysnomia. Recently obtained bounds on the mass of Dysnomia demonstrate that Eris must be unexpectedly dissipative for it to have despun over the age of the solar system. Here, we show that Eris must have differentiated into an ice shell and rocky core to explain the dissipation. We further demonstrate that Eris's ice shell must be convecting to be sufficiently dissipative, which distinguishes it from Pluto's conductive shell. The difference is likely due to Eris's apparent depletion in volatiles compared with Pluto, perhaps as the result of a more energetic impact.

Origin and Evolution of Enceladus's Tidal Dissipation.

Enceladus possesses a subsurface ocean beneath a conductive ice shell. Based on shell thickness models, the estimated total conductive heat loss from Enceladus is 25-40 GW; the measured heat output from the South Polar Terrain (SPT) is 4-19 GW. The present-day SPT heat flux is of order 100 mWm-2, comparable to estimated paleo-heat fluxes for other regions of Enceladus. These regions have nominal ages of about 2 Ga, but the estimates are uncertain because the impactor flux in the Saturnian system may not resemble that elsewhere. Enceladus's measured rate of orbital expansion implies a low dissipation factor Qp for Saturn, with Qp≈3×103 (neglecting the role of Dione). This value implies that Enceladus's present-day equilibrium tidal heat production (roughly 50 GW, but with large uncertainties) is in approximate balance with its heat loss. If Qp is constant, Enceladus cannot be older than 1.5 Gyr (because otherwise it would have migrated more than is permissible). However, Saturn's dissipation may be better described by the "resonance-locking" theory, in which case Enceladus's orbit may have only evolved outwards by about 35% over the age of the Solar System. In the constant-Qp scenario, any ancient tidal heating events would have been too energetic to be consistent with the observations. Because resonance-locking makes capture into earlier mean-motion orbital resonances less likely, the inferred ancient heating episodes probably took place when the current orbital resonance was already established. In the resonance-locking scenario, tidal heating did not change significantly over time, allowing for a long-lived ocean and a relatively stable ice shell. If so, Enceladus is an attractive target for future exploration from a habitability standpoint.

Exploring the Interior of Europa with the Europa Clipper.

The Galileo mission to Jupiter revealed that Europa is an ocean world. The Galileo magnetometer experiment in particular provided strong evidence for a salty subsurface ocean beneath the ice shell, likely in contact with the rocky core. Within the ice shell and ocean, a number of tectonic and geodynamic processes may operate today or have operated at some point in the past, including solid ice convection, diapirism, subsumption, and interstitial lake formation. The science objectives of the Europa Clipper mission include the characterization of Europa's interior; confirmation of the presence of a subsurface ocean; identification of constraints on the depth to this ocean, and on its salinity and thickness; and determination of processes of material exchange between the surface, ice shell, and ocean. Three broad categories of investigation are planned to interrogate different aspects of the subsurface structure and properties of the ice shell and ocean: magnetic induction, subsurface radar sounding, and tidal deformation. These investigations are supplemented by several auxiliary measurements. Alone, each of these investigations will reveal unique information. Together, the synergy between these investigations will expose the secrets of the Europan interior in unprecedented detail, an essential step in evaluating the habitability of this ocean world.

Hadaean to Palaeoarchaean stagnant-lid tectonics revealed by zircon magnetism.

Plate tectonics is a fundamental factor in the sustained habitability of Earth, but its time of onset is unknown, with ages ranging from the Hadaean to Proterozoic eons1-3. Plate motion is a key diagnostic to distinguish between plate and stagnant-lid tectonics, but palaeomagnetic tests have been thwarted because the planet's oldest extant rocks have been metamorphosed and/or deformed4. Herein, we report palaeointensity data from Hadaean-age to Mesoarchaean-age single detrital zircons bearing primary magnetite inclusions from the Barberton Greenstone Belt of South Africa5. These reveal a pattern of palaeointensities from the Eoarchaean (about 3.9 billion years ago (Ga)) to Mesoarchaean (about 3.3 Ga) eras that is nearly identical to that defined by primary magnetizations from the Jack Hills (JH; Western Australia)6,7, further demonstrating the recording fidelity of select detrital zircons. Moreover, palaeofield values are nearly constant between about 3.9 Ga and about 3.4 Ga. This indicates unvarying latitudes, an observation distinct from plate tectonics of the past 600 million years (Myr) but predicted by stagnant-lid convection. If life originated by the Eoarchaean8, and persisted to the occurrence of stromatolites half a billion years later9, it did so when Earth was in a stagnant-lid regime, without plate-tectonics-driven geochemical cycling.

Viscous relaxation as a probe of heat flux and crustal plateau composition on Venus.

It has recently been suggested that deformed crustal plateaus on Venus may be composed of felsic (silica-rich) rocks, possibly supporting the idea of an ancient ocean there. However, these plateaus have a tendency to collapse owing to flow of the viscous lower crust. Felsic minerals, especially water-bearing ones, are much weaker and thus lead to more rapid collapse, than more mafic minerals. We model plateau topographic evolution using a non-Newtonian viscous relaxation code. Despite uncertainties in the likely crustal thickness and surface heat flux, we find that quartz-dominated rheologies relax too rapidly to be plausible plateau-forming material. For plateaus dominated by a dry anorthite rheology, survival is possible only if the background crustal thickness is less than 29 km, unless the heat flux on Venus is less than the radiogenic lower bound of 34 [Formula: see text]. Future spacecraft determinations of plateau crustal thickness and mineralogy will place firmer constraints on Venus's heat flux.

Loss of a satellite could explain Saturn's obliquity and young rings.

The origin of Saturn's ~26.7° obliquity and ~100-million-year-old rings is unknown. The observed rapid outward migration of Saturn's largest satellite, Titan, could have raised Saturn's obliquity through a spin-orbit precession resonance with Neptune. We use Cassini data to refine estimates of Saturn's moment of inertia, finding that it is just outside the range required for the resonance. We propose that Saturn previously had an additional satellite, which we name Chrysalis, that caused Saturn's obliquity to increase through the Neptune resonance. Destabilization of Chrysalis's orbit ~100 million years ago can then explain the proximity of the system to the resonance and the formation of the rings through a grazing encounter with Saturn.

The Tides of Enceladus' Porous Core.

The inferred density of Enceladus' core, together with evidence of hydrothermal activity within the moon, suggests that the core is porous. Tidal dissipation in an unconsolidated core has been proposed as the main source of Enceladus' geological activity. However, the tidal response of its core has generally been modeled assuming it behaves viscoelastically rather than poroviscoelastically. In this work, we analyze the poroviscoelastic response to better constrain the distribution of tidal dissipation within Enceladus. A poroviscoelastic body has a different tidal response than a viscoelastic one; pressure within the pores alters the stress field and induces a Darcian porous flow. This flow represents an additional pathway for energy dissipation. Using Biot's theory of poroviscoelasticity, we develop a new framework to obtain the tidal response of a spherically symmetric, self-gravitating moon with porous layers and apply it to Enceladus. We show that the boundary conditions at the interface of the core and overlying ocean play a key role in the tidal response. The ocean hinders the development of a large-amplitude Darcian flow, making negligible the Darcian contribution to the dissipation budget. We therefore infer that Enceladus' core can be the source of its geological activity only if it has a low rigidity and a very low viscosity. A future mission to Enceladus could test this hypothesis by measuring the phase lags of tidally induced changes of gravitational potential and surface displacements.

Early Cambrian renewal of the geodynamo and the origin of inner core structure.

Paleomagnetism can elucidate the origin of inner core structure by establishing when crystallization started. The salient signal is an ultralow field strength, associated with waning thermal energy to power the geodynamo from core-mantle heat flux, followed by a sharp intensity increase as new thermal and compositional sources of buoyancy become available once inner core nucleation (ICN) commences. Ultralow fields have been reported from Ediacaran (~565 Ma) rocks, but the transition to stronger strengths has been unclear. Herein, we present single crystal paleointensity results from early Cambrian (~532 Ma) anorthosites of Oklahoma. These yield a time-averaged dipole moment 5 times greater than that of the Ediacaran Period. This rapid renewal of the field, together with data defining ultralow strengths, constrains ICN to ~550 Ma. Thermal modeling using this onset age suggests the inner core had grown to 50% of its current radius, where seismic anisotropy changes, by ~450 Ma. We propose the seismic anisotropy of the outermost inner core reflects development of a global spherical harmonic degree-2 deep mantle structure at this time that has persisted to the present day. The imprint of an older degree-1 pattern is preserved in the innermost inner core.

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 (>104 km3) 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.

Timing of Martian Core Formation from Models of Hf-W Evolution Coupled with N-body Simulations.

Determining how and when Mars formed has been a long-standing challenge for planetary scientists. The size and orbit of Mars are difficult to reproduce in classical simulations of planetary accretion, and this has inspired models of inner solar system evolution that are tuned to produce Mars-like planets. However, such models are not always coupled to geochemical constraints. Analyses of Martian meteorites using the extinct hafnium-tungsten (Hf-W) radioisotopic system, which is sensitive to the timing of core formation, have indicated that the Martian core formed within a few million years of the start of the solar system itself. This has been interpreted to suggest that, unlike Earth's protracted accretion, Mars grew to its modern size very rapidly. These arguments, however, generally rely on simplified growth histories for Mars. Here, we combine likely accretionary histories from a large number of N-body simulations with calculations of metal-silicate partitioning and Hf-W isotopic evolution during core formation to constrain the range of conditions that could have produced Mars. We find that there is no strong correlation between the final masses or orbits of simulated Martian analogs and their 182W anomalies, and that it is readily possible to produce Mars-like Hf-W isotopic compositions for a variety of accretionary conditions. The Hf-W signature of Mars is very sensitive to the oxygen fugacity (fO2) of accreted material because the metal-silicate partitioning behavior of W is strongly dependent on redox conditions. The average fO2 of Martian building blocks must fall in the range of 1.3-1.6 log units below the iron-wüstite buffer to produce a Martian mantle with the observed Hf/W ratio. Other geochemical properties (such as sulfur content) also influence Martian 182W signatures, but the timing of accretion is a more important control. We find that while Mars must have accreted most of its mass within ~5 million years of solar system formation to reproduce the Hf-W isotopic constraints, it may have continued growing afterwards for over 50 million years. There is a high probability of simultaneously matching the orbit, mass, and Hf-W signature of Mars even in cases of prolonged accretion if giant impactor cores were poorly equilibrated and merged directly with the proto-Martian core.

The origin of the Moon's Earth-like tungsten isotopic composition from dynamical and geochemical modeling.

The Earth and Moon have identical or very similar isotopic compositions for many elements, including tungsten. However, canonical models of the Moon-forming impact predict that the Moon should be made mostly of material from the impactor, Theia. Here we evaluate the probability of the Moon inheriting its Earth-like tungsten isotopes from Theia in the canonical giant impact scenario, using 242 N-body models of planetary accretion and tracking tungsten isotopic evolution, and find that this probability is <1.6-4.7%. Mixing in up to 30% terrestrial materials increases this probability, but it remains <10%. Achieving similarity in stable isotopes is also a low-probability outcome, and is controlled by different mechanisms than tungsten. The Moon's stable isotopes and tungsten isotopic composition are anticorrelated due to redox effects, lowering the joint probability to significantly less than 0.08-0.4%. We therefore conclude that alternate explanations for the Moon's isotopic composition are likely more plausible.

Astronomical context of Solar System formation from molybdenum isotopes in meteorite inclusions.

Calcium-aluminum-rich inclusions (CAIs) in meteorites are the first solids to have formed in the Solar System, defining the epoch of its birth on an absolute time scale. This provides a link between astronomical observations of star formation and cosmochemical studies of Solar System formation. We show that the distinct molybdenum isotopic compositions of CAIs cover almost the entire compositional range of material that formed in the protoplanetary disk. We propose that CAIs formed while the Sun was in transition from the protostellar to pre-main sequence (T Tauri) phase of star formation, placing Solar System formation within an astronomical context. Our results imply that the bulk of the material that formed the Sun and Solar System accreted within the CAI-forming epoch, which lasted less than 200,000 years.

Paleomagnetism indicates that primary magnetite in zircon records a strong Hadean geodynamo.

Determining the age of the geomagnetic field is of paramount importance for understanding the evolution of the planet because the field shields the atmosphere from erosion by the solar wind. The absence or presence of the geomagnetic field also provides a unique gauge of early core conditions. Evidence for a geomagnetic field 4.2 billion-year (Gy) old, just a few hundred million years after the lunar-forming giant impact, has come from paleomagnetic analyses of zircons of the Jack Hills (Western Australia). Herein, we provide new paleomagnetic and electron microscope analyses that attest to the presence of a primary magnetic remanence carried by magnetite in these zircons and new geochemical data indicating that select Hadean zircons have escaped magnetic resetting since their formation. New paleointensity and Pb-Pb radiometric age data from additional zircons meeting robust selection criteria provide further evidence for the fidelity of the magnetic record and suggest a period of high geomagnetic field strength at 4.1 to 4.0 billion years ago (Ga) that may represent efficient convection related to chemical precipitation in Earth's Hadean liquid iron core.

Explaining the Galilean Satellites' Density Gradient by Hydrodynamic Escape.

The Galilean satellites exhibit a monotonic decrease in density (and increase in ice mass fraction) with distance from Jupiter (Pollack & Fanale 1982). Whether this is because of the background conditions when they formed (Lunine & Stevenson 1982; Canup & Ward 2002; Mosqueira & Estrada 2003a; Ronnet et al. 2017), the process of accretion itself (Dwyer et al. 2013), or later loss due to tidal heating (Canup & Ward 2009), has been in dispute for forty years. We find that a hitherto largely neglected process - vapor loss driven by accretional heating (Kuramoto & Matsui 1994) - can reproduce the observed density trend for accretion timescales ≳300 kyr, consistent with gas-starved satellite formation models (Canup & Ward 2002, 2006). In this model both Io and Europa develop an early surface liquid water ocean. Vapor escape from this ocean causes the water inventories of Io and Europa to be completely and mostly lost, respectively. Isotopic fractionation arising from vapor loss means that Europa will develop a higher D/H ratio compared with Ganymede and Callisto. We make predictions that can be tested with in situ measurements of D/H of potential Europa plumes (Roth et al. 2014) by the Europa Clipper spacecraft, or infrared spectroscopic determinations (Clark et al. 2019) of D/H at all three bodies.

The top-down solidification of iron asteroids driving dynamo evolution.

The cores of some small planetesimals, such as asteroid (16) Psyche, are thought to have been exposed through collisions during the early solar system that removed their mantles. These small bodies likely solidified from the top down representing a fundamentally different solidification regime to that of Earth's core. Here we derive simplified models of the downwards solidification of the metallic crust, and consider thermal convection and the potential for viscous delamination of the weak, warm base of the crust to provide a buoyancy flux sufficient to drive a dynamo. Thermal buoyancy is very short lived (~1000 years), and therefore cannot be the source of measured paleomagnetic remanence. In contrast, viscous delamination is found to provide a long-lasting buoyancy flux sufficient to generate an intense, multipolar magnetic field, while not greatly affecting the crustal solidification time. Our results suggest that a Psyche-sized (150 km radius) body solidified in roughly 6.7 - 20 Myr, and that delamination produced a strong magnetic field over much of this time. Finally, including light, insoluble impurities, such as sulfur, results in a partially solid mushy zone at the base of the crust. This further weakens the base of the crust and results in smaller scale delamination events. Despite a significant change in the dynamics of delamination, the time to total solidification and the predicted properties of the magnetic field are broadly comparable to the sulfur-free case, though we argue this may result in observable compositional stratification of the body.

Constraints on terrestrial planet formation timescales and equilibration processes in the Grand Tack scenario from Hf-W isotopic evolution.

We examine 141 N-body simulations of terrestrial planet late-stage accretion that use the Grand Tack scenario, coupling the collisional results with a hafnium-tungsten (Hf-W) isotopic evolution model. Accretion in the Grand Tack scenario results in faster planet formation than classical accretion models because of higher planetesimal surface density induced by a migrating Jupiter. Planetary embryos which grow rapidly experience radiogenic ingrowth of mantle tungsten which is inconsistent with the measured terrestrial value, unless much of the tungsten is removed by an impactor core that mixes thoroughly with the target mantle. For physically Earth-like surviving planets, we find that the fraction of equilibrating impactor core kcore ≥ 0.6 is required to produce results agreeing with observed terrestrial tungsten anomalies (assuming equilibration with relatively large volumes of target mantle material; smaller equilibrating mantle volumes would require even larger kcore ). This requirement of substantial core re-equilibration may be difficult to reconcile with fluid dynamical predictions and hydrocode simulations of mixing during large impacts, and hence this result disfavors the rapid planet building of Grand Tack accretion.

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.