Next Article in Journal
The Relocation Problem of Field Calibrated Low-Cost Sensor Systems in Air Quality Monitoring: A Sampling Bias
Previous Article in Journal
Kinematic ME-MAFA for Pseudolite Carrier-Phase Ambiguity Resolution in Precise Single-Point Positioning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Study of the Seasonal Variations in TEC Depletion and the ROTI during 2013–2019 over Hong Kong

1
School of Electronics and Information Engineering, Beihang University, Beijing 100191, China
2
Aviation Data Communication Corporation, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(21), 6200; https://doi.org/10.3390/s20216200
Submission received: 10 October 2020 / Revised: 29 October 2020 / Accepted: 29 October 2020 / Published: 30 October 2020
(This article belongs to the Section Remote Sensors)

Abstract

:
Equatorial plasma bubbles (EPBs) can cause large total electron content (TEC) gradient magnitudes and significant density irregularities. In this paper, depletions and irregularities due to EPBs are identified by using the Global Positioning System (GPS)-TEC time series extracted from nine Global Navigation Satellite System (GNSS) stations over Hong Kong near the equatorial ionization anomaly (EIA) crest region from 2013 to 2019. The correlation analyses between the daily variation in the rate of TEC change index (ROTI) and that of the EPB occurrence rate, depth, and duration are presented. The monthly EPB occurrence rate, depth, duration, and ROTI show strong seasonal variations, with maxima during equinoctial seasons, especially during the moderate-to-high solar activity years of 2013–2016. Furthermore, two seasonal asymmetries can be clearly seen for these parameters from 2013 to 2016. The EPB occurrences rate, depth, and duration vary annually with the solar radio flux at 10.7 cm (F10.7) index. The correlation analyses of the EPB occurrence rate, depth, and duration are found to be much more strongly correlated with the F10.7 index on an annual basis than on a monthly basis. The correlation analysis of monthly variations shows the impacts of solar activity on EPB occurrence, depth, and duration are seasonally dependent, which is significantly greater in the equinoctial seasons and summer than in winter.

1. Introduction

Equatorial plasma bubbles (EPBs) refer to plasma depletions frequently observed in the evening and are the major sources of electron density irregularities in equatorial and low-latitude regions [1,2]. The plasma density irregularities within the bubbles can scatter and diffract radio waves, leading to rapid fluctuations in the signal amplitude and phase in a received electromagnetic wave [3]. The significant magnitude of total electron content (TEC) gradients associated with the side walls of bubbles or with plasma density irregularities could potentially cause very large Global Navigation Satellite System (GNSS) positioning errors [4,5].
Physically, the occurrence of an EPB largely depends on the mechanism for producing initial density perturbations and the conditions affecting the Rayleigh–Taylor (R-T) growth rate [6]. Through the R-T instability mechanism acting on the steep density gradient of the F region bottom side, equatorial spread F (ESF) develops shortly after sunset and propagates both upward and toward off-equatorial latitudes [7,8,9]. As the lower-density ionosphere moves upward into the higher-density topside ionosphere, a plasma bubble is created. Most ionospheric density irregularities are created by velocity shear mixing within the bubble gradients along with the formation of drift waves [1,3]. The irregularities start to decay first with the shortest scale length as the night progresses, and then they disappear when solar ionization refills the flux tube at sunrise [3].
Nighttime equatorial irregularities have been observed in ionograms as ESF for a long time in history [1,9]. The EPB/ESF occurrence was found to be dependent on the solar cycle, season, local time (LT), and longitudinal sector based on many studies using observational techniques such as ionosondes [8], very-high-frequency (VHF) radar systems [10], and in situ satellites [11]. Since EPBs cause depletions and irregularities in total electron content (TEC) time series when satellite-to-receiver ray paths pass through the region of plasma bubbles, studies of EPB occurrence based on TEC or rate of TEC change index (ROTI) have also been conducted in recent years [12].
Previous studies have demonstrated that the easily accessible ROTI is an appropriate indicator to describe the strength of ionospheric irregularities [13]. Monthly, seasonal, semiannual, and latitudinal variations in EPB occurrence based on the Global Positioning System (GPS) ROTI over five global regions between 2000 and 2006 were analyzed [14]. The regional diurnal and seasonal variability in the GPS-TEC depletions over South America in 2008 were studied [15]. The occurrence rate, time duration, and depth of TEC depletions and their dependence on the LT and season from 2002 to 2014 were analyzed at the global scale [16]. The characteristics and causes of seasonal variations in EPB over Hong Kong from 2014 to 2017 were studied [17].
Through previous studies, the characteristics of EPB occurrence were observed to be seasonally and longitudinally dependent [14,18]. Based on ROTI observations, Nishioka [14] showed that the two maxima were seen in the annual occurrence variations in African, Asian, and Atlantic regions. Furthermore, the maxima of occurrence rates in the Asian region were in March/April and September/October. The analysis of EPB occurrences over Hong Kong found a maximum during equinoctial months and a minimum during December from 2001 to 2012, except during the solar minimum in 2007–2009 [19]. Although a geometric relationship between the geomagnetic field line and the sunset terminator line was proposed to successfully explain the observed maxima in the occurrence rates during spring and autumn [18], two asymmetries of EPB occurrences observed in previous studies reveal that other factors need to be studied [14,17,19].
An EPB could affect different aspects of received GPS signals. In addition to temporal and spatial variations in EPB occurrences based on solitary observations, potential correlations among multiple observations, such as scintillation indices, ROTIs, and TEC depletions, have been studied. Occurrences of deep GPS-TEC depletions over China’s low-latitude region during the equinoctial seasons from 2011 to 2012 were found to coincide qualitatively with the scintillation index and ROTI [20]. The correlation coefficients between the ROTIs and scintillation indices over Hong Kong from 2012 to 2013 were found to be 0.6 if the data from all GPS satellites were used together [21]. The statistical relationship between the intensity scintillation index S4 and the ROTI was analyzed over the Indian equatorial region [22].
Correlation studies are appealing not only because the very large number of conventional geodetic receivers provides much more data coverage for research opportunities than the limited number of ionospheric scintillation monitoring receivers (ISMRs) [23] but also because there exists a straightforward relationship between the ROTI, irregularity strength, and intensity scintillation index [24]. For instance, the degree to which ionospheric irregularities produce scintillations is dependent on the frequency of the signal compared to the plasma frequency, propagation geometry, and strength of the irregularities [3,24].
In this paper, depletions and irregularities due to EPBs are identified by using the GPS-TEC time series extracted from nine GNSS stations over the Hong Kong region from 2013 to 2019. The temporal variations in the occurrence rate, depletion depth, duration, and ROTI are analyzed on daily, monthly, and semiannual time scales. The relation between TEC depletion and irregularities is studied by conducting correlation analyses between the ROTI and three characteristic parameters of EPBs: the occurrence rate, depth, and duration. The seasonal and yearly variations in the EPB occurrence rates, depths, and durations and their solar activity dependences are presented. The paper shows that the seasonal dependence and asymmetries are manifested not only in the monthly variations of occurrence rate and ROTI but also in the variations of depth and duration. In addition, EPB occurrence rate, depth, and duration are found to be impacted greatly by solar activity and the impact degree is more significant on an annual than on a monthly basis.

