A case study of a high latitude ionospheric electron density depletion

Incoherent scatter measurements inside and outside an ionospheric electron density depletion are described. The density depletion is probably caused by an enhancement of NO⁺-ions and subsequent dissociative recombination. The NO⁺-ions are increased because high electric fields present at the geographical location of the density depletion speed up the reaction O⁺ + N₂ → NO⁺ + N. The electron as well as the ion temperature within the density depletion are strongly enhanced, the latter due to Joule heating, also caused by the electric field.     

A study of ionospheric electron density deviations during two great storms

A large foF2 data base has been collected for the two great storms occurring in March and October 1989. Plots of foF2 deviations versus time and latitude show the large depressions cannot be caused by a composition change alone; the storm induced dynamo electric field must also play an important role, especially for low latitudes. The observed ionospheric response is hemispherically asymmetric, with the autumnal hemisphere suffering a longer lasting, more latitudinally extended, and deeper depression.     

HF produced ionospheric electron density irregularities diagnosed by UHF radio star scintillations

High frequency waves incident on an overdense ionosphere (i.e. HF < penetration frequency) are known to produce large-scale irregularities with scale sizes of several hundred meters in the F-region of the ionosphere. Three observations of radio star intensity fluctuations at UHF are reported for HF ionospheric modification experiments performed at the Arecibo Observatory. Two observations at 430 MHz and one observation at 1400 MHz indicate that the thin phase screen theory is a good approximation to the observed power spectra. However, the theory has to be extended to include antenna filtering. Such filtering is important for UHF radio star scintillations since the antenna usually has a narrow beam width. HF power densities of less than 37 μW m⁻² incident on the ionosphere produce electron density irregularities larger than 13% of the ambient density (at 260 km) having scale sizes of ~510 m perpendicular to the geomagnetic field. The irregularities form within 20–25 s after the HF power is turned on. From the observed power spectra driftvelocities of the irregularities can be estimated.     

The effects of ionospheric horizontal electron density gradients on whistler mode signals

Whistler mode group delays observed at Faraday, Antarctica (65° S, 64° W) and Dunedin, New Zealand (46° S, 171° E) show sudden increases of the order of hundreds of milliseconds within 15 minutes. These events (‘discontinuities’) are observed during sunrise or sunset at the duct entry regions, close to the receiver's conjugate point. The sudden increase in group delay can be explained as a tilting of the up-going wave towards the sun by horizontal electron density gradients associated with the passage of the dawn/dusk terminator. The waves become trapped into higher L-shell ducts. The majority of the events are seen during June-August and can be understood in terms of the orientation of the terminator with respect to the field aligned ducts. The position of the source VLF transmitter relative to the duct entry region is found to be important in determining the contribution of ionospheric electron density gradients to the L-shell distribution of the whistler mode signals.     

The global distribution of ionospheric small-scale irregularities from topside sounding data

This paper examines the global distribution of electron density irregularities with scales of the order of several tens to hundreds of meters in the ionosphere by using topside sounder data from the COSMOS-1809 satellite obtained in May–June and December 1987. The diffuse traces of Z-waves on topside ionograms in a frequency band just below the upper hybrid resonance are used for diagnostics. These traces are attributed to the scattering of sounder-generated ordinary and slow extraordinary mode waves.     

Estimation of the true ionospheric height profile,with a continuous gradient,from oblique sounding data

A technique to calculate the true height ionospheric profile from multifrequency oblique soundings is proposed. In the classic isotropic statement of the problem, it avoids a discontinuity in the gradient of electron density with height at nodal points and. as shown in test examples, it considerably increases the accuracy of the reconstructed profile as compared with well known, currently employed methods. The influence of the geomagnetic field is taken into account in the technique by constructing the equivalent isotropic ionogram. Both numerical and experimental results are presented and discussed.     

Coherent-array HF Doppler sounding of traveling ionospheric disturbances: I. Basic technique

