skip to main content

Title: 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 ± 40 km by modeling the triplicated waveforms. From this depth range, we infer a mantle potential temperature of 1,605 ± 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/m 2 .  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In recent years there have been several attempts to make the link between mineral properties and seismic anisotropy in the D’’ region but have yet to reach consensus with regards to the dynamics in lower mantle minerals that could give rise to the observed seismic anisotropy. Here, we aim to provide further constraints on the observed long wavelength shear velocity patterns seen in seismic tomography studies. We introduce a forward model of deformation in a subducting slab as it impacts the core mantle boundary (D’’ layer) and proceeds to upwelling at the edge of a simulated LLSVP. By implementing the most recent results from atomistic modeling and high-pressure deformation experiments coupled with a 3-dimensional geodynamic model, we compare the microstructural evolution of an aggregate with a pyrolytic composition to the macroscopically observed seismic anisotropy of the lowermost mantle. We account for topotaxial relations in the forward and reverse phase transitions of MgSiO3-perovskite (Pv) to post-perovskite (pPv) within the slab as well as explore the effects introduced by partial melting near the CMB. Comparisons in the two leading candidate deformation mechanisms in the post-perovskite phase, (001) and (010), are compared. In this study we find that the reverse transition (pPv to Pv) occurs at a depth which is ~ 150 km deeper than that of the forward transition due to increasing temperature near the CMB providing a varying topography of the D’’ discontinuity. Our model also produces good fits with the isotropic velocities of PREM for the bulk lower mantle. When coupled with temperature and pressure dependent forward and reverse phase transitions, a pPv system with dominant (001) slip provides good correlation with the currently observed VSH fast horizontal (~ 1 – 6%) in D’’ and with VSV consistently fast in upwelling areas. Azimuthal variations along the streamline are also investigated showing a symmetry lower than that of the assumed VTI in D’’ introduced by ‘rolling’ effects near the slab’s edge. The addition of 1% partial melting at the CMB is shown to increase S and P wave anisotropy beneath the slab at the base of upwelling with up to ~2.5 & 4.0% P and S wave reductions respectively compared to the global reference. 
    more » « less
  2. SUMMARY The detailed structure near the 410-km discontinuity provides key constraints of the dynamic interactions between the upper mantle and the lower mantle through the mantle transition zone (MTZ) via mass and heat exchange. Meanwhile, the temperature of the subducting slab, which can be derived from its fast wave speed perturbation, is critical for understanding the mantle dynamics in subduction zones where the slab enters the MTZ. Multipathing, i.e. triplicated, body waves that bottom near the MTZ carry rich information of the 410-km discontinuity structure and can be used to constrain the discontinuity depth and radial variations of wave speeds across it. In this study, we systematically analysed the trade-off between model parameters in triplication studies using synthetic examples. Specifically, we illustrated the necessity of using array-normalized amplitude. Two 1-D depth profiles of the wave speed below the Tatar Strait of Russia in the Kuril subduction zone are obtained. We have observed triplications due to both the 410-km discontinuity and the slab upper surface. And, seismic structures for these two interfaces are simultaneously inverted. Our derived 410-km discontinuity depths for the northern and southern regions are at 420$\pm $15 and 425$\pm $15 km, respectively, with no observable uplift. The slab upper surface is inverted to be located about 50–70 km below the 410-km discontinuity. This location is between the depths of the 1 and 2 per cent P-wave speed perturbation contours of a regional 3-D full-waveform inversion (FWI) model, but we found twice the wave speed perturbation amplitude. A wave speed increase of 3.9–4.6 per cent within the slab, compared to 2.0–2.4 per cent from the 3-D FWI model, is necessary to fit the waveforms with the shortest period of 2 s, indicating that high-frequency waves are required to accurately resolve the detailed structures near the MTZ. 
    more » « less
  3. Comparing compositional models of the terrestrial planets provides insights into physicochemical processes that produced planet-scale similarities and differences. The widely accepted compositional model for Mars assumes Mn and more refractory elements are in CI chondrite proportions in the planet, including Fe, Mg, and Si, which along with O make up >90% of the mass of Mars. However, recent improvements in our understandings on the composition of the solar photosphere and meteorites challenge the use of CI chondrite as an analog of Mars. Here we present an alternative model composition for Mars that avoids such an assumption and is based on data from Martian meteorites and spacecraft observations. Our modeling method was previously applied to predict the Earth’s composition. The model establishes the absolute abundances of refractory lithophile elements in the bulk silicate Mars (BSM) at 2.26 times higher than that in CI carbonaceous chondrites. Relative to this chondritic composition, Mars has a systematic depletion in moderately volatile lithophile elements as a function of their condensation temperatures. Given this finding, we constrain the abundances of siderophile and chalcophile elements in the bulkMars and its core. The Martian volatility trend is consistent with <7 wt% S in its core, which is significantly lower than that assumed in most core models (i.e., >10 wt% S). Furthermore, the occurrence of ringwoodite at the Martian core-mantle boundary might have contributed to the partitioning of O and H into the Martian core. 
    more » « less
  4. Abstract

    To quantitatively convert upper mantle seismic wave speeds measured into temperature, density, composition, and corresponding and uncertainty, we introduce theWhole‐rockInterpretativeSeismicToolboxForUltramaficLithologies (WISTFUL). WISTFUL is underpinned by a database of 4,485 ultramafic whole‐rock compositions, their calculated mineral modes, elastic moduli, and seismic wave speeds over a range of pressure (P) and temperature (T) (P = 0.5–6 GPa,T = 200–1,600°C) using the Gibbs free energy minimization routine Perple_X. These data are interpreted with a toolbox of MATLAB® functions, scripts, and three general user interfaces:WISTFUL_relations, which plots relationships between calculated parameters and/or composition;WISTFUL_geotherms, which calculates seismic wave speeds along geotherms; andWISTFUL_inversion, which inverts seismic wave speeds for best‐fit temperature, composition, and density. To evaluate our methodology and quantify the forward calculation error, we estimate two dominant sources of uncertainty: (a) the predicted mineral modes and compositions, and (b) the elastic properties and mixing equations. To constrain the first source of uncertainty, we compiled 122 well‐studied ultramafic xenoliths with known whole‐rock compositions, mineral modes, and estimatedPTconditions. We compared the observed mineral modes with modes predicted using five different thermodynamic solid solution models. The Holland et al. (2018, solution models best reproduce phase assemblages (∼12 vol. % phase root‐mean‐square error [RMSE]) and estimated wave speeds. To assess the second source of uncertainty, we compared wave speed measurements of 40 ultramafic rocks with calculated wave speeds, finding excellent agreement (VpRMSE = 0.11 km/s). WISTFUL easily analyzes seismic datasets, integrates into modeling, and acts as an educational tool.

    more » « less
  5. Abstract

    We present ab initio (LDA + Usc) studies of high‐temperature and high‐pressure elastic properties of pure as well as iron‐bearing (ferrous, Fe2+, and ferric, Fe3+) and aluminum‐bearing MgSiO3postperovskite, the likely dominant phase in the deep lower mantle of the Earth. Thermal effects are addressed within the quasiharmonic approximation by combining vibrational density of states and static elastic coefficients. Aggregate elastic moduli and sound velocities for the Mg end members are successfully compared with scarce experimental data available. Effects of iron (Fe) and aluminum (Al) substitutions on elastic properties and their pressure and temperature dependence have been thoroughly investigated. At the observed perovskite to postperovskite transition (P = 125 GPa andT = 2,500 K), compressional and shear velocities increase by 0–1% and 1.5–3.75%, respectively. This observation is consistent with some seismic studies of the Ddiscontinuity beneath the Caribbean, which suggests that our robust estimates of elastic properties of the postperovskite phase will be very helpful to understand lateral velocity variations in the deep lower mantle region and to constrain its composition and thermal structure.

    more » « less