2. Data and Methods

2.1. Dataset

The paper analyzes the EPBs based on ionospheric slant TEC (STEC) values extracted from GNSS observations. The GNSS observations were recorded from continuously operating reference stations affiliated with the Hong Kong satellite positioning reference station network (SatRef). The distributions of the nine GNSS stations used in the study are shown in Figure 1a. The nine stations were chosen according to the data availability from 2013 to 2019. The station codes, geographic coordinates and geomagnetic latitudes of these stations are shown in Table 1.
Figure 1b shows the solar cycle variations from 1975 to 2020, corresponding to solar cycles (SCs) 21 to 24. The blue dots denote the solar radio flux at 10.7 cm (F10.7) observations, and the red line denotes the 81 day moving averages. The solar activities became weaker from SC 22 to SC 24, the latter of which had the weakest solar activity during the period. The studied period from 2013 to 2019 covers the solar maximum and the solar minimum in SC 24. The solar activity, i.e., the F10.7 index, increased from 2013 to 2014 and then decreased year by year from 2014 to 2019.

2.2. Detection Method of Plasma Bubbles

When a GNSS signal propagates through the ionosphere, because of the dispersive nature of the ionosphere the free electrons induce a group delay and a phase advance, which is inversely proportional to the square of the frequency of the signal. This feature could be used to estimate the ionospheric delay by making geometry-free signal combinations from the GNSS dual-frequency signals [25]. Specifically, ionospheric delays from a satellite to a receiver are obtained by differencing the dual-frequency observations between the pseudoranges ( P 1 and P 2 ) and between the phases ( L 1 and L 2 ) [25,26,27]. The equations are as follows:
I p = def P 2 P 1 γ 12 = STEC + B r + B s + ε 2 ε 1 / γ 12  
I L = def L 1 L 2 γ 12 = STEC + λ 1 N 1 λ 2 N 2 + υ 1 υ 2 / γ 12
I p 1 _ L 1 = def P 1 L 1 / 80.6 / f 1 2 = STEC +   ε 1 λ 1 N 1 + υ 1 / 80.6 / f 1 2
γ 12 = 40.3 / f 2 2 40.3 / f 1 2
where P i and L i are the code and carrier phase measurements, respectively, on frequency f i , and γ 12 is a frequency-dependent parameter. B r and B s are unknown hardware biases from the receiver and satellite, respectively, which are relatively stable over several hours to a number of days [28,29]. N i is the carrier phase ambiguity parameter at a wavelength of λ i . ε i and υ i are multipath and noise terms of the code and carrier, respectively, on frequency f i .
In the above equations, all the ionospheric observables are converted to ionospheric delays expressed in TEC units (TECU), which are convenient for comparing the ionospheric delays from different frequency combinations; 1 TECU is equivalent to 1016 electrons/m2 column. I p contains the absolute ionospheric delay plus hardware biases. I L and I p 1 _ L 1 are impacted by unknown ambiguities in the phase measurement. I p 1 _ L 1 relies only on measurements on the first GPS frequency.
For various reasons, a number of cycle slips may occur and lead to sudden discontinuities in carrier phase observations. The carrier phase discontinuity could be detected and repaired by using the method from TurboEdit [26]. Then, to remove effects of unknown phase ambiguities, the leveled carrier-derived ionospheric observable I ϕ and leveled code-minus-carrier I cmc are derived as follows [30,31]:
I ϕ = def I L + I p I L arc = STEC + B r + B s + ε 2 ε 1 arc / γ 12
I cmc = def I p 1 _ L 1 + I p I p 1 _ L 1 arc = STEC + B r + B s + ε 2 ε 1 arc γ 12 + ε 1 ε 1 arc 80.6 / f 1 2
where arc means averaging along a continuous arc. The leveled carrier-derived ionospheric observable I ϕ has a higher precision than I p and I cmc , and the latter will be used for validation purposes. Although unknown hardware biases are contained in ionospheric observables, they will be canceled out as we are only concerned with the relative value of the STEC in the study. An elevation mask of 15° is used for each satellite and receiver pair.
The temporal rate of ionospheric delay (rate of TEC, ROT) is computed from the difference in the ionospheric delay I ϕ divided by the time interval, usually expressed in TECU per minute as follows [12]:
ROT = I ϕ t 2 I ϕ t 1 t 2 t 1
The ROTI is usually used to investigate ionospheric irregularities. The ROTI is computed as the standard deviation of the ROT within 5 min windows using the following equation [12]:
ROTI = ROT 2 ROT 2
When the satellite–receiver ray path passes through a plasma bubble, a sudden decrease in the STEC value and then a recovery appears in the GPS-TEC time series. The bubble can be detected using STEC data using the method proposed by Pradipta et al. [32]. In this method, the TEC trend variation is extracted from the original STEC ( I ϕ ) time series by using mechanical rotary motion-implicit terrain (MRMIT) filtering [32]. Then, the detrended STEC (dTEC) time series containing the plasma bubble is computed as follows:
dTEC = STEC STEC f
where STEC f is the TEC trend variation using MRMIT filtering. To identify plasma bubbles, the threshold of minimum TEC depletion is set to 4 TECU. The depletion due to an actual plasma bubble is usually larger than 4 TECU according to previous studies [15,16]. The threshold is chosen to exclude the impact of other disturbances such as traveling ionospheric disturbance.
Three parameters can be identified using the dTEC time series. The depth of an EPB is defined as the maximum magnitude of the TEC depletion derived from the dTEC. The occurrence time of an EPB is defined as the epoch when the dTEC reaches the maximum depletion. The duration of an EPB is computed as the period from the start time of the dTEC decrease to the end time of dTEC recovery.
Figure 2 shows a typical case of an EPB observed by satellite GPS 19 (PRN) from the Hong Kong station HKWS on 9 March 2014 (quiet day with Kp < 1 and Dst ≥ −6 nT). The plasma bubble observed by GPS 19 has a duration of 88 min and a depth of 77 TECU. The large ROTI value corresponds well to the EPB occurrence period, which indicates that severe electron density irregularities inside plasma bubbles cause significant fluctuations in the TEC [14].

3. Results

3.1. Diurnal Variation in Plasma Bubbles over the HKWS Station

