Showing papers by "Mark A. Wieczorek published in 2022"

Surface waves and crustal structure on Mars

Abstract: We detected surface waves from two meteorite impacts on Mars. By measuring group velocity dispersion along the impact-lander path, we obtained a direct constraint on crustal structure away from the InSight lander. The crust north of the equatorial dichotomy had a shear wave velocity of approximately 3.2 kilometers per second in the 5- to 30-kilometer depth range, with little depth variation. This implies a higher crustal density than inferred beneath the lander, suggesting either compositional differences or reduced porosity in the volcanic areas traversed by the surface waves. The lower velocities and the crustal layering observed beneath the landing site down to a 10-kilometer depth are not a global feature. Structural variations revealed by surface waves hold implications for models of the formation and thickness of the martian crust. Description An insightful impact On 24 December 2021, the seismometer for the InSight mission on Mars detected a large seismic event with a distinct signature. Posiolova et al. discovered that the event was caused by a meteor impact on the surface of Mars, which was confirmed by satellite observations of a newly formed 150-kilometer crater. The surface nature and size of the impact allowed Kim et al. to detect surface waves from the event, which have yet to be observed on Mars. These surface waves help to untangle the structure of the Martian crust, which has various amounts of volcanic and sedimentary rock, along with subsurface ice, in different regions of the planet (see the Perspective by Yang and Chen). The characteristics of the impact itself are important because they provide a seismic fingerprint of an impact event that is different from the marsquakes observed so far. —BG A new crater formed on the surface of Mars was detected with the seismometer on the InSight mission.

InSight Constraints on the Global Character of the Martian Crust

Abstract: Analyses of seismic data from the InSight mission have provided the first in situ constraints on the thickness of the crust of Mars. These crustal thickness constraints are currently limited to beneath the lander that is located in the northern lowlands, and we use gravity and topography data to construct global crustal thickness models that satisfy the seismic data. These models consider a range of possible mantle and core density profiles, a range of crustal densities, a low‐density surface layer, and the possibility that the crustal density of the northern lowlands is greater than that of the southern highlands. Using the preferred InSight three‐layer seismic model of the crust, the average crustal thickness of the planet is found to lie between 30 and 72 km. Depending on the choice of the upper mantle density, the maximum permissible density of the northern lowlands and southern highlands crust is constrained to be between 2,850 and 3,100 kg m−3. These crustal densities are lower than typical Martian basaltic materials and are consistent with a crust that is on average more felsic than the materials found at the surface. We argue that a substantial portion of the crust of Mars is a primary crust that formed during the initial differentiation of the planet. Various hypotheses for the origin of the observed intracrustal seisimic layers are assessed, with our preferred interpretation including thick volcanic deposits, ejecta from the Utopia basin, porosity closure, and differentiation products of a Borealis impact melt sheet.

Marsquake Locations and 1‐D Seismic Models for Mars From InSight Data

Abstract: We present inversions for the structure of Mars using the first Martian seismic record collected by the InSight lander. We identified and used arrival times of direct, multiples, and depth phases of body waves, for 17 marsquakes to constrain the quake locations and the one‐dimensional average interior structure of Mars. We found the marsquake hypocenters to be shallower than 40 km depth, most of them being located in the Cerberus Fossae graben system, which could be a source of marsquakes. Our results show a significant velocity jump between the upper and the lower part of the crust, interpreted as the transition between intrusive and extrusive rocks. The lower crust makes up a significant fraction of the crust, with seismic velocities compatible with those of mafic to ultramafic rocks. Additional constraints on the crustal thickness from previous seismic analyses, combined with modeling relying on gravity and topography measurements, yield constraints on the present‐day thermochemical state of Mars and on its long‐term history. Our most constrained inversion results indicate a present‐day surface heat flux of 22 ± 1 mW/m2, a relatively hot mantle (potential temperature: 1740 ± 90 K) and a thick lithosphere (540 ± 120 km), associated with a lithospheric thermal gradient of 1.9 ± 0.3 K/km. These results are compatible with recent seismic studies using a reduced data set and different inversion approaches, confirming that Mars' potential mantle temperature was initially relatively cold (1780 ± 50 K) compared to that of its present‐day state, and that its crust contains 10–12 times more heat‐producing elements than the primitive mantle.

