Next Article in Journal
Assessment of New Satellite Missions within the Framework of Numerical Weather Prediction
Next Article in Special Issue
The First Pi2 Pulsation Observed by China Seismo-Electromagnetic Satellite
Previous Article in Journal
An Evaluation of Citizen Science Smartphone Apps for Inland Water Quality Assessment
Previous Article in Special Issue
Nadir-Dependent GNSS Code Biases and Their Effect on 2D and 3D Ionosphere Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Small-Scale Ionospheric Irregularities of Auroral Origin at Mid-latitudes during the 22 June 2015 Magnetic Storm and Their Effect on GPS Positioning

1
Institute of Solar-Terrestrial Physics of Siberian Branch of Russian Academy of Sciences, 664033 Irkutsk, Russia
2
Lomonosov Moscow State University, Faculty of Physics, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(10), 1579; https://doi.org/10.3390/rs12101579
Submission received: 16 April 2020 / Revised: 12 May 2020 / Accepted: 14 May 2020 / Published: 15 May 2020
(This article belongs to the Special Issue Remote Sensing of Ionosphere Observation and Investigation)

Abstract

:
Small-scale ionospheric irregularities affect navigation and radio telecommunications. We studied small-scale irregularities observed during the 22 June 2015 geomagnetic storm and used experimental facilities at the Institute of Solar-Terrestrial Physics of the Siberian Branch of the Russian Academy of Sciences (ISTP SB RAS) located near Irkutsk, Russia (~52°N, 104°E). The facilities used were the DPS-4 ionosonde (spread-F width), receivers of the Irkutsk Incoherent Scatter Radar (Cygnus A signal amplitude scintillations), and GPS/GLONASS receivers (amplitude and phase scintillations), while 150 MHz Cygnus A signal recording provides a unique data set on ionosphere small-scale structure. We observed increased spread-F, Cygnus A signal amplitude scintillations, and GPS phase scintillations near 20 UT on 22 June 2015 at mid-latitudes. GPS/GLONASS amplitude scintillations were at a quiet time level. By using global total electron content (TEC) maps, we conclude that small-scale irregularities are most likely caused by the auroral oval expansion. In the small-scale irregularity region, we recorded an increase in the precise point positioning (PPP) error. Even at mid-latitudes, the mean PPP error is at least five times that of the quiet level and reaches 0.5 m.

Graphical Abstract

1. Introduction

Small-scale ionospheric irregularities affect radio telecommunications, because those irregularities can scatter electromagnetic signals [1]. The most pronounced effect is observed when the irregularity scale is of the radiowave’s first Fresnel zone ~ √λz (where λ is the signal wavelength, z is the distance to the irregularity). The scattering leads to radio signal amplitude scintillations that deteriorate the quality of global navigation satellite system (GNSS) operation. In the extreme case, irregularities may even cause losses-of-lock [2] and interrupt positioning services.
Intensive amplitude scintillations and GNSS signal losses-of-lock are usually observed in the high-latitude and equatorial ionosphere [3,4,5,6,7,8], where irregularities are effectively generated. Equatorial small-scale irregularities are usually seen as spread-F on ionograms. Hence, a long-term history of irregularity observations is provided by ionosondes [9,10,11]. At the same time, there are studies showing that such irregularities could exist in the mid-latitude ionosphere [12,13,14]. The auroral oval broadening is a reason for an irregularity to be generated at the middle latitudes [15,16]. Another reason is a large-scale ionospheric structure (including traveling ionospheric disturbances (TIDs) [17]) that appeared outside the mid-latitude region and traveled down to mid- latitudes. For example, an equatorial super plasma bubble can penetrate into the mid-latitude ionosphere [18] and produce small-scale irregularities that manifest themselves as the GPS losses of phase lock [19]. However, there is no large amount of amplitude scintillation data at mid-latitude. This is due to a small number of facilities at the mid-latitudes and insufficient coverage of the radio observations.
Several approaches have appeared recently to overcome this problem. One of them concerns indirect measurements. The Low Frequency Array (LOFAR) project [20] reveals a promising capability. Another approach is measuring scintillations of discrete space radio sources by incoherent scatter radars. The latter approach has been developed at the Institute of Solar-Terrestrial Physics of the Siberian Branch of the Russian Academy of Sciences (ISTP SB RAS) [21] and enables researchers to obtain unique data.
The small-scale irregularities (the first Fresnel zone scale) almost do not affect the carrier phase. It is due to power-law dependence of ionospheric irregularities amplitudes. Yeh and Liu [1] showed that only irregularities with the size 10–102 times greater than those of the first Fresnel zone have a noticeable effect on the phase. It is conventional to study small-scale irregularities with instruments operating at wavelength within the meter-decameter range. However, such an approach provides the information only on a small part of the irregularity spectrum. To get a comprehensive study, one should employ instruments capable of sensing various part of the irregularity spectrum simultaneously.
The 22-23 June 2015 geomagnetic storm was “an anomalously intense event with large magnetic fields” [22]. It was one of the most intense storms of solar cycle 24 [23]. A number of significant papers have described and analyzed this event. Astafyeva et al. [23] observed an extreme enhancement of the electron density Ne and vertical total electron content VTEC in the Northern (summer) Hemisphere. They concluded that the extreme topside response was the result of a combination of the prompt penetration electric fields, disturbance dynamo, and the storm-time thermospheric circulation. Astafyeva et al. [24] studied global the ionospheric and thermospheric effects of the storm. They found that dayside neutral mass density increased by three to five times as compared to the quiet time level. Astafyeva et al. [25] studied the equatorial electrodynamics. They found an increase in the eastward electric field and in the upward E×B drift on the dayside (which led to rapid strong enhancement of the VTEC and electron/ion density), as well as downward vertical drift with a decrease in the VTEC and in the plasma density on the nightside. Singh and Sripathi [26] studied small-scale structures in the equatorial region. They suggested that the westward penetration of the electric field into the equatorial sector during the local midnight caused an abrupt decrease in the virtual height to ~200 km and suppressed plasma bubbles in the Indian sector. Comprehensive analysis was carried out by Piersanti et al. [27] who showed an increase in the ionospheric polar convection and an expansion of the lower boundary of the convection region toward low latitudes (from 60°N to 50°N).
Geomagnetic storms and TIDs are also known as sources for GPS positioning quality deterioration. GPS scintillations lead to a range error in GNSS due to diffraction [28]. Deep signal fades that appear during small-scale irregularities effect result in navigation outages [29]. Luo et al. [30] studied several strong storms during solar cycle 24, including the 17 March 2015 St. Patrick storm. They reported that the precise point positioning (PPP) error was less than 0.32 m at mid-latitudes.
In this paper, we study small-scale structures at mid-latitudes during the 22 June 2015 magnetic storm. The mid-latitude irregularity formation mechanism is not quite clear. We have tried to fill the knowledge gap by recording new scintillation data in different frequency bands at mid-latitudes. The data help to reveal physical mechanism for the irregularity generation in the mid-latitude ionosphere. Another point is to reveal how the irregularities affect GPS precision during the storm.

2. Experimental Facilities and Data

