Analysis of regolith properties using seismic signals generated by InSight’s HP3 Penetrator

International audienceInSight’s Seismic Experiment for Interior Structure (SEIS) provides a unique and unprecedented opportunity to conduct the first geotechnical survey of the Martian soil by taking advantage of the repeated seismic signals that will be generated by the mole of the Heat Flow and Physical Properties Package (HP3). Knowledge of the elastic properties of the Martian regolith have implications to material strength and can constrain models of water content, and provide context to geological processes and history that have acted on the landing site in western Elysium Planitia. Moreover, it will help to reduce travel-time errors introduced into the analysis of seismic data due to poor knowledge of the shallow subsurface. The challenge faced by the InSight team is to overcome the limited temporal resolution of the sharp hammer signals, which have significantly higher frequency content than the SEIS 100 Hz sampling rate. Fortunately, since the mole propagates at a rate of ∼1 mm per stroke down to 5 m depth, we anticipate thousands of seismic signals, which will vary very gradually as the mole travels.Using a combination of field measurements and modeling we simulate a seismic data set that mimics the InSight HP3-SEIS scenario, and the resolution of the InSight seismometer data. We demonstrate that the direct signal, and more importantly an anticipated reflected signal from the interface between the bottom of the regolith layer and an underlying lava flow, are likely to be observed both by Insight’s Very Broad Band (VBB) seismometer and Short Period (SP) seismometer. We have outlined several strategies to increase the signal temporal resolution using the multitude of hammer stroke and internal timing information to stack and interpolate multiple signals, and demonstrated that in spite of the low resolution, the key parameters—seismic velocities and regolith depth—can be retrieved with a high degree of confidence

Planned Products of the Mars Structure Service for the InSight Mission to Mars

Abstract The InSight lander will deliver geophysical instruments to Mars in 2018, including seismometers installed directly on the surface (Seismic Experiment for Interior Structure, SEIS). Routine operations will be split into two services, the Mars Structure Service(MSS) and Marsquake Service (MQS), which will be responsible, respectively, for defining the structure models and seismicity catalogs from the mission. The MSS will deliver a series of products before the landing, during the operations, and finally to the Planetary Data System (PDS) archive. Prior to the mission, we assembled a suite of a priori models of Mars, based on estimates of bulk composition and thermal profiles. Initial models during the mission will rely on modeling surface waves and impact-generated body waves independent of prior knowledge of structure. Later modeling will include simultaneous inversion of seismic observations for source and structural parameters. We use Bayesian inversion techniques to obtain robust probability distribution functions of interior structure parameters. Shallow structure will be characterized using the hammering of the heatflow probe mole, as well as measurements of surface wave ellipticity. Crustal scale structure will be constrained by measurements of receiver function and broadband Rayleigh wave ellipticity measurements. Core interacting body wave phases should be observable above modeled martian noise levels, allowing us to constrain deep structure. Normal modes of Mars should also be observable and can be used to estimate the globally averaged 1D structure, while combination with results from the InSight radio science mission and orbital observations will allow for constraint of deeper structure

Double-difference measurements in global full-waveform inversions

International audienceWe demonstrate double-difference (DD) tomography, a method involving differential measurements between stations, for 2-D and 3-D adjoint inversions based on realistic source–receiver distributions, with a focus on the global scale. We first present 2-D synthetic inversion results using cross-correlation traveltime and L2 waveform difference objective functions. Introducing a weighting term to DD objective functions based on the number of measurement pairs per station speeds up convergence and reduces bias in the final inverted model due to uneven data coverage. We next demonstrate frequency-dependent multitaper DD measurements in a 3-D experiment with real earthquake data by computing global-scale gradients. At the global scale, careful selection of station pairs is required for differential measurements in terms of geographical distance or geological context. In our technique, if no suitable pairs are found for a particular station, the DD measurement reduces to a classical misfit measurement. Furthermore, we compare 2-D and 3-D DD results with those from corresponding conventional misfits. By exploiting previously unused information in the recorded wavefield, DD tomography shows promise for balancing the gradient and speeding up convergence, especially around dense regional seismic networks

Balancing unevenly distributed data in seismic tomography: a global adjoint tomography example

International audienceSUMMARY The uneven distribution of earthquakes and stations in seismic tomography leads to slower convergence of nonlinear inversions and spatial bias in inversion results. Including dense regional arrays, such as USArray or Hi-Net, in global tomography causes severe convergence and spatial bias problems, against which conventional pre-conditioning schemes are ineffective. To save computational cost and reduce model bias, we propose a new strategy based on a geographical weighting of sources and receivers. Unlike approaches based on ray density or the Voronoi tessellation, this method scales to large full-waveform inversion problems and avoids instabilities at the edges of dense receiver or source clusters. We validate our strategy using a 2-D global waveform inversion test and show that the new weighting scheme leads to a nearly twofold reduction in model error and much faster convergence relative to a conventionally pre-conditioned inversion. We implement this geographical weighting strategy for global adjoint tomography

Crustal and Uppermost Mantle Heterogeneities Across the Ailaoshan Red River Shear Zone, SE Tibet: Implications for Cenozoic Magmatic Activity

