Modulation of the auroral proton spectrum in the upper atmosphere

An energetic auroral proton entering the atmosphere will alternate between being a proton and a neutral hydrogen atom by charge-exchange collisions with atmospheric constituents. This study uses a simple procedure to evaluate the energy degradation of the penetrating protons/hydrogen atoms by using semi-empirical range relations in air, and derives the particle energy variation as a function of altitude, starting from proton spectra observed from rockets above the main collision region. The main assumptions are that the geomagnetic field is homogeneous and vertical and that the pitch angle of the proton/hydrogen atom is preserved in collisions with atmospheric constituents before being thermalized. The calculations show that the incoming particle flux first loses the low energy particles at the highest pitch angles, even if the beam itself widens as it penetrates the atmosphere. The largest energy loss for particles with initial energy between 10 and 1000 keV occurs in the height interval between 100 and 125 km.     

Calculation of auroral Bahner volume emission height profiles in the upper atmosphere

Energetic protons entering the atmosphere will either travel as auroral protons or as neutral hydrogen atoms due to charge-exchange and excitation interactions with atmospheric constituents. Our objective is to develop a simple procedure to evaluate the Balmer excitation rates of H_α and H_β, and produce the corresponding volume emission rates vs height, using semi-empirical range relations in air, starting from proton spectra observed from rockets above the main collision region as measured by Reasoneret al. [(1968) J. geophys. Res.73, 4185] and Søbraaset al. [(1974) J. geophys. Res.79, 1851]. The main assumptions are that the geomagnetic field is parallel and vertical, and that the pitch angle of the proton/hydrogen atom is preserved in collisions with atmospheric constituents before being thermalized. Calculations show that the largest energy losses occur in the height interval between 100 and 125 km, and the corresponding volume emission rate vs height profiles have maximum values in this height interval. The calculted volume emission rate height profile of H_β compares favorably with that measured with a rocket-borne photometer.     

The solar eclipse of 26 February 1979: analysis and modelling of primary pair production,electron density and ionic composition

The solar eclipse of 26 February 1979 was observed from Red Lake, Canada, (52 °N, 91 °W) where totality occurred at about 1053 local time. Several research groups and government agencies participated in an extensive ground- and rocket-based observational program directed at the middle atmosphere. At the time of the eclipse, an extensive geomagnetic storm was in progress and the middle atmosphere was undergoing temperature and circulation changes associated with a stratospheric warming. Concurrent observations of atmospheric constituents, solar radiation, electron flux and other middle atmosphere parameters were obtained as inputs for a D-region predictive chemical computer code, DAIRCHEM, tailored to eclipse conditions. Ion pair production rates were computed by an E-region infrared radiance model and were used as necessary source function input values for DAIRCHEM computations. The computations yielded predictions of electron and total positive ion densities about totality. The positive ion measurements of a supersonic Gerdien condenser and a subsonic blunt probe during the eclipse were in agreement with the model computations and provided normalizing summations of total positive ions for the interpretation of mass spectrometer measurements. The chemical computer code identified principal routes for increase and removal of key species such as O₂⁺, NO⁺, hydrated clusters and negative ions. The dominant precursor ion for pair production hydrates was O₂⁺ and the chemistry was characteristic of the disturbed D-region.     

Stratospheric positive ion composition measurements,ion abundances and related trace gas detection

Positive ion spectra obtained from measurements with a balloon-borne mass spectrometer during three balloon flights are critically investigated and compared with other data. Ion abundances for proton hydrates [H⁺(H₂O)_n ions] at different stratospheric temperatures are compared, as well as the abundances of non proton hydrates H⁺X_itl(H₂O)_m, X being most likely CH₃CN.The detection of trace gases from ion composition measurements is discussed and an upper limit for the number densities of minor constituents such as NH₃ and CH₃OH is estimated at 35 km. Although sodium compounds cannot be responsible for the major positive ions, a closer investigation of high resolution daytime spectra suggests a small contribution of sodium in daytime ion chemistry.     

