Radar Sounding Data Research Articles

Evidence of Ice‐Rich Layered Deposits in the Medusae Fossae Formation of Mars

AbstractSubsurface reflectors in radar sounder data from the Mars Advanced Radar for Subsurface and Ionospheric Sounding instrument aboard the Mars Express spacecraft indicate significant dielectric contrasts between layers in the Martian Medusae Fossae Formation (MFF). Large density changes that create dielectric contrasts are less likely in deposits of volcanic ash, eolian sediments, and dust, and compaction models show that homogeneous fine‐grained material cannot readily account for the inferred density and dielectric constant where the deposits are more than a kilometer thick. The presence of subsurface reflectors is consistent with a multi‐layer structure of an ice‐poor cap above an ice‐rich unit analogous to the Martian Polar Layered Deposits. The volume of an ice‐rich component across the entire MFF below a 300–600 m dry cover corresponds to a global equivalent layer of water of ∼1.5 to ∼2.7 m or ∼30%–50% of the total estimated in the North Polar cap.

Producing 3D radargrams from orbital radar sounding data at Mars: History, results, methods, lessons and plans

The Shallow Radar (SHARAD) instrument aboard Mars Reconnaissance Orbiter (MRO) has been conducting sounding operations of Mars' near surface for more than 16 years. Results of the SHARAD investigation have been well documented in the literature, with relatively high coverage density in the polar regions leading to some of the earliest scientific discoveries there. Data collection and ensuing surprises of SHARAD's primary and first extended science phases were the impetus for choosing the north polar region for producing a first 3D radar image. This early attempt was successful, producing the first 3D radar image of Mars' north polar layered deposits in Planum Boreum, and shortly thereafter a companion 3D radar image of the south polar layered deposits in Planum Australe. Subsequent work refined the 3D methodology used to produce these images, leading to a significantly higher quality 3D image of the former as well as the first 3D radar image of the debris-covered glaciers in the east-central portion of Deuteronilus Mensae (DM) in the mid-latitudes. Prior experience in terrestrial seismic data processing and analysis methods forms the basis of the 3D methodology used with Martian radar data, and has been indispensable in ongoing efforts to improve this methodology to further clarify the 3D images of the targets. The purposes of this article are to review 1) the SHARAD 3D work history and results, 2) the methodology developed and challenges encountered thus far in producing the SHARAD 3D images, 3) the broader impact on the processing and analysis of SHARAD data, and 4) current efforts and plans for producing follow-on SHARAD 3D images. Similarities and differences between orbital radar and seismic sounding and data processing are sprinkled throughout the article as reminders that the outcomes in this case are very much the product of cross-disciplinary knowledge and experience.

Ice Content of Mantling Materials in Deuteronilus Mensae, Mars

AbstractLayers of ice and dust (i.e., “mantle”) are found throughout the mid‐latitudes of Mars, where they drape and mute topography and record the recent climate history. Many of the youngest Late Amazonian “latitude‐dependent mantles” have been characterized, but the full range of mantle characteristics, including variations in ice content and age, has not been well documented. To advance our understanding of these ice‐dust units, we use SHAllow RADar (SHARAD) radar sounding data to constrain the subsurface physical properties of a widespread Middle Amazonian mantle unit within Deuteronilus Mensae. This region hosts a high density of glacial landforms, which allows assessment of the stratigraphic relationship between mantle deposits and glacial ice. SHARAD reflectors at the base of the mantle units are used to determine moderate dielectric constants (∼5.75) and high loss tangents (∼0.038). These radar properties imply relatively limited ice content and high‐loss mantling materials with possible enhanced roughness or subsurface scattering. Further, we interpret deeper reflectors that are continuous with glacier basal reflectors and that suggest materials with lower dielectric constants (∼3.2–4.5) and loss tangents (0.0071–0.0120) to be due to extensions of stagnant, possibly lithic‐rich, glacier ice that are buried by a mantle layer. This study demonstrates the high variability in the origins and evolution of mid‐latitude mantling units on Mars. While the most recent latitude‐dependent mantling deposits may still be ice‐rich, older deposits, such as those preserved in Deuteronilus Mensae and elsewhere in the mid‐latitudes of Mars, may be largely desiccated or contain a more heterogeneous and limited distribution of ice.

