An overview of radar soundings of the martian ionosphere from the Mars Express spacecraft

The Mars Express spacecraft carries a low-frequency radar called MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) that is designed to study the subsurface and ionosphere of Mars. In this paper, we give an overview of the ionospheric sounding results after approximately one year of operation in orbit around Mars. Several types of ionospheric echoes are commonly observed. These include vertical echoes caused by specular reflection from the horizontally stratified ionosphere; echoes from a second layer in the topside ionosphere, possibly associated with O⁺ ions; oblique echoes from upward bulges in the ionosphere; and a variety of other echoes that are poorly understood. The vertical echoes provide electron density profiles that are in reasonable agreement with the Chapman photo-equilibrium model of planetary ionospheres. On the dayside of Mars the maximum electron density is approximately 2 × 10⁵ cm⁻³. On the nightside the echoes are often very diffuse and highly irregular, with maximum electron densities less than 10⁴ cm⁻³. Surface reflections are sometimes observed in the same frequency range as the diffuse echoes, suggesting that small isolated holes exist in the nightside ionosphere, possibly similar to those that occur on the nightside of Venus. The oblique echoes arise from upward bulges in the ionosphere in regions where the crustal magnetic field of Mars is strong and nearly vertical. The bulges tend to be elongated in the horizontal direction and located in regions between oppositely directed arch-like structures in the crustal magnetic field. The nearly vertical magnetic field lines in the region between the arches are thought to connect into the solar wind, thereby allowing solar wind electrons to heat the lower levels of the ionosphere, with an attendant increase in the scale height and electron density.     

Temporal and spatial variations observed in the ionospheric composition of Venus: Implications for empirical modelling

In situ measurements of the thermal ion composition of the ionosphere of Venus have been obtained for a period of two Venus years from the Bennett rf ion mass spectrometer on the Pioneer Venus Orbiter. Ion measurements within an altitude interval of 160 to 300 kilometers, corresponding to an overall latitude interval of about ?4° to 34°N, are assembled from the interval December 1978 to March 1980. This time interval corresponds to two revolutions of Venus about the Sun, designated as two “diurnal cycles”. The distributions of several ion species in this data base have been sorted to identify temporal and spatial variations, and to determine the feasibility of an analytical representation of the experimental results. The first results from the sorting of several prominent ions including O⁺, O₂⁺, and H⁺ and several minor ions including CO₂⁺, C⁺, and H₂⁺ reveal significant diurnal variations, with superimposed modulation associated with solar activity and solar wind variations. The diurnal variation consists of strong day to night contrast in the ion concentrations, with differences of one to two orders of magnitude, depending upon ion mass and altitude. The concentrations of O₂⁺, O⁺, CO₂⁺ and C⁺ peak throughout the dayside decreasing sharply at the terminators to nightside levels, lower by one to two orders of magnitude relative to the dayside. The diurnal variations of the light ions H⁺ and H₂⁺ peak during the night, exhibiting asymmetric nightside bulges favoring the pre-dawn sector, near 0400 solar hour angle. Superimposed upon the diurnal distributions are modulation signatures which correlate well with modulation in the F_10.7 index, indicating a strong influence of solar variability on the ion production and distribution. The influence of solar wind perturbations upon the ion distributions are also indicated, by a significant increase in the scatter of the observations with increasing altitude as higher altitudes, approaching 300 kilometers, are sampled. Together, these temporal and spatial variations make the task of modelling the ionosphere of Venus both very interesting and challenging.     

Solar cycle effect on temporal and spatial variation of topside ion density measured by SROSS C2 and ROCSAT 1 over the Indian longitude sector

