Contemporary Ionospheric Scintillation Studies: Statistics, 2D Analytical and 3D Numerical Inversion

The propagation of radiowaves through ionospheric irregularities can lead to random amplitude and phase fluctuations of the signal, otherwise known as scintillation, which can severely impact the performance of Global Navigation Satellite System (GNSS) and communication systems. Research into high latitude scintillation, through statistical analysis and inverse modeling, was completed to provide insight into the temporal and spatial distribution, and irregularity parameters, which can ultimately support the development of impact mitigation techniques, and deepen our understanding of the underlying physics. The work in this dissertation focused on the statistical analysis of Global Positioning System (GPS) scintillation data, data inversion, two-dimensional (2D) and three-dimensional (3D) scintillation modeling. The statistical analysis revealed distinct trends in the distribution of scintillation, while demonstrating that for GPS signals, phase scintillation occurs most frequently and can be treated as stochastic Total Electron Content (TEC); findings which have significant implications for impact mitigation. For the first of two inversion studies, scintillation data associated with a series of Polar Cap Patches (PCPs), which are common large-scale high latitude structures, was inverted to gain insight into the composition of the underlying irregularities. The results of this study suggest that the irregularities can be modeled as rods interbedded with sheets, which is knowledge that is crucial for the anchoring of models used to develop system mitigation techniques. The final study presents the results of modeling and inversion work to identify the conditions under which a 2D analytic version of the 3D numerical Satellite-beacon Ionospheric-scintillation global model of the upper atmosphere (SIGMA) model can be used to perform modeling in high latitude regions. During the study, it was found that the analytic model tends to diverge for electron density variance times irregularity layer thickness values exceeding 2, matched reasonably well for correlation length to thickness ratios up to 0.2, and was incompatible when ratios approached 0.35. An elevation angle limitation was also identified for the 2D model, and inflated values for the electron density variance were observed overall, which are thought to result from the weak scatter limits of the analytic model. These inflated values were particularly acute in the auroral zone during elevated conditions and suggest that the analytic model used in the study is not well suited for modeling the highly elongated irregularities associated with auroral precipitation.