Search

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.

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.

SUMMARY: The Earth's mantle transition zone (MTZ) plays a key role in the thermal and compositional interactions between the upper and lower mantle. Seismic anisotropy provides useful information about mantle deformation and dynamics across the MTZ. However, seismic anisotropy in the MTZ is difficult to constrain from surface wave or shear wave splitting measurements. Here, we investigate the sensitivity to anisotropy of a body wave method, SS precursors, through 3-D synthetic modelling and apply it to real data. Our study shows that the SS precursors can distinguish the anisotropy originating from three depths: shallow upper mantle (80–220 km), deep upper mantle above 410 km, and MTZ (410–660 km). Synthetic resolution tests indicate that SS precursors can resolve $\ge $3 per cent azimuthal anisotropy where data have an average signal-to-noise ratio (SNR = 7) and sufficient azimuthal coverage. To investigate regional sensitivity, we apply the stacking and inversion methods to two densely sampled areas: the Japan subduction zone and a central Pacific region around the Hawaiian hotspot. We find evidence for significant VS anisotropy (15.3 ± 9.2 per cent) with a trench-perpendicular fast direction (93° ± 5°) in the MTZ near the Japan subduction zone. We attribute the azimuthal anisotropy to the grain-scale shape-preferred orientation of basaltic materials induced by the shear deformation within the subducting slab beneath NE China. In the central Pacific study region, there is a non-detection of MTZ anisotropy, although modelling suggests the data coverage should allow us to resolve at least 3 per cent anisotropy. Therefore, the Hawaiian mantle plume has not produced detectable azimuthal anisotropy in the MTZ.

The Seismic Experiment for Interior Structure (SEIS) of the InSight mission to Mars, has been providing direct information on Martian interior structure and dynamics of that planet since it landed. Compared to seismic recordings on Earth, ground motion measurements acquired by SEIS on Mars are made under dramatically different ambient noise conditions, but include idiosyncratic signals that arise from coupling between different InSight sensors and spacecraft components. This work is to synthesize what is known about these signal types, illustrate how they can manifest in waveforms and noise correlations, and present pitfalls in structural interpretations based on standard seismic analysis methods. We show that glitches, a type of prominent transient signal, can produce artifacts in ambient noise correlations. Sustained signals that vary in frequency, such as lander modes which are affected by variations in temperature and wind conditions over the course of the Martian Sol, can also contaminate ambient noise results. Therefore, both types of signals have the potential to bias interpretation in terms of subsurface layering. We illustrate that signal processing in the presence of identified nonseismic signals must be informed by an understanding of the underlying physical processes in order for high fidelity waveforms of ground motion to be extracted. While the origins of most idiosyncratic signals are well understood, the 2.4 Hz resonance remains debated and the literature does not contain an explanation of its fine spectral structure. Even though the selection of idiosyncratic signal types discussed in this paper may not be exhaustive, we provide guidance on best practices for enhancing the robustness of structural interpretations.

SUMMARY: We investigated the likelihood of radial anisotropy in the shallow and deep upper mantle, including the mantle transition zone (MTZ) under the Indian Ocean. Seismic anisotropy can be an indicator of mantle deformation through lattice preferred orientation of anisotropic crystals in the mantle. It has thus the potential to illuminate Earth’s dynamic interior, but previous seismic tomography studies have not achieved consensus on the existence of radial anisotropy below ∼250 km depth. We developed a fully nonlinear transdimensional hierarchical Bayesian Markov Chain Monte Carlo approach to invert fundamental and higher mode surface wave dispersion data and applied it to a subset of a global Love and Rayleigh wave data set. We obtained posterior model parameter distributions for shear wave velocity (VS) and radial anisotropy ξ under the Indian Ocean. These posterior model distributions were used to calculate the probability of having radial anisotropy at different depths. We demonstrated that separate inversions of Love and Rayleigh waves yield models compatible with the results of joint inversions within uncertainties. The obtained pattern of VS anomalies agrees with most previous studies. They display negative anomalies along ridges in the uppermost mantle, but those are stronger than for regularized inversions. The Central Indian Ridge and the Southeastern Indian Ridge present velocity anomalies that extend to ∼200 km depth, whereas the Southwestern Indian Ridge seems to have a shallower origin. Weaker, laterally variable velocity perturbations were found at larger depths. The anisotropy models differ more strongly from regularized inversion results and their uncertainties were rather large. We found that anisotropy models from regularized inversions also depend on the chosen parametrization, which is consistent with the existence of a large model null-space. Apart from a fast horizontally polarized shear wave signal in the top 100 km, likely reflecting the horizontal plate motion due to asthenospheric deformation, no clear relation to surface geology was detected in the anisotropy models. We found that, although the anisotropy model uncertainties are rather large, and lateral variations are present, the data generally prefer at least 1 per cent anisotropy in the MTZ with fast vertically polarized shear waves, within errors. Incorporating group velocity data did not help better constrain deep structure by reducing parameter trade-offs. We also tested the effect of prior constraints on the 410- and 660-km topography and found that the undulations of these discontinuities had little effect on the resulting models in our study region.

