First composition measurements of positive ions in the upper troposphere

First mass-spectrometric composition measurements of atmospheric ions between 3250 and 11700 m altitude are reported. They reveal the presence of very massive cluster ions, the majority of which cannot be attributed to a single hydrated ion family like, for example H⁺(H₂O)_n. The observed fraction of very massive ions increases with decreasing altitude. Masses as large as about 540 amu were observed at 8200 m altitude. Implications of the observations for ion and nucleation processes are discussed.     

The upper ionospheres of Jupiter and Saturn

We use a 1-D chemical diffusive model, in conjunction with the measured neutral atmospheric structure, to analyze the Voyager RSS electron density, n_e, profiles for the ionospheres of Jupiter and Saturn. As with previous studies we find serious difficulties in explaining the n_e measurements. The model calculates ionospheres for both Jupiter and Saturn with n_e peaks of 10 times the measured peaks at altitudes which are 900–1000 km lower than the altitude of peaks in the RSS electron densities. Based on our knowledge of neutral atmospheric structure, ionization sources, and known recombination mechanisms it seems that, vibrational excitation of H₂ must play some role in the conversion of slowly radiatively recombining H⁺ ions to the relatively more rapidly recombining H₂⁺ and H₃⁺ ions. In addition, vertical ion flow induced by horizontal neutral winds or electric fields probably also play some role in maintaining the plasma peaks observed both for Jupiter and Saturn to be at high altitudes. For the ionosphere of Saturn, the electron densities are affected by a putative influx of H₂O molecules, Φ_H2O, from the rings. To reproduce the RSS V2 exit n_e results model requires an influx of Φ_H2O 2 × 10⁷ molecules cm⁻² s⁻¹ without invoking H_2f vibrational excitation. To maintain the model n_e peak at the measured altitude vertical plasma drift maintained by meridional winds or vertical electric fields is required. The amounts of H₂O are consistent with earlier estimates of Connerney and Waite (1984) and do not violate any observational constraints.     

Response of mesopheric ozone to particle precipitation

In the mesosphere, water vapor photolysis is the major source of odd hydrogen (H, OH and HO₂) under normal conditions. The odd hydrogen produced may then be converted to H₂ by the reaction H + HO₂→ H₂ + O₂. This process is responsible for the calculated decrease in the H₂O mixing ratio and accompanying increase in the H₂ mixing ratio with altitude in the upper mesosphere and lower thermosphere. Charged particle precipitation events are calculated to produce the same effect, particularly in the 70–85 km region, thus temporarily resulting in enhanced conversion of H₂O to H₂ following such an event. Since odd hydrogen is produced predominantly by water vapor photolysis at these altitudes, decreased odd hydrogen concentrations are also anticipated. Odd hydrogen processes dominate ozone destruction in this region, and so an increase in ozone may occur if odd hydrogen concentrations decrease. We have examined the calculated time behavior of these processes in a numerical model using the August 1972 solar proton event as an example, and we present calculations indicating what might be observed in future events.     

Properties of energetic water-group ions in the extended pick-up region surrounding comet Giacobini-Zinner

The properties of energetic (65–95 keV) cometary water-group ions in the extended solar wind pick-up region surrounding comet Giacobini-Zinner are examined using data from the EPAS instrument on the ICE spacecraft. In the outer part of this region, extending from cometocentric distances of several hundred thousand to a few million kilometres (the limit of pick-up ion detectability), it is found that large modulations of the ion flux occur (with J_MAX/J_MIN 10²-10³) which are related to the direction of the magnetic field. It is also found that the ions stream in a direction which is intermediate between the directions of the solar wind flow and the E × B drift, and that ions are present at energies somewhat above the local pick-up energy. These properties indicate that the waves which are excited by the unstable “ring-beam” pick-up ion velocity distributions do result in significant scattering of the ions in this region, both in pitch angle and in energy, but that they have insufficient amplitude to scatter the ions into near isotropy in the solar wind frame. Closer to the comet (but still upstream from the bow shock), the ion flux modulations are considerably reduced in amplitude and the ions respond less to the E × B drift, indicating that the ions are scattered nearer to isotropy in this region. Inbound, this transition takes place relatively abruptly at a distance of 4 × 10⁵ km in association with an increase in the solar wind speed, after which the ion flux increases, and ceases to be modulated by the field direction, while the streaming direction is continuously antisolar and unmodulated by the direction of the E × B drift. Outbound, weak vestiges of the ring-beam ion anisotropy are present in the region immediately upstream from the bow shock (at −1 × 10⁵ km), but these become more marked at distances in excess of t4 × 10⁵ km, increasing gradually with increasing distance from the comet. It is shown that the evolution of the ion properties is qualitatively consistent with expectations based on quasi-linear diffusion of the ions by the magnetosonic waves observed during the encounter.     

