Single-frequency Global Positioning System (GPS) receivers do not accurately compensate for the ionospheric delays imposed upon GPS signals. This can lead to significant errors and single-frequency systems rely upon models to compensate. This investigation applies 4D (four-dimensional) ionospheric tomography to GPS timing for the first time. The tomographic algorithm, MIDAS (Multi-Instrument Data Analysis System), is used to correct for the ionospheric delay and the results are compared to existing single and dual-frequency techniques. Days during the solar maximum years 2002, 2003 and 2004 have been chosen to display results when the ionospheric delays are large and variable. Maps of the ionospheric electron density, across Europe, are produced by using data collected from a fixed network of dual-frequency GPS receivers. Results that improve upon the use of existing ionospheric models are achieved for fixed (static) and mobile (moving) GPS receiver scenarios. The effects of excluding all of the GPS satellites below various elevation masks, ranging from 5° to 40°, on timing solutions for fixed and mobile situations are also presented. The greatest timing accuracies when using the fixed GPS receiver technique are obtained by using the highest mask. The mobile GPS timing solutions are most accurate when satellites at lower elevations continue to be included. Furthermore, timing comparisons are made across baselines up to ~4000 km and the ionospheric errors are shown to increase with increasing baseline. GPS time transfer is then investigated and MIDAS is shown to improve the time transfer stabilities of a single-frequency GPS system. The results are comparable to the dual-frequency time transfer after ~2 hours averaging time. Overall, the MIDAS technique provides the most accurate and most stable results (comparable to dual-frequency) for a single-frequency based GPS system. Ionospheric corrections (via MIDAS) may be broadcast to users nationally or via the internet for example, opening up the possibility of improving the accuracy and stability of single-frequency GPS systems in real-time.

Global Navigation Satellite Systems can provide position, velocity and time information to users using receiver hardware. The United States developed Global Positioning System (GPS) is the only current fully operational system; however further systems are in development. The GPS has shown considerable success for navigation, but it still has a number of problems that limit its accuracy. The two main problems are the ionosphere and local environment of the receiver. The ionosphere causes a delay and random rapid shifts in phase and amplitude (scintillation) to the signal. The local environment can provide the signal with multiple routes (multi-path) to the receiver. In this project a GPS signal simulator is developed, which models the effects of the ionosphere and multi-path on the modulated signals. The focus is made on the GPS system as the simulator measurements can be compared to the real measurements; however other systems will be considered in the future. A number of experiments investigating multi-path and ionospheric effects on a receiver’s ability to track the signals have been completed. The simulator has been used to replicate a real local multi-path environment and the results have been compared. Further investigations of the multi-path have shown a unique multi-path signature in the receiver power output. The later part of the thesis describes a case study investigating a short but rapid period of scintillation observed on three receivers based in Norway. An analysis of the multi-path environment was completed, but was found not to be the cause. The ionosphere was investigated using equipment based across Scandinavia. The equipment showed that geomagnetic conditions were disturbed at the time of the event. The GPS measurements were compared with all-sky camera data to show that the scintillation can be attributed to the GPS signal path crossing electron density structures associated with the aurora.

This thesis aims to answer the question, “What physical process dominates the formation of plasma irregularities, capable of directly or indirectly causing GPS L1 band scintillation, in polar cap plasma patches during magnetic storm conditions?.” A novel modelling technique utilising an ionospheric imaging algorithm is developed and used to elucidate the relative importance of the two most commonly discussed processes. These are the Gradient Drift Instability (GDI) and turbulence induced by electric field mapping to the ionosphere from the magnetosphere. The results show that in magnetic storm conditions, at times the GDI process is dominant, but that at other times turbulence may be as significant as the GDI in determining how the plasma within a polar cap patch behaves, possibly more so. This in turn suggests that further study of the turbulence process is necessary in order to fully understand how big a role it plays in causing GPS L1 band scintillation in the polar cap. The success of the modelling technique developed here shows the utility of ionospheric imaging as a tool for understanding physical problems of the ionosphere; efforts to improve it and to apply it in other contexts would be worthwhile.

