Planetary Interiors

Early models of the composition of the Earth relied heavily on meteorites. In all these models Earth had different layers, each layer corresponded to a different type of meteorite or meteorite component. Later, more realistic models based on analyses of samples from Earth began with Ringwood’s pyrolite composition in the 1960s. Further improvement came with the analyses of rare MgO rich peridotites from a variety of occurrences all over the Earth, as xenoliths enclosed in melts from the upper mantle or as ultramafic massifs, tectonically emplaced on the Earth’s surface. Chemical systematics of these rocks allow the determination of the major element composition of the primitive upper mantle (PUM), the upper mantle after core formation and before extraction of basalts ultimately leading to the formation of the crust. Trace element analyses of upper mantle rocks confirmed their primitive nature. Geochemical and geophysical evidence argue for a bulk Earth mantle of uniform composition, identical to the PUM, also designated as “bulk silicate Earth” (BSE). The formation of a metal core was accompanied by the removal of siderophile and chalcophile elements into the core. Detailed modeling suggests that core formation was an ongoing process parallel to the accretion of Earth. The composition of the core is model dependent and thus uncertain and makes reliable estimates for siderophile and chalcophile element concentrations of bulk Earth difficult. Improved stable isotope analyses show isotopic similarities with noncarbonaceous chondrites (NCC), while the chemical composition of the mantle of the Earth indicates similarities with carbonaceous chondrites (CC). In detail, however, it can be shown that no single known meteorite group, nor any mixture of meteorite groups can match the chemical and isotopic composition of Earth. This conclusion is extremely important for any formation model of the Earth.

Element partitioning—at its most basic—is the distribution of an element of interest between two constituent phases as a function of some process. Major constituent elements generally affect the thermodynamic environment (chemical equilibrium) and therefore trace element partitioning is often considered, as trace elements are present in minute quantities and their equilibrium exchange reactions do not impart significant changes to the larger system. Trace elements are responsive to thermodynamic conditions, and thus they act as passive tracers of chemical reactions without appreciably influencing the bulk reactions themselves. In planetary sciences, the phase pairs typically considered are mineral-melt, metal-silicate, and sulfide-silicate, owing largely to the ubiquity of their coexistence in planetary materials across scales and context, from the micrometer-sized components of meteorites up to the size of planets (thousands of kilometers). It is common to speak of trace elements in terms of their tendency toward forming metallic, sulfidic, or oxide phases, and the terms “siderophile,” “chalcophile,” and “lithophile” (respectively) are used to define these tendencies under what is known as the Goldschmidt Classification scheme. The metric of an element’s tendency to concentrate into one phase relative to another is expressed as the ratio of its concentration (as a weight or molar fraction) in one phase over another, where convention dictates the reference frame as solid over liquid, and metal or sulfide over silicate; this mathematical term is the element’s partition coefficient, or distribution coefficient, between the two respective phases, D M Phase B Phase A (where M is the element of interest, most often reported as molar fraction), or simply D M . In general, trace elements obey Henry’s Law, where the element’s activity and concentration are linearly proportional. Practically speaking, this means that the element is sufficiently dilute in the system such that its atoms interact negligibly with one another compared to their interactions with major element phases, and thus the trace element’s partition coefficient in most settings is not appreciably affected by its concentration. The radius and charge of an element’s ionized species (its ionic radius and valence state)—in relation to either the major element ion for which it is substituting or the lattice site vacancy or interstitial space it is filling—generally determine the likelihood of trace element substitution or vacancy/interstitial fill (along with the net charge of the lattice space). The key energy consideration that underlies an element’s partitioning is the Gibbs free energy of reaction between the phases involved. Gibbs free energy is the change in internal energy associated with a chemical reaction (at a given temperature and pressure) that can be used to do work, and is denoted as Δ G rxn . Reactions with negative Δ G rxn values are spontaneous, and the magnitude of this negative value for a given phase, for example, a metal oxide, denotes the relative affinity of the metal toward forming oxides. That is to say, an element with a highly negative Δ G rxn for its oxide species at relevant pressure-temperature conditions will tend to be found in oxide and silicate minerals, that is, it will be lithophile (and vice versa for siderophile elements). Trace element partitioning systematics in mineral-melt and metal-/sulfide-silicate systems have boundless applications in planetary science. A growing collective understanding of the partition coefficients of elements has been built on decades of physical chemistry, deterministic theory, petrology, experimental petrology, and natural observations. Leveraging this immense intellectual, technical, and methodological foundation, modern trace element partitioning research is particularly aimed at constraining the evolution of plate tectonics on Earth (conditions and timing of onset), understanding the formation history of planetary materials such as chondrite meteorites and their constituents (e.g., chondrules), and de-convolving the multiply operating processes at play during the accretion and differentiation of Earth and other terrestrial planets.

