Strong turbulent flow in the sub-auroral region in the Antarctic can deteriorate satellite-based navigation signals

doi:10.21203/rs.3.rs-5004088/v1

Download PDF

Article

Strong turbulent flow in the sub-auroral region in the Antarctic can deteriorate satellite-based navigation signals

https://doi.org/10.21203/rs.3.rs-5004088/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

In the subauroral zone at the boundary of the auroral oval in the evening and night hours during geomagnetic disturbances, a narrow (about 1º-2º) and extended structure (several hours in longitude) is formed. It is known as a polarization jet (PJ) or the sub-auroral ion drift (SAID). The PJ/SAID is a fast westward ion drift and is one of the main signatures of a geomagnetic disturbance in the subauroral ionosphere at the altitudes of the F-layer, when the geomagnetic AE index reaches more than 500 nT. Plasma speed in the PJ/SAID can reach several kilometres per second, and the size of plasma irregularities inside it can reach scales from tens of meters to several hundred meters. Such high velocities and structured plasma can affect trans-ionospheric radio signals and lead to scintillations in the received signal. We show that at the moment of auroral activity intensification, an increase in the magnitude of phase scintillation index (σ_ϕ) as well as loss of satellite signals lock were observed in the region of the PJ/SAID equatorward of the auroral oval over Dronning Maud Land (Queen Maud Land) in Antarctica. We find that fluctuations inside the PJ/SAID can lead to serious deterioration of radio communication or navigational services. We emphasize the importance of considering the geometry of the beam passing from the GNSS satellite to the receiver on the ground. We highlight the mutual contribution of the PJ/SAID and the diffuse aurora boundary, which are almost impossible to separate in practice. Our results demonstrate the importance of considering the subauroral zone, where very dynamic plasma formations can occur with a strong flow and various-scale irregularities inside that lead to serious interference in satellite communications.

Physical sciences/Astronomy and planetary science/Space physics/Magnetospheric physics

Earth and environmental sciences/Space physics

Earth and environmental sciences/Space physics/Aurora

Physical sciences/Physics/Space physics

The subauroral zone is located between the projection of the position of the plasmapause along the magnetic field lines of the geomagnetic field and the equatorial boundary of the auroral oval. During increasing geomagnetic activity, the position of the subauroral zone changes rapidly. Phenomena observed in the subauroral zone are unique. For example, it is precisely in the subauroral region where narrow streams of fast westward subauroral ion drifts occur near the projection of the plasmapause at altitudes within the F-layer of the ionosphere. These streams are known as polarization jets (PJ) or Sub-Auroral Ion Drifts (SAID) [Galperin et al., 1974; Spiro et al., 1979] and are most noticeable during storms/substorms against the background of large-scale plasma convection. The decrease in plasma density in the ionospheric F-layer within PJ/SAID significantly affects the conditions for the propagation of shortwave radio emissions, indicating the practical importance of studying PJ/SAID.

Various studies involving different space missions, numerical modelling, and radar research have allowed for the identification of key large-scale features in the formation and development of PJ/SAID [e.g., Khalipov et al., 2001; Anderson et al., 1991, 1993; Karlsson et al., 1998; Figueiredo et al., 2004; Mishin and Puhl-Quinn, 2007; Mishin, 2023; Horvath and Lovell, 2023a]. Previous studies have also helped to determine the time and location where subauroral ion drifts are most likely to occur and establish their connection with geomagnetic disturbances. In recent studies using data from the Norwegian microsatellite NorSat-1, equipped with a high-frequency Langmuir probe system, the small-scale structure of PJ/SAID has been investigated the first time [Sinevich et al., 2021, 2022, 2023]. Thus, it is now well-known that within PJ/SAID, there exist plasma irregularities of various scales that can significantly impact the propagation of radio waves, radio communication, and cause scintillations in navigation signals.

Since it is now known that PJ/SAID can be a source of various inhomogeneous structures in the ionosphere, it became important to understand whether it is possible to trace their effects on the quality of trans-ionospheric radio signals, and hence positioning with the Global Navigation Satellite Systems (GNSS), and to separate these effects from those originating from irregularities caused by the energetic particle precipitation during auroras, especially in the southern hemisphere. The polar and subpolar regions in the southern hemisphere, particularly around the Antarctic, are the least studied regions to date. This is mainly because the influence of various ionospheric effects on the propagation of radio signals, including those from navigation satellites (GPS, GLONASS, Galileo, BeiDou), has been studied primarily for the northern hemisphere, where there is a dense network of ground-based receivers, radars, ionosondes, and other instruments. Recently, the pace of development of international and regional programs for the exploration of Antarctica has been increasing, along with the need to ensure the reliable functioning of radio-technical systems in the Antarctic region. It is important to note that the subauroral region is located approximately near the coast of Antarctica, where there are a significant number of polar stations (some with airstrips for delivering supplies and personnel to Antarctica) and where there is relatively dense ship traffic. Therefore, it is necessary to understand the processes occurring in the ionosphere in order to avoid or predict potential issues with radio communication and the degradation of positioning accuracy in navigation systems.

In this study, we explore for the first time how strong turbulent flows occurring in the subauroral ionosphere of the southern hemisphere during the PJ/SAID development affect slips (i.e., block of radio signals from the GNSS satellites) in satellite-based navigation systems and lead to scintillations in the received signal on the ground. To achieve this, a unique multi-instrumental approach is employed involving a large number of instruments and data sources. The uniqueness of the considered case is that at that moment of time all instruments were working, and we were able to comprehensively study the geomagnetic disturbance and subsequent scintillations and loss of signals from the GNSS satellites. We used data from three space missions (NorSat-1, DMSP, Swarm), data from two receivers installed in Antarctica, and data from magnetometers. Such a rich set of tools made it possible to study in detail the effects in the subauroral region during the presence of PJ/SAID and its impact on GNSS operations, as well as obtaining new results relevant not only to the Antarctic but also to the Arctic region since many processes occur simultaneously in the polar zones during varying geomagnetic activity.

