Helium escape from the Earth''s atmosphere: The charge exchange mechanism revisited

We have studied the escape of neutral helium from the terrestrial atmosphere through exothermic charge exchange reactions between He⁺ ions and the major atmospheric constituents N₂, O₂, and O. Elastic collisions with the neutral background particles were treated quantitatively using a recently developed kinetic theory approach. An interhemispheric plasma transport model was employed to provide a global distribution of He⁺ ions as a function of altitude, latitude and local solar time and for different levels of solar ionization. Combining these ion densities with neutral densities from an MSIS model and best estimates for the reaction rate coefficients of the charge exchange reactions, we computed the global distribution of the neutral He escape flux. The escape rates show large diurnal and latitudinal variations, while the global average does not vary by more than a factor of three over a solar cycle. We find that this escape mechanism is potentially important for the overall balance of helium in the Earth's atmosphere. However, more accurate values for the reaction rate coefficients of the charge exchange reactions are required to make a definitive assessment of its importance.     

Measurements of radio noise at 0.700 mc and 2.200 mc from a high-altitude rocket probe

Radio noise observations at frequencies of 0·700 Mc and 2·200 Mc were made at altitudes between 3000 and 11,000 km from a Blue Scout Jr. high-altitude rocket probe on 30 July 1963. A steady background flux of (7·5₋₃⁺⁶) × 10⁻¹⁹ W m⁻²)(c/s)⁻¹ at 0·700 Mc and (1·8^+1.0_−0.5 × 10⁻¹⁹ W m⁻² (c/s)⁻¹ at 2·200 Mc was observed. Assuming a galactic origin of the observed fluxes at both frequencies, the averaged sky brightnesses are b(0·700 Mc) = (6₋₃⁺⁵) × 10⁻²⁰ W m⁻² (c/s)⁻¹ sr⁻¹b(2·200 Mc) = (1.4^+1.0_−0.5 × 10⁻²⁰ W m⁻² (c/s)⁻¹ sr⁻¹ The observed brightness at 2·200 Mc is in reasonable agreement with the results of other observers. The apparent brightness at 0·700 Mc is, however, greater than was expected from previous observations. An alternative source of the 0·700 Mc flux in the terrestrial exosphere, as well as characteristics of several noise bursts observed during the flight, is briefly discussed.     

Mapping of the cusp plasma precipitation on the surface of Mercury

The presence of a magnetosphere around Mercury plays a fundamental role on the way the solar wind plasma interacts with the planet. Since the observations suggest that Mercury should occupy a large fraction of its magnetosphere and because of lack of an atmosphere, significant differences in solar wind-magnetosphere coupling are expected to exist with respect to the Earth case. On the basis of a modified Tsyganenko T96 model we describe the geometry of the magnetic field that could characterize Mercury, and its response to the variations of the impinging solar wind and of the interplanetary magnetic field. The investigation is focused on the shape and dimension of the open magnetic field regions (cusps) that allow the direct penetration of magnetosheath plasma through the exosphere of Mercury, down to its surface. The precipitating particle flux and energy are evaluated as a function of the open field line position, according to different solar wind conditions. A target of this study is the evaluation of the sputtered particles from the crust of the planet, and their contribution to the exospheric neutral particle populations. Such estimates are valuable in the frame of a neutral particle analyser to be proposed on board of the ESA/BepiColombo mission.     

Rocket measurements of the 6300 Å and 3914 Å dayglow features

Rocket results are presented on the OI 6300 Å line and on the N₂⁺ 3914 Å band in the dayglow. An altitude range of 78–335 km is covered. Theoretical interpretations are given, using results of simultaneous measurements of electron density and electron temperature. The apparent brightness of the 6300 Å line at the base of the emitting region is found to be 13 kR, of which 5.5 kR are ascribed to excitation through the Schumann-Runge dissociation of O₂ by the solar UV radiations, 0.55 kR to the dissociative recombination of O₂⁺ and NO⁺ ions, and 0.03 kR to the excitation of O by thermal electrons. An additional source of excitation above 280 km is suggested. The deactivation of O(¹D) by O₂(X³Σ_g⁻) is found to be appreciable below 200 km, and its rate coefficient is estimated to be 2 × 10⁻¹⁰ cm³/sec. The apparent brightness of the 3914 Å band at the base of the emitting region is found to be 6.5 kR, decreasing to 3.2 kR at 330 km. Assuming that fluorescent scattering of solar radiation is the mechanism involved the distribution of N₂⁺ ions is calculated. The rate coefficients for the loss of these ions are hence calculated.     

Mercury's sodium exosphere

Mercury's neutral sodium exosphere is simulated using a comprehensive 3D Monte Carlo model following sodium atoms ejected from Mercury's surface by thermal desorption, photon stimulated desorption, micro-meteoroid vaporization and solar wind sputtering. The evolution of the sodium surface density with respect to Mercury's rotation and its motion around the Sun is taken into account by considering enrichment processes due to surface trapping of neutrals and ions and depletion of the sodium available for ejection from the surfaces of grains. The change in the sodium exosphere is calculated during one Mercury year taking into account the variations in the solar radiation pressure, the photo-ionization frequency, the solar wind density, the photon and meteoroid flux intensities, and the surface temperature. Line-of-sight column densities at different phase angles, the supply rate of new sodium, average neutral and ion losses over a Mercury year, surface density distribution and the importance of the different processes of ejection are discussed in this paper. The sodium surface density distribution is found to become significantly nonuniform from day to night sides, from low to high latitudes and from morning to afternoon because of rapid depletion of sodium atoms in the surfaces of grains mainly driven by thermal depletion. The shape of the exosphere, as it would be seen from the Earth, changes drastically with respect to Mercury's heliocentric position. High latitude column density maxima are related to maxima in the sodium surface concentration at high latitudes in Mercury's surface and are not necessarily due to solar wind sputtering. The ratio between the sodium column density on the morning side of Mercury's exosphere and the sodium column density on the afternoon side is consistent with the conclusions of Sprague et al. (1997, Icarus 129, 506-527). The model, which has no fitting parameters, shows surprisingly good agreement with recent observations of Potter et al. (2002, Meteor. Planet. Sci. 8, 3357-3374) successfully explaining their velocity and column density profiles vs. heliocentric distance. Comparison with this data allows us to constrain the supply rate of new sodium atoms to the surface. We also discuss the possible origins of the strong high latitude emissions (Potter and Morgan, 1990, Science 248, 835-838; 1997a, Adv. Space Res. 19, 1571-1576; 1997b, Planet. Space Sci. 45, 95-100; Sprague et al., 1998, Icarus 135, 60-68) and the strong variations of the total content of the sodium exosphere on short (Potter et al., 1999, Planet. Space Sci. 47, 1441-1449) and long time scales (Sprague et al., 1997, Icarus 129, 506-527).     

Cyclotron amplification of geomagnetic micropulsations Pc1 in the magnetosphere

A numerical analysis of cyclotron instabilities is carried out by computing the dispersion relation for a three component cold plasma-beam system. Rates of growth and damping for various values of the stream density are calculated from the dispersion relation. The rates of growth and damping increase monotonically as the number density of the proton stream increases. It is found that the frequencies at the rates of maximum growth and the damping decrease slightly to lower frequencies and a sharp peak at these frequencies becomes blunt. The minimum e-folding times of an ion cyclotron wave for (a) σ_s = 10⁻⁴, σ_i = 10⁻² and (b) σ_s = 10⁻¹, σ_i = 10⁻² are about 3·84 and 0·16 sec respectively in the vicinity of the equatorial plane at 6 Re, where σ_s and σ_i are the ratios of the beam density N_s and the helium ion (H₆⁺) density N_i to the total positive ions in the plasma-beam system.     

Electron transport and precipitation at Mercury during the MESSENGER flybys: Implications for electron-stimulated desorption

To examine electron transport, energization, and precipitation in Mercury's magnetosphere, a hybrid simulation study has been carried out that follows electron trajectories within the global magnetospheric electric and magnetic field configuration of Mercury. We report analysis for two solar-wind parameter conditions corresponding to the first two MESSENGER Mercury flybys on January 14, 2008, and October 6, 2008, which occurred for similar solar wind speed and density but contrasting interplanetary magnetic field (IMF) directions. During the first flyby the IMF had a northward component, while during the second flyby the IMF was southward. Electron trajectories are traced in the fields of global hybrid simulations for the two flybys. Some solar wind electrons follow complex trajectories at or near where dayside reconnection occurs and enter the magnetosphere at these locations. The entry locations depend on the IMF orientation (north or south). As the electrons move through the entry regions they can be energized as they execute non-adiabatic (demagnetized) motion. Some electrons become magnetically trapped and drift around the planet with energies on the order of 1–10 keV. The highest energy of electrons anywhere in the magnetosphere is about 25 keV, consistent with the absence of high-energy (>35 keV) electrons observed during either MESSENGER flyby. Once within the magnetosphere, a fraction of the electrons precipitates at the planetary surface with fluxes on the order of 10⁹ cm⁻² s⁻¹ and with energies of hundreds of eV. This finding has important implications for the viability of electron-stimulated desorption (ESD) as a mechanism for contributing to the formation of the exosphere and heavy ion cloud around Mercury. From laboratory estimates of ESD ion yields, a calculated ion production rate due to ESD at Mercury is found to be on par with ion sputtering yields.     

Helium ions in the mid-latitude plasmasphere

An initial study of the behaviour of He⁺ ions in the mid-latitude plasmasphere is carried out by solving the time-dependent equations of continuity and momentum. Starting with a low He⁺ tube content, results are obtained for a period of 8 days. In the topside ionosphere there is an upflow of He⁺ during the day and downflow at night, for the sunspot maximum conditions considered. The downflow at night differs from the behaviour of H⁺ for these atmospheric conditions. However, little of the He⁺ produced in the daytime is lost by recombination at night; it is suggested that the supply of He⁺ to the mid-latitude plasmasphere is, in effect, an escape process for neutral helium.     

Minor ions in the solar wind

Ions heavier than ⁴He are treated as “minors” in the solar wind. This is justified for many applications since minor ions have no significant influence on the dynamics of the interplanetary plasma. However, minor ions carry information on many aspects of the formation, on the acceleration and on the transfer of solar plasma from the corona into the interplanetary space. This review concentrates on various aspects of minor ions as diagnostic tracers. The elemental abundance patterns of the solar wind are shaped in the chromosphere and in the lower transition region by processes, which are not fully understood at this moment. Despite this lack of detailed understanding, observed abundance patterns have been classified and are now commonly used to characterize the sources, and to trace back solar-wind flows to their origins in the solar atmosphere. Furthermore, the solar wind is the most important source of information for solar isotopic abundances and for solar abundances of volatile elements. In order to fully exploit this information, a comprehensive understanding of elemental and isotopic fractionation processes is required. We provide observational clues to distinguish different processes at work.     

Helioseismological constraint on solar axion emission

Helioseismological sound-speed profiles severely constrain possible deviations from standard solar models, allowing us to derive new limits on anomalous solar energy losses by the Primakoff emission of axions. For an axion-photon coupling g_ay 5 × 10⁻¹⁰ GeV⁻¹, the solar model is almost indistinguishable from the standard case, while g_ay 10 × 10⁻¹⁰ GeV⁻¹ is probably excluded, corresponding to an axion luminosity of about 0.20 L. This constraint on g_ay is much weaker than the well-known globular-cluster limit, but about a factor of 3 more restrictive than previous solar limits. Our result is primarily of interest to the large number of current or proposed search experiments for solar axions because our limit defines the maximum g_ay for which it is self-consistent to use a standard solar model to calculate the axion luminosity.     

AXIOM: Advanced X‐ray imaging of the magnetosheath

AXIOM (Advanced X‐ray Imaging Of the Magnetosphere) is a concept mission which aims to explain how the Earth's magnetosphere responds to the changing impact of the solar wind using a unique method never attempted before; performing wide‐field soft X‐ray imaging and spectroscopy of the magnetosheath, magnetopause and bow shock at high spatial and temporal resolution. Global imaging of these regions is possible because of the solar wind charge exchange (SWCX) process which produces elevated soft X‐ray emission from the interaction of high charge‐state solar wind ions with primarily neutral hydrogen in the Earth's exosphere and near‐interplanetary space (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Estimation prospects of the source number density of ultra-high-energy cosmic rays

We discuss the possibility of accurately estimating the source number density of ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their arrival distribution. The arrival distribution has information on their source and source distribution. We calculate the propagation of UHE protons in a structured extragalactic magnetic field (EGMF) and simulate their arrival distribution at the Earth using our previously developed method. The source number density that can best reproduce observational results by Akeno Giant Air Shower Array is estimated at about 10⁻⁵ Mpc⁻³ in a simple source model. Despite having large uncertainties of about one order of magnitude, due to small number of observed events in current status, we find that more detection of UHECRs in the Auger era can sufficiently decrease this so that the source number density can be more robustly estimated. Two hundred event observation above 4 × 10¹⁹ eV in a hemisphere can discriminate between 10⁻⁵ and 10⁻⁶ Mpc⁻³. Number of events to discriminate between 10⁻⁴ and 10⁻⁵ Mpc⁻³ is dependent on EGMF strength. We also discuss the same in another source model in this paper.     

Upper limits for the Martian exospheric number density during the Planet B/Nozomi mission

Dissociative recombination (DR) of ionospheric O₂⁺ ions is an important source of suprathermal atomic oxygen in the exosphere as previous studies about the Martian upper atmosphere have shown. Because of the weaker gravitational attraction a hot oxygen corona on Mars should be denser than that observed on Venus. Since the most important mechanism for the production of the hot oxygen atoms in the Martian exosphere is DR, we investigated the variability of this production mechanism depending of solar activity. The Japanese Nozomi spacecraft will have the possibility to detect with the neutral mass spectrometer (NMS) for the first time in-situ the theoretically predicted hot oxygen corona on Mars, if the corona number density above the cold background atmosphere is of the order of 10,000 cm⁻³. Due to a problem in the propulsion system Nozomi failed its planned arrival rendevouzs with Mars in October 1999 and will, therefore, arrive at the red planet not before January 2004. Solar activity will reach its maximum in 2001, so the related production rate of hot oxygen atoms will be in the medium range during the new arrival date of Nozomi. We used the ionospheric profiles from the Viking mission for low solar activity conditions (F_10.7≈70) and the Mariner 9 mission with a solar activity of about 120 for medium solar wind activity. The latter is comparable to the level we expect for the Mars arrival of Nozomi. The resulting influence of the hot oxygen corona number density distribution was calculated with a Monte Carlo technique. This technique is used to compute a hot particle density distribution function. We studied the atomic diffusion process in the Martian atmosphere by simulating the collision probability, particle direction and energy loss after collisions by generating random numbers. Compared to previous studies we have improved the Monte Carlo model by using more and smaller altitude steps and more detailed treatment of particles with a temporary downward motion. This has resulted in an increased amount of collisions and a shift to lower energies in the energy spectrum. Our results show that the hot oxygen component should begin to dominate above the cold background atmosphere at an altitude of about 500 km above the Martian surface. The NMS instrument on board of Nozomi should detect the hot oxygen component after its arrival at Mars in January 2004, at an altitude of about 600 km above the Martian surface. Since the solar activity will decrease during the mission the measurements during the first orbits will be the most significant ones. The first in-situ measurements of the hot oxygen number density would be very important for adjusting atmospheric escape models by separating ballistic, satellite and escape trajectories of the hot oxygen atoms, which are significant for studies of the evolution and solar wind interaction of the Martian atmosphere.     

Interaction of interstellar pick-up ions with the solar wind

The interaction of interstellar pick-up ions with the solar wind is studied by comparing a model for the velocity distribution function of pick-up ions with actual measurements of He⁺ ions in the solar wind. The model includes the effects of pitch-ang'e diffusion due to interplanetary Alfvén waves, adiabatic deceleration in the expanding solar wind and the radial variation of the source function. It is demonstrated that the scattering mean free path is in the range 0.1 AU and that energy diffusion can be neglected as compared with adiabatic deceleration. The effects of adiabatic focusing, of the radial variation of the neutral density and of a variation of the solar wind velocity with distance from the Sun are investigated. With the correct choice of these parameters we can model the measured energy spectra of the pick-up ions reasonably well. It is shown that the measured differential energy density of the pick-up ions does not vary with the solar wind velocity and the direction of the interplanetary magnetic field for a given local neutral gas density and ionization rate. Therefore, the comparison of the model distributions with the measurements leads to a quantitative determination of the local interstellar gas density.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

Comparative Study of MHD Modeling of the Background Solar Wind

Knowledge about the background solar wind plays a crucial role in the framework of space-weather forecasting. In-situ measurements of the background solar wind are only available for a few points in the heliosphere where spacecraft are located, therefore we have to rely on heliospheric models to derive the distribution of solar-wind parameters in interplanetary space. We test the performance of different solar-wind models, namely Magnetohydrodynamic Algorithm outside a Sphere/ENLIL (MAS/ENLIL), Wang–Sheeley–Arge/ENLIL (WSA/ENLIL), and MAS/MAS, by comparing model results with in-situ measurements from spacecraft located at 1 AU distance to the Sun (ACE, Wind). To exclude the influence of interplanetary coronal mass ejections (ICMEs), we chose the year 2007 as a time period with low solar activity for our comparison. We found that the general structure of the background solar wind is well reproduced by all models. The best model results were obtained for the parameter solar-wind speed. However, the predicted arrival times of high-speed solar-wind streams have typical uncertainties of the order of about one day. Comparison of model runs with synoptic magnetic maps from different observatories revealed that the choice of the synoptic map significantly affects the model performance.     

Nightside ionosphere of Venus

Venera 9, 10 measurements of the nightside ionospheric profile and the night airglow were used for investigating ionosphere formation processes. The upper ionospheric layer may be formed by HeI 584 Å radiation; the lower layer by meteorite ionization. Upper limits on the electron energy flux, <4 × 10⁸eV cm⁻² s⁻¹, the helium ion flux <10⁷ cm⁻² s⁻¹, the nitric oxide mixing ratio, <1.5 × 10⁻⁴ and the atomic sulphur mixing ratio, <10⁻⁶, are deduced for ionospheric altitudes.     

Hydrogen atoms and ions in the thermosphere and exosphere

The distribution of atomic hydrogen in the thermosphere and exosphere is computed taking into account the upward flow which balances the escape flux. Because of the upward flow the number-density gradient is much steeper than it would be in a static atmosphere. Attention is drawn to the fact that the ratio of the amount of hydrogen above the 100 or 110km levels to the amount of hydrogen above the 200 or 300 km levels is a sensitive measure of the temperature of the exosphere. The evidence on the absolute abundance of atomic hydrogen is examined. It is concluded that the number density at the 120km level is probably about 5 × 10⁵/cm³. The Ly. absorption line at this level is beyond the linear part of the curve of growth.

Consideration is also given to the steady-state distributions of O⁺ and H⁺ ions. In the lower part of the exosphere the number density of O⁺ ions falls with increase in altitude (the associated scale height being twice that of the O atoms) and the number density of H⁺ ions rises at the same rate (as was first pointed out by Dungey). The altitude at which the number densities of O⁺ and H⁺ ions become equal is calculated on various assumptions regarding the temperature and hydrogen content of the exosphere. It is found to be about 1200 km when the temperature is 1250° K and the hydrogen content corresponds to the number density cited near the end of the preceding paragraph. The gradient of the predicted electrondensity distribution at several Earth radii is much less than that deduced from whistler studies.

The passage from charge transfer to diffusive equilibrium is discussed in an Appendix.     

The shape of the power spectrum of cosmic ray at ground level up to 7 × 10 Hz

The power spectral density of cosmic ray fluctuations observed at ground level during the years 1966–1968 has been calculated. In order to obtain the correct shape of the spectrum, the Fast Fourier Transform method with a triangular data window was used and corrections were made for uncorrelated errors and aliasing effects. When ignoring the Earth rotation peaks, the spectral index , for a sample of polar, middle latitude and equatorial stations, is − 1.96 in the frequency range 3 × 10⁻⁷–10⁻⁴ Hz. A possible break around 10⁻⁵ Hz, if existing, would be, on the whole, barely significant as a would change from − 1.96 to − 2.10. There are indications that beyond 10⁻⁴ Hz up to 7 × 10⁻³ Hz the spectrum continues with − 2.     

Magnetosphere,exosphere, and surface of mercury

The discovery of an atomic sodium exosphere at Mercury raises the question of whether Mercury, like Io at Jupiter, can maintain a heavy ion magnetosphere. We suggest that it does, and that heavy ions (mainly Na⁺) from the exosphere are typically accelerated to keV energies and make important or dominant contributions to the mass (∼300 g sec⁻¹) and energy (∼3 × 10⁹W) budgets of the magnetosphere. The sodium supply to the exosphere is largely from within Mercury itself, with external sources like meteroid infall and the solar wind being relatively unimportant. Therefore Mercury is in the process of losing its semivolatiles. Photosputtering dominates charged particle sputtering and can maintain an adequate rate of Na ejection from the surface.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.