A&A 465, 641-645 (2007)
DOI: 10.1051/0004-6361:20065991
G. Gronoff1 - J. Lilensten1 - C. Simon1 - O. Witasse2 - R. Thissen3 - O. Dutuit3 - C. Alcaraz3
1 - Laboratoire de Planétologie de Grenoble, Université Joseph Fourier - CNRS, France
2 -
Research and Scientific Support Department of ESA, ESTEC, Noordwijk, The Netherlands
3 -
Laboratoire de Chimie Physique - CNRS - UMR 8000, Univ. Paris Sud, 91405 Orsay, France
Received 7 July 2006 / Accepted 9 December 2006
Abstract
Context. Previous studies have shown that doubly-charged positive ions, including molecular ions, can have detectable densities in the ionospheres of Mars, Titan, and the Earth.
Aims. In the present paper, we continue our approach by modelling the Venusian ionosphere. The neutral atmosphere of Venus has a similar composition to the ones encountered in the other terrestrial planets, with carbon dioxide (as on Mars), nitrogen (as on Earth and Titan), and atomic oxygen (as on Earth).
Methods. We computed the doubly-charged positive ion production through photoionisation (primary production) and electron impact ionisation (secondary production). We computed the densities under the photochemical equilibrium assumption.
Results. We predict a
layer centred around
altitude with a density reaching 30 dications per
in active solar conditions and a
layer centred around a
altitude with density reaching
.
We find good agreement between the modelled
densities and Pioneer Venus Orbiter measurements. Finally, we investigate the problem of possible detection of these ions by the European spacecraft Venus Express.
Conclusions. Although molecular doubly-charged positive ion densities are low, these ions cannot be considered as negligible compared to other minor ion species. Their detection should be possible in the future through remote sensing or in-situ methods including ion mass spectrometry and UV spectroscopy.
Key words: planets and satellites: individual: Venus - atmospheric effects - Sun: UV radiation - space vehicles: instruments
The study of planetary ionospheres is among the major scientific goals of many planetary missions. Venus Express (VEX), a European Space Agency mission, is the first dedicated mission to Venus since Magellan in the 1990s.
The VEX spacecraft entered into the Venus orbit on 11 April 2006, and it has since then started acquiring its planetary data, thereby providing a global view of Venus' atmosphere. The scientific goals cover the following themes: atmosphere dynamics, structure, composition, clouds and haze, radiative balance, surface properties, and - most important in the frame of this paper - plasma/escape processes. The study of the ionosphere of Venus is clearly included in the last theme, and the following experiments are expected to contribute: VERA (Häusler et al. 2006), SPICAV (Bertaux et al. 2006), and ASPERA (Barabash et al. 2006). A review of our current knowledge of certain key topics associated with the aeronomy of Venus can be found in Witasse & Nagy (2006). We report here on a study of the modelling of doubly-charged positive ions in the diurnal ionosphere.
Plasma densities depend on production
and loss mechanisms, as well as on the dynamics.
Production can be of either a chemical or physical
origin. In the latter case we consider
two sources: photoabsorption (called primary photoproduction) and
electron impact ionisation (called secondary production).
In optically thin upper atmospheres, as in the case for Venus in the EUV range (<
)
above
,
the primary production is computed through a Beer-Lambert law.
The secondary production computation
requires solving a kinetic transport equation (Lilensten & Blelly 2002).
In Witasse et al. (2002,2003), this approach was adapted to Mars,
allowing the authors to predict the existence of a
layer.
The same approach was developed by Simon et al. (2005) and Lilensten et al. (2005) for the Earth and Titan cases, respectively.
The present paper aims at presenting the calculation of the three main stable dications in the Venusian ionosphere:
,
,
and
.
The
dication has already been detected and quantified in the Venusian ionosphere by Pioneer
Venus Orbiter (PVO) (Taylor et al. 1980), eventually allowing comparisons to be made with models such as the one of Fox & Victor (1981).
The
and
molecular ions have not yet been identified and are modelled here for the first time in the ionosphere of Venus.
Several input parameters must be taken into account for modelling the contribution doubly-charged ions in terms of productions and densities: the solar flux and the neutral atmosphere.
We use the semi-empirical solar EUV/XUV model called solar 2000 (Tobiska & Eparvier 1998; Tobiska 1991,1993). Like most models currently used in planetary aeronomy (Torr & Torr 1985,1979; Richards et al. 1994), this model is based on measurements recorded onboard the Dynamics Explorer satellites (Hinteregger et al. 1973) and onboard rockets. It also takes data into account from other sources, such as the SME, OSO, and AEROS spacecrafts, as well as ground-based facilities. The UV spectrum is divided into 39 energy boxes.
The neutral atmosphere is yielded by the VTS3 model (Hedin et al. 1983), restricted to the main species (
,
,
,
,
).
This model is run for the conditions that match the orbit in terms of solar activity, location, time, and solar-zenith angle.
The
molecule is the major constituent below
altitude,
below
,
while
and
become dominant above.
The
density is about
above
altitude
and does not affect the results below
(Fox & Sung 2001).
However, as will be shown later, the evolution of hydrogen density has an influence on the
density above this altitude. In order to take its variations into account, we use the results of Hartle et al. (1996) who examined the solar cycle variations of the hydrogen in the thermosphere of Venus.
The altitude range for our computations is
altitude.
The dication density computed for a neutral species A is given by
We consider 32 ionisation processes, summarised in Table 1.
The cross sections for the single ionisation of neutral species are explained in Witasse et al. (2002,2003), Lilensten et al. (2005), Simon et al. (2005), and references therein.
In Table 2, we give the references for the double-ionisation cross sections of the different neutral species.
The
double-photoionisation cross section and its evolution with photon energy was estimated on the basis of a method proposed in Samson (1990).
Samson shows that the double-photoionisation cross section
of a species A (
)
is proportional to the electron-impact ionisation
of the monocation (A
), multiplied by the total photoabsorption
of the neutral species (
.
The proportionality constant q is estimated as
in the case of
.
Thus, we have:
.
Table 1: Summary of the ion productions taken into account here and their corresponding ionisation thresholds.
Table 2: Sources for double-ionisation cross sections, either from theoretical estimates or from laboratory measurements.
The Pioneer Venus Orbiter carried 17 experiments, including a neutral mass spectrometer and an ion-mass spectrometer (Russell 1992).
It allowed Hedin et al. (1983) to build a model of the Venusian thermosphere. However, although the mission lasted 14 years, no low-altitude data were being recorded during low solar activity, which put some constraints and uncertainties on the thermospheric model.
The Ion and Neutral Mass Spectrometer (INMS) allowed the detection of atomic doubly-charged ions such as
(Taylor et al. 1980). However, the spectrometer was designed to detect only specific, preselected ions excluding
.
The sensitivity threshold of INMS was 5 particles per
.
The
dication cannot be distinguished from
because both ions have the same
ratio (where m is the mass and q the charge).
We have selected three orbits in the PVO data base, which represent three different geophysical conditions:
![]() |
Figure 1:
The
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The results are shown in Fig. 1.
For orbit 185, the
density was modelled by Fox & Victor (1981).
Using a photoionisation model, Fox computed a reaction rate constant of
for the reaction
.
This model took neither the electron impact ionisation nor the electron recombination of
into account. More recently, the Simon et al. (2005) calculations on Earth showed a better agreement considering a constant rate of
.
We considered both values here, and found much better agreement with the second one.
We did not need to adjust the other input parameters, like the input solar flux, the chemical reaction rates, the absorption and collision cross sections, and the neutral atmosphere.
The agreement is already very good in orbits 185 and 441 when the
density used in the model comes from PVO measurements with
above
.
This
density value is compatible with the solar conditions during orbit 185 as seen from Hartle et al. (1996). When this value is enhanced to
,
the model and measurements diverge above
.
When neglecting the single ionisation of
(Fig. 1) a worse divergence occurs above the same altitude.
At low altitudes, the model slightly overestimates the
density. According to Fox & Sung (2001), the atomic
mixing ratio can be too high at high solar activity in the VTS3 model (Hedin et al. 1983), explaining this small discrepancy.
The 14th orbit constitutes a boundary condition. The latitude is about
,
while the solar zenith angle is
.
In that case, we are at the planet's terminator. Moreover, the solar activity is lower (
). In these conditions, the
density is likely to be
(Hartle et al. 1996). With this value, our model fits the measurements perfectly, while a lower
density value (
)
results in a degradation above
.
From these three cases, we have highlighted the usability range in which our model is valid.
![]() |
Figure 2:
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The Venus Express SPICAV ultraviolet spectrometer has a range of
with a resolution of
;
therefore, the detection of such a fluorescence may be difficult but not impossible.
The
dication has no significant known fluorescence line. Therefore, the most appropriate method of detection is the ion-mass
spectrometry technique. For
the mass-to-charge ratio equals 22. Therefore this ion is a very good
candidate for detection, since no other species in the Cytherean
atmosphere has the same ratio. The closest ratios are 16 for
and 28 for both
and
.
The only species with the same
mass-to-charge ratio is the
ion. However the main isotope of
has an m/q ratio of 20, and its density in the Cytherean ionosphere is known to be extremely low, well below the density of
predicted here.
Thus no contamination from this component is likely to occur.
Finally, the Analyser of Space Plasma and Energetic Atoms (ASPERA-4) onboard Venus Express was not designed to detect
.
For the first time, the density of molecular dications has been
modelled in the diurnal upper atmosphere of Venus.
The
density was also computed and successfully compared to PVO measurements. These densities have been calculated
by using a kinetic code for the production rates and a simple chemical
scheme for the chemical losses. Doubly-charged ions were produced on the dayside through
primary photoionisation and secondary electron impact ionisation.
In the subsolar case, an additional source of ionisation may be solar-wind electrons and protons precipitating in the ionosphere and colliding with the thermosphere. These ions are essentially lost by dissociative recombination
with thermal electrons and by chemical reactions with the neutral atmosphere.
Consequently, a molecular dication layer can be created. The
dication layer is centred on an altitude of
with a density reaching
in active solar conditions.
The density of the
layer centred at
altitude can reach
.
The
dication is most abundant at about
,
with typical values of
.
In the near future, we plan to extend our approach to electron precipitation from the solar-wind, a process which can play a major role in the upper atmosphere, depending on the solar activity.
The measurement of molecular doubly-charged ions is a challenge for Venus Express and any planetary orbiter. The detection of such ions is of prime interest for the community involved in aeronomy studies that seek to understand their role in planetary ionospheres.
Acknowledgements
The authors wish to thank Mathieu Barthélemy (LPG, France) for helpful discussions and Héloise Soldi-Lose (Institut fur Chemie, Berlin, Germany) for useful advices.