GNSS
coordinates


Observation of travelling ionospheric disturbances over Morocco

Sep 2024 | No Comment

We investigate and present results on the propagation characteristics of AGWs associated with the Godzilla SDS event that occurred between 15th and 26th June 2020 over Morocco

Uluma Edward

Department of Physics, Masinde Muliro University of Science & Technology, Kakamega, Kenya

Chali Idosa Uga

Department of Space Science, University of Alabama in Huntsville, Alabama, USA Department of Physics, Jimma University, Jimma, Oromia, Ethiopia

Solomon Otoo Lomotey

Department of Smart Technologies, University of Environment and Sustainable Development, Somanya, Eastern Region, Ghana

Athwart Davis Odhiambo

Department of Physics, Masinde Muliro University of Science & Technology, Kakamega, Kenya

Fashae Joshua Bankole

Physics Programme, Bowen University, Iwo, Osun State, Nigeria

Kouassi Nguessan

Département de Physique, Université Nangui Abrogoua, Abidjan, Côte d’Ivoire

Muniafu Wilberforce

Department of Physics, Masinde Muliro University of Science & Technology, Kakamega, Kenya

Boniface Ndinya

Department of Physics, Masinde Muliro University of Science & Technology, Kakamega, Kenya

Omondi George

Department of Physics and Materials Science, Maseno University, Maseno, Kenya

Abstract

Using Vertical Total Electron Content (VTEC) data from four GNSS receiver stations: IFR1 (Ifrane Seismic), MELI (Melilla), TETN (Tetouan) and OUCA (Ouca) in Morocco, we investigate the Godzilla sand and dust storm (SDS) event of June 2020 as a source of energetics for generation of travelling ionospheric disturbances (TIDs). Godzilla SDS event began on 5th June 2020 in Algeria following a decrease in pressure and spread to other areas across the Sahara between 6th and 28th June 2020. It was tracked using the Sentinel- 5P Satellite mission. Solar wind parameters: Horizontal component of interplanetary magnetic field (IMF Bz), Y-component of interplanetary electric field (IE-Ey), and solar wind speed (Vz) and geomagnetic indices: disturbance storm time (Dst) and planetary K (kp) indices were examined and showed very minimal geomagnetic influence during the period. The study showed that TEC exhibited wave-like structures having distinct troughs and peaks over IFR1-MELITETN and OUCA-IFR1-TETN arrays which were clear indicators of generation of TIDs. The arrays and de-trended TEC plots showed that the TIDs propagated poleward. Neutral winds were seen to play a very important role in the propagation of the atmospheric gravity waves (AGWs) which are manifestations of TIDs. 1.

Introduction

Travelling ionospheric disturbances (TIDs) are known to be signatures of atmospheric gravity waves (AGWs) in the ionosphere [1,2]. These wave-like or TID perturbations in the F-layer depend on ionospheric parameters such as wind, electron density and ion temperature [3]. These fluctuations propagate as waves through the ionosphere with wide range of velocities and frequencies. TIDs play an important role in the exchange of energy and momentum between various regions of the upper atmosphere and ionosphere [4]. TIDs and AGWs are related in such a way that TIDs are manifestations of AGWs in the atmosphere [3]. This relationship has been highlighted by various notable studies by Hines, [1,2], Richmond, [5], Hunsucker, [6], Jing & Hunsucker, [7], Hocke & Schlegel, [3]. TIDs are triggered by the ions being forced along the field lines by the neutral air winds driven by the pressure wave [8]. They can also be triggered by geomagnetic or auroral activity, tropospheric activities associated with earthquakes, tsunami, volcanic eruptions, underground nuclear explosions, seismic events, and sand and dust storm (SDS) [4].

