Geophysical evidence for an enriched molten silicate layer above Mars's core.

The detection of deep reflected S waves on Mars inferred a core size of 1,830 ± 40 km (ref. 1), requiring light-element contents that are incompatible with experimental petrological constraints. This estimate assumes a compositionally homogeneous Martian mantle, at odds with recent measurements of anomalously slow propagating P waves diffracted along the core-mantle boundary2. An alternative hypothesis is that Mars's mantle is heterogeneous as a consequence of an early magma ocean that solidified to form a basal layer enriched in iron and heat-producing elements. Such enrichment results in the formation of a molten silicate layer above the core, overlain by a partially molten layer3. Here we show that this structure is compatible with all geophysical data, notably (1) deep reflected and diffracted mantle seismic phases, (2) weak shear attenuation at seismic frequency and (3) Mars's dissipative nature at Phobos tides. The core size in this scenario is 1,650 ± 20 km, implying a density of 6.5 g cm-3, 5-8% larger than previous seismic estimates, and can be explained by fewer, and less abundant, alloying light elements than previously required, in amounts compatible with experimental and cosmochemical constraints. Finally, the layered mantle structure requires external sources to generate the magnetic signatures recorded in Mars's crust.

First observations of core-transiting seismic phases on Mars.

We present the first observations of seismic waves propagating through the core of Mars. These observations, made using seismic data collected by the InSight geophysical mission, have allowed us to construct the first seismically constrained models for the elastic properties of Mars' core. We observe core-transiting seismic phase SKS from two farside seismic events detected on Mars and measure the travel times of SKS relative to mantle traversing body waves. SKS travels through the core as a compressional wave, providing information about bulk modulus and density. We perform probabilistic inversions using the core-sensitive relative travel times together with gross geophysical data and travel times from other, more proximal, seismic events to seek the equation of state parameters that best describe the liquid iron-alloy core. Our inversions provide constraints on the velocities in Mars' core and are used to develop the first seismically based estimates of its composition. We show that models informed by our SKS data favor a somewhat smaller (median core radius = 1,780 to 1,810 km) and denser (core density = 6.2 to 6.3 g/cm3) core compared to previous estimates, with a P-wave velocity of 4.9 to 5.0 km/s at the core-mantle boundary, with the composition and structure of the mantle as a dominant source of uncertainty. We infer from our models that Mars' core contains a median of 20 to 22 wt% light alloying elements when we consider sulfur, oxygen, carbon, and hydrogen. These data can be used to inform models of planetary accretion, composition, and evolution.

Seismic evidence for a 1000 km mantle discontinuity under the Pacific.

Seismic discontinuities in the mantle are indicators of its thermo-chemical state and offer clues to its dynamics. Ray-based seismic methods, though limited by the approximations made, have mapped mantle transition zone discontinuities in detail, but have yet to offer definitive conclusions on the presence and nature of mid-mantle discontinuities. Here, we show how to use a wave-equation-based imaging method, reverse-time migration of precursors to surface-reflected seismic body waves, to uncover both mantle transition zone and mid-mantle discontinuities, and interpret their physical nature. We observe a thinned mantle transition zone southeast of Hawaii, and a reduction in impedance contrast around 410 km depth in the same area, suggesting a hotter-than-average mantle in the region. Here, we furthermore reveal a 4000-5000 km-wide reflector in new images of the mid mantle below the central Pacific, at 950-1050 km depth. This deep discontinuity exhibits strong topography and generates reflections with polarity opposite to those originating at the 660 km discontinuity, implying an impedance reversal near 1000 km. We link this mid-mantle discontinuity to the upper reaches of deflected mantle plumes upwelling in the region. Reverse-time migration full-waveform imaging is a powerful approach to imaging Earth's interior, capable of broadening our understanding of its structure and dynamics and shrinking modeling uncertainties.

Seismic detection of a deep mantle discontinuity within Mars by InSight.

Constraining the thermal and compositional state of the mantle is crucial for deciphering the formation and evolution of Mars. Mineral physics predicts that Mars' deep mantle is demarcated by a seismic discontinuity arising from the pressure-induced phase transformation of the mineral olivine to its higher-pressure polymorphs, making the depth of this boundary sensitive to both mantle temperature and composition. Here, we report on the seismic detection of a midmantle discontinuity using the data collected by NASA's InSight Mission to Mars that matches the expected depth and sharpness of the postolivine transition. In five teleseismic events, we observed triplicated P and S waves and constrained the depth of this discontinuity to be 1,006 [Formula: see text] 40 km by modeling the triplicated waveforms. From this depth range, we infer a mantle potential temperature of 1,605 [Formula: see text] 100 K, a result consistent with a crust that is 10 to 15 times more enriched in heat-producing elements than the underlying mantle. Our waveform fits to the data indicate a broad gradient across the boundary, implying that the Martian mantle is more enriched in iron compared to Earth. Through modeling of thermochemical evolution of Mars, we observe that only two out of the five proposed composition models are compatible with the observed boundary depth. Our geodynamic simulations suggest that the Martian mantle was relatively cold 4.5 Gyr ago (1,720 to 1,860 K) and are consistent with a present-day surface heat flow of 21 to 24 mW/m2.

Seismic detection of the martian core.