We chose conditions close to the spring equinox on March 18, 2018, during a period of low solar activity, when a moderate geomagnetic disturbance was observed with a maximum of Kp = 6, SYM-H = − 53 nT and Auroral Electrojet Index (AE) was up to 1600 nT. In the evening local sector there was an outstanding combination of instruments available that could identify the presence of PJ/SAID and this was the determining factor for focusing on this event. A significant feature of this work is the rich amount of satellite data available, which has shed light on important details associated with structures in the subauroral region. Figure 1 presents UV snapshots of auroral emissions from the SSUSI instrument installed onboard DMSP satellite F-18 when it entered a new orbit from the south pole at 20:45 UT towards the equator. In Fig. 1, the trajectory of the F-18 satellite is shown with black line and the part with PJ/SAID is marked with a thick line. Figure 1 also shows in a similar manner the overflights of other LEO satellites: Swarm B, C and DMSP F-17, as well as NorSat-1 (colour coded). More details on how the PJ/SAID has been identified from the measurements is given in the Methods section. In the middle of the figure, we indicate locations of two GNSS Ionospheric Scintillation and TEC Monitor (GISTM) receivers at the Troll (TRL) and SANAE IV (SAN) research stations and two magnetometers at the SANAE IV and Neumayer stations (shown with the red, black and white stars respectively). Phase scintillation indices σ_ϕ observed with the GISTM receivers are presented at locations corresponding to ionospheric piercing points at the height of 350 km (circles for TRL receiver and squares for the SNA receiver), and σ_ϕ intensity is colour-coded according to the colour bar scale. The values of the phase scintillation index are given for the moment of time 21:25 UT, which is four minutes after the PJ/SAID was detected according to data from the NorSat-1 satellite (more images, as well as videos showring the dynamics of changes of the indices are given in the Supplementary Information). Since the SSUSI images are taken dynamically (about 20 minutes of shooting is required), the image of the boundaries of the precipitating particles in the UV range shows the average location of the diffuse and discrete auroras. However, even now it is clear from this figure that the area of increased values of the phase scintillation index up to 0.6 rad coincides with the area of the intense discrete aurora, as well as on the equatorial boundary of the diffuse aurora, where, according to satellite data and theory, the PJ/SAID is located. During detailed analysis of plots and results, our attention was focused on selected GNSS satellites (Troll station, GLONASS system: PRN 45 at 21:23 UT, PRN 49 at 21:50 and 21:52 UT, PRN 46 at 22:39, 22:46 and 22:50 UT; SANAE IV station, GPS system: PRN 1 at 21:30 UT, PRN 14 at 22:59 and 23:01 UT and the Galileo system at TRL, PRN 73 two hours later at 01:19 UT), that showed extreme values of the phase scintillation index. In the usual practice of scintillation data analysis such values should be removed from the data set because the corresponding lock time values are less than 240s. However for us, that means that at some particular moment, the receiver cannot track the satellite. In our case, such satellites were united by the fact that they were: 1) lost by the receiver and the positioning was disrupted, such that a new tracking count started; 2) in the area where the PJ/SAID was located; 3) lost at low elevation angles (less than 25°). Note that in order to separate different satellite systems on one plot we use PRN + 0 for GPS, PRN + 40 for GLONASS, and PRN + 70 for Galileo systems.

The analysis of other times and days showed that satellite PRN 1 of the SANAE IV receiver has a constant phase tracking failure in the same place. Which, apparently, may be associated with local features of the conditions near the receiver. The same is true for the PRN 49 satellite of the Troll receiver, where approximately in the same place and at times shifted due to the orbits of the satellites, a tracking failure was also observed. This again indicates interference associated with the area around the satellite, rather than with the state of the environment along the signal path from the GNSS satellite to the receiver on the ground. The remaining four satellites (PRN 45, 46, 14 and 73) indicate a disruption associated with changes in the ionosphere. In addition, in the area where the tracking of satellite PRN 45 is disrupted, after 1.5 hours the tracking of satellite PRN 14 also fails.

When we look at the dynamics in detail, we see that at 21:08 UT the phase scintillation indices σ_ϕ begin to increase in the line parallel to the geomagnetic isoline approximately at 62–63° MLAT above the receivers (see details in Supplementary Information). The enhancement of σ_ϕ above stations and more poleward is associated with the flow shear instabilities that are modulated by the auroral particle precipitation, as shown by Skjæveland et al. (2021) who also considered this event. It is noteworthy that the amplitude scintillation index S4 is more sensitive in the zone of discrete aurora, but we also observe increased values S4 from satellites passing through the zone of the equatorial boundary of a diffuse aurora where PJ/SAID is usually observed. At 21:17 UT, according to magnetometer data at SANAE IV and Neumayer (see Fig. 6 in Methods), a sharp increase in the northern component N of the magnetic field begins, and at the same time the values of σ_ϕ above the receivers increase by more than 0.8 rad (PRN satellites 14 = 1.67 rad, 3 = 1.49 rad, 32 = 0.88 rad, 31 = 1.63 rad, 47 = 1.72 rad, 46 = 1.77 rad, 75 = 2.03 rad, 12 = 0.82 rad). In the next few minutes, even more polar satellites show values σ_ϕ > 1 rad. A subsequent decrease in the N component at 21:20 UT at both magnetometers indicates an increase in the westward electrojet associated with increased electric field or conductivity. The behaviour of the S4 index is not so strongly influenced by the expansion phase of the substorm. However, we observe an increase in the S4 value from the satellites passing through the edges of the auroral oval. This generally corresponds to the statistical study of Spogli et al. [2009], where it was shown that the amplitude scintillations index gets enhanced in correspondence with the boundaries of the auroral oval. Thus, satellite PRN 45 already had high values σ_ϕ and S4 even before the loss of lock at 21:23 UT (elevation angle 𝛼 = 13.4°). The GPS PRN 1 satellite data obtained with TRL and SNA receivers also show an increase in S4 in the position of the PJ/SAID (𝛼 is about 13°). The PRN 46 satellite had low σ_ϕ index values of approximately 0.25 rad and a bit higher S4 index of about 0.27 before entering the expected PJ/SAID zone. The latest satellite measurements confirming the presence of a PJ/SAID are by Swarm B at 22:20 UT. Then, in the PJ/SAID passage zone, communication with satellite 46 is lost 19 minutes later at 22:39 UT (𝛼 = 23.49°). The loss of lock is occurring again at 22:46 UT (𝛼 = 19.83°), and in the period from 22:50 (𝛼 = 17.76°) to 22:54 UT (𝛼 = 15.7°). Afterwards, although the minimum cut-off elevation angle of 10° chosen in the analysis had not yet been reached, communication with the satellite was no longer established. GPS satellite 14 shows an increase in the values of the S4 and σ_ϕ indices after 22:35 UT as the satellite approaches the area where the PJ/SAID was detected. After passing through a slightly disturbed region at 22:59 UT there was a loss of lock for the L1CA signal from the SANAE IV receiver (𝛼 = 12.42°). Later at night Galileo satellite 73 experienced three times the loss of signal lock at 01:19, 01:23 and 01:28 UT, corresponding to approximately − 57°-58° MLAT (𝛼 = 15.85°, 14.40°, 12.59°). Measurements from the driftmeter as well as UV images from SSUSI onboard DMSP F-17 and F-18 satellites confirm the still stable presence of the PJ/SAID structure around − 60° MLAT at 01:56 and 02:00 UT, respectively. This indicates that the structure is long-lasting, dynamic (with poleward and equatorward movements), and also varying in thickness from 1°-3°. All satellites considered experienced an increase in the phase and amplitude scintillation indices, as well as loss of signal lock in areas where their signal passed through the PJ/SAID structure. This again indicates the strong influence of the PJ/SAID on GNSS signals.