The Earth’s core formed by multiple collisions with differentiated protoplanets. The Hf-W (hafnium-tungsten) isotopic system reveals that these collisions took place over a timescale of tens of megayears (Myr), in agreement with accretion simulations. The degree to which the iron and silicates re-equilibrated during each collision is uncertain and affects the apparent core age derived from tungsten isotopic measurements. Seismological data reveal that the core contains light elements in addition to Fe-Ni, and the outer core is more enriched in such elements than the inner core. Because O is excluded efficiently from solid iron, O is almost certainly an important constituent of the outer core. The identity of other elements is less certain, despite intensive measurements of their effects on seismic velocities, densities, and partitioning behavior at appropriate pressures and temperatures. Si and O are very likely present, with perhaps some S; C and H are less likely. Si and Mg may have exsolved over time, potentially helping to drive the geodynamo and producing a low-density layer at the top of the core. Radioactive elements (U, Th, K) are unlikely to be present in important concentrations. The cooling of the core is controlled by the mantle’s ability to extract heat. The geodynamo has existed for at least 3.5 gigayears (Gyr), placing a lower bound on the heat flow out of the core. Because the thermal conductivity of the core is uncertain by a factor of ~3, the lower bound on this heat flow is similarly uncertain. Once the inner core started to crystallize, additional sources of energy were available to power the geodynamo. Inner core crystallization likely started in the time range 0.5 to 2.0 Gyr Before Present (BP); paleomagnetic arguments have been advanced for inner core growth starting at several different epochs within this time range.

This article consists of three sections. The first discusses how we determine satellite internal structures and what we know about them. The primary probes of internal structure are measurements of magnetic induction, gravity, and topography, as well as rotation state and orientation. Enceladus, Europa, Ganymede, Callisto, Titan, and (perhaps) Pluto all have subsurface oceans; Callisto and Titan may be only incompletely differentiated. The second section describes dynamical processes that affect satellite interiors and surfaces: tidal and radioactive heating, flexure and relaxation, convection, cryovolcanism, true polar wander, non-synchronous rotation, orbital evolution, and impacts. The final section discusses how the satellites formed and evolved. Ancient tidal heating episodes and subsequent refreezing of a subsurface ocean are the likeliest explanation for the deformation observed at Ganymede, Tethys, Dione, Rhea, Miranda, Ariel, and Titania. The high heat output of Enceladus is a consequence of Saturn’s highly dissipative interior, but the dissipation rate is strongly frequency-dependent and does not necessarily imply that Saturn’s moons are young. Major remaining questions include the origins of Titan’s atmosphere and high eccentricity, the regular density progression in the Galilean satellites, and the orbital evolution of the Saturnian and Uranian moons.

Probing the interiors of the gaseous giant planets in our solar system is not an easy task. It requires a set of accurate measurements combined with theoretical models that are used to infer the planetary composition and its depth dependence. The masses of Jupiter and Saturn are 317.83 and 95.16 Earth masses (M ⊕ ), respectively, and since a few decades, it has been known that they mostly consist of hydrogen and helium. The mass of heavy elements (all elements heavier than helium) is not well determined, nor are their distribution within the planets. While the heavy elements are not the dominating materials inside Jupiter and Saturn, they are the key to understanding the planets’ formation and evolutionary histories. The planetary internal structure is inferred from theoretical models that fit the available observational constraints by using theoretical equations of states (EOSs) for hydrogen, helium, their mixtures, and heavier elements (typically rocks and/or ices). However, there is no unique solution for determining the planetary structure and the results depend on the used EOSs as well as the model assumptions imposed by the modeler. Major model assumptions that can affect the derived internal structure include the number of layers, the heat transport mechanism within the planet (and its entropy), the nature of the core (compact vs. diluted), and the location (pressure) of separation between the two envelopes. Alternative structure models assume a less distinct division between the layers and /or a non-homogenous distribution of the heavy elements. The fact that the behavior of hydrogen at high pressures and temperatures is not perfectly known and that helium may separate from hydrogen at the deep interior add sources of uncertainty to structure models. In the 21st century, with accurate measurements of the gravitational fields of Jupiter and Saturn from the Juno and Cassini missions, structure models can be further constrained. At the same time, these measurements introduce new challenges for planetary modelers.

