火星电离层无线电掩星探测仿真研究

对中俄联合火星星-星电离层掩星技术体制进行了分析和介绍,采用三维射线追踪方法对电离层掩星事件的电波观测值进行了模拟计算,并利用模拟的掩星观测数据进行了电子密度廓线反演,结果说明仿真算法可靠.利用仿真的方法,分别对掩星电波相位观测误差和卫星轨道误差等带来的反演误差进行了个例计算和分析,结果得到:5%周的相位测量误差对白天电离层掩星探测结果的影响可以忽略,而夜间电子密度测量的绝对误差小于4×108 m-3;卫星轨道误差对掩星的主要影响是导致电离层高度抬升或下降.结果表明,中俄联合火星电离层掩星探测技术体制先进,可望获得高精度的电子密度廓线;其技术体制也可以用于月球电离层环境的探测.     

电离层掩星反演的正则最大熵法

电离层掩星数据现已成为电离层观测数据的重要来源,对掩星数据的反演研究一直是掩星研究的热点.传统采用的改正TEC(1btal Electron Content)的Abel变换反演法为线性反演法,它会把测量误差带入反演结果中.为改善反演效果受测量误差的影响,引入两种非线性的反演方法一正则化和正则最大熵反演法.随后设计模拟试验,对3种方法进行验证、比较,得到正则最大熵反演法可以很好地减小测量误差的影响.     

Inversions of satellite to satellite electron content: Simulation studies with NeUoG-plas