Lithospheric Structure of Venusian Crustal Plateaus

Abstract: Crustal plateaus are Venusian highlands characterized by tectonized terrains. It is commonly interpreted that their topography is isostatically supported and that they represent fossils of an extinct tectonic regime. Using gravity and topography we perform a comprehensive investigation of the lithospheric structure of six crustal plateaus. We computed the admittance (gravity to topography wavelength‐dependent ratio) for each region and compared them to modeled admittances. Three compensation scenarios were tested: Airy isostasy, a surface‐loading flexural model, and a flexural model with surface and subsurface loads. Our results show that the topography of most plateaus is supported by crustal thickening and that the addition of a mantle support component is not necessary at the investigated wavelengths. The elastic thickness was constrained to be less than 35 km with a best‐fitting average of 15 km, confirming that these regions are consistent with an isostatic regime. The average crustal thickness of the plateaus ranges from 15 to 34 km, and if they are in Airy isostasy, this implies that the global average crustal thickness of Venus is about 20 km. Phoebe Regio is the sole exception of our analysis in that crustal thicknesses that are compatible with the other plateaus are obtained only when a buoyant layer is included. Heat flow estimations computed from the elastic thickness indicate that the plateaus formed under higher heat flow conditions compared to the current global average and could have caused localized melting. Present‐day heat flow predictions suggest that eclogitization could occur where the crust is thickest.

Distinguishing the Origin of Asteroid (16) Psyche

Abstract: The asteroid (16) Psyche may be the metal-rich remnant of a differentiated planetesimal, or it may be a highly reduced, metal-rich asteroidal material that never differentiated. The NASA Psyche mission aims to determine Psyche's provenance. Here we describe the possible solar system regions of origin for Psyche, prior to its likely implantation into the asteroid belt, the physical and chemical processes that can enrich metal in an asteroid, and possible meteoritic analogs. The spacecraft payload is designed to be able to discriminate among possible formation theories. The project will determine Psyche's origin and formation by measuring any strong remanent magnetic fields, which would imply it was the core of a differentiated body; the scale of metal to silicate mixing will be determined by both the neutron spectrometers and the filtered images; the degree of disruption between metal and rock may be determined by the correlation of gravity with composition; some mineralogy (e.g., modeled silicate/metal ratio, and inferred existence of low-calcium pyroxene or olivine, for example) will be detected using filtered images; and the nickel content of Psyche's metal phase will be measured using the GRNS.

Distinguishing the Origin of Asteroid (16) Psyche

Abstract: The asteroid (16) Psyche may be the metal-rich remnant of a differentiated planetesimal, or it may be a highly reduced, metal-rich asteroidal material that never differentiated. The NASA Psyche mission aims to determine Psyche's provenance. Here we describe the possible solar system regions of origin for Psyche, prior to its likely implantation into the asteroid belt, the physical and chemical processes that can enrich metal in an asteroid, and possible meteoritic analogs. The spacecraft payload is designed to be able to discriminate among possible formation theories. The project will determine Psyche's origin and formation by measuring any strong remanent magnetic fields, which would imply it was the core of a differentiated body; the scale of metal to silicate mixing will be determined by both the neutron spectrometers and the filtered images; the degree of disruption between metal and rock may be determined by the correlation of gravity with composition; some mineralogy (e.g., modeled silicate/metal ratio, and inferred existence of low-calcium pyroxene or olivine, for example) will be detected using filtered images; and the nickel content of Psyche's metal phase will be measured using the GRNS.

Different Martian Crustal Seismic Velocities Across the Dichotomy Boundary From Multi‐Orbiting Surface Waves

Abstract: We have observed both minor‐arc (R1) and major‐arc (R2) Rayleigh waves for the largest marsquake (magnitude of 4.7 ± 0.2) ever recorded. Along the R1 path (in the lowlands), inversion results show that a simple, two‐layer model with an interface located at 21–29 km and an upper crustal shear‐wave velocity of 3.05–3.17 km/s can fit the group velocity measurements. Along the R2 path, observations can be explained by upper crustal thickness models constrained from gravity data and upper crustal shear‐wave velocities of 2.61–3.27 and 3.28–3.52 km/s in the lowlands and highlands, respectively. The shear‐wave velocity being faster in the highlands than in the lowlands indicates the possible existence of sedimentary rocks, and relatively higher porosity in the lowlands.

Structure Along the Martian Dichotomy Constrained by Rayleigh and Love Waves and Their Overtones

Abstract: Using seismic recordings of event S1222a, we measure dispersion curves of Rayleigh and Love waves, including their first overtones, and invert these for shear velocity (VS) and radial anisotropic structure of the Martian crust. The crustal structure along the topographic dichotomy is characterized by a fairly uniform vertically polarized shear velocity (VSV) of 3.17 km/s between ∼5 and 30 km depth, compatible with the previous study by Kim et al. (2022), https://doi.org/10.1126/science.abq7157. Radial anisotropy as large as 12% (VSH > VSV) is required in the crust between 5 and 40 km depth. At greater depths, we observe a large discontinuity near 63 ± 10 km, below which VSV reaches 4.1 km/s. We interpret this velocity increase as the crust‐mantle boundary along the path. Combined gravimetric modeling suggests that the observed average crustal thickness favors the absence of large‐scale density differences across the topographic dichotomy.

An autonomous lunar geophysical experiment package (ALGEP) for future space missions

Abstract: Abstract Geophysical observations will provide key information about the inner structure of the planets and satellites and understanding the internal structure is a strong constraint on the bulk composition and thermal evolution of these bodies. Thus, geophysical observations are a key to uncovering the origin and evolution of the Moon. In this article, we propose the development of an autonomous lunar geophysical experiment package, composed of a suite of instruments and a central station with standardized interface, which can be installed on various future lunar missions. By fixing the interface between instruments and the central station, it would be possible to easily configure an appropriate experiment package for different missions. We describe here a series of geophysical instruments that may be included as part of the geophysical package: a seismometer, a magnetometer, a heat flow probe, and a laser reflector. These instruments will provide mechanical, thermal, and geodetic parameters of the Moon that are strongly related to the internal structure. We discuss the functionality required for future geophysical observations of the Moon, including the development of the central station that will be used commonly by different payloads.

Interior Dynamics and Thermal Evolution of Mars -- a Geodynamic Perspective

Abstract: Over the past decades, global geodynamic models have been used to investigate the thermal evolution of terrestrial planets. With the increase of computational power and improvement of numerical techniques, these models have become more complex, and simulations are now able to use a high resolution 3D spherical shell geometry and to account for strongly varying viscosity, as appropriate for mantle materials. In this study we review global 3D geodynamic models that have been used to study the thermal evolution and interior dynamics of Mars. We discuss how these models can be combined with local and global observations to constrain the planet’s thermal history. In particular, we use the recent InSight estimates of the crustal thickness, upper mantle structure, and core size to show how these constraints can be combined with 3D geodynamic models to improve our understanding of the interior dynamics, present-day thermal state and temperature variations in the interior of Mars. Our results show that the crustal thickness variations control the surface heat flow and the elastic thickness pattern, as well as the location of melting zones in the present-day martian mantle. The lithospheric temperature and the seismic velocities pattern in the shallow mantle reflect the crustal thickness pattern. The large size of the martian core leads to a smaller scale convection pattern in the mantle than previously suggested. Strong mantle plumes that produce melt up to recent times become focused in Tharsis and Elysium, while weaker plumes are distributed throughout the mantle. The thickness of the seismogenic layer, where seismic events can occur, can be used to discriminate between geodynamic models, if the source depth and location of seismic events is known. Furthermore model predictions of present-day martian seismicity can be compared to the values measured by InSight. Future models need to consider recent estimates from the present-day elastic lithosphere thickness at the north pole of Mars, the effects of lateral variations of seismic velocities on waves propagation through the mantle and lithosphere, and to test the spatial distribution of seismicity by comparing model predictions to observations.

Crustal Structure Constraints From the Detection of the SsPp Phase on Mars