Comprehensive Radar Mapping of Malaspina Glacier (Sít' Tlein), Alaska—The World's Largest Piedmont Glacier—Reveals Potential for Instability

AbstractMalaspina Glacier, located on the coast of southern Alaska, is the world's largest piedmont glacier. A narrow ice‐cored foreland zone undergoing rapid thermokarst erosion separates the glacier from the relatively warm waters of the Gulf of Alaska. Glacier‐wide thinning rates for Malaspina are greater than 1 m/yr, and previous geophysical investigations indicated that bed elevation exceeds 300 m below sea level in some places. These observations together give rise to the question of glacial stability. To address this question, glacier evolution models are dependent upon detailed observations of Malaspina's subglacial topography. Here, we map 2,000 line‐km of the glacier's bed using airborne radar sounding data collected by NASA's Operation IceBridge. When compared to gridded radar measurements, we find that glaciological models overestimate Malaspina's volume by more than 30%. While we report a mean bed elevation 100 m greater than previous models, we find that Malaspina inhabits a broad basin largely grounded below sea level. Several subglacial channels dissect the glacier's bed: the most prominent of these channels extends at least 35 km up‐glacier from the terminus toward the throat of Seward Glacier. Provided continued foreland erosion, an ice‐ocean connection may promote rapid retreat along these overdeepened subglacial channels, with a global sea‐level rise potential of 1.4 mm.

A Dictionary-Based Integrated Simulation Approach to Model Large- and Small-Scale Coherent Surface Scattering Phenomena in Radar Sounder Data

A Dictionary-Based Integrated Simulation Approach to Model Large- and Small-Scale Coherent Surface Scattering Phenomena in Radar Sounder Data

Analysis of Lava Tubes’ Roughness and Radar Near-Nadir Regime Backscattering Properties

Lava tubes are terrestrial tunnel-like natural subsurface caves. Mounting evidence suggests their presence on the Moon and Mars. Planetary radar sounders are nadir-looking instruments operating in the high-frequency (HF)/very-HF (VHF) part of the spectrum with subsurface penetration capabilities. Recently, several studies either proposed future mission concepts for lava tubes’ detection or attempted to locate them on the Moon and Mars with the available radar-sounding data. Lava tubes are typically modeled as quasi-cylindrical structures but their actual geometry and their influence on the radar backscattering in near-nadir regime have never been investigated in the literature. These are crucial information for understanding the feasibility of detecting lava tubes by current and future planetary radar sounding systems. Accordingly, in this article: 1) we assess whether lava tubes are self-affine fractal surfaces at horizontal scales relevant to radio and microwave scattering and 2) we evaluate the effect of lava tube topography on the radar backscattering response in the near-nadir regime. Our experimental results, which are inferred from 3-D terrestrial laser scanning (TLS) data of planetary lava tube analogs, show that lava tubes: 1) are self-affine fractals at horizontal scales relevant to radar sounding and 2) they are electromagnetically rough surfaces, especially in the VHF band. We provide quantitative values on the lava tube fractal parameters and radar roughness losses along with a discussion on both: 1) the implication of our results on current radar sounding systems’ ability to detect lava tubes and 2) the planning of future missions devoted to lava tube detection and characterization.

Bistatic Radar Tomography of Shear Margins: Simulated Temperature and Basal Material Inversions

A significant portion of the sea level contributions of Antarctica and Greenland comes from ice streams, but the physical processes controlling ice stream width are poorly understood, especially when topographic controls are absent. Recent modeling studies have indicated that ice stream width may be controlled by elevated temperatures inside ice stream shear margins. While radio echo sounders can provide measurements of englacial water storage and subglacial conditions, existing radar-sounding techniques cannot measure temperature profiles at the scale required to test this hypothesis. We propose using a wide-angle radar survey and tomographic inversion to resolve temperature profiles, gradients, and anomalies at the scale required to study the thermophysical controls on shear margins. Recent work produced a bistatic radar system capable of obtaining the long offsets required for well-constrained inversions; however, shear-margin-specific temperature inversion techniques have not been developed for this system. In this article, we develop Newton’s method and alternating direction method of multipliers’ inversions for estimating temperature distribution and basal material across ice stream shear margins. We evaluate the performance of these inversion techniques on simulated bistatic radar-sounding data. Our results suggest that bistatic radar tomography experiments should be able to produce temperature maps on 50 m <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\times50$ </tex-math></inline-formula> m grids with <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$0.83~^{\circ} \text{C}~\pm ~0.084~^{\circ} \text{C}$ </tex-math></inline-formula> mean temperature error, <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$3.58~^{\circ} \text{C}~\pm ~0.20~^{\circ} \text{C}$ </tex-math></inline-formula> maximum temperature error, and an error in relative basal permittivity of <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$0.63~\pm ~0.08$ </tex-math></inline-formula> for a 4-km transect.