In this work, we studied how the PJ/SAIDs affect slips, phase and amplitude scintillations of the received radio signal from GNSS satellites in the subauroral region in the southern hemisphere near and over Dronning Maud Land (Queen Maud Land), Antarctica. To explicitly determine whether the PJ/SAID was present at a given location, we used data from several satellites. Data received from satellites DMSP F-17, DMSP F-18, Swarm C, Swarm B and NorSat-1 confirmed the presence of the PJ/SAID in the subauroral region during the geomagnetic substorm of March 18, 2018. According to the analysis, PJ/SAID existed during a considered event at least from 19 UT March 18, 2018 to 02 UT March 19, 2018. Based on their research Anderson et al. [1991] concluded that a typical PJ/SAID lifetime is between 30 minutes and 3 hours. The mechanisms behind the PJ/SAID formation allow the existence of these phenomena for a longer period during geomagnetic storms due to a continuous energy input from the magnetosphere/ionosphere [Anderson et al., 1993; Southwood and Wolf, 1978; Mishin and Puhl-Quinn, 2007]. It is also likely that the PJ/SAID can fade and reappear during one geomagnetic event, depending on geomagnetic activity variations. According to various satellite data, the PJ/SAID considered in this study is located between 53° and 65° MLAT and from 16 to 02 MLT, which is in good agreement with Spiro at al. [1979] and Karlsson et al. [1998] who, based on statistics, claim that typical invariant latitudes of the PJ/SAID are 55°-65° and that the PJ/SAID can be stretched along latitude for more than ~ 6 hours of MLT in the midnight and evening sectors. In Fig. 1 one can see that PJ/SAID detected by DMSP F-18, DMSP F-17 and Swarm C satellites is located at higher latitudes (~ 59°, ~ 60° and ~ 65° MLAT respectively) than PJ/SAID detected by Swarm B (~ 55° MLAT at 22:20 UT and ~ 57° MLAT at 20:48 UT) and NorSat-1 (~ 54° MLAT). This is because the mean magnetic latitude of PJ/SAID lowers with an increase of MLT moving from the evening to post-midnight sector according to Plate 1 in Karsson et al. [1998] and Foster, Vo [2002] for SAPS (Subauroral Polarization Stream), which most likely has similar features to PJ/SAID. The latitudinal width of the considered PJ/SAID is around 1° which is as well the characteristic of PJ/SAID and most researchers agree that typical PJ/SAID width is 1°-2° [Spiro et al., 1979; Anderson et al., 1991; Karlsson et al., 1998].

We decided to check whether our case was perhaps only an isolated event. We can confirm that it is a regular phenomenon. According to the DMSP satellite data, the PJ/SAID was also present during days close to the selected event. One of the cases on the 25th of March is shown in Fig. 2. Here the geomagnetic disturbance was weaker, and the equatorial boundary of the diffuse aurora was located above the receivers. The thick areas along the satellite's flight path show how dynamic the PJ/SAID structure is over time. The 10-minute difference between the two flybys results in the PJ/SAID moving more poleward about 1° MLAT. The poleward shift is associated with the general dynamics of displacing boundaries of the auroral oval with a decrease in geomagnetic activity [Zossi et al., 2020]. However, at the boundary of the diffuse aurora, where the narrow region of the PJ/SAID begins, we also see an increase in the phase and amplitude scintillation index (see Fig. 2 and Supplementary Information). There is also a loss of signal lock of some satellites observed with the Troll receiver, GPS PRN 3 at 23:04 UT (𝛼 = 15.13°, -57.7° MLAT); GLONASS PRN 55 from 19:44 − 19:55 UT (𝛼 = 19.77° and 25.6°, -59.6° and − 60.4° MLAT), PRN 46 at 23:29 UT (𝛼 = 25.84°, -60.0° MLAT) and the receivers stops tracking the satellite after that, PRN 58 at 23:42 UT (𝛼 = 11.3°, -63.8° MLAT). These satellites are characterized by a higher elevation angle during loss of lock than the considered March 18th event. This corresponds to the PJ/SAID position which is more poleward on March 25th than on March 18th. Thus, it is clear that the increase in scintillation indices during the development of the PJ/SAID during the geomagnetic event of March 18th, 2018, is not an “isolated” event, but it is a regular phenomenon in the high-latitude region. It is of course worth noting that not all GNSS satellites are affected.

Plasma structuring, high-speed flow and irregularities in the ionosphere can cause degradation of communication and precise positioning performance, cycle slips and loss of lock [Kintner and Ledvina, 2005; Oskavik et al., 2015; Jin and Oskavik, 2018; Chernyshov et el., 2020b]. It is known that rapid changes in the ionosphere can lead to scintillations and fluctuation in phase and amplitude measurements [e.g., Kintner and Ledvina, 2005]. At high latitudes, both phase and amplitude scintillations of radio signals interfere with communication and navigation as ionospheriс irregularities usually intensify during geomagnetic activity periods [Spogli et a., 2009]. The dynamics of the subauroral ionosphere during a substorm is driven by the interaction of the ring current and the plasmasphere/plasmapause, resulting in increased fluxes of the westward convection, that is, to the emergence of the PJ/SAID where areas of irregular plasma density troughs and enhanced electron temperature are located, leading to scintillations at subauroral latitudes [Mishin and Blaunstein 2008]. These are among the fastest flows observed anywhere in the ionosphere, so scintillations observed in this region should have different properties [Wang et al., 2018]. In the region of the PJ/SAID, we observe an increase in the indices of both phase and amplitude scintillations. The phase scintillation index is sensitive to large (larger than the Fresnel scale) and small-scale irregularities along the signal path, while the S4 index is an indicator of small-scale irregularities below Fresnel scales. The Fresnel scale at n = 1 is about 365 m at a height of 350 km for L1 frequency) [Kintner et al., 2007]. Variations in these two indices indicate the presence of plasma irregularities of different scales.