Clues to a planet's geologic history are contained in its interior structure, particularly its core. We detected reflections of seismic waves from the core-mantle boundary of Mars using InSight seismic data and inverted these together with geodetic data to constrain the radius of the liquid metal core to 1830 ± 40 kilometers. The large core implies a martian mantle mineralogically similar to the terrestrial upper mantle and transition zone but differing from Earth by not having a bridgmanite-dominated lower mantle. We inferred a mean core density of 5.7 to 6.3 grams per cubic centimeter, which requires a substantial complement of light elements dissolved in the iron-nickel core. The seismic core shadow as seen from InSight's location covers half the surface of Mars, including the majority of potentially active regions-e.g., Tharsis-possibly limiting the number of detectable marsquakes.

Core Formation and Geophysical Properties of Mars.

The chemical and physical properties of the interiors of terrestrial planets are largely determined during their formation and differentiation. Modeling a planet's formation provides important insights into the properties of its core and mantle, and conversely, knowledge of those properties may constrain formational narratives. Here, we present a multi-stage model of Martian core formation in which we calculate core-mantle equilibration using parameterizations from high pressure-temperature metal-silicate partitioning experiments. We account for changing core-mantle boundary (CMB) conditions, composition-dependent partitioning, and partial equilibration of metal and silicate, and we evolve oxygen fugacity (fO2) self-consistently. The model successfully reproduces published meteorite-based estimates of most elemental abundances in the bulk silicate Mars, which can be used to estimate core formation conditions and core composition. This composition implies that the primordial material that formed Mars was significantly more oxidized (0.9-1.4 log units below the iron-wüstite buffer) than that of the Earth, and that core-mantle equilibration in Mars occurred at 42-60% of the evolving CMB pressure. On average, at least 84% of accreted metal and at least 40% of the mantle were equilibrated in each impact, a significantly higher degree of metal equilibration than previously reported for the Earth. In agreement with previous studies, the modeled Martian core is rich in sulfur (18-19 wt%), with less than one weight percent O and negligible Si. We have used these core and mantle compositions to produce physical models of the present-day Martian interior and evaluate the sensitivity of core radius to crustal thickness, mantle temperature, core composition, core temperature, and density of the core alloy. Trade-offs in how these properties affect observable physical parameters like planetary mass, radius, moment of inertia, and tidal Love number k2 define a range of likely core radii: 1620-1870 km. Seismic velocity profiles for several combinations of model parameters have been used to predict seismic body-wave travel times and planetary normal mode frequencies. These results may be compared to forthcoming Martian seismic data to further constrain core formation conditions and geophysical properties.

Global quieting of high-frequency seismic noise due to COVID-19 pandemic lockdown measures.

Human activity causes vibrations that propagate into the ground as high-frequency seismic waves. Measures to mitigate the coronavirus disease 2019 (COVID-19) pandemic caused widespread changes in human activity, leading to a months-long reduction in seismic noise of up to 50%. The 2020 seismic noise quiet period is the longest and most prominent global anthropogenic seismic noise reduction on record. Although the reduction is strongest at surface seismometers in populated areas, this seismic quiescence extends for many kilometers radially and hundreds of meters in depth. This quiet period provides an opportunity to detect subtle signals from subsurface seismic sources that would have been concealed in noisier times and to benchmark sources of anthropogenic noise. A strong correlation between seismic noise and independent measurements of human mobility suggests that seismology provides an absolute, real-time estimate of human activities.

Inferring Earth's discontinuous chemical layering from the 660-kilometer boundary topography.

Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth's upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle.

Seismically determined elastic parameters for Earth's outer core.

Turbulent convection of the liquid iron alloy outer core generates Earth's magnetic field and supplies heat to the mantle. The exact composition of the iron alloy is fundamentally linked to the processes powering the convection and can be constrained by its seismic properties. Discrepancies between seismic models determined using body waves and normal modes show that these properties are not yet fully agreed upon. In addition, technical challenges in experimentally measuring the equation-of-state (EoS) parameters of liquid iron alloys at high pressures and temperatures further complicate compositional inferences. We directly infer EoS parameters describing Earth's outer core from normal mode center frequency observations and present the resulting Elastic Parameters of the Outer Core (EPOC) seismic model. Unlike alternative seismic models, ours requires only three parameters and guarantees physically realistic behavior with increasing pressure for a well-mixed homogeneous material along an isentrope, consistent with the outer core's condition. We show that EPOC predicts available normal mode frequencies better than the Preliminary Reference Earth Model (PREM) while also being more consistent with body wave-derived models, eliminating a long-standing discrepancy. The velocity at the top of the outer core is lower, and increases with depth more steeply, in EPOC than in PREM, while the density in EPOC is higher than that in PREM across the outer core. The steeper profiles and higher density imply that the outer core comprises a lighter but more compressible alloy than that inferred for PREM. Furthermore, EPOC's steeper velocity gradient explains differential SmKS body wave travel times better than previous one-dimensional global models, without requiring an anomalously slow ~90- to 450-km-thick layer at the top of the outer core.

Regional variation of inner core anisotropy from seismic normal mode observations.

Earth's solid inner core is surrounded by a convecting liquid outer core, creating the geodynamo driving the planet's magnetic field. Seismic studies using compressional body waves suggest hemispherical variation in the anisotropic structure of the inner core, but are poorly constrained because of limited earthquake and receiver distribution. Here, using normal mode splitting function measurements from large earthquakes, based on extended cross-coupling theory, we observe both regional variations and eastern versus western hemispherical anisotropy in the inner core. The similarity of this pattern with Earth's magnetic field suggests freezing-in of crystal alignment during solidification or texturing by Maxwell stress as origins of the anisotropy. These observations limit the amount of inner core super rotation, but would be consistent with oscillation.