Inversion of Meteor Rayleigh Waves on Earth and Modeling of Air Coupled Rayleigh Waves on Mars

Meteor impacts and/or meteor events generate body and surface seismic waves on the surface of a planet. When meteoroids burst in the atmosphere, they generate shock waves that subsequently convert into acoustic waves in the atmosphere and seismic waves in the ground. This effect can be modeled as the amplitude of Rayleigh and other Spheroidal modes excitation, due to atmospheric/ground coupling effects.First, an inversion of the seismic source of Chelyabinsk superbolide is performed. We develop an approach in order to model a line source in the atmosphere, corresponding to the consecutive generation of shock waves by the interaction with the atmosphere. The model is based on the known trajectory. We calculate the synthetic seismograms of Rayleigh waves associated with the event by the summation of normal modes of a model of the solid part and the atmosphere of the planet. Through an inversion technique based on singular value decomposition, we perform a full Rayleigh wave inversion and we provide solutions for the moment magnitude.SEIS will likely detect seismic waves generated by impacts and the later might be further located by remote sensing differential processing. In the case of Mars, we use the same method to obtain waveforms associated with impacts on the planetary surface or in low altitudes in the Martian atmosphere. We show that the contribution of the fundamental spheroidal solid mode is dominating the waveforms, compared to that of the first two overtones. We perform an amplitude comparison and we show that small impactors (diameter of 0.5 to 2 m), can be detected by the SEIS VBB seismometer of InSight mission, even in short epicentral distances, in the higher frequencies of the Rayleigh waves. We perform an analysis based on impact rate estimations and we calculate the number of detectable events of 1 meter diameter meteor impacts to be 6.7 to 13.4 per 1 Martian year for a \(Q=500\).     

Experimental Investigation of InSight HP<Superscript>3</Superscript> Mole Interaction with Martian Regolith Simulant

The InSight mission launches in 2018 to characterize several geophysical quantities on Mars, including the heat flow from the planetary interior. This quantity will be calculated by utilizing measurements of the thermal conductivity and the thermal gradient down to 5 meters below the Martian surface. One of the components of InSight is the Mole, which hammers into the Martian regolith to facilitate these thermal property measurements. In this paper, we experimentally investigated the effect of the Mole’s penetrating action on regolith compaction and mechanical properties. Quasi-static and dynamic experiments were run with a 2D model of the 3D cylindrical mole. Force resistance data was captured with load cells. Deformation information was captured in images and analyzed using Digitial Image Correlation (DIC). Additionally, we used existing approximations of Martian regolith thermal conductivity to estimate the change in the surrounding granular material’s thermal conductivity due to the Mole’s penetration. We found that the Mole has the potential to cause a high degree of densification, especially if the initial granular material is relatively loose. The effect on the thermal conductivity from this densification was found to be relatively small in first-order calculations though more complete thermal models incorporating this densification should be a subject of further investigation. The results obtained provide an initial estimate of the Mole’s impact on Martian regolith thermal properties.     

Atmospheric Science with InSight

In November 2018, for the first time a dedicated geophysical station, the InSight lander, will be deployed on the surface of Mars. Along with the two main geophysical packages, the Seismic Experiment for Interior Structure (SEIS) and the Heat-Flow and Physical Properties Package (HP³), the InSight lander holds a highly sensitive pressure sensor (PS) and the Temperature and Winds for InSight (TWINS) instrument, both of which (along with the InSight FluxGate (IFG) Magnetometer) form the Auxiliary Sensor Payload Suite (APSS). Associated with the RADiometer (RAD) instrument which will measure the surface brightness temperature, and the Instrument Deployment Camera (IDC) which will be used to quantify atmospheric opacity, this will make InSight capable to act as a meteorological station at the surface of Mars. While probing the internal structure of Mars is the primary scientific goal of the mission, atmospheric science remains a key science objective for InSight. InSight has the potential to provide a more continuous and higher-frequency record of pressure, air temperature and winds at the surface of Mars than previous in situ missions. In the paper, key results from multiscale meteorological modeling, from Global Climate Models to Large-Eddy Simulations, are described as a reference for future studies based on the InSight measurements during operations. We summarize the capabilities of InSight for atmospheric observations, from profiling during Entry, Descent and Landing to surface measurements (pressure, temperature, winds, angular momentum), and the plans for how InSight’s sensors will be used during operations, as well as possible synergies with orbital observations. In a dedicated section, we describe the seismic impact of atmospheric phenomena (from the point of view of both “noise” to be decorrelated from the seismic signal and “signal” to provide information on atmospheric processes). We discuss in this framework Planetary Boundary Layer turbulence, with a focus on convective vortices and dust devils, gravity waves (with idealized modeling), and large-scale circulations. Our paper also presents possible new, exploratory, studies with the InSight instrumentation: surface layer scaling and exploration of the Monin-Obukhov model, aeolian surface changes and saltation / lifing studies, and monitoring of secular pressure changes. The InSight mission will be instrumental in broadening the knowledge of the Martian atmosphere, with a unique set of measurements from the surface of Mars.     