We used the ISTP SB RAS facilities to comprehensively study small-scale irregularities. We employed the S4 index [31] to study amplitude scintillations and σφ [32] to study phase scintillations. We used the spread-F [33] width to study irregularities within the ionosonde wavelength range. The S4 index is the signal intensity E standard deviation normalized to the average signal intensity:
S 4 = E 2 E 2 E .
The σφ index is simply the carrier phase φ standard deviation:
σ ϕ = ϕ 2 ϕ 2 .
The spread-F width is the measure of trace broadening on an ionogram near the ionosphere peak.
Below, we describe the instruments and data:
(1) Ionosonde DPS-4 [34] enabled us to measure the spread-F width. We used the ionosonde located in Irkutsk – IR352 (52.3°N, 104.3°E; red star in Figure 1), at ~48° magnetic latitude. We manually scaled the ionograms to obtain the spread-F width. The ionosonde antenna is not typical. This did not enable us to measure the directions of signal receiving. During the geomagnetic storm, the ionosonde in Irkutsk was operating with a 15-min cadence. Additionally, we analyzed data from the ionosondes in Norilsk—NO369 (69.4°N, 88.1°E), Yakutsk—YA462 (62.0°N, 129.7°E), Mohe—MH453 (52.0°N, 122.5°E), and Moscow—MO155 (55.47°N, 37.3°E). As the ionosonde records the spread-F within ~3–10 Mhz at 200–300 km, the irregularities that primarily affect the ionosonde signal have sizes 2.5–5 km.
(2) Irkutsk Incoherent Scatter Radar (IISR) [35]. The IISR is located next to Irkutsk (52.9°N, 103.3°E; red star in Figure 1), at ~49° magnetic latitude. Its radiation pattern center is at 52.6°N, 102.9°E. The IISR receivers enabled us to measure scintillations of discrete space radio sources, when they crossed the radar sector of view [21]. In our experiments, we recorded the Cygnus A signal. Cygnus A is the most powerful radio galaxy in our part of the Universe. Cygnus A emits in a broad frequency band including the IISR frequency band (154–162 Mhz). We should note that the radiation pattern is shown for the chosen frequency. By changing the frequency, we changed the radiation pattern, and thus could track the radio source [21].
Unfortunately, during the storm, the IISR was operating in passive mode. So, we have no corresponding electron density profiles or ion velocities. The highest temporal resolution for measurements in the passive mode is of the order of seconds. We used the 4.5-s resolution data. The integration time for the S4 index calculation is 1 minute. Since the IISR operates at ~150 MHz, at 100–300 km, the size of the irregularities primarily affecting the IISR signal amplitude is 400–800 m.
(3) GPStation-6 specialized GPS/GLONASS receiver. The receiver is located in Irkutsk (52.3°N, 104.3°E; red star in Figure 1), at ~48° magnetic latitude. GPStation-6 receiver enabled us to measure the amplitude scintillations intensity, S4, and phase variations, σφ [36]. The data sampling rate used for the S4 and σφ calculation is 50 Hz (0.02-s resolution). The integration time for the S4 and σφ is 1 minute. The S4 and σφ measured by GPStation-6 were transformed into equivalent vertical values [37]. The GPStation-6 receives signals within the L1 and L2 bands for all the GPS/GLONASS satellites. There is an option with some satellites to receive a third frequency band, L5 and L3 for GPS and GLONASS, respectively. As receivers operate at ~1.6/1.2 GHz (L1/L2) at 100–300 km, the size of the irregularities primarily affecting the navigation signal amplitude is 100–300 m.
(4) Additionally, we used the total electron content (TEC) variations obtained by using the dual frequency GPS/GLONASS phase measurement. The technique for the TEC calculation is well-known [38,39,40].
I = 1 40.308 f 1 2 f 2 2 f 1 2 f 2 2 [ L 1 λ 1 L 2 λ 2 + c o n s t + σ ϕ ] ,
where L1 λ 1 and L2 λ 2 are the radio signal carrier phase ranges including ionospheric delay (m), L1 and L2 are a number of phase rotations at f1 and f2 frequencies; λ1 and λ 2 are corresponding wavelengths; const is an unknown initial phase ambiguity (m); and σϕ is a term corresponding to noise in phase measurements (m). By convention, 1 TEC Unit (TECU) = 1016 electrons/m2.
We used the data from 5172 receivers of the global and regional GPS/GLONASS receiver networks (see Acknowledgments), including the IGS [41], UNAVCO, SONEL, etc., as well as the regional data, like SibNet [42], CHAIN [43], EUREF [44] HIVE (https://hive.geosystems.aero/), etc. Blue dots in Figure 1 show the location of the used receivers.
Single cycle slips in the TEC data were eliminated by the formula:
TEC(i) = TEC(−1) + TEC(i) − TEC(0) + [TEC(−1) − TEC(−2)]/2 + [TEC(1) − TEC(0)]/2.
where i is the time instance after a cycle slip, TEC(0) is the TEC value, when the cycle slip is recorded, TEC(1) corresponds to TEC records after cycle slip, TEC(−2) and TEC(−1) correspond to two TEC records before the cycle slip occurred. We used an approximation that the time derivative equals the mean derivative before the cycle slip and after it. A similar approach was used to eliminate loss-of-lock effects.
Variations in TEC within 15 min–3 h correspond to TIDs and other wave-like ionospheric disturbances. As TEC variations, we considered detrended and filtered TEC values. We performed the following procedure [45]. We detrended TEC data based on smoothing spline. After this, we used the running average with a 10–20 min (20–60 min) window to filter variations corresponding to medium-scale travelling ionospheric disturbances. Thus, we obtained the TEC variations series dI. We used a single layer model [46] as a mapping function to reduce TEC variations’ dependence on the elevation angle. We mapped the variations on the globe given the ionospheric pierce point location and obtain the global spatial distribution of TEC variations [40]. We marked each ionospheric pierce point (at 300 km height) with the color corresponding to the TEC variation value. All calculations of TEC variations were performed by SIMuRG (https://simurg.iszf.irk.ru) [47].
(5) To estimate the GPS precision, we calculated the station coordinates in the kinematic GPS dual-frequency precise point positioning (PPP) mode. For that purpose, we used the GAMP open-source software [48]. Receiver and satellite clock offsets were considered in GAMP PPP solution by applying IGS precise satellite orbit and clock products. In addition, using linearized equations of original pseudorange and carrier phase observations, GAMP considers line-of-sight (LOS) ionospheric delay, receiver non-calibrated code delay, receiver and satellite non-calibrated phase delays, and zenith wet delay. We used a dual-frequency PPP model, in which the receiver non-calibrated code delays are absorbed by both receiver clock offset and LOS ionospheric delay parameters. The 24-hour median values of X, Y, Z coordinates for a station were regarded as reference positions. The positioning error was calculated as the difference between the reference and the instant position
σ = Δ X 2 + Δ Y 2 + Δ Z 2 .
By using the worldwide GNSS network, we can obtain the accuracy of the positioning distribution similar to that in [49].

3. Results

3.1. 22 June 2015 Geomagnetic Storm

The storm is related to active region (AR) 12,371 on the Sun according to the SOHO LASCO CME catalog (https://cdaw.gsfc.nasa.gov/CME_list/; SOHO—the Solar and Heliospheric Observatory, LASCO—the Large Angle and Spectrometric Coronagraph, CME—Coronal Mass Ejection). The AR produced several CMEs. The 22 June 2015 magnetic storm is mainly due to ~1350 km/s Halo CME that appeared at 2:36. The CME with kinetic energy of ~9.0∙1031 erg (http://cdaw.gsfc.nasa.gov/CME_list/) is associated with an M2.7 flare with maximum at 02:34 UT. On 22 June, a cluster of shocks passed the Advanced Composition Explorer (ACE) spacecraft at 17:58 UT [22]. Figure 2 presents time series of the data to help tracking storm evolution. We see a sharp increase in the proton density (Figure 2c) and a sudden storm commencement (SSC) (Figure 2a) at ~05:45 and ~18:30 UT on 22 June, when the shocks hit the Earth’s magnetosphere.
The data presented in Figure 2 (Bz component of the interplanetary magnetic field (IMF), H-SYM index, AE-index, and proton density) were retrieved from the OMNI database (ftp://spdf.gsfc.nasa.gov) [50]. The spread-F (Figure 2c, red line) was measured based on the Irkutsk ionosonde data. As the shock hits the 22 June 2015 magnetosphere at 18:30 UT, we see the storm onset as a SSC, and the H-SYM index starts decreasing. The proton density changed from 12 cm−3 (before the shock) to 52 cm−3 (after the shock). The H-SYM had a local minimum (135 nT) at 20:15 UT. After 20:15 UT, the magnetosphere was recovering. The recovery of the magnetosphere lasted until 00:50 UT, when Bz turned south again, and the proton density decreased to ~9 cm−3. The minimal H-SYM was observed at 04:25 UT with the value of −207 nT. We should note that the substorm occurred before the storm. The AE-index shows an increase after 14 UT on 22 June. H-SYM dropped to −50 nT at ~16:15 UT.

3.2. Spread-F and Scintillations

Figure 3 shows the experiment geometry. The left panel shows the GPS ionospheric pierce points (IPP) at 300 km as black lines with 10–20 minute TEC variations superimposed, red curves, at 19:30–20:30 on 22 June 2015. The vertical segment in the right bottom part of panel (a) shows 0.5 TECU variations. The data shown are measured by the ISTP SB RAS station within the SibNet. The TEC variations are shown only for four satellites (PRN15, PRN20, PRN21, PRN22). The right panel shows the geometry in details. Black lines show the GPS IPP (the same as on the left panel). Gray lines show the IPP of Cygnus A at 100, 250, and 400 km for the IISR location that is 60 km away from the ISTP SB RAS station. Filled gray region and filled blue region are approximate radiation patterns (at 250 km) of the ionosonde and of the IISR, respectively. The ionosonde was not capable of measuring the direction of the received signals during the interval of interest. Thus, we only can attribute the spread-F data to somewhere in the gray region. Figure 3 allows us to see that the measurements are almost in the same region.
We observed the spread-F on 22 June 2015 (Figure 2b) rather than on 23 June 2015. Note that, when the 22 June 2015 H-SYM is minimal at 04:25 UT, there was no noticeable spread-F (Figure 2b). There were occasional black-outs that prevented reliable processing of the ionograms. However, when there is no black-out, there is still no spread-F on the ionograms. The solar wind proton density at this time was even smaller than that of the quiet time. This leads to a smaller energy input into the auroral region; hence, no large-scale disturbances exist and no small-scale irregularities are generated.
Figure 4 shows the spread-F width time series from the ionosondes in Irkutsk (a) as well as in Yakutsk, Mohe, and Moscow (d). Panel b shows phase scintillations σφ calculated by using observations from all the visible GPS satellites (all the GPS satellite data are mixed here), and Cygnus A radio signal amplitude scintillations (c). The 22 June 2015 data are plotted for 15:00 – 24:00 UT. On the panel c, we also show the 24 June 2015 data (grey) as a quiet time reference. Dashed-dotted line shows the angle γ between the “IISR – Cygnus A” line of sight and the Earth magnetic field calculated based on the IGRF−12 model [51].
We see the spread-F enhancement 19:45 UT through 21:15 UT in Irkutsk. After the magnetic storm onset, the spread-F was at the background level for an extra hour. The background was ~0.08 MHz, while the maximal spread-F width was 0.46 MHz at 20:15 UT. We also observed a σφ increase for ~19:30–19:50 UT. The σφ mean value was 0.07 rad, with ~0.1–0.2 rad during the enhancement. We did not observe any significant increase in the GPS/GLONASS navigation signal amplitude scintillation. During the time interval presented in Figure 4, the S4 index did not exceed the background level, and, hence, we do not show the GPS/GLONASS S4 here.
Unfortunately, there are no observations from the Norilsk ionosonde during the storm. Black-out caused by strong absorption in the D-region started shortly after the SSC and prevented collecting useful information from the ionograms. However, we were able to retrieve the data for the ionosondes (see Figure 4d) in Yakutsk (blue line), Mohe (black line), and Moscow (red line).
The spread-F is observed at all the stations. A strong spread-F appears at the similar time at Mohe and Moscow. Also, the time was close to that when we observed the spread-F at Irkutsk (20:15 UT, 22 June 2015). Small deviations for the Moscow station are explained by its slightly higher latitude than that of Irkutsk and Mohe. There is no noticeable spread-F appearance dependence on the longitude. At Yakutsk, the spread-F is observed early, which we attribute to a higher latitude of that station as compared to others. The above (Section 3.1) substorm results in spread-F at Yakutsk latitude before the main storm onset.
The Cygnus A was observed within the radar field of view 17:00 through 22:00 UT. The Cygnus A signal scintillation index S4 was calculated by using formula (1). Figure 4c shows the result. The black continuous line shows the 22 June 2015 data. On 22 June 2015, the S4 first increase was observed between 18:30 and 19:30 UT. A similar increase was recorded on 24 June 2015, the day that was chosen as a reference day (grey line). The increase corresponds to the angle γ minimum, i.e., it is the magnetic zenith effect. Such an effect is regularly observed in the IISR scintillation data and fades quickly, when the angle between the LOS and the magnetic field exceeds 10° [52]. The highest S4 for Cygnus A was observed around 20:00 UT. It coincides with the time, when the enhanced spread-F was observed. The scintillation intensity dropped at 20:30 UT, while the spread-F lasted for an extra hour.

3.3. Precise Point Positioning

To estimate the GPS positioning accuracy, we calculated the PPP error (5) for the worldwide network. The dynamics for the PPP error spatial distribution can be found in the Supplementary Materials (Video S1). Figure 5 shows the average PPP error vs. latitude. The latitude bin size is 2.5°. PPP errors are averaged over −0°…180° (European-Asian sector, upper panels a-b) and −180°…0° (American-Atlantic sector, bottom panels c-d). The dashed line on panels b and d marks 18:30 UT. Panels a and c show the data for the reference day (21 June 2015, 172nd Day of the Year, DOY), while panels b and d show the data for the 22 June 2015 storm (173 DOY).
The results for the reference day show that the PPP error increases due to satellite geometry or data processing peculiarities. Such effects have similar patterns for both the storm and the reference days. For example, there is a PPP error increase at ~10UT in the Northern Hemisphere. Note that, over 0–1(2) UT, the quality of the PPP solution obtained by the GAMP software is significantly worse than that at other times. So, we do not consider this time. In Figure 5, the navy blue regions (zero values) mark lack of the data. Data gaps can be seen on both the reference and the storm days. One can see such patterns at the North Pole region, as well as at ~60°S, where there are no stations for calculations.
The typical average errors are less than 0.5 m (and even less than 0.1–0.3 m). From Figure 5, in the Northern Hemisphere after 18:30 UT (dashed line), one can see a sharp decrease in the PPP quality at 65–75°N. The error exceeded 0.5 m. So, we recorded a precision deterioration, which results in fivefold worse positioning accuracy. There is an equatorward “propagation” of the error with time. In the American-Atlantic sector (Figure 5d), we can observe such a “propagation” as far as ~50°N. The PPP error maximum appears at ~20 UT at this latitude. In the European-Asian sector, we observed a similar pattern with an increased error at ~20 UT, but not in such details due to fewer stations. There is a signature that propagated to the equatorial region in the Northern Hemisphere, but the effect is not so pronounced, and it is observed during the reference day.

4. Discussion

Using the data from Figure 4, we can conclude that, during the geomagnetic storm, 500–5500 m small-scale irregularities were present in the ionosphere. The absence of the GPS/GLONASS amplitude signal scintillations means that there were no intensive irregularities with corresponding scales. The ionospheric irregularity amplitude exponentially decreased as the irregularity size decreased [1,53]. That is why there were no significant 100–300 m irregularities.
We also examined mechanisms for the existence of such intensive small-scale irregularities. We had two hypotheses. The first possible mechanism is an expansion of the auroral oval toward mid-latitudes, so we would observe auroral phenomena at mid-latitude stations. Particle precipitation results in local increases and decreases in the plasma density more than 10 km wide. Such structures can be cascaded into smaller scales through instability processes and turbulent diffusion [54]. The 0.1–30 km irregularities are generated at high latitudes mainly due to ExB instability [55]. Also, the thermomagnetic instability can produce F-region 1 km-scale irregularities at high-latitudes [56].
The second mechanism suggests generation of small-scale irregularities due to passage of a medium-scale (MS) TIDs [57]. For example, generating medium-scale irregularities due to the large-scale disturbance passage was mentioned by Astafyeva et al. [16]. Astafyeva et al. [16] found medium-scale 2–10 min irregularities, when a well-distinguished strong disturbance of the 40-min time scale passed by. There are a number of papers by G.G. Bowman focusing on the spread-F caused by TIDs. Bowman and Mortimer [58] showed that, during the passage of a large scale TID, the average ionospheric F2 layer height rises. Bowman [17] showed that the spread-F observed in equatorial regions several hours after increased geomagnetic activity is directly related to ionosphere height rise. Studies by Farley et al. [59] indicated that a favorable condition for the instability growth is induced rise of the ionosphere height. The “rise-in-the-ionospheric-layer-height” mechanism, which makes the spread-F observable at equatorial stations, may also produce the spread-F at mid-latitudes [58]. In fact, gravity waves are the prevailing source for both MS TIDs and small-scale irregularities [60]. Gravity waves induce currents that drive plasma instability.
To find a proper mechanism, we used, first, the OVATION model data and the SupernDARN data to estimate the oval boundary dynamics and its correspondence with the observed small-scale irregularity appearance. Second, we studied the global distribution of TEC variations. On the one hand, the global TEC variations maps can reveal TID origin and propagation. On the other hand, the maps provide some evidence of small-scale structure development. Third, we analyzed the changes in PPP error over the globe. An increase in the PPP error and its localization can provide additional information about the source of the recorded small-scale irregularities.
To study oval expansion we used the OVATION prime energy flux data (not shown here) (http://iswa.gsfc.nasa.gov; https://sourceforge.net/projects/ovation-prime/) [61]. For the time interval of interest, the oval boundary did not reach 52°N for the Asian sector. However, the OVATION is a model and cannot represent the actual precipitation (especially during such intensive storms). Based on SuperDARN data, Piersanti et al. [27] studied the polar convection for this very storm and concluded that the auroral oval expanded to ~50°N (we should note that it is magnetic latitude). We also used the SuperDARN data to study this expansion in detail. The data are available at http://vt.superdarn.org/ [62]. Figure 6 shows the maps for the electric potential from SuperDARN for 18:30 UT (a), 18:40 UT (b), 19:00 UT (c), 19:30 UT (d), 20:00 UT (e), 20:30 UT (f), 21:00 UT (g), 21:30 UT (h). Green dashed curve shows the Heppner–Maynard boundary.
From Figure 6 one can see that, after 18:30 UT, the expansion of auroral oval started. In the Irkutsk sector (look for Lake Baikal at the map bottom), the Heppner–Maynard boundary reached 50° magnetic latitude at ~19:00 UT and persisted at this position until ~21:30 UT. Over that time, the boundary approached and then moved away from the facilities without reaching them. However, the boundary was very close to the radiation patterns of the instruments that we used (Figure 3). The ionosonde has wider radiation patterns. So, it was able to record more reflections from small-scale structures due to the expanded oval than the Irkutsk incoherent scatter radar was. The difference in the observation time seems to be related to the radiation pattern difference. This point proves the hypothesis that the oval expansion is the reason for generating small-scale irregularities.
Another hypothesis is that small-scale irregularities are caused by the passage of a large-scale disturbance. We assume that a small-scale irregularity source is originated in the auroral area and then travels down to mid-latitudes. To prove this, we used TEC variation maps (see Experimental Facilities). On these maps, the TEC variation amplitude is plotted. The point positions correspond to the ionospheric pierce points (300 km height). Figure 7 shows the 10–20 min TEC variations maps for 18:30 UT (a), 19:00 UT (b), 19:30 UT (c), 20:00 UT (d), 20:30 UT (e), 21:00 UT (f). Figure 8 provides similar maps, but for 20–60 min TEC variations. For reference, we provide the movies that illustrate the global spatial-temporal distribution of the 10–20 min and 20–60 min TEC variations for 17:00–24:00 UT at 1-min time resolution (see Supplementary Materials, Videos S2 and S3).
There are always sporadic disturbances in the auroral oval that do not form any pronounced structures. We can clearly see this before the storm from Videos S2 and S3. After the storm onset at 18:30 UT, we see a large-scale well-defined negative disturbance that covers all longitudes. Later, we see a positive disturbance. We can see a long negative front at 19:30 UT reaching the northern part of the Irkutsk region. Maps of 20–60 min TEC variations show an extended sharp “edge” in the European-Asian region. The “edge” dynamics corresponds to that for the Heppner–Maynard boundary from the SuperDARN data.
The 10–20 min TEC variation dynamics (Video S2) show that, in the European-Asian region, the irregularities of the auroral origin (“irregularity cloud“) approached 50°N. The same spatial limit is seen from maps of 20–60 min TEC variations. This agrees with the results in Figure 2 (left panel), where TEC variations exceeding 0.1 TECU were recorded only for the IPPs located northward 52°N. These auroral-origin irregularities, presumably, are manifestation of delayed instability development in the auroral oval. The recorded TEC variations (Figure 2b, Figure 7, Figure 8) can be related to the movement of the boundary with an increased ionization level. We should also note that the polar cap features smaller TEC variations. Along with the propagation of the auroral–origin boundary, we can see the polar cap expansion.
After ~19:20 UT, we can observe an appearance of large-scale TID, close to the position of the “edge” related to the oval expansion. Thus, there were observed both auroral oval expansion and TID generation on the auroral oval boundary. We should note that it is difficult to definitely separate these two phenomena due to the complexity of the considered event.
We recorded a TID in the Northern Hemisphere with an extended wave front. In Figure 7f, we see large positive TEC variations crossing the equator in the American region. According to the movie (see Supplementary Materials, Videos S2 and S3), the disturbance originated in the Northern Hemisphere, crossed the equator, and reached ~40°S. The disturbance, hence, travels more than 10,000 km from its source. The disturbance can be seen in the American and Australian longitude sectors, which implies that it is ~200° wide in longitude. It can be clearly identified in the American-Atlantic region, especially, when it was reaching the equatorial region (see Figure 7e, Figure 7f). TID is also seen over the European region, and a part of its structure can be found over Northern Africa. It is difficult to clearly identify TID southward over the Northern African and Asian regions due to low data coverage.
Finally, we can compare the TEC variations maps with the maps for the PPP quality. The observed PPP error does not coincide with the sharp oval extension (Figure 6). The oval reached 50° magnetic latitude in the Asian sector quickly, while the PPP deterioration at such latitudes was recorded only at ~20 UT. However, TEC variation maps (Figure 7, Figure 8 and Movies S1, S2) also show a gradual propagation of irregularities from the oval. Moreover, in the American sector, we do not observe any “propagation” of the PPP error to the Southern Hemisphere, as we observe for TIDs. If the TID was the source for scintillations resulting in enhancement in PPP errors at mid-latitudes, we would expect to see an increase in PPP error along the TID propagation trajectory. However, we did not observe any “propagating” pattern (up to the Southern Hemisphere) in PPP error distributions. While this is not an incontrovertible conformation, it indirectly indicates that TID propagation is not a source for the increase in PPP error (and scintillations). In contrast, this is additional evidence in favor of the auroral oval expansion mechanism for the revealed small-scale irregularities.
Thereby, the auroral oval expansion and the development of small-scale structure seem to be a more preferable mechanism to produce the observed small-scale irregularities at mid-latitudes. The irregularities lead to increased phase scintillations of navigation signal. The increase is more pronounced at the start of the irregularity formation. Phase scintillations decay very quickly. We observed spread-F, Cygnus A signal amplitude scintillations, and GPS phase scintillations simultaneously. The spread-F width maximum was recorded after the GPS phase scintillations recover their background level. The ~2.5–5.5 km irregularities (ionosonde) seem to persist over the entire interval, when the spread-F is observed. The ~400–800 m irregularities (IISR) emerged after the GPS phase scintillations started and ended before the enhanced spread-F ended. We refer the absence of the GPS signal amplitude scintillations to the exponential shape of the irregularity spectrum. Consequently, the irregularities with the size proper for a significant effect on the navigation signal had too low amplitudes to affect the latter.
The results on the PPP error agree, in general, with those obtained by Luo et al. [30]. However, while Luo et al. [30] reported a <0.32 m PPP error, we recorded a higher error (>0.5 m) even at mid-latitudes. We recorded, at least, the GPS positioning error being five times increased at mid-latitudes during several hours of the storm.

5. Conclusions

We observed small-scale irregularities as the 22 June 2015 geomagnetic storm evolved. We used spread-F data from ionosonde, and the Cygnus A signal amplitude scintillation data from IISR. The most pronounced effect is observed when the irregularity scale is of the radiowave first Fresnel zone ~ √λz (where λ is the signal wavelength, z is the distance to the irregularity). An approximate distance z is ~100(ISSR)/200(ionosonde)−300 km. As the ionosonde recorded F-spread within the 3–10 MHz band and the Cygnus A signal frequency is 150 MHz, the irregularity size is within ~0.5 km (IISR) and ~5.5 km (ionosonde) range.
The exponential behavior of the irregularity spectrum explains why we did not observe noticeable GPS/GLONASS signal amplitude scintillations. It also explains why the ionosonde signal scintillations were more pronounced than those of Cygnus A (received by IISR). In addition, a longer observation time of small-scale irregularities by the ionosonde may be connected to the radiation pattern difference (ionosonde has a wider radiation pattern than the IISR).
The spatial distribution of the TEC variations and the maps for the electric potential from the SuperDARN data revealed that small-scale irregularities are most likely caused by an the auroral oval expansion. Scintillations and the PPP errors at mid-latitudes increased, when the irregularity region with an increased level of the TEC variations approached the measuring facilities. At mid-latitudes, the PPP errors exceeded the background level in the region of small-scale irregularities appearance over several hours by at least fivefold.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/12/10/1579/s1, https://www.mdpi.com/2072-4292/12/10/1579/s2, https://www.mdpi.com/2072-4292/12/10/1579/s3; Video S1: Global dynamics for the GPS PPP error as derived from a ground-based GPS-receiver for 0 through 24 UT on 22 June 2015. The animation cadence is 1 min. The time (in UT) is shown on the top; Video S2: Global dynamics for the 10–20 min TEC variations as derived from a ground-based GNSS-receiver for 17 through 24 UT on 22 June 2015. The animation cadence is 1 min. The scale is on the right. The corresponding UT is shown in the animation title; Video S3: Global dynamics for the 20–60 min TEC variations as derived from a ground-based GNSS-receiver for 17 through 24 UT on 22 June 2015. The animation cadence is 1 min. The scale is on the right. The corresponding UT is shown in the animation title.

Author Contributions

Conceptualization and Methodology, Y.Y., R.V. and K.R. Data curation and formal analysis: IISR—R.V., A.S. and M.G.; Ionosonde—K.R. and A.V.; GNSS—Y.Y., A.Y., A.K. and A.V. PPP calculation and formal analysis, S.S. Investigation, Y.Y, R.V., K.R. and A.V. Writing—Original Draft Preparation, Y.Y., A.V., R.V. and K.R. Writing—Review & Editing, Y.Y., A.V. and A.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by Russian Federation President Grant No. MK-3265.2019.5 and by Russian Foundation for Basic Research Grant No. 19-05-00889. The study is based on the data recorded by using the ISTP SB RAS Angara Multi-access Center facilities (http://ckp-rf.ru/ckp/3056/), including the Unique Research Facility Irkutsk Incoherent Scatter Radar (http://ckp-rf.ru/usu/77733/), under budgetary funding from Basic Research Program II.16.

Acknowledgments

The authors thank Alexandra Kustavinova for the graphical abstract made for the paper. We acknowledge the Scripps Orbit and Permanent Array Center, IGS [41], Bundesamt für Kartographie und Geodäsie (BKG) Data Center, WuHan University, Crustal Dynamics Data Information System (CDDIS), Korea Astronomy and Space Science Institute, National Geodetic Survey (Data from NOAA's National Geodetic Survey (NGS) Continuously Operating Reference Station (CORS) network of Global Navigation Satellite System (GNSS) (https://www.ngs.noaa.gov/CORS/data.shtml), the UNAVCO Facility (supported by the National Science Foundation (NSF) and National Aeronautics and Space Administration (NASA) under NSF Cooperative Agreement No. EAR-0735156), Système d'Observation du Niveau des Eaux Littorales (SONEL), the EUREF Permanent Network Services [44], Geoscience Australia, the New Zealand GeoNet project (supported by EQC, GNS Science and LINZ), the State GPS network of the Republic of Bulgaria, Institut Geographique National, Instituto Geográfico Nacional, Instituto Tecnológico Agrario de Castilla y León, Geodetic Data Archiving Facility, Instituto Brasileiro de Geografia e Estatística, Canadian High Arctic Ionospheric Network [43] for the GPS/GLONASS data. We are grateful to the Industrial Geodetic Systems (HIVE) (https://hive.geosystems.aero/) and, especially, to S. Sorokin for the access to the GPS/GLONASS network data in Russia, which allowed us to significantly increase the TEC and PPP visibility in the addressed region. The authors acknowledge the use of AE, H-SYM, Bz, and the proton density data from the OMNI database; the ionosonde data in Yakutsk, Moscow, Mohe from GIRO database; the SuperDARN data, which is a collection of radars funded by national scientific funding agencies of Australia, Canada, China, France, Italy, Japan, Norway, South Africa, United Kingdom, and the United States of America; and the SuperDARN software and web tools available at Virginia Tech.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yeh, K.C.; Liu, C.-H. Radio wave scintillations in the ionosphere. Proc. IEEE 1982, 70, 324–360. [Google Scholar] [CrossRef]
  2. Kintner, P.; Ledvina, B.M.; De Paula, E.R. GPS and ionospheric scintillations. Space Weather 2007, 5. [Google Scholar] [CrossRef]
  3. Aarons, J. Global morphology of ionospheric scintillations. Proc. IEEE 1982, 70, 360–378. [Google Scholar] [CrossRef]
  4. Hysell, D. An overview and synthesis of plasma irregularities in equatorial spread F. J. Atmos. Sol. Terr. Phys. 2000, 62, 1037–1056. [Google Scholar] [CrossRef]
  5. Hlubek, N.; Berdermann, J.; Wilken, V.; Gewies, S.; Jakowski, N.; Wassaie, M.; Damtie, B. Scintillations of the GPS, GLONASS, and Galileo signals at equatorial latitude. J. Space Weather Space Clim. 2014, 4. [Google Scholar] [CrossRef] [Green Version]
  6. Jiao, Y.; Morton, Y.; Taylor, S. Comparative studies of high-latitude and equatorial ionospheric scintillation characteristics of GPS signals. In Proceedings of the Position, Location and Navigation Symposium (PLANS 2014), Monterey, CA, USA, 5–8 May 2014; pp. 37–42. [Google Scholar] [CrossRef]
  7. Jiao, Y.; Morton, Y.J. Comparison of the effect of high-latitude and equatorial ionospheric scintillation on GPS signals during the maximum of solar cycle 24. Radio Sci. 2015, 50, 886–903. [Google Scholar] [CrossRef]
  8. Prikryl, P.; Jayachandran, P.T.; Chadwick, R.; Kelly, T.D. Climatology of GPS phase scintillation at northern high latitudes for the period from 2008 to 2013. Ann. Geophys. 2015, 33, 531–545. [Google Scholar] [CrossRef] [Green Version]
  9. Booker, H.G.; Wells, H.W. Scattering of radio waves by theF-region of the ionosphere. J. Geophys. Res. Space Phys. 1938, 43, 249. [Google Scholar] [CrossRef]
  10. Basu, S.; Kelley, M.C. A review of recent observations of equatorial scintillations and their relationship to current theories of Fregion irregularity generation. Radio Sci. 1979, 14, 471–485. [Google Scholar] [CrossRef]
  11. Woodman, R.F. Spread F—An old equatorial aeronomy problem finally resolved? Ann. Geophys. 2009, 27, 1915–1934. [Google Scholar] [CrossRef]
  12. Skone, S.H. The impact of magnetic storms on GPS receiver performance. J. Geod. 2001, 75, 457–468. [Google Scholar] [CrossRef]
  13. Afraimovich, E.; Astafyeva, E.; Demyanov, V.; Gamayunov, I. Mid-latitude amplitude scintillation of GPS signals and GPS performance slips. Adv. Space Res. 2009, 43, 964–972. [Google Scholar] [CrossRef]
  14. Astafyeva, E.; Yasyukevich, Y.; Maksikov, A.; Zhivetiev, I.V. Geomagnetic storms, super-storms, and their impacts on GPS-based navigation systems. Space Weather 2014, 12, 508–525. [Google Scholar] [CrossRef]
  15. Fridman, S. The formation of small-scale irregularities as a result of ionospheric plasma mixing by large-scale drifts. Planet. Space Sci. 1990, 38, 961–972. [Google Scholar] [CrossRef]
  16. Astafyeva, E.; Afraimovich, E.; Voeykov, S. Generation of secondary waves due to intensive large-scale AGW traveling. Adv. Space Res. 2008, 41, 1459–1462. [Google Scholar] [CrossRef]
  17. Bowman, G. A relationship between polar magnetic substorms, ionospheric height rises and the occurrence of spread-F. J. Atmos. Terr. Phys. 1978, 40, 713–722. [Google Scholar] [CrossRef]
  18. Ma, G.; Maruyama, T. A super bubble detected by dense GPS network at east Asian longitudes. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef]
  19. Demyanov, V.; Yasyukevich, Y.; Ishin, A.; Astafyeva, E.I. Ionospheric super-bubble effects on the GPS positioning relative to the orientation of signal path and geomagnetic field direction. GPS Solut. 2011, 16, 181–189. [Google Scholar] [CrossRef]
  20. Mevius, M.; Tol, S.; Pandey, V.N.; Vedantham, H.K.; Brentjens, M.A.; Bruyn, A.G.; Abdalla, F.; Asad, K.M.B.; Bregman, J.D.; Brouw, W.N.; et al. Probing ionospheric structures using the LOFAR radio telescope. Radio Sci. 2016, 51, 927–941. [Google Scholar] [CrossRef] [Green Version]
  21. Vasilyev, R.V.; Kushnarev, D.; Lebedev, V.; Medvedev, A.; Nevidimov, N.; Ratovsky, K. Perspectives of usage of Irkutsk incoherent scatter radar (IISR) as an imaging riometer and radio-heliograph. J. Atmos. Sol. Terr. Phys. 2013, 105, 273–280. [Google Scholar] [CrossRef]
  22. Reiff, P.; Daou, A.G.; Sazykin, S.; Nakamura, R.; Hairston, M.R.; Coffey, V.; Chandler, M.O.; Anderson, B.J.; Russell, C.T.; Welling, D.; et al. Multispacecraft observations and modeling of the 22/23 June 2015 geomagnetic storm. Geophys. Res. Lett. 2016, 43, 7311–7318. [Google Scholar] [CrossRef] [Green Version]
  23. Astafyeva, E.; Zakharenkova, I.; Alken, P. Prompt penetration electric fields and the extreme topside ionospheric response to the 22–23 June 2015 geomagnetic storm as seen by the Swarm constellation. Earth Planets Space 2016, 68, 1309. [Google Scholar] [CrossRef]
  24. Astafyeva, E.; Zakharenkova, I.; Huba, J.D.; Doornbos, E.; Ijssel, J.V.D. Global Ionospheric and Thermospheric Effects of the June 2015 Geomagnetic Disturbances: Multi-Instrumental Observations and Modeling. J. Geophys. Res. Space Phys. 2017, 122, 11716–11742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Astafyeva, E.; Zakharenkova, I.; Hozumi, K.; Alken, P.; Coïsson, P.; Hairston, M.R.; Coley, W.R. Study of the Equatorial and Low-Latitude Electrodynamic and Ionospheric Disturbances During the 22–23 June 2015 Geomagnetic Storm Using Ground-Based and Spaceborne Techniques. J. Geophys. Res. Space Phys. 2018, 123, 2424–2440. [Google Scholar] [CrossRef] [PubMed]
  26. Singh, R.; Sripathi, S. Ionospheric Response to 22–23 June 2015 Storm as Investigated Using Ground-Based Ionosondes and GPS Receivers Over India. J. Geophys. Res. Space Phys. 2017, 122, 645. [Google Scholar] [CrossRef]
  27. Piersanti, M.; Alberti, T.; Bemporad, A.; Berrilli, F.; Bruno, R.; Capparelli, V.; Carbone, V.; Cesaroni, C.; Consolini, G.; Cristaldi, A.; et al. Comprehensive Analysis of the Geoeffective Solar Event of 21 June 2015: Effects on the Magnetosphere, Plasmasphere, and Ionosphere Systems. Sol. Phys. 2017, 292, 5–169. [Google Scholar] [CrossRef]
  28. Gherm, V.E.; Zernov, N.; Strangeways, H. Effects of diffraction by ionospheric electron density irregularities on the range error in GNSS dual-frequency positioning and phase decorrelation. Radio Sci. 2011, 46. [Google Scholar] [CrossRef]
  29. Andrew, O.-O.A.; Doherty, P.H.; Carrano, C.S.; Valladares, C.E.; Groves, K.M. Impacts of ionospheric scintillations on GPS receivers intended for equatorial aviation applications. Radio Sci. 2012, 47. [Google Scholar] [CrossRef]
  30. Luo, X.; Gu, S.; Lou, Y.; Xiong, C.; Chen, B.; Jin, X. Assessing the Performance of GPS Precise Point Positioning Under Different Geomagnetic Storm Conditions during Solar Cycle 24. Sensors 2018, 18, 1784. [Google Scholar] [CrossRef] [Green Version]
  31. Fremouw, E.J.; Bates, H.F. Worldwide Behavior of Average VHF-UHF Scintillation. Radio Sci. 1971, 6, 863–869. [Google Scholar] [CrossRef]
  32. Fremouw, E.J.; Leadabrand, R.L.; Livingston, R.C.; Cousins, M.D.; Rino, C.L.; Fair, B.C.; Lang, R.A. Early results from the DNA Wideband satellite experiment—Complex-signal scintillation. Radio Sci. 1978, 13, 167–187. [Google Scholar] [CrossRef]
  33. Berkner, L.V.; Wells, H.W. F-region ionosphere-investigations at low latitudes. J. Geophys. Res. Space Phys. 1934, 39, 215. [Google Scholar] [CrossRef]
  34. Reinisch, B.; Haines, D.M.; Bibl, K.; Galkin, I.; Kitrosser, D.F.; Sales, G.S.; Scali, J.L.; Huang, X. Ionospheric sounding in support of over-the-horizon radar. Radio Sci. 1997, 32, 1681–1694. [Google Scholar] [CrossRef]
  35. Potekhin, A.P.; Medvedev, A.; Zavorin, A.V.; Kushnarev, D.S.; Lebedev, V.P.; Shpynev, B.G. Development of diagnostic capabilities of the Irkutsk incoherent scattering radar. Cosmo Res. 2008, 46, 347–353. [Google Scholar] [CrossRef]
  36. Shanmugam, S.; Jones, J.; Macaulay, A.; Van Dierendonck, A. Evolution to modernized GNSS ionoshperic scintillation and TEC monitoring. In Proceedings of the 2012 IEEE/ION Position, Location and Navigation Symposium; Institute of Electrical and Electronics Engineers (IEEE), Myrtle Beach, South Carolina, USA, 23–26 April 2012; pp. 265–273. [Google Scholar]
  37. Spogli, L.; Alfonsi, L.; De Franceschi, G.; Romano, V.; Aquino, M.H.O.; Dodson, A. Climatology of GPS ionospheric scintillations over high and mid-latitude European regions. Ann. Geophys. 2009, 27, 3429–3437. [Google Scholar] [CrossRef]
  38. Hofmann-Wellenhof, B.; Lichtenegger, H.; Collins, J. Global Positioning System; Springer Science and Business Media LLC: Berlin, Germany, 2001. [Google Scholar]
  39. Jin, S.; Feng, G.; Gleason, S. Remote sensing using GNSS signals: Current status and future directions. Adv. Space Res. 2011, 47, 1645–1653. [Google Scholar] [CrossRef]
  40. Afraimovich, E.L.; Astafyeva, E.I.; Demyanov, V.; Edemskiy, I.K.; Gavrilyuk, N.S.; Ishin, A.; Kosogorov, E.A.; Leonovich, L.A.; Lesyuta, O.S.; Palamartchouk, K.; et al. A review of GPS/GLONASS studies of the ionospheric response to natural and anthropogenic processes and phenomena. J. Space Weather Space Clim. 2013, 3, A27. [Google Scholar] [CrossRef] [Green Version]
  41. Dow, J.M.; Neilan, R.E.; Rizos, C. The International GNSS Service in a changing landscape of Global Navigation Satellite Systems. J. Geod. 2009, 83, 191–198. [Google Scholar] [CrossRef]
  42. Яcюкeвич, Ю.; Yasyukevich, Y.; Пeрeвaлoвa, Н.; Perevalova, N.; Becнин, A.; Vesnin, A. SibNet—Siberian Global Navigation Satellite System Network: Current state. Sol. Terr. Phys. 2018, 4, 63–72. [Google Scholar] [CrossRef]
  43. Jayachandran, P.T.; Langley, R.B.; MacDougall, J.W.; Mushini, S.C.; Pokhotelov, D.; Hamza, A.M.; Mann, I.R.; Milling, D.K.; Kale, Z.C.; Chadwick, R.; et al. Canadian High Arctic Ionospheric Network (CHAIN). Radio Sci. 2009, 44. [Google Scholar] [CrossRef]
  44. Bruyninx, C.; Habrich, H.; Söhne, W.; Kenyeres, A.; Stangl, G.; Völksen, C. Enhancement of the EUREF Permanent Network Services and Products. In International Association of Geodesy Symposia; Springer Science and Business Media LLC: Berlin, Germany, 2011; Volume 136, pp. 27–34. [Google Scholar]
  45. Maletckii, B.; Yasyukevich, Y.; Vesnin, A. Wave Signatures in Total Electron Content Variations: Filtering Problems. Remote Sens. 2020, 12, 1340. [Google Scholar] [CrossRef] [Green Version]
  46. Schaer, S. Mapping and Predicting the Earth’s Ionosphere Using the Global Positioning System. Ph.D. Dissertation, Astronomical Institute, University of Berne, Berne, Switzerland, 25 March 1999. [Google Scholar]
  47. Yasyukevich, Y.; Kiselev, A.V.; Zhivetiev, I.V.; Edemskiy, I.K.; Syrovatskii, S.V.; Maletckii, B.M.; Vesnin, A.M. SIMuRG: System for Ionosphere Monitoring and Research from GNSS. GPS Solut. 2020, 24, 1–12. [Google Scholar] [CrossRef]
  48. Zhou, F.; Dong, D.; Li, W.; Jiang, X.; Wickert, J.; Schuh, H. GAMP: An open-source software of multi-GNSS precise point positioning using undifferenced and uncombined observations. GPS Solut. 2018, 22, 33. [Google Scholar] [CrossRef]
  49. Yasyukevich, Y.; Astafyeva, E.; Padokhin, A.; Ivanova, V.A.; Syrovatskii, S.; Podlesnyi, A. The 6 September 2017 X-Class Solar Flares and Their Impacts on the Ionosphere, GNSS, and HF Radio Wave Propagation. Space Weather 2018, 16, 1013–1027. [Google Scholar] [CrossRef]
  50. Mathews, G.J.; Towheed, S.S. NSSDC OMNIWeb: The first space physics WWW-based data browsing and retrieval system. Comput. Netw. ISDN Syst. 1995, 27, 801–808. [Google Scholar] [CrossRef]
  51. Thebault, E.; Finlay, C.C.; Beggan, C.D.; Alken, P.; Aubert, J.; Barrois, O.; Bertrand, F.; Bondar, T.; Boness, A.; Brocco, L.; et al. International Geomagnetic Reference Field: The 12th generation. Earth Planets Space 2015, 67, 68. [Google Scholar] [CrossRef]
  52. Vasilyev, R.V.; Globa, M.; Kushnarev, D.; Medvedev, A.; Ratovsky, K.; Globa, M. Spectral characteristics of ionospheric scintillations of VHF radiosignal near magnetic zenith. J. Atmos. Sol. Terr. Phys. 2017, 160, 48–55. [Google Scholar] [CrossRef]
  53. Crane, R.K. Spectra of ionospheric scintillation. J. Geophys. Res. Space Phys. 1976, 81, 2041–2050. [Google Scholar] [CrossRef]
  54. Kelley, M.C.; Vickrey, J.F.; Carlson, C.W.; Torbert, R. On the origin and spatial extent of high-latitudeFregion irregularities. J. Geophys. Res. Space Phys. 1982, 87, 4469. [Google Scholar] [CrossRef]
  55. Kelley, M.C. 10 Instabilities and Structure in the High-Latitude Ionosphere. Int. Geophys. 2009, 96, 469–544. [Google Scholar] [CrossRef]
  56. Kagan, L.M.; Myasnikov, E.N.; Kosolapenko, V.I.; Kryazhev, V.A.; Cheremnyj, V.A.; Persson, M.A. F-layer irregularities’ formation at auroral latitudes: Radio wave scintillation and EISCAT observations. J. Atmos. Terr. Phys. 1995, 57, 917–928. [Google Scholar] [CrossRef]
  57. Bowman, G.G. Periodicities and spread-F associated with largescale travelling ionospheric disturbances. Indian J. Radio Space Phys. 1996, 25, 93–100. Available online: http://nopr.niscair.res.in/handle/123456789/35645 (accessed on 16 April 2020).
  58. Bowman, G.G.; Mortimer, I.K. Some characteristics of large-scale travelling ionospheric disturbances and a relationship between the F2 layer height rises of these disturbances and equatorial pre-sunrise events. Ann. Geophys. 2010, 28, 1419–1430. [Google Scholar] [CrossRef] [Green Version]
  59. Farley, D.T.; Balsey, B.B.; Woodman, R.F.; McClure, J.P. Equatorial spreadF: Implications of VHF radar observations. J. Geophys. Res. Space Phys. 1970, 75, 7199–7216. [Google Scholar] [CrossRef]
  60. Hysell, D.; Larsen, M.; Fritts, D.; Laughman, B.; Sulzer, M. Major upwelling and overturning in the mid-latitude F region ionosphere. Nat. Commun. 2018, 9, 3326. [Google Scholar] [CrossRef]
  61. Newell, P.T.; Sotirelis, T.; Wing, S. Diffuse, monoenergetic, and broadband aurora: The global precipitation budget. J. Geophys. Res. Space Phys. 2009, 114. [Google Scholar] [CrossRef]
  62. Chisham, G.; Lester, M.; Milan, S.E.; Freeman, M.P.; Bristow, W.A.; Grocott, A.; McWilliams, K.A.; Ruohoniemi, J.M.; Yeoman, T.K.; Dyson, P.L.; et al. A decade of the Super Dual Auroral Radar Network (SuperDARN): Scientific achievements, new techniques and future directions. Surv. Geophys. 2007, 28, 33–109. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Locations of Irkutsk Incoherent Scatter Radar/Irkutsk ionosonde DPS-4 (red star) and GPS/GLONASS receivers (blue dots).
Figure 1. Locations of Irkutsk Incoherent Scatter Radar/Irkutsk ionosonde DPS-4 (red star) and GPS/GLONASS receivers (blue dots).
Remotesensing 12 01579 g001
Figure 2. Evolution of the 22 June 2015 geomagnetic storm from the OMNI data. Panel a: H-SYM index (black line) and the interplanetary magnetic field (IMF) Bz component (blue line). Panel b: AE-index. Panel c: solar wind proton density (black line) and spread-F width (red line). Spread-F data are provided for the Irkutsk ionosonde (52°N, 104°E).
Figure 2. Evolution of the 22 June 2015 geomagnetic storm from the OMNI data. Panel a: H-SYM index (black line) and the interplanetary magnetic field (IMF) Bz component (blue line). Panel b: AE-index. Panel c: solar wind proton density (black line) and spread-F width (red line). Spread-F data are provided for the Irkutsk ionosonde (52°N, 104°E).
Remotesensing 12 01579 g002
Figure 3. Experiment geometry. Panel (a) shows the GPS IPP (black lines) at 300 km and 10–20 min total electron content (TEC) variations (superimposed red lines) within 19:30–20:30 UT on 22 June 2015. Panel (b) shows the geometry in details. Black lines show the GPS IPP at 19:30–20:30 UT. Gray lines show the Cygnus A ionospheric pierce point at 100, 250, and 400 km. On both panels, the filled grey region is the ionosonde approximate radiation pattern, the filled blue region is the IISR approximate radiation pattern. Vertical segment in the right bottom part of the panel (a) shows 0.5 TECU variations.
Figure 3. Experiment geometry. Panel (a) shows the GPS IPP (black lines) at 300 km and 10–20 min total electron content (TEC) variations (superimposed red lines) within 19:30–20:30 UT on 22 June 2015. Panel (b) shows the geometry in details. Black lines show the GPS IPP at 19:30–20:30 UT. Gray lines show the Cygnus A ionospheric pierce point at 100, 250, and 400 km. On both panels, the filled grey region is the ionosonde approximate radiation pattern, the filled blue region is the IISR approximate radiation pattern. Vertical segment in the right bottom part of the panel (a) shows 0.5 TECU variations.
Remotesensing 12 01579 g003
Figure 4. Effects of ionospheric irregularities on different-band radio signals. Panel a presents the spread-F at the Irkutsk ionosonde; panel b displays the σφ for the L1 GPS signal; panel c reflects the Cygnus A radio amplitude scintillations from the Irkutsk Incoherent Scatter Radar (IISR) data; panel d shows spread-F. Panel d shows data for the ionosondes in Yakutsk (YA462, 62.0°N, 129.7°E; blue line), Mohe (MH453; 52.0°N, 122.5°E; black line), and Moscow (MO155; 55.47°N, 37.3°E; red line) during the 22 June 2015 magnetic storm. The horizontal line on panel d shows 0.5 MHz. The black lines and the histogram on panels a-c show the 22 June 2015 data, the gray line on the panel c shows the 24 June 2015 data (reference day). The dashed-dotted line in panel c shows the angle γ between the “IISR—Cygnus A” LOS and the Earth magnetic field. The dashed lines mark 19:30 UT, 20:15 UT, 21:30 UT.
Figure 4. Effects of ionospheric irregularities on different-band radio signals. Panel a presents the spread-F at the Irkutsk ionosonde; panel b displays the σφ for the L1 GPS signal; panel c reflects the Cygnus A radio amplitude scintillations from the Irkutsk Incoherent Scatter Radar (IISR) data; panel d shows spread-F. Panel d shows data for the ionosondes in Yakutsk (YA462, 62.0°N, 129.7°E; blue line), Mohe (MH453; 52.0°N, 122.5°E; black line), and Moscow (MO155; 55.47°N, 37.3°E; red line) during the 22 June 2015 magnetic storm. The horizontal line on panel d shows 0.5 MHz. The black lines and the histogram on panels a-c show the 22 June 2015 data, the gray line on the panel c shows the 24 June 2015 data (reference day). The dashed-dotted line in panel c shows the angle γ between the “IISR—Cygnus A” LOS and the Earth magnetic field. The dashed lines mark 19:30 UT, 20:15 UT, 21:30 UT.
Remotesensing 12 01579 g004
Figure 5. Average precise point positioning (PPP) error (obtained by GAMP software) vs. latitude. The latitude bin size is 2.5°. The longitude region was 0°…180° (European-Asian sector, upper panels a-b), −180°…0° (American-Atlantic sector, bottom panels c-d). The dashed line on panels b and d marks 18:30 UT. Panels a and c show the data for the reference day (172 DOY), while panels b and d show the data for the 22 June 2015 storm (173 DOY).
Figure 5. Average precise point positioning (PPP) error (obtained by GAMP software) vs. latitude. The latitude bin size is 2.5°. The longitude region was 0°…180° (European-Asian sector, upper panels a-b), −180°…0° (American-Atlantic sector, bottom panels c-d). The dashed line on panels b and d marks 18:30 UT. Panels a and c show the data for the reference day (172 DOY), while panels b and d show the data for the 22 June 2015 storm (173 DOY).
Remotesensing 12 01579 g005
Figure 6. Maps for the electric potential from the SuperDARN data for 18:30 UT (a), 18:40 UT (b), 19:00 UT (c), 19:30 UT (d), 20:00 UT (e), 20:30 UT (f), 21:00 UT (g), 21:30 UT (h). Green dashed curve shows the Heppner–Maynard boundary.
Figure 6. Maps for the electric potential from the SuperDARN data for 18:30 UT (a), 18:40 UT (b), 19:00 UT (c), 19:30 UT (d), 20:00 UT (e), 20:30 UT (f), 21:00 UT (g), 21:30 UT (h). Green dashed curve shows the Heppner–Maynard boundary.
Remotesensing 12 01579 g006
Figure 7. Maps of 10–20 min TEC variations for 18:30 UT (a), 19:00 UT (b), 19:30 UT (c), 20:00 UT (d), 20:30 UT (e), 21:00 UT (f). Scale is shown at the bottom in TECU units.
Figure 7. Maps of 10–20 min TEC variations for 18:30 UT (a), 19:00 UT (b), 19:30 UT (c), 20:00 UT (d), 20:30 UT (e), 21:00 UT (f). Scale is shown at the bottom in TECU units.
Remotesensing 12 01579 g007
Figure 8. Maps of 20–60 min TEC variations for 18:30 UT (a), 19:00 UT (b), 19:30 UT (c), 20:00 UT (d), 20:30 UT (e), 21:00 UT (f). Scale is shown at the bottom in TECU units.
Figure 8. Maps of 20–60 min TEC variations for 18:30 UT (a), 19:00 UT (b), 19:30 UT (c), 20:00 UT (d), 20:30 UT (e), 21:00 UT (f). Scale is shown at the bottom in TECU units.
Remotesensing 12 01579 g008

Share and Cite

MDPI and ACS Style

Yasyukevich, Y.; Vasilyev, R.; Ratovsky, K.; Setov, A.; Globa, M.; Syrovatskii, S.; Yasyukevich, A.; Kiselev, A.; Vesnin, A. Small-Scale Ionospheric Irregularities of Auroral Origin at Mid-latitudes during the 22 June 2015 Magnetic Storm and Their Effect on GPS Positioning. Remote Sens. 2020, 12, 1579. https://doi.org/10.3390/rs12101579

AMA Style

Yasyukevich Y, Vasilyev R, Ratovsky K, Setov A, Globa M, Syrovatskii S, Yasyukevich A, Kiselev A, Vesnin A. Small-Scale Ionospheric Irregularities of Auroral Origin at Mid-latitudes during the 22 June 2015 Magnetic Storm and Their Effect on GPS Positioning. Remote Sensing. 2020; 12(10):1579. https://doi.org/10.3390/rs12101579

Chicago/Turabian Style

Yasyukevich, Yury, Roman Vasilyev, Konstantin Ratovsky, Artem Setov, Maria Globa, Semen Syrovatskii, Anna Yasyukevich, Alexander Kiselev, and Artem Vesnin. 2020. "Small-Scale Ionospheric Irregularities of Auroral Origin at Mid-latitudes during the 22 June 2015 Magnetic Storm and Their Effect on GPS Positioning" Remote Sensing 12, no. 10: 1579. https://doi.org/10.3390/rs12101579

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop