PRIKRYL_#6219_OK_corrected_Layout 6 ANNALS OF GEOPHYSICS, 56, 2, 2013, R0222; doi:10.4401/ag-6219 R0222 Probabilistic forecasting of ionospheric scintillation and GNSS receiver signal tracking performance at high latitudes Paul Prikryl1,*, Veettil Sreeja2, Marcio Aquino2, Periyadan T. Jayachandran3 1 Communications Research Centre, Ottawa, ON, Canada 2 University of Nottingham, Nottingham Geospatial Institute, Nottingham, United Kingdom 3 University of New Brunswick, Physics Department, Fredericton, NB, Canada ABSTRACT At high latitudes, phase scintillation occurs predominantly on the day- side in the ionospheric footprint of the magnetospheric cusp, and in the nightside auroral oval. A new technique of probabilistic forecasting of phase scintillation occurrence relative to the arrival time of high-speed solar wind from coronal holes and interplanetary coronal mass ejections has recently been proposed [Prikryl et al. 2012]. Cumulative probability distribution functions for the phase-scintillation occurrence that are ob- tained can be specified for low (below-median) and high (above-median) values of various solar wind plasma parameters. Recent advances in mod- eling of high-speed solar wind and coronal mass ejections, combined with the probabilistic forecasting of scintillation, will lead to improved opera- tional space weather forecasting applications. Scintillation forecasting and mitigation techniques need to be developed to avoid potential costly failures of technology-based Global Navigation Satellite Systems in the near future, in particular during the upcoming solar maximum. The Global Navigation Satellite Systems receiver-tracking performance during severe scintillation conditions can be assessed by the analysis of receiver phase-locked-loop jitter. Tracking jitter maps [Sreeja et al. 2011] offers a potentially useful tool to provide users with expected tracking conditions, if based on scintillation predictions as proposed above. Scintillation in- dices are obtained from L1 GPS data collected with the Canadian High Arctic Ionospheric Network. Combined with high rate amplitude and phase data, they can be used as input to receiver tracking models to de- velop scintillation mitigation techniques. 1. Introduction Global Navigation Satellite Systems (GNSS), such as the US Global Positioning System (GPS), the Russian GNSS (GLONASS) and the EU Galileo, have become indispensable for precise positioning and time keeping for the support of a wide range of civilian and military applications. The main threat to reliable and safe oper- ation of modern technology that relies on GNSS is the variable propagation conditions that are encountered by radio signals as they pass through the ionosphere. Variable ionospheric delays and signal scintillation can affect the performance of radio communications and navigation systems, and can severely degrade GPS po- sitional accuracy. Scintillation, rapid signal amplitude fading, and phase fluctuation, are phenomena of great concern at low and high latitudes. As pointed out by Conker et al. [2003], to determine the effects of scintil- lation on GNSS receivers, including the GPS/Satellite- Based Augmentation System (SBAS), "both predictions of scintillation activity in the ionosphere and residual effects of this activity after processing by a receiver" must be considered. These authors estimated the scin- tillation effects on the availability of GPS and SBAS for L1 (1575.42° MHz) C/A and L2 (1227.60° MHz) semi- codeless receivers. In the present study, the benefits of a probabilistic forecasting of phase scintillation to ap- plications relying on robust and accurate phase signal tracking are addressed. Ionospheric scintillation, which is seen as fluctua- tions of radio signal amplitude and phase that are caused by ionospheric irregularities [Aarons 1982, Basu et al. 1998], is most intense at low latitudes during post sunset hours [Basu et al. 2002], and at high latitudes in the nightside auroral oval and on the dayside in the re- gion called the cusp [Basu et al. 1998, Spogli et al. 2009, Prikryl et al. 2011a]. At high latitudes, ionospheric scin- tillation is largely controlled by the process of coupling between the solar wind, including the interplanetary magnetic field, and the Earth magnetic field. In the nightside auroral oval, intense electric currents (auroral electrojets), various plasma instabilities, and energetic particle precipitation result in steep plasma density gra- dients and irregularities. It has been shown that scintil- lation strongly affects GPS receiver performance in the Article history Received December 20, 2011; accepted January 13, 2012. Subject classification: Scintillation, Space weather forecasting, GNSS tracking performance. auroral zone [e.g., Skone and de Jong 2000]. In the cusp, just poleward of the dayside auroral oval, where intense ionospheric convection and forma- tion of patches of enhanced ionization occur, scintilla- tion occurrence is particularly enhanced. Scintillation occurrence is more seasonally dependent in the cusp than in the nightside auroral oval [Prikryl et al. 2011a], which is characterized by annual and semi-annual vari- ations, respectively. It is also strongly enhanced during intense solar wind events, as either solar wind co-rotat- ing interaction regions (CIRs), or interplanetary coro- nal mass ejections (ICMEs). ICMEs and CIRs on the leading edge of high-speed streams (HSSs) from coro- nal holes are major causes of ionospheric disturbances. The strongest of these solar wind events result in major geomagnetic storms [e.g., Borovsky and Denton 2006, Richardson et al. 2006] that have a global impact on the ionosphere. Based on a statistical study of scintillation occurrence during such events from the years 2008 to 2010, a new method of probabilistic forecasting of phase scintillation was proposed [Prikryl et al. 2012]. A summary of this method using a larger number of events is presented in the present study. Scintillation indices, S4 and vU, are widely used as the measures of amplitude and phase scintillation, re- spectively. The S4 index is the standard deviation of the received power normalized by its mean value, and the phase scintillation index vU is the standard deviation of the detrended phase using a filter in the receiver with a 0.1 Hz cut-off frequency. However, as pointed out by Aquino et al. [2007], these scintillation indices do not fully characterize the actual instantaneous values of phase and amplitude fluctuations that affect GNSS re- ceiver performance. They combined scintillation in- dices and high rate GNSS phase and amplitude data with state-of-the-art receiver tracking models, to study receiver tracking performance under scintillation con- ditions. This strategy was also applied for an intense scintillation event at low latitudes [Sreeja et al. 2011], where the variance of the phase-locked-loop (PLL) tracking error (tracking jitter) of a GPS receiver was used to construct maps that can provide an indication of tracking conditions for that particular receiver. In this study, tracking-jitter maps are obtained for a high- latitude scintillation event. 2. The Canadian High Arctic Ionospheric Network At present, the Canadian High Arctic Ionospheric Network (CHAIN) [Jayachandran et al. 2009] consists of ten specialized GPS receivers and six Canadian Ad- vanced Digital Ionosondes that are distributed in the auroral oval and the polar cap (Figure 1). The latitudes and longitudes of the stations can be found elsewhere [e.g., Prikryl et al. 2011a, their Table 1]. The receivers are GPS ionospheric scintillation and total electron content (TEC) monitors (GISTMs) GSV 4004B [Van Dierendonck and Arbesser-Rastburg 2004]. The GISTMs consist of a NovAtel OEM4 dual frequency re- ceiver with special firmware specifically configured to measure and log the power and phase of the GPS L1 signal at a high sampling rate (50 Hz). The receiver computes ionospheric TEC from the GPS L1 and L2 signals. The GSV 4004B can also automatically com- pute and record S4 and vU indices; 1-min values are used in the present study. Each receiver can track and report scintillation and TEC measurements simultane- ously from up to 10 GPS satellites in view. 3. Occurrence of GPS phase scintillation at high lat- itudes A scintillation climatology study for the solar min- imum years of 2008 and 2009 showed that phase scin- tillation, as a function of magnetic local time and geomagnetic latitude, primarily occurs in the night-side auroral oval and on the dayside in the ionospheric foot- print of the cusp [Prikryl et al. 2011a]. As expected, these scintillation regions shift in latitude in response to varying geomagnetic activity. We have extended the climatology data set to 4 years. Figure 2a, b shows the maps of the percentage occurrence of phase scintilla- tion with vU that exceeds 0.1 radians as a function of magnetic local time (MLT) and corrected geomagnetic (CGM) latitude for 1,272 quiet and 189 disturbed days selected using the threshold of 3-h occurrence of Kp >2 PRIKRYL ET AL. 2 Figure 1. Canadian High Arctic Ionospheric Network (CHAIN). The loci of IPPs at 350 km altitude for all GPS satellites, as seen from Iqaluit at elevation angles above 30° during one day. The CGM latitudes 70° and 80° (yellow) are superposed over the geographic grid in the dotted line. 3 for more than 60% of the day. The boundaries of the Feldstein statistical auroral oval for quiet and moder- ately disturbed conditions are superposed on the maps [Holtzworth and Meng 1975], indicated by the value of the IQ parameter. A dotted concentric circle in Figure 2b shows the latitude of Cambridge Bay. The Cam- bridge Bay dataset samples the poleward edge of the auroral oval in the night [Prikryl et al. 2010], but it is representative of the cusp during the day. The Cam- bridge Bay data time series are used to obtain cumula- tive probability distribution functions for the phase scintillation occurrence for the days before and after the arrivals of CIR/HSSs and ICMEs. Figure 3a, b shows scintillation occurrence maps for the first two days of HSS/CIRs and for ICME days, respectively. In total, 88 HSS/CIRs with maximum solar wind velocity exceeding 500 km/s are included. The updated catalog of near-Earth ICMEs [Richardson and Cane 2010] includes 57 ICME events from 2008 to 2011 (101 days in total). However, most of these ICMEs did not cause significant geomagnetic storms [e.g., Prikryl et al. 2011b]. Scintillation occurred predominantly in the cusp for both datasets, but the occurrence was higher for ICME days than for HSS/CIR days. This is opposite to the 2008 to 2010 results [Prikryl et al. 2012] and it is a sign that many more significant ICMEs oc- curred in 2011. Also, comparing the total numbers of HSS/CIRs and ICMEs, they increased by 114% and 33%, respectively, from 2010 to 2011. 4. The superposed-epoch analysis of geomagnetic and phase scintillation indices Arrivals of HSSs originating from coronal holes are followed by significantly enhanced phase scintillation occurrence, particularly in the cusp, which usually con- tinues for a few days. To obtain the mean response of the scintillation relative to a new CIR arrival time, a su- perposed-epoch analysis on time series of scintillation occurrence (Figure 4a) and solar wind parameters was applied, using the key times of 88 CIRs. Figure 4a shows that the mean phase scintillation occurrence steeply rises at the CIR arrival day and tapers off a few days later, which is similar to the mean Kp index and the quasi-invariant index, QI [Osherovich et al. 1999]. Also, the average number of GPS cycle slips, which is nearly zero a few days before the CIR key day, increases at the PROBABILISTIC FORECASTING OF SCINTILLATION Figure 2. Phase scintillation occurrence maps for geomagnetically quiet (a) and moderately disturbed (b) days. Figure 3. Phase scintillation occurrence maps for days of CIRs/HSSs (a) and ICMEs (b). PRIKRYL ET AL. 4 CIR interface and remains elevated during the HSS. The one-day periodicity in scintillation is due to the Cam- bridge Bay cusp location and is only partly smoothed out in this relatively small statistical sample of 88 CIRs. Figure 4b shows the results of the superposed- epoch analysis for 39 ICMEs for which the Dst index dropped below −10 nT, eliminating the weakest ICMEs. The QI increases more gradually but to higher values than for CIRs. This is because many ICMEs are actually magnetic clouds, and the ratio of the magnetic pressure to ram pressure that the QI represents is higher inside magnetic clouds. The phase scintillation occur- rence peaks at ~15% and the mean number of cycle slips peaks at more than 5 /h on the day of the ICME arrival. A cycle slip [Horvath and Crozier 2007] is de- fined here as a jump in differential phase TEC of more than or equal to 1.5 TEC Units in 1 s (1 TEC Unit cor- responds to 1016 electrons/m2). If the superposed- epoch analysis is further limited to only 15 of the strongest ICMEs (Dst < −30 nT), the peak mean scin- tillation occurrence, the mean number of cycle slips, and the Kp index are 22%, 15 and 4.3, respectively. 5. Probabilistic forecasting of phase scintillation in the cusp based on the superposed-epoch analysis of solar wind CIRs and ICMEs McPherron and Siscoe [2004] and McPherron and Weygand [2006] demonstrated a new promising tech- nique to forecast geomagnetic disturbance indices rel- ative to the stream interface between the fast and slow solar wind. In particular, they showed geoeffectiveness of what they called co-rotating compression ridges, their reference to CIRs. This method is readily applica- ble to other manifestations of space weather in the geo- space environment that are sensitive to arrivals of HSSs, including ionospheric scintillation. Figure 5a shows the cumulative probability distri- bution function (PDF) for the phase scintillation oc- currence at Cambridge Bay for the days before and after the major CIR/HSSs arrivals. The probability that phase scintillation will exceed a given value is shown, as plotted on the abscissa. The PDFs are stepped by 2-day intervals from 6 days before and 6 days after CIRs. Fig- ure 5b shows the PDFs for days before and after ICMEs. The probabilities are clearly highest for the first two days after CIR and ICME arrivals. These PDFs can be used to develop high-latitude scintillation forecasting techniques [Prikryl et al. 2012]; e.g., scintillation occurrence in the cusp relative to the stream interface between the fast and slow solar wind and/or shocks associated with CIRs and ICMEs. To gain further practical application of the probabilistic forecasting of scintillation occurrence, tracking-jitter maps can be constructed to assess GNSS receiver sig- nal-tracking performance under anticipated scintilla- tion conditions [Sreeja et al. 2011]. 6. Tracking-jitter maps In a GNSS receiver, the PLL aims to minimize the error between the input phase and its estimated phase output. It is the magnitude of this error that determines the ability of the loop to remain locked. The variance of the error at the output of the PLL (the tracking-jit- ter variance) increases during scintillation, and hence is a good measure of the effects of scintillation on the re- ceiver. The tracking-jitter variance at the output of the PLL can be calculated using the scintillation sensitive tracking model of Conker et al. [2003]. The formula for the GPS L1 carrier PLL accounts for the effects of scin- tillation on the input phase and computes the tracking- jitter variance , in radians2, as [Conker et al. 2003]: (1) where are, respectively, the error vari- ance components relating to the phase scintillation, the Figure 4. Superposed epoch analysis of QI, Kp, number of cycle slips Ncys and phase scintillation occurrence. The time series are keyed by the start times of CIRs (a) and ICMEs (b). 2vz^ h , andS T OSC 2 2 2v v v 2 2 2 2 S T OSC v v v v= + +z z z z 5 PROBABILISTIC FORECASTING OF SCINTILLATION thermal noise (amplitude scintillation), and the oscilla- tor noise (assumed to be 0.01 radians2 in the receiver). In Equation (1), the amplitude scintillation is mod- eled as an increase in the thermal noise, related to the decrease in the received signal power, in radians2, as [Conker et al. 2003]: (2) where Bn is the L1 third-order PLL one-sided bandwidth (equal to 10 Hz in this study), (c/n0)L1−C/A is the frac- tional form of the signal-to-noise density ratio, equal to 100.1C/N0, h is the predetection integration time (equal to 0.02s in this study), and S4(L1) is the ampli- tude scintillation index at L1 frequency. The phase scintillation component is modeled, in radians2, as [Conker et al. 2003]: (3) where T is the spectral strength of the phase power spectrum density (PDS) at 1 Hz, p is the spectral slope of the phase PSD, k is the order of the PLL (equal to 3 in this study), and fn is the loop natural frequency (equal to 1.91 Hz). For this study, the one minute scintillation indices and the high data rate (50-Hz phase and amplitude fluc- tuations) for the GPS L1 C/A signal recorded by the re- ceivers in the CHAIN network on November 1, 2011, were considered. A satellite elevation angle cut-off of 30° was applied to the data to reduce the impact of non-scintillation-related effects on the tracking errors (e.g., induced by multipath). The latitude of the ionos- pheric pierce point (IPP) for the different satellite-to-re- ceiver links was calculated every 1 min, assuming a single-shell ionospheric model at an altitude of 350 km. The PLL jitter for the different satellites in view with elevation angles greater than 30° at each epoch were calculated using the equations above. These PLL jitter relate to the PLL tracking error assumed in the slant di- rection of the line of sight between receiver and satel- lite. For practical use of the proposed PLL jitter maps it is necessary that a verticalised version of this error is represented, so that the user can estimate the error for any specific line of sight by applying an appropriate mapping function. This mapping function, as explained in Sreeja et al. [2011], is the subject of ongoing research funded by the UK Engineering and Physical Sciences Research Council (EPSRC). For the construction of the map presented in the present study, however, an ap- proximation was used to convert the slant PLL jitter to vertical PLL jitter, by assuming a standard mapping function cos |, where | is the zenith angle at the IPP. The values of the verticalised PLL jitter were then grid- ded in bins with a resolution of 1 min in time and 0.5° in IPP latitude, to produce the PLL jitter map of Fig- ure 7 in the next section. Figure 5. Cumulative probability distribution functions for the phase scintillation occurrence for Cambridge Bay, stepped by 2-day intervals before and after CIRs (n = 88) (a) and ICMEs (Dst < −10 nT; n = 39) (b). B / 2 / c n S L c n S L 1 2 1 1 1 1 2 1 / / L C A n L C A2 0 1 4 2 0 1 4 2 T v h = - + - z - - ^ ^^ ^ ^^ h hh h hh; E Tr k k p 2 2 1 n p 2 1 S v r = + -z - sinkf 6 @c m PRIKRYL ET AL. 6 7. The event of November 1, 2011 Figure 6 shows 5-min averages of the solar wind velocity, Vsw, the interplanetary magnetic field (IMF) BZ component in geocentric-solar-magnetospheric (GSM) coordinates, and the magnetic field magnitude, BT, from the OMNI dataset projected to subsolar bow shock (http://omniweb.gsfc.nasa.gov/). The arrival of an ICME on November 1 (shown by a vertical dotted line) resulted in an AE index of about 1000 nT and the Dst index dipped to −61 nT, intensifying a geomagnetic storm that started earlier after the arrival of another ICME and a moderate HSS stream on October 30. The hourly percent occurrence of phase scintillation index vU >0.1 rad was enhanced between 10:00 and 22:00 UT, as shown for the Cambridge Bay GPS receiver. Figure 7 shows the maps of vU values (top panel) and PLL jitter (bottom panel) as a function of UT and IPP geographic latitude for the CHAIN GPS L1 signals, as described in the previous section. These maps were constructed using the data from all of the stations of Figure 1, except station Qikiqtarjuaq. It is also worth noting that for the construction of the maps in Figure 7, not only the PLL jitter have been verticalised, but also the vU values. The latter were verticalised using the mapping function described in Spogli et al. [2009]. It can be observed from Figure 7 that the regions of en- hanced PLL jitter generally coincide with enhance- ments in vU. The phase scintillation events at high latitudes occurred primarily in the cusp and dayside polar cap between ~10:00 and 20:00 UT, which also co- incided with the occurrence of PLL jitter enhance- ments. This indicates increased likelihood for the occurrence of cycle slips and loss of lock, which de- grade the positioning accuracy. In Figure 8 the percentage occurrence of scintilla- tion and PLL jitter are mapped as functions of mag- netic latitude and MLT, to show the main regions of enhanced scintillation relative to the auroral oval. Fig- ure 8 shows the occurrence maps for S4 and the vU in- dices, as well as for PLL jitter, on a coarser grid Figure 6. Solar wind magnetic field and plasma parameters from OMNIWeb dataset combining available solar wind monitor data projected to the nose of the Earth bow shock. The hourly percent occurrence of vU > 0.1 rad from Cambridge Bay is highlighted in blue. Figure 7. Phase scintillation index vU (a) and receiver PLL jitter maps (b) as functions of UT and IPP geographic latitude for the CHAIN GPS L1 signals. 7 PROBABILISTIC FORECASTING OF SCINTILLATION (compared to Figure 7) of 1° × 0.25 h of MLT. The po- sition of the statistical auroral oval for very disturbed conditions (IQ = 6) [Feldstein 1967, Holzworth and Meng 1975] is superposed. The phase, and to a lesser degree the amplitude, scintillation occurrence exceeds given threshold values maximized in the cusp and the polar cap, although significant scintillation also oc- curred at auroral and even at subauroral latitudes. The percentage occurrence of PLL jitter exceeding 0.15 rad was significantly lower, although this generally mapped to the same regions as the phase scintillation; i.e., the cusp, the polar cap, and the auroral oval. Figure 9 shows the 60-s values of vU and PLL jitter obtained for the CHAIN GPS receiver in Iqaluit on No- vember 1, 2011. The scatter plot (Figure 9c) shows that vU is moderately correlated with the PLL jitter (Pear- son's correlation coefficient, 0.68). 8. Discussion The climatology of GPS phase scintillation at high latitudes [Prikryl et al. 2011a] has been extended by two more years to a period from 2008 to 2011. This larger dataset has confirmed that scintillation predom- inantly occurs in the ionospheric footprint of the cusp. The primary driver of the cusp dynamics is the solar wind coupling to the dayside magnetopause, which peaks with arrivals of ICMEs and HSS/CIRs. Based on the statistics of these solar wind events, using a super- posed epoch analysis technique, a probabilistic fore- casting method that focuses on the phase scintillation in the cusp has been proposed [Prikryl et al. 2012]. Solar activity has been gradually rising since 2010, and there was a significant increase in the number of mod- erate ICMEs in addition to HSS/CIRs. The number of ICMEs that resulted in Dst < −10 nT doubled in 2011 alone, when compared with the period of extended solar minimum (2008-2010). The superposed-epoch analysis based on improved statistics of events shows consistently higher probability of scintillation occur- rence in the cusp on the days of arrival of ICMEs and HSS/CIRs, in contrast to the low probability 1-2 days earlier. This pattern of significant increase in scintilla- tion occurrence makes it possible to specify phase scin- tillation forecasts for anticipated arrival days of ICMEs and HSS/CIRs [Prikryl et al. 2012]. While the scintillation activity in the polar cap was generally quite low from 2008 to 2010, during moderate to strong events, irregularities in a tongue of ionization or in polar patches that were drawn into the polar cap caused significantly enhanced scintillation occurrence in polar regions [Prikryl et al. 2011b, Prikryl et al. 2013a, this issue]. This depended on the phase of the season in either the southern or northern hemisphere. In partic- ular, the event discussed in Section 7 was associated with strong ionospheric convection across the polar cap and enhanced TEC in a tongue of ionization [Prikryl et al. 2013b]. Fast drifting ionospheric irregularities within a structured tongue of ionization resulted in high occur- rence of strong phase scintillation (Figure 8b). Under scintillation, maps of tracking jitter [Sreeja et al 2011] can be used to assess the tracking conditions, the likelihood for the occurrence of cycle slips, and the loss of lock, which are known to degrade the GNSS po- sitioning accuracy. The PLL jitter is correlated with vU and the enhancements of PLL jitter appear to coincide, in general, with the occurrence of phase scintillation. Figure 8. The occurrence of S4 index (a), phase scintillation (b) and receiver PLL jitter (c) as functions of MLT and CGM latitude for November 1, 2011. The statistical auroral oval for very disturbed conditions is superposed. PRIKRYL ET AL. 8 Hence, applications relying on robust and accurate phase signal tracking should benefit from the proba- bilistic forecasting of phase scintillation. 9. Conclusions Phase scintillation occurrence maps from 2008 to 2011 for geomagnetically quiet and disturbed days, and more specifically for the prime days of HSS/CIRs or ICMEs, all show predominance of scintillation occurrence in the cusp. Cumulative probability distribution functions for the phase scintillation occurrence that are a by-product of superposed epoch analysis of growing databases of CIRs and ICMEs can be used in probabilistic forecasting of scintillation activity relative to the arrival times of these major solar wind disturbances. Phase scintillation occur- rence can impact on the GNSS receiver tracking per- formance, which can be assessed by the tracking jitter maps. The PLL jitter is correlated with the phase scintil- lation index and the regions of enhanced PLL jitter ap- proximately coincide with phase scintillation occurrence. Hence the probabilistic forecasting of scintillation activity should be beneficial in assessing prevailing GNSS tracking conditions at high latitudes, and particularly in the cusp. Acknowledgements. Infrastructure funding for CHAIN was provided by the Canada Foundation for Innovation and the New Brunswick Innovation Foundation. CHAIN operation is conducted in collaboration with the Canadian Space Agency. Research activi- ties at Nottingham Geospatial Institute, University of Nottingham, are funded by the UK Engineering and Physical Sciences Research Council project Polaris (http://www.bath.ac.uk/elec-eng/polaris/) and the FP7 funded project TRANSMIT (http://www.nottingham. ac.uk/transmit/index.aspx). The solar wind data were obtained from Goddard Space Flight Center Space Physics Data Facility OMNIWeb (http://omniweb.gsfc.nasa.gov/). References Aarons, J. (1982). Global morphology of ionospheric scintillations, Proc. IEE, 70, 360-378. Aquino, M., M. Andreotti, A. Dodson and H. Strange- ways (2007). On the use of ionospheric scintillation indices as input to receiver tracking models, Adv. Space Res., 40, 426-435. Basu, S., E.J. Weber, T.W. Bullett, M.J. Keskinen, E. MacKenzie, P. Doherty, R. Sheehan, H. Kuenzler, P. Ning and J. Bongiolatti (1998). Characteristics of plasma structuring in the cusp/cleft region at Svalbard, Radio Sci., 33, 1885-1899; doi:10.1029/98RS01597. Basu, S., K.M. Groves, Su. Basu and P.J. Sultan (2002). Specification and forecasting of scintillations in com- munication/navigation links: current status and fu- ture plans, J. Atmos. Sol.-Terr. Phys., 64, 1745-1754. Borovsky, J.E., and M.H. Denton (2006). Differences be- tween CME-driven storms and CIR-driven storms, J. Geophys. Res., 111, A07S08; doi:10.1029/2005JA0 11447. Conker, R.S., M.B. El Arini, C.J. Hegarty and T. Hsiao (2003). Modeling the effects of ionospheric scintil- lation on GPS/SBAS availability, Radio Sci., 38, 1001; doi:10.1029/2000RS002604. Feldstein, Y.I., and G.V. Starkov (1967). Dynamics of au- roral belt and polar geomagnetic disturbances, Planet. Space Sci., 15, 209-230. Holzworth, R.H., and C.-I. Meng (1975). Mathematical representation of the auroral oval, Geophys. Res. Lett., 2, 377-380. Horvath, I., and S. Crozier (2007). Software developed for obtaining GPS derived total electron content val- ues, Radio Sci., 42, RS2002; doi:10.1029/2006RS00 3452. Jayachandran, P.T., et al. (2009). Canadian High Arctic Ionospheric Network (CHAIN), Radio Sci., 44 (1), RS0A03; doi:10.1029/2008RS004046. McPherron, R.L., and G. Siscoe (2004). Probabilistic forecasting of geomagnetic indices using solar wind air mass analysis, Space Weather, 2, S01001; doi:10. 1029/2003SW000003. McPherron, R.L., and J. Weygand (2006). The solar wind and geomagnetic activity as a function of time rela- tive to corotating interaction regions, In: B. Tsuru- tani, R.L. McPherron, W. Gonzalez, G. Lu, J.H.A. Sobral and N. Gopalswamy (eds.), Recurrent Mag- netic Storms: Corotating Solar Wind Streams, AGU Monograph, 167, 125-137. Osherovich, V.A., J. Fainbergand and R.G. Stone (1999). Solar wind quasi-invariant as a new index of solar activity, Geophys. Res. Lett., 26, 2597-2600. Prikryl, P., P.T. Jayachandran, S.C. Mushini, D. Pokhot- elov, J.W. MacDougall, E. Donovan, E. Spanswick Figure 9. Phase scintillation index vU (a), receiver PLL jitter (b), and correlation between vU and the PLL jitter (c). 9 PROBABILISTIC FORECASTING OF SCINTILLATION and J.-P. St.-Maurice (2010). GPS TEC, scintillation and cycle slips observed at high latitudes during solar minimum, Annales Geophysicae, 28, 1307-1316. Prikryl, P., P.T. Jayachandran, S.C. Mushini and R. Chad- wick (2011a). Climatology of GPS phase scintillation and HF radar backscatter for the high-latitude iono- sphere under solar minimum conditions, Annales Geophysicae, 29, 377-392; doi:10.5194/angeo-29-377. Prikryl, P., et al. (2011b). Interhemispheric comparison of GPS phase scintillation at high latitudes during the magnetic-cloud-induced geomagnetic storm of 5-7 April 2010, Annales Geophysicae, 29, 2287-2304; doi:10.5194/angeo-29-2287. Prikryl, P., P.T. Jayachandran, S.C. Mushini and I.G. Richardson (2012). Toward the probabilistic forecast- ing of high-latitude GPS phase scintillation, Space Weather, 10, S08005; doi:10.1029/2012SW000800. Prikryl, P., Y. Zhang, Y. Ebihara, R. Ghoddousi-Fard, P.T. Jayachandran, J. Kinrade, C.N. Mitchell, A.T. Weath- erwax, G. Bust, P. J. Cilliers, L. Spogli, L. Alfonsi, V. Romano, B. Ning, G. Li, M.J. Jarvis, D.W. Danskin, E. Spanswick, E. Donovan and M. Terkildsen (2013a). An interhemispheric comparison of GPS phase scin- tillation with auroral emission observed at the South Pole and from the DMSP satellite, Annals of Geo- physics, 56 (2), R0216; doi:10.4401/ag-6227. Prikryl, P., R. Ghoddousi-Fard, B.S.R. Kunduri, E.G. Thomas, A.J. Coster, P.T. Jayachandran, E. Span- swick and D.W. Danskin (2013b). GPS phase scintil- lation and proxy indices observed at high latitudes during a moderate geomagnetic storm, Annales Geophysicae, 31, 805-816; doi:10.5194/angeo-31- 805-2013. Richardson, I.G., D.F. Webb, J. Zhang, D.B. Berdichevsky, D.A. Biesecker, J.C. Kasper, R. Kataoka, J.T. Stein- berg, B.J. Thompson, C. Wu and A.N. Zhukov (2006). Major geomagnetic storms (Dst ≤ −100 nT) gener- ated by co-rotating interaction regions, J. Geophys. Res., 111, A07S09; doi:10.1029/2005JA011476. Richardson, I.G., and H.V. Cane (2010). Near-Earth in- terplanetary coronal mass ejections during solar cycle 23 (1996-2009): catalog and summary of prop- erties, Solar Phys., 264, 189-237; doi:10.1007/s11207- 010-9568-6. Skone, S., and M. de Jong (2000). The impact of geo- magnetic substorms on GPS receiver performance, Earth Planets Space, 52, 1067-1071. Spogli, L., Lu. Alfonsi, G. De Franceschi, V. Romano, M.H.O. Aquino and A. Dodson (2009). Climatology of GPS ionospheric scintillations over high and mid- latitude European regions, Annales Geophysicae, 27, 3429-3437. Sreeja, V., M. Aquino and Z.G. Elmas (2011). Impact of ionospheric scintillation on GNSS receiver tracking performance over Latin America: introducing the concept of tracking jitter variance maps, Space Weather, 9, S10002; doi:10.1029/2011SW000707. Van Dierendonck, A.J., and B. Arbesser-Rastburg (2004). Measuring ionospheric scintillation in the equato- rial region over Africa, including measurements from SBAS geostationary satellite signals, Proc. ION GNSS 17th technical meeting of the satellite divi- sion, Long Beach, CA, USA, 316. *Corresponding author: Paul Prikryl, Communications Research Centre, Ottawa, Ontario, Canada; email: paul.prikryl@crc.gc.ca. © 2013 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice