Trimpi events on low latitude paths: An investigation of gyroresonance interactions at low L-values

We report on Trimpi events observed at Durban (L = 1.69, 29°53′S, 31°00′E) and investigate the efficacy of gyroresonance scattering in precipitating electrons into the atmosphere at low L (<2). The rate of occurrence of Trimpis at Durban is less than one per day. Our observations include a number of daytime events on OMEGA signals from La Reunion. Using the full relativistic equations of motion, a test particle simulation is employed to find the region in parameter space where large pitch angle scattering occurs. We find that at low L the conditions for pitch angle scattering are less favourable than at higher L (L ∼ 4). Resonant electrons have high (relativistic) energies, interaction times are of the order of milliseconds (T_i ∼ 5 ms) and large wave amplitudes (B_w ∼ 200 pT) are required at whistler frequencies to produce pitch angle changes of greater than 1°. Large pitch angle scattering is needed near Durban since particles near the loss cone will have been lost in the South Atlantic Geomagnetic Anomaly. We note that the radio frequencies transmitted into the magnetosphere from lightning are too low to give effective electron scattering at low L. We suggest an explanation for the low rate of occurrence of Trimpis at Durban.     

Broad-band auroral VLF hiss and inverted-V electron precipitation in the polar magnetosphere

A polar map of the occurrence rate of broad-band auroral VLF hiss in the topside ionosphere was made by a criterion of simultaneous intensity increases more than 5 dB above the quiet level at 5, 8, 16 and 20 kHz bands, using narrow-band intensity data processed from VLF electric field (50 Hz–30 kHz) tapes of 347 ISIS passes received at Syowa Station, Antarctica, between June 1976 and January 1983.The low-latitude contour of occurrence rate of 0.3 is approximately symmetric with respect to the 10–22 MLT (geomagnetic local time) meridian. It lies at 74° around 10 MLT, and extends down to 67° around 22 MLT. The high-latitude contour of 0.3 lies at invariant latitude of about 82° for all geomagnetic local times. The polar occurrence map of broad-band auroral VLF hiss is qualitatively similar to that of inverted-V electron precipitation observed by Atmospheric Explorer.(AE-D) (Huffman and Lin, 1981, American Geophys. Union, Geophysics Monograph, No. 25, p. 80), especially concerning the low-latitude boundary and axial symmetry of the 10–22 h MLT meridian.The frequency range of the broad-band auroral VLF hiss is discussed in terms of whistler Aode Cerenkov radiation by inverted-V electrons (1–30 keV) precipitated from the boundary plasma sheet. High-frequency components, above 12 kHz of whistler mode Cerenkov radiation from inverted-V electrons with energy below 40 keV, may be generated at altitudes below 3200 km along geomagnetic field lines at invariant latitudes between 70 and 77°. Low-frequency components below 2 kHz may be generated over a wide region at altitudes below 6400 km along the same field lines. Thus, the frequency range of the downgoing broad-band auroral hiss seems to be explained by the whistler mode Cerenkov radiation generated from inverted-V electrons at geocentric distances below about 2 R_E (Earth's radius) along polar geomagnetic field lines of invariant latitude from 70 to 77°, since the whistler mode condition for all frequencies above 1 kHz of the downgoing hiss is not satisfied at geocentric distance of 3 r_e on the same field lines.     

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.     

Double structure of ionospheric conductivity in the midnight auroral oval during a substorm

Measurements of precipitating particles on board DMSP F7 spacecraft are used to analyze the distribution of ionospheric conductance in the midnight auroral zone during substorms. The distribution is compared with the meridional profile of ionospheric currents calculated from magnetic data from the Kara meridional chain. Two regions of high Hall conductance are found; one of them is the traditional auroral zone, at latitudes 64–68°, and the other is a narrow band at latitudes 70–73°. The position of high conductance zones is in agreement with the location of the intense westward currents. The accelerated particle population is typical of electrons E_e > 5 keV in the high conductance region.     

Solar-induced variation of proton precipitation near the equator

This paper reports the solar condition dependences of the quasi-trapped component (low energy) of the proton population of energy 0.65–35 MeV which peaks in the equatorial zone centered on the minimum magnetic field equator in the altitude range 170–850 km. The proton populations compared pertain to AZUR observation in 1969–1970, S81-1 mission observation in 1982 and EXOS-C observation during 1984–1986. In the equatorial zone, the dependence of the flux normalization constant, which represents the absolute proton population, upon factors like L (1.1–1.3), B (0.29–0.32 gauss), latitude ( ± 20°), longitude (0–360°) and anisotropy index q (~6–12) of the pitch angle distribution function is not so significant in the given range of these factors as it is upon the solar epoch. It is found that the absolute proton flux in 1982 was, at least, forty times that in 1984–1986 and, almost, three times that in 1969–1970, possibly, due to, varying solar conditions in those epochs.     

A study of equatorial waves in the Indian zone

An equatorial wave campaign was conducted at Trivandrum (8.5°N, 77°E), Minicoy (8.3°N, 73°E) and Port Blair (11.7°N, 92.7°E) during June-July 1988. The campaign provided balloon-measured daily wind profiles at all the three stations for 48 days in the 0–30 km altitude range and rocket-measured daily wind profiles at Trivandrum for 42 days in the 31–60 km altitude range. Using these daily wind data a study was made on different equatorial wave modes present in this region. The study revealed evidence of Kelvin waves with period 12–16 days and vertical wavelength ∼ 10 km in the lower stratosphere, with period 6–9.6 days and vertical wavelength of ∼ 10–15 km in the stratospheric-lower mesospheric region and MRG waves with periods 4–4.4 days and vertical wavelength of 10 km in the upper troposphere and lower stratosphere.     

Quasi-periodic particle precipitation and Trimpi activity at Halley,Antarctica

The relationship between quasi-periodic VLF emissions and micropulsations is briefly reviewed, and then discussed with reference to an event recorded at Halley, Antarctica, on day 257 in 1986. VLF emissions at 2 kHz with a quasi-period of 9 s were observed simultaneously with Pi1 and Pe1 micropulsations. Also observed was a quasi-periodic Trimpi event on the amplitudes and phases of the VLF transmitters NAA and NSS. It is deduced that the VLF emissions are modulated in the generation region by a hydromagnetic wave, giving rise to particle precipitation. The emissions are also modulated by the bounce period of the generating particles. The Trimpi effect is due to 120 keV electrons being precipitated into the lower ionosphere by the interaction with the VLF emissions. This event shows that the Trimpi effect can be used to detect particle precipitation occurring during the ULF/VLF interaction, and can give information which helps to define the mechanism reponsible for the interaction.     

Experimental daytime VLF ionospheric parameters

Measured field strengths from VLF transmitters are used to determine improved daytime values of ionospheric parameters to enable improved VLF propagation predictions. These parameters are the traditional H′ (height in km) and β (sharpness in km⁻¹) as used by Wait and by NOSC in their Earthionosphere waveguide computer program. They are found by comparing the predictions of the NOSC program with the observed VLF field strengths over both long and short paths.Experimental observations from two nearly north-south paths are used to determine the solar zenith angle dependence of both H′ and β for low latitude (or summer mid-latitude) conditions. These results are then used to predict the daytime variations in VLF field strengths with solar zenith angle (and hence time) on other suitable paths and good agreement is found with measurements made on these paths.The absolute value of β for overhead Sun is found to be 0.45 km⁻¹ and is principally determined by the attenuation on the very long, west to east, fully sunlit, 14.4 Mm path from NWC (Australia, 22°S) to San Francisco (37°N), after applying small corrections for the solar zenith angle variations along the path at midday. Further support is obtained from results from the 8.6 Mm path NDT (Japan) to San Francisco, an 8.2 Mm path NPM (Hawaii) to New Zealand, and an east to west 7.5 Mm path from NPM to Townsville, Australia. The conditions studied are solar maximum. The frequencies studied are 15–30 kHz.     

Tropical behavior of mesospheric ozone as observed by SMM

The seasonal behavior of low latitude mesospheric ozone, as observed by the SMM satellite solar occultation experiment, is detailed for the 1985–1989 period. Annual as well as semi-annual waves are observed in the 50–70 km altitude region. In the latitude range of ±30 the ozone phase and amplitude are functions of temperature and seasonal changes in solar flux. Temperature is the controlling factor for the equatorial region and seasonal changes in solar flux become more dominant at latitudes outside the equatorial zone (greater than ±15). There is a hemispheric asymmetry in the ozone annual wave in the 20 30 region, with northern hemispheric ozone having a larger amplitude than southern hemispheric ozone. In this region temperature is nearly in phase with ozone in both hemispheres and is reduced in amplitude in the northern hemisphere. The equatorial region is characterized by a strong semi-annual wave in addition to the annual variation, and temperature is nearly out of phase with ozone. At all latitudes there is a larger ozone concentration at sunrise than at sunset. The sunrise sunset difference increases with increasing altitude     

Study of equatorial plasma bubble dynamics using GHz scintillation observations in the Indian sector

To study equatorial plasma bubble dynamics, telemetry signals (4 GHz) were recorded simultaneously from two geostationary satellites. INSAT-1B (74°E) and INSAT-1C (94°E) at Sikandarabad satellite Earth station (dip 42.0°) from January to December 1989 and at the Chenglepet satellite Earth station (dip 10.5°) during September–October 1989 along the same geomagnetic meridian. The characteristics and occurrence pattern of the scintillations suggest that these are equatorial plasma bubble induced events. Observations from the two satellites recorded simultaneously at each of these locations were utilized to estimate the east-west plasma bubble irregularity motion. Plasma bubble rise velocities over the magnetic equator were calculated from the systematic onset time differences observed between an equatorial and a low latitude station. The east-west plasma bubble velocity estimated at Sikandarabad, corresponding to 1200 km altitude in the equatorial plane, shows a night time variation pattern with a peak at around 2100 LT. The mean values over Chenglepet, which correspond to 400 km altitude, start decreasing right from 1900 LT and seem to be influenced by the plasma bubble rise velocities. The differences in magnitude observed between the present results and those reported elsewhere by other techniques are interpreted in terms of vertical shears in the plasma zonal flow over the equator. The near alignment of the two observing stations along a common geomagnetic meridian and the simultaneous use of two satellites located twenty degrees apart in longitude provided an excellent data base to study plasma bubble dynamics.     

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.     

Parachute measurements of positive ion density of the middle atmosphere over the dip equator by spherical probe

A parachute-borne gridded spherical probe has been used to measure the total positive ion density. Two launches were made, using Soviet M-100 rockets, on 22 and 29 April 1987, at 1200 UT, from an equatorial station, Thumba (8°N, 76°E) India. Data were obtained for the altitude region 10 to 80 km. A broad maximum around 15 km and a broad minimum around 60 km have been noticed in the ion current profiles obtained in both flights. The theory of the operation of the probe has been given. A detailed discussion of the results obtained has also been included.     

A model for the electron temperature variation along geomagnetic field lines and its effect on electron density profiles and VLF paths

A realistic model for the temperature variation along geomagnetic field lines is described. For high altitudes (>1500 km) the temperature is taken to increase as the nth power of radial distance (n−2), giving temperatures consistent with those measured in situ by high altitude satellites. For realistic temperatures at low altitude an extra term is included. The temperature gradient along the field line is then 0.9–1.6° km⁻¹ during the day and 0.5–0.7° km⁻¹ during the night at 1000 km, reducing to about half these values at 2000 km, for the latitude range 35–50°. This is consistent with calculations made from nearly simultaneous satellite measurements at 1000 and 2500 km. It is shown that assuming diffusive equilibrium, including the new temperature model, more realistic equatorial electron density profiles result than for isothermal field lines.The temperature gradient model is also purposely formulated to be of a form that enables the temperature modified geopotential height to be obtained without numerical integration. This renders the model particularly suitable for ray-tracing calculations. A ray-tracing model is developed and it is shown that unducted ray paths are significantly altered from the corresponding paths in an equivalent isothermal model; there is greater refraction and magnetospheric reflection takes place at lower altitudes. For summer day conditions, an inter-hemispheric unducted ray path becomes possible from 26° latitude that can reach the ground at the conjugate.     

Numerical simulation of the penetration and reflection of a whistler beam incident on the lower ionosphere at very low latitude

By the full-wave algorithm with Fourier synthesis, 3-D propagation of a whistler beam incident on the pre-dawn lower ionosphere at very low latitude is numerically investigated. Processes of transmission, reflection, and coupling with the Earth-ionosphere waveguide are discussed via the wave energy and polarisation distributions and their dependence on the wave parameters and the ionospheric profile (such as the E_s-layer). It is shown that the dominant wave above 90 km altitude has the propagation characteristics of the magneto-ionic whistler mode, and absorption, spreading, reflection and mode conversion mainly occur at, and are greatly affected by, the bottom of the ionosphere. It is found that the transmitted energy density along the Earth's surface is reduced by 20 dB or more. Beam transmission loss varies asymmetrically with the incident angle, but changes little with the frequency. In the region 150 km (for 5 kHz) away from the ‘exit area’ where whistlers emerge, the bearing measurements using ground-based VLF direction-finders may be in error because direction-finding algorithms assume plane wave propagation. Only a small portion (about −25 dB at 5 kHz) of the incident energy is reflected up to an altitude of 150 km, and major reflection takes place in a small range of altitude at the bottom of the ionosphere with little spreading and lateral shift with respect to the incident beam. Reflection is enhanced considerably at lower frequency. Our results also suggest that an E_s-layer or an ionospheric gradient refracting waves to higher latitudes would be favorable factors for multi-hop echoes to be received on the ground.     

Simultaneous observations of travelling ionospheric disturbances by two-dimensional radio interferometry and the differential Doppler technique applied to satellite signals

The main object of the campaign reported here was to compare TID characteristics obtained from two essentially different observation techniques: (1) observation of the apparent angular position shifts of Virgo A by the Nançay radioheliograph (47.33°N, 2.15°E) gave azimuths and periods of travelling ionospheric disturbances (TIDs); (2) differential Doppler shifts of signals from NNSS-satellites recorded simultaneously at Tours (47.35°N, 0.70°E), Nançay and Besançon (47.32°N, 5.99°E) provided azimuths and latitudinal wavelengths. Observations were made during the period 10–30 November 1987, between 6 and 12 h UT. It is found that azimuths obtained from the two techniques are consistent if sufficient averaging over wave trains is performed: averaging over several hours for radio interferometry and averaging over the whole satellite trace for the differential Doppler technique. Averaging is necessary because of (1) the intrinsic dispersion in wave azimuth, (2) the broadness of observed wave spectra and the dispersive properties of gravity waves, and (3) the spatial separation of ionospheric points for the two techniques. Good agreement between the azimuths was achieved by setting the altitude of the TIDs, which is used in the differential Doppler analysis, to about 250 km, appreciably lower than the maximum in electron density (about 350 km). The mean azimuth of observed TIDs was 12° East from South with a standard deviation of about 30°. The dominant period and horizontal wavelength of the observed TIDs were 40 min and 450 km. The East-West coherence length of the TIDs was found to be only of the order of 200 km.     

Upper stratospheric and mesospheric temperatures derived from lidar observations at Aberystwyth

From lidar observations of relative atmospheric density above Aberystwyth (52.4°N, 4.1°W) upper stratospheric and mesospheric temperatures have been derived for a total of 93 nights between December 1982 and February 1985. Excellent agreement was found between radiances synthesised from these temperatures and those measured by satellite-borne instruments. Summer temperatures showed a smooth and regular variation with altitude and reasonably good agreement with the CIRA (1972) model atmosphere. By contrast, winter temperatures showed a much greater variability with altitude and greater changes from night to night, with the frequent occurrence of a large amplitude wave-like perturbation in the mesosphere with about 15 km vertical wavelength and amplitude about 20K between 60 and 80 km.Pronounced warmings of the stratosphere were observed during the three winters of observation. During the warming event occurring in early February 1983 the stratopause temperature increased to 303K at 43 km, while the major warming event of late December 1984/early January 1985 produced a stratospheric temperature gradient of 16K km⁻¹ between 34 and 36 km. During the latter event a distinct local temperature minimum at 32.6 km was observed on New Year's Eve, this descending to 29 km by the following night and being accompanied by a lowering of the stratopause from 43 to 38.5 km in the same period. These results demonstrate the ability of the present technique to resolve the high stratopause temperatures and steep stratospheric temperature gradients which occur during stratospheric warmings, in marked contrast to the limited resolution achieved by satellite experiments.     

Balloon measurements of stratospheric ion conductivities over the tropics

Conductivity measurements of negative and positive ions were made from about 20 to 35 km by two identical balloon-borne spherical probes at Hyderabad (17.5°N, 78.6°E), India on 22 April 1989 and 22 December 1990. One balloon was launched at 0158 h IST (Indian Standard Time) which reached its ceiling around 0330 h IST. After that time, it floated for about 3 h, 1.5 h before sunrise and 1.5 h after sunrise. Thus it gave data for both day- and night-time conditions at float altitude. The other balloon was launched at 0535 h IST. It gave data for daytime only. Several interesting results have been obtained at the float altitudes. During the night, in the flight of 22 April 1989 the conductivity values of positive ions were found to be about 1.5 times those of negative ions at the float altitude. During the day, in the flight of 22 April 1989, the positive ion conductivity values were found to increase with the increase of solar elevation angle at around 37.5 km altitude. The negative ion conductivity values, however, did not show any day-night variation. In the flight of 22 December 1990, these features were not seen. Instead, a pocket was found where conductivity values were very high (of the order of 10⁻¹¹ mho m⁻¹) at an altitude of about 32.5 km. Also in this flight, the positive ion conductivity was always found to be approximately equal to that of the negative ion conductivity.     

Neutral density cells in the high latitude thermosphere—1. Solar maximum cell morphology and data analysis

A new high latitude thermospheric neutral density structure has been revealed in NCAR-TIGCM simulations at 120–350 km altitude. The structure consists of density cells above 50° latitude with radii of approximately 1000 km. There are between two to four cells present depending on the altitude and magnetic activity. For example, at 200 km under magnetically active conditions, the density structure consists of four cells: low density cells are located near dawn and dusk and high density cells are located near noon and midnight. Density variations among the cells range from 5 to 50% for magnetically quiet and active conditions respectively. The cells are present at all seasons, for a wide range of magnetic activity levels, and at solar minimum and solar maximum. The density cell morphology is established for equinox solar maximum as a function of altitude and magnetic activity. Departures of the cell structure from this morphology due to seasonal and solar cycles are discussed. The cell morphology provides a new framework in which to interpret lower thermospheric density data. Data to test and confirm the model predictions were provided by the SETA-1 satellite.     

Sudden neutral sodium layers: a strong link to sporadic E layers

A sodium LIDAR instrument located at Andenes, Norway (69°N; 16°E) observed several sudden developments of narrow sodium layers in the 90–100 km altitude region. These layers grow with typical time constants of 5 min and have a width of 1 km in altitude. We present the temporal and spatial properties for a number of these events. In a first step towards identifying the processes which create these layers we study the correlation of the growth phase of sudden sodium layers and of sporadic E layers. The latter have been recorded by an ionosonde located 129 km east of the LIDAR site. Within the mutual altitude and time resolution available in our common records a strong correlation of simultaneous occurrence of sudden sodium layers and E_s−l layers is observed, which establishes a strong link between the formation of the two types of layers. We further discuss processes which potentially could give rise to the formation of sudden sodium layers.     

A simple theoretical model of the diffraction of VLF electromagnetic waves by the Antarctic icecap

The Antarctic continent has been modelled as a spherical cap whose pole is coincident with that of the South Pole, which totally absorbs VLF radio waves attempting to propagate over it. The propagation of Omega navigation signals around this model icecap has then been computed using Kirchhoff diffraction theory. Spherical caps extending to 66.5 and 75.5°S have been found to accurately model the signals from Omega La Reunion and Argentina, respectively, received on flights between Christchurch, New Zealand and Scott Base in Antarctica, up to the boundary of the theoretical icecap. These model icecaps were found to be good fits to the boundary of the Antarctic continent, when measured at the 1–1.5 km contour of ice thickness, in the region where the VLF waves diffracted around the icecap. The good agreement obtained between the experimental field strength data and those computed theoretically, using only simple diffraction theory, suggests that coastal refraction plays at most only a secondary role in circumpolar propagation.