We present an introduction to the use of phase-coherent, multi-receiver HF Doppier sounding arrays for measuring the horizontal velocity of traveling ionospheric disturbances (TID's). The point of departure is the theorem of Pfister (1971, J. atmos. terr. Phys. 33, 999) relating ray Doppler to ray zenith angle for a monostatic full reflection sounder. Retaining the simple model of a specular, smooth ionospheric reflector which is deformed by a propagating undulation, we first generalize the theorem to bistatic sounding geometry and then include the effects of amplitude in addition to phase. Next, these results are cast into an algorithm for treating multi-receiver phase sounders containing many diverse baselines, in order to obtain an accurate and unambiguous solution in the plane of wave slowness (inverse of velocity). The point spread function of this solution is controlled by process bandwidth and by array geometry. We illustrate the coherent-array approach using data from an eight receiver array during passage of a TID.     

Coherent-array HF Doppier sounding of traveling ionospheric disturbances: II. Assessment of errors

We examine errors which are likely to occur in the compact-array HF sounding technique for determining the velocity of traveling ionospheric disturbances. The errors arise from ionospheric irregularities (either large-scale specular facets, or fine-scale scatterers) and from the presence of multiple disturbance undulations propagating simultaneously with different velocities. The errors are either systematic (that is, not causing internal inconsistencies within the array velocimetry data) or random (that is, causing quantifiable inconsistencies between alternative velocity estimates furnished by the array). The random errors can be directly assessed using array data, as is illustrated for actual data from a 16-channel array.     

Modeling the low latitude ionospheric total electron content

The undisturbed ambient total electron content of the ionosphere in the equatorial region exhibits two characteristic features:

• 1.(i) a longitudinal behavior of the post-sunset variation of the ionization near the crests of the equatorial anomaly
• 2.(ii) an enhancement at lower latitudes following the post-sunset decay. During high solar activity periods the southern crest of the equatorial anomaly in the African longitude sector is characterized by a post-sunset maximum often exceeding the afternoon maximum. In the Indian and other longitude zones, the post-sunset peak is not so prominent. Instead, a ledge is obtained in the corresponding local time period. At lower magnetic latitudes, the ionization decays very rapidly around sunset, but an enhancement lasting 2–4 h is observed afterwards.

Numerical solution of the plasma continuity equation, including the effects of ionization production by solar ultraviolet radiation, loss through charge exchange and transport by diffusion, electrodynamic drift and neutral wind, has been used to investigate the above two features. It is found that the pre-reversal peak of the E × B drift at the magnetic equator around sunset is the dominant mechanism responsible for the post-sunset behavior near the crests of the equatorial anomaly. The zonal wind causes an asymmetry of the total content in the northern and southern hemispheres. In African longitudes, where the magnetic declination is about 20°W, the southern crest is more developed at the expense of the northern counterpart. The north-south asymmetry is practically absent in the Asian sector, with its low (< 5°) declination angle. In the Pacific area, an easterly declination (about 9°E) results in a higher post-sunset ionization at the northern crest, although the asymmetry is less pronounced than that in the African zone. The night-time enhancement at lower latitudes has been found to be controlled by the post-sunset increase in the vertical drift, possibly also modulated by the neutral wind.     

Short term variabilities of ionospheric electron content (IEC) and peak electron density (NP) during solar cycles 20 and 21 for a low latitude station

Hourly values of IEC and of f₀F₂ (critical frequency) for a low latitude station, Hawaii (21.2°N, 157.7°W), during the solar maxima (1969 and 1981) and minima (1965 and 1985) years of two consecutive solar cycles, 20 and 21, are used to study the day to day variabilities of the ionospheric parameters IEC and NP. It is found that there is good correspondence in the day to day variations of IEC and NP from one solar cycle to the other for both solar maximum and minimum years in the two solar cycles. Depending on solar phase and season, while the mean daytime IEC and NP variations range from about 20% to 35%, the mean night time values vary from about 25% to 60%. The mean daytime variations in NP for the solar minimum phase are remarkably higher in all the three seasons compared to the solar maximum phase. However, no such increase is observed in the mean daytime IEC variations, indicating the highly variable nature of the daytime ionospheric F region compared to the topside during solar minimum for this low latitude station. The winter night time IEC also seems to be a relatively stable parameter during the solar minimum. The short term day to day variabilities of the day time peak values of IEC and NP (ie IEC_max and NP_max) are not closely associated with the variations in F_10.7 solar flux. Contrary to the common expectation, the variabilities in both the parameters, particularly in NP_max, are somewhat reduced during the solar maximum (when the variability in F_10.7 solar flux is much higher compared to the solar minimum) which is more evident in the stronger 21 solar cycle. A larger number of significant components are seen in the spectra of the percentage variation of both IEC_max and NP_max during both solar phases of the two solar cycles compared to the corresponding F_10.7 solar flux spectra. The number of additional components for both the parameters with periods less than 15 days are more for the low solar activity years than for the solar maximum years.     

GPS satellite measurements: ionospheric slab thickness and total electron content

A preliminary analysis was made of ionospheric slab thickness, τ, and total electron content, TEC, for southern Australia using GPS satellite measurements. It was found that at mid-latitudes τ has similar overall diurnal, seasonal and latitudinal variations in the southern hemisphere as in the northern hemisphere. However, there are appreciable differences between τ in the two hemispheres which would justify appropriate modifications to ionospheric models based on northern hemisphere data before being applied confidently to the southern hemisphere. The usefulness of GPS satellites together with ionosondes over a spread of latitudes was demonstrated in determining long-term variations of TEC and τ over a large area. It was concluded that as few as four GPS receivers could provide TEC for the whole of Australia in real-time, though approximately six receivers in convenient locations would be required in practice.     

Combining TID observations: NNSS and radio interferometry data

The determination of the kinematic parameters of travelling ionospheric disturbances (TIDs) is an essential basis for their physical understanding and investigation of their possible excitation mechanisms. This paper describes a method to determine these parameters on the basis of combined radio interferometric observations and measurements of differential Doppler shift observed for NNSS satellites.The analysis of Dutch interferometric and differential Doppler observations, based on this method, shows that the average direction of propagation of a sample of medium scale TIDs observed during the first three months of 1982 and 1983 is towards about 10° west of south; the speeds and horizontal ‘wavelengths’ are typical for medium scale TIDs. The results are in satisfactory agreement with those of similar research performed at Nançay (France) more than 4 yr later and using a different methodology, but not with results obtained at Nançay during the same winter months as used in the present study. Differences may indicate that TID parameters vary with geographic location and time.     

A study of gravity waves in ionospheric electron content at L = 4

Using satellite radio beacon transmissions, travelling ionospheric disturbances have been observed in the electron content at L = 4. Waves are a common feature at this latitude, present for at least 98% of all daylight hours. The amplitude is usually 1–4% of the mean electron content and periods range between 15 and 90 minutes. Simultaneous observation of two satellite beacons, giving an effective east-west separation of 350 km, indicated apparent east-to-west velocities of 200–700 m/s.A search was made for a likely source of the waves, using data from magnetometers and riometers, from incoherent scatter radar measurements of Joule heating, and from orbiting satellite measurements of electron influx, but no definite source could be established.It is also shown that travelling disturbances are closely related to occurrences of spread-F on ionograms at high latitudes.     

Changes of the electron concentration profile during local heating of the ionospheric plasma

Using numerical simulation of a non-stationary problem of thermodiffusion and diffusive spreading of the electron component of the dense cold ionospheric plasma, the processes of formation and relaxation of strong disturbances of the electron temperature and concentration in the E- and F-regions of the middle-latitude ionosphere are examined, taking into account the altitudinal distribution of the electron transport coefficients. The cases of local heating and heating at separated altitudes of the ionospheric plasma by powerful radio waves generated from ground-based HF-facilities are numerically investigated. The numerical simulations of the non-stationary problem are compared with the analytical evaluations carried out for the stationary and quasi-stationary heating models. Results obtained from numerical experiments give good explanations of the experimentally observed deformation of the altitudinal ionospheric plasma density profile and the creation of negative cavities in the upper ionosphere and positive cavities in the lower ionosphere during the process of plasma heating.     

On medium scale travelling ionospheric disturbances as seen through total electron content power spectra

Modulation phase of 140 MHz with respect to 360 MHz signals from ATS 6 satellite recorded at Slough (51.5°N, 0.6°W geographic latitude and longitude; 54.3° geomagnetic latitude) at one-minute interval are power spectral analysed to derive dominant periodicities corresponding to Travelling Ionospheric Disturbances. From the significant peaks in the spectra, an occurrence peak at periods between 10 and 15 min and a secondary peak at 60–65 min are seen. From cross-spectral analysis of the same records from three stations separated by a few hundred kilometers for a short period, the speed and azimuth of the propagating disturbances are determined. During the day-time, most of the waves in the period range 30–100 min are seen to propagate at azimuths of 90–160°. At night-time they propagate poleward. Theoretical computations of the azimuth response of TEC to typical gravity waves, including the effect of neutral winds, show that the observed azimuths of propagation are in reasonably good agreement with theoretical predictions.     

Comparison of incoherent scatter observations of electron density,and electron and ion temperature at Millstone Hill with the International Reference Ionosphere

Millstone Hill incoherent scatter (IS) observations of electron density (N_e, electron temperature (T_e) and ion temperature (T_i) are compared with the International Reference Ionosphere (IRI-86) for both noon and midnight, for summer, equinox and winter, at both solar maximum (1979–1980) and solar minimum (1985–1986). The largest difference inN_e is found in the topside, where values of N_e given by IRI-86 are generally larger than those obtained from IS measurements, by a factor which increases with increasing height, and which has a mean value near two at 600 km. Apart from the bottom of the profile, which is tied to the CIRA neutral temperature, the IRI-86 T_e model has no solar cycle variation. However, the IS measurements during the summer reveal larger T_e at solar maximum than at solar minimum. At other seasons higher T_e at solar maximum occurs only during the daytime at the greater heights. Nighttime T_e is shown by the IS radar to be generally larger in winter than in summer, an effect not included in the IRI. This is apparently due to photoelectron heating during winter from the sunlit ionosphere conjugate to Millstone Hill. The day-night difference in T_i given by IRI-86 above 600km is not as large in the IS measurements.     

The relationship between electron temperature and electron density under quiet conditions in the auroral zone

A theoretical model is described which predicts electron temperature in the day-time F-region above EISCAT on geomagnetically quiet days, given the observed values of electron concentration, ion temperature and heat conduction, the daily average of solar radiation at 10.7cm and the MSIS-86 model of the neutral atmosphere. The values predicted by the model agree very closely with the observed temperatures, both for average conditions and for individual days. On two occasions the onset of a geomagnetic disturbance after a period of quiet conditions was accompanied by a growing divergence between the predicted and observed values, which corresponds to an additional source of electron heating such as would be provided by low energy particle precipitation.     

Electric fields and electron density irregularities in the equatorial electrojet

Two Centaure rockets were launched from Thumba (0 47′S dip). India, with a new arrangement of double probe sensors for the simultaneous measurements of the irregularities in the electron density and the electric field along and perpendicular to the spin axis of the rocket. These experiments were carried out during the period when type I irregularities were observed with the VHF backscatter radar at Thumba. Irregularities with scale sizes ranging from a few meters to a few kilometers in the electron density and in the electric field components both in the east-west and the vertical direction could be studied with these experiments. Irregularities in the electric field in the medium scale size range (30–300 m) were observed with peak to peak amplitudes up to 20 mV m⁻¹ and in the small scale (⩽ 15 m) with peak to peak amplitudes up to 5 mV m⁻¹. Horizontally propagating waves with horizontal scale sizes up to 2.5 km were observed in the region below 105.5 km. Using linear theory for the electrojet irregularities, it was found that for 5 % perturbations in the electron density, the amplitude of the electric field can be as large as 20–30 mV m⁻¹. The spectrum of the irregularities in the vertical electric field in the rocket frame of reference was calculated and it was found that for the range of scale sizes between 10 and 70 m, the mean spectral index was −2.7 and −2.6. while in the scale size range 2–10 m it was −4.0 and −5.1 for the flights C-77 and C-73, respectively.