We investigated the diurnal, seasonal and latitudinal variations of ion density Ni over the Indian low and equatorial topside ionosphere within 17.5°S to 17.5°N magnetic latitudes by combining the data from SROSS C2 and ROCSAT 1 for the 9 year period from 1995 to 2003 during solar cycle 23. The diurnal maximum density is found in the local noon or in the afternoon hours and the minimum occurs in the pre sunrise hours. The density is higher during the equinoxes as compared to that in the June and December solstice. The local time spread of the daytime maximum ion density increases with increase in solar activity. A north south asymmetry with higher ion density over northern hemisphere in the June solstice and over southern hemisphere in December solstice has been observed in moderate and high solar activity years. The crest to crest distance increases with increase in solar flux. Ion density bears a nonlinear relationship with F_10.7 cm solar flux and EUV flux in general. The density increases linearly with solar flux up to ∼150 sfu (1 sfu = 10⁻²²Wm⁻²Hz⁻¹) and EUV flux up to ∼50 units (10⁹ photons cm⁻² s⁻¹). But beyond this the density saturates. Inverse saturation and linear relationship have been observed in some season or latitude also. Inter-comparison of the three solar activity indices F_10.7 cm flux, EUV flux and F_10.7P (= (F_10.7 + F_10.7A)/2, where F_10.7A is the 81 day running average value of F_10.7) shows that the ion density correlates better with F_10.7P and F_10.7 cm fluxes. The annual average daytime total ion density from 1995 to 2003 follows a hysteresis loop as the solar cycle reverses. The ion density at 500 km over the Indian longitude sector as obtained by the international reference ionosphere is in general lower than the measured densities during moderate and high solar activity years. In low solar activity years the model densities are equal or higher than measured densities. The IRI EIA peaks are symmetric (±10°) in equinox while densities are higher at 10°N in June solstice and at 10°S in the December solstice. The model density follows F_10.7 linearly up to about F_10.7 > ∼150 sfu and then saturates.     

Magnetosphere response to the 2005 and 2006 extreme solar events as observed by the Cluster and Double Star spacecraft

The four identical Cluster spacecraft, launched in 2000, orbit the Earth in a tetrahedral configuration and on a highly eccentric polar orbit (4–19.6 R_E). This allows the crossing of critical layers that develop as a result of the interaction between the solar wind and the Earth’s magnetosphere. Since 2004 the Chinese Double Star TC-1 and TC-2 spacecraft, whose payload comprise also backup models of instruments developed by European scientists for Cluster, provided two additional points of measurement, on a larger scale: the Cluster and Double Star orbits are such that the spacecraft are almost in the same meridian, allowing conjugate studies. The Cluster and Double Star observations during the 2005 and 2006 extreme solar events are presented, showing uncommon plasma parameters values in the near-Earth solar wind and in the magnetosheath. These include solar wind velocities up to ∼900 km s⁻¹ during an ICME shock arrival, accompanied by a sudden increase in the density by a factor of ∼5 and followed by an enrichment in He⁺⁺ in the secondary front of the ICME. In the magnetosheath ion density values as high as 130 cm⁻³ were observed, and the plasma flow velocity there reached values even higher than the typical solar wind velocity. These resulted in unusual dayside magnetosphere compression, detection of penetrating high-energy particles in the magnetotail, and ring current development following several successive injections of energetic particles in the inner magnetosphere, which “washed out” the previously formed nose-like ion structures.     

In situ results on the variation of neutral atmospheric hydrogen at venus

The concentrations of neutral hydrogen within the atmosphere of Venus are investigated for the period 1979–1980. During this period, the planet made nearly three orbits about the Sun, so that nearly three complete diurnal cycles were observed from the Pioneer Venus Orbiter (PVO). Values of n(H) are derived from in-situ ion and neutral composition measurements from the Orbiter Ion Mass Spectrometer (OIMS) and the Orbiter Neutral Mass Spectrometer (ONMS) using a charge exchange relationship involving O⁺, H⁺, O and CO₂. The dawn bulge in the diurnal distribution of n(H), reported from the first diurnal cycle by Brinton et al., is found to persist with n(H) peaking at levels near 2 - 5 × 10⁷/cm³ at altitudes below 165 km. At peak levels, the bulge exhibits a concentration ratio up to 400/1 relative to dayside values. Large day to day variations of up to a factor of five in n(H) are frequently encountered, and are attributed to perturbations induced by the solar wind interaction. These short term variations, plus a suggestion of some local time variation in the bulk location, make precise assessment of interannual variations in the n(H) difficult. Between the first diurnal cycle in early 1979 and the third in mid 1980, the decline in solar euv flux was of the order of 10% or less. Allowing for uncertainties due to short term variations, no clear evidence is found for an interannual variation in the hydrogen concentrations.     