Many large volcanic channel systems are recognized at the surfaces of rocky bodies of the inner solar system. The more than 200 channels known for the Moon mainly have simple sinuous forms with widths of up to several kilometers and lengths of up to hundreds of kilometers, typically commencing at topographic depressions and extending downslope until they fade into associated volcanic units. The Rima Hadley system was a key target of the Apollo 15 mission and was confirmed as a product of volcanic processes related to the emplacement of lavas in the Palus Putredinis region of Mare Imbrium. The more than 200 channels known for Venus are in many cases morphologically similar to sinuous lunar rilles, but some systems are especially large and complex, with widths of up to tens of kilometers and lengths that can exceed 1,000 km. Such systems typically commence at structural features or in regions of disturbed terrain and possess anastomosing reaches associated with prominent streamlined uplands. In contrast, Venusian canali typically maintain sinuous forms with widths of only a few kilometers but remarkably can be characterized by lengths of thousands of kilometers. Some Venusian channels were involved in the emplacement of fluidized ejecta in the vicinities of impact craters whereas others may have formed in such environments as a result of later volcanic events. The 10 large volcanic channels that are recognized on Mercury have lengths no greater than ~161 km but can have widths of up to several tens of kilometers. These systems developed as conduits for voluminous lavas that extend across adjacent impact basins. Terrestrial komatiitic channels of Archean and Proterozoic ages can have sizes that are comparable to those of lunar rilles, and the formation of these systems is likely to have played an important role in the development of associated Ni-Cu-(PGE) ores. The outflow channels of Mars have widths of up to tens of kilometers and lengths of up to thousands of kilometers and are widely interpreted as aqueous systems formed by catastrophic discharges from aquifers, but the properties and geological associations of these features and numerous other large Martian channel systems are arguably well aligned with those expected of volcanic origins. Overall, large volcanic channels of the inner solar system are mainly ancient products of the emplacement of low-viscosity lava flows of mafic or ultramafic composition, involving eruptions that were characterized by extraordinarily high effusion rates and total lava volumes that in some cases are likely to have been as great as those that characterized some Large Igneous Provinces on Earth. The deeply rooted igneous plumbing systems most favorable to the development of large volcanic channels would have been especially common in the earlier history of the solar system, when the interior temperatures of rocky bodies were greater than today. The early development of large volcanic channel systems is likely typical of the geological histories of large rocky bodies in the universe.

Since 1973 space missions carrying vector magnetometers have shown that most, but not all, solar system planets have a global magnetic field of internal origin. They have also revealed a surprising diversity in terms of field strength and morphology. While Jupiter’s field, like that of Earth, is dominated by a dipole moderately tilted relative to the planet’s spin axis, the fields of Uranus and Neptune are multipole-dominated, whereas those of Saturn and Mercury are highly symmetric relative to the rotation axis. Planetary magnetism originates from a dynamo process, which requires a fluid and electrically conducting region in the interior with sufficiently rapid and complex flow. The magnetic fields are of interest for three reasons: (i) they provide ground truth for dynamo theory, (ii) the magnetic field controls how the planet interacts with its space environment, for example, the solar wind, and (iii) the existence or nonexistence and the properties of the field enable us to draw inferences on the constitution, dynamics, and thermal evolution of the planet’s interior. Numerical simulations of the geodynamo, in which convective flow in a rapidly rotating spherical shell representing the outer liquid iron core of the Earth leads to induction of electric currents, have successfully reproduced many observed properties of the geomagnetic field. They have also provided guidelines on the factors controlling magnetic field strength and morphology. For numerical reasons the simulations must employ viscosities far greater than those inside planets and it is debatable whether they capture the correct physics of planetary dynamo processes. Nonetheless, such models have been adapted to test concepts for explaining magnetic field properties of other planets. For example, they show that a stable stratified conducting layer above the dynamo region is a plausible cause for the strongly axisymmetric magnetic fields of Mercury or Saturn.

