The effect of rapid changes in ionospheric flow on velocity vectors deduced from radar beam-swinging experiments

The effects on the horizontal ionospheric velocity vectors deduced from radar beam-swinging experiments, which occur when changes in the flow take place on short time scales compared with the experiment cycle time, are analysed in detail. The further complications which arise in the interpretation of beam-swinging data, due to longitudinal gradients in the flow and to field-aligned flows, are also considered. It is concluded that these effects are unlikely to seriously compromise statistical determinations of the response time of the flow, e.g. to changes in the north-south component of the IMF, such as have been recently reported by Etemadiet al. (1988, Planet. Space Sci.36, 471), using EISCAT ‘Polar’ data.

The influence of anisotropic F region ion velocity distributions on ionospheric ion outflows into the magnetosphere

The contribution to the field-aligned ionospheric ion momentum equation, due to coupling between pressure anisotropy and the inhomogeneous geomagnetic field, is investigated. We term this contribution the “hydrodynamic mirror force” and investigate its dependence on the ion drift and the resulting deformations of the ion velocity distribution function from an isotropic form. It is shown that this extra upforce increases rapidly with ion drift relative to the neutral gas but is not highly dependent on the ion-neutral collision model employed. An example of a burst of flow observed by EISCAT, thought to be the ionospheric signature of a flux transfer event at the magnetopause, is studied in detail and it is shown that the nonthermal plasma which results is subject to a hydrodynamic mirror force which is roughly 10% of the gravitational downforce. In addition, predictions by the coupled University College London-Sheffield University model of the ionosphere and thermosphere show that the hydrodynamic mirror force in the auroral oval is up to 3% of the gravitational force for Kp of about 3, rising to 10% following a sudden increase in cross-cap potential. The spatial distribution of the upforce shows peaks in the cusp region and in the post-midnight auroral oval, similar to that of observed low-energy heavy ion flows from the ionosphere into the magnetosphere. We suggest the hydrodynamic mirror force may modulate these outflows by controlling the supply of heavy ions to regions of ion acceleration and that future simulations of the effects of Joule heating on ion outflows should make allowance for it.

The dependence of high-latitude dayside ionospheric flows on the North-South component of the IMF: A high time resolution correlation analysis using EISCAT “Polar” and AMPTE UKS and IRM data

In 1984 and 1985 a series of experiments was undertaken in which dayside ionospheric flows were measured by the EISCAT “Polar” experiment, while observations of the solar wind and interplanetary magnetic field (IMF) were made by the AMPTE UKS and IRM spacecraft upstream from the Earth's bow shock. As a result, 40 h of simultaneous data were acquired, which are analysed in this paper to investigate the relationship between the ionospheric flow and the North-South (Bz) component of the IMF. The ionospheric flow data have 2.5 min resolution, and cover the dayside local time sector from ∼ 09:30 to ∼ 18:30 M.L.T. and the latitude range from 70.8° to 74.3°. Using cross-correlation analysis it is shown that clear relationships do exist between the ionospheric flow and IMF Bz, but that the form of the relations depends strongly on latitude and local time. These dependencies are readily interpreted in terms of a twinvortex flow pattern in which the magnitude and latitudinal extent of the flows become successively larger as Bz becomes successively more negative. Detailed maps of the flow are derived for a range of Bz values (between ± 4 nT) which clearly demonstrate the presence of these effects in the data. The data also suggest that the morning reversal in the East-West component of flow moves to earlier local times as Bz, declines in value and becomes negative. The correlation analysis also provides information on the ionospheric response time to changes in IMF Bz, it being found that the response is very rapid indeed. The most rapid response occurs in the noon to mid-afternoon sector, where the westward flows of the dusk cell respond with a delay of 3.9 ± 2.2 min to changes in the North-South field at the subsolar magnetopause. The flows appear to evolve in form over the subsequent ~ 5 min interval, however, as indicated by the longer response times found for the northward component of flow in this sector (6.7 ±2.2 min), and in data from earlier and later local times. No evidence is found for a latitudinal gradient in response time; changes in flow take place coherently in time across the entire radar field-of-view.

Ionospheric ion upwelling in the wake of flux transfer events at the dayside magnetopause