Based on their intrinsic properties, TIDs can be categorized either as Large Scale Traveling Ionospheric Disturbances (LSTIDs) or Medium Scale Traveling Ionospheric Disturbances (MSTIDs). LSTIDs have propagation period of 30 minutes to an hour while MSTIDs have propagation period of 15 to 60 minutes. LSTIDs have a wavelength greater than 1000 km and travel with velocities between 300-1000 m/s. MSTIDs have a wavelength of less than 1000 km and travel at a speed of between 100-300 m/s [1,2,5,6,7]. Geomagnetic or aurora activities generate LSTIDs while the lower atmosphere phenomena are a source of MSTIDs [9,10]. When the AGWs propagate vertically upwards from the lower atmosphere carrying energy, they become unstable and start losing energy. However, their amplitude increases due to reduced air density and mechanical friction [11]. MSTIDs have no correlation with geomagnetic activity. They are known to occur during both daytime and nighttime even though the mechanisms of generations causing the daytime and nighttime occurrences vary [12].

Daytime MSTIDs are known to originate from lower or middle atmosphere [13] while nighttime MSTIDs originate from electrodynamical processes which include the Perkins instability, that is associated with the E and F region coupling [14]. The climatology of MSTIDs has a strong dependence on longitude [4]. Tsagouri et al., [4] analyzed LSTIDs based on Digisonde observations over European region over Athens (380N, 23.50E), Dourbes (50.10N, 4.60N) and Ebre (40.80N, 0.50E). They observed a significant uplifting of the F2 layer which corresponded with an oscillation pattern in the foF2 during the LSTID activity. Variations in the height of the peak electron density hmF2 and the corresponding scale height Hm were also observed. However, it was difficult to make conclusion on any sensitivity of the method’s results to the electric-field-related disturbances which are common in mid-latitudes zone during geomagnetic storms. Habarulema et al, [14] investigated the storm time equatorward and poleward LSTIDs on global scale. They used total electron content (TEC) observations from GNSS receivers over South American, African and Asian sectors. The results showed the American and African sectors exhibiting predominantly strong poleward TIDs while the Asian sector exhibiting equatorward TIDs which crossed the geomagnetic equator on 9th March 2012. Katamzi & Habarulema, [15] also investigated TIDs over South African latitudes during the 29th to 31st October 2003 which was a geomagnetically disturbed period, using Global Positioning and System (GPS) observations. The results revealed presence of LSTIDs on the days of the geomagnetic storm using line-of-sight vertical TEC measurements from individual satellites. The wave-like structures were seen to propagate equatorward with velocities between 587.04 and 1635.09 m/s.

Sand and dust storm (SDS) events can be sources of energetics that can lead to generation of AGWs which would induce MSTIDs in the ionosphere. The MSTIDs give rise to plasma density gradients as they propagate through the ionosphere leading to electron density instabilities which manifests as TEC variation [16]. An interesting SDS event known as “Godzilla” SDS which occurred between 1st June and 30th June 2020 across the Sahara Desert was tracked by the Copernicus Sentinel-5P satellite which stores a variety of satellite images used to detect changes in landscape [17] and the offline Satellite Absorbing Index (SAI), which qualitatively shows elevated aerosol layers in atmosphere [18]. The National Oceanic and Atmospheric Administration (NOAA) approximated the Godzilla SDS to be between 60% to 70% larger than the Sand and dust storms experienced in the past [19]. The greatest impact of the Godzilla SDS was felt between 15th and 20th June 2020 [20,21,22].

In the present paper, we investigate and present results on the propagation characteristics of AGWs associated with the Godzilla SDS event that occurred between 15th and 26th June 2020 over Morocco. The study is investigated within geomagnetically quiet period, hence geomagnetic storm-related phenomena is eliminated. We analyze GPS-TEC data from four Global Navigation Satellite System (GNSS) receiver stations situated in Morocco.

2. Materials and methods

2.1 Geomagnetic Indices Data and Solar Wind Parameters

In this study, the periods between 15th and 26th June 2020 were considered due to greatest intensity of the Godzilla SDS event over Morocco [20]. The geophysical conditions such as Kp index showing the level of the geomagnetic activity, the hourly values of Z-component of the interplanetary magnetic field (IMF-Bz), Y-component of the interplanetary electric field (IEFy), solar wind speed (Vz), Planetary K (Kp) and disturbance storm time (Dst) indices were considered during these periods. These parameters were obtained from the archive of OmniWeb website via https:// omniweb.gsfc.nasa.gov/form/dx1.