The effects of differential ion flows on EISCAT observations in the auroral ionosphere

During geomagnetic storms different partial pressure gradients in the auroral ionosphere may result in H⁺, He⁺, O⁺ and molecular ions drifting with different velocities along the Earth's magnetic field line. For relative drift velocities ⪡ 400 m s⁻¹ it is shown that differential ion flows may be identified by two signatures in the autocorrelation function (ACF) measured by EISCAT. For larger relative drifts numerical simulations show that these signatures still exist and may result in an asymmetry in the incoherent scatter spectrum for O⁺ and molecular ions. It is demonstrated that UHF data can be reliably analysed for k²λ_D² ≲ 1, but at high altitudes, where O⁺–H⁺ flows are expected, UHF observations will be restricted by large Debye lengths (k²λ_D² > 1). Examples of ACFs based on polar wind theory are presented and discussed for the VHF system and finally it is shown that large ion temperature ratios (T_i(H⁺) >T_i(O⁺)) can significantly affect the velocity determination.     

Can the high latitude ionosphere support large field-aligned ion drifts?

Ion velocities perpendicular and parallel to the geomagnetic field have recently been deduced by Smith et al. from bistatic measurements at 71° geomagnetic latitude in the afternoon sector. The results of this experiment include large (>400 m s⁻¹) downward ion velocities parallel to the magnetic field that persist for hours, small (100 m s⁻¹) ion velocities perpendicular to the magnetic field and electron density profiles with extremely narrow full-width at half-maximum. The explanation of these results was that the ionospheric flux tubes observed were near the terminator, and thus, sunlit at the top and in darkness at the bottom. The difference in production between the top and bottom of the flux tube creates an excess of ions at the top, which rapidly diffuse downwards. A three-dimensional, time-dependent model of the ionosphere has been used to test this explanation. Numerical experiments were performed to determine upper limits for the downward ion velocity. Assuming reasonable vertically-induced ion drifts due to either neutral winds or plasma convection, these upper limits were substantially smaller than the measurements. The location of the terminator was found to contribute a maximum of about 60 m s⁻¹ to the vertical ion velocity due to diffusion in a partially illuminated flux tube. In an attempt to explain the narrow density profiles without invoking an additional ionization source, the downward force in the model was arbitrarily increased, as would occur due to parallel electric fields in the ionosphere. Since the interpretation of these measurements as large field-aligned flows seems untenable by a model thought to be consistent with the currently accepted physics of the atmosphere, an alternate hypothesis is presented. If the common volume measurement is made in a region of O⁺ precipitation, then the line profile would not be Doppler shifted when viewed off-zenith. Therefore, the field-aligned velocities would be small, and the narrow width of the profiles would be due to enhanced electron densities in an O⁺ arc.     

The trouble with thermospheric vertical winds: geomagnetic,seasonal and solar cycle dependence at high latitudes

The vertical wind component is frequently used to determine the zero-velocity baseline for measurements of thermospheric winds by Fabry-Perot and other interferometers. For many of the upper atmospheric emission lines from which Doppler shifts are determined, for example for the OI 630 nm emission, available laboratory sources are not convenient for long-term use at remote automatic observatories. Therefore, the assumption that the long-term average vertical wind is zero is frequently used to create a baseline from which the Doppler shifts corresponding with the line-of-sight wind from other observing directions can then be calculated. A data base consisting of 1242 nights of thermospheric wind measurements from Kiruna (68°N, 20°E), a high-latitude site, has been analysed. There are many interesting short-term fluctuations of the vertical wind which will be discussed in future papers. However, the mean vertical wind at Kiruna also has a systematic variation dependent on geomagnetic activity, season and solar cycle. This means that the assumption that the average value of the vertical wind is zero over the observing period cannot be used in isolation to determine the instrument reference or baseline. Despite this note of caution, even within the auroral oval, the assumption of a zero mean vertical wind can be used to derive a baseline which is probably valid within 5 ms⁻¹ during periods of quiet geomagnetic activity (K_p < 2), near winter solstice. During other seasons, and during periods of elevated geomagnetic activity, a systematic error in excess of 10 ms⁻¹ may occur.     