Electric probes have been the primary instruments for the in situ investigation of plasma parameters in the Earth’s ionosphere. This dissertation is a compendium of three papers, each dealing with a separate spacecraft that carried one or more instruments based on the electric probe technique. The first paper presents data from the Sudden Atom Layer sounding rocket that carried an RF Impedance Probe, a DC fixed-bias Langmuir Probe (DCP), and an Electric Field Probe. The combined dataset indicates a case of payload surface charging, the causes of which are investigated within the paper. A generic circuit model is developed to analyze payload charging and behavior of Langmuir-type instruments. Our analysis indicates that the anomalous charging event was an outcome of triboelectrification of the payload surface from neutral dust particles present in the Earth’s mesosphere. These results suggest caution in interpreting observations from the Langmuir class of instrumentation within dusty environments. The second paper presents data from the Floating Potential Measurement Unit (FPMU) that is deployed on the International Space Station. The FPMU instrument suite consists of three different Langmuir-type probes and a Plasma Impedance Probe (PIP). We first give a brief overview of the instrumentation, and then describe the algorithm used to reduce Langmuir probe I-V curves to plasma parameters. It is shown that the derived temperatures agree well with International Reference Ionosphere (IRI) model, while the derived density matches better with the USU-Global Assimilation of Ionospheric Measurement model. The third paper presents the dataset from the EQUIS II sounding rocket campaign. The rocket payloads carried a PIP, a DCP, and an internally heated Sweeping Langmuir Probe. The ratio of the payload surface area to the cumulative area of the instrument and its guard was about 250. We show that on small sounding rocket payloads the DCP technique of relative electron density measurement is not very accurate. We further show that the ion saturation region analysis of the I-V curve produces absolute ion density that matches very well with the absolute electron density derived from the PIP, and the derived temperatures agree reasonably well with the IRI model.

Aurorae often break down into elongated filaments parallel to the geomagnetic field lines (B) with cylindrically symmetric structures. The object of this thesis is to study the ion distribution function and transport properties in response to the sudden introduction of a radial electric field (E) in such a cylindrical geometry. Both collision-free and collisional situations are considered.<p> The thesis starts by solving a collision-free problem where the electric field is constant in time but increases linearly with radius, while the initial ion density is uniform in space. The attendant Boltzmann equation is solved by tracking the ions back in time, thereby using the temporal link between the initial position and velocity of an ion and its position and velocity at an arbitrary time and place. Complete analytical solutions show that the ion distribution function is a pulsating Maxwellian in time, and all transport parameters (e.g., bulk speed, temperature, etc.) oscillate in time but independent of radius. If the ion-neutral collisions are taken into account by employing a simple relaxation model, analytical solutions are also obtained. In this case, the ion distribution function can be driven to horseshoe shapes which are symmetric with respect to the ExB direction. The bulk parameters evolve in a transition period of the order of one collision time as they go from oscillating to the non-oscillating steady state.<p> In more realistic electric field structures which are spatially inhomogeneous but still constant in time, a generalized semi-numerical code is developed under collision-free conditions. This code uses a backmapping approach to calculate the ion velocity distribution and bulk parameters. With arbitrarily selected electric field rofiles, calculations reveal various shapes of ion velocity distribution functions (e.g., tear-drop, core-halo, ear-donut, etc). The associated transport properties are also obtained and discussed.<p> Under both collision-free and collisional conditions, the effect of the density inhomogeneities at the initial time is studied in an electric field which is proportional to radius and constant in time. With two profiles of the initial ion density for the collision-free case, and one profile for the collisional case, complete analytical solutions are obtained. The results reveal that the distribution function and the bulk properties are now strongly dependent on radial position.<p> If the radial electric field is unable to stay constant with time but modulated by in-coming charged particles, a fluid formalism is used to study the excitation of several plasma waves under different kinds of initial conditions. These identified waves include the ion cyclotron oscillation, the ion and electron upper-hybrid oscillations, and the lower-hybrid oscillation.<p> The results of this thesis are expected to be applicable to high-resolution observations. Future work should also include the mirror effect and the formation of conics in velocity space. Finally, the velocity distributions obtained in this thesis could trigger various plasma instabilities, and this topic should also be looked at in the future.

In 2011, the Canadian enhanced Polar Outflow Probe (ePOP) satellite will be launched. The ePOP satellite is equipped with several scientific Earth observation instruments, including a Radio Receiver Instrument (RRI) which will be used to detect High Frequency (HF) radio waves transmitted from a ground-based transmitter. The ground-based instrument will be one of the Super Dual Auroral Radar Network (SuperDARN) array of radars. A radio wave transmitted from the SuperDARN radar will propagate through the ionosphere and be detected by the RRI on ePOP. Analysis of the characteristics of the signal received by the RRI will provide information about the plasma density in the ionosphere between the transmitter and receiver. As the ePOP satellite is not yet operational, extensive ray path modelling has been performed to simulate the expected signal at the RRI for various ionospheric conditions.<p> The other major objective of this research was to examine the effect of the variable refractive index in the ionosphere on SuperDARN drift velocity measurements. Past comparisons between velocities measured by SuperDARN and other instruments have found that velocities measured by SuperDARN typically were about 20-30% lower. This research has shown that underestimation of drift velocities by SuperDARN is a consequence of not including the refractive index when these velocities are calculated. As refractive index measurements are not readily available, this research has involved developing and implementing various methods to estimate the refractive index in the ionosphere. These methods have demonstrated that plasma density values within the SuperDARN scattering volume are appreciably higher than background plasma densities in the ionosphere. Application of these methods, which has resulted in a much better understanding of the physics of the coherent scattering process, has resulted in agreement between velocities measured by SuperDARN and other instruments.