2.2 Godzilla SDS Data

The movement of the Godzilla SDS over Morocco was tracked using the Copernicus Sentinel-5P satellite mission which was launched by the European Space agency on 13th October 2017. This satellite uses the Google Earth Engine (GEE), which is a cloud-based platform that stores a variety of satellite images used to detect changes in landscape [23] and the Sentinel-5P Aerosol Index (SAI), which qualitatively shows elevated aerosol layers in atmosphere and the Sentinel-5P Aerosol Index (SAI), which qualitatively shows elevated aerosol layers in atmosphere [18].

2.3 GPS/GNSS Data

The GPS-TEC data were obtained from four GNSS receiver stations in Morocco as shown in the map displayed in Fig. 1 with the station’s geographical locations. Table 1 shows the station codes, the station names, geographical and geomagnetic coordinates, and their respective local time. The GPS-TEC data were accessed from the University of NAVSTAR Consortium dual frequency website: http://unavco.org/data/gps-gnssdata/.

The GPS-TEC data in the GPS receivers is usually saved in the zipped Receiver Independent Exchange (RINEX) Format which is converted to GPS observable files using appropriate software [24]. In this study, the Gopi software developed by Prof Gopi Krishna Seemala [25] was used to convert the RINEX to observation files. GPS measurements are either code pseudoranges (P) or carrier phases (ϕ). The GPS receiver receives the code time delay and carrier phase difference by cross-correlating the f 1 and f 2 modulated carrier signals, which are normally considered to travel along the same path through the ionosphere [26]. Estimates of GPS-derived ionospheric TEC can be obtained using dual frequency GPS measurements [27,28]. GPS receiver data is critical for estimating the electron density along a ray path between a GPS satellite and a ground receiver [29,30]. Dual-frequency GPS receivers may offer integral information on the ionosphere and plasmasphere. This is done by computing the differential of the code and carrier phase measurements, in addition to removing ionospheric inaccuracies in TEC estimates [31,32]. As a result, the GPSTEC computed by the dual-frequency receivers is offered as an input to an ionosphere assimilation model [33]. For the present study, GPS-TEC data collected in dual-frequency receivers was used. The GPS-TEC data was obtained using the pseudo-range and carrier phase measurements. The TEC calculated from the pseudo-range measurement (slant TEC) is given by Equation 1:

Where α is the satellite’s elevation angle, RE is the Earth’s mean radius, and h is the height of the ionospheric layer, which is considered to be 400 km.

To reduce multipath effects, the data selected was for elevation angles of 300and above [35]. This was to eliminate multipath errors.

Information for the four GNSS receiver stations: MELI, IFR1, OUCA and TETN is given in Table 1.

The average of VTEC for all PRNs for IFR1, MELI, OUCA and TETN during study period were obtained by averaging the VTEC values for all identical pseudorandom numbers (PRNs) within a 24 hour period [36] and used to plot VTEC against Universal Time (UT) for each day and station. The VTEC against universal time (UT) plots for each day were analyzed.

2.4 Methodology

2.4.1 Extracting travelling Ionospheric disturbances (TIDs) from GPS-TEC data

Background or unperturbed ionospheric TEC was determined by fitting each satellite’s time series of the VTEC data with a fourth polynomial across all the four GNSS receiver stations. This was done using equation 6.

Detrended TEC from GPS-TEC data was created using the Savitzky-Golay filter with a 60-minute sliding window for all the four GNSS receiver stations and visible satellites [38]. The Savitzky-Golay filter uses quadratic polynomial that is fitted over each sliding window. The detrended TEC over all geographic latitudes against Universal time (UT) were directly constructed in a temporal resolution of 10 minutes. The derived detrended TEC plots along geographical latitudes were then used to investigate propagation behavior of the Godzilla SDS induced TIDs [39].

3. Results and discussion

3.1 The Movement of the Godzilla SDS between 15th and 26th June 2020 over Sahara