ИССЛЕДОВАНИЕ ИЗМЕНЧИВОСТИ РАДИОРЕФРАКЦИИ СВЧ-ИЗЛУЧЕНИЯ В ТРОПОСФЕРЕ АРКТИЧЕСКОГО РЕГИОНА РОССИЙСКОЙ ФЕДЕРАЦИИ

Using the dataset of aerological sounding carried out at five Arctic stations, the profiles of the main characteristics that determine radio refraction of the weather radar sounding beam were calculated and analyzed. These characteristics are the refractive index and its gradient. Seasonal parameters of the exponential model of tropospheric refraction were obtained over Naryan-Mar, Salekhard, Norilsk, Tiksi, and Chokurdakh stations. The frequency of occurrence of various refraction types with account of seasonal variability was investigated. The obtained data can be used for software processing of atmospheric radar sounding data over the Arctic region.

Double ridge formation over shallow water sills on Jupiter\u2019s moon Europa

Jupiter’s moon Europa is a prime candidate for extraterrestrial habitability in our solar system. The surface landforms of its ice shell express the subsurface structure, dynamics, and exchange governing this potential. Double ridges are the most common surface feature on Europa and occur across every sector of the moon, but their formation is poorly understood, with current hypotheses providing competing and incomplete mechanisms for the development of their distinct morphology. Here we present the discovery and analysis of a double ridge in Northwest Greenland with the same gravity-scaled geometry as those found on Europa. Using surface elevation and radar sounding data, we show that this double ridge was formed by successive refreezing, pressurization, and fracture of a shallow water sill within the ice sheet. If the same process is responsible for Europa’s double ridges, our results suggest that shallow liquid water is spatially and temporally ubiquitous across Europa’s ice shell.

Noise Removal and Feature Extraction in Airborne Radar Sounding Data of Ice Sheets

The airborne ice-penetrating radar (IPR) is an effective method used for ice sheet exploration and is widely applied for detecting the internal structures of ice sheets and for understanding the mechanism of ice flow and the characteristics of the bottom of ice sheets. However, because of the ambient influence and the limitations of the instruments, IPR data are frequently overlaid with noise and interference, which further impedes the extraction of layer features and the interpretation of the physical characteristics of the ice sheet. In this paper, we first applied conventional filtering methods to remove the feature noise and interference in IPR data. Furthermore, machine learning methods were introduced in IPR data processing for noise removal and feature extraction. Inspired by a comparison of the filtering methods and machine learning methods, we propose a fusion method combining both filtering methods and machine-learning-based methods to optimize the feature extraction in IPR data. Field data tests indicated that, under different conditions of IPR data, the application of different methods and strategies can improve the layer feature extraction.

A Novel Approach to the Detection and Imaging of Candidate Martian Subglacial Water Bodies by Radar Sounder Data

Research based on Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) data detected unusual radar bright basal reflections located at about 1.5 km depth in a Mars region denoted as Ultimi Scopuli. These reflections were interpreted as a signature of subglacial liquid water even though this interpretation is still being debated in the literature. In this article, we propose a novel approach to the detection and imaging of candidate subglacial liquid water from radar sounding data. The approach combines the radar echo power traces with a suitable digital elevation model to provide a bidimensional representation of the surface. Even if the imaging method reconstructs a representation of the surface, we prove that it can be used to identify subsurface bright reflections in icy regions. Imaging is feasible even if the basal interface is not directly included in the processed data for image generation. To support this experimental evidence, we show that a relationship exists between the value of the reflected echo power originating from the englacial layers and the basal-to-surface-echo-power ratio. The observed relationship holds on both Ultimi Scopuli radar sounding data acquired on Mars and Lake Vostok data acquired on Earth. Our results show that the 2-D imaging provides an alternative way for locating candidate subglacial liquid water bodies on Mars over large areas also where the basal interface is not directly measured. The proposed approach complements previous research for further evaluation of the actual presence of liquid water on Mars.

Altimetry Measurements From Planetary Radar Sounders and Application to SHARAD on Mars