The Ailaoshan Red River shear zone (ARSZ) was formed in the Mesozoic as a suture zone between the Indochina block and the Yangtze craton. Since the Cenozoic, block extrusion due to the Indo-Asian collision has reactivated the fault zone and caused large-scale shearing. Affected by the Cenozoic orogeny, a large volume of magmatic and metamorphic rocks developed in the ARSZ, forming many orogenic gold deposits. However, the source and the geodynamic process of these magmatic activities are still unclear. To gain a basic understanding of the subsurface magmatic activity, we deployed a dense array of 24 broadband seismic stations across the Daping and Chang'an gold deposits at the southern end of the ARSZ. Receiver function analysis, common conversion point stacking, and a joint inversion of receiver functions and surface wave dispersions are performed to image the detailed structure of the crust and uppermost mantle. Low-velocity zones in the mid-lower crust and thinned lithosphere (∼70 km) are imaged under the ARSZ. The observed subsurface structures are verified by 3D numerical modeling with the SEM-FK method. We speculate that the mantle upwelling caused by lithospheric delamination has provided the main source of the mantle component in the magmatic rocks since ∼35 Ma; afterward, high temperatures produced partial melting in the lower crust, which was emplaced along active shear zones.Published versionThis study was supported by the National Key Research and Development Program of China (2016YFC0600302), the National Natural Science Foundations of China (42174056), the Postgraduate Research & Practice Innovation Program of Jiangsu Province, China (KYCX20_0059)

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity.

The compositional and thermal state of Earths mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earths early history

Global adjoint tomography—model GLAD-M25

Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones

Data and Workflow Management for Exascale Global Adjoint Tomography: Chapter 13

International audienc

Simulations of Seismic Wave Propagation on Mars

International audienceWe present global and regional synthetic seismograms computed for 1D and 3D Mars models based on the spectral-element method. For global simulations, we implemented a radially-symmetric Mars model with a 110 km thick crust (Sohl and Spohn in J. Geophys. Res., Planets 102(E1):1613–1635, 1997). For this 1D model, we successfully benchmarked the 3D seismic wave propagation solver SPECFEM3D_GLOBE (Komatitsch and Tromp in Geophys. J. Int. 149(2):390–412, 2002a; 150(1):303–318, 2002b) against the 2D axisymmetric wave propagation solver AxiSEM (Nissen-Meyer et al. in Solid Earth 5(1):425–445, 2014) at periods down to 10 s. We also present higher-resolution body-wave simulations with AxiSEM down to 1 s in a model with a more complex 1D crust, revealing wave propagation effects that would have been difficult to interpret based on ray theory. For 3D global simulations based on SPECFEM3D_GLOBE, we superimposed 3D crustal thickness variations capturing the distinct crustal dichotomy between Mars’ northern and southern hemispheres, as well as topography, ellipticity, gravity, and rotation. The global simulations clearly indicate that the 3D crust speeds up body waves compared to the reference 1D model, whereas it significantly changes surface waveforms and their dispersive character depending on its thickness. We also perform regional simulations with the solver SES3D (Fichtner et al. Geophys. J. Int. 179:1703–1725, 2009) based on 3D crustal models derived from surface composition, thereby addressing the effects of various distinct crustal features down to 2 s. The regional simulations confirm the strong effects of crustal variations on waveforms. We conclude that the numerical tools are ready for examining more scenarios, including various other seismic models and sources

Simulations of Seismic Wave Propagation on Mars

International audienceWe present global and regional synthetic seismograms computed for 1D and 3D Mars models based on the spectral-element method. For global simulations, we implemented a radially-symmetric Mars model with a 110 km thick crust (Sohl and Spohn in J. Geophys. Res., Planets 102(E1):1613–1635, 1997). For this 1D model, we successfully benchmarked the 3D seismic wave propagation solver SPECFEM3D_GLOBE (Komatitsch and Tromp in Geophys. J. Int. 149(2):390–412, 2002a; 150(1):303–318, 2002b) against the 2D axisymmetric wave propagation solver AxiSEM (Nissen-Meyer et al. in Solid Earth 5(1):425–445, 2014) at periods down to 10 s. We also present higher-resolution body-wave simulations with AxiSEM down to 1 s in a model with a more complex 1D crust, revealing wave propagation effects that would have been difficult to interpret based on ray theory. For 3D global simulations based on SPECFEM3D_GLOBE, we superimposed 3D crustal thickness variations capturing the distinct crustal dichotomy between Mars’ northern and southern hemispheres, as well as topography, ellipticity, gravity, and rotation. The global simulations clearly indicate that the 3D crust speeds up body waves compared to the reference 1D model, whereas it significantly changes surface waveforms and their dispersive character depending on its thickness. We also perform regional simulations with the solver SES3D (Fichtner et al. Geophys. J. Int. 179:1703–1725, 2009) based on 3D crustal models derived from surface composition, thereby addressing the effects of various distinct crustal features down to 2 s. The regional simulations confirm the strong effects of crustal variations on waveforms. We conclude that the numerical tools are ready for examining more scenarios, including various other seismic models and sources