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.     

Generalized limits to the number of light particle degrees of freedom from big bang nucleosynthesis

We compute the big bang nucleosynthesis limit on the number of light neutrino degrees of freedom in a model-independent likelihood analysis based on the abundances of ⁴He and ⁷Li. We use the two-dimensional likelihood functions to simultaneously constrain the baryon-to-photon ratio and the number of light neutrinos for a range of ⁴He abundances Y_p = 0.225–0.250, as well as a range in primordial ⁷Li abundances from (1.6 to 4.1) ×10⁻¹⁰. For (⁷Li/H)_p = 1.6 × 10⁻¹⁰, as can be inferred from the ⁷Li data from Population II halo stars, the upper limit to N_ν based on the current best estimate of the primordial ⁴He abundance of Y_p = 0.238 is N_ν < 4.3 and varies from N_ν < 3.3 (at 95% C.L.) when Y_p = 0.225 to N_ν < 5.3 when Y_p = 0.250. If ⁷Li is depleted in these stars the upper limit to N_ν is relaxed. Taking (⁷Li/H)_p = 4.1 × 10⁻¹⁰, the limit varies from N_ν < 3.9 when Y_p = 0.225 to N_ν 6 when Y_p = 0.250. We also consider the consequences on the upper limit to N_ν if recent observations of deuterium in high-redshift quasar absorption-line systems are confirmed.     

Spectrum of the relic neutrino background from past supernovae and cosmological models