Model calculations of the dayside ionosphere of Venus

Model calculations of the dayside ionosphere of Venus are presented. The coupled continuity and momentum equations were solved for O₂⁺, O⁺, CO₂⁺, C⁺, N⁺, He⁺, and H⁺ density distributions, which are compared with measurements from the Pioneer Venus ion mass spectrometer. The agreement between the model results and the measurements is good for some species, such as O⁺, and rather poor for others, such as N⁺, indicating that our understanding of the dayside ion composition of Venus is incomplete. The coupled heat conduction equations for ions and electrons were solved and the calculated temperatures compared with Pioneer Venus measurements. It is shown that fluctuations in the magnetic field have a significant effect on the energy balance of the ionosphere.     

Quantifying the relationship between the measurement precision and specifications of a UV/visible sensor on a geostationary satellite

To investigate the feasibility of new satellite observations, including air quality (AQ) observations from geostationary (GEO) orbit, it is essential to link the measurement precision (ε) with sensor specifications in advance. The present study attempts to formulate the linkage between ε and specifications of a UV/visible sensor (signal-to-noise ratio (SNR), full width at half maximum (FWHM) of the slit function, and sampling ratio (SR)) on a GEO satellite. A sophisticated radiative transfer model (JACOSPAR) is used to calculate synthetic radiance spectra that would be measured by a UV/visible sensor observing the atmosphere over Tokyo (35.7°N, 139.7°E) from GEO orbit at 120°E longitude. The spectra, modified according to given sensor specifications, are analyzed by the differential optical absorption spectroscopy technique to estimate the ε for slant column densities of O₃ and NO₂. We find clear relationships: for example, the ε of the O₃ slant column density (molecules cm⁻²) and SNR at 330 nm are linked by the equation log(ε) = −1.06 · log(SNR) + 20.71 in the UV region, and the ε of the NO₂ slant column density and SNR at 450 nm are linked by log(ε) = −0.98 · log(SNR) + 18.00, at a FWHM = 0.6 nm (for the Gaussian slit function) and SR = 4. The relationships are mostly independent of other specifications (e.g., horizontal and temporal resolutions), as they affect ε primarily through SNR, providing constraints in determining the optimal SNR (and alternatively FWHM and SR) for similar UV/visible sensors dedicated for AQ studies.     

Hybrid modelling the Pioneer Venus Orbiter magnetic field observations

This study presents comparisons between the Pioneer Venus Orbiter (PVO) magnetometer (OMAG) observations and the HYB-Venus hybrid simulation code. The comparisons are made near periapsides of four PVO orbits using the full resolution PVO/OMAG data. Also, the statistics of the solar wind and interplanetary magnetic field (IMF) conditions at Venus are studied using the PVO interplanetary dataset. The statistics include the histograms and the probability density maps of the selected upstream parameters. The confidence intervals derived from the upstream statistics demonstrate the nominal simulation input parameter space. Moreover, the probability density maps give the dependencies between the upstream parameters. The comparisons between the simulation code and the data along the spacecraft trajectory show that the basic, large scale, trends seen in the magnetic field can be understood by the current simulation runs. The discrepancies between the simulation and the data were found to arise at low altitudes close to the planetary ionosphere in the region which cannot be resolved in detail by the grid size of the runs.     

Corotating interplanetary streams and associated ionospheric disturbances at Venus and Earth