The PJ/SAIDs themselves are very structured, having within them irregularities of different scales, which are clearly visible from data from the NorSat-1 satellite, equipped with a high-resolution Langmuir probe system. In this work, inside the considered PJ/SAID small-scale irregularities in plasma parameters are found (see Methods, Fig. 5). The spatial size of irregularities of electron density is up to several hundreds of meters which does not exceed the background level of subauroral plasma during disturbed geomagnetic conditions. On the contrary, irregularities in the electron temperature inside the PJ/SAID have significantly smaller minimum sizes (dozens of meters) and are more intensive inside than outside of PJ/SAID.

Although the PJ/SAID has a rather narrow structure (1°-2° in latitude), strong plasma flows and large gradients of plasma parameters can significantly impact the passage of radio waves. In the considered example of a geomagnetic event, this is manifested in increased values of both scintillation indices (but without a significant increase in the S4 values). However, it is necessary to consider the geometry of the study case. Since we are considering the Antarctic region, the navigation signal arrives at the receiver at a lower elevation angle compered to mid-latitudes, as the inclination of the GNSS satellites is relatively small. Therefore, the radio signal from the "equatorward" satellites (as illustrated in Fig. 3) passes not only through the PJ/SAID region with its manifestation in the ionospheric F-region but also through the region of diffuse auroras.

Therefore, it is quite challenging to separate one effect from another or to consider them in context from each other. Our results show that it is a combination of two effects together from the PJ/SAID and the diffuse aurora boundary. Note that satellite signals passing only through areas of diffuse aurora, as well as the boundary between diffuse and discrete aurora, were not that affected and there was no loss of tracking of the GNSS satellites. A large network of receivers in future is necessary to identify how scintillation is related to subauroral density structures.

In this study, we were the first to use a multi-instrumental approach to study the influence of the PJ/SAID on the occurrence of scintillations and inhomogeneous structures of various scales (from dozens of meters up to several hundreds of meters) in the ionosphere in the Antarctic subauroral region. Our multi-instrumental approach is based on the use of various satellites (NorSat-1, DMSP, Swarm) and ground-based data (GISTM receivers, magnetometers). This technique allowed for a comprehensive consideration of geomagnetic events, which made it possible to draw firm and reasonable conclusions. A multi-instrumental complex study confirmed the presence of a PJ/SAID in the region of two receivers under consideration. We presented the influence of the PJ/SAID effect on the propagation of the trans-ionospheric radio signals during two geomagnetic events, where we observed an increase in the phase and amplitude scintillation indices, as well as loss of the signal lock.

Our study shows that in practical applications (such as maritime navigation) during geomagnetic disturbances, it is necessary to consider not only the precipitation of charged particles in the auroral oval but also the development of a PJ/SAID in the subauroral region in the evening and night local hours, which can lead to loss of communication with the satellite. Further studies should be carried out also for the northern hemisphere. A distinctive feature of such studies for the northern hemisphere is since most instruments at high latitudes are concentrated in the auroral zone and the polar cap since this is where the main interest is directed. We demonstrate the importance of taking into account the subauroral zone, which is where fairly stable plasma formations occur, with a strong flow and various-scale irregularities inside, which lead to serious interference in satellite communications.

This section describes instrumentations and methods used to process the observation data. To conduct this study, we had at our disposal an unprecedented set of data from different instruments in this region. Two scintillation receivers at SANAE IV (SAN, geographic coordinates ɸ = 71.7ºS; ℓ = 2.8ºW) and Troll (TRL, ɸ = 72.0ºS; ℓ = 2.5ºE) stations, which are located on the equatorial side and boundary of the quiet auroral oval. NorSat-1 satellite with a very high-resolution sampling data, flying in the evening sector in the field of view of the receivers, which allows us to evaluate the small-scale irregularities. DMSP satellites show the occurrence of particle precipitation and measure drift. Swarm satellites show the electron density variations. All data from the LEO satellites allow us to trace the dynamics of the position and development of PJ/SAID over time (represented in Fig. 4 by thick sections of the flyby trajectories of various LEO satellites). Magnetometers allow us to check the position and direction of PJ/SAID above the stations.

In our study, we used data from three space missions in low Earth orbits: Norsat-1, Swarm and DMSP. The Norwegian microsatellite NorSat-1 (∼16 kg, 23×39×44 cm) was launched on July 14, 2017 from the cosmodrome of Baikonur into a circular sun-synchronous orbit with an inclination of 98, altitude of ∼600 km, and orbital period of 95 min. The satellite has a multi-needle Langmuir Probe (m-NLP) developed at the University of Oslo. It should be noted that the classic Langmuir probe makes measurements over a range of voltages, building a current-voltage characteristic, from which the plasma parameters are calculated. Since the coverage of the entire range takes time, the time resolution of the obtained plasma parameters is usually low, which makes it difficult to study small-scale plasma structures. The m-NLP consists of four cylindrical probes (needles) at different voltage biases within the electron saturation region of the current-voltage characteristic [Hoang et al., 2018]. The absence of the need to cover the voltage range makes it possible to measure the electron current at a much higher frequency, which allows obtaining high-frequency characteristics of the plasma. To determine the electron density Ne in measurements by the system of LPs, at least two cylindrical probes operating at different fixed voltages are used [Jacobsen et al., 2010]. Also, a main feature of the m-NLP approach is that it does not need information about plasma potential and electron temperature to measure electron density. While the absolute electron temperature cannot be derived without knowing the satellite potential in accordance with the Langmuir probe theory, we can estimate the change of temperature with time and latitude [Chernyshov et al., 2020]. The high resolution of m-NLP instrument made it possible to study for the first time the small-scale structure of PJ/SAID in the subauroral region [Sinevich et al., 2021, 2022], as well as to discover a new phenomenon called Stratified Subauroral Ion Drift (SSAID) [Sinevich et al., 2023].