Low-frequency radar sounders have the potential to generate altimetric profiles, but the feasibility of utilizing planetary radar sounding data as an alternative to laser altimetry has not been assessed using existing data to date. Therefore, we have developed, implemented, and evaluated an algorithm to process SHAllow RADar sounder (SHARAD) data on Mars (Experiment Data Records as available on the planetary data system) first into altimetry profiles and ultimately into digital terrain models (DTMs). The minimally processed data are pulse compressed, corrected for ionospheric distortion, zero-Doppler filtered, and incoherently summed. We then apply pulse re-tracking techniques adapted from terrestrial ocean altimetry to identify the surface return. From the surface return we compute the time-of-flight and hence the range from the spacecraft to the surface of the planet. The altimetry groundtracks are then co-registrated with Mars Orbiter Laser Altimeter (MOLA) to remove any biases resulting from residual ionospheric effects or timing issues. The altimetric profiles are finally used to create DTMs based on SHARAD data. While the SHARAD altimetry data have coarser inherent resolution than laser altimeters or imaging radars, we demonstrate that radar sounding data is still a viable source for satellite-based altimetry measurements. This is particularly important for future planetary missions not carrying laser altimeters but radar sounders, such as the upcoming Europa Clipper mission.

Empirical Models of Ion Density Distribution in the Dayside Martian Ionosphere

AbstractWith the aid of the recent Neutral Gas and Ion Mass Spectrometer (NGIMS) measurements made onboard the Mars Atmosphere and Volatile Evolution spacecraft, we construct empirical models for the vertical density distributions of 12 species in the dayside Martian ionosphere over the altitude range of 150–450 km, from abundant species such as , , and O+ to trace species such as and OCOH+. Two different sets of formulism are proposed to parameterize the ion density distribution characterized either by a distinctive layer structure or by a near exponential decay above 150 km. Both the solar control and magnetic control of ion distribution are taken into account in our modelings. In general, our empirical models are in good agreement with the NGIMS measurements. The full set of parameters constrained by data‐model comparison is tabulated for reference. This allows the easy computation of any individual ion density profile for a unique combination of solar zenith angle, solar extreme ultraviolet flux, and ambient magnetic elevation, which in turn facilitates comparisons with photochemical model results. The empirical models reported in this study are an extension of previous empirical models constructed for the electron distribution in the dayside Martian ionosphere using both radar sounding and radio occultation data.

An Unsupervised Fuzzy System for the Automatic Detection of Candidate Lava Tubes in Radar Sounder Data

Lava tubes are buried channels that transport thermally insulated lava. Nowadays, lava tubes on the Moon are believed to be empty and thus indicated as potential habitats for humankind. In recent years, several studies investigated possible lava tube locations, considering the gravity anomaly distribution and surficial volcanic features. This article proposes a novel and unsupervised method to map candidate buried empty lava tubes in radar sounder data (radargrams) and extract their physical properties. The approach relies on a model that describes the geometrical and electromagnetic (EM) properties of lava tubes in radargrams. According to this model, reflections in radargrams are automatically detected and analyzed with a fuzzy system to identify those associated with lava tube boundaries and reject the others. The fuzzy rules consider the EM and geometrical properties of lava tubes, and thus, their appearance in radargrams. The proposed method can address the complex task of identifying candidate lava tubes on a large number of radargrams in an automatic, fast, and objective way. The final decision on candidate lava tubes should be taken in postprocessing by expert planetologists. The proposed method is tested on both a real and a simulated data set of radargrams acquired on the Moon by the Lunar Radar Sounder (LRS). Identified candidate lava tubes are processed to extract geometrical parameters, such as the depth and the thickness of the crust (roof).

Measurement of sea wave height by radar sounding with polarizationmodulated signals

Methods for determining the height of sea waves, which are relevant for maritime navigation, have been investigated. The measurement of the height of sea waves by the method of radar sounding by polarization-modulated signals is considered. The energy and spectral characteristics of the radar signal reflected by the sea surface at low angles of incidence are described. An algorithm for calculating the average wave height based on the results of measuring the polarization anisotropy of the specific effective scattering surface of the sea area is considered. The relation of the average frequency of the amplitude fluctuations of the reflected radar signals with the sea waves is given. The spectral characteristics of radar signals reflected by the sea surface when irradiated with polarization-modulated signals are analyzed. The procedure for forming an effective statistical estimate of the average wave height based on the results of spectral measurements is considered. The errors of the generated estimates of the wave height associated with the measurement errors of the polarization anisotropy of the effective scattering surface and the average frequency of the envelope of the amplitude fluctuations caused by the Doppler effect are analyzed. The results of experimental measurements of the spectra of the amplitude fluctuations of the envelope of the reflected signals at different states of the sea surface are presented. The analytical dependence of the polarization anisotropy of the effective sea scattering surface on the irradiation angle relative to the wave propagation direction is given. A comparison of estimates of the height of sea waves obtained from radar sounding and meteorological observation data is carried out.