Aurorae often break down into elongated filaments parallel to the geomagnetic field lines (B) with cylindrically symmetric structures. The object of this thesis is to study the ion distribution function and transport properties in response to the sudden introduction of a radial electric field (E) in such a cylindrical geometry. Both collision-free and collisional situations are considered.<p> The thesis starts by solving a collision-free problem where the electric field is constant in time but increases linearly with radius, while the initial ion density is uniform in space. The attendant Boltzmann equation is solved by tracking the ions back in time, thereby using the temporal link between the initial position and velocity of an ion and its position and velocity at an arbitrary time and place. Complete analytical solutions show that the ion distribution function is a pulsating Maxwellian in time, and all transport parameters (e.g., bulk speed, temperature, etc.) oscillate in time but independent of radius. If the ion-neutral collisions are taken into account by employing a simple relaxation model, analytical solutions are also obtained. In this case, the ion distribution function can be driven to horseshoe shapes which are symmetric with respect to the ExB direction. The bulk parameters evolve in a transition period of the order of one collision time as they go from oscillating to the non-oscillating steady state.<p> In more realistic electric field structures which are spatially inhomogeneous but still constant in time, a generalized semi-numerical code is developed under collision-free conditions. This code uses a backmapping approach to calculate the ion velocity distribution and bulk parameters. With arbitrarily selected electric field rofiles, calculations reveal various shapes of ion velocity distribution functions (e.g., tear-drop, core-halo, ear-donut, etc). The associated transport properties are also obtained and discussed.<p> Under both collision-free and collisional conditions, the effect of the density inhomogeneities at the initial time is studied in an electric field which is proportional to radius and constant in time. With two profiles of the initial ion density for the collision-free case, and one profile for the collisional case, complete analytical solutions are obtained. The results reveal that the distribution function and the bulk properties are now strongly dependent on radial position.<p> If the radial electric field is unable to stay constant with time but modulated by in-coming charged particles, a fluid formalism is used to study the excitation of several plasma waves under different kinds of initial conditions. These identified waves include the ion cyclotron oscillation, the ion and electron upper-hybrid oscillations, and the lower-hybrid oscillation.<p> The results of this thesis are expected to be applicable to high-resolution observations. Future work should also include the mirror effect and the formation of conics in velocity space. Finally, the velocity distributions obtained in this thesis could trigger various plasma instabilities, and this topic should also be looked at in the future.

In 2011, the Canadian enhanced Polar Outflow Probe (ePOP) satellite will be launched. The ePOP satellite is equipped with several scientific Earth observation instruments, including a Radio Receiver Instrument (RRI) which will be used to detect High Frequency (HF) radio waves transmitted from a ground-based transmitter. The ground-based instrument will be one of the Super Dual Auroral Radar Network (SuperDARN) array of radars. A radio wave transmitted from the SuperDARN radar will propagate through the ionosphere and be detected by the RRI on ePOP. Analysis of the characteristics of the signal received by the RRI will provide information about the plasma density in the ionosphere between the transmitter and receiver. As the ePOP satellite is not yet operational, extensive ray path modelling has been performed to simulate the expected signal at the RRI for various ionospheric conditions.<p> The other major objective of this research was to examine the effect of the variable refractive index in the ionosphere on SuperDARN drift velocity measurements. Past comparisons between velocities measured by SuperDARN and other instruments have found that velocities measured by SuperDARN typically were about 20-30% lower. This research has shown that underestimation of drift velocities by SuperDARN is a consequence of not including the refractive index when these velocities are calculated. As refractive index measurements are not readily available, this research has involved developing and implementing various methods to estimate the refractive index in the ionosphere. These methods have demonstrated that plasma density values within the SuperDARN scattering volume are appreciably higher than background plasma densities in the ionosphere. Application of these methods, which has resulted in a much better understanding of the physics of the coherent scattering process, has resulted in agreement between velocities measured by SuperDARN and other instruments.

A range of shortwave frequencies will be returned to earth by the ionosphere such that long-distance shortwave communication is accomplished. Because of ionospheric variations, operation parameters are not the same in different time duration. We must select optimized parameters with the aim of utilizing the ionosphere effectively for shortwave radio communications. Shortwave propagation prediction programs can provide such information. Empirical models are often used to represent the variations of the electron density with height in the ionosphere, since radio propagations are concerned with electron density. We develop a graphical user interface based on VOACAP that is one of the most widely used prediction programs. We also improve the electron density model in VOACAP. A gradient discontinuity occurs in the original parabolic ionospheric model. We incorporate an electron density model that does not contain any gradient discontinuities. In this model the electron density profile is composed of three segments: (1) a cosine F2 layer; (2) an E-F layer comprising a secant function, and optional cosine function to provide an F1 ledge (3) a parabolic E layer. Compared with experimental data, our modified model results in more accurate predictions.

Thesis (M. Sc.)--University of Hong Kong, 1968. / Mimeographed.