The movement of the Godzilla SDS over Morocco between 15th and 26th June 2020 was tracked using the Copernicus Sentinel- 5P satellite. During the study period, the Western part of Africa such as Mauritania, Mali, Morocco and Western Sahara felt the greatest impact of the Godzilla SDS [2] as in 2 and 3. By 20th June 2020, the first traces of dust had reached the Caribbean [40]. Large SDS traces were observed over the upper parts of the Sahara between 21st June and 26th June 2020 as shown in Fig. 3.

In Figs. 2 and 3, the cream colour on the SAI indicates low aerosol concentrations in the atmosphere while red colour shows presence of high aerosol concentrations in the atmosphere during the Godzilla SDS of June 2020 as tracked by the Copernicus Sentinel-5P satellite.

3.2 The Variation of Solar Wind Parameters between 15th and 26th June 2020

Fig. 4 depicts the changes in IMF Bz, IEF Ey, the solar wind speed, kp and Dst indices between 15th and 26th June 2020.

Noted that from Fig. 4, on 15th June 2020, there was a southward turning of the IMF-Bz of -3 nT with a corresponding rise of IEF-Ey of 1. The solar wind speed rose to maximum, of 330 km/s. The highest kp index value of 1.5 and lowest Dst value of -8 nT were attained. On 16th June 2020, the southward turning of -2.5 nT was attained with a corresponding rise in IEF-Ey of 1.5. The lowest solar wind speed of 280 km/s was attained. The maximum kp index value attained was 2.2 with minimum Dst value of -11 nT. On 17th and 18th June 2020 there was a southward turning of -5 nT with a corresponding IEF-Ey of 1. The solar wind speed dropped from 320 km/s to 280 km/s on 17th June and rose to a maximum of 320 km/s on 18th June 2020. The maximum Dst value of 1.8 nT was attained on 17th and 18th June with a corresponding Dst value of -10 nT. On 20th June 2020, the IMF-Bz southern turning reached -6 nT with a corresponding IEF-Ey value of 2. This corresponded well with a solar wind speed of 350km/s. Between 15th and 26th June 2020, the highest Kp value of 2.2 was attained with a corresponding Dst index value of -11 nT on 19th June 2020. The The solar wind parameter shows a minimum solar wind speed of 270 km/s on 16th June 2020 and maximum solar wind of 360 km/s on 20th June 2020. During solar minimum, solar wind of speeds between 250 km/s and 400 km/s originates from regions close to the heliospheric current sheet at the heliomagnetic equator while the fast solar wind having speeds of between 400 km/s and 800 km/s which originates from coronal holes when the magnetic field is open [41]. The solar wind parameters and geomagnetic indices in Fig. 4 indicate a very minimal geomagnetic field influence on the ionosphere between 15th and 26th June 2020.

3.3 Variations of TEC and VTECfit against UT

Figs. 5 and 6 show the VTEC and VTECfit against UT for OUCA-IFR1-TETN receiver array while Figs. 7 and 8 shows the VTEC and VTECfit against UT for IFR1- MELI-TETN receiver array between 15th and 26th June 2020 respectively.

3.3.1 OUCA-IFR1-TETN receiver array

Based on Figs. 5 and 6, on 15th, 23rd and 24th June 2020 triple peak structures (indicated by black arrows) were observed between 7:00 and 21:00 UT with the first peak appearing approximately at 08:00 UT, the second peak appearing approximately at 16:00 UT, while the third peak appearing approximately at 20:00 UT. From 16th to 20th June 2020, double peak structures were observed between 07:00 and 12:00 UT with the first peak appearing approximately at 08:00 UT and the second peak appearing approximately at 20:00 UT. On 16th, 17th, 18th, 19th , 21st , 22nd, 25th and 26th June 2020, double peak structures however were observed between 7:00 and 12:00 UT with the first peak appearing approximately at 08:00 UT and the second peak appearing approximately at 20:00 UT.

3.3.2 IFR1-MELI-TETN receiver array