It is greatly expected that the relic neutrino background from past supernovae will be detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, Λ, on the spectrum, since some recent cosmological observations strongly suggest the existence of Λ. We find following results. (1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimates (6–10 MeV). (2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae are dominant, is about 25h₅₀²(R_SN/0.1 yr⁻¹)(n_Gh₅₀⁻³/0.02 Mpc⁻³) events per year, where R_SN is the supernova rate in a galaxy, n_G is the number density of galaxies, and h₅₀ = H₀/(50 km/s Mpc), where H₀ is the Hubble constant. (3) The event rate is almost insensitive to Λ. The flux increases in the low energy side (< 10 MeV) with increasing Λ, but decreases in the high energy side (> 10 MeV) in models in which the integrated number of supernovae in one galaxy is fixed.     

A theoretical study of the time-dependent behaviour of O in the mid-latitude plasmasphere

The behaviour of O²⁺ at L = 3 in the plasmasphere is studied. Starting with a low O²⁺ flux-tube content to characterize post-magnetic-storm conditions the time-dependent equations of continuity and momentum for O²⁺ are solved to give densities and fluxes for a period of several days using both sunspotmaximum and sunspot-minimum parameters. Our results show large amounts of O²⁺ near the equator at sunspot maximum but relatively little at sunspot minimum, and emphasize the key role of the collisional process between O²⁺ and O⁺. It is the combined effects of O²⁺---O ⁺ collisions and thermal diffusion that lead to the large O²⁺ densities near the equator at sunspot maximum. Both of these mechanisms have less influence at sunspot minimum. At sunspot maximum the O⁺ layer acts as a collisional barrier below the O²⁺ production region preventing O²⁺ from sinking towards regions of high recombination rate. In this production region the effects of thermal diffusion are small and upward flow of O²⁺ results from the action of the O²⁺ pressure gradient and the polarization electric field. When the upward flowing O²⁺ reaches regions in which thermal diffusion has a strong influence it is accelerated to even higher altitudes. The O ⁺ barrier is so effective that the diurnal variation of the O⁺ layer is reflected in the diurnal variation of O²⁺ near the equator at sunspot maximum. Our sunspot maximum results also indicate that certain types of temperature profiles are more likely to enhance equatorial O²⁺ densities. The existence of large temperature gradients below 1000 km altitude does not help the flow of O²⁺ towards the equator. The associated changes in the O⁺ layer lead to more O²⁺-O ⁺collisions and a smaller O²⁺ thermal-diffusion coefficient, the latter being sensitive to the ratio n(H⁺)/n(O⁺).     

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.     

A study of F2 region night-time vertical ionization fluxes at millstone hill

Vertical fluxes of ionization in the F2 region have been measured by the incoherent scatter technique over Millstone Hill in 1969. The results obtained near midnight for the region above h_maxF2 have been examined to determine whether there is a significant flux of ionization from the magnetosphere to the ionosphere that serves to maintain the F-layer. It is found that H⁺ ions are a minor constituent over the altitude range in which useful measurements can be made, so that any conclusion must rest upon properly interpreting the observed O⁺ fluxes. By selecting periods when the layer did not appear to be decaying rapidly it was hoped to find cases where the O⁺ flux did not vary with altitude in the range 500 h 800 km (i.e. where losses are unimportant), since this would imply that the flux is of magnetospheric origin.

While three cases exhibited this behaviour, the majority exhibited a decrease in the O⁺ flux with height, indicating that the layer was descending. Attempts to correct for this were made, and the average flux from the magnetosphere was estimated as 3 × 10⁷ el/cm²/sec. This is in fair agreement with other recent estimates, and implies that at this latitude the ionosphere is not maintained solely by the magnetospheric flux. Moreover, large increases in flux that could give rise to nocturnal increases in the total content of the layer do not appear to have been seen.     

Satellite beacon measurements of protonospheric fluxes

The ionospheric and protonospheric regions of the plasmasphere, which are dominated by the O⁺ and H⁺ ionic species, respectively, interact by means of proton fluxes within tubes of magnetic force. The present study is concerned with the determination of these fluxes by the beacon satellite technique as used in the ATS-6 experiment in relation to three observing sites: Boulder, Colorado; Lancaster, U.K.; and Fairbanks, Alaska. From plasmasphere models based on solutions of the time dependent O⁺ and H⁺ momentum and continuity equations, it is shown that the time differential of the “residual content” as measured at Lancaster, provides a good estimate of the protonospheric flux at 4000km altitude in the L = 1.8 magnetic shell under quite geomagnetic conditions. The effect of the neutral thermospheric wind on the protonospheric flux is also investigated. Fluxes determined by the beacon technique for the period from September 1975 to July 1976 are shown, and these are compared with typical results derived from other techniques.     

Interpreting vertical plasma drift in the mid-lattitude ionosphere using ionosonde measurements

Measurements of the density at the F₂ peak (N_mF₂) were obtained by the Boulder, Colorado, ionosonde as part of the SUNDIAL-86 campaign. The measurements were made during a period of low to moderate geomagnetic activity following a “disturbed” day. These measurements were then used to estimate the height of the F₂ peak (h_mF₂). A three-dimensional time-dependent model of Earth's ionosphere was used to calculate N_mF₂ and h_mF₂ using the vertical plasma drift as a free parameter. Since the plasmasphere-ionosphere exchange flux can remain upward during the night for these conditions, different feasible flux scenarios were inputed to the ionospheric model. These different flux scenarios had a large effect on the “induced” vertical plasma drifts required to match the measurements (i.e. at times greater than a factor of 2 in speed or a difference in direction). Futhermore, uncertainty in the O⁺---O collision frequency changes the required vertical plasma drift at night. Despite knowledge of h_mF₂, interpretation of the vertical plasma drifts as meridional neutral winds is compromised by a lack of knowledge of the plasmasphere-ionosphere exchange flux following disturbed days.     

The spectrum of electron content fluctuations in the ionosphere

Continuous records of the electron content of the ionosphere, from 1965 to 1970, are used to obtain power spectra covering periods from 30 sec to 2 yr at latitudes of 34°S and 42°S. At periods up to 5 min, amplitudes were less than 0.2 per cent of the total electron content. Variations produced by gravity waves were very common in the range 20–80 min, with no preferred periods. The r.m.s. amplitude per octave A₀ was about 10¹⁵ electrons/m², or 0.6 per cent of the mean electron content. The amplitude increased during the day, particularly in winter when periodic components predominated. The cut-off at about 17 min was sharply defined, giving a mean scale height for the neutral atmosphere (at 300 km) of about 43 km in summer, 47 km on winter days and 42 km on winter nights.

From 12 hr to 1 month A₀ was about 12 per cent of the mean electron content in both summer and winter at 34°S, and 10 per cent at 42°S. The 24 hr and 27 day peaks were largest just before sunspot maximum, and almost disappeared near sunspot minimum. Variations between 1 and 27 days reflect the random occurrence of ionospheric storms and show no consistent peaks. Day to day and night to night variations were both about 10 per cent of the background content for periods from 2 days to 2 yr, apart from a slight decrease between 1 and 6 months.     

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.     

A measurement of the O2 ultra-violet nightglow emission

New measurements of the Herzberg I emission height profile in the night airglow are reported and indicate a peak emission height near 96 km in agreement with previous measurements. Using an atomic oxygen concentration profile determined from the oxygen green line profile measured on the same rocket it is concluded that the O₂(A³Σ_u⁺) state is not excited in the direct three body recombination of atomic oxygen. It is suggested that the excitation mechanism is a two step process, similar to the Barth mechanism for the atomic oxygen green lineand that the excited intermediate state is C³Δ_u.     

Kp index correlations with solar-wind parameters during the first and second stages of a recurrent geomagnetic storm

The dependence of geomagnetic activity during a recurrent magnetic storm on the solar-wind magnetic field and plasma parameters has been studied. According to variations of solar-wind magnetic field strength B, a recurrent magnetic storm is divided into two stages: the first proceeding during the peak of B, and the second proceeding after the return of B to quiet level. The K_p index vs solar-wind parameters scattering diagrams for stages I and II differ significantly. In particular, the random scattering for stage I is much larger than for stage II. It was found that for stage I the K_p index correlates with B, with the sign and value of northsouth field component B_z and with the magnitude ΔB of field fluctuations, the situation being similar to that during sporadic magnetic storms, though the scale of the event is smaller. For stage II, the K_p index does not correlate with B, but strongly correlates with ΔB and weaker—with B_z. So geomagnetic activity at stage II is supported mainly by solar-wind magnetic field fluctuations. The dependence of the K_p index on plasma parameters (concentration of protons n, bulk velocity v and temperature T) is weak for both stages.     

The relation between low-latitude neutral density variations near 400 km and magnetic activity indices

Neutral density data were obtained near 400km (1600 LT) from a microphone density gauge on OGO-6 from 0°G to 40°N magnetic latitude for 25 September–3 October 1969. Several geomagnetic storms occurred during this period (a_p varied from 0 to 207). Least-squares fits were made to data points on density-a_p and density-D_st scatter diagrams, where the density values selected were delayed in time behind a_p and D_st. An equation representing the least-squares fit was computed for each delay time. The equation of best fit (and the corresponding time delay between the density and the magnetic index which resulted in this best fit) was found by choosing the equation that gave the minimum standard error. For example, the best fit at 10°N geomagnetic latitude occurred for a_p at t — 3 hr, where t is the time of the density values. The implications of the time differences associated with the best fits at various latitudes and longitudes are discussed with regard to the time delays involved in geomagnetic heating of the neutral upper atmosphere.

A low-latitude density bulge has been found between 0°N and 40°N whose magnitude varies with a_p. DeVries (1972b) has independently discovered this daytime phenomenon. If the bulge is a semi-permanent feature near the equinoxes because of the enhanced geomagnetic activity, this may help explain the semi-annual effect in density, which was uncovered first in the drag data from low inclination satellites.     

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.     

Auroral N2(c′4Σu → aΠg emission

Inspection of recent spectra presented by Sivjee (1983) show evidence of the 0–4 and 0–5 bands of the N₂(c′₄¹Σ_u⁺ → a¹Π_g) Gaydon-Herman system. In conjunction with earlier spectra, it is now possible that this band system is a significant auroral component, with an intensity approx. 7% that of the N2 2P system. The absence in aurorae of the potentially far stronger N₂(c′₄¹Σ_u⁺ → X¹Π_g) system is discussed. It is that the O₂(A³Σ_u⁺ → X³Σ_g⁻) band system is indiscernible in Sivjee's auroral spectra, under conditio the foreground nightglow is expected to be clearly visible. On the other hand, at least one relatively strong O₂(A′³Δ_u → a¹Δ_g) band appears to be present in these spectra.     

Convection and field-aligned currents, related to polar cap arcs, during strongly northward IMF (11 January 1983)

Electric and magnetic fields and auroral emissions have been measured by the Intercosmos-Bulgaria-1300 satellite on 10–11 January 1983. The measured distributions of the plasma drift velocity show that viscous convection is diminished in the evening sector under IMF B_y < 0 and in the morning sector if IMF B_y > 0. A number of sun-aligned polar cap arcs were observed at the beginning of the period of strongly northward IMF and after a few hours a θ-aurora appeared. The intensity of ionized oxygen emission [O⁺(²P), 7320 Å] increased significantly reaching up to several kilo-Rayleighs in the polar cap arc. A complicated pattern of convection and field-aligned currents existed in the nightside polar cap which differed from the four-cell model of convection and NBZ field-aligned current system. This pattern was observed during 12 h and could be interpreted as six large scale field-aligned current sheets and three convective vortices inside the polar cap. Sun-aligned polar cap arcs may be located in regions both of sunward and anti-sunward convection. Structures of smaller spatial scale correspond to the boundaries of hot plasma regions related to polar cap arcs. Obviously these structures are due to S-shaped distributions of electric potential. Parallel electric fields in these S-structures provide electron acceleration up to 1 keV at the boundaries of polar cap arcs. The pairs of field-aligned currents correspond to those S-structures: a downward current at the external side of the boundary and an upward current at the internal side of it.     

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.     

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.     

Small-scale auroral arc distortions

Small-scale spatially periodic distortions of auroral forms have been studied utilizing low-light level television observations made at various locations in the Northern and Southern Hemispheres. The most commonly observed features were folds and vortex-like curl formations. The curls, identified here with the Kelvin-Helmholtz instability due to fluid shear, invariably had a counterclockwise rotational shape and motion when viewed in a direction anti-parallel to the Earth's magnetic field. The typical measured wavelength (5 km) and measured growth rate (4.2 sec⁻¹) were used to evaluate the Kelvin-Helmholtz dispersion relation for the apparent shear ω_s = ∂ ν_x/ ∂y (28 sec⁻¹). The apparent horizontal velocities of both folds (0–5 km/sec) and curls (0–22 km/sec) were invariably observed to be counterclockwise with respect to the multiple arc centre when viewed antiparallel to B. Consistent agreement between rotational shape and rotational motion suggests that the apparent growth rate and the apparent horizontal velocities closely approximate the actual values. If the shear results from E×B drifts in a space charge field, the calculated value for ω_s, implies an unneutralized electron density 0–1 cm⁻³ and a ΔE across the arc element 500mV/m. The velocity measurements indicate that the ΔE values for individual elements can combine to produce transient electric fields at the edges of multiple arcs as high as 1000 mV/m.     

Twilight effects of solar ionizing radiation

The absorption of solar ionizing radiation during twilight is investigated. Ion production rates are obtained as a function of altitude and twilight intensities and altitude profiles of emissions arising from the fluorescence of solar ionizing radiation are calculated for various solar depression angles. For an atmosphere with an exospheric temperature of 750°K, the predicted overhead intensity from fluorescence of the O⁺(²P — ²D) lines at 7319–7330 diminishes from 175 R at dusk to 10 R at a solar depression angle of 10°. The predicted overhead intensities from fluorescence of the N₂⁺ Meinel and first negative systems are respectively about 175 R and 20 R at dusk diminishing to respectively 1.5 R and 0.1 R at a solar depression angle of 10°.

It is suggested that a charge transfer reaction of O⁺²D in N₂ is a significant source of N₂⁺ ions. This reaction offers a possible explanation for the high apparent rotational temperatures in the first negative system observed by Broadfoot and Hunten. Other excitation and ionization mechanisms are briefly discussed.