Radar Sounding of Subsurface Structure in Eastern Coprates and Capri Chasmata, Mars

AbstractWe surveyed the subsurface structure in eastern Coprates and Capri Chasmata in the equatorial region using high‐resolution visible images, digital terrain models, and radar sounding data. We identified subsurface reflectors in four areas of the chasmata. At the stratigraphic exposure on the chasmata walls, the corresponding depth of the reflector is ∼60 m. The bulk dielectric constants of the layers above the reflectors are calculated as 3.4–4.0, suggesting a rock‐air mixture with ∼39.3% and 46.1% porosity or a rock‐air‐ice mixture with &lt;21.2% water ice fraction. This high porosity corresponds to nonwelded and unconsolidated sediments emplaced by aeolian, fluvial, and volcanic activities. If water ice actually exists, further studies and discussions are required for the mechanism to maintain it within low latitudes.

Mars Surface Imaging by Exploiting Off-Nadir Radar Sounding Data

Radar sounder surface imaging is a rather unexplored approach to the analysis of planetary bodies. While a radar sounder is an instrument specifically designed for subsurface investigations, a particular set of power measurements (denoted as off-nadir surface echoes) can be exploited together with an external digital elevation model to produce images of the investigated surface at meters wavelength. The use of the off-nadir data may also reveal the presence of previously undetected subsurface features. In this article, we present a method for producing surface roughness images by high-frequency (HF) radar sounder data. The study of surface roughness in the HF band is particularly useful for both geologic studies and landing-zone reconnaissance as it is evaluated at meters to hundreds of meters horizontal scale. The proposed method combines off-nadir data of the Shallow Radar Sounder (SHARAD) with the Mars Orbiter Laser Altimeter (MOLA) digital elevation model. The produced roughness images at 20 MHz (15-m wavelength) of the Martian surface provide higher coverage and resolution of the surface roughness characterization at a 10-100-m horizontal scale than previous SHARAD work. By comparing the experimental roughness image with the one obtained by radar backscattering simulations, it is possible to identify subsurface features. In our experiments, we were able to produce a bidimensional image of a previously undetected large buried crater (10 km ×12 km) located in the Nili Fossae. This finding opens up new possibilities in exploiting radar sounding data for better detecting shallow subsurface features.

Automated detection and characterization of Antarctic basal units using radar sounding data: demonstration in Institute Ice Stream, West Antarctica

AbstractBasal units – visibly distinct englacial structures near the ice-bed interface – warrant investigation for a number of reasons. Many are of unknown composition and origin, characteristics that could provide substantial insight into subglacial processes and ice-sheet history. Their significance, moreover, is not limited to near-bed depths; these units appear to dramatically influence the flow of surrounding ice. In order to enable improved characterization of these features, we develop and apply an algorithm that allows for the automatic detection of basal units. We use a tunable layer-optimized SAR processor to distinguish these structures from the bed, isochronous englacial layers and the ice-sheet surface, presenting a conceptual framework for the use of radio-echo character in the identification of ice-sheet features. We also outline a method by which our processor could be used to place observational constraints on basal units’ configuration, composition and provenance.

Interferometric discrimination of cross-track bed clutter in ice-penetrating radar sounding data

AbstractThe interpretations of relevant interfaces (i.e. the surface and bed) in radar sounding datasets over glaciers and ice sheets are primary boundary conditions in a variety of climate studies and particularly subglacial water routing models. It is therefore necessary to ensure these interpretations are consistent and not affected by cross-track clutter. For the surface interface, interferometry and a family of methods relying on digital elevation models have been used to successfully discriminate cross-track surface clutter. Here we present how interferometry can be applied to the problem of basal clutter from cross-track bed topography. Our approach is based on a comparison of the differential phases of ambiguous reflectors that may represent bed clutter and the differential phase of a reflector in an adjacent area that appears unaffected by basal clutter. The reflector yielding the smallest interferometric phase difference relative to the unambiguous bed reflector is considered to represent its consistent continuation. We successfully demonstrate our approach using 60 MHz center frequency MARFA data collected over Devon Ice Cap in the Canadian Arctic. Finally, we investigate the effects of clutter-affected and interferometry-corrected bed interpretations on ice layer thickness estimates, basal hydraulic head gradients and the potential extent of inferred subglacial water bodies.

Five decades of radioglaciology

AbstractRadar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of theAnnals of Glaciologyon ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciologicalAnnals of Glaciologyissue of 2014, the context of their history and future prospects.