The effects of flux transfer events (FTE) on the dayside auroral ionosphere are studied, using a simple twin-vortex model of induced ionospheric plasma flow. It is shown that the predicted and observed velocities of these flows are sufficient to drive nonthermal plasma in the F region, not only within the newly opened flux tube of the FTE, but also on the closed, or "old" open, field lines around it. In fact, with the expected poleward neutral wind, the plasma is more highly nonthermal on the flanks of, but outside, the open flux tube: EISCAT observations indicate that plasma is indeed driven into nonthermal distributions in these regions. The nonthermal plasma is thereby subject to additional upforce due to the resulting ion temperature anisotropy and transient expansion due to Joule heating and also to ion accelerations associated with the FTE field aligned current system. Any upflows produced on closed field lines in the vicinity of the FTE are effectively bunched-up in the "wake" of the FTE. Observations from the AMPTE-UKS satellite at the magnetopause reveal ion upflows of energy ∼100 eV flowing out from the ionosphere on closed field lines which are only found in the wake of the FTE. Such flows are also only found shortly after two, out of all the FTEs observed by AMPTE-UKS. The outflow from the ionosphere is two orders of magnitude greater than predicted for the "classical" polar wind. It is shown that such ionospheric ion flows are only expected in association with FTEs on the magnetopause which are well removed from the sub-solar point-either towards dusk or, as in the UKS example discussed here, towards dawn. It is suggested that such ionospheric ions will only be observed if the center of the FTE open flux tube passes very close to the satellite. Consequently, we conclude the ion upflows presented here are probably driven by the second of two possible source FTEs and are observed at the satellite with a lag after the FTE which is less than their time-of-flight.

The modelled occurrence of non-thermal plasma in the ionospheric F-region and the possible consequences for ion outflows into the magnetosphere

A global, time-dependent, three-dimensional, coupled ionosphere-thermosphere model is used to predict the spatial distribution of non-thermal plasma in the F-layer. It is shown that, even for steady-state conditions with Kp as low as 3, the difference between the ion and neutral velocities often exceeds the neutral thermal speed by a factor, D', which can be as large as 4. Theoretically, highly non-Maxwellian, and probably toroidal, ion velocity distributions are expected when D' exceeds about 1.5. The lack of response of the neutral winds to sunward ion drifts in the dawn sector of the auroral oval cause this to be the region most likely to contain toroidal distributions. The maximum in D' is found in the throat region of the convection pattern, where the strong neutral winds of the afternoon sector meet the eastward ion flows of the morning sector. These predictions are of interest, not only to radar scientists searching for non-thermal ionospheric plasma, but also as one possible explanation of the initial heating and upward flows of ions in the cleft ion fountain and nightside auroral oval, both of which are a major source of plasma for the magnetosphere.

The geomagnetic mass spectrometer— mass and energy dispersions of ionospheric ion flows into the magnetosphere

NASA's Dynamics Explorer (DE) mission was designed to study the coupling between the Earth's magnetosphere, ionosphere and neutral thermosphere1. One area of major interest is the outflow of ionospheric plasma into the magnetosphere, the scale and significance of which is only now becoming apparent with the advent of mass-resolving, low-energy ion detectors. Here we compare observations of ion flows in the polar magnetosphere, made by the retarding ion mass spectrometer (RIMS)2 on DE1, with those made simultaneously in the topside ionosphere by the ion drift meter (IDM)3 on the lower-altitude DE2 spacecraft. The results show the dayside auroral ionosphere to be a significant and highly persistent source of plasma for the magnetosphere. The upwelling ionospheric ions are spatially dispersed, according to both their energy and mass, by the combined actions of the geomagnetic field and the dawn-to-dusk convection electric field, in an effect analogous to the operation of an ion mass spectrometer.

Simplified estimation of ray-path mirroring height for HF radiowaves reflected from the ionospheric F-region

From Milsom's equations, which describe the geometry of ray-path hops reflected from the ionospheric F-layer, algorithms for the simplified estimation of mirror-reflection height are developed. These allow for hop length and the effects of variations in underlying ionisation (via the ratio of the F2- and E-layer critical frequencies) and F2-layer peak height (via the M(3000)F2-factor). Separate algorithms are presented which are applicable to a range of signal frequencies about the FOT and to propagation at the MUF. The accuracies and complexities of the algorithms are compared with those inherent in the use of a procedure based on an equation developed by Shimazaki.