Positive and negative ion composition measurements in the D- and E-regions during the 26 February 1979 solar eclipse

Two rockets bearing quadrupole mass spectrometers capable of measuring both positive and negative ion composition were launched from Red Lake, Canada, during the solar eclipse. Both instruments had liquid helium cryopumps and shock-attaching conical samplers. The payloads also contained two Gerdien condensers to measure total positive and negative ion concentrations and ion mobilities. Attitude control systems aligned the payloads with the velocity vector throughout ascent and descent. The first rocket was launched so that the D-region was in darkness 35 ± 8 s on the upleg and about 150 ± 15 s on the downleg for the study of ionospheric decay processes. The second rocket was fired after totality into 75% solar illumination for the study of ionospheric recovery. The positive ion composition above 105 km exhibited a strongly increasing NO⁺/O₂⁺ ratio with time after second contact due to O₂⁺ charge transfer with NO and a sharply diminished ionization rate. However, in both nights, the ionization below 105 km was created mainly by energetic particle deposition as exemplified by the increased ion concentrations and the composition signatures of a particle event: asignificant enhancement of O₂⁺ below 105 km and large amounts of H₅O₂⁺ ions in the D-region which result from the O₂⁺ clustering scheme. H₅O₂ was the major ion in the upper D-region while H₇O⁺₃, H₉O₄⁺ and H₅O₂⁺ were dominant ions at lower altitudes. Numerous minor species were also detected. The negative ion distributions in both flights exhibited a distinct shelf at 83 ± 2 km, decreasing by more than an order of magnitude by 90 km and with minima near 75 km. In the 75–90 km range, a significant percentage of the negative ions had masses exceeding 160 a.m.u. Comparisons are made with prior negative ion measurements during similar daytime auroral zone absorption (AZA) events. Two striking characteristics of the precipitating particles were apparent from these and past observations in daytime AZA events: there is a near absence of low energy electrons capable of ionizing above about 105 km and there is'a significant spatial and/or temporal variability in the electron flux. This paper is devoted principally to a presentation of the ion composition measurements and associated uncertainties.     

Topside and intelhemispheric ion flows in the mid-latitude plasmasphere

A modelling study has been carried out of field-aligned ion flows in the topside ionospheres of conjugate hemispheres under solstice conditions at mid to low latitudes. In the model calculations coupled time-dependent O⁺, H⁺ and electron continuity, momentum and heat balance equations are solved along dipole magnetic field lines at L = 1.5 and 3.0 Sunspot medium and sunspot minimum atmospheric conditions are considered.It has been found that thermal coupling between conjugate hemispheres gives rise to strong flows of O⁺ in the topside ionosphere of the summer hemisphere that are directed upwards at conjugate sunrise and directed downwards at conjugate sunset. At conjugate sunrise in the winter hemisphere there is a small upward-directed signature in the O⁺ field-aligned flux; there is no observable signature in the O⁺ field-aligned flux in the winter hemisphere at conjugate sunset. There are strong upward and downward flows of O⁺ at local sunrise and local sunset, respectively, in both the summer and winter hemispheres.At both L = 1.5 and 3.0 the 24 h time-integrated interhemispheric H⁺ flux is in the direction summer hemisphere to winter hemisphere. At L = 1.5 its magnitude is in good agreement with the magnitude of the 24 h time-integrated plasma (O⁺ + H⁺) field-aligned flux at 1000 km altitude; there are no such agreements at L = 3.0.A study of the roles played by the individual terms of the O⁺ momentum equation has demonstrated the complex structure of momentum balance. Certain of the terms may be orders of magnitude greater than the combined total of the individual terms, i.e. the O⁺ field-aligned flux.     