Ionospheric conditions

A general analysis of ionospheric conditions has been made in the light of possible ionic reactions occurring in the upper atmosphere. Data obtained on various parameters, such as ionic production and recombination, show that precise knowledge of the spectral distribution of solar radiation is needed and that other experimental determinations on dissociative recombinations are required.

The ionic complexity of the ionosphere is underlined by describing how the atomic ions O⁺ and N⁺ react with N₂, O₂ and NO molecules. The behavior of the molecular ions N⁺₂, O⁺₂and NO⁺depends on a group of simultaneous processes involving charge transfers and ionatom interchanges which are more important than dissociative recombinations. The altitude distribution of ions is exemplified by discussing the relative importance of various loss coefficients in the D-, E- and F-regions. It is seen that molecular nitrogen ions are subject to important charge transfer processes, that nitric oxide ions are always final products destroyed only by dissociative recombination. Additionally, the entire production of atomic oxygen ions is related to the photoionization of molecular nitrogen. Some information is also given on possible anomalies in the ratio of O⁺₂ and NO⁺ densities in the lower ionosphere. From the lack of sufficient experimental information on ionic processes it is shown that a precise analysis of ionospheric behavior remains highly speculative.     

C-type shocks in the interstellar medium: profiles of CH+ and CH absorption lines

We have computed optical absorption-line profiles of CH⁺ and CH, as predicted by a model of a C-type shock propagating in a diffuse interstellar cloud. Both these species are produced in the shock wave in the reaction sequence that is initiated by C⁺(H₂, H)CH⁺. Whilst CH⁺ flows at the ion speed, CH, which forms in the dissociative recombination reaction CH⁺₃(e⁻, H₂)CH, flows at a speed which is intermediate between those of the ions and the neutrals. The predicted velocity shift between the CH⁺ and CH line profiles is found to be no more than approximately 2 km s⁻¹, which is smaller than has previously been assumed. We also investigate OH and HCO⁺, finding that the correlation between their column densities, recently observed in the diffuse interstellar medium, can be reproduced by the model.     

Electron-impact rotational excitation of linear molecular ions

Molecular R -matrix calculations are performed to give rotational excitation rates for electron collisions with linear molecular ions. Results are presented for CO⁺, HCO⁺, NO⁺ and H₂⁺ up to electron temperatures of 10 000 K. De-excitation rates and critical electron densities are also given. It is shown that the widely used Coulomb–Born approximation is valid for Δ j =1 transitions when the molecular ion has a dipole greater than about 2D, but otherwise is not reliable for studying electron-impact rotational excitation. In particular, transitions with Δ j >1 are found to have appreciable rates and are found to be entirely dominated by short-range effects.     

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.     

Dynamics of the O H2O charge-transfer reaction and its implications for space-borne measurements

