火星轨道器次表层探测雷达数据处理技术与现状研究

火星是人类重要的地外天体探测目标之一,对火星表面进行的探测和研究表明,火星表面曾经存在液态水,水是生命存在的基础,因此,在次表层寻找不同形式的水是目前火星探测的重要科学目标之一.近17年来,欧洲火星快车(Mars Express)上搭载的火星次表层和电离层探测先进雷达(Mars Advanced Radar for S...     

Correction of the ionospheric distortion on the MARSIS surface sounding echoes

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft has made numerous measurements of the Martian surface and subsurface. However, all of these measurements are distorted by the ionosphere and must be compensated before any analysis. We have developed a technique to compensate for the ionospheric distortions. This technique provides a powerful tool to derive the total electron content (TEC) and other higher-order terms of the limited expansion of the plasma dispersion function that are related to overall shape of the electron column profile. The derived parameters are fitted by using a Chapman model to derive ionospheric parameters like n₀, electron density primary peak (maximum for solar zenith angle (SZA) equal 0), and the neutral height scale H.

Our estimated ionospheric parameters are in good agreement with Mars Global Surveyor (MGS) radio-occultation data. However, since MARSIS does not have the observation geometry limitations of the radio occultation measurements, our derived parameters extend over a large range of SZA for each MEX orbit.

The first results from our technique have been discussed by Safaeinili et al. [2007, Estimation of the total electron content of the Martian ionosphere using radar sounder surface echoes. Geophys. Res. Lett. 34, L23204, doi:10.1029/2007GL032154].     

Radar absorption due to a corotating interaction region encounter with Mars detected by MARSIS

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a subsurface and topside ionosphere radar sounder aboard the European Space Agency spacecraft Mars Express, in orbit at Mars since 25 December 2003, and in operation since 17 June 2005. The ionospheric sounding mode of MARSIS is capable of detecting the reflection of the sounding wave from the martian surface. This ability has been used in previous work to show that the surface reflection is absorbed and disappears during periods when high fluxes of energetic particles are incident on the ionosphere of Mars. These absorption events are believed to be the result of increased collisional damping of the sounding wave, caused by increased electron density below the spacecraft, in turn caused by impact ionization from the impinging particles. In this work we identify two absorption events that were isolated during periods when the surface reflection is consistently visible and when Mars is nearly at opposition. The visibility of the surface reflection is viewed in conjunction with particle and photon measurements taken at both Mars and Earth. Both absorption events are found to coincide with Earth passing through solar wind speed and ion flux signatures indicative of a corotating interaction region (CIR). The two events are separated by an interval of approximately 27 days, corresponding to one solar rotation. The first of the two events coincides with abruptly enhanced particle fluxes seen in situ at Mars. Simultaneous with the particle enhancement there are an abrupt decrease in the intensity of electron oscillations, typically seen by the Mars Express particle instrument ASPERA-3 between the magnetic pileup boundary and the martian bow shock, and a sharp drop in the solar wind pressure, seen in the proxy quantity based on MGS magnetometer observations. The decrease in oscillation intensity is therefore the probable effect of a relaxation of the martian bow shock. The second absorption event does not show a particle enhancement and complete ASPERA-3 data during that time are unavailable. Other absorption events are the apparent result of solar X-ray and XUV enhancements. We conclude that surface reflection absorption events are sometimes caused by enhanced ionospheric ionization from high energy particles accelerated by the shocks associated with a CIR. A full statistical analysis of CIRs in relation to observed absorption events in conjunction with a quantitative analysis of the deposition of ionization during space weather events is needed for a complete understanding of this phenomenon. If such analyses can be carried out, radar sensing of the martian ionosphere might be useful as a space weather probe.     

MARSIS surface reflectivity of the south residual cap of Mars

The south residual cap of Mars is commonly described as a thin and bright layer of CO₂-ice. The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a low-frequency radar on board Mars Express operating at the wavelength between 55 and 230 m in vacuum. The reflection of the radar wave on a stratified medium like the residual cap can generate interferences, causing weaker surface reflections compared to reflections from a pure water ice surface. In order to understand this anomalous low reflectivity, we propose a stratified medium model, which allows us to estimate both the thickness and the dielectric constant of the optically thin slab. First, we consider the residual cap as single unit and show that the decrease in the reflected echo strength is well explained by a mean thickness of 11 m and a mean dielectric constant of 2.2. This value of dielectric constant is close to the experimental value 2.12 for pure CO₂-ice. Second, we study the spatial variability of the radar surface reflectivity. We observe that the reflectivity is not homogeneous over the residual cap. This heterogeneity can be modeled either by variable thickness or variable dielectric constant. The surface reflectivity shows that two different units comprise the residual cap, one central unit with high reflectivity and surrounding, less reflective units.     

The influence of the ionosphere in subsurface Martian soil sounding experiments and a method of its correction

The problem of radar sounding of the Martian soil from an orbiting spacecraft is analyzed. The influence of the ionosphere on the results of remote sounding is revealed through the results of numerical modeling and actual subsurface soil sounding experiments carried out using the MARSIS radar sounder aboard the Mars-Express spacecraft. A method of correction of the dispersion distortion is proposed and tested through radiolocation of the polar cap on Mars.     

Reconstruction of nonmonotonic electron density profiles of the Martian topside ionosphere

One of the problems in reconstructing the real ionosphere from an ionogram is the occurrence of a ‘valley,’ where electron density decreases with altitude and make a non-monotonic profile. For the case of the Earth ionosphere, the ordinary and extraordinary ray data, accompanied with an empirical model, based on the observations, are necessary to obtain a mathematical solution for a ‘valley,’ such as the region between the E and F layers. MARSIS/MEX is a topside sounder designed to observe the ionosphere of Mars. Some ‘valley’ structures were found in the ionograms measured by MARSIS. The echoes of the extraordinary ray are not available owing to the absence of the strong magnetic field on Mars. Therefore, it is difficult to have a mathematical solution for the valleys in the Martian ionosphere. In this paper, a least square method with a simple model is presented to solve the ‘valley’ problem in the topside ionosphere of Mars. The electron density profiles with ‘valleys’ observed by the Radio Occultation experiment onboard MGS are used to rebuild the virtual depths at MARSIS frequencies. The reconstructed electron density profile by the least square method with a simple model from the rebuilt virtual depth curve is compared with the original electron density profile. It is proved that this method can reproduce small valleys in the profile of the Martian ionosphere well.     

Evolution of the Martian atmosphere and hydrosphere: Solar wind erosion studied by ASPERA-3 on Mars Express

The evolution of the Martian atmosphere and the potential existence of a past hydrosphere is a scientific issue of great interest in planetary research. Although the first missions to Mars had a focus on surface features and atmospheric properties, some of the missions (e.g., The Soviet Mars 2, 3 and 5) also carried instruments addressing the solar wind interaction with the Martian atmosphere and ionosphere and the potential existence of an intrinsic magnetic field on Mars. However, it took until 1989 before a spacecraft, Phobos-2, was able to carry out a more detailed investigation of the solar wind interaction with Mars. Phobos-2 gave valuable data on the Solar wind interaction with Mars during about 2 months of operations, leading to a better understanding of the solar wind impact on a weakly magnetized planet. However, Phobos-2 also raised a number of critical issues that has left science without adequate data since 1989.Investigations planned for Mars Express will cast new light on important aspects of the solar wind interaction with Mars. ASPERA-3 (Analyzer of Space Plasma and Energetic Atoms) on Mars Express will focus on the overall plasma outflow and monitor remotely the outflow and inflow of energetic neutral atoms produced by charge exchange processes. This report will discuss some of the unsolved issues about the solar wind interaction with Mars and how we plan to address these issues with Mars Express.     

Mars: Dual-polarization radar observations with extended coverage

Thirteen-centimeter-wavelength radar observations of Mars made in 1982 at Arecibo Observatory yield accurate measurements of the full backscatter spectrum in two orthogonal polarizations. The data, which were obtained for several widely separated subradar longitudes at 24°N latitude, provide the first global view of the distribution of small-scale surface roughness on Mars. The diffuse component of the echo exhibits strong spatial variations. Areas of maximum depolarization correlate well with volcanic regions (Tharsis and Elysium), while the heavily cratered upland terrain yields relatively low depolarization. Parts of Tharsis give near-complete depolarization (polaziation ratio $\text{μ}{}_{\mathrm{<mtext>c</mtext>}}\text{? 1}$ when viewed at oblique angles of incidence). Northern Martian plains regions (Tharsis, Elysium, and Amazonis) may comprise the most extensive area of severe decimeter-scale surface roughness in the inner Solar System. On the average, the northern Martian tropics yield higher diffuse radar cross sections (σ^D = 0.05–0.12) and a higher of degree disk-integrated depolarization (μ_c = 0.1–0.4) than is found for the Moon, Mercury, and Venus. Comparisons between the Moon and Mars using radar data, ground truth, and simple scattering models suggest that Mars possesses a relatively high average coverage by decimeter-scale rocks. Also discussed are several of the more interesting quasispecular scattering results, the most unsual of which were obtained over the Olympus Mons aureole region.     

A plasma flow velocity boundary at Mars from the disappearance of electron plasma oscillations

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft is capable of measuring ionospheric electron density by the use of two main methods: remote radar sounding and from the excitation of local plasma oscillations. The frequency of the locally excited electron plasma oscillations is used to measure the local electron density. However, plasma oscillations are not observed when the plasma flow velocity is higher than about 160 km/s, which occurs mainly in the solar wind and magnetosheath. As a consequence, in many passes, there is a sudden disappearance of the plasma oscillations as the spacecraft enters into the magnetosheath. This fact allows us to identify a flow velocity boundary on the dayside, between the ionosphere of Mars and the shocked solar wind. This paper summarizes the results of the measurement of 552 orbits mostly over a period from August 4, 2005 to August 17, 2007. The boundary points found using MARSIS have been verified by measurements from the Analyzer of Space Plasma and Energetic Atoms (ASPERA-3) Electron Spectrometer (ELS) instrument on Mars Express. The average position of the flow velocity boundary is compared to flow velocity simulations computed using hybrid model and other boundaries. The boundary altitude is slightly lower than the magnetic pile-up boundary determined using Phobos 2 and Mars Global Surveyor (MGS) crossings, but it is in good agreement with the induced magnetospheric boundary determined by ASPERA-3. Investigation of the effect of the crustal magnetic field revealed that the flow velocity boundary is raised at the locations with strong crustal magnetic fields.     

The 3–5 MHz global reflectivity map of Mars by MARSIS/Mars Express: Implications for the current inventory of subsurface H2O

We extracted the surface echo power from 2 years of MARSIS measurements. The retrieved values are calibrated to compensate for changes in the distance of the spacecraft to the surface and for the attenuation of the signal by the ionosphere. The results are used to build the first global map of surface echo power at 3–5 MHz. The surface echo power variations are primarily caused by kilometer-scale surface roughness. Then, we derive the values of dielectric constant of the shallow subsurface materials by normalizing the surface echo power map using a simulation of MARSIS signal from the MOLA topography. As a result, we obtain a map that characterizes the dielectric properties of the materials down to a few decameters below the surface. Dielectric properties vary with latitude, with high values in mid-latitudes belts (20–40°) and lower values at both equatorial and high latitudes. From the comparison of MARSIS reflectivity map to GRS observations, we conclude that the reflectivity decrease observed poleward of 50–60° corresponds to the onset of water-ice occurrence within the regolith. Assuming homogenous ground composition and texture at the scale of the MARSIS resolution cell, our inferred volume of ground water ice is of 10⁶ km³, equivalent to a polar cap. Low reflectivity areas are also observed in equatorial regions. From radar studies alone, equatorial low dielectric constant values could have different interpretations but the correlation with GRS hydrogen distribution rather points toward a water-related explanation.     

Determination of asteroid masses from their close encounters with Mars

An obstacle to the asteroid mass determination lies in the difficulty in isolating the gravitational perturbation exerted by a single asteroid on the planets, being strongly correlated and mixed up with those of many other asteroids. This hindrance may be avoided by the method of analysis presented here: an asteroid mass is estimated in correspondence with its close encounters with Mars where the acceleration it induces on the planet can be sufficiently disentangled from those generated by the remaining asteroid masses to calculate. We test this technique in the analysis of range observations to Mars Global Surveyor and Mars Express performed from 1999 to 2007. For this purpose, we adopt the dynamical model of the planetary ephemeris INPOP06 (Fienga et al., 2008), which includes the gravitational influences of the 300 most perturbing asteroids of the Martian orbit. We obtain the solutions of 10 asteroid masses that have the largest effects on this orbit over the period examined: they are generally in good agreement with determinations recently published.     

MARSIS radar sounder observations in the vicinity of Ma'adim Vallis, Mars

The MARSIS radar experiment aboard the ESA Mars Express satellite has recorded several unusual reflections in the Ma'adim Vallis region of Mars. These reflections display a wide variety of morphologies which are very different from those of reflections seen beneath the Polar Layered Deposits, Medusae Fossae Formation and Dorsa Argentea Formation. Their morphologies are sometimes very laterally extensive, parabolic or hyperbolic, and apparently deep, but they can also appear horizontal and shallow. Aided by a geological map of the Ma'adim Vallis region, the morphological, locational and temporal characteristics of the reflections have been studied individually in an attempt to constrain their origin. While some may be subsurface reflections based on their shallow morphologies and correlation with the Eridania Planitia basin network, all of the reflections are ambiguous to some degree, displaying characteristics that do not allow a definite subsurface- or possibly ionospheric-sourced mechanism to be proposed for their creation. Those with more exaggerated morphologies are regarded as being much more likely to result from ionospheric distortion rather than subsurface inhomogeneity.     

Development and testing of subsurface sampling devices for the Beagle 2 lander

For planetary landing missions, the capability to acquire samples of soil and rock is of high importance whenever complex analyses (e.g. isotopic studies) on these materials are to be carried out, or when samples are to be returned to Earth. Not only surface samples are of relevance, but in recent concepts at least for Mars landing missions also subsurface samples are required. Subsurface material on Mars is believed to have been protected from the inferred oxidants at the immediate surface while also being protected from the UV influx. Therefore, there is considerable hope that in subsurface soil samples on Mars, at least organic matter delivered by meteorites may be detected, and possibly also relics of earlier simple microbial life on the planet. Likewise, samples from the inside of Martian surface rocks promise to have been protected from weathering and for the same reason they are important for organic chemistry studies. In this paper, an overview is given of the development and science of two different subsurface sampling devices for the Beagle 2 lander of ESA's Mars Express mission, being a “Mole” subsurface soil sampler and a small rock coring and sampling mechanism. Besides their sampling function, both the Mole and the Corer/Grinder will provide data on physical properties of Martian soils and rock, respectively, through the way they interact with the sampled materials. Details of the Mole and Corer/Grinder design are presented, along with results of recent tests with prototypes in the laboratory on physically analogous sample materials.     

Sublimation of the Martian CO2 Seasonal South Polar Cap

The polar condensation/sublimation of CO₂, that involve about one fourth of the atmosphere mass, is the major Martian climatic cycle. Early observations in visible and thermal infrared have shown that the sublimation of the Seasonal South Polar Cap (SSPC) is not symmetric around the geographic South Pole.Here we use observations by OMEGA/Mars Express in the near-infrared to detect unambiguously the presence of CO₂ at the surface, and to estimate albedo. Second, we estimate the sublimation of CO₂ released in the atmosphere and show that there is a two-step process. From Ls=180° to 220°, the sublimation is nearly symmetric with a slight advantage for the cryptic region. After Ls=220° the anti-cryptic region sublimation is stronger. Those two phases are not balanced such that there is 22% ± 9 more mass the anti-cryptic region, arguing for more snow precipitation. We compare those results with the MOLA height measurements. Finally we discuss implications for the Martian atmosphere about general circulation and gas tracers, e.g. Ar.     

Dayside induced magnetic field in the ionosphere of Mars

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) onboard the Mars Express spacecraft has occasionally displayed surprising features. One such feature is the occurrence of a series of broadband, low-frequency echoes at equally spaced delay times after the sounder transmitter pulse. The interval between the echoes has been shown to be at the cyclotron period of electrons orbiting in the local magnetic field. The electrons are believed to be accelerated by the large voltages applied to the antenna by the sounder transmitter. Measurements of the period of these “electron cyclotron echoes” provide a simple technique for determining the magnitude of the magnetic field near the spacecraft. These measurements are particularly useful because Mars Express carries no magnetometer, so this is the only method available for measuring the magnetic field magnitude. Using this technique, results are presented showing the large scale structure of the draped field inside the magnetic pile-up boundary. The magnitude of the draped field is shown to vary from about 40 nT at a solar zenith angle of about 25°, to about 25 nT at a solar zenith angle of 90°. The results compare favorably with similar results from the Mars Global Surveyor spacecraft. A fitting technique is developed to derive the vector direction and magnitude of the draped magnetic field in cases where the spacecraft passes through regions with significant variation in the crustal field. The magnetic field directions are consistent with current knowledge of the draping geometry of the magnetic field around Mars.     