The GPS-STEC from station HKWS was analyzed to detect plasma bubbles during the 7 years from 2013 to 2019. As a GPS satellite–receiver ray path can be impacted by plasma bubbles several times during an observation arc, multiple EPB occurrence events can be observed for a GPS-TEC time series. Figure 3 presents the daily variation in EPB occurrence event numbers and the number of GPS satellites with EPB observations. Overall, the daily variations in both numbers showed obvious consistent tendencies from 2013 to 2019. They both vary greatly with solar cycle and season. The correlation coefficient of both numbers was found to be 0.95. Statistical testing for the significance of the relationships resulted in a small p-value (p < 0.0001) between EPB occurrence event number and the number of GPS satellites with EPB observations, which indicates a strong correlation.
In total, 2668 EPB occurrence events were observed on the 431 quiet days and 230 EPB occurrence events were observed on the 39 active days (Dst < −50 nT). Figure 4 presents the observed distribution of the EPB depths, ROTI, LT at maximum depth, and durations from 2013 to 2019. It can be seen that the maximum depth can reach as large as 100 TECU (in the STEC domain), and most of the ROTI values are larger than 0.5 TECU/min, the LTs of the EPB occurrences range from 19:00 to 05:00 LT, and the durations can reach as long as 120 min.
Most EPB occurrences (2668 out of 2898 events) were observed during the quiet period (431 out of 470 days), and furthermore, EPB occurrences can be enhanced or suppressed during the storm period through complicated mechanisms; thus, in the following analysis, only EPBs observed on the geomagnetically quiet days (Dst ≥ −50 nT) were used.
Since we care more about how many satellites would be affected by EPBs, we computed the daily EPB occurrence rate for a particular day using the following equation:
d a i l y   o c c u r r e n c e   r a t e % = N u m b e r   o f   G P S   s a t e l l i t e s   w i t h   E P B s   o n   t h e   d a y T o t a l   n u m b e r   o f   v i s i b l e   G P S   s a t e l l i t e s   o n   t h e   d a y   × 100
We studied the relationship between the variation in the daily occurrence rate based on STEC depletions and the variation in the maximum nighttime ROTI. Figure 5 presents diurnal variations in the daily occurrence rates and ROTI over the HKWS station during 2013–2019. The ROTI values are expressed in units of 0.1 TECU/min for better comparison with the occurrence rate. The occurrence rate and ROTI increased or decreased in a similar manner. The correlation between the daily occurrence rate and the ROTI was found to be strong (R = 0.85).
Similarly, the variations in the daily maximum EPB depth and duration were compared with the ROTI. Diurnal variations in the EPB depth and ROTI are presented in Figure 6. The variation in EPB depths was strongly correlated with the ROTI, and the correlation coefficient was 0.85. As shown in Figure 7, the variations in EPB duration and ROTI had a modest correlation coefficient (R = 0.72) and the durations under higher ROTI were slightly more scattered.

3.2. Monthly and Seasonal Variation in EPBs

Monthly EPB variations from 2013 to 2019 were analyzed for the nine Hong Kong stations. As shown in Figure 8a, the cumulative numbers of days with EPB occurrence in each month for the nine stations were in agreement with each other. The monthly nighttime (from 20:00 to 02:00 LT) ROTI maxima were consistent for the nine stations, as shown in Figure 8b.
To study the monthly and seasonal variations in the EPB occurrences and ROTI, the monthly EPB occurrence rate per station was computed using the following equation:
M o n t h l y   o c c u r r e n c e   r a t e   % = N u m b e r   o f   d a y s   w i t h   E P B s   i n   t h e   m o n t h T o t a l   n u m b e r   o f   d a y s   i n   t h e   m o n t h   × 100
The monthly EPB occurrence rate over the Hong Kong region was then obtained by averaging the monthly EPB occurrence rate per station across the nine stations. Similarly, by averaging the monthly maximum nighttime ROTI per station, the monthly ROTI was derived over the region.
As shown in Figure 9, there were good correspondences between the EPB occurrence rate and ROTI. The correlation coefficient of the monthly variation was high (R = 0.83), slightly lower than that (R = 0.85) of the daily variation. The values of the occurrence rate and ROTI were obviously larger in the two equinoctial seasons than in the two solstitial seasons during the years from 2013 to 2019 except 2018. During 2018, the difference among the four seasons was not significant. Furthermore, there are two seasonal asymmetries: (1) solstitial asymmetry and (2) equinoctial asymmetry. For the solstitial asymmetry, the occurrence rate was obviously larger in June than in December during 2013–2016. In a similar manner, the ROTI was larger in June than in December during 2014–2016. However, in 2013, the occurrence rate in June was larger than that in December, whereas the ROTI in December was larger than that in June. This phenomenon means that solstitial asymmetry manifested differently in the EPB occurrence rate and ROTI. For the equinoctial asymmetry, both the occurrence rate and ROTI in spring months (March/April) were larger than those in the autumn months (September/October) during 2013–2016.
The monthly variation in the maximum EPB depth was analyzed and compared with the ROTI as shown in Figure 10. Generally, EPB depths showed obvious equinoctial maxima throughout 2013–2019 except in 2018. The EPB depths also presented solstitial and equinoctial asymmetry similar to that observed in the ROTI during 2013–2016. The correlation in the monthly variations between the EPB depths and the ROTI was very high (R = 0.94), which is even larger than that based on diurnal variations (R = 0.85).
The correlation analysis of the monthly variations between the maximum EPB duration and the ROTI is shown in Figure 11. The correlation coefficient is 0.85, which is larger than that based on daily variations (R = 0.72). Two maxima can be seen in the equinoctial seasons from 2013 to 2019. There exists solstitial asymmetry during 2014–2016, whereas there is no obvious difference in 2013. Equinoctial asymmetry similar to that of the ROTI is clearly seen during 2014 to 2016. However, the maximum duration in autumn was larger than that in spring in 2013.

3.3. Effect of Solar Activity