Based on Figs. 7 and 8, it was also noted that on 15th , 23rd and 24th June 2020, triple peak structures were observed between 7:00 and 21:00 UT with the first peak appearing approximately at 08:00 UT, the second peak appearing approximately at 16:00 UT while the third peak appearing approximately at 20:00 UT Double peaks were observed on 16th, 17th, 18th, 19th, 20th, 21st, 22nd, 25th and 26th June 2020 with the first peak appearing approximately at 10:00 UT , the third peak appearing approximately at 18:00 UT. A close analysis of Figs. 5, 6, 7 and 8 showed that the two receiver arrays exhibited similar wave-like structures. The peaks were observed to all appear between 7:00 and 12:00UT after sunrise and between 18:00 and 21:00UT after sunset.

It was also noted that the two receiver arrays show peak structures for the lower latitude stations (OUCA and IFR1) being observed first as compared to the other receiver stations. This indicates that there was a pole-ward propagation of the TIDs.

 

3.4 Changes in Detrended TEC (∆TEC)

Figs. 9 to 14 show plots of detrended TEC along geographic latitude against UT for the four GNSS receiver stations between 15th and 26th June 2020.

In these Figures, band structures were observed over the 24- hour period for all the days between 15th and 26th June 2020. In Figs. 9(a), 9(b), 11(b), 12(a) and 12(b) the band structures for DTEC ranged between 0.5 and 1.0 TECU. In Figs. 10(a), 10(b), 11(a), 13(a), 13(b) the band structures for DTEC ranged between 0.6 and 1.5 TECU. In Figs. 14(a) and 14(b), the band structures for DTEC ranged between 0.5 and 1.2 TECU.

In Figs. 9 – 14, we can observe evident poleward propagating band structures during the 24-hour period for all the days between 15th and 26th June 2020. The band structures depicts present of TIDs. We also observe important values of detrended DTEC at the time of VTEC peaks. For example, Figs. 10(a) and 11(a) show DTEC values of above 1 when VTEC has the large peaks around 18:00 UT. The band structures depict presence of TIDs that is important when VTEC is maximum.

The present study successfully demonstrates that SDS events are a source of energetics for generation of TIDs which are signatures of AGWs. The Dst and kp indices in Fig. 4 rule out the influence of geomagnetic activity in the generation of TIDs. Therefore, the TEC results from this study exclude the effect of geomagnetic disturbances but are fully attributed to the effect of the SDS. During SDS events, internal waves are continuously generated and breaking throughout the atmosphere [11]. Wave breaking effectively mixes sand and dust aerosols through the atmosphere and contributes to driving some larger-scale flows. The gravity waves (GWs) generated are affected by the changing stratification as they propagate. Each of these interactions, in addition to gravity wave dissipation, may contribute to the vertical flux of horizontal momentum and the universal frequency spectrum in the middle atmosphere [42,43,44,45]. The possible cause of these types of disturbances is the neutral winds [46]. Neutral winds play an important role in the wave breaking/dissipation in the mesosphere and lower thermosphere to global redistribution of energy and momentum deposited at high latitudes by the magnetosphere [47]. They are usually slower than the actual wind in stable conditions and faster in unstable conditions i.e. land surface or sea stress. When GWs propagate into the ionosphere, the measured TEC exhibits clear wavelike-structures, as shown in Figs. 5, 6, 7 and 8. Their peaks and troughs are imaged on the ionospheric pierce points (IPP) trajectory at different moments and different propagation distances. This makes observation of GWs propagation possible.

4. Conclusion

We have investigated the generation of TIDs during the Godzilla SDS event of June 2020 over Morocco using GPSTEC data. The obtained results showed TEC exhibited wavelike structures having distinct troughs and peaks over IFR1- MELI-TETN and OUCA-IFR1-TETN arrays which were clear indicators of generation of TIDs. The arrays showed that the TIDs propagated poleward (along latitude). This was supported by the evident poleward propagating band structures on the detrended TEC plots during the 24-hour period for all the days between 15th and 26th June 2020. The band structures on the plots also depicted the presence of TIDs. Neutral winds were seen to play a very important role in the propagation of the AGWs which are manifestations of TIDs. In conclusion, the study confirms that SDS events can be a source of energetics for the generation of MSTIDs.

Disclaimer (Artificial Intelligence)

Author(s) hereby declare that NO generative AI technologies such as Large Language Models (ChatGPT, COPILOT, etc) and text-to-image generators have been used during writing or editing of manuscripts.

Acknowledgements

The authors thank the University of NAVSTAR Consortium: http://unavco.org/data/gps-gnssdata/ for the GNSS data; Kyoto: wdc.kugi. kyoto.u.ac.jp/dst/index.html and Kyoto: www. kugi.kyoto-ua.ac.jp/kp for the geomagnetic activity data; https://omniweb.gsfc.nasa.gov/ form/dx1.html for the daily values of solar wind parameters and the European Space Agency for the Copernicus Sentinel-5P satellite images of June 2020 Godzilla SDS over the Sahara. They also thank Prof Gopi Seemala of the Indian Institute of Geomagnetism for the GPS-TEC analysis software.

Competing interests

Authors have declared that no competing interests exist.

References

1. Hines CO. Internal atmospheric gravity waves at ionospheric heights, Can. J Phys. 1960;38:1414.

2. Hines CO. The upper atmosphere in motion, Geophysical Monograph, American Geophysical Union; 1974.

3. Hocke K, Schlegel K. A review of atmospheric gravity waves and travelling ionospheric disturbances: 1982–1995, Annales Geophysicae.1996;14:917.

4. Tsagouri I, Belehaki A, Koutroumbas K, Herekakis T. Identification of large scale ionospheric travelling disturbances (Lstids) based on digisonde observations. Atmosphere. 2023;14:331. Available:https://doi.org/10.3390/ atmos140 20331

5. Richmond A. Gravity wave generation, propagation, and dissipation in the thermosphere, J. Geophys. Res. 1978;83:4131.

6. Hunsucker RD. Atmospheric gravity waves generated in the high latitude ionosphere: A review, Rev. Geophys. Space Phys. 1982;20:293.

7. Jing N, Hunsucker RD. A theoretical investigation of sources of large and medium scale atmospheric gravity waves in the auroral oval, J. Atmos. Terr. Phys. 1993;55:1667.

8. Balthazor RL, Moffett RJ. A study of atmospheric gravity waves and travelling ionospheric disturbances at equatorial latitudes. Annales Geophysicae. 1997;15:1048-1056.

9. Ding F, Wan W, Ning B, Zhao B, Li Q, Zhang R, Xiong B, Song B. Two dimensional imaging of large scale travelling ionospheric disturbances over China based on GPS data. J. Geophys. Res. 2012;117:A08318.

10. Ferreira AA, Borries C, Xiong C, Borges RA, Mielich J, Kouba D. Identification of potential precursors for the occurrence of large scale travelling Ionospheric disturbances in a case study during September 2017. J. Space weather. Space Clim. 2020;2022:10, 32.

11. Fritts DC, Yuan L. Stability analysis of inertia-wave structure in the middle atmosphere. Journal of Atmospheric Sciences. 1989;46(12):1738-1745. Available:https://doi.org/10. 1175/1520- 0469(1989)046<1738 : saoiws>2.0.co;2

12. Kotake N, Otsuka Y, Tsugawa T, Ogawa T, Saito A. Climatological study of GPS total electron content variations caused by medium-scale travelling ionospheric disturbances. J. Geophys. Res. 2006;111:1-8.

13. Frissell NA, Baker JBH, Ruohoniemi JM, Greenwald RA, Gerrard AJ, Miller ES, West ML. Sources and characteristics of mediumscale travelling ionospheric disturbances observed by high frequency radars in the North America Sector. J. Geophys. Res. Space Physics. 2016;121: 3722-3739.

14. Habarulema JB, Katamzi ZT, Yingezaw E, Yamazaki Y, Seemala G. Simultaneous storm time equatorward and poleward large scale TIDs on a global scale. Geophysical Research Letters. 2016; 43(13):6678-6686.

15. Katamzi ZT, Habarulema JB. Travelling ionospheric disturbances observed at south african mid-latitudes during the 29th – 31st October 2003 geomagnetically disturbed period. Advances in space Research. 2003;53(2014):48-62.

16. McNamara L. The ionosphere: Communications, Surveillance and direction finding Orbit, a foundation series. Krieger Pub. Co; 1991.

17. Gorelick N, Hancher M, Dixon M, Ilyushchenko S, Thau D, Moore R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sensing and Environment. 2017;202:18-27.

18. De Graaf M, Stammes P, Torres O, Koelemeijer RBA. Absorbing aerosol index: Sensitivity analysis, application to GOME, and comparison to TOMS. Journal of Geophysical Research: Atmospheres. 2005;110(D1).

19. Zhongming Z, Linong L, Xiaona Y, Wangqiang Z, Wei L. Satellites and Unusual Sahara dust plume; 2020.

20. Asare-Ansah AB, Twumasi YA, Ning ZH, Ansah PB, Frimpong DB, Owusu F, Apraku CY, Anokye M, Loh M, Armah RDN, Opong J. Tracking The Godzilla Dust plume using the Google Earth Engine Platform. The international Archives of the photogrammetry, Remote sensing and spatial information sciences, Vol. XLVI-M-2-2022, ASPRS 2022 Annual conference, 6-8 February and 21-25 March, Denver, Colorado, USA and Virtual; 2022. Available:https://doi.org/10.5194/ isprs-Archives-XLVI-M-2-2022-33-2022

21. Althaf P, Shaeb KHB, Kumar KR. Hotspot analysis and long-term trends of Absorbing aerosol index from dust emissions were measured by the Ozone measuring Instrument at different Urban locations in India from 2005 to 2008. Atmospheric Environment. 2022;118933.

22. Jethva H, Satheesh SK, Srinivasan J.Seasonal variability of aerosols over theIndoGangetic basin. Journal of Geophysical Research: Atmospheres. 2005;110(D21).

23. Gorelick N, Hancher M, Dixon M,Ilyushchenko S, Thau D, Moore R. GoogleEarth Engine: Planetary-scale geospatialanalysis for everyone. Remote sensingand Environment. 2017;202:18-27.

24. Gurtner W, Estey L. Rinex-the receiverindependent exchange formatversion3.00. Astronomical Institute, University ofBern and UNAVCO, Bolulder, Colorado;2007.

25. Seemala GK. GPS-TEC AnalysisApplication Read Me. Institute for ScientificResearch, Boston College, USA; 2011.

26. Tariku YA. Patterns of GPS-TEC variationover low-latitude regions (African sector)during the deep solar minimum (2008 to2009) and solar maximum (2012 to 2013) phases. Earth, Planets and Space. 2015;67(1):1–9.

27. Horvath I, Crozier S. Software developed for obtaining GPS-derived total electron content values. Radio Science. 2007; 42(02):1–20.

28. Cepni MS, Potts LV, Miima JB. Highresolution station- based diurnal ionospheric total electron content (TEC) from dual-frequency GPS observation Space Weather. 2013;11(9): 520–528.

29. Ma X, Maruyama T, Ma G, et al. Threedimensional ionospheric tomography using observation data of gps ground receivers and ionosonde by neural network. Journal of Geophysical Research: Space Physics. 2005;110(A5).

30. Jin S, Park J, Wang J, Choi B, Park P. Electron density profiles derived from ground-based observations. The Journal of Navigation. 2006;59(3):395-401.

31. Heise S, Jakowski N, Wehrenpfennig A, et al. Sounding of the topside ionosphere/ plasmasphere based on GPS measurements from champ: Initial results. Geophysical Research Letters. 2002; 29(14):44–1.

32. Jawoski N, Mayer C, Hoque M, Wilken V. Total electron content models and their use in ionosphere monitoring. Radio Science. 2011;46(06):1-11.

33. Ciraolo L, Azpilicueta F, Brunini C, Meza A, Radicella SM. Calibration errors on experimental slant total electron content determined by the GPS. Journal of Geodesy. 2007;81:111-120.

34. Kassa T, Damtie B. Ionospheric irregularities over Bahir Dar, Ethiopia during selected geomagnetic storms, Advances in Space Research; 2017 Available:http://dx.doi. org/10.1016/j.asr.2017.03.036