The signals of Global Navigation Satellites have found a large number of uses in atmospheric and ionospheric research. Reception of the signals from a satellite in a Low Earth Orbit (LEO) leads regularly to occultation of the signals by the surface of the Earth. Before an occultation the signals traverse the ionosphere with rays with decreasing height of their perigees. Satellite electron content observed prior to ‘setting’ occultations or after `rising' occultations can be used as input data for inversion. The inversion procedure gives horizontally averaged height profiles of electron density.Assessment studies are needed to find out under which conditions the profiles from inversions are representative for ‘true’ electron density profiles above the Earth occultation point.A great number of such studies have been carried out using the ionosphere/plasmasphere model NeUOG-plas for forward and backward modelling. Different transmitter-receiver scenarios have been investigated.We describe the assessment procedure and report on results showing the most interesting cases and statistics.     

TEC Compensated Inversion Method of Ionospheric Occultation

For the inversion of ionospheric occultation data, the Abel inversion method of calibrated TEC is usually adopted, but the spherical asymmetry of the real ionosphere will bring errors to the inversion result of electron densities. This paper studies a TEC compensated inversion method. By combining with a background field, this method may eliminate the effect of the horizontal variability of electron density on the TEC, and improve the applicability of the assumption of spherical symmetry as well as the inversion accuracy. Applying this method to the inversion of simulative occultation events, the result demonstrates that compared with the traditional Abel transform inversion, the TEC compensated inversion method can reduce the inversion error. Comparing the various background fields obtained in different ways, it is shown that the higher the coincidence level between the background field and the real field, the better the inversion result.     

The solar wind interaction with Venus through the eyes of the Pioneer Venus Orbiter

Venus has no internal magnetic dynamo and thus its ionosphere and hot oxygen exosphere dominate the interaction with the solar wind. The solar wind at 0.72 AU has a dynamic pressure that ranges from 4.5 nPa (at solar max) to 6.6 nPa (at solar min), and its flow past the planet produces a shock of typical magnetosonic Mach number 5 at the subsolar point. At solar maximum the pressure in the ionospheric plasma is sufficient to hold off the solar wind at an altitude of 400 km above the surface at the subsolar point, and 1000 km above the terminators. The deflection of the solar wind occurs through the formation of a magnetic barrier on the inner edge of the magnetosheath, or shocked solar wind. Under typical solar wind conditions the time scale for diffusion of the magnetic field into the ionosphere is so long that the ionosphere remains field free and the barrier deflects almost all the incoming solar wind. Any neutral atoms of the hot oxygen exosphere that reach the altitude of the magnetosheath are accelerated by the electric field of the flowing magnetized plasma and swept along cycloidal paths in the antisolar direction. This pickup process, while important for the loss of the Venus atmosphere, plays a minor role in the deceleration and deflection of the solar wind. Like at magnetized planets, the Venus shock and magnetosheath generate hot electrons and ions that flow back along magnetic field lines into the solar wind to form a foreshock. A magnetic tail is created by the magnetic flux that is slowed in the interaction and becomes mass-loaded with thermal ions.The structure of the ionosphere is very much dependent on solar activity and the dynamic pressure of the solar wind. At solar maximum under typical solar wind conditions, the ionosphere is unmagnetized except for the presence of thin magnetic flux ropes. The ionospheric plasma flows freely to the nightside forming a well-developed night ionosphere. When the solar wind pressure dominates over the ionospheric pressure the ionosphere becomes completely magnetized, the flow to the nightside diminishes, and the night ionosphere weakens. Even at solar maximum the night ionosphere has a very irregular density structure. The electromagnetic environment of Venus has not been well surveyed. At ELF and VLF frequencies there is noise generated in the foreshock and shock. At low altitude in the night ionosphere noise, presumably generated by lightning, can be detected. This paper reviews the plasma environment at Venus and the physics of the solar wind interaction on the threshold of a new series of Venus exploration missions.     

Maximum Entropy Regularization Method of Ionospheric Occultation

With the increasing of applications of Global Positioning System (GPS), the research on the factors affecting the radio signals is becoming more and more important. One of the most significant effects on ionosphere monitoring is the ionospheric radio occultation (IRO), and the IRO data is now serving as one of the most important monitoring sources for the ionospheric measurements. Meanwhile the inversion of the occultation data is a popular topic. The traditional Abel inversion process through compensated total electron content (TEC) is a linear inversion method, thus it would transfer the measurement errors directly to the inversion results. In order to improve the occultation results, we introduce in this paper two nonlinear methods, namely the regularization method and the maximum entropy regularization method. Through designed simulative experiments, we verify and compare these three methods, and conclude that the maximum entropy regularization method can reduce significantly the influence of measurement errors.     

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].     

SPICAV on Venus Express: Three spectrometers to study the global structure and composition of the Venus atmosphere

Spectroscopy for the investigation of the characteristics of the atmosphere of Venus (SPICAV) is a suite of three spectrometers in the UV and IR range with a total mass of 13.9 kg flying on the Venus Express (VEX) orbiter, dedicated to the study of the atmosphere of Venus from ground level to the outermost hydrogen corona at more than 40,000 km. It is derived from the SPICAM instrument already flying on board Mars Express (MEX) with great success, with the addition of a new IR high-resolution spectrometer, solar occultation IR (SOIR), working in the solar occultation mode. The instrument consists of three spectrometers and a simple data processing unit providing the interface of these channels with the spacecraft.A UV spectrometer (118–320 nm, resolution 1.5 nm) is identical to the MEX version. It is dedicated to nadir viewing, limb viewing and vertical profiling by stellar and solar occultation. In nadir orientation, SPICAV UV will analyse the albedo spectrum (solar light scattered back from the clouds) to retrieve SO₂, and the distribution of the UV-blue absorber (of still unknown origin) on the dayside with implications for cloud structure and atmospheric dynamics. On the nightside, γ and δ bands of NO will be studied, as well as emissions produced by electron precipitations. In the stellar occultation mode the UV sensor will measure the vertical profiles of CO₂, temperature, SO₂, SO, clouds and aerosols. The density/temperature profiles obtained with SPICAV will constrain and aid in the development of dynamical atmospheric models, from cloud top (∼60 km) to 160 km in the atmosphere. This is essential for future missions that would rely on aerocapture and aerobraking. UV observations of the upper atmosphere will allow studies of the ionosphere through the emissions of CO, CO⁺, and CO₂⁺, and its direct interaction with the solar wind. It will study the H corona, with its two different scale heights, and it will allow a better understanding of escape mechanisms and estimates of their magnitude, crucial for insight into the long-term evolution of the atmosphere.The SPICAV VIS-IR sensor (0.7–1.7 μm, resolution 0.5–1.2 nm) employs a pioneering technology: an acousto-optical tunable filter (AOTF). On the nightside, it will study the thermal emission peeping through the clouds, complementing the observations of both VIRTIS and Planetary Fourier Spectrometer (PFS) on VEX. In solar occultation mode this channel will study the vertical structure of H₂O, CO₂, and aerosols.The SOIR spectrometer is a new solar occultation IR spectrometer in the range λ=2.2–4.3 μm, with a spectral resolution λ/Δλ>15,000, the highest on board VEX. This new concept includes a combination of an echelle grating and an AOTF crystal to sort out one order at a time. The main objective is to measure HDO and H₂O in solar occultation, in order to characterize the escape of D atoms from the upper atmosphere and give more insight about the evolution of water on Venus. It will also study isotopes of CO₂ and minor species, and provides a sensitive search for new species in the upper atmosphere of Venus. It will attempt to measure also the nightside emission, which would allow a sensitive measurement of HDO in the lower atmosphere, to be compared to the ratio in the upper atmosphere, and possibly discover new minor atmospheric constituents.     

GPS/LEO掩星方法中的电离层延迟量对太阳射电辐射流量的响应

对无线电信号在电离层中的传播路径进行了数值模拟;针对掩星观测中的五组实际卫星轨道数据(包括GPS和LEO卫星),给出了在太阳射电辐射流量(FLUX)分别为0,70,160和240等几种情况下的电离层延迟量对10.7cm波长的太阳辐射流量的响应结果;这分别对应着无太阳辐射、太阳射电辐射处于活动低谷、平静期和高峰期等几种情形。从中可以发现,掩星观测中电离层延迟量对10.7cm波长的太阳射电辐射流量的响应表现为如下特性:电离层延迟同参数FLUX的大小有明显的对应关系,即FLUX越大,则掩星观测中的电离层延迟越大;对于上升掩星情况而言,掩星观测中电离层延迟量起先逐渐增加,然后达到某一峰值,其后逐渐下降;在掩星观测的初期和中期,太阳射电辐射处于活动低谷、平静期和高峰期等几种情形之间的电离层延迟量差异都较为显著,而在掩星观测的后期,几种情形相互间电离层延迟量的差异都比较小。     

Non-thermal water loss of the early Mars: 3D multi-ion hybrid simulations

In this study we analyze the non-thermal loss rates of O⁺, O₂⁺ and CO₂⁺ ions over the last 4.5 billion years (Gyr) in the Martian history by using a 3D hybrid model. For this reason we derived the past solar wind conditions in detail. We take into account the intensified particle flux of the early Sun as well as an Martian atmosphere, which was exposed to a sun's extreme ultraviolet (EUV) radiation flux 4.5 Gyr ago that was 100 times stronger than today. Furthermore, we model the evolution of the interplanetary magnetic field by a Weber & Davis solar wind model. The ‘external’ influences of the Sun's radiation flux and solar wind flux lead to the formation of an ionospheric obstacle by photoionization, charge exchange and electron impact. For the early Martian conditions we could show that charge exchange was the dominant ionization mechanism. Several hybrid simulations for different stages in the evolution of the Martian atmosphere, at 1, 2, 5, 10, 30 and 100 EUV, were performed to analyze the non-thermal escape processes by ion pick-up, momentum transfer from the solar wind to the ionosphere and detached ionospheric plasma clouds. Our results show a non-linear evolution of the loss rates. Using mean solar wind parameters the simulations result in an oxygen loss equivalent to the depth of a global Martian ocean of about 2.6 m over the last 4.5 Gyr. The induced magnetic field strength could be increased up to about 2000 nT. A simulation run with high solar wind density results in an oxygen loss of a Martian ocean up to 205 m depth during 150 million years after the sun reached the zero age mean sequence (ZAMS).     

Comparisons of Cassini flybys of the Titan magnetospheric interaction with an MHD model: Evidence for organized behavior at high altitudes

Recent papers suggest the significant variability of conditions in Saturn’s magnetosphere at the orbit of Titan. Because of this variability, it was expected that models would generally have a difficult time regularly comparing to data from the Titan flybys. However, we find that in contrast to this expectation, it appears that there is underlying organization of the interaction features roughly above ～1800 km (1.7 Rt) altitude by the average external field due to Saturn’s dipole moment. In this study, we analyze Cassini’s plasma and magnetic field data collected at 9 Titan encounters during which the external field is close to the ideal southward direction and compare these observations to the results from a 2-fluid (1 ion, 1 electron) 7-species MHD model simulations obtained under noon SLT conditions. Our comparative analysis shows that under noon SLT conditions the Titan plasma interaction can be viewed in two layers: an outer layer between 6400 and 1800 km where interaction features observed in the magnetic field are in basic agreement with a purely southward external field interaction and an inner layer below 1800 km where the magnetic field measurements show strong variations and deviate from the model predictions. Thus the basic features inferred from the Voyager 1 flyby seem to be generally present above ～1800 km in spite of the ongoing external variations from SLT excursions, time variability and magnetospheric current systems as long as a significant southward external field component is present. At around ～1800 km kinetic effects (such as mass loading and heavy ion pickup) and below 1800 km ionospheric effects (such as drag of ionospheric plasma due to coupling with neutral winds and/or magnetic memory of Titan’s ionosphere) complicate what is observed.     

Trapped Ar isotopes in meteorite ALH 84001 indicate Mars did not have a thick ancient atmosphere

Water is not currently stable in liquid form on the martian surface due to the present mean atmospheric pressure of ～7 mbar and mean global temperature of ～220 K. However, geomorphic features and hydrated mineral assemblages suggest that Mars’ climate was once warmer and liquid water flowed on the surface. These observations may indicate a substantially more massive atmosphere in the past, but there have been few observational constraints on paleoatmospheric pressures. Here we show how the ⁴⁰Ar/³⁶Ar ratios of trapped gases within martian meteorite ALH 84001 constrain paleoatmospheric pressure on Mars during the Noachian era [～4.56–3.8 billion years (Ga)]. Our model indicates that atmospheric pressures did not exceed ～1.5 bar during the first 400 million years (Ma) of the Noachian era, and were <400 mbar by 4.16 Ga. Such pressures of CO₂ are only sufficient to stabilize liquid water on Mars’ surface at low latitudes during seasonally warm periods. Other greenhouse gases like SO₂ and water vapor may have played an important role in intermittently stabilizing liquid water at higher latitudes following major volcanic eruptions or impact events.     

Mars mesospheric zonal wind around northern spring equinox from infrared heterodyne observations of CO2

We present direct observations of Mars zonal wind velocities around northern spring equinox (L_S = 336°, L_S = 355°, L_S = 42°) during martian year 27 and 29. Data was acquired by means of infrared heterodyne spectroscopy of CO₂ features at 959.3917 cm^?1 (10.4232 μm) and 957.8005 cm^?1 (10.4405 μm) using the Cologne Tuneable Heterodyne Infrared Spectrometer (THIS) at the McMath–Pierce telescope of the National Solar Observatory on Kitt Peak in Arizona and the NASA Infrared Telescope Facility on Mauna Kea, Hawaii between 2005 and 2008. Winds were measured on the dayside of Mars with an unprecedented spatial resolution allowing sampling of up to nine independent latitudes over the martian disk. Retrieved wind velocities depend strongly on latitude and season with values ranging from 180 m/s prograde to ?94 m/s retrograde. A comparison of the observational results to predicted values from the Mars Climate Database yield a reasonable agreement between modeling and observation.     

Is the primordial crust of Mars magnetized?

A meteorite impact capable of creating a 200 km diameter crater can demagnetize the entire crust beneath, and produce an appreciable magnetic anomaly at satellite altitudes of ～400 km in case the pre-existing crust is magnetized. In this study we examine the magnetic field over all of the craters and impact-related Quasi-Circular Depressions (QCDs) with diameters larger than 200 km that are located on the highlands of Mars, excluding the Tharsis bulge, in order to estimate the mean magnetization of the highland crust. Using the surface topography and the gravity of Mars we first identify those QCDs that are likely produced by impacts. The magnetic map of a given crater or impact-related QCD is derived using the Mars Global Surveyor high-altitude nighttime radial magnetic data. Two extended ancient areas are identified on the highlands, the South Province and the Tempe Terra, which have large number of craters and impact-related QCDs but none of them has an appreciable magnetic signature. The primordial crust of these areas is not magnetized, or is very weakly magnetized at most. We examine some plausible scenarios to explain the weak magnetization of these areas, and conclude that no strong dynamo existed in the first ～100 Myr of Mars’ history when the newly formed primordial crust was cooling below the magnetic blocking temperatures of its minerals.     

Local plasma processes and enhanced electron densities in the lower ionosphere in magnetic cusp regions on Mars

Both the MARSIS ionospheric sounder and the charged particle instrument package ASPERA-3 are experiments on board the Mars Express spacecraft. Joint observations have shown that events of intense ionospheric electron density enhancements occur in the lower ionosphere of magnetic cusp regions, and that these enhancements are not associated with precipitation of charged particles above a few hundred electron volts (<300 eV). To account for the enhancement by particle precipitation, electron fluxes are required with mean energy between 1 and 10 keV. No ionizing radiation, neither energetic particles nor X-rays, could be identified, which could produce the observed density enhancement only in the spatially limited cusp regions. Actually, no increase in ionizing radiation, localized or not, was observed during these events. It is argued that the process causing the increase in density is controlled mainly by convection of ionosphere plasma driven by the interaction between the solar wind and crustal magnetic field lines leading to excitation of two-stream plasma waves in the cusp ionosphere. The result is to heat the plasma, reduce the electron–ion recombination coefficient and thereby increase the equilibrium electron density.     

Small-scale variability in Saturn’s lower ionosphere

We perform and present a wavelet analysis on all 31 Cassini electron density profiles published to date (Nagy, A.F. et al. [2006]. J. Geophys. Res. 111 (A6), CiteID A06310; Kliore, A.J. et al. [2009]. J. Geophys. Res. 114 (A4), CiteID A04315). We detect several discrete scales of variability present in the observations. Small-scale variability (S < 700 km) is observed in almost all data sets at different latitudes, both at dawn and dusk conditions. The most typical scale of variability is 300 km with scales between 200 km and 450 km being commonly present in the vast majority of the profiles. A low latitude dawn/dusk asymmetry is noted in the prevalent scales with the spectrum peaking sharply at the 300 km scale at dusk conditions and being broader at dawn conditions. Compared to dawn conditions the dusk ionosphere also shows more significant variability at the 100 km scale. The 300 km vertical scale is also present in the few available profiles from the northern hemisphere. Early observations from 2005 show a dominant scale at 350 km whereas later in 2007–2008 the spectrum shifts to the shorter scales with the most prominent scale being 300 km. The performed wavelet analysis and the obtained results are independent of assumptions about the nature of the layers and do not require a definition for a “background” electron density profile.In the second part of the paper we present a gravity wave propagation/dissipation model for Saturn’s upper atmosphere and compare the wave properties to the characteristics of the observed electron density variability at different scales. The general features observed in the data are consistent with gravity waves being present in the lower ionosphere and causing layering of the ions and the electrons. The wave-driving mechanism provides a simultaneous explanation for several of the properties of the observed variability: (i) lack of variability in the electron density above the predicted region of wave dissipation; (ii) in most cases the peak amplitude of variability occurs within the altitude range for dissipation of gravity waves or below; (iii) shorter scales have smaller amplitudes than the longer scales; (iv) shorter scales are present at lower altitudes whereas longer scales persist to higher altitudes; and (v) several layers often form a system of equally spaced maxima and minima that can be traced over a large altitude range.     

A New look at the ionosphere of Jupiter in light of the UVS occultation results

Simple photochemical models cannot reconcile Jupiter's ionospheric electron density profiles with the observed neutral atmosphere. The location of the peak electron density predicted when the neutral atmosphere determined by theVoyager Ultraviolet Spectrometer is combined with simple models falls about 1000km lower than the peak determined by radio occultation. The locations and magnitudes of the peaks in electron density can be accounted for by including the effects of vertical transport of ions in the ionospheric models. This vertical transport may be induced by meridional winds in the neutral atmosphere or external electric fields. It is probable that precipitating particles and an altitude-variable H₂ vibrational temperature play important roles in determining the character of the iono?phere. In view of the complex relationship between the ionosphere and neutral atmosphere, an attempt to infer one from the other cannot succeed. However, combining independent information on the two leads to new insights into the coupling of the neutral atmosphere, the ionosphere and the magnetosphere.     

Dynamical evolution of the Hungaria asteroids

The Hungarias are a stable asteroid group orbiting between Mars and the main asteroid belt, with high inclinations (16–30°), low eccentricities (e < 0.18), and a narrow range of semi-major axes (1.78–2.06 AU). In order to explore the significance of thermally-induced Yarkovsky drift on the population, we conducted three orbital simulations of a 1000-particle grid in Hungaria a–e–i space. The three simulations included asteroid radii of 0.2, 1.0, and 5.0 km, respectively, with run times of 200 Myr. The results show that mean motion resonances—martian ones in particular—play a significant role in the destabilization of asteroids in the region. We conclude that either the initial Hungaria population was enormous, or, more likely, Hungarias must be replenished through collisional or dynamical means. To test the latter possibility, we conducted three more simulations of the same radii, this time in nearby Mars-crossing space. We find that certain Mars crossers can be trapped in martian resonances, and by a combination of chaotic diffusion and the Yarkovsky effect, can be stabilized by them. Therefore, some Hungarias (around 5% of non-family members with absolute magnitudes H < 15.5 and 10% for H < 17) may represent previously transient Mars crossers that have been adopted in this manner.     

Low-magnetic crust underlying South Province of Mars

The lack of distinct magnetic signatures observed by Mars Global Surveyor (MGS) over the impact craters and impact-related Quasi-Circular-Depressions (QCDs) with diameters greater than 200 km located on South Province, south of 30S and from almost the west of Hellas to Argyre basins, implies a weakly magnetized crust. Using MOLA topography and the recent JPL gravity model of Mars we determine the structure of the crust beneath the craters and impact-related QCDs, and show that the impacts that have created these features were capable of strongly disturbing the crust directly beneath. On the basis of theoretical magnetic anomaly modeling and shock demagnetization models, we demonstrate that the impacts are capable of demagnetizing the entire crust beneath and creating distinct magnetic anomalies at the satellite altitude of 400 km in case the crust was appreciably magnetized prior to the impacts. We derive the magnetic anomalies of these features using the radial component of the high-altitude nighttime MGS data. An upper limit of <2 × 10⁴ A for the bulk magnetization of the crust beneath South Province is estimated, which is about 30 times less than that underlying Terra Cimmeria and Terra Sirenum. Similar weak bulk magnetization is obtained for part of the crust surrounding Hellas, Isidis, and Argyre basins.     

Evidence for topographic effects on the Martian ionosphere

The Radio occultation experiment on board Mariner 9 has been used to demonstrate that the altitude of the main electron density peak in the Martian ionosphere is closely related to the height of Mars’ surface at the occultation point. This is direct evidence for topographic effects on the Martian ionosphere. Modeling indicates that topographic-induced modulations of the neutral density in the upper atmosphere can account for the observed ionospheric effects. The neutral density modulation is likely to be caused by nonmigrating tides in the Martian thermosphere.