Modeling of Ground Deformation and Shallow Surface Waves Generated by Martian Dust Devils and Perspectives for Near-Surface Structure Inversion

We investigated the possible seismic signatures of dust devils on Mars, both at long and short period, based on the analysis of Earth data and on forward modeling for Mars. Seismic and meteorological data collected in the Mojave Desert, California, recorded the signals generated by dust devils. In the 10–100 s band, the quasi-static surface deformation triggered by pressure fluctuations resulted in detectable ground-tilt effects: these are in good agreement with our modeling based on Sorrells’ theory. In addition, high-frequency records also exhibit a significant excitation in correspondence to dust devil episodes. Besides wind noise, this signal includes shallow surface waves due to the atmosphere-surface coupling and is used for a preliminary inversion of the near-surface S-wave profile down to 50 m depth. In the case of Mars, we modeled the long-period signals generated by the pressure field resulting from turbulence-resolving Large-Eddy Simulations. For typical dust-devil-like vortices with pressure drops of a couple Pascals, the corresponding horizontal acceleration is of a few nm/s² for rocky subsurface models and reaches 10–20 nm/s² for weak regolith models. In both cases, this signal can be detected by the Very-Broad Band seismometers of the InSight/SEIS experiment up to a distance of a few hundred meters from the vortex, the amplitude of the signal decreasing as the inverse of the distance. Atmospheric vortices are thus expected to be detected at the InSight landing site; the analysis of their seismic and atmospheric signals could lead to additional constraints on the near-surface structure, more precisely on the ground compliance and possibly on the seismic velocities.     

Michelson Interferometer for Global High-Resolution Thermospheric Imaging (MIGHTI): Instrument Design and Calibration

The Michelson Interferometer for Global High-resolution Thermospheric Imaging (MIGHTI) instrument was built for launch and operation on the NASA Ionospheric Connection Explorer (ICON) mission. The instrument was designed to measure thermospheric horizontal wind velocity profiles and thermospheric temperature in altitude regions between 90 km and 300 km, during day and night. For the wind measurements it uses two perpendicular fields of view pointed at the Earth’s limb, observing the Doppler shift of the atomic oxygen red and green lines at 630.0 nm and 557.7 nm wavelength. The wavelength shift is measured using field-widened, temperature compensated Doppler Asymmetric Spatial Heterodyne (DASH) spectrometers, employing low order échelle gratings operating at two different orders for the different atmospheric lines. The temperature measurement is accomplished by a multichannel photometric measurement of the spectral shape of the molecular oxygen A-band around 762 nm wavelength. For each field of view, the signals of the two oxygen lines and the A-band are detected on different regions of a single, cooled, frame transfer charge coupled device (CCD) detector. On-board calibration sources are used to periodically quantify thermal drifts, simultaneously with observing the atmosphere. The MIGHTI requirements, the resulting instrument design and the calibration are described.     

Towards a Global Unified Model of Europa’s Tenuous Atmosphere