Earth’s moon, hereafter referred to as “the Moon,” has been an object of intense study since before the time of the Apollo and Luna missions to the lunar surface and associated sample returns. As a differentiated rocky body and as Earth’s companion in the solar system, much study has been given to aspects such as the Moon’s surface characteristics, composition, interior, geologic history, origin, and what it records about the early history of the Earth-Moon system and the evolution of differentiated rocky bodies in the solar system. Much of the Apollo and post-Apollo knowledge came from surface geologic exploration, remote sensing, and extensive studies of the lunar samples. After a hiatus of nearly two decades following the end of Apollo and Luna missions, a new era of lunar exploration began with a series of orbital missions, including missions designed to prepare the way for longer duration human use and further exploration of the Moon. Participation in these missions has become international. The more recent missions have provided global context and have investigated composition, mineralogy, topography, gravity, tectonics, thermal evolution of the interior, thermal and radiation environments at the surface, exosphere composition and phenomena, and characteristics of the poles with their permanently shaded cold-trap environments. New samples were recognized as a class of achondrite meteorites, shown through geochemical and mineralogical similarities to have originated on the Moon. New sample-based studies with ever-improving analytical techniques and approaches have also led to significant discoveries such as the determination of volatile contents, including intrinsic H contents of lunar minerals and glasses. The Moon preserves a record of the impact history of the solar system, and new developments in timing of events, sample based and model based, are leading to a new reckoning of planetary chronology and the events that occurred in the early solar system. The new data provide the grist to test models of formation of the Moon and its early differentiation, and its thermal and volcanic evolution. Thought to have been born of a giant impact into early Earth, new data are providing key constraints on timing and process. The new data are also being used to test hypotheses and work out details such as for the magma ocean concept, the possible existence of an early magnetic field generated by a core dynamo, the effects of intense asteroidal and cometary bombardment during the first 500 million–600 million years, sequestration of volatile compounds at the poles, volcanism through time, including new information about the youngest volcanism on the Moon, and the formation and degradation processes of impact craters, so well preserved on the Moon. The Moon is a natural laboratory and cornerstone for understanding many processes operating in the space environment of the Earth and Moon, now and in the past, and of the geologic processes that have affected the planets through time. The Moon is a destination for further human exploration and activity, including use of valuable resources in space. It behooves humanity to learn as much about Earth’s nearest neighbor in space as possible.

Mercury, like its inner Solar System planetary neighbors Venus, Mars, and the Moon, shows no evidence of having ever undergone plate tectonics. Nonetheless, the innermost planet boasts a long record of tectonic deformation. The most prominent manifestation of this history is a population of large scarps that occurs throughout the planet’s cratered terrains; some of these scarps rise kilometers above the surrounding landscape. Mercury’s smooth plains, the majority of which are volcanic and occupy over a quarter of the planet, abound with low-relief ridges. The scarps and ridges are underlain by thrust faults and point to a tectonic history dominated by crustal shortening. At least some of the shortening strain recorded by the ridges may reflect subsidence of the lavas in which they formed, but the widespread distribution of scarps attests to a planetwide process of global contraction, wherein Mercury experienced a reduction in volume as its interior cooled through time. The onset of this phenomenon placed the lithosphere into a net state of horizontal compression, and accounts for why Mercury hosts only a few instances of extensional structures. These landforms, shallow troughs that form complex networks, occur almost wholly in volcanically flooded impact craters and basins and developed as those lavas cooled and thermally contracted. Tellingly, widespread volcanism on Mercury ended at around the same time the population of scarps began to form. Explosive volcanism endured beyond this point, but almost exclusively at sites of lithospheric weakness, where large faults penetrate deep into the interior. These observations are consistent with decades-old predictions that global contraction would shut off major volcanic activity, and illustrate how closely Mercury’s tectonic and volcanic histories are intertwined. The tectonic character of Mercury is thus one of sustained crustal shortening with only localized extension, which started almost four billion years ago and extends into the geologically recent past. This character somewhat resembles that of the Moon, but differs substantially from those of Earth, Venus, or Mars. Mercury may represent how small rocky planets tectonically evolve and could provide a basis for understanding the geological properties of similarly small worlds in orbit around other stars.