On 22 November 2013, the European Space Agency launched three Swarm satellites into nearly polar orbits. These satellites eventually reached heights of 460 km (Alpha and Charlie) and 510 km (Bravo). The main purpose of the Swarm mission is to carry out accurate multipoint measurements of low-frequency magnetic and electric fields in the ionosphere of the Earth. Two spherical LPs that were placed on the satellites were used to measure the values of electron density and temperature [Knudsen et al. 2017]. Swarm satellite data has also been used to study PJ/SAID (often in conjunction with other satellite data), as well as such new phenomena as STEVE (strong thermal emission velocity enhancement), the physical mechanism of which is closely related to PJ/SAID [e.g., Archer et al., 2019; MacDonald et al., 2018; Sinevich et al., 2022; Nishimura et al., 2019, 2023; Mishin and Streltsov, 2023].

A large number of studies conducted by various scientific teams in different years on the investigation of PJ/SAID are associated with the use of data from the Defense Meteorological Satellite Program (DMSP) satellite series. The DMSP satellites have been collecting weather data for U.S. military operations for over five decades and provide assured, secure global weather and space weather data. The DMSP satellites are polar-orbiting spacecraft at about 835–860 km altitude with 101 min orbital period. Many studies of PJ/SAID are based on data from the DMSP satellites [e.g., Anderson et al., 2001; Khalipov et al., 2016; Puhl-Quinn et al., 2007; Mishin, Puhl-Quinn, 2007; Mishin et al., 2010; Mishin, 2013; He et al., 2014; Horvath and Lovell, 2023b among others]. Also, using the data from these satellites, interesting cases were found when two peaks of PJ/SAID velocities were observed side by side on one satellite pass, that is, the so-called Double peak SAID (DSAID) [He et al., 2016; Wei et al., 2019; Horvath and Lovell, 2017]. Since NorSat-1 and Swarm are not equipped with ion drift measurements, data from the DMSP mission is used to indicate the existence of PJ/SAID explicitly in our work.

The Special Sensor Ultraviolet Spectrographic (SSUSI) instrument [Paxton et al., 1992] measures ultraviolet (UV) emissions in five different wavelengths and these instruments are flown on board the DMSP satellites. SSUSI is mounted on the nadir-looking panel of the satellite. The multicolour images from SSUSI cover the visible Earth disk from horizon to horizon and the anti-sunward limb up to an altitude of approximately 520 km. The SSUSI on board the DMSP satellites provide nearly hourly, 3000 km wide high-resolution (10km×10km) UV snapshots of auroral emissions. These UV data have been converted to average energies and energy fluxes of precipitating electrons. Thus, using SSUSI it is possible to determine the boundaries of auroral precipitation and thereby clearly understand where the boundary between the auroral and subauroral zones is in the ionosphere. As is known from various works (for example by Galperin [2002]), PJ/SAID is located inside the main ionospheric trough in the subauroral region and equatorward of the soft electron precipitation boundary, therefore, it is important to know these boundaries for a more accurate determination of the location of PJ/SAID.

The auroral energy flux derived from SSUSI images taken onboard DMSP is organised in Altitude Adjusted Corrected Geomagnetic (AACGM) coordinates for the nightside (18:00–06:00 MLT) of the southern hemisphere [Luan et al., 2018]. Therefore, the AACGM coordinate system is used everywhere in this study for uniformity. The AACGM coordinates [Shepherd, 2014] are an extension of corrected geomagnetic (CGM) coordinates that more accurately represent the actual magnetic field. In the AACGM coordinates, points along a given magnetic field line are given the same coordinates and thus better reflect the magnetic conjugacy.

Figure 5a, b shows position of the PJ/SAID at ~ 20:56:30 UT at geomagnetic latitudes from − 58.1° to -57.5° and at 18 MLT according to DMSP F-18 data. Despite the gap in the data, which is probably associated with instrument overscale, inside the electron density dip, which lies within the geomagnetic latitudes between − 58.5° and − 57.4°, an increase in the horizontal westward drift velocity of ions up to 2400 m/s is noticeable. According to Fig. 5b, the increase in the ion temperature by approximately twice the average level does not spatially coincide with the peak of the vertical and horizontal drift velocities and is located more poleward by approximately 2° in latitude. This may be due to a significantly lower spatial frequency of ion and electron temperature measurements (~ 0.25 Hz) than the spatial frequency of drift velocity measurements (1 Hz) on the DMSP satellite. The horizontal drift velocity of ions above 1000 m/s, coinciding with the electron density drop in the subauroral region of the ionosphere in the evening and midnight MLT sectors, indicates explicitly the presence of PJ/SAID in the considered region of the geomagnetic latitudes [Anderson et al., 1991]. As already mentioned, the position of this PJ/SAID is marked in Figs. 1–2, 4 with thick sections of the fly-by lines.

Figure 5c illustrates PJ/SAID detected from the high-frequency electron density and electron temperature data measured by the m-NLP onboard NorSat-1 satellite. As can be seen from Fig. 5c in an interval of the geomagnetic latitudes from − 53.5° to -54.2°, there is an electron density dip that spatially coincides with an increase in electron temperature value and fluctuations. Since NorSat-1 is not equipped with an ion drift velocity instrument, we cannot directly determine that the PJ/SAID is present in the considered time interval, however, the well-known properties of PJ/SAID (increase in temperature and decrease in ion density) indicate the presence of a PJ/SAID. In addition, thanks to DMSP F-18 horizontal ion drift velocity data it is already known that PJ/SAID existed at 20:56 UT at 18 MLT. According to the fact that average PJ/SAID spans over ~ 6 MLT [Karlsson et al., 1998] we conclude that PJ/SAID also existed at 21:20 UT at 00 MLT where NorSat-1 took measurements. Location inside the electron density decrease that coincides with an elevation in electron temperature at the subauroral latitudes in the evening and midnight MLT sector is characteristic of PJ/SAID [Anderson et al., 1991]. Besides that, the increase in fluctuations of the electron density and electron temperature are also characteristic of PJ/SAID [Sinevich et al., 2021, 2022; Mishin and Blaunstein, 2008]. Thus, we conclude that the considered electron density dip coinciding with an electron temperature increase detected from NorSat-1 data is PJ/SAID. Similarly, we conclude that PJ/SAID is present during the Swarm B and Swarm C passages since the Swarm Spherical LPs electron density and electron temperature data at 19:50, 20:48 and 22:20 UT show similar behaviour to the NorSat-1 data (Swarm satellites data are not shown in the figures).

