skip to main content


Title: Autocorrelation Infrasound Interferometry on Mars
Abstract

A high‐sensitivity pressure sensor was deployed as part of the Mars Interior Exploration using Seismic Investigations, Geodesy and Heat Transport lander on Elysium Planitia in November 2018. We use pressure records from 1 October to 31 December 2019 (Sol 301–389) for frequencies between 0.1 and 0.5 Hz to infer relative sound‐speed changes in the Martian atmosphere using the autocorrelation infrasound interferometry method. We find that relative sound‐speed changes are up to ±15%, follow a similar pattern to Martian‐daily variations of atmospheric temperature and horizontal wind velocity, and are similar to those inferred from in‐situ observations and Martian climatology. The relative sound‐speed changes and horizontal wind speed variations are synchronous, while temperature peaks ∼1.88 hr after these time series. The strong and continuous emergence of coherent phases in the autocorrelation codas suggests the presence of continuous infrasound on Mars.

 
more » « less
Award ID(s):
1847736
NSF-PAR ID:
10368352
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
8
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Seismic and infrasound multistation ambient‐noise interferometry has been widely used to infer ground and atmospheric properties, and single‐station and autocorrelation seismic interferometry has also shown potential for characterizing Earth structure at multiple scales. We extend autocorrelation seismic interferometry to ambient atmospheric infrasound recordings that contain persistent local noise from waterfalls and rivers. Across a range of geographic settings, we retrieve relative sound‐speed changes that exhibit clear diurnal oscillations consistent with temperature and wind variations. We estimate ambient air temperatures from variations in relative sound speeds. The frequency band from 1 to 2 Hz appears most suitable to retrieve weather parameters as nearby waterfalls and rivers may act as continuous and vigorous sources of infrasound that help achieve convergence of coherent phases in the autocorrelation codas. This frequency band is also appropriate for local sound‐speed variations because it has infrasound with wavelengths of ∼170–340 m, corresponding to a typical atmospheric boundary layer height. After applying array analysis to autocorrelation functions derived from a three‐element infrasound array, we find that autocorrelation codas are composed of waves reflected off nearby topographic features, such as caldera walls. Lastly, we demonstrate that autocorrelation infrasound interferometry offers the potential to study the atmosphere over at least several months and with a fine time resolution.

     
    more » « less
  2. SUMMARY

    The SEIS (Seismic Experiment for Interior Structure) seismometer deployed at the surface of Mars in the framework of the NASA-InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has been continuously recording the ground motion at Elysium Planitia for more than one martian year. In this work, we investigate the seasonal variation of the near-surface properties using both background vibrations and a particular class of high-frequency seismic events. We present measurements of relative velocity changes over one martian year and show that they can be modelled by a thermoelastic response of the Martian regolith. Several families of high-frequency seismic multiplets have been observed at various periods of the martian year. These events exhibit complex, repeatable waveforms with an emergent character and a coda that is likely composed of scattered waves. Taking advantage of these properties, we use coda wave interferometry (CWI) to measure relative traveltime changes as a function of the date of occurrence of the quakes. While in some families a stretching of the coda waveform is clearly observed, in other families we observe either no variation or a clear contraction of the waveform. These various behaviors correspond to different conditions of illumination at the InSight landing site, depending on the season. Measurements of velocity changes from the analysis of background vibrations above 5 Hz are consistent with the results from CWI. We identify a frequency band structure in the power spectral density (PSD) that can be tracked over hundreds of days. This band structure is the equivalent in the frequency domain of an autocorrelogram and can be efficiently used to measure relative traveltime changes as a function of frequency. We explain how the PSD analysis allows us to circumvent the contamination of the measurements by the Lander mode excitation which is inevitable in the time domain. The observed velocity changes can be adequately modelled by the thermoelastic response of the regolith to the time-dependent incident solar flux at the seasonal scale. In particular, the model captures the time delay between the surface temperature variations and the velocity changes in the subsurface. Our observations could serve as a basis for a joint inversion of the seismic and thermal properties in the first 20 m below InSight.

     
    more » « less
  3. Abstract

    Data assimilation is an increasingly popular technique in Mars atmospheric science, but its effect on the mean states of the underlying atmosphere models has not been thoroughly examined. The robustness of results to the choice of model and assimilation algorithm also warrants further study. We investigate these issues using two Mars general circulation models (MGCMs), with particular emphasis on zonal wind and temperature fields. When temperature retrievals from the Mars Global Surveyor Thermal Emission Spectrometer (TES) are assimilated into the U.K.‐Laboratoire de Météorologie Dynamique MGCM to create the Mars Analysis Correction Data Assimilation (MACDA) reanalysis, low‐level zonal jets in the winter northern hemisphere shift equatorward and weaken relative to a free‐running control simulation from the same MGCM. The Ensemble Mars Atmosphere Reanalysis System (EMARS) reanalysis, which is also based on TES temperature retrievals, also shows jet weakening (but less if any shifting) relative to a control simulation performed with the underlying Geophysical Fluid Dynamics Laboratory MGCM. Examining higher levels of the atmosphere, monthly mean three‐dimensional temperature and zonal wind fields are in generally better agreement between the two reanalyses than between the two control simulations. In conjunction with information about the MGCMs' physical parameterizations, intercomparisons between the various reanalyses and control simulations suggest that overall the EMARS control run is plausibly less biased (relative to the true state of the Martian atmosphere) than the MACDA control run. Implications for future observational studies are discussed.

     
    more » « less
  4. Abstract

    The extent to which the eddy statistics of the Martian atmosphere can be inferred from the mean state and highly simplified assumptions about diabatic and frictional processes is investigated using an idealized general circulation model (GCM) with Newtonian relaxation thermal forcing. An iterative technique, adapted from previous terrestrial studies, is used to generate radiative equilibrium temperatures such that the three-dimensional time-mean temperature fields of the idealized model match means computed from the Mars Analysis Correction Data Assimilation (MACDA). Focusing on a period of strong Northern Hemisphere eddy activity prior to winter solstice, it is found that the idealized model reproduces some key features of the spatial patterns of the MACDA eddy temperature variance and kinetic energy fields. The idealized model can also simulate aspects of MACDA’s seasonal cycle of spatial patterns of low-level eddy meridional wind and temperature variances. The most notable weakness of the model is its eddy amplitudes—both their absolute values and seasonal variations are quite unrealistic, for reasons unclear. The idealized model was also run with a mean flow based on output from the Geophysical Fluid Dynamics Laboratory (GFDL) full-physics Mars GCM. The idealized model captures the difference in mean flows between MACDA and the GFDL Mars GCM and reproduces a bias in the more complex model’s eddy zonal wavenumber distribution. This implies that the mean flow is an important influence on transient eddy wavenumbers and that improving the GFDL Mars GCM’s mean flow would make its eddy scales more realistic.

     
    more » « less
  5. Abstract

    The martian North Polar layered deposits (NPLD) are composed of alternating water‐ice and dust‐rich layers resulting from atmospheric deposition and are key to understanding Mars climate cycles. Within these deposits are spiral troughs whose migration affects deposition signals. To understand the relationship between NPLD stratigraphy and Martian climate, we must identify modern‐day drivers of NPLD ice migration. Prevailing theory posits migration driven by upstream‐migrating bed undulations bounded by hydraulic jumps, caused by katabatic winds flowing over trough walls with asymmetric cross‐sectional relief. This is supported by trough‐parallel clouds, whose formation has been attributed to hydraulic jumps. We present a cloud atlas across the Martian north pole using ∼13,800 THEMIS images spanning ∼18 Earth years. We find trough‐parallel clouds in ∼400 images, with regions nearer to the pole having higher cloud frequency. We compare spiral trough geometry to our cloud atlas and find regions with trough‐parallel clouds often correlate with metrics associated with modern‐day sublimation‐deposition cycles (i.e., relief and asymmetry), but not always. In some regions, troughs with morphologies conducive to cloud formation have no clouds. Overall, trough geometry varies greatly across the deposits, both within and between troughs, suggesting localized differences in deposition relative to migration, varying katabatic wind intensities, differing past climatic states influencing the troughs, varying trough initiation properties, or the possibility of additional mechanisms for trough initiation and migration (e.g., in situ trough erosion). Understanding what controls trough shape variability across the NPLD and how these controls change through time and space is key when interpreting Martian paleoclimate.

     
    more » « less