Thinning of the near-earth (10∼15 RE) plasma sheet preceding the substorm expansion phase

The timing of the plasma-sheet thinning relative to the onset of the expansion phase of substorms is examined by the analysis of the OGO 5 electron (79 ± 23 keV) and proton (100～150 keV) data with the aid of simultaneous magnetic field observations. It is found that the timing of the thinning is significantly dependent on the distance. At $\text{x}{}^2\text{+ y}{}^2\text{? 15 R}{}_{\mathrm{E}}$ the thinning often starts before the onset, while at $\text{x}{}^2\text{+ y}{}^2\text{? 15 R}{}_{\mathrm{E}}$ it tends to occur after the onset, where x and y refer to solar magnetospheric coordinates. The thinning that precedes the expansion-phase onset has been found to reduce the thickness to ～1 R_E, and further thinning may occur in a spatially limited region. Hence it is conceivable that the formation of the neutral line characterizing the substorm expansion phase is the consequence of the thinning of the plasma sheet in the near-Earth region.     

Theory for the critical ionization velocity phenomenon

A theoretical consideration is given to the critical ionization velocity. A macroscopic theory is constructed, assuming that (1) the newly born ions have an unstable velocity distribution and (2) the resultant fluctuating electric field heats the electrons. For the critical ionization velocity to operate in the case of $\text{m}{}_{\mathrm{a}}\text{V}{}^2\text{2}\text{> eφ}{}_{\mathrm{<mtext>ion</mtext>}}$, the initial values of high speed ion (or newly ionized ion) density and suprathermal electron energy density should be above a certain threshold. This threshold depends upon the velocity diffusion co-efficient and the collisional loss rate of the newly ionized ion.     

On the simultaneous ionization-excitation of the OII(λ834 Å) resonance transition by electron impact on atomic oxygen

Laboratory cross-section data on the excitation of the OII(2s 2p⁴⁴P → 2s² 2p³⁴S; λ834 Å) resonance transition and on the production of O⁺ and O²⁺ ions by electron impact on atomic oxygen are used to show that the ratio $\text{σ(λ834}\text{A}\text{?}\text{)}\text{σ(O}{}^{\mathrm{+}}\text{+ O}{}^{\mathrm{2+}}\text{)}$ is nearly constant for incident electron energies > 50 eV. Under auroral conditions, the total electron-ion pair production rate from electron impact on O can be inferred from λ834 Å volume emission rate measurements using the result that $\text{η(O}{}^{\mathrm{+}}\text{+ O}{}^{\mathrm{2+}}\text{)}\text{\$}\text{?}\text{8.4η(λ834}\text{A}\text{?}\text{)}$. These findings, along with earlier work on the simultaneous ionization-excitation of the 1 Neg (0,0) band of N₂⁺ and the 1 Neg (1, 0) band of O⁺₂, allow the specific ionization rates for the principal atmospheric constituents (O⁺, O⁺₂, N⁺₂), for the multiply-ionized species (O²⁺, O²⁺₂, N²⁺₂), and for the dissociatively produced atomic ions to be inferred in aurora from remote satellite observations.     

Electron runaway in turbulent astrophysical plasmas

An astrophysical electron acceleration process is described which involves turbulent plasma effects: the acceleration mechanism will operate in ‘collision free’ magnetoactive astrophysical plasmas when ion-acoustic turbulence is generated by an electric field which acts parallel to the ambient magnetic lines of force. The role of ‘anomalous’ (ion-sound) resistivity is crucial in maintaining the parallel electric field. It is shown that, in spite of the turbulence, a small fraction of the electron population can accelerate freely, i.e. runaway, in the high parallel electric potential. The number density n^(B) of the runaway electron component is of order $\text{n}{}^{\mathrm{(B)}}\text{?}\text{n}\text{2}\text{(}\text{c}{}_{\mathrm{s}}\text{U}{}_{\mathrm{\parallel }}{}^{\mathrm{?}}\text{)}{}^2$, where n = background electron number density, c_s = ion-sound speed and U_∥^? = relative drift velocity between the electron and ion populations. The runaway mechanism and the number density n^(B) do not depend critically on the details of the non-linear saturation of the ion-sound instability.     