Loss cone fluxes and pitch angle diffusion at the equatorial plane during auroral radio absorption events

The flux and pitch angle distribution of energetic electrons near the loss cone have been investigated over the energy range 15–300 keV, using measurements on the geosynchronous satellite GEOS-2 at the times of auroral radio absorption events detected by riometers in Scandinavia. It is shown that conditions of strong pitch angle diffusion apply only during the most intense absorption events ( 6 dB at 30 MHz) which are relatively infrequent. During most events the loss cone is partially depleted, with the degree of depletion increasing as the absorption becomes weaker. The variation of the pitch angle diffusion coefficient with the observed radio absorption is estimated. A consequence of loss cone depletion is a tendency to overestimate the smaller events when computing the radio absorption from flux measurements in the 0°–5° range of detector pointing angles. An empirical law is derived which enables the computation of radio absorption consistent with measurements. D-region recombination laws are discussed and limits are set on the height profile of the effective recombination coefficient.     

Vertical motions in the thermosphere over Mawson,Antarctica

High time resolution measurements of Doppler shift and broadening of the (OI) >1630 nm emission in the night airglow and aurora have provided determinations of vertical velocities and temperatures in the neutral thermosphere over Mawson, Antarctica. The vertical wind exhibits a large, rapid and complex response to geomagnetic energy input. Upward winds greater than 50 m s⁻¹ are frequently associated with the expansion phase of auroral substorms. Following the disturbance, prolonged periods of downward winds produce temperature enhancements of 200K outside the source region, thus providing a mechanism for the redistribution of geomagnetic energy. Oscillatory behaviour consistent with thermospheric gravity waves is observed during both quiet and disturbed conditions.     

D-region observations of polar cap absorption events during the EISCAT operation in 1981–1989

During the years 1981–1989, 71 solar proton events altogether were observed. Dividing the events into strong, p.f.u. > 1000 (p.f.u.—proton flux measured at geosynchronous satellite orbit in units of (cm² s sr)⁻¹), medium, 100 < p.f.u. < 1000 and weak events, p.f.u. < 100, only the strong and medium events have a considerable effect on the lower ionosphere. The mean daily absorption at 30 MHz (A), measured in the auroral zone, is >2 dB during strong events, <2 dB during medium events and < l dB during weak events. The most active year during the EISCAT operation was 1989 when 23 solar proton events were observed including six strong events. Diurnal variation of the electron density in the D-region during PCA is a function of the solar zenith angle. However, south of L = 5 a minimum in absorption is observed during the noon hours. During sunrise the absorption increases simultaneously with solar elevation angle, but during sunset there is about 2 h delay between the decrease of absorption and the solar elevation angle.     

Soft X-rays from the sunlit earth's atmosphere

The HEAO-1 A-2 experiment low energy proportional counters have been used to measure the X-ray spectrum of the sunlit earth in the energy range 0.2–0.8 keV. The X-rays arise by coherent scattering of, or fluorescence of atmospheric constituents by, solar coronal X-rays incident on the atmosphere. Although the relative spectral contributions of the two processes depend upon the sun-earth-satellite geometry, fluorescent oxygen and nitrogen K X-ray emission is always important. The observed spectra were compared with calculations in order to derive the coronal temperature and emission measure, parameters that characterize the incident solar spectrum. These derived parameters agree well with the expected values for the nonflaring sun, and good agreement was obtained between measurements closely spaced in time but having a wide range of geometries and counting rates. Thus X-ray observations of the sunlit earth's atmosphere can be a useful monitor of solar activity for satellite-borne instrumentation unable to view the sun directly. The total measured fluorescent line flux agreed well with calculations, but the N : O line ratio did not. This disagreement is attributed to several causes, including the relative weakness of N emission at high altitudes where fluorescence is particularly important, the presence of line emission in the solar spectrum, and possible small calibration errors. Since present detectors cannot resolve the oxygen and nitrogen K X-rays, observation of X-ray fluorescence from the sunlit atmosphere will be useful in monitoring atmospheric constituents only to the extent that total line counting rates depend upon composition. X-rays from the sunlit earth are briefly examined as a source of background in auroral X-ray observations. During nonflare periods this background should be unimportant above about 2 keV.     

