Simultaneous Estimation of Dust Content and Layer Thickness in Martian Polar Layered Deposits from SHARAD Radargrams Using Hidden Markov Model

Xiaofeng Liu,Wenzhe Fa
DOI: https://doi.org/10.1109/tgrs.2024.3370809
IF: 8.2
2024-01-01
IEEE Transactions on Geoscience and Remote Sensing
Abstract:The Martian polar layered deposits (PLD) are extensive, stratified sheets composed primarily of water ice with a small fraction of dust, which is believed to have formed due to climatic variations in the late Amazonian period. The thickness and dust content of layers in the PLD are crucial for understanding the paleoclimate they potentially record. The SHAllow RADar (SHARAD) onboard Mars Reconnaissance Orbiter (MRO) is currently the highest resolution radar sounding instrument for acquiring internal information of PLD. However, previous studies based on high-resolution optical and topographic data show that the thicknesses of the layers are usually smaller than the vertical resolution of SHARAD (similar to 8.5 m in water ice, epsilon(ice)' = 3.15). In this case, a bright radar echo is caused by a thin dusty ice layer and its reflectivity is determined by multiple reflections between the upper and lower boundaries of the layer. This leads to the nonunique solution for layer thickness and dust content from the observed reflectivity. In previously published studies, the dust content of a layer was determined by assuming an equal probability distribution of layer thickness, which is usually not realistic. To resolve this problem, we propose a novel algorithm for layer parameter inversion based on the hidden Markov model (HMM) that utilizes the spatial correlation of layer parameters. In this algorithm, the layer parameters and SHARAD measured reflection strength of a dusty layer correspond to the state and observation sequences, respectively. The spatial continuity and correlation of the layer parameter (i.e., layer thickness and dust content) are utilized to construct probability distributions of HMM, and the layer parameter with the highest probability is obtained as the inversion result. Compared to the previous method, our algorithm does not rely on prior knowledge of the layer thickness distribution. The algorithm is validated using both simulated reflectivity and SHARAD data. In the experiment with simulated reflectivity, the average inversion errors of the dust content and thickness obtained are 0.048 and 1.66 m, respectively. The error of the inverted dust content is reduced by 23.8% compared to the previous method. The layer thickness inverted from SHARAD data is comparable to that obtained by the high-resolution topographic data in the nearby outcrop region. In the future, our algorithm can be utilized to analyze radargrams collected by the orbital sounding radar onboard China's Tianwen-1 mission, providing a more detailed view of the internal structure of PLD.
What problem does this paper attempt to address?