The microwave spectrum of Mars: An analysis

A weighted least squares fit to the best available data on the Martian microwave spectrum indicates that the brightness temperature decreases from long to short wavelengths, rather than increasing as expected from the solution of the one-dimensional equation of heat conduction. Reasonable assumptions on the ratio of electrical to thermal skin depths, on internal heat sources, on ferromagnetic materials, on radiative conduction, on compaction with depth, and on surface rpughness all fail in reproducing the deduced spectrum. A thin near-surface layer of a material with high dielectric constant and high millimeter wave absorption is needed. Since Mars exhibits marked surface overturn, a condensible material, namely liquid water, seems indicated. A layer of liquid water some tens of microns thick, on the average, localized in the top few millimeters of a Martian epilith with refractive index ? 1.6 fits the microwave spectrum, and the infrared and radar data as well. The origin of such a layer of liquid water and its possible exobiological significance are discussed. The distribution of water should be nonuniform over the disk and may help explain discordant microwave observations and the anomalous variation of infrared brigthness temperature with latitude. Further millimeter wave radio and radar studies of Mars are needed.     

Inversion of relief on Mars

Relief inversion has been invoked to explain a number of geomorphic features of the martian surface. Terrestrial relief inversion occurs when former depressions become elevated because their fill is more resistant to erosion than the surrounding terrain. It is a common product of long-term landscape evolution on Earth, especially in relatively stable intra-cratonic settings and flat, or near flat lying successions. The inverted relief will preserve relicts of former land surfaces and is therefore older than the surrounding terrain. Relief inversion can occur by a range of processes, including infill of depressions by intrinsically resistant material, selective secondary cementation via diagenesis and weathering, or surface armouring. We examine a number of possible cases of inverted relief on Mars that appear to have formed by these three processes. We suggest that the most likely cementing agents for surface induration are iron oxides, silica, and sulfates. Possible cementation mechanisms include fluid mixing during regional groundwater flow, cooling of hydrothermal or basinal fluids as they near the surface, and evaporation and sublimation of near surface water. Wind action appears the most common erosive process on Mars capable of the regional landscape lowering necessary for relief inversion to occur, unlike on Earth where both deflation and runoff are important. Preliminary crater densities of selected features show that the tops of the proposed inverted relief have considerably more craters than the surrounding plains, as is predicted by the inversion hypothesis. More accurate dating of inverted surfaces and the adjacent areas may provide a simple way of measuring the degree of erosion over time in at least some areas of Mars.     

