We present analysis methods for studying the structuring and motion of ionospheric irregularities at the subkilometer scale sizes that produce L band scintillations. Spaced‐receiver methods are used for Global Navigation Satellite System (GNSS) receivers' phase measurements over approximately subkilometer to kilometer length baselines for the first time. The quantities estimated by these techniques are plasma drift velocity, diffraction anisotropy magnitude and orientation, and characteristic velocity. Uncertainties are quantified by ensemble simulation of noise on the phase signals carried through to the observations of the spaced‐receiver linear system. These covariances are then propagated through to uncertainties on drifts through linearization about the estimated values of the state. Five receivers of SAGA, the Scintillation Auroral Global Positioning System (GPS) Array, provide 100 Hz power and phase data for each channel at L1 frequency. The array is sited in the auroral zone at Poker Flat Research Range, Alaska. A case study of a single scintillating satellite observed by the array is used to demonstrate the spaced‐receiver and uncertainty estimation process. A second case study estimates drifts as measured by multiple scintillating channels. These scintillations are correlated with auroral activity, based on all‐sky camera images. Measurements and uncertainty estimates made over a 30 min period are compared to a collocated incoherent scatter radar and show good agreement in horizontal drift speed and direction during periods of scintillation for which the characteristic velocity is less than the drift velocity.

We introduce and analyze the first data from an array of closely spaced Global Positioning System (GPS) scintillation receivers established in the auroral zone in late 2013 to measure spatial and temporal variations in L band signals at 100–1000 m and subsecond scales. The seven receivers of the Scintillation Auroral GPS Array (SAGA) are sited at Poker Flat Research Range, Alaska. The receivers produce 100 s scintillation indices and 100 Hz carrier phase and raw in‐phase and quadrature‐phase samples. SAGA is the largest existing array with baseline lengths of the ionospheric diffractive Fresnel scale at L band. With an initial array of five receivers, we identify a period of simultaneous amplitude and phase scintillation. We compare SAGA power and phase data with collocated 630.0 nm all‐sky images of an auroral arc and incoherent scatter radar electron precipitation measurements, to illustrate how SAGA can be used in multi‐instrument observations for subkilometer‐scale studies.

We describe a method to detect and classify global positioning system (GPS) scintillation, and then hypothesize a possible ionospheric layer scattering the signal. The objective is to routinely identify events of interest to investigate in detail in future work. Scintillation types include amplitude, phase, or both amplitude and phase. A scintillation event is one for which a scintillation index remains above a threshold across a majority of closely spaced receivers viewing a single satellite. Events are categorized by signal frequency and scintillation type. An event is then hypothesized to be due to the E or F layer using an independent data source. Data from the scintillation auroral GPS array located in Poker Flat Research Range, Alaska, are used to analyze L1 and L2C frequencies in 2014 and 2015. The irregularity layer associated with each scintillation event is hypothesized to be due to the activity in the E layer of the ionosphere (below 150 km) or in the F layer (above 195 km) using collocated Poker Flat Incoherent Scatter Radar electron density measurements. Events in a transition layer (150-195 km) and inconclusive results are also recorded. We find that nearly all of the over 4,000 events are phase scintillations. The majority of the events are hypothesized to occur when the peak density is at E-layer altitudes. This indicates that E-layer-related scintillation may be quite common at auroral latitudes and that GPS receivers are sensitive to the irregularities occurring both there and at F layer altitudes.Plain Language Summary As global positioning system (GPS) signals pass through Earth's atmosphere, they may "twinkle" when received at the ground due to variations in the number of charged particles present in the ionized layer of the Earth's atmosphere, the ionosphere. The twinkling or rapid fluctuation in the GPS signals is called ionospheric scintillation. We develop a way, using data from GPS receivers and a nearby radar, to detect and classify scintillations measured in Alaska in 2014 and 2015, and then make an educated guess about whether they are in the upper or lower layer of the ionosphere. This method will be useful for handling large data sets of scintillation, which are now becoming more common. We find that, even though most of the ionosphere's charged particles are usually above 195 km height, GPS is sensitive to scintillations at high latitudes due to variations in a lower layer below 150 km, and these happen a majority of the time.

This work develops a method by which a kilometer-spaced array of Global Navigation Satellite System (GNSS) scintillation receivers can be used to estimate the ionospheric irregularity layer height and thickness and associated uncertainties on those estimates. Spectra of filtered signal power and phase data are used to estimate these quantities by comparing the observed ratio of the log of the power spectrum to the phase spectrum with the Rytov weak scatter theoretical ratio. A Monte Carlo simulation of noise on the input signal and the irregularity drift velocity is used to quantify the error in estimates of height and thickness. The method is tested using data from the Scintillation Auroral Global Positioning System (GPS) Array (SAGA) sited in the auroral zone at Poker Flat Research Range, Alaska. For the 30-min scintillation period studied, the technique identifies ionospheric scattering from a thick F layer, which correlates well with on-site incoherent scatter radar measurements of peak electron density, for an event previously identified in the literature as likely due to F layer.