Annual variations in the EPB occurrence rates, depths, and durations were derived for the 7 years from 2013 to 2019 in this study. Figure 12 shows the semiannual variations during the 7 years. The semiannual occurrence rate was derived as the ratio of the number of days with EPB occurrences to the number of observation days in a half year from January to June and from July to December. The semiannual F10.7 index was derived as the average daily value across a half year. The semiannual EPB depth and duration were obtained by averaging the monthly maximum values across a half year.
Figure 12 shows that the semiannual occurrence rate, EPB depth, and duration vary with the mean F10.7 index. A clear SC pattern was observed from 2013 to 2019. In the ascending phase of SC 24 from 2013 to 2014, a small increase was seen in the EPB occurrence rate and depth, whereas a more obvious increase was observed in the EPB duration. These values reached a maximum in the high solar activity year of 2014, and then they decreased gradually from 2014 to the low solar activity year of 2019. Furthermore, the correlation coefficients between the semiannual EPB occurrence rate, depth, and duration and the F10.7 index were 0.96, 0.94, and 0.95, respectively, indicating that they are closely correlated with the F10.7 index on an annual basis.
To study the seasonal impact of solar activity on EPB occurrences, depths, and durations, a correlation analysis was conducted for the monthly occurrence rates grouped into 3 seasons, namely, summer months (May to August), equinoctial months (March, April, September, and October), and winter months (January, February, November, and December).
The monthly variation in the occurrence rates of EPBs and the monthly mean F10.7 index are shown in Figure 13. A clear seasonal variation was seen in the occurrence rates of EPBs, with the maximum occurrence found in March 2014. The correlation analysis of monthly EPB occurrences against the monthly mean F10.7 index was carried out for the grouped months. The seasonal analysis shows that EPB occurrences were influenced more by solar activity during summer (R = 0.89) and the equinoctial seasons (R = 0.83) than during winter (R = 0.40). Furthermore, the impact of solar activity on EPB occurrences varied more obviously on an annual basis (R = 0.96) than on a monthly basis (R = 0.59).
Similarly, the monthly maximum EPB depths along with the F10.7 index are presented in Figure 14. The maximum depths occurred in February 2014. The seasonal analysis shows that EPB depths were influenced more by solar activity during equinoctial seasons (R = 0.88) and summer (R = 0.78) than during winter (R = 0.57). Furthermore, the impacts of solar activity on EPB depths varied more obviously on an annual (R = 0.94) basis than on a monthly (R = 0.68) basis.
The monthly maximum durations of EPBs along with the F10.7 index are shown in Figure 15. The maximum duration occurred in April 2014. The seasonal analysis shows that EPB durations were influenced more by solar activity during the equinoctial seasons (R = 0.89) and summer (R = 0.83) than during winter (R = 0.45). Furthermore, the impact of solar activity on EPB duration varied more obviously on an annual basis (R = 0.95) than on a monthly basis (R = 0.63).

4. Discussion

4.1. Diurnal Variation in EPB Occurrence, Depth, and Duration

Based on the analysis from a single station (HKWS) during 2013–2019, the diurnal variations in the EPB occurrences, depths, and durations corresponded well with that of ROTI. They showed similar patterns in maximum and minimum in seasonal scale and solar cycle. On the diurnal scale, the EPB occurrence rate, depth, and duration were highly correlated with the ROTI, with correlation coefficients equal to 0.85, 0.85, and 0.72, respectively.
The magnitude of EPB depth determines the degree of decrease in plasma density of the bubbles. The EPB duration reflects the degree of spatial or temporal scale that a plasma bubble would impact. The ROTI is the metric used to measure the strength of ionospheric irregularity inside or around the region of plasma bubbles with a much larger gradient [5,12]. Ionospheric irregularities occurring inside or around plasma bubbles usually lead to large ROTIs [13]. The relationship between the three EPB parameters and the ROTIs reflects the degree of correlation between electron density depletions and density irregularities. The significant correlation coefficients indicate that greater value in the ROTI is usually closely associated with larger occurrence rate, depth, and duration of EPBs.

4.2. Monthly and Seasonal Variation in EPBs

The paper analyzed monthly EPB variations from 2013 to 2019 for the nine Hong Kong stations. The monthly variations in EPBs occurrence rates, depths, and durations were generally consistent with those of the ROTI from 2013 to 2019, and they all showed strong seasonal variations, especially during moderate-to-high solar activity years 2013–2016. The maximum observations in EPB occurrence and ROTI during spring and autumn in this paper were consistent with those from previous reports [14,17,19]. Additionally, the paper showed that the maximum observations in EPB depth and duration also occurred in the spring and autumn.
The correlation coefficients between these parameters and the ROTI were 0.83, 0.94, and 0.85, respectively. Although the correlation coefficient of the monthly variation in EPB occurrences was slightly smaller than that of the daily variation, the correlation coefficients of the monthly variations between the EPB depth and duration and the ROTI were greater than those of the daily variations. This finding indicates that daily variations are more scattered than monthly variations, which might be because subtle variations in various factors that contribute to the nocturnal variations in the EPBs and ROTI are smoothed out during monthly correlation analysis.
Using twenty-three GPS receivers around the magnetic dip equator, Nishioka et at. [14] showed that occurrence rates based on ROTI observations achieved maxima in March/April and September/October in the Asian region. Using GPS-TEC data from 2001 to 2012 over Hong Kong, Kumar et al. [19] reported two clear seasonal maxima in the EPB occurrences corresponding to the spring and autumn months. Tang et al. [17] reported that the seasonal occurrence rate was significantly larger in the two equinoctial seasons than in the solstitial seasons in 2014 and 2015. In the low solar activity year of 2017, the difference in the occurrence rate between summer and the two equinoctial seasons was not obvious. The correlation coefficients between the ROTIs and scintillation indices over Hong Kong from 2012 to 2013 were found to be 0.6 if the data from all GPS satellites were used together [21].
The seasonal variation in plasma bubbles is controlled by the geometry between the geomagnetic field line and the sunset terminator line [18]. In the equinoctial seasons, the E region conductivity is reduced most rapidly due to the closest alignment of the solar terminator with the magnetic meridian, and the rapid reduction is responsible for the enhancement of the vertical plasma drift, which creates a favorable condition for the formation of plasma bubbles.
The paper presented that the seasonal asymmetries were not only manifested in the monthly variation of occurrence rate but also in the variation of ROTI, depth, and duration of EPB during 2014–2019, especially during the moderate-to-high solar activity years of 2013–2016. The patterns of the two seasonal asymmetries were generally consistent with those from previous studies [17]. However, during 2013, the maximum ROTI, depth, and duration were larger in December than in June, and the maximum duration in September was larger than that in March.
The seasonal asymmetries have been studied by previous researchers. Tang et al. [17] reported that the seasonal occurrence rate was greater in spring than in autumn, especially during 2014–2015. Furthermore, Tang et al. [17] observed that seasonal occurrences were more frequent in summer than in winter during 2014 and 2015. Kumar et al. [19] found that the EPB occurrence rate had a higher value in September than in March during the high solar activity year of 2001.
The slight discrepancy between the result in the paper and previous studies indicated complicated mechanisms were involved in interpreting the seasonal asymmetries. The equinoctial asymmetry was assumed to be related to neutral wind in the F region [14] or asymmetry in the electron density distribution [17,18]. The solstitial asymmetry was interpreted by monthly variations in flux-tube-integrated conductivities in the F regions [14] or the seeding mechanism of thunderstorm-driven atmospheric gravity waves [17]. Previous studies [14] also suggested that the discrepancy in the observations of the occurrence rates could be due to the difference in measurement techniques and/or observational periods and/or the location.
Further analysis (not shown here) based on 90 percent instead of the maximum monthly ROTI, occurrence rate, depth, and duration indicates that the values are consistently larger in June than in December during 2013–2016, and they are larger in spring than in autumn. Therefore, the discrepancy in the asymmetries could be attributed to the difference in observation techniques.

4.3. Effect of Solar Activity

Annual variations in the EPB occurrence rates, depths, and durations were derived for the 7 years from 2013 to 2019 in this study. The solar cycle dependence was clearly observed from 2013 to 2019. They were highest in the high solar activity year of 2014, and then they decreased gradually from 2014 to the low solar activity year of 2019. The correlation coefficients between the semiannual EPB occurrence rate, depth, and duration and the F10.7 index were found to be 0.96, 0.94, and 0.95. This finding indicates that the EPB occurrence rate, depth, and duration are strongly dependent on solar activity on an annual basis.
Furthermore, the seasonal impacts of solar activity on EPB occurrence, depth, and duration indicated that the impacts of solar activity were most pronounced in the equinoctial months than during winter. Moreover, the differences on impact degree of solar activity were compared respectively between annual and monthly bases for EPB occurrences (annual R = 0.96 vs. monthly R = 0.59), depths (annual R = 0.94 vs. monthly R = 0.68), and durations basis (annual R = 0.95 vs. monthly R = 0.63). The results show that EPBs are influenced more obviously on an annual basis than on a monthly basis.
Tang et al. [17] showed that EPB occurrence from 2014 to 2017 was highest in 2014. Kumar et al. [19] indicated that the EPB occurrence rate varied with solar activity from 2001 to 2012 and reached a maximum during the high solar activity year of 2002. They reported that the correlation coefficient between the annual total number of EPB occurrences and the annual mean F10.7 index was 0.92. Moreover, they observed that the effects of solar activity on EPB occurrences were more pronounced on an annual basis than on a monthly basis. They found that the annual correlation coefficient of 0.92 between the annual number of EPB occurrences and the annual mean F10.7 index was significantly larger than the monthly correlation coefficient of 0.63 between the monthly occurrences and the monthly mean F10.7 index. For the seasonal impacts of solar activity on the number of EPB occurrences, their seasonal analysis showed that EPB occurrences were influenced more during spring and autumn (R = 0.80) than during summer (R = 0.62) and winter (R = 0.68).
The paper shows that the results of solar impact on EPB occurrences are consistent with previous studies. Additionally, the analysis presents that EPB depths and durations are also impacted greatly by solar activity, and the impact degree is more significant on an annual than on a monthly basis. The solar cycle dependence of ionospheric irregularities on the F10.7 index is explained in terms of the enhancement of the pre-reversal electric field amplitude with an increase in the solar UV radiation intensity [33]. The TEC depletion depth and duration characterize different aspects of EPBs, and they are controlled by the same electron dynamical process as the EPBs, and thus are indirectly related to solar activity in a similar manner.

5. Conclusions

The paper analyzed the features of plasma bubbles based on ionospheric STEC data extracted from GNSS observations over Hong Kong from 2013 to 2019 covering the solar maximum and the solar minimum in SC 24. By using MRMIT filtering, the detrended STEC time series containing plasma bubbles and related parameters were computed.
The diurnal analysis from a single station (HKWS) showed that the EPB occurrence rate, depth, and duration were highly correlated with the ROTI. The large correlation coefficient between the diurnal variation in the ROTI and that of the occurrence rate (R = 0.85), depth (R = 0.85), and duration (R = 0.72) indicated that TEC depletions and ionospheric irregularities are usually closely associated.
The monthly EPB occurrence rate, depth, duration, and ROTI over the region were derived from 2013 to 2019. They all showed strong seasonal variations, with maxima during equinoctial seasons, especially during the moderate-to-high solar activity years of 2013–2016. The correlation coefficients between these three parameters and the ROTI were 0.83, 0.94, and 0.85, respectively, on a monthly time scale. Although the correlation coefficient of the monthly variation in EPB occurrences was slightly smaller than that of the daily variation, the correlation coefficients between the monthly variations in EPB depths and durations and the ROTI were greater than those of the daily variations. Furthermore, two seasonal asymmetries could be clearly seen for these parameters during the moderate-to-high solar activity period from 2013 to 2016, although they had slightly different patterns, attributed to the difference in observation techniques.
The annual variations in EPB occurrence rate, depth, and duration were presented, which varied with solar activity exhibited by the F10.7 index. In the ascending phase of SC 24 from 2013 to 2014, a small increase was seen, and then they decreased gradually from 2014 to the low solar activity year of 2019. The correlation analysis of the EPB occurrence rate, depth, and duration were found to be more strongly correlated with the F10.7 index on an annual basis than on a monthly basis. The correlation analysis of the monthly variations for three seasonal categories showed that the impacts of solar activity on EPB occurrence, depth, and duration were seasonally dependent, which was obviously more significant in the equinoctial seasons and summer than in winter.
It should be noted that the stations are located near the region at the northern crest of the EIA, therefore the observed features in EPB and ROTI are mostly related to patterns over the specific region. Because of high variability in EPB morphology at the EIA, it is probable that distinct features of EPB may be observed in other different regions.

Author Contributions

Q.L. and Y.Z. conceived and designed the experiments; Q.L. developed the method; Q.L. and J.F. gathered the data and developed the software; K.F. validated the data analysis; Q.L. drafted the manuscript; Q.L. and K.F. reviewed and revised the manuscript; Y.Z. offered guidance and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (project no. 2017YFB0503404), the National Natural Science Foundation of China (grant no. 61871012 and U1833125), the Open Fund Project of the Intelligent Operation Key Laboratory of Civil Aviation Airport Group (grant no. KLAGIO20180405), the National Key Research and Development Program of China (grant no. 2018YFB0505105), and the Beijing Nova Program of Science and Technology (grant no. Z191100001119134).

Acknowledgments