The O(1S) quantum yield from O2+ dissociative recombination

Recent laboratory studies show that the $\text{O}\text{(}{}^1\text{S}\text{)}$ quantum yield, $\text{f}\text{(}{}^1\text{S}\text{)}$, from O₂⁺ dissociative recombination varies considerably with the degree r of vibrational excitation. However, the suggestion that the high values for $\text{f}\text{(}{}^1\text{S}\text{)}$ deduced from airglow and auroral observations can be explained by invoking vibrational excitation, creates a number of problems. Firstly, the rapid vibrational deactivation of O₂⁺ ions by collisions with O atoms will keep r too low to account for the magnitude of $\text{f}\text{(}{}^1\text{S}\text{)}$; secondly, r varies considerably from one atmospheric source to another but its relative values (which should be reliable) do not co-vary with those of $\text{f}\text{(}{}^1\text{S}\text{)}$; thirdly, because r increases markedly above the peak of the $\text{X5577}\text{A}\text{?}$ dissociative recombination layer, the fits which theorists have obtained to the observed volume emission rate profiles would have to be regarded as fortuitious. It is tentatively suggested that $\text{f}\text{(}{}^1\text{S}\text{)}$ is higher in the airglow and aurora than in the laboratory plasma studied by Zipf (1980) because of the electron temperature dependence of the $\text{O}\text{(}{}^1\text{S}\text{)}$ specific recombination coefficient for O₂⁺(v^' ? 3) ions.The repulsive ${}^1\text{Σ}{}_{\mathrm{<mtext>u</mtext>}}\text{[}{}^1\text{D}\text{+}{}^1\text{s}\text{]}$ state of O₂ does not provide a suitable channel for the dissociative recombination. A possible alternative is the bound ${}^3\text{Π}{}_{\mathrm{<mtext>u</mtext>}}\text{[}{}^5\text{S}\text{+}{}^3\text{s}\text{]}$ state with predissociation to the repulsive ${}^3\text{Π}{}_{\mathrm{<mtext>u</mtext>}}\text{[}{}^3\text{P}\text{+}{}^1\text{s}\text{]}$ state.     

Non-adiabatic processes in a two-fluid plasma and energization of the plasma sheet plasma

The change of energy of a collisionless, two-fluid plasma consists of the adiabatic gain or loss of energy, which is due to the work done by the electromagnetic forces, and of the non-adiabatic change associated with the presence of the “rest” field $\text{E}{}^1\text{=}\text{E}\text{+ (}\text{1}\text{c}\text{)}\text{V}\text{×}\text{B}$. The non-adiabatic gain or loss of energy per unit ti may be expressed by the relation

$\text{Q=}\text{E}{}_{\mathrm{\parallel }}\text{·}\text{i}{}_{\mathrm{\parallel }}\text{+}\text{c}\text{eNB}{}^2\text{B·}\text{f}\text{?}\text{×f}$

where i is the density of conductive current, N the ion number-density, and f (f?) the sum of inertia and pressure divergence of ions (electrons). Symbols of parallelism refer to the direction of B.A special case of non-adiabatic energization of a slowly convecting plasma sheet plasma is discussed in some detail. Regardless of the value of V, the non-adiabatic energization may significantly exceed any conceivable energization associated with the electric field $\text{?(}\text{1}\text{c}\text{)}\text{V}\text{×}\text{B}$.     

Comment on the group velocity of surface waves on the plasma sheet boundary

In the recent estimation by Maltsev and Lyatsky (1984) of the group velocity of surface waves on the inner boundary of the plasma sheet, the effect of the curvature of the field lines of the ambient magnetic field of the Earth on the spectrum has been assessed. The authors have not accounted for the fact, however, that the group velocity of the compressional surface magnetohydrodynamic waves itself is nonzero transverse to the magnetic field—a characteristic which has been omitted in the spectrum of Chen and Hasegawa (1974), being used by Maltsev and Lyatsky.This characteristic of compressional surface MHD waves is inherent for the spectrum $\text{ω = (}\text{k}{}_6\text{k}\text{)V}{}_{\mathrm{A}}\text{(k}{}^2{}_6\text{+ 2k}{}^2{}_{\mathrm{\perp }}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, obtained by Nenovski (1978) in the cold plasma limit V_A ? V_S(V_A is Alfvén velocity, and V_S, sound velocity). A comment has been made on the restrictions, proceeding from the approximation, used by Maltsev and Lyatsky. The estimation of the velocities for movements of auroral riometer absorption bays have been reviewed.     

Characteristics of localized resonance coupling oscillations of the slow magnetosonic wave in a non-uniform plasma

We analyze linear resonance oscillations in a non-uniform one-fluid finite-β plasma, which is oversimplified to understand easily fundamental characteristics of the resonance oscillations. A linear resonance oscillation of localized slow magnetosonic mode $\text{[ω}{}^2{}_{\mathrm{s}}\text{=}\text{ω}{}^2{}_{\mathrm{A}}\text{(1 +}\text{V}{}^2{}_{\mathrm{A}}\text{V}{}^2{}_{\mathrm{s}}\text{)}\text{]}$, which has the diamagnetic property in a uniform plasma, is newly found to be excited in the radially non-uniform plasma. The localized slow resonance indicates a radially polarized compressional oscillation (δB_∥ ? δB_H ? δB_D). The sense of the Alfvénic polarizations in the H-D plane near the resonant point is a function of both the propagation in the azimuthal direction and the slope of wave amplitude in the radial direction, whereas the sense of the resonant slow magnetosonic polarizations changes in accordance only with the switch in the azimuthal propagation direction. Further multi-satellite studies are necessary to establish the resonant structures of the slow magnetosonic waves in the magnetosphere.     

Linewidth studies on the NI(4S−4P) resonance multiplet

We have measured the linewidths of the NI multiplets $\text{[2p}{}^2\text{3p}{}^4\text{D}{}^0\text{?2p}{}^2\text{3s}{}^4\text{P, λ8691}\text{A}\text{?}\text{; 2p}{}^2\text{3p}{}^4\text{P}{}^0\text{?2p}{}^2\text{3s}{}^4\text{P, λ8212}\text{A}\text{?}\text{; 2p}{}^2\text{3s}{}^4\text{P?2p}{}^3{}^4\text{S}{}^0\text{, λ1200}\text{A}\text{?}\text{]}$ produced in the dissociative excitation of N₂ by energy electrons. The infrared transitions excite the N(⁴P) resonance state by cascade and they account for > 50% of the total N(⁴P) cross section at 100 eV. Both the i.r. and v.u.v. lines are found to be highly Doppler broadened ( ～ 25 times the thermal Doppler line width). These results indicate that dissociative excitation of N₂ produces N (⁴P) atoms with sufficient kinetic energy so that the $\text{λ1200}\text{A}\text{?}$ resonance radiation $\text{[2p}{}^2\text{3s}{}^4\text{P ?2p}{}^3{}^4\text{S}{}^0\text{]}$ emitted by these excited atoms would be optically thin in the Earth's upper atmosphere. We also found that the line strength ratios for the resolved components of the $\text{λ1200}\text{A}\text{?}$ triplet excited by dissociative excitation differ from those predicted by the multiplicities of the states involved and used in current entrapment models; the intensity ratios also vary with the energy of the incident electron. These developments introduce new complications into the analysis of the terrestrial ultraviolet dayglow.     

Numerical studies of the transport of solar protons in interplanetary space

Numerical solutions of the Fokker-Planck equation governing the transport of solar protons are obtained using the Crank-Nicholson technique with the diffusion coefficient represented by K_r=K₀r^b where r is radial distance from the Sun and b can take on positive or negative values. As b ranges from +1 to ?3, the time to the observation of peak flux decreases by a factor of 5 for 1 MeV protons when $\text{V}\text{K}{}_0\text{=}\text{3 AU}{}^{\mathrm{b?1}}$ where V is the solar wind speed. The time to peak flux is found to be very insensitive to assumptions concerning the solar and outer scattering boundary conditions and the presence of exponential time decay in the flux does not depend on the existence of an outer boundary. At $\text{V}\text{K}{}_0\text{? 15 AU}{}^{\mathrm{b?1}}$, 1 MeV particles come from the Sun by an almost entirely convective process and suffer large adiabatic deceleration at b?0 but for b=+1, large Fermi acceleration is possible at all reasonable $\text{V}\text{K}{}_0$ values. Implications of this result for the calculation and measurement of particle diffusion coefficients is discussed. At $\text{b?0}$, the pure diffusion approximation to transport overestimates by a factor 2 or more the time to peak flux but as b becomes more negative, the additional effects of convection and energy loss become less important.     

Pitch angle dependence of the charge-exchange lifetime of ring current ions

Previous work has parameterized the pitch angle dependence of the charge-exchange lifetime τ of ring current ions in terms of γ, the power of the cosine of the mirror latitude λ_m of the particles, such that $\text{τ(λ}{}_{\mathrm{m}}\text{)}\text{τ(0) ≌}\text{cos}{}^{\mathrm{\gamma }}\text{λ}{}_{\mathrm{m}}$ at given L. Using the atomic hydrogen density model of Johnson and Fish, previous authors have suggested values of γ = 5 or 6. We here evaluate γ as a function of λ_m and L using the more recent Chamberlain density models, and show that γ = 3?4 is more appropriate over most of the pitch angle and L range. Consequently, ion distributions in the ring current decay phase are expected to become rather less anisotropic in pitch angle due to chargeexchange than previously believed. We have also investigated the use of several other simple approximate analytic forms for $\text{τ(λ}{}_{\mathrm{m}}\text{)}\text{τ(0)}$, one of which gives far better agreement with the numerical results than the cos^γ λ_m, variation, and should hence be used in future studies.     

Thermospheric neutral composition from auroral ion ratios

For nighttime auroras, we find that positive ion ratios are only a function of the neutral atmospheric composition and of the pertinent ionic processes if the ions are depleted mainly by ion-molecule reactions. Ionic ratios calculated for $\text{[}\text{N}{}^{\mathrm{+}}\text{]}\text{[}\text{N}{}_2{}^{\mathrm{+}}\text{]}$ using the 1976 U.S. Standard Atmosphere and laboratory rate coefficients (with one exception) rise smoothly with altitude: 0.1 (120 km), 0.3 (140 km), 0.6 (160 km), 1.0 (180 km) and 1.5 (200 km). These values compare favorably with experimental ratios from three different auroral experiments. The exception refers to our use of a larger rate coefficient for N₂⁺ + O → NO⁺ + N than found in the laboratory. We also determine an $\text{[}\text{N}{}_2{}^{\mathrm{+}}\text{]}\text{[}\text{O}{}^{\mathrm{+}}\text{]}$ ratio with altitude: 0.36 (120 km), 0.078 (140 km), 0.030 (160 km), 0.014 (180 km) and 0.0075 (200 km). These values compare favorably with results from the same three auroral experiments. However, the match with a fourth auroral experiment is poor. Except for this last case, we conclude that the neutral composition at auroral latitudes in late winter is similar to the U.S. Standard for the altitudes examined.     

Non-linear motion of a weak inhomogeneity in a magneto-active plasma

The equations of motion of a magnetic field aligned column of ionization embedded in an unbounded weakly ionized plasma immersed in an external electric field E are solved to second order terms in the strength of the initial perturbation, which is assumed to be weak. It is shown that when the Hall currents induced by the internal electrostatic field are taken into account, the centre of mass of the column will, in general, move along a curved path before acquiring the familiar $\text{E ×}\text{B}\text{B}{}^2$ drift velocity.     

Temperature variation of the plasma sheet during substorms

The temperature and density of the plasma in the Earth's distant plasma sheet at the downstream distances of about 20–25 R_e are examined during a high geomagnetic disturbance period. It is shown that the plasma sheet cools when magnetospheric substorm expansion is indicated by the AE index. During cooling, the plasma sheet temperature, T, and the number density, N, are related by $\text{T ∝ N}{}^{\mathrm{<mtext>2</mtext><mtext>3</mtext>}}$ (adiabatic process) in some instances, while by T ∝ N^?1 (isobaric process) in other cases. The total plasma and magnetic pressure decreases when $\text{T ∝ N}{}^{\mathrm{<mtext>2</mtext><mtext>3</mtext>}}$ and increases when T ∝ N^?1. Observation also indicates that the dawn-dusk component of plasma flow is frequently large and comparable to the sunward-tailward flow component near the central plasma sheet during substorms.     

Perturbed particle disks

The Boltzmann moment equations are solved to determine the velocity ellipsoid in a particle disk near an isolated satellite resonance. In a coordinate frame which rotates with the pattern speed of the perturbation potential, the solutions are stationary functions of the azimuthal angle. From the velocity ellipsoid we obtain the stress tensor due to particle collisions and consequently, the viscous angular momentum flux. We show that the magnitude of the rate of deformation tensor in a perturbed particle disk is bounded from above by $\text{KΩ(1 + τ}{}^2\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{,}\text{where}\text{Ω}$ is the orbital angular velocity, τ is the optical depth, and K is a dimensionless constant of order unity. It is also found that in sufficiently perturbed regions there are ranges of azimuthal angle within which the radial component of the angular momentum flux is negative. It is even possible for the angular momentum luminosity, the radial flux integrated over azimuth, to be negative. These results are important for understanding sharp edges and the decay of density waves in planetary rings. They are also relevant to the damping of differential precession and eccentricity in narrow ringlets.     

Total current to cylindrical collectors in collisionless plasma flow

A theory is presented for charged-particle collection by a cylindrical conducting object, such as a spacecraft or an electrostatic probe, which is moving transversely through a collisionless plasma, such as those in the upper atmosphere and space. The calculation is approximate, using symmetric potential profiles which are exact for the infinite-cylinder stationary case. Theoretical current predictions are presented for ratios of collector potential to electron thermal energy eφ_c/kT_e from 0 to ?25, for ion-to-electron temperature ratios T_i/T_c = 1 and 0.5, ratio of collector radius to electron Debye length r_c/λ_D from 0 to 100, and ratio of flow speed to ion thermal speed $\text{S}{}_{\mathrm{i}}\text{= U/(2kT}{}_{\mathrm{i}}\text{/m}{}_{\mathrm{i}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}$ from 0 to 10. Comparisons with existing exact calculations by other authors show that none of these fulfil all of the requirements for nontrivial comparison. Appropriate parameter ranges for future exact calculations are thereby suggested. These are as follows: (a) r_c/λ_D should be large enough that the collector not be in or near orbit-limited conditions; (b) the ratio $\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,\; i}}\text{¦}$ of ion directed energy to potential energy change in the sheath, should be close to unity or if

$\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,i}}\text{¦? 1,}\text{then}\text{S}{}_{\mathrm{i}}\text{? 1}$