A theoretical study of the distribution of ionization in the high-latitude ionosphere and the plasmasphere: first results on the mid-latitude trough and the light-ion trough

A model of the O⁺ and H⁺ distributions in the plasmasphere and high-latitude ionosphere is described and first results are presented. The O⁺ and H⁺ continuity and momentum equations are solved from the F-region to the equatorial plane in the inner plasmasphere, and to an altitude of 1400 km in the outer plasmasphere and high-latitude ionosphere. Account is taken of high-latitude convection, departure from corotation inside the plasmasphere, and neutral air winds. The neutral air winds are consistent with the assumed convection pattern. For equinox and magnetically quiet conditions the calculations show that a mid-latitude trough in F-layer electron concentration is present from 1600 to 0600 LT and the trough may occur either inside or outside the plasmasphere. The movement of the trough in this period is from higher to lower latitudes and is in qualitative agreement with AE-C and ESRO-4 data. A light-ion trough feature is apparent in the H⁺ distribution in the topside ionosphere at all local times. During the day the upward H⁺ flow increases with latitude to produce the light-ion trough. At night the H⁺ trough may be directly produced by the occurrence of the mid-latitude O⁺ trough. The relationships between the position of the plasmapause and the trough are discussed. Also discussed are the influence of particle ionization in the auroral zone and the effect of the neutral air wind.     

Optical measurements of a nightside poleward expanding aurora

Coordinated optical observations were performed from the poleward side of the midnight auroral oval. Height measurements of the auroral emissions at 4278, 5577 and 6300 Å, as well as their intensity ratios in the poleward expanded auroral substorm, have been carried out. The findings indicate a significantly softened electron spectrum compared with similar data from the equatorward part of this substorm. Typical values for the poleward expanded aurora are 300 eV and lower, while keV particles dominate the auroras at 10° lower latitudes. Emission altitudes and spectral characteristics are comparable to the transient burst emissions frequently observed from the same site in the post-noon sector, i.e. within the cusp.The 6300 Å atomic oxygen emission is used as a tracer of F-region wind and temperature. Interferometer observations show that there exists a prevailing crosspolar antisunward wind, increasing with geomagnetic activity to several hundred m s⁻¹. The temperature shows an increase of 150 K associated with high geomagnetic activity.     

A comparative study of H+ and He+ at sunspot minimum and sunspot maximum

The time-dependent equations of continuity and momentum for O⁺, H⁺ and He⁺ are solved for the section of a mid-latitude flux tube from the lower F2-region to the equator. H⁺ and He⁺ behaviour is compared for both sunspot minimum and sunspot maximum conditions at equinox concentrating on light ion replenishment following a magnetic storm and light ion fluxes.It is found that the He⁺ fluxes vary little from day to day throughout the replenishment period following a magnetic storm. In contrast, the H⁺ fluxes can vary considerably throughout the replenishment period, particularly at sunspot minimum. In addition, our results show that temperatures above 1000 km have an important influence on the variation and magnitude of H⁺ fluxes and consequently also influence the replenishment of H⁺. At sunspot minimum the He⁺ content reaches its maximum value after eight days of replenishment and then starts to decrease whilst the H⁺ content continues increasing. This is explained in terms of night-time H⁺ return flows dragging the He⁺ to regions where it can readily recombine. Consequently lightion depletion following a magnetic storm provides a loss process for neutral helium at sunspot maximum but not necessarily at sunspot minimum.     

Observation and analysis of NI 520.0 nm auroral emissions