During the summer of 1979, solar coronal structure was such that a sequence of recurrent regions produced a corresponding sequence of corotating solar wind streams, with pronounced downstream signatures. One of these stream events passed Earth on July 3, and was observed later at Venus late on July 11th, with similar characteristics. Corresponding in-situ measurements at Earth from the Atmospheric Explorer-E satellite and at Venus from the Pioneer Venus Orbiter are examined for evidence of comparable perturbations of the planetary ionospheres. The passage of the stream shock front is found to be associated with pronounced fluctuations in n(0⁺) which appear as pronounced local depletion of ion concentrations in both ionospheres. The ionosphere disturbances appear to be closely associated with large variations in the solar wind momentum flux. The implied local ionospheric depletions observed at each planet are interpreted to be the consequence of plasma redistribution, rather than actual depletions of plasma.     

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.     

Relative contributions of terrestrial and solar wind ions in the plasma sheet

A major uncertainty concerning the origins of plasma sheet ions is due to the fact that terrestrial H⁺ can have similar fluxes and energies as H⁺ from the solar wind. The situation is especially ambiguous during magnetically quiet conditions (AE < 60γ) when H⁺ typically contributes more than 90% of the plasma sheet ion population. In this study we examine that problem using a large data set obtained by the ISEE-1 Plasma Composition Experiment. The data suggest that one component of the H⁺ increases in energy with increasing activity, roughly in proportion to $\text{1}\text{4}$ the energy of the He⁺⁺, whereas the other H⁺ component has about the same energy at all activity levels, as do the O⁺ and the He⁺. If we can assume that the H⁺ of solar wind origin on the average has about the same energy-per-nucleon as the He⁺⁺, which is presumably almost entirely from the solar wind, then the data imply that as much as 20–30% of the H⁺ can be of terrestrial origin even during quiet conditions.     

A new approach to solar wind monitoring

The monitoring of solar wind parameters is a key problem of the space weather program. We are presenting a new solution of plasma parameter determination suitable for small and fast solar wind monitors. The first version will be launched during the SPECTR-R project into a highly elongated orbit with apogee ∼350,000 km. The method is based on simultaneous measurements of the total ion flux and ion integral energy spectrum by six identical Faraday cups. Three of them are dedicated to determination of the ion flow direction, the other three (equipped with control grids supplied by a retarding potential) are used for determination of the density, temperature, and speed of the plasma flow. The version under development is primarily designed for the measurements in the solar wind and tail magnetosheath, thus for velocities range from 270 to 750 km/s, temperatures from 1 to 30 eV, and densities up to 200 cm⁻³. However, the instrument design can be simply modified for measurements in other regions with a substantial portion of low-energy plasma as a subsolar magnetosheath, cusp or low-latitude boundary layer. Testing of the engineering model shows that the proposed method can provide reliable plasma parameters with a high time resolution (up to 8 Hz). The paper presents not only the method and its technical realization but it documents all advantages and peculiarities of the suggested approach.     

Heating heavy ions in the polar corona by collisionless shocks: A one-dimensional simulation

Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T_⊥ ? T_∥, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O⁵⁺, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O⁵⁺ ions are easily reflected and gain more than mass proportional energy with respect to protons.     

Solar wind-Venus interaction affecting ionospheric electron density profile

The equilibrium electron density profile has been computed and compared with measured profiles by Venera 9 and Mariner 5 and 10. The computed electron density profile is seen to show discrepancies with measured data. The contribution of solar wind interaction induced convection to equilibrium electron density profile has been estimated. It is found that the convective processes are less important at lower altitudes, whereas at higher altitudes its contribution becomes dominant. The night side Venus ionosphere is formed due to the transport of O⁺ and impact ionization of neutral gases by suprathermal electrons. The discrepancies in theoretical and measured electron density profiles provide clear indication of additional energy source of solar wind origin.     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Changes in sea-level pressure over South Korea associated with high-speed solar wind events