Power spectrograms of the discrete Fourier transform of electron density and electron temperature are presented in Fig. 5d and Fig. 5e respectively. The spectrogram in Fig. 5d shows that fluctuations of electron density are present inside the PJ/SAID (an increase in spectral power at frequencies up to ~ 50 Hz which corresponds to a spatial size of plasma irregularities of several hundreds of meters) yet do not exceed the background level. On the contrary, the spectrogram in Fig. 5e illustrates a significant rise in fluctuations of electron temperature inside the PJ/SAID relative to the background level. The spectral power of electron temperature rises at frequencies up to 300–400 Hz which corresponds to the spatial size of plasma irregularities of several dozens of meters. Thus, we can conclude that PJ/SAID includes fluctuations and plasma irregularities of various scales, which can significantly affect radio communication, radio wave propagation, and GNSS positioning accuracy. This emphasises the relevance and importance of this study.

Together with the high-speed flow (that is usually above 1 km/s for PJ/SAID), such irregularities can lead to rapid fluctuations in the phase and amplitude measurements observed by receivers on the ground. To check how strong the influence of plasma structuring in PJ/SAID can be, we used data from two GNSS Ionospheric Scintillation and TEC Monitor (GISTM) receivers. The station SANAE IV has a NovAtel GSV4004B receiver, operated by the South African National Space Agency (SANSA). In December 2017, a NovAtel GPStation-6 receiver was installed at the Norwegian research station Troll as a part of Troll Ionospheric Observatory managed by the University of Oslo. GISTM receiver records signals from the GPS, GLONASS, and Galileo satellites, and it outputs raw observational data every second, including 50 Hz phase and amplitude measurements and 1 Hz TEC measurements. The raw phase measurements detrended with a 6th-order Butterworth high-pass filter with a 0.1 Hz cutoff frequency. NovAtel GSV4004B receiver records signals only from GPS satellites. We consider two scintillation indices: S4 defined as the normalized variance of the signal intensity [Yeh & Liu, 1982] and is an indicator for small-scale irregularities at and below the Fresnel radius, and σ_ϕ - the standard deviation of the measured phase (e.g., more in Rino, 1979) which is a proxy of plasma structuring on large and small scales. We focus also on the lock time parameter that indicates how long the receiver follows the satellite and has been locked to the carrier phase on the signal. The thin ionospheric layer approach and the AACGM system were used to recalculate the coordinates of ionospheric piercing points (IPPs) at 350 km [Li et al., 2018]. We chose an altitude of 350 km since PJ/SAID is most pronounced in the F region of the ionosphere. For both receivers, we used 1 Hz data (for the scintillation indices) from the L1 signal type (to comply with data resolution from SANAE IV) and thresholds of 10° elevation angle to collect more data that can cover regions located equatorward of the auroral oval.

The cut-off elevation mask to a GNSS satellite (elevation angle) is an essential parameter for studying the ionosphere. The receiver starts to track the satellite radio signal after the satellite rises above the horizon at a certain angle. Usually, the cut-off elevation mask value is set between 10° and 30° for ionospheric research. The advantages of using lower elevation angles are the following: 1) Lower cut-off angle would provide better spatial data coverage which is important because in this work we are considering the Antarctic region, where GNSS satellites are visible only from a certain angle due to the low inclination of the orbits of GNSS satellites. 2) Even at elevation angles larger than 30°, the multipath effects (due to signal propagation along several paths due to scattering and reflections in a medium with inhomogeneities) can result in phase slips and pseudorange distortions due to the signal interference [Chernyshov et al 2020b]. Therefore, in this study, we choose the cut-off angle of 10°. The Antarctic region is a less noisy area (lack of buildings, trees, etc), which makes it possible to consider lower cut-off elevation angles.

In addition, data from magnetometers at SANAE IV and Neumayer (ɸ = 70.7ºS; ℓ = 8.3ºW) stations were used to infer the ionospheric current system (Fig. 6). Neumayer data is available on the SuperMAG database. The baselines were subtracted from the measurements of the magnetic fields.

Author Contribution

D.K., A.S., A.Ch., D.Ch. and W.M. conceived the idea of the study. D.K. and A.S. conducted the data analysis and prepared figures. A.Ch. and D.K. wrote the main part of the manuscript. All authors discussed, analysed and interpreted the results, and worked on the manuscript.

Acknowledgement

DK, YJ, and WM acknowledge funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (ERC Consolidator Grant agreement No. 866357, POLAR-4DSpace).

Data Availability

The Swarm data is available at ftp://swarm-diss.eo.esa.int. Data from NorSat-1 is available at http://tid.uio.no/plasma/norsat. Data from the DMSP mission is available at the NOAA website https://satdat.ngdc.noaa.gov/dmsp/data. The SSUSI measurements data is available at https://ssusi.jhuapl.edu/data_products. Magnetometer data from the Neumayer station is available at the SuperMAG collaboration website https://supermag.jhuapl.edu/indices. Data from the SANAE IV research base is made available through the South African National Space Agency (SANSA) http://www.sansa.org.za. The Troll station datasets analysed during the current study are available in the zenodo repository, https://doi.org/10.5281/zenodo.12759371.