.     

Associative ionization of N(2D) and O

Recent laboratory measurements have shown that N(²P) atoms, and thus probably hot N(²D) atoms, will recombine with atomic oxygen via an associative ionization process at the gas kinetic rate. While the reaction is endothermic, it has been suggested that this has interesting implications for the upper atmosphere in that N(²D) atoms in the tail of the velocity distribution could provide an additional source of NO⁺ through the reverse of the dissociative recombination reaction

$\text{NO}{}^{\mathrm{+}}\text{+ e ?}\text{N}\text{(}{}^2\text{D}\text{) +}\text{O}$

. It has also been suggested that this process might account for the difference between a laboratory determination of the rate coefficient and that determined from the Atmospheric Explorer Satellite data. In this paper we investigate further the likelihood of the associative ionization of N(²D) and O playing a significant role in the normal ionosphere, in the light of several recent relevant studies. We conclude that the associative ionization process is not an important factor and that a more probable cause for disagreements in the various determinations of the recombination coefficient, is the difference in excited states of the ions in the various experiments.     

Direct measurement of conjugate photoelectrons and predawn 630 nm airglow

A sounding rocket was flown during the predawn on 17 January, 1976 from Uchinoura, Japan, to measure directly the behaviour of the conjugate photoelectrons at magnetically low latitudes. On board the rocket were an electron energy analyzer, 630 nm airglow photometer, and plasma probes to measure electron density and temperature. The incoming flux of the photoelectrons was measured in the altitude range between 210 and 340 km. The differential flux at the top of the atmosphere was determined to be $\text{F = (1.3 ± 0.4) × 10}{}^{11}\text{exp}\text{[?}\text{E(}\text{eV}\text{)}\text{12}\text{]}$ electron · m^?2 · sr^?1 · s^?1 in the energy range 10 ? E ? 50 eV. The emission rate of the 630 nm airglow was observed in the altitude range between 90 and 360 km. The apparent emission rate observed at 80 km was 32 ± 5 R. From a theoretical calculation of the optical excitation rate using the observed electron flux data along with a model distribution of atomic oxygen, it was estimated that more than 65% of the emission could be produced by direct impact of the photoelectrons with atomic oxygen in the thermosphere between 200 and 360 km. Using the observed electron density and the model distribution of oxygen molecules the residual of the emission was ascribed to the excitation of O(¹D) through dissociative recombination, $\text{O}{}_2{}^{\mathrm{+}}\text{+}\text{e}\text{→}\text{O}{}^1\text{+}\text{O}{}^7$. The direct collisional excitation by ambient electrons is estimated to be negligibly small at the level of observed electron temperature.     