O⁺(⁴S) + H₂O charge-transfer cross-section and product-ion time-of-flight measurements are presented over a center-of-mass collision energy range of 0.2–15 eV. The results are obtained with a newly constructed guided-ion beam experiment. The measured energy dependence of the charge-transfer crosssection agrees well with previous measurements and ADO model predictions of Bates [Chem. Phys. Lett. 82, 396 (1981) Proc. R. Soc. Lond. A384, 289 (1982) Chem. Phys. Lett. 111, 428 (1984)] at low collision energies, but exhibits significant differences at collision energies above 2 eV. The product kinetic energy analysis shows that most of the product ions are produced with near-thermal kinetic energies. At low collision energies, a significant backscattered channel is observed that is associated with complex formation. This channel exhibits a cross-section of approximately 1.4 Ųat a collision energy of 2.65 eV, corresponding to a laboratory ion energy of 5 eV.     

The winter anomaly in ionospheric absorption and the D-region ion chemistry

A detailed analysis of the D-region ion composition measurements performed by Zbinden et al. (1975), during a winter day of high ionospheric absorption, has been carried out. The study examines the interactive mesosphere-D-region processes which occur in such anomalous conditions and their implication for water cluster ion chemistry. Two clustering regimes for NO⁺ have been observed in the data. Association with N₂ is identified as the dominant process below 76 km. Between 76 and 78 km altitude the effective loss rate of NO⁺ drops by two orders of magnitude. Above 77 km, the three-body reaction NO⁺ + CO₂+M→NO⁺CO₂+M seems to be the main NO⁺ loss. A mesospheric temperature profile could be derived from the ion composition data. This indicates the presence of a strong inversion above 76 km altitude. The wavelike structure obtained, is shown to be consistent with in situ winter temperature measurements. The sharp suppression of the N₂ association reaction could, thus, be explained by an increase in the collisional break-up of the NO⁺N₂ ion because of the enhanced temperature. In conclusion, our study indicates that, besides the increase in the production of NO⁺ and O₂⁺, due to an enhancement in the minor ionizable constituents, an additional thermal mesosphere-D-region interaction seems necessary to explain this winter anomalous ion composition data.     

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.     

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.     

Resonant interactions between cometary ions and low frequency electromagnetic waves

We explore the conditions for resonance between cometary pick-up ions and parallel propagating electromagnetic waves. A model ring—beam distribution for the pick-up H₂O⁺ ions is adopted which allows a direct comparison of the source of free energy for growth from either the beam or the gyrating ring in the limit near marginal stability. Under average solar wind conditions in the inner solar system, the gyrating ring provides the dominant contribution to wave growth. The presence of a field-aligned beam is only important to allow resonance with R-mode waves which occur in two distinct frequency bands either well above or below the pick-up ion gyrofrequency. The most unstable mode is the low frequency R-mode or fast MHD wave, though higher frequency whistlers or low frequency L-mode waves may also be excited by the same source of free energy. The nature of the unstable waves is strongly influenced by the inclination of the interplanetary field. For 3° the rate of the low frequency R-mode growth is dramatically reduced and resonant L-mode waves should experience net ion beam damping. Conversely for 75°, the ion beam velocity will be insufficient to allow resonant R-mode instability; L-mode waves should therefore predominate. The low frequency fast MHD mode should experience the most rapid amplification for intermediate inclination; 30° 75°. In the frame of the solar wind such waves must propagate along the field in the direction upstream towards the Sun with a phase speed lower than the beaming velocity of the pick-up ions. The waves are consequently blown back away from the Sun and would thus be detected with a left-hand polarization by an observer in the cometary frame. We consider this the most likely mechanism to account for the interior MHD waves observed by satellites over an extended spatial region surrounding comets Giacobini-Zinner and Halley.     

An auroral F-region study using in situ measurements by the Atmosphere Explorer-C satellite