Archer, W. et al. The optical signature of intense subauroral ion drifts. Geophys. Res. Lett. 46, 6279–6286. https://doi.org/10.1029/2019GL082687 (2019).
Anderson, P. C., Heelis, R. A. & Hanson, W. B. The ionospheric signatures of rapid subauroral ion drifts. J. Geophys. Res. Space Phys. 96 (A4), 5785–5792. https://doi.org/10.1029/90JA02651 (1991).
Anderson, P. C. et al. A proposed production model of rapid subauroral ion drifts and their relationship to substorm evolution. J. Geophys. Res. Space Phys. 98 (A4), 6069–6078. https://doi.org/10.1029/92JA01975 (1993).
Anderson, P. C., Carpenter, D. L., Tsuruda, K., Mukai, T. & Rich, F. J. Multisatellite observations of rapid subauroral ion drifts (SAID). J. Geophys. Res. Space Phys. 106 (A12), 29585–29599. https://doi.org/10.1029/2001JA000128 (2001).
Figueiredo, S., Karlsson, T. & Marklund, G. Investigation of subauroral ion drifts and related field-aligned currents and ionospheric Pedersen conductivity distribution. Ann. Geophys. 22 (3), 923–934. https://doi.org/10.5194/angeo-22-923-2004 (2004).
Foster, J. C. & Vo, H. B. Average characteristics and activity dependence of the subauroral polarization stream. J. Geophys. Res. Space Phys. 107 (A12), 1475–1481. https://doi.org/10.1029/2002JA009409 (2002).
Galperin, Y. I., Ponomarev, V. N. & Zosimova, A. G. Plasma convection in the polar ionosphere. Ann. Geophys. 30 (1), 1–7 (1974).
Galperin, Y. I. Polarization jet: Characteristics and a model. Ann. Geophys. 20 (3), 391–404. https://doi.org/10.5194/angeo-20-391-2002 (2002).
Chernyshov, A. A. et al. In situ observations of ionospheric heating effects: First results from a joint SURA and NorSat-1 experiment. Geophys. Res. Lett. 47(13), e2020GL088462. (2020a). https://doi.org/10.1029/2020GL088462
Chernyshov, A. A., Miloch, W. J., Jin, Y. & Zakharov, V. I. Relationship between TEC jumps and auroral substorm in the high-latitude ionosphere. Sci. Rep. 10 https://doi.org/10.1038/s41598-020-63422-9 (2020b).
He, F., Zhang, X. X. & Chen, B. Solar cycle, seasonal, and diurnal variations of subauroral ion drifts: Statistical results. J. Geophys. Res. Space Phys. 119, 5076–5086. https://doi.org/10.1002/2014JA019807 (2014).
He, F., Zhang, X. X., Wang, W. & Chen, B. Double-peak subauroral ion drifts (DSAIDs). Geophys. Res. Lett. 43 (11), 5554–5562. https://doi.org/10.1002/2016GL069133 (2016).
Hoang, H. et al. The multi-needle Langmuir probe system on board NorSat-1. Space Sci. Rev. 214 (4), 75. https://doi.org/10.1007/s11214-018-0509-2 (2018).
Horvath, I. & Lovell, B. C. Investigating the development of double-peak subauroral ion drift (DSAID). J. Geophys. Res. Space Phys. 122 (4), 4526–4542. https://doi.org/10.1002/2016JA023506 (2017).
Horvath, I., Lovell, B. C. & JA031677. Antisunward streaming westward sub-auroral ion drifts (SAID) developed in the postmidnight (1–4) magnetic local time sector during 2013. J. Geophys. Res. Space Phys. 128, e (2023a). (2023). https://doi.org/10.1029/2023JA031677
Horvath, I. & Lovell, B. C. Abnormal Sub-Auroral Ion Drifts (ASAID) developed in the postmidnight (1–4) magnetic local time sector in 2013. J. Geophys. Res. Space Phys. 128 https://doi.org/10.1029/2023JA031808 (2023b). e2023JA031808.
Jacobsen, K. S., Pedersen, A., Moen, J. I. & Bekkeng, T. A. A new Langmuir probe concept for rapid sampling of space plasma electron density. Meas. Sci. Technol. 21 (8), 085902. https://doi.org/10.1088/0957-0233/21/8/085902 (2010).
Jin, Y. & Oksavik, K. GPS scintillations and losses of signal lock at high latitudes during the 2015 St. Patrick’s Day storm. J. Geophys. Res. Space Phys. 123 (9), 7943–7957. https://doi.org/10.1029/2018JA025933 (2018).
Karlsson, T., Marklund, G. T., Blomberg, L. G. & Mälkki, A. Subauroral electric fields observed by the Freja satellite: A statistical study. J. Geophys. Res. Space Phys. 103 (A3), 4327–4341. https://doi.org/10.1029/97JA00333 (1998).
Kintner, P. & Ledvina, B. The ionosphere, radio navigation,and global navigation satellite systems. Adv. Space Res. 35, 788. https://doi.org/10.1016/j.asr.2004.12.076 (2005).
Kintner, P. M., Ledvina, B. M. & de Paula, E. R. GPS and ionospheric scintillations. Space Weather. 5 (9), S09003. 10.1029/2006SW000260 (2007).
Khalipov, V. L., Galperin, Y. I., Stepanov, A. E. & Shestakova, L. V. Formation of a polarization jet at the outburst stage of a substorm: Results of ground-based measurements. Cosm. Res. 39 (3), 244–253. https://doi.org/10.1023/a:1017573319665 (2001).
Khalipov, V. L., Stepanov, A. E., Kotova, G. A. & Bondar’, E. D. Formation of polarization jet during injection of ions into the inner magnetosphere. Geomag. Aeron. 56 (2), 187–193. https://doi.org/10.1134/s0016793216020080 (2016).
Knudsen, D. J. et al. Thermal ion imagers and Langmuir probes in the Swarm electric field instruments. J. Geophys. Res. Space Phys. 122 (2), 2655–2673. https://doi.org/10.1002/2016JA022571 (2017).
Li, M. et al. Determination of the optimized single-layer ionospheric height for electron content measurements over China. J. Geod. 92, 169–183. https://doi.org/10.1007/s00190-017-1054-6 (2018).
Luan, X., Zhou, S. & Dou, X. Auroral energy flux distribution over the nightside auroral oval observed by the DMSP F16/SSUSI: Seasonal, geomagnetic, and solar activity dependences. J. Geophys. Res. Space Phys. 123, 4457–4466. https://doi.org/10.1029/2017JA023970 (2018).
MacDonald, E. et al. New science in plain sight: Citizen scientists lead to the discovery of optical structure in the upper atmosphere. Sci. Adv. 4 (3), eaaq0030. https://doi.org/10.1126/sciadv.aaq0030 (2018).
Mishin, E. V., Puhl-Quinn, P. A. & Said Plasmaspheric short circuit of substorm injections. Geophys. Res. Lett. 34 (24), L24101. https://doi.org/10.1029/2007GL031925 (2007).
Mishin, E. V. & Blaunstein, N. Irregularities within subauroral polarization stream-related troughs and GPS radio interference at midlatitudes. Am. Geophys. Union Geophys. Monogr. Ser. 181, 291–295. https://doi.org/10.1029/181GM26 (2008).
Mishin, E. V. Interaction of substorm injections with the subauroral geospace: 1. Multispacecraft observations of SAID. J. Geophys. Res. Space Phys. 118 (9), 5782–5796. https://doi.org/10.1002/jgra.50548 (2013).
Mishin, E. V. The evolving paradigm of the subauroral geospace. Front. Astronomy Space Sci. 10, 1118758. https://doi.org/10.3389/fspas.2023.1118758 (2023).
Mishin, E. V. & Streltsov, A. V. The inner structure of STEVE-linked SAID. Geophys. Res. Lett. 50 https://doi.org/10.1029/2023GL102956 (2023). e2023GL102956.
Nishimura, Y. et al. Magnetospheric signatures of STEVE: Implications for the magnetospheric energy source and interhemispheric conjugacy. Geophys. Res. Lett. 46, 5637–5644. https://doi.org/10.1029/2019GL082460 (2019).
Nishimura, Y., Dyer, A., Kangas, L., Donovan, E. & Angelopoulos, V. Unsolved problems in strong thermal emission velocity Enhance-ment (STEVE) and the picket fence. Front. Astronomy Space Sci. 10, 1087974. https://doi.org/10.3389/fspas.2023.1087974 (2023).
Oksavik, K., van der Meeren, C., Lorentzen, D. A., Baddeley, L. J. & Moen, J. Scintillation and loss of signal lock from poleward moving auroral forms in the cusp ionosphere. J. Geophys. Res. Space Phys. 120 (10), 9161–9175. https://doi.org/10.1002/2015JA021528 (2015).
Paxton, L. J. et al. Special sensor ultraviolet spectrographic imager: an instrument description, Proc. SPIE 1745, Instrumentation for Planetary and Terrestrial Atmospheric Remote Sensing, (1992). https://doi.org/10.1117/12.60595
Puhl-Quinn, P. A. et al. Cluster and DMSP observations of SAID electric fields. J. Geophys. Res. Space Phys. 112 (A5), A05219. https://doi.org/10.1029/2006JA012065 (2007).
Rino, C. L. A power law phase screen model for ionospheric scintillation: 1. Weak scatter. Radio Sci. 14 (6), 1135–1145. https://doi.org/10.1029/RS014i006p01135 (1979).
Shepherd, S. G. Altitude-adjusted corrected geomagnetic coordinates: Definition and functional approximations. J. Geophys. Res. Space Phys. 119, 7501–7521. https://doi.org/10.1002/2014JA020264 (2014).
Sinevich, A. A., Chernyshov, A. A., Chugunin, D. V., Miloch, W. J. & Mogilevsky, M. M. Studying the small-scale structure of a polarization jet during the April 20, 2018 geomagnetic storm. Sol Terr. Phys. 7 (1), 17–26. https://doi.org/10.12737/stp-71202103 (2021).
Sinevich, A. A. et al. Small-scale irregularities within polarization jet/SAID during geomagnetic activity. Geophys. Res. Lett. 49, e2021GL097107. (2022). https://doi.org/10.1029/2021GL097107
Sinevich, A. A. et al. Stratified Subauroral Ion Drift (SSAID). J. Geophys. Res. Space Phys. 128, e2022JA031109. (2023). https://doi.org/10.1029/2022JA031109
Skjæveland, A. H., Kotova, D. S. & Miloch, W. J. Case studies of ionospheric plasma irregularities over Queen Maud Land, Antarctica. J. Geophys. Res. Space Phys. 126 https://doi.org/10.1029/2021JA029963 (2021). e2021JA029963.
Southwood, D. J. & Wolf, R. A. An assessment of the role of precipitation in magnetospheric convection. J. Geophys. Res. Space Phys. 83, 5227–5232, https://doi.org/10.1029/JA083iA11p05227
Spiro, R. W., Heelis, R. A. & Hanson, W. B. Rapid subauroral ion drifts observed by Atmosphere Explorer C. Geophys. Res. Lett. 6(8), 657–660. (1978). (1979). https://doi.org/10.1029/GL006i008p00657
Spogli, L. et al. Climatology of GPS ionospheric scintillations over high and mid-latitude European regions. Ann. Geophys. 27, 3429–3437. https://doi.org/10.5194/angeo-27-3429-2009 (2009).
Wang, Y. et al. Experimental evidence on the dependence of the standard GPS phase scintillation index on the ionospheric plasma drift around noon sector of the polar ionosphere. J. Geophys. Res. Space Phys. 123, 2370–2378. https://doi.org/10.1002/2017JA024805 (2018).
Wei, D., Yu, Y., Ridley, A. J., Cao, J. & Dunlop, M. W. Multi-point observations and modeling of subauroral polarization streams (SAPS) and double-peak subauroral ion drifts (DSAIDs): A case study. Adv. Space Res. 63 (11), 3522–3535. https://doi.org/10.1016/j.asr.2019.02.004 (2019).
Zossi, B., Fagre, M., Amit, H. & Elias, A. G. Geomagnetic field model indicates shrinking northern auroral oval. J. Geophys. Res. Space Phys. 125 https://doi.org/10.1029/2019JA027434 (2020). e2019JA027434.
Yeh, K. C. & Liu, C. H. Radio wave scintillations in the ionosphere. Proceedings of the IEEE 70(4), 324–360. (1982). https://doi.org/10.1109/PROC.1982.12313

No competing interests reported.

Download PDF

Editorial decision: Revision requested
15 Oct, 2024
Reviews received at journal
13 Oct, 2024
Reviews received at journal
07 Oct, 2024
Reviewers agreed at journal
25 Sep, 2024
Reviewers agreed at journal
25 Sep, 2024
Reviewers invited by journal
23 Sep, 2024
Editor assigned by journal
23 Sep, 2024
Editor invited by journal
12 Sep, 2024
Submission checks completed at journal
09 Sep, 2024
First submitted to journal
30 Aug, 2024

You are reading this latest preprint version

Strong turbulent flow in the sub-auroral region in the Antarctic can deteriorate satellite-based navigation signals

Status:

Version 1

Abstract

Figures

Introduction

Main

Discussion

Conclusions

Methods

Declarations

Author Contribution

Acknowledgement

Data Availability

References

Additional Declarations

Supplementary Files

Status:

Version 1