35. Otsuka Y, Ogawa T, Saito A, Tsugawa T, Fukao S, Miyazaki S. A new technique for mapping of total electron content using GPS network in Japan. Earth, Planets and Space. 2002;54(1):63-70.

36. Sardon E, Zarraoa N. Estimation of total electron using GPS data: How stable are the differential satellite and receiver instrumental biases? Radio Science. 1997;32(5):1899-1910.

37. Aryal S, Geddes G, Finn S, et al. Multispectral and Multi-instrumentation Observation of TIDs following the solar eclipse of 21st August 2017. Journal of Geophysical Research: Space Physics. 2019;124.

38. Zhang SR, Coaster AJ, Erickson PJ, Goncharenko LP, Rideout W, Vierinen J. Travelling ionospheric disturbances and ionospheric perturbations associated with solar flares in September 2017. Journal of Geographical Research: Space Physics. 2019;124(7):5894-5917.

39. Zakharenkova I, Astafyeava E, Cherniak I. GPS and GLONAS observations of laregscale travelling ionopspheric disturbances during the 2015 St. Patrick’s day storm. Journal of Geophysical Research: Space Physics. 2016;121(12):12-138. Available:https://doi. org/10.1002/2016ja032332

40. Warren C. Godzilla dust storm traced to shaky northern jet stream; 2023. Available:https://www.science.org/ content/article/godzilla-dust-storm-Tracedshaky-northern-jet-stream, 2022

41. Kamide Y, Chian A. Hand Book of the SolarTerrestrial environment. Springer Berlin Heidelberg, New York. ISBN 978-3-540- 46314-6; 2007 DOI: 10.007/6104478

42. VanZandt TE. A universal spectrum of buoyancy waves in the atmosphere. Geophys. Res. Lett. 1982;9:575- 578.

43. Balsley BB, Carter DA. The spectrum of atmospheric velocity fluctuations at 8 and 86 km. Geophys. Res. Lett. 1982;9:465–468.

44. Nakamura T, Tsuda T, Fukao S, Kato S, Manson AH, Meek CE. Comparative observations of short-period gravity waves (10-100 min) in the mesosphere in 1989 by Saskatoon MF radar (52on), Canada and the MU radar (35on), Japan. Radio Sci. 1993;28:729–746.

45. Collins RL, Nomura A, Gardner CS. Gravity waves in the upper mesosphere over Antarctica: Lidar observations at the South Pole and Syowa. J. Geophys. Res. 1994;99:5475–5485.

46. Jonah OF, Zhang S, Coster AJ, Goncharenko LP, Erickson PJ, Rideout W, De Paula ER, De Jesus R. Understanding Inter-Hemispheric Traveling Ionospheric Disturbances and Their Mechanisms. Remote Sens. 2020;12:228.

47. Dhadly M, Sassi F, Emmert J, Drob D, Conde M, Wu Q, Makela J, Budzien S, Nicholas A. Neutral winds from mesosphere to thermosphere—past, present, and future outlook. Front. Astron. Space Sci. 9: 1050586. DOI: 10.3389/fspas.2022.1050586

The paper originally published in International Astronomy and Astrophysics Research Journal may be cited as: Edward, Uluma, Chali Idosa Uga, Solomon Otoo Lomotey, Athwart Davis Odhiambo, Fashae Joshua Bankole, Kouassi Nguessan, Muniafu Wilberforce, Boniface Ndinya, and Omondi George. 2024. “Observation of Travelling Ionospheric Disturbances over Morocco During the Godzilla Sand and Dust Storm of 15th to 26th June 2020 Using GNSS”. International Astronomy and Astrophysics Research Journal 6 (1):18-39. https://www.journaliaarj. com/index.php/IAARJ/article/view/100.

© Copyright (2024): Author(s). The licensee is the journal publisher. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The ariicle is republished with author’s permission

 

Leave your response!

Add your comment below, or trackback from your own site. You can also subscribe to these comments via RSS.

Be nice. Keep it clean. Stay on topic. No spam.