Simple M-factor algorithm for improved estimation of the basic maximum usable frequency of radio waves reflected from the ionospheric F-region

The equations of Milsom are evaluated, giving the ground range and group delay of radio waves propagated via the horizontally stratified model ionosphere proposed by Bradley and Dudeney. Expressions for the ground range which allow for the effects of the underlying E- and F1-regions are used to evaluate the basic maximum usable frequency or M-factors for single F-layer hops. An algorithm for the rapid calculation of the M-factor at a given range is developed, and shown to be accurate to within 5%. The results reveal that the M(3000)F2-factor scaled from vertical-incidence ionograms using the standard URSI procedure can be up to 7.5% in error. A simple addition to the algorithm effects a correction to ionogram values to make these accurate to 0.5%.

Ionospheric origin of magnetospheric O+ ions

Flows of thermal atomic oxygen (O+) ions are deduced from topside ionospheric plasma density profiles. The mean flux within most of the polar cap is of the order of 10^12 m^{−2} s^{−1}, a figure which is consistent with both theoretical and experimental estimates of the light ion polar wind at greater altitudes. Larger flows (up to 6 × 10^12 m^{−2} s^{−1}) are observed near the poleward edge of the night-side statistical auroral oval, a feature not reproduced in the light ion flux. The implication is one of a low altitude acceleration mechanism, acting upon the O+ ions at these latitudes and at heights above that at which the fluxes are observed. Such a process would enable ions to escape from the ionosphere because they do not exchange charge with neutral hydrogen. The observations are in general agreement with energetic O+ ions as previously observed in various parts of the magnetosphere.

The bottomside mid-latitude ionospheric trough

The measured power losses and Doppler shifts of h.f. radio waves propagated over a long, west-east, sub-auroral path are found to exhibit features which cannot be explained by simple predictions and models. Both theory and the limited available data indicate that a bottomside F2-layer depletion should be present below the topside mid-latitude trough. Introducing this into the models (using the mean statistical positions of the trough deduced from Alouette I and II soundings) is shown to explain many of these features. From the Doppler shifts and a simple ray-tracing model the height of the depleted F2-peak inside the trough is deduced to be greater than its value outside the trough by an amount of the order of only 30–80 km.

The wavelet method as an alternative for reducing ionospheric effects from single-frequency GPS receivers

The ionospheric effect is one of the major errors in GPS data processing over long baselines. As a dispersive medium, it is possible to compute its influence on the GPS signal with the ionosphere-free linear combination of L1 and L2 observables, requiring dual-frequency receivers. In the case of single-frequency receivers, ionospheric effects are either neglected or reduced by using a model. In this paper, an alternative for single-frequency users is proposed. It involves multiresolution analysis (MRA) using a wavelet analysis of the double-difference observations to remove the short- and medium-scale ionosphere variations and disturbances, as well as some minor tropospheric effects. Experiments were carried out over three baseline lengths from 50 to 450 km, and the results provided by the proposed method were better than those from dual-frequency receivers. The horizontal root mean square was of about 0.28 m (1 sigma).

Application of ionospheric corrections in the equatorial region for L1 GPS users

In the absence of the selective availability, which was turned off on May 1, 2000, the ionosphere can be the largest source of error in GPS positioning and navigation. Its effects on GPS observable cause a code delays and phase advances. The magnitude of this error is affected by the local time of the day, season, solar cycle, geographical location of the receiver and Earth's magnetic field. As it is well known, the ionosphere is the main drawback for high accuracy positioning, when using single frequency receivers, either for point positioning or relative positioning of medium and long baselines. The ionosphere effects were investigated in the determination of point positioning and relative positioning using single frequency data. A model represented by a Fourier series type was implemented and the parameters were estimated from data collected at the active stations of RBMC (Brazilian Network for Continuous Monitoring of GPS satellites). The data input were the pseudorange observables filtered by the carrier phase. Quality control was implemented in order to analyse the adjustment and to validate the significance of the estimated parameters. Experiments were carried out in the equatorial region, using data collected from dual frequency receivers. In order to validate the model, the estimated values were compared with ground truth. For point and relative positioning of baselines of approximately 100 km, the values of the discrepancies indicated an error reduction better than 80% and 50% respectively, compared to the processing without the ionospheric model. These results give an indication that more research has to be done in order to provide support to the L1 GPS users in the Equatorial region.

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth's magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.