Accurate Mars Express orbits to improve the determination of the mass and ephemeris of the Martian moons

The determination of the ephemeris of the Martian moons has benefited from observations of their plane-of-sky positions derived from images taken by cameras onboard spacecraft orbiting Mars. Images obtained by the Super Resolution Camera (SRC) onboard Mars Express (MEX) have been used to derive moon positions relative to Mars on the basis of a fit of a complete dynamical model of their motion around Mars. Since, these positions are computed from the relative position of the spacecraft when the images are taken, those positions need to be known as accurately as possible. An accurate MEX orbit is obtained by fitting two years of tracking data of the Mars Express Radio Science (MaRS) experiment onboard MEX. The average accuracy of the orbits has been estimated to be around 20–25 m. From these orbits, we have re-derived the positions of Phobos and Deimos at the epoch of the SRC observations and compared them with the positions derived by using the MEX orbits provided by the ESOC navigation team. After fit of the orbital model of Phobos and Deimos, the gain in precision in the Phobos position is roughly 30 m, corresponding to the estimated gain of accuracy of the MEX orbits. A new solution of the GM of the Martian moons has also been obtained from the accurate MEX orbits, which is consistent with previous solutions and, for Phobos, is more precise than the solution from the Mars Global Surveyor (MGS) and Mars Odyssey (ODY) tracking data. It will be further improved with data from MEX-Phobos closer encounters (at a distance less than 300 km). This study also demonstrates the advantage of combining observations of the moon positions from a spacecraft and from the Earth to assess the real accuracy of the spacecraft orbit. In turn, the natural satellite ephemerides can be improved and participate to a better knowledge of the origin and evolution of the Martian moons.     