We explore a possibility that the daily sea-level pressure (SLP) over South Korea responds to the high-speed solar wind event. This is of interest in two aspects: first, if there is a statistical association this can be another piece of evidence showing that various meteorological observables indeed respond to variations in the interplanetary environment. Second, this can be a very crucial observational constraint since most models proposed so far are expected to preferentially work in higher latitude regions than the low latitude region studied here. We have examined daily solar wind speed V, daily SLP difference ΔSLP, and daily log(BV²) using the superposed epoch analysis in which the key date is set such that the daily solar wind speed exceeds 800 km s⁻¹. We find that the daily ΔSLP averaged out of 12 events reaches its peak at day +1 and gradually decreases back to its normal level. The amount of positive deviation of ΔSLP is +2.5 hPa. The duration of deviation is a few days. We also find that ΔSLP is well correlated with both the speed of solar wind and log(BV²). The obtained linear correlation coefficients and chance probabilities with one-day lag for two cases are r ? 0.81 with P > 99.9%, and r ? 0.84 with P > 99.9%, respectively. We conclude by briefly discussing future direction to pursue.     

Self-similar stochastic processes in solar wind turbulence

Solar wind data is used to estimate the autocorrelation function for the stochastic process x(τ) = y(t + τ) − y(t), considered as a function of τ, where y(t) is any one of the quantities B²(t), n_p(t)V²(t), or n_p(t). This process has stationary increments and a variance that increases like a power law τ^2γ where γ is the scaling exponent. For the kinetic energy density and the proton density the scaling exponent is close to the Kolmogorov value γ = 1/3, for the magnetic energy density it is slightly larger. In all three cases, it is shown that the autocorrelation function estimated from the data agrees with the theoretical autocorrelation function for a self-similar stochastic process with stationary increments and finite variance. This is far from proof, but it suggests that these stochastic processes may be self-similar for time scales in the small scale inertial range of the turbulence, that is, from approximately 10 to 10³ s.     

Modeling and data analysis of a Forbush decrease

We study the Forbush decrease of the galactic cosmic ray intensity observed in 9–25 September 2005 using the experimental data and a newly developed time-dependent three dimensional modeling. We analyze neutron monitors and muon telescopes, and the interplanetary magnetic field data. We demonstrate a clear relationship between the rigidity (R) spectrum exponent (γ) of the Forbush decrease and the exponent (ν) of the power spectral density of the components of the interplanetary magnetic field in the frequency range of ∼ 10⁻⁶–10 ⁻⁵ Hz. We confirm that an inclusion of the time-dependent changes of the exponent ν makes the newly developed nonstationary three dimensional model of the Forbush decrease compatible with the experimental data. Also, we show that the changes of the rigidity spectrum exponent γ does not depend on the level of convection of the galactic cosmic rays stream by solar wind; depending on the changes of the exponent ν, i.e. on the state of the turbulence of the interplanetary magnetic field.     

Statistical properties of the IMF clock angle in the solar wind with northward and southward interplanetary magnetic field based on ACE observation from 1998 to 2009: Dependence on the temporal scale of the solar wind

Utilizing ACE satellite observations from 1998 to 2009, we performed the elaborate study on the properties of the clock angle θ_CA (arctan(By/Bz) (?90° to 90°) of the interplanetary magnetic field (IMF) in the solar wind at 1?AU. The solar wind with northward IMF (NW-IMF) and southward IMF (SW-IMF) are analyzed, independently. Statistical analysis shows that the solar wind with SW-IMF and NW-IMF has similar properties in general, including their durations, the IMF Bz and By components, and the IMF θ_CA. Then, the solar wind with NW-IMF (SW-IMF) is classified into five different temporal scales according to the duration of the NW-IMF (SW-IMF), i.e., very-short wind of 10–30?min, short-scale wind of 0.5–1?h, moderate-scale wind of 1–3?h, long-scale wind of 3–5?h, and super-long wind >5?h. Our analysis reveals that the IMF θ_CA has a distinct decrease with increase of the temporal scale of the solar wind. Next, the solar wind is classified into two groups, i.e., the high-speed solar wind (>450?km/s) and the low-speed solar wind (<450?km/s). Our analysis indicates that the IMF θ_CA depends highly on the solar wind speed. Statistically, high-speed solar wind tends to have larger IMF θ_CA than low-speed solar wind. The evolutions of the solar wind and IMF with the solar activity are further studied, revealing no clear solar variation of the IMF θ_CA. Finally, we analyze the monthly variation of the IMF θ_CA. Superposed epoch result strongly suggests the seasonal variation of the IMF θ_CA.