On 14 July 1974 the Atmosphere Explorer-C satellite flew through an aurora at F-region altitudes just after local midnight. The effects of the particle influx are clearly evident in the ion densities, the 6300 Å airglow, and the electron and ion temperatures. This event provided an opportunity to study the agreement between the observed ion densities and those calculated from photochemical theory using in situ measurements of such atmospheric parameters as the neutral densities and the differential electron energy spectra obtained along the satellite track. Good agreement is obtained for the ions O₂⁺, NO⁺ and N₂⁺ using photochemical theory and measured rate constants and electron impact cross sections. Atomic nitrogen densities are calculated from the observed [NO⁺]/[O₂⁺] ratio. In the region of most intense electron fluxes (20 erg cm⁻² sec⁻¹) at 280 km, the N density is found to be between 2 and 7 × 10⁷ cm⁻³. The resulting N densities are found to account for approx. 60% of the production of N⁺ through electron impact on N and the resonant charge exchange of O⁺(²P) with N(⁴S). This reaction also provides a significant source of O(¹S) in the aurora at F-region altitudes. In the region of intense fast electron influx, the reaction with atomic nitrogen is found to be the main loss of O⁺(²P).     

First model of negative ion composition in the troposphere

A chemical model of negative ions in the troposhere (0–15 km) is presented. This model is an extension of the negative ion composition model in the lower stratosphere (Kawamoto and Ogawa, 1984, Planet. Space Sci. 32, 1223) with some modifications. The computed result shows that the predominant ions are NO₃⁻HNO₃H₂O below 10km and NO⁻₃(HNO₃)₂ above 10km, and that the fractional abundance of cluster ions having a HSO⁻₄ core increases with height below 12km and decreases with height above it. The ions having CO⁻₃ cores are at most 2% in fractional abundance. The other kinds of negative ions are far smaller in fractional abundance than the NO⁻₃, HSO⁻₄ and CO⁻₃ core ions. The result is compared with the two mass spectrometric observed results (Heitmann and Arnold, 1983, Nature, Lond. 306, 747; Perkins and Eisele, 1984, J. geophys. Res. 89, 9649). The problems on the tropospheric negative ions which arose are discussed.     

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.     

Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

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.     

Atmospheric transmission at submillimetre wavelengths from Mauna Kea

The submillimetre atmospheric transmission spectrum above Mauna Kea has been measured at a resolution of 0.005 cm⁻¹ (150 MHz) with a Fourier transform spectrometer at the James Clerk Maxwell Telescope, using the Sun as a source. Column abundances of O₂, H₂O and O₃ determined from these spectra are found to be in excellent agreement with independent measurements. The derived column abundances have been used as inputs to the atmospheric spectral modelling program fascod . The synthetic transmission spectrum is found to be in excellent agreement with the measured spectrum, and provides a template for submillimetre observations from the JCMT.     

Chemistry of the nightside ionosphere of Venus

We have constructed a one-dimensional model of the nightside ionosphere of Venus in which it is assumed that the ionization is maintained by day-to-night transport of atomic ions. Downward fluxes of O⁺, C⁺ and N⁺ in the ratios measured on the dayside at high altitudes are imposed at the upper boundary of the model (about 235 km). We discuss the resulting sources and sinks of the molecular ions NO⁺,CO⁺,N₂⁺,CO₂⁺ and O₂⁺. As the O⁺ flux is increased, the peak density of O⁺ increases proportionally and the altitude of the peak decreases. The O₂⁺ peak density is approximately proportional to the square root of the O⁺ flux and the peak rises as the O⁺ flux increases. NO⁺ densities near the peak are relatively unaffected by changes in the O⁺ flux. If the ionosphere is maintained mostly by transport, the ratio of the peak densities of O⁺ and O₂⁺ indicates the downward flux ofO⁺, independent of the absolute magnitudes of the densities. The densities of mass-28 ions are, however, still considered to be the most sensitive indicator of the importance of electron precipitation. We examine here the inbound and outbound portions of six early nightside orbits with low periapsis and use data from the Pioneer Venus orbiter ion mass spectrometer, the retarding potential analyzer and the electron temperature probe to determine the relative importance of ion transport and electron precipitation. For most of the orbits, precipitation is inferred to be of low to moderate importance. Only for orbit 65, which was the first nightside orbit published by Taylor et al. [J. geophys. Res. 85, 7765 (1980)] and for the inbound portion of orbit 73 does the ionization structure appear to be greatly affected by electron precipitation.