Abstract: The shallowest intracrustal layer (extending to 8 ± 2 km depth) beneath the Mars InSight Lander site exhibits low seismic wave velocity, which is likely related to a combination of high porosity and other lithological factors. The SsPp phase, an SV‐ to P‐wave reflection on the receiver side, is naturally suited for constraining the seismic structure of this top crustal layer since its prominent signal makes it observable with a single station without the need for stacking. We have analyzed six broadband and low‐frequency seismic events recorded on Mars and made the first coherent detection of the SsPp phase on the red planet. The timing and amplitude of SsPp confirm the existence of the ∼8 km interface in the crust and the large wave speed (or impedance) contrast across it. With our new constraints from the SsPp phase, we determined that the average P‐wave speed in the top crustal layer is between 2.5 and 3.2 km/s, which is a more precise and robust estimate than the previous range of 2.0–3.5 km/s obtained by receiver function analysis. The low velocity of Layer 1 likely results from the presence of relatively low‐density lithified sedimentary rocks and/or aqueously altered igneous rocks that also have a significant amount of porosity, possibly as much as 22%–30% by volume (assuming an aspect ratio of 0.1 for the pore space). These porosities and average P‐wave speeds are compatible with our current understanding of the upper crustal stratigraphy beneath the InSight Lander site.

The Psyche Gravity Investigation

Abstract: Abstract The objective of the NASA Psyche mission gravity science investigation is to map the mass distribution within asteroid (16) Psyche to elucidate interior structure and to resolve the question of whether this metal-rich asteroid represents a remnant metal core or whether it is a primordial body that never melted. Measurements of gravity will be obtained via the X-band telecommunication system on the Psyche spacecraft, collected from progressively lower mapping altitudes. Orbital gravity will allow an estimate of $GM$ G M to better than 0.001 km 3 s −2 . A spherical harmonic model of gravity to degree and order 10 will be achievable and, in concert with spherical harmonic data sets from topography and magnetometry, as well as surface composition data, will provide information regarding the spatial and radial distribution of mass that will be used to constrain the origin and evolution of (16) Psyche.

The Psyche Topography and Geomorphology Investigation

Abstract: Abstract Detailed mapping of topography is crucial for the understanding of processes shaping the surfaces of planetary bodies. In particular, stereoscopic imagery makes a major contribution to topographic mapping and especially supports the geologic characterization of planetary surfaces. Image data provide the basis for extensive studies of the surface structure and morphology on local, regional and global scales using photogeologic information from images, the topographic information from stereo-derived digital terrain models and co-registered spectral terrain information from color images. The objective of the Psyche topography and geomorphology investigation is to derive the detailed shape of (16) Psyche to generate orthorectified image mosaics, which are needed to study the asteroids’ landforms, interior structure, and the processes that have modified the surface over geologic time. In this paper we describe our approaches for producing shape models, and our plans for acquiring requested image data to quantify the expected accuracy of the results. Multi-angle images obtained by Psyche’s camera will be used to create topographic models with about 15 m/pixel horizontal resolution and better than 10 m height accuracy on a global scale. This is slightly better as global imaging obtained during the Dawn mission, however, both missions yield resolutions of a few m/pixel locally. Two different techniques, stereophotogrammetry and stereophotoclinometry, are used to model the shape; these models will be merged with the gravity fields obtained by the Psyche spacecraft to produce geodetically controlled topographic models. The resulting digital topography models, together with the gravity data, will reveal the tectonic, volcanic, impact, and gradational history of Psyche, and enable co-registration of data sets to determine Psyche’s geologic history.

A Positive Feedback Between Crustal Thickness and Melt Extraction for the Origin of the Martian Dichotomy

Abstract: A North/South difference in crustal thickness is likely at the origin of the Martian dichotomy in topography. Recent crustal thickness maps were obtained by inversion of topography and gravity data seismically anchored at the InSight station. On average, the Martian crust is 51–71 km thick with a southern crust thicker by 18–28 km than the northern one. The origin of this crustal dichotomy is still debated although the hypothesis of a large impact is at present very popular. Here, we propose a new mechanism for the formation of this dichotomy that involves a positive feedback between crustal growth and mantle melting. As the crust is enriched in heat‐producing elements, the lid of a one‐plate planet is hotter and thinner where the crust is thicker, inducing a larger amount of partial melt below the lid and hence a larger rate of melt extraction and crustal growth. We first demonstrate analytically that larger wavelength perturbations, that is, hemispherical perturbations, grow faster because smaller wavelengths are more attenuated by thermal diffusion. We then use a parameterized thermal evolution model with a well‐mixed mantle topped by two different lids characterized by their thermal structures and thicknesses to study the growth of the crust in the two hemispheres. Our results demonstrate that this positive feedback can generate a significant crustal dichotomy.