All the observation data used in this study were made available free of charge by the International GNSS Service (IGS), the Survey and Mapping Office of the Lands Department of the Hong Kong Special Administrative Region, the World Data Center at Kyoto University in Japan (http://wdc.kugi.kyoto-u.ac.jp/), and National Oceanic and Atmospheric Administration (NOAA)′s National Centers for Environmental Information (http://www.ngdc.noaa.gov). The authors are grateful to these institutions and the station operators for their valuable service.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kelley, M.C. Equatorial spread-F: Recent results and outstanding problems. J. Atmos. Terr. Phys. 1985, 47, 745–752. [Google Scholar] [CrossRef]
  2. Kil, H. The morphology of equatorial plasma bubbles—A review. J. Astron. Space Sci. 2015, 32, 13–19. [Google Scholar] [CrossRef] [Green Version]
  3. Kintner, P.M.; Ledvina, B.M.; de Paula, E.R. GPS and ionospheric scintillations. Space Weather 2007, 5, S09003. [Google Scholar] [CrossRef]
  4. Pullen, S.; Park, Y.S.; Enge, P. Impact and mitigation of ionospheric anomalies on ground-based augmentation of GNSS. Radio Sci. 2009, 44, RS0A21. [Google Scholar] [CrossRef]
  5. Pradipta, R.; Doherty, P.H. Assessing the occurrence pattern of large ionospheric TEC gradients over the brazilian airspace. Navig. J. Inst. Navig. 2016, 60, 335–343. [Google Scholar] [CrossRef]
  6. Sultan, P.J. Linear theory and modeling of the Rayleigh-Taylor instability leading to the occurrence of equatorial spread F. J. Geophys. Res. 1996, 101, 26875–26891. [Google Scholar] [CrossRef]
  7. Dungey, J.W. Convective diffusion in the equatorial F region. J. Atmos. Terr. Phys. 1956, 9, 304–310. [Google Scholar] [CrossRef]
  8. Basu, S.; Kelley, M.C. A review of recent observations of equatorial scintillations and their relationship to current theories of F region irregularity generation. Radio Sci. 1979, 14, 471–485. [Google Scholar] [CrossRef]
  9. Fejer, B.G.; Kelley, M.C. Ionospheric irregularities. Rev. Geophys. Space Phys. 1980, 18, 401–454. [Google Scholar] [CrossRef]
  10. McNamara, L.F.; Caton, R.G.; Parris, R.T.; Pedersen, T.R.; Thompson, D.C.; Wiens, K.C.; Groves, K.M. Signatures of equatorial plasma bubbles in VHF satellite scintillations and equatorial ionograms. Radio Sci. 2013, 48, 89–101. [Google Scholar] [CrossRef]
  11. Smith, J.; Heelis, R.A. Equatorial plasma bubbles: Variations of occurrence and spatial scale in local time, longitude, season, and solar activity. J. Geophys. Res. Space Phys. 2017, 122, 5743–5755. [Google Scholar] [CrossRef]
  12. Pi, X.; Mannucci, A.J.; Lindqwister, U.J.; Ho, C.M. Monitoring of global ionospheric irregularities using the Worldwide GPS Network. Geophys. Res. Lett. 1997, 24, 2283–2286. [Google Scholar] [CrossRef]
  13. Beach, T.L.; Kintner, P.M. Simultaneous global positioning system observations of equatorial scintillations and total electron content fluctuations. J. Geophys. Res. Space Phys. 1999, 104, 22553–22565. [Google Scholar] [CrossRef]
  14. Nishioka, M.; Saito, A.; Tsugawa, T. Occurrence characteristics of plasma bubble derived from global ground-based GPS receiver networks. J. Geophys. Res. Space Phys. 2008, 113, A05301. [Google Scholar] [CrossRef]
  15. Seemala, G.K.; Valladares, C.E. Statistics of total electron content depletions observed over the South American continent for the year 2008. Radio Sci. 2011, 46, RS5019. [Google Scholar] [CrossRef]
  16. Blanch, E.; Altadill, D.; Juan, J.M.; Camps, A.; Barbosa, J.; González-Casado, G.; Riba, J.; Sanz, J.; Vazquez, G.; Orús-Pérez, R. Improved characterization and modeling of equatorial plasma depletions. J. Space Weather Space Clim. 2018, 8, A38. [Google Scholar] [CrossRef] [Green Version]
  17. Tang, L.; Chen, W.; Louis, O.-P.; Chen, M. Study on seasonal variations of plasma bubble occurrence over Hong Kong area using GNSS observations. Remote Sens. 2020, 12, 2423. [Google Scholar] [CrossRef]
  18. Tsunoda, R.T. Control of the seasonal and longitudinal occurrence of equatorial scintillations by the longitudinal gradient in integrated E region pedersen conductivity. J. Geophys. Res. 1985, 90, 447–456. [Google Scholar] [CrossRef]
  19. Kumar, S.; Chen, W.; Liu, Z.Z.; Ji, S. Effects of solar and geomagnetic activity on the occurrence of equatorial plasma bubbles over Hong Kong. J. Geophys. Res. Space Phys. 2016, 121, 9164–9178. [Google Scholar] [CrossRef]
  20. Deng, B.; Huang, J.; Kong, D.; Xu, J.; Wan, D.; Lin, G. Temporal and spatial distributions of TEC depletions with scintillations and ROTI over south China. Adv. Space Res. 2015, 55, 259–268. [Google Scholar] [CrossRef]
  21. Yang, Z.; Liu, Z.Z. Correlation between ROTI and ionospheric scintillation indices using Hong Kong low-latitude GPS data. GPS Solut. 2016, 20, 815–824. [Google Scholar] [CrossRef]
  22. Acharya, R.; Majumdar, S. Statistical relation of scintillation index S4 with ionospheric irregularity index ROTI over Indian equatorial region. Adv. Space Res. 2019, 64, 1019–1033. [Google Scholar] [CrossRef]
  23. Nguyen, V.K.; Rovira-Garcia, A.; Juan, J.M.; Sanz, J.; González-Casado, G.; Ta, T.H. Measuring phase scintillation at different frequencies with conventional GNSS receivers operating at 1 Hz. J. Geod. 2019, 93, 1985–2001. [Google Scholar] [CrossRef] [Green Version]
  24. Carrano, C.S.; Groves, K.M.; Rino, C.L. On the relationship between the rate of change of total electron content index (ROTI), irregularity strength (CkL), and the scintillation index (S4). J. Geophys. Res. Space Phys. 2019, 124, 2099–2112. [Google Scholar] [CrossRef]
  25. Ma, G.; Maruyama, T. Derivation of TEC and estimation of instrumental biases from GEONET in Japan. Ann. Geophys. 2003, 49, 121–133. [Google Scholar] [CrossRef] [Green Version]
  26. Blewitt, G. An automatic editing algorithm for GPS data. Geophys. Res. Lett. 1990, 17, 199–202. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, N.; Yuan, Y.; Li, Z.; Montenbruck, O.; Tan, B. Determination of differential code biases with multi-GNSS observations. J. Geod. 2016, 90, 209–228. [Google Scholar] [CrossRef]
  28. Coco, D.S.; Coker, C.; Dahlke, S.R.; Clynch, J.R. Variability of GPS satellite differential group delay biases. IEEE Trans. Aerosp. Electron. Syst. 1991, 27, 931–938. [Google Scholar] [CrossRef]
  29. Zhang, B.; Teunissen, P.J.; Yuan, Y.; Zhang, X.; Li, M. A modified carrier-to-code leveling method for retrieving ionospheric observables and detecting short-term temporal variability of receiver differential code biases. J. Geod. 2019, 93, 19–28. [Google Scholar] [CrossRef] [Green Version]
  30. Hernández-Pajares, M.; Juan, J.M.; Sanz, J.; Aragón-Àngel, À.; García-Rigo, A.; Salazar, D.; Escudero, M. The ionosphere: Effects, GPS modeling and the benefits for space geodetic techniques. J. Geod. 2011, 85, 887–907. [Google Scholar] [CrossRef]
  31. Datta-Barua, S.; Walter, T.; Pullen, S.; Luo, M.; Blanch, J.; Enge, P. Using WAAS ionospheric data to estimate LAAS short baseline gradients. In Proceedings of the 2002 National Technical Meeting of the Institute of Navigation, San Diego, CA, USA, 28–30 January 2002; pp. 523–530. [Google Scholar]
  32. Pradipta, R.; Valladares, C.E.; Doherty, P.H. An effective TEC data detrending method for the study of equatorial plasma bubbles and traveling ionospheric disturbances. J. Geophys. Res. Space Phys. 2015, 120, 11048–11055. [Google Scholar] [CrossRef]
  33. Abdu, M.A. Equatorial ionosphere–thermosphere system: Electrodynamics and irregularities. Adv. Space Res. 2005, 35, 771–787. [Google Scholar] [CrossRef]
Figure 1. (a) Distributions of the nine Global Navigation Satellite System (GNSS) stations in Hong Kong used in the analysis. (b) Solar radio flux at 10.7 cm (F10.7) index from 1975 to 2019. The period denoted by two black dashed lines indicates 2013–2019, which is the period studied in the paper.
Figure 1. (a) Distributions of the nine Global Navigation Satellite System (GNSS) stations in Hong Kong used in the analysis. (b) Solar radio flux at 10.7 cm (F10.7) index from 1975 to 2019. The period denoted by two black dashed lines indicates 2013–2019, which is the period studied in the paper.
Sensors 20 06200 g001
Figure 2. An example of a plasma bubble detected over station HKWS and GPS satellite PRN 19 on 9 March 2014. (a) Time series of slant total electron content (STEC) I cmc (black line), I phi (green line) and filtered trend using mechanical rotary motion-implicit terrain (MRMIT) filtering (blue line); (b) Detrended STEC time series (green line); (c) rate of TEC (ROT) (black line) and rate of TEC change index (ROTI) (blue line) time series. The red dashed-line rectangle denotes the period of STEC depletion. The horizontal axis represents the local time (LT, UT + 8).
Figure 2. An example of a plasma bubble detected over station HKWS and GPS satellite PRN 19 on 9 March 2014. (a) Time series of slant total electron content (STEC) I cmc (black line), I phi (green line) and filtered trend using mechanical rotary motion-implicit terrain (MRMIT) filtering (blue line); (b) Detrended STEC time series (green line); (c) rate of TEC (ROT) (black line) and rate of TEC change index (ROTI) (blue line) time series. The red dashed-line rectangle denotes the period of STEC depletion. The horizontal axis represents the local time (LT, UT + 8).
Sensors 20 06200 g002
Figure 3. Daily equatorial plasma bubble (EPB) occurrences over HKWS from 2013 to 2019. (a) Diurnal variation in the number of EPB events and number of GPS satellites with EPB observations. The total numbers are shown. (b) Scatterplot showing the correlation between the daily number of EPB events and number of GPS satellites with EPB observations during a geomagnetically quiet period (Dst ≥ −50 nT). The red solid line is obtained by linear regression. The correlation coefficient (R) and p-value (P) are shown.
Figure 3. Daily equatorial plasma bubble (EPB) occurrences over HKWS from 2013 to 2019. (a) Diurnal variation in the number of EPB events and number of GPS satellites with EPB observations. The total numbers are shown. (b) Scatterplot showing the correlation between the daily number of EPB events and number of GPS satellites with EPB observations during a geomagnetically quiet period (Dst ≥ −50 nT). The red solid line is obtained by linear regression. The correlation coefficient (R) and p-value (P) are shown.
Sensors 20 06200 g003
Figure 4. EPB histograms over the HKWS station from 2013 to 2019 of (a) the depletion depth in TEC units (TECU); (b) ROTI in TECU/min; (c) LT at maximum depletion; (d) depletion duration in minutes.
Figure 4. EPB histograms over the HKWS station from 2013 to 2019 of (a) the depletion depth in TEC units (TECU); (b) ROTI in TECU/min; (c) LT at maximum depletion; (d) depletion duration in minutes.
Sensors 20 06200 g004aSensors 20 06200 g004b
Figure 5. Diurnal occurrence rate over the HKWS station during the geomagnetically quiet period (Dst ≥ −50 nT) from 2013 to 2019. (a) Variations in the daily occurrence rate (red bars) and ROTI in 0.1 TECU/min (black bars). (b) Correlation between the daily occurrence rate and the ROTI in TECU/min. The red solid line is obtained via linear regression.
Figure 5. Diurnal occurrence rate over the HKWS station during the geomagnetically quiet period (Dst ≥ −50 nT) from 2013 to 2019. (a) Variations in the daily occurrence rate (red bars) and ROTI in 0.1 TECU/min (black bars). (b) Correlation between the daily occurrence rate and the ROTI in TECU/min. The red solid line is obtained via linear regression.
Sensors 20 06200 g005
Figure 6. Diurnal EPB depth in TECU over the HKWS station during the geomagnetically quiet period (Dst ≥ −50 nT) from 2013 to 2019. (a) Variations in the daily EPB depth (red bars) and ROTI in 0.1 TECU/min (black bars). (b) Correlation between the daily EPB depth and the ROTI in TECU/min. The red solid line is obtained via linear regression.
Figure 6. Diurnal EPB depth in TECU over the HKWS station during the geomagnetically quiet period (Dst ≥ −50 nT) from 2013 to 2019. (a) Variations in the daily EPB depth (red bars) and ROTI in 0.1 TECU/min (black bars). (b) Correlation between the daily EPB depth and the ROTI in TECU/min. The red solid line is obtained via linear regression.
Sensors 20 06200 g006
Figure 7. Diurnal EPB duration in minutes over the HKWS station during the geomagnetically quiet period (Dst ≥ −50 nT) from 2013 to 2019. (a) Variations in the daily EPB duration (red bars) and ROTI in 0.1 TECU/min (black bars). (b) Correlation between the daily EPB duration and the ROTI in TECU/min. The red solid line is obtained via linear regression.
Figure 7. Diurnal EPB duration in minutes over the HKWS station during the geomagnetically quiet period (Dst ≥ −50 nT) from 2013 to 2019. (a) Variations in the daily EPB duration (red bars) and ROTI in 0.1 TECU/min (black bars). (b) Correlation between the daily EPB duration and the ROTI in TECU/min. The red solid line is obtained via linear regression.
Sensors 20 06200 g007
Figure 8. Variation in the monthly EPB occurrences over nine stations from 2013 to 2019. (a) Cumulative number of days with EPB occurrences in each month; (b) monthly maximum ROTI at nighttime.
Figure 8. Variation in the monthly EPB occurrences over nine stations from 2013 to 2019. (a) Cumulative number of days with EPB occurrences in each month; (b) monthly maximum ROTI at nighttime.
Sensors 20 06200 g008
Figure 9. Monthly EPB occurrence rate over Hong Kong from 2013 to 2019. (a) Monthly occurrence rate and ROTI in 0.1 TECU/min; (b) scatterplot of the monthly EPB occurrence rate vs. the ROTI in TECU/min. The red line is obtained via linear regression.
Figure 9. Monthly EPB occurrence rate over Hong Kong from 2013 to 2019. (a) Monthly occurrence rate and ROTI in 0.1 TECU/min; (b) scatterplot of the monthly EPB occurrence rate vs. the ROTI in TECU/min. The red line is obtained via linear regression.
Sensors 20 06200 g009
Figure 10. Monthly EPB depth over Hong Kong from 2013 to 2019. (a) Monthly EPB depth and ROTI in 0.1 TECU/min; (b) scatterplot of the monthly maximum EPB depth vs. the ROTI in TECU/min. The red line is obtained via linear regression.
Figure 10. Monthly EPB depth over Hong Kong from 2013 to 2019. (a) Monthly EPB depth and ROTI in 0.1 TECU/min; (b) scatterplot of the monthly maximum EPB depth vs. the ROTI in TECU/min. The red line is obtained via linear regression.
Sensors 20 06200 g010
Figure 11. Monthly EPB duration over Hong Kong from 2013 to 2019. (a) Monthly EPB duration and ROTI in 0.1 TECU/min; (b) Scatterplot of the monthly maximum EPB duration vs. the ROTI in TECU/min. The red line is obtained via linear regression.
Figure 11. Monthly EPB duration over Hong Kong from 2013 to 2019. (a) Monthly EPB duration and ROTI in 0.1 TECU/min; (b) Scatterplot of the monthly maximum EPB duration vs. the ROTI in TECU/min. The red line is obtained via linear regression.
Sensors 20 06200 g011
Figure 12. Semiannual variation in EPBs over Hong Kong from 2013 to 2019. (a) Semiannual mean F10.7 index, EPB occurrence rate, depth, and duration; (b) scatterplot of the semiannual EPB occurrence rate vs. the mean F10.7 index; (c) scatterplot of the semiannual EPB depth vs. the mean F10.7 index; (d) scatterplot of the semiannual EPB duration vs. the mean F10.7 index. The red lines are obtained via linear regression.
Figure 12. Semiannual variation in EPBs over Hong Kong from 2013 to 2019. (a) Semiannual mean F10.7 index, EPB occurrence rate, depth, and duration; (b) scatterplot of the semiannual EPB occurrence rate vs. the mean F10.7 index; (c) scatterplot of the semiannual EPB depth vs. the mean F10.7 index; (d) scatterplot of the semiannual EPB duration vs. the mean F10.7 index. The red lines are obtained via linear regression.
Sensors 20 06200 g012
Figure 13. (a) Monthly variations in the occurrence rate and F10.7 index; (b) scatterplot showing the correlation between the monthly occurrence rate and the mean F10.7 index. The months have been grouped into summer (red circles), winter (black triangles), and equinoctial (blue squares) months. The lines are obtained via linear regression for summer (red line), winter (black line), and the equinoctial seasons (blue line).
Figure 13. (a) Monthly variations in the occurrence rate and F10.7 index; (b) scatterplot showing the correlation between the monthly occurrence rate and the mean F10.7 index. The months have been grouped into summer (red circles), winter (black triangles), and equinoctial (blue squares) months. The lines are obtained via linear regression for summer (red line), winter (black line), and the equinoctial seasons (blue line).
Sensors 20 06200 g013
Figure 14. (a) Monthly variations in the EPB depth in TECU and the F10.7 index; (b) scatterplot showing the correlation between the monthly EPB depth and the mean F10.7 index. The months have been grouped into summer (red circles), winter (black triangles), and equinoctial (blue squares) months. The lines are obtained via linear regression in summer (red line), winter (black line), and equinoctial seasons (blue line).
Figure 14. (a) Monthly variations in the EPB depth in TECU and the F10.7 index; (b) scatterplot showing the correlation between the monthly EPB depth and the mean F10.7 index. The months have been grouped into summer (red circles), winter (black triangles), and equinoctial (blue squares) months. The lines are obtained via linear regression in summer (red line), winter (black line), and equinoctial seasons (blue line).
Sensors 20 06200 g014
Figure 15. (a) Monthly variations in the EPB duration in minutes and the F10.7 index; (b) scatterplot showing the correlation between the monthly EPB duration and the mean F10.7 index. The months have been grouped into summer (red circles), winter (black triangles), and equinoctial (blue squares) months. The lines are obtained from the linear regression in summer (red line), winter (black line), and the equinoctial seasons (blue line).
Figure 15. (a) Monthly variations in the EPB duration in minutes and the F10.7 index; (b) scatterplot showing the correlation between the monthly EPB duration and the mean F10.7 index. The months have been grouped into summer (red circles), winter (black triangles), and equinoctial (blue squares) months. The lines are obtained from the linear regression in summer (red line), winter (black line), and the equinoctial seasons (blue line).
Sensors 20 06200 g015
Table 1. List of GNSS stations in Hong Kong with their geographic and geomagnetic coordinates.
Table 1. List of GNSS stations in Hong Kong with their geographic and geomagnetic coordinates.
Station CodeGeographic LongitudeGeographic LatitudeGeomagnetic Latitude
HKKT114°03′ E22°26′ N12°40′ N
HKNP113°53′ E22°14′ N12°28′ N
HKOH114°13′ E22°14′ N12°29′ N
HKPC114°02′ E22°17′ N12°31′ N
HKSC114°08′ E22°19′ N12°33′ N
HKSL113°55′ E22°22′ N12°54′ N
HKSS114°16′ E22°25′ N12°41′ N
HKST114°18′ E22°23′ N12°38′ N
HKWS114°20′ E22°26′ N12°40′ N
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Q.; Zhu, Y.; Fang, K.; Fang, J. Statistical Study of the Seasonal Variations in TEC Depletion and the ROTI during 2013–2019 over Hong Kong. Sensors 2020, 20, 6200. https://doi.org/10.3390/s20216200

AMA Style

Li Q, Zhu Y, Fang K, Fang J. Statistical Study of the Seasonal Variations in TEC Depletion and the ROTI during 2013–2019 over Hong Kong. Sensors. 2020; 20(21):6200. https://doi.org/10.3390/s20216200

Chicago/Turabian Style

Li, Qiang, Yanbo Zhu, Kun Fang, and Jisi Fang. 2020. "Statistical Study of the Seasonal Variations in TEC Depletion and the ROTI during 2013–2019 over Hong Kong" Sensors 20, no. 21: 6200. https://doi.org/10.3390/s20216200

APA Style

Li, Q., Zhu, Y., Fang, K., & Fang, J. (2020). Statistical Study of the Seasonal Variations in TEC Depletion and the ROTI during 2013–2019 over Hong Kong. Sensors, 20(21), 6200. https://doi.org/10.3390/s20216200

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

Article Metrics

Back to TopTop