Initial analysis of ion fluxes in the magnetotail of Mars based on simultaneous measurements on Mars Express and Maven

Simultaneous operation of two Mars satellites, equipped with instruments for the study of the plasma environment close to Mars, the European satellite Mars Express and American satellite MAVEN, allows one to investigate the influence of the interplanetary environment on the Martian magnetosphere and atmospheric losses, induced by the solar wind, for the first time, with a sufficient degree of confidence. In this paper, the data from measurements on the Mars Express satellite (MEX) of heavy ion losses are analyzed in comparison with the solar wind and magnetic field measurements on the MAVEN satellite. The main issue is the spatial structure of the escaping ion flux and the influence of the nonstationarity of the solar wind flux on the escape rate.     

Water distribution in Martian permafrost regions from joint analysis of HEND (Mars Odyssey) and MOLA (Mars Global Surveyor) data

We jointly analyze data from the High-Energy Neutron Detector (HEND) onboard the NASA Mars Odyssey spacecraft and data from the Mars Orbiter Laser Altimeter (MOLA) onboard the Mars Global Surveyor spacecraft. The former instrument measures the content of hydrogen (in the form of H₂O or OH) in the subsurface layer of soil and the latter instrument measures the surface albedo with respect to the flux of solar energy. We have checked the presence of a correlation between these two data sets in various Martian latitude bands. A significant correlation has been found between these data at latitudes poleward of 40° in the northern hemisphere and at latitudes 40°–60° in the southern hemisphere. This correlation is interpreted as evidence for the presence of stable water ice in these regions under a dry layer of soil whose thickness is determined by the condition for equilibrium between the condensation of water from the atmosphere and its sublimation when heated by solar radiation. For these regions, we have derived an empirical relation between the flux of absorbed solar radiation and the thickness of the top dry layer. It allows the burial depth of the water ice table to be predicted with a sub-kilometer resolution based on near-infrared albedo measurements. We have found no correlation in the southern hemisphere at latitudes >60°, although neutron data also suggest that water ice is present in this region under a layer of dry soil. We conclude that the thickness of the dry layer in this region does not correspond to the equilibrium condition between the water ice table and the atmosphere.