The results from the analysis of simultaneous auroral ground-based optical measurements of the N(²D) 520.0 nm, N₂⁺ 1NG 470.9 nm, O(³P) 844.6 nm and the O(¹D) 630.0 nm emission intensities are presented. The data were obtained during auroral observations at Gillam (56.35°N, 265.32°E) over an observation period of about 8 hours, from UT 2:33 hrs to UT 10:06 hrs, on 20 March 1985. The soft electron flux measurements on board the DMSP satellite for the time of the experiment have also been considered in the analysis. The N(²D) density and the N(²D) 520.0 nm integral emission rate I(520.0) were calculated employing a one- and two-dimensional time-dependent ion-chemistry model and the model predictions have been compared with the experimental I(520.0) nm emission rates. It was found that the model predictions of the NI I(520.0) nm intensity based on the electron energy fluxes inferred from the experimental I(844.6)/I(427.8) emission rate ratios are smaller in magnitude than the experimental values by a factor of 5–8 after allowing for horizontal transport of [N(²D)] by neutral winds. Assuming soft electron precipitation, suggested by the OI I(630.0) nm emission measurements and the DMSP satellite electron flux data, provided good agreement between the model and experimental results. Based on the results obtained it was concluded that horizontal transport played a minor role and that the observed N(²D) I(520.0) nm emissions were mostly produced by precipitating soft electron fluxes with energies below about 100 eV.     

The suppressing effect of field-aligned currents on VLF emissions in the magnetosphere

The growth rate of whistler-mode waves is calculated analytically for a bi-Maxwellian plasma in the presence of a beam of cool electrons. This beam is moving in the same direction as the gyroresonant electrons and in the opposite direction to the waves which are considered to propagate parallel (or anti-parallel) to the imposed geomagnetic field. A somewhat surprising result is found. This is that even if the anisotropy is greater than a critical value, which is strongly frequency dependent, the beam reduces the growth of the waves near half the electron gyrofrequency. For a field-aligned current density ~ 1 μA m⁻², this mechanism can explain the lack of signals near 1.4 kHz on auroral (return current) flux tubes. It can also explain the observed absorption of signals at half the electron gyrofrequency, around 7 kHz on L = 4 flux tubes, near the equatorial plane and just outside the plasmapause.     

Comparison of ionospheric electrical conductances inferred from coincident radar and spacecraft measurements and photoionization models

Height-integrated electrical conductivities (conductances) inferred from coincident Sondrestrom incoherent scatter radar and DMSP-F7 observations in the high-latitude ionosphere during solar minimum are compared with results from photoionization models. We use radar and spacecraft measurements in combination with atmospheric and ionospheric models to distinguish between the contributions of the two main sources of ionization of the thermosphere, namely, solar UV/EUV radiation and auroral electron precipitation. The model of Robinsonet al. (1987, J. geophys. Res.89, 3951) of Pedersen and Hall conductances resulting from electron precipitation appears to be in accordance with radar measurements. Published models of the conductances resulting from photoionization that use the solar zenith angle and the solar 10.7-cm radio flux as scaling parameters are, however, in discrepancy with radar observations. At solar zenith angles of less than 90°, the solar radiation components of the Pedersen and Hall conductances are systematically overestimated by most of these models. Geophysical conditions that have some bearing on the state of the high-latitude thermosphere (e.g. geomagnetic and substorm activity and a seasonal variation of the neutral gas distribution) seem to influence the conductivity distribution but are to our knowledge not yet sufficiently well modelled.     

Observation of low energy particle precipitation at low altitude in the equatorial zone

Precipitation of protons (~ 1 MeV) in the equatorial zone was investigated by the Phoenix-1 experiment on board the S81-1 mission from May–November, 1982. The protons show a precipitation peak along the line of minimum magnetic field strength with a full width at half maximum (FWHM) of 13°. The index of equatorial pitch angle distribution is q ~ 19. The peak proton flux shows a fifth-power altitude dependence, and the proton flux shows approximately a factor of 3 times increase in 1982 compared to that in 1969 due, possibly, to the stronger (~ 1.2 times) solar maximum conditions of 10.7cm radiation in 1982.