Emission line shapes produced by dissociative excitation of atmospheric gases

Special line shapes are derived fro the λ 1356 Å (${}^5\text{S}{}^0\text{-}{}^3\text{P}$) transition of atomic oxygen from metastable (${}^5\text{S}{}^0\text{-}{}^3\text{P}$) time-of-flight spectra produced by electron impact dissociative excitation of O₂, CO₂, CO, and NO, and they are compared with the broadened λ 1304 A resonance line shapes deduced by Poland and Lawrence (1973) from atomic oxygen absorption studies. The non-thermal line shapes for both airglow emission features are shown to have an effective width comparable to a 60,000 K thermal doppler line shape for an electron impact energy of 100eV. The variation of the effective line width with electron-impact energy from threshold to 300 eV is given. Since the effective line width of the resonance radiation produced by dissociative excitation is very large compared with the doppler absorption widths of the ambient O atoms at normal exospheric temperatures, the anomalously broadened resonance lines will propagate through a planetary atmosphere as though they were optically thin. Thus, electron-impact dissociation of CO and CO₂ will contribute to the observed optically thin component of the λ 1304 Å emission in the upper atmospheres of Venus and Mars. However, the process cannot account for more than 10% of the observed optically thin emission because of the small magnitude of the excitation cross-section and the comparatively high-energy threshold for the process. The possibility that the source of the kinetically energetic $\text{O}\text{(}{}^3\text{S}$) atoms is the dissociative recombination of vibrationally excited CO₂⁺ ions is discussed.