Microsoft Word - ETASR_V13_N3_pp10936-10940 Engineering, Technology & Applied Science Research Vol. 13, No. 3, 2023, 10936-10940 10936 www.etasr.com Nimmakayala & Vemuri: Analysis of Ionospheric Scintillations using GPS and NavIC Combined … Analysis of Ionospheric Scintillations using GPS and NavIC Combined Constellation Shiva Kumar Nimmakayala Department of EECE, GITAM School of Technology, GITAM Deemed to be University, India 121960402502@gitam.in (corresponding author) Bala Sai Srilatha Indira Dutt Vemuri Department of EECE, GITAM School of Technology, GITAM Deemed to be University, India svemuri@gitam.in Received: 17 March 2023 | Revised: 19 April 2023 | Accepted: 21 April 2023 Licensed under a CC-BY 4.0 license | Copyright (c) by the authors | DOI: https://doi.org/10.48084/etasr.5863 ABSTRACT The disturbances and irregularities in the ionosphere are the primarily recognized ramifications of space weather called scintillations. Irregularities in the electron densities are the source of the ionospheric scintillations. This article investigates the ionospheric scintillations, which are predominant in the trans- equatorial and equatorial regions. Based on the data from a multi-constellation Global Navigation Satellite Systems (GNSS) receiver at the Chaitanya Bharathi Institute of Technology Hyderabad, the relationship between the amplitude scintillation index S4 and the rate of change of total electron content (ROTI) is examined. The correlation coefficient between S4 and ROTI is demonstrated in this article. The outcome validates the usefulness of the ROTI in identifying the scintillations. Keywords-ionospheric scintillations; rate of change of total electron content; amplitude scintillation index; Global Navigation Satellite System (GNSS) I. INTRODUCTION The advancement of the Global Navigation Satellite System (GNSS), which encompasses GLONASS, GPS, BeiDO, Galileo, and numerous regional navigation systems, directed advanced steps in the earth’s atmosphere research [1]. The space-based investigation of the GNSS signals is one of best enforced ways to investigate the ionosphere [2]. The space weather trans-ionospheric signals can undergo expeditious variations and disturbances both in amplitude and phase when they pass through indiscriminate electron density irregularities which are present in the ionosphere [3-4]. This anomaly is generally referred to as ionospheric scintillations. The temporal and spatial status of the ionospheric scintillations depends on geographical and geophysical locations, the time of the day, and season [5-6]. The ionospheric scintillations are the major sources of degradation of the VHF band and C band signals, immensely effecting navigation and communication systems [7]. Hence, it is vital to investigate and analyze them. The GNSS signals are considered a robust tool to investigate the Total Electron Content (TEC) and the existence of ionospheric scintillations [8]. Phase and amplitude scintillations are both ionosphere irregularities. II. GNSS DATA COLLECTION For this analysis, the data of the multi constellation GNSS receiver of Chaitanya Bharati Institute of Technology, which is located in Hyderabad, during the period from 2020-06-07 to 2020-06-13 were taken into account. The considered parameters are time, PRN (Pseudo Random Number) of the satellite, carrier to noise values of the received signal, and the ionospheric delay of the received signal. To reduce complexity, only the results of 2020-06-07 are shown in this paper. III. METHODOLY The ionospheric irregularities and scintillations can be investigated in two ways when using GNSS signals. The first of these methods uses carrier to noise (� ��⁄ ) ratio to calculate the amplitude scintillation index S� [9] and the second one uses ionospheric delay to calculate the rate of change of TEC index (ROTI) which is the replacement index of S� [10]. The current work is based on the multi-step interpretation and processing of the L (1575.42MHz), L� (1176.45MHz), and S (2492.028MHz) GNSS signals corresponding to GPS and NavIC satellite constellations. In the first method, the amplitude scintillation index S� is calculated as the ratio of the standard deviation to the mean value of the averaged � ��⁄ [9]. The effectiveness of the amplitude scintillations index is divided into 3 categories which are week scintillations when S� is less than 0.3, moderate scintillations when S� ranges between 0.3 and 0.6, and strong scintillations when S� is above 0.6 [11]. The amplitude scintillations are calculated as follows: �� = 〈�� �〉�〈��〉� 〈��〉� (1) Engineering, Technology & Applied Science Research Vol. 13, No. 3, 2023, 10936-10940 10937 www.etasr.com Nimmakayala & Vemuri: Analysis of Ionospheric Scintillations using GPS and NavIC Combined … where SI is the intensity of the received trans-ionospheric signal. The S� value is computed using the signals received on the L (1575.42MHz), L� (1176.45MHz) and S (2492.028MHz) frequencies for a week-long data from 2020- 06-07 to 2020-06-13. The second method is to calculate the ROTI using the delay due to ionosphere in the received signals [12]. The ionospheric delay is computed by: �����= ��.� �� *TEC (2) where TEC indicates the Total Electron Content on the frequency of the received trans-ionospheric signal and ����� is the ionospheric delay of the received signal. TEC is used to describe all the electrons that are present along the path from the satellite to the receiver [12-13]. The TEC is expressed below and is measured in TECU: 1 TECU = 10 � e/m! (3) TEC = �"#$% ∗�� ��.� (4) The rate of change of TEC (ROT) is calculated at regular intervals by (5): ROT = '()*+�'()*,-+ .*�.*,- (5) where b specifies the epoch time and a specifies the noticeable satellite. ROTI is calculated by the standard deviation of ROT: ROTI = /〈ROT!〉 3 〈ROT〉! (6) The effectiveness of ROTI is divided into three categories: It is fragile when the ROTI index ranges between 0.25 and 0.5, modest when the ROTI index ranges between 0.5 and 1, and vigorous when the ROTI index is above 1 [14-15]. ROTI was calculated for signals received on L , L�, and S frequencies for a weak long data from 2020-06-07 to 2020-06-13. To observe the correlation between ROTI and S� , the covariance and standard deviation of S� and ROTI are computed [16] by: 41�5 = 〈S�!〉 3 〈S�〉! (7) 416�'� = /〈ROTI!〉 3 〈ROTI〉! (8) The covariance of S� and ROTI is: � 6�'��5 = ∑ 9:5�:5 ;;;;)9<=>?�<=>?;;;;;;;;) @ @� (9) where A�B is the mean of �� and ROTI;;;;;;; is the mean of ROTI. The correlation coefficient is a measure of the correlation between ROTI and S� methods which is expressed by C 6�'��5 , and it is given as: C 6�'� �5 = Ϲ-EFGHI5 J EFGHJ I5 (10) where Ϲ16�'��5 specifies the covariance of S� and ROTI . The flowchart of the followed methodology is shown in Figure 1. Fig. 1. Flowchart of the methodology. IV. RESULTS AND DISCUSSION The pseudo range measurements data of multi-constellation GNSS receiver, which were sampled every 1s, are considered for computation and analysis of amplitude scintillations, ROT, and ROTI. Figure 2 specifies S� for all the available satellites on 2020-06-07 for the L1 frequency. The highest scintillations are observed from the satellite PRN 9. Therefore, the values of S� and ROTI are compared to the PRN 9 values. Fig. 2. S� index on 2020-06-07 for the L1 frequency. Figure 3 shows the comparison plot for S� and ROTI for the PRN-9, which exhibits maximum scintillations. Figure 4 shows the correlation between ROTI and S�. The x axis specifies the epoch time for every hour and the y axis the correlation coefficient between S� and ROTI. The highest correlation coefficient is 1 and the minimum is higher than the cutoff level. As a result, the correlation coefficient is robust due to the persistent link between ROTI and S�. Figure 5 depicts the comparison of S� for all the available NavIC satellites on 2020-06-07 for the L5 frequency. The maximum observed S� index is 0.0205 for IR-02, which is an indication of no scintillations on 2020-06-07 for the L5 frequency. The highest S� values are observed in the IR-02 satellite. Therefore, the values of S� and ROTI are compared to the IR-02 values. Engineering, Technology & Applied Science Research Vol. 13, No. 3, 2023, 10936-10940 10938 www.etasr.com Nimmakayala & Vemuri: Analysis of Ionospheric Scintillations using GPS and NavIC Combined … Fig. 3. Comparison of �� and ROTI on 2020-06-07 for PRN-9 (L1 frequency). Fig. 4. Correlation of of S� and ROTI on 2020-06-07 for PRN-9 (L1 frequency). Fig. 5. S4 amplitute variations on 2020-06-07 for L5 frequency. Fig. 6. Comparison of S� and ROTI on 2020-06-07 for IR-02 (L5 frequency). In Figure 6, the maximum S� index value of 0.0078 and the maximum ROTI index value of 0.1911 are observed. From these, it can be concluded that no scintillations are observed on the L5 frequency for IR-02. Figure 7 shows the correlation between S� and ROTI. The maximum correlation of 1 is observed many times, which is an indicative of the similarity between S� and ROTI. Figure 8 shows the S� variations of all visible NavIC satellites on 2020-06-07 for the � frequency. The maximum observed S� index is 0.0294 for IR-02, which is an indication of no scintillations there on 2020-06-07 for the S frequency. The highest S4 values are observed in the IR-02 satellite. Therefore, all S� and ROTI values are compared with the ones of the IR-02 satellite. Fig. 7. Correlation of S� and ROTI on 2020-06-07 for IR-02 (L5 frequency). Fig. 8. S4 amplitute variations on 2020-06-07 for S1 frequency. Fig. 9. Comparison of S4 and ROTI on 2020-06-07 for IR-02 (S1 frequency). Figure 9 shows that the maximum S� index value is 0.0078 and the maximum ROTI index value is 0.1911. From this it can Engineering, Technology & Applied Science Research Vol. 13, No. 3, 2023, 10936-10940 10939 www.etasr.com Nimmakayala & Vemuri: Analysis of Ionospheric Scintillations using GPS and NavIC Combined … be concluded that no scintillations are observed on the S frequency for IR-02. Figure 10 shows the correlation between S4 and ROTI for NavIC satellite IR-02 on 2020-06-07 on the S frequency. Maximum correlation is observed many times, which indicates that ROTI can be considered as an alternative metric for the S� index. Figure 11 shows the skyplot of all the visible GPS and NavIC satellites at 15:00hrs of 2020-06-07 for the IGS (International GNSS Service) station which is located at Hyderabad (Lat/Long:17.417 0 N/78.551 0 E). Fig. 10. Correlation of S4 and ROTI on 2020-06-07 for IR-02 (S1 frequency). Fig. 11. Sky plot of visible satellites at 15:00 hrs on 2020-06-07. V. CONCLUSIONS In this article, the analysis of GNSS trans-ionospheric signals for a typical day (2020-06-07) is presented. In order to identify the ionospheric perturbations, the amplitude of the scintillations S�, the rate of the change of total electron content (ROT), and its index (ROTI) are computed on L1, L5, and S frequencies using GPS and NavIC combined constellation. Based on the correlation coefficient, maximum similarity between the S� and ROTI is observed, and, thus, ROTI can be considered as an alternative index for the investigation of the amplitude scintillations. It is also observed that the ionospheric scintillations are robust only on L1 frequency of the GNSS trans-ionospheric signals and no ionospheric scintillations are observed on the L5 and S frequencies. As the L5 and S are the frequencies for NavIC constellation, it is concluded that the ionospheric scintillation effects are minimum on the NavIC signal. Hence, it is concluded that the ionosphere scintillation effects can be minimized with higher frequencies, like L5 and S1. ACKNOWLEDGMENT The authors are grateful to the Navigation and Communication Centre (NCRC), Chaitanya Bharathi Institute of Technology (CBIT), Hyderabad, Space Applications Centre (SAC), and ISRO Ahmadabad for allowing the use of the IGS receiver data for this research project. REFERENCES [1] S. K. Pagoti and S. I. D. Vemuri, "Development and performance evaluation of Correntropy Kalman Filter for improved accuracy of GPS position estimation," International Journal of Intelligent Networks, vol. 3, pp. 1–8, Jan. 2022, https://doi.org/10.1016/j.ijin.2022.01.002. [2] Ch. Rajasekhar, V. B. S. Srilatha Indira Dutt, and G. Sasibhushana Rao, "Weighted GDoP for improved position accuracy using NavIC and GPS hybrid constellation over Indian sub-continent," International Journal of Intelligent Networks, vol. 2, pp. 42–45, Jan. 2021, https://doi.org/10.1016/j.ijin.2021.06.001. [3] O. M. Mubarak, "The Effect of Carrier Phase on GPS Multipath Tracking Error," Engineering, Technology & Applied Science Research, vol. 10, no. 5, pp. 6237–6241, Oct. 2020, https://doi.org/10.48084/etasr.3578. [4] T. L. Dammalage, "The Effect of Multipath on Single Frequency C/A Code Based GPS Positioning," Engineering, Technology & Applied Science Research, vol. 8, no. 4, pp. 3270–3275, Aug. 2018, https://doi.org/10.48084/etasr.2206. [5] A. Hussain, F. Akhtar, Z. H. Khand, A. Rajput, and Z. Shaukat, "Complexity and Limitations of GNSS Signal Reception in Highly Obstructed Enviroments," Engineering, Technology & Applied Science Research, vol. 11, no. 2, pp. 6864–6868, Apr. 2021, https://doi.org/10.48084/etasr.3908. [6] P. Sirish Kumar, V. B. S. Srilatha Indira Dutt, and L. Ganesh, "Implementation of new navigation algorithm based on cross- correntropy for precise positioning in low latitude regions of South India," International Journal of Speech Technology, vol. 23, no. 4, pp. 747–756, Dec. 2020, https://doi.org/10.1007/s10772-020-09727-6. [7] O. E. Abe, Y. O. Migoya-Orué, and S. M. Radicella, "Correlation analysis of normalized ROTI and S4 as observed during different geomagnetic conditions of the 2013 September equinox over the stations within EIA African sector," Advances in Space Research, Oct. 2022, https://doi.org/10.1016/j.asr.2022.09.058. [8] S. Mulugeta and T. Kassa, "Nighttime ionospheric irregularities inferred from Rate of Total Electron Content Index (ROTI) values over Bahir Dar, Ethiopia," Advances in Space Research, vol. 67, no. 4, pp. 1261– 1266, Feb. 2021, https://doi.org/10.1016/j.asr.2020.11.030. [9] T. K. Gogie, "The Rate of Ionospheric Total Electron Content Index (ROTI) as a Proxy for Nighttime Ionospheric Irregularity Using Ethiopian Low-Latitude GPS Data," Geomagnetism and Aeronomy, vol. 61, no. 3, pp. 464–475, May 2021, https://doi.org/10.1134/S0016793221030051. [10] C. Li, C. M. Hancock, N. A. S. Hamm, S. V. Veettil, and C. You, "Analysis of the Relationship between Scintillation Parameters, Multipath and ROTI," Sensors, vol. 20, no. 10, Jan. 2020, Art. no. 2877, https://doi.org/10.3390/s20102877. [11] K. Kotulak, I. Zakharenkova, A. Krankowski, I. Cherniak, N. Wang, and A. Fron, "Climatology Characteristics of Ionospheric Irregularities Described with GNSS ROTI," Remote Sensing, vol. 12, no. 16, Jan. 2020, Art. no. 2634, https://doi.org/10.3390/rs12162634. [12] W. Li, S. Song, and X. Jin, "Ionospheric Scintillation Monitoring With ROTI From Geodetic Receiver: Limitations and Performance Evaluation," Radio Science, vol. 57, no. 5, 2022, Art. no. e2021RS007420, https://doi.org/10.1029/2021RS007420. [13] D. Zhao, W. Li, C. Li, C. M. Hancock, G. W. Roberts, and Q. Wang, "Analysis on the ionospheric scintillation monitoring performance of ROTI extracted from GNSS observations in high-latitude regions," Engineering, Technology & Applied Science Research Vol. 13, No. 3, 2023, 10936-10940 10940 www.etasr.com Nimmakayala & Vemuri: Analysis of Ionospheric Scintillations using GPS and NavIC Combined … Advances in Space Research, vol. 69, no. 1, pp. 142–158, Jan. 2022, https://doi.org/10.1016/j.asr.2021.09.026. [14] Q. Li, Y. Zhu, K. Fang, and J. Fang, "Statistical Study of the Seasonal Variations in TEC Depletion and the ROTI during 2013–2019 over Hong Kong," Sensors, vol. 20, no. 21, Jan. 2020, Art. no. 6200, https://doi.org/10.3390/s20216200. [15] Ming O. U., Jiayan W. U., Longjiang C., and Weimin Z., "Global Ionospheric TEC and ROTI Variations during a Moderate Geomagnetic Storm," Chinese Journal of Space Science, vol. 41, no. 6, pp. 887–897, Nov. 2021, https://doi.org/10.11728/cjss2021.06.887. [16] X. Luo, S. Gu, Y. Lou, L. Cai, and Z. Liu, "Amplitude scintillation index derived from C/N0 measurements released by common geodetic GNSS receivers operating at 1 Hz," Journal of Geodesy, vol. 94, no. 2, Feb. 2020, Art. no. 27, https://doi.org/10.1007/s00190-020-01359-7.