Despite the numerous modeling efforts of the past, our knowledge on the radiation-induced physical and chemical processes in Europa’s tenuous atmosphere and on the exchange of material between the moon’s surface and Jupiter’s magnetosphere remains limited. In lack of an adequate number of in situ observations, the existence of a wide variety of models based on different scenarios and considerations has resulted in a fragmentary understanding of the interactions of the magnetospheric ion population with both the moon’s icy surface and neutral gas envelope. Models show large discrepancy in the source and loss rates of the different constituents as well as in the determination of the spatial distribution of the atmosphere and its variation with time. The existence of several models based on very different approaches highlights the need of a detailed comparison among them with the final goal of developing a unified model of Europa’s tenuous atmosphere. The availability to the science community of such a model could be of particular interest in view of the planning of the future mission observations (e.g., ESA’s JUpiter ICy moons Explorer (JUICE) mission, and NASA’s Europa Clipper mission). We review the existing models of Europa’s tenuous atmosphere and discuss each of their derived characteristics of the neutral environment. We also discuss discrepancies among different models and the assumptions of the plasma environment in the vicinity of Europa. A summary of the existing observations of both the neutral and the plasma environments at Europa is also presented. The characteristics of a global unified model of the tenuous atmosphere are, then, discussed. Finally, we identify needed future experimental work in laboratories and propose some suitable observation strategies for upcoming missions.     

Daytime Ionosphere Retrieval Algorithm for the Ionospheric Connection Explorer (ICON)

The NASA Ionospheric Connection Explorer Extreme Ultraviolet spectrograph, ICON EUV, will measure altitude profiles of the daytime extreme-ultraviolet (EUV) OII emission near 83.4 and 61.7 nm that are used to determine density profiles and state parameters of the ionosphere. This paper describes the algorithm concept and approach to inverting these measured OII emission profiles to derive the associated \(\mathrm{O}^{+}\) density profile from 150–450 km as a proxy for the electron content in the F-region of the ionosphere. The algorithm incorporates a bias evaluation and feedback step, developed at the U.S. Naval Research Laboratory using data from the Special Sensor Ultraviolet Limb Imager (SSULI) and the Remote Atmospheric and Ionospheric Detection System (RAIDS) missions, that is able to effectively mitigate the effects of systematic instrument calibration errors and inaccuracies in the original photon source within the forward model. Results are presented from end-to-end simulations that convolved simulated airglow profiles with the expected instrument measurement response to produce profiles that were inverted with the algorithm to return data products for comparison to truth. Simulations of measurements over a representative ICON orbit show the algorithm is able to reproduce hmF2 values to better than 5 km accuracy, and NmF2 to better than 12% accuracy over a 12-second integration, and demonstrate that the ICON EUV instrument and daytime ionosphere algorithm can meet the ICON science objectives which require 20 km vertical resolution in hmF2 and 18% precision in NmF2.     

Crustal S-Wave Velocity from Apparent Incidence Angles: A Case Study in Preparation for InSight

Retrieval of crustal structure and thickness of Mars is among the main goals of InSight. Here we investigate which constraints on the crust at the landing site can be provided by apparent P-wave incidence angles derived from P-receiver functions. We consider receiver functions for six different Mars models, calculated from synthetic seismograms generated via Instaseis from the Green’s function databases of the Marsquake Service, in detail. To allow for a larger range of crustal thicknesses and structures, we additionally analyze data from five broad-band stations across Central Europe. We find that the likely usable epicentral distance range for P-wave receiver functions on Mars lies between \(35^{\circ}\) and the core shadow, and can be extended to more than \(150^{\circ}\) by also using the PP-phase. Comparison to models for the spatial distribution of Martian seismicity indicates that sufficient seismicity should occur within the P-wave distance range around InSight within the nominal mission duration to allow for the application of our method. Apparent P-wave incidence angles are derived from the amplitudes of vertical and radial receiver functions at the P-wave onset within a range of period bands, up to 120 s. The apparent incidence angles are directly related to apparent S-wave velocities, which are inverted for the subsurface S-wave velocity structure via a grid search. The veracity of the forward calculated receiver functions and apparent S-wave velocities is ensured by benchmarking various algorithms against the Instaseis synthetics. Results indicate that apparent S-wave velocity curves provide valuable constraints on crustal thickness and structure, even without any additional constraints, and considering the location uncertainty and limited data quantity of InSight. S-wave velocities in the upper half of the crust are constrained best, but if reliable measurements at long periods are available, the curves also provide constraints down to the uppermost mantle. Besides, it is demonstrated that the apparent velocity curves can differentiate between crustal velocity models that are indistinguishable by other methods.