The geophysics of extrasolar planets is a scientific topic often regarded asstanding largely beyond the reach of near-term observations. This reality in noway diminishes the central role of geophysical phenomena in shaping planetaryoutcomes, from formation, to thermal and chemical evolution, to numerous issuesof surface and near-surface habitability. We emphasize that for a balancedunderstanding of extrasolar planets, it is important to look beyond the naturalbiases of current observing tools, and actively seek unique pathways tounderstand exoplanet interiors as best as possible during the long interimprior to a time when internal components are more directly accessible. Suchpathways include but are not limited to: (a) enhanced theoretical and numericalmodeling, (b) laboratory research on critical material properties, (c)measurement of geophysical properties by indirect inference from imprints lefton atmospheric and orbital properties, and (d) the purpose-driven use of SolarSystem object exploration expressly for its value in comparative planetologytoward exoplanet-analogs. Breaking down barriers that envision local SolarSystem exploration, including the study of Earth's own deep interior, asseparate from and in financial competition with extrasolar planet research, maygreatly improve the rate of needed scientific progress for exoplanetgeophysics. As the number of known rocky and icy exoplanets grows in the yearsahead, we expect demand for expertise in 'exogeoscience' will expand at acommensurately intense pace. We highlight key topics, including: how wateroceans below ice shells may dominate the total habitability of our galaxy byvolume, how free-floating nomad planets may often attain habitable subsurfaceoceans supported by radionuclide decay, and how deep interiors may criticallyinteract with atmospheric mass loss via dynamo-driven magnetic fields.

The boundary between the lithosphere and asthenosphere is associated with a plate-wide high seismic velocity “lid” overlying lowered velocities, consistent with thermal models. Seismic body waves also intermittently detect a sharp velocity reduction at similar depths, the Gutenberg (G) discontinuity, which cannot be explained by temperature alone. We compared an anisotropic tomography model with detections of the G to evaluate their context and relation to the lithosphere-asthenosphere boundary (LAB). We find that the G is primarily associated with vertical changes in azimuthal anisotropy and lies above a thermally controlled LAB, implying the two are not equivalent interfaces. The origin of the G is a result of frozen-in lithospheric structures, regional compositional variations of the mantle, or dynamically perturbed LAB.

The boundary between the lithosphere and asthenosphere is associated with a platewide high-seismic velocity "lid" overlying lowered velocities, consistent with thermal models. Seismic body waves also intermittently detect a sharp velocity reduction at similar depths, the Gutenberg (G) discontinuity, which cannot be explained by temperature alone. We compared an anisotropic tomography model with detections of the G to evaluate their context and relation to the lithosphere-asthenosphere boundary (LAB). We find that the G is primarily associated with vertical changes in azimuthal anisotropy and lies above a thermally controlled LAB, implying that the two are not equivalent interfaces. The origin of the G is a result of frozen-in lithospheric structures, regional compositional variations of the mantle, or dynamically perturbed LAB.

A conceptual model, based on field observations and assumed physics of a perennial firn aquifer near Helheim Glacier (southeast Greenland), is evaluated via steady-state 2-D simulation of liquid water flow and energy transport with phase change. The simulation approach allows natural representation of flow and energy advection and conduction that occur in vertical meltwater recharge through the unsaturated zone and in lateral flow within the saturated aquifer. Agreement between measured and simulated aquifer geometry, temperature, and recharge and discharge rates confirms that the conceptual field-data-based description of the aquifer is consistent with the primary physical processes of groundwater flow, energy transport and phase change. Factors that are found to control simulated aquifer configuration include surface temperature, meltwater recharge rate, residual total-water saturation and capillary fringe thickness. Simulation analyses indicate that the size of perennial firn aquifers depends primarily on recharge rates from surface snowmelt. Results also imply that the recent aquifer expansion, likely due to a warming climate, may eventually produce lakes on the ice-sheet surface that would affect the surface energy balance. [ABSTRACT FROM AUTHOR]