磁暴期间中国中低纬电离层不规则体与扰动分析

2017年9月8日发生了一次强磁暴，Kp指数最大值达到8.利用区域电离层格网模型（Regional Ionosphere Map，RIM）和区域ROTI（Rate of TEC Index）地图，分析了磁暴期间中国及其周边地区电离层TEC扰动特征和低纬地区电离层不规则体的产生与发展情况，同时利用不同纬度IGS（International GNSS Service）测站BJFS（39.6°N，115.9°E），JFNG（30.5°N，114.5°E）和HKWS（22.4°N，114.3°E）的GPS双频观测值，获取各测站的ROTI和DROT（Standard Deviation of Differential ROT）指数变化趋势.结果表明:此次磁暴发生期间电离层扰动先以正相扰动为主，主要发生在中低纬区域，dTEC（differential TEC）最大值达到14.9TECU，随后电离层正相扰动逐渐衰减，在低纬区域发生电离层负相扰动，dTEC最小值达到-7.2TECU；在12:30UT-13:30UT时段，中国南部低纬地区发生明显的电离层不规则体事件；相比BJFS和JFNG两个测站，位于低纬的HKWS测站的ROTI和DROT指数变化更为剧烈，这表明电离层不规则体结构存在纬度差异.      

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.     

Diagnostics of equatorial and low latitude ionosphere by TEC mapping over Brazil

The total electron content (TEC) in the equatorial and low-latitude ionosphere over Brazil was monitored in two dimensions by using 2011 data from the ground-based global navigation satellite system (GNSS) receiver network operated by the Brazilian Institute for Geography and Statistics. It was possible to monitor the spatial and temporal variations in TEC over Brazil continuously during both day and night with a temporal interval of 10 min and a spatial resolution of about 400 km. The daytime equatorial ionization anomaly (EIA) and post-sunset plasma enhancement (PS-EIA) were monitored over an area corresponding to a longitudinal extension of 4000 km in South America. Considerable day-to-day variation was observed in EIA and PS-EIA. A large latitudinal and longitudinal gradient of TEC indicated a significant ionospheric range error in application of the GNSS positioning system. Large-scale plasma bubbles after sunset were also mapped over a wide range. Depletions with longitudinally separated by more than 800 km were observed. They were extended by more than 2000 km along the magnetic field lines and drifted eastward. It is expected that 2-dimensional TEC mapping can serve as a useful tool for diagnosing ionospheric weather, such as temporal and spatial variation in the equatorial plasma trough and crest, and particularly for monitoring the dynamics of plasma bubbles.     

基于陆态网的区域电离层TEC空间变化特性分析

为研究中国陆态网区域电离层TEC在空间小尺度、高分辨率情况下的变化特性及适用精度范围，利用陆态网260个GNSS连续运行观测站数据，解算并生成2016-2017年731天陆态网区域电离层RIM格网，并进行精度验证.在同一RIM格网中，分别在经度和纬度方向上对间隔不同经纬度的TEC格网点作差分析.结果表明:陆态网区域内经度方向上TEC最大变化率和平均变化率分别为0.30TECU·（°）-1和0.11TECU·（°）-1；经度间隔1°时，TEC差值小于2TECU，且随着经度间隔的增大，其TEC差值也随之增大，并表现出一定的半年和周年变化规律；纬度方向上TEC最大变化率和平均变化率分别为1.7TECU·（°）-1和0.46TECU·（°）-1；陆态网区域内电离层TEC随纬度减小而增大，纬度间隔1°时，99.4%的TEC差值小于4TECU，且随着纬度间隔的增大，其TEC差值也随之增大，并表现出一定的半年和周年变化规律；间隔相同情况下，纬度方向上TEC的变化比经度方向上大.      

Solar and magnetic control on night-time enhancement in TEC near the crest of the Equatorial Ionization Anomaly

The ionospheric Total Electron Content (TECs), derived by dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Bhopal (23.2°N, 77.4°E, Geomagnetic 14.2°N) were analyzed for the period of January, 2005 to February, 2008. The work deals with monthly, diurnal, solar and magnetic activity variations on night-time enhancement in TEC. From a total of 157 night-time enhancements, 75 occur during pre-midnight and 82 post-midnight hours. The occurrence of night-time enhancement in TEC is utmost during summer months, followed by equinox and winter months. The occurrence of night-time enhancement in TEC decreases with increase in solar and magnetic activities. We observed that peak size and half amplitude duration are positively correlated, while time of occurrence of night-time enhancement in TEC and time of peak enhancement are negatively correlated with solar activity. The peak size, half amplitude duration, time of peak enhancement and time of occurrence of night-time enhancement in TEC shows negative correlation with magnetic activity. The results have been compared with the earlier ones and discussed in terms of possible source mechanism responsible for the enhancement at anomaly crest region.     

Ionospheric anomalies observed over South Korea preceding the Great Tohoku earthquake of 2011

We investigated the ionospheric anomalies observed before the Tohoku earthquake, which occurred near the northeast coast of Honshu, Japan on 11 March, 2011. Based on data from a ground-based Global Positioning System (GPS) network on the Korean Peninsula, ionospheric anomalies were detected in the total electron content (TEC) during the daytime a few days before earthquake. Ionospheric TEC anomalies appeared on 5, 8 and 11 March. In particular, the ionospheric disturbances on 8 March evidenced a remarkable increase in TEC. The GPS TEC variation associated with the Tohoku earthquake was an increase of approximately 20 total electron content units (TECU), observed simultaneously in local and global TEC measurements. To investigate these pre-earthquake ionospheric anomalies, space weather conditions such as the solar activity index (F10.7) and geomagnetic activity indices (the Kp and Dst indices) were examined. We also created two-dimensional TEC maps to visual the spatial variations in the ionospheric anomalies preceding the earthquake.     

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

Investigation of the characteristics of wavelike oscillations of post-sunset equatorial ionospheric irregularity by decomposing fluctuating TEC

Post-sunset ionospheric irregularities are common features of the equatorial ionosphere that affect radio communication and navigation systems; their triggering physical mechanism is not yet fully understood. Atmospheric gravity wave is considered as a seeding mechanism for the occurrence of ionospheric irregularities (Abdu et al., 2009). To understand the effects of atmospheric waves, characteristics of wavelike oscillation from ionospheric total electron content (TEC) fluctuation that can be obtained from superposition of different oscillation modes have been investigated. Decomposing fluctuating TEC into different oscillation modes and investigating oscillation characteristics of each component is also important to get insight about the characteristics of individual atmospheric waves that may cause TEC fluctuation. In this paper we have investigated characteristics of components of fluctuating TEC obtained from SCINDA GPS receiver installed at Bahir Dar, (geographic coordinate, 11.5°N, 37.6° E, and dip latitude of 2.5°N) Ethiopia during April 2012. First Empirical Mode Decomposition (EMD) has been applied to decompose TEC fluctuation into different oscillation modes that are known as Intrinsic Mode Function (IMF). Hilbert-Huang Transform (HHT) and Continuous Wavelet Transform (CWT) have been applied to investigate the characteristics of wave-like oscillations. Applying EMD on fluctuating vTEC corresponding to a GPS satellite, five components are found. Results from HHT and CWT have shown excellent agreement. In addition, it is found out that the median periods of oscillation of those five components are 9, 17, 47, 78, and 118 min. Of these periods, 17 and 47 min respectively are oscillation periods of components of TEC fluctuation with occurrence frequency of 92% and 91% that may be interpreted as the manifestation of two frequently occurring components of atmospheric gravity waves that are likely generated by the motion of solar terminator.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

Demonstration of the use of the Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS) measurements to validate GPS ionospheric imaging

There is a lack of independent ionospheric data that can be used to validate GPS imaging results at mid latitudes over severe storm times. Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS), a global network of dual-frequency ground to satellite observations, provides this missing data and here is employed as verification to show the accuracy of the ionospheric GPS images in terms of the total electron content (TEC). In this paper, the large-scale ionospheric structures that appeared during the strong geomagnetic storm of 20 November 2003 are reconstructed with a GPS tomographic algorithm, known as MIDAS, and validated with DORIS TEC measurements. The main trough shown in an extreme equatorward position in the ionospheric imaging over mainland Europe is confirmed by DORIS satellite measurements. Throughout the disturbed day, the variations of relative slant TECs between DORIS data and MIDAS results agree quite well, with the average of the mean differences about 2 TECu. We conclude that as a valuable supplement to GPS data, DORIS ionospheric measurements can be used to analyse TEC variations with a relatively high resolution, ∼10 s in time and tens of kilometres in space. This will be very helpful for identification of some highly dynamic structures in the ionosphere found at mid-latitudes, such as the main trough, TID (Travelling Ionospheric Disturbances) and SED (Storm Enhanced Density), and could be used as a valuable auxiliary data source in ionospheric imaging.     

On the estimation of regional covariance functions of TEC variations over Canada

Spatial and temporal variations of Total Electron Content (TEC) can affect GNSS high accuracy positioning. Enhanced estimation of ionospheric variations and their de-correlation can benefit differential and point positioning rapid solutions. Global and regional TEC maps can provide the overall state of ionopsheric variations in space and time domains within their accuracy limits. In this paper, these maps are exploited to retrieve ionospheric variations by means of variograms and their associated covariance functions of TEC residuals over Canadian region during geomagnetically quiet and disturbed conditions. A number of theoretical variogram functions are reviewed for modeling covariance of TEC residuals. The variogram modeling of residuals during a strong geomagnetic storm revealed variances of one order of magnitude larger compared to a rather quiet condition. Variogram models are also used in regional and local kriging interpolation experiments and their performances are evaluated. Global maps of TEC RMS by International GNSS Service and two of its analysis centres are also compared over the Canadian region during a two-year period. Realistic representation of regional variances using estimated variograms when compared to global ionospheric RMS maps are also presented.     

Observations of TEC depletions over Indian low latitudes

Results pertaining to the latitudinal extent of the ionospheric irregularities in terms of TEC depletions have been presented for the two nights namely; 28 October 2004 and 7 February 2005. This study has been carried out using the GPS–TEC over the Indian low latitude stations, at Udaipur, Hyderabad and Bengaluru. This is probably the first report of simultaneous GPS observation of TEC depletions over different latitudes from the Indian sector. The results show that the amplitude of TEC depletions due to the equatorial spread F may vary with time and the location of the observation. The maximum amplitude of the TEC depletion has been found to be about 30 TECU over Hyderabad. The depletions in TEC are found to be field aligned.     

一种基于GNSS数据的全球电离层不规则结构活动表征方法

针对如何利用GNSS(Global Navigation Satellite System)数据进行电离层扰动监测的问题,提出了一种基于GNSS数据表征全球电离层扰动的方法.利用大约400个GNSS地面站点的观测数据,计算总电子含量(Total Electron Content,TEC)变化率的标准差——ROTI(Ra...     

A nonlinear background removal method for seismo-ionospheric anomaly analysis under a complex solar activity scenario: A case study of the M9.0 Tohoku earthquake

A precise determination of ionospheric total electron content (TEC) anomaly variations that are likely associated with large earthquakes as observed by global positioning system (GPS) requires the elimination of the ionospheric effect from irregular solar electromagnetic radiation. In particular, revealing the seismo-ionospheric anomalies when earthquakes occurred during periods of high solar activity is of utmost importance. To overcome this constraint, a multiresolution time series processing technique based on wavelet transform applicable to global ionosphere map (GIM) TEC data was used to remove the nonlinear effect from solar radiation for the earthquake that struck Tohoku, Japan, on 11 March, 2011. As a result, it was found that the extracted TEC have a good correlation with the measured solar extreme ultraviolet flux in 26–34 nm (EUV_26–34) and the 10.7 cm solar radio flux (F_10.7). After removing the influence of solar radiation origin in GIM TEC, the analysis results show that the TEC around the forthcoming epicenter and its conjugate were significantly enhanced in the afternoon period of 8 March 2011, 3 days before the earthquake. The spatial distributions of the TEC anomalous and extreme enhancements indicate that the earthquake preparation process had brought with a TEC anomaly area of size approximately 1650 and 5700 km in the latitudinal and longitudinal directions, respectively.     

Statistical characteristics of locally generated ESF during equinoctial months over Sanya

Understanding the local generation rate of equatorial spread-F (ESF) is important for forecasting ionospheric scintillation. Using the GPS ionospheric scintillation/TEC and VHF radar data during March-April and September-October from 2010 to 2014, the occurrence of ionospheric scintillation, TEC fast fluctuation, and backscatter plume were studied. Through analyzing the simultaneous occurrence of ionospheric scintillation, TEC fast fluctuation and backscatter plume, the local generation rate of ESF over Sanya was investigated. The results show that the monthly generation rate varies between 0% and 68%. A significant equinoctial asymmetry of local generation rate of ESF can be found in 2010, 2013 and 2014. The local generation rate of ESF increases from 2010 to 2014 during March-April, while it does not have similar trend during September-October. The plasma vertical drift influenced by solar activity has a significant impact on the monthly generation rate. The equinoctial asymmetry of plasma vertical drift may contribute a lot to the equinoctial asymmetry of the generation rate of ESF.     

A study of L band scintillations during the initial phase of rising solar activity at an Indian low latitude station

The ionospheric scintillation and TEC (Total Electron Content) variations are studied using GPS (Global Positioning System) measurements at an Indian low latitude station Surat (21.16°N, 72.78°E; Geomagnetic: 12.90°N, 147.35°E), situated near the northern crest of the equatorial anomaly region. The results are presented for data collected during the initial phase of current rising solar activity (low to moderate solar activity) period between January 2009 and December 2011. The results show that within a total number of 656 night-time scintillation events, 340 events are observed with TEC depletions, Rate of change of TEC (ROT) fluctuations and enhancement of Rate of change of TEC Index (ROTI). A comparison of night-time scintillation events from the considered period reveal strong correlation amongst the duration of scintillation activity in S₄ index_, TEC depletion, ROT fluctuations and ROTI enhancement in the year 2011, followed by the year 2010 and least in 2009. The statistical analyses of scintillation activity with enhancement of ROTI also show that about 70–96% scintillation activity took place in equinox and winter months. Moreover, from a nocturnal variation in occurrence of scintillation with (S₄ ? 0.2) and enhancement of ROTI with (ROTI ? 0.5), a general trend of higher occurrence in pre-midnight hours of equinox and winter seasons is observed in both indices during the year 2011 and 2010, while no significant trend is observed in the year 2009. The results suggest the presence of F-region ionospheric irregularities with scale sizes of few kilometers and few hundred meters over Surat and are found to be influenced by solar and magnetic activity.     

基于因果卷积与LSTM网络的电离层总电子含量预报

电离层总电子含量（TEC）不仅是分析电离层形态的关键参数之一,同时为导航及定位等空间应用系统消除电离层附加时延提供重要支撑。由于电离层TEC的时空变化特征,本文融合因果卷积和长短时记忆网络,以太阳活动指数F_10.7、地磁活动指数Dst和电离层TEC历史数据作为特征输入,构建深度学习模型,实现提前24 h预报电离层TEC。进一步利用2005－2013年连续9年的CODE TEC数据,全面评估了模型在北京站（40°N,115°E）、武汉站（30.53°N,114.36°E）和海口站（20.02°N,110.38°E）的预报性能。结果显示不同太阳活动条件下三个站的TEC值与真实测量值的相关系数都大于0.87,均方根误差大都集中在0～1 TECU以内,且模型预报精度与纬度、太阳、地磁活动程度、季节变化相关。与仅由长短时记忆网络构成的预报模型相比,本实验模型均方根误差降低了15%,为电离层TEC预报模型的实际应用提供了参考。      

Ionospheric anomalies of local earthquakes detected by GPS TEC measurements using data from Tashkent and Kitab stations

This paper reports the ionospheric anomalies observed during strong local earthquakes (M?5.0) which occurred mostly in and around Uzbekistan in seismically active zones, during years 2006 to 2009 within approximately 1000 km distance from the observing GPS stations located in Tashkent and Kitab, Uzbekistan. The solar and geomagnetic conditions were quiet during occurrence of the selected strong earthquakes. We produce Total Electron Content (TEC) time series over both sites and apply them to detect anomalous TEC signals preceding or accompanying the local earthquakes. The results show anomalous increase or decrease of TEC before or during the earthquakes. In general the anomalies occurred 1–7 days before the earthquakes as ionospheric electromagnetic precursors. To identify the anomalous values of TEC we calculated differential TEC (dTEC). dTEC is obtained by subtracting monthly averaged diurnal vTEC from the values of observed vTEC at each epoch. This procedure removes normal diurnal variations of vTEC. The present results are in good agreement with the previous observations on ionospheric earthquake precursors reported by various researchers.     

Relative amplitude of the total electron content variations depending on geomagnetic activity

We developed a method of estimation of a relative amplitude dI/I of the total electron content (TEC) variations in the ionosphere as deduced from the data of the global GPS receivers network. To obtain statistically significant results we picked out three latitudinal belts provided in the Internet by the maximum number of GPS sites. They are high-latitudinal belt (50–80°N, 200–300°E; 59 sites), mid latitude belt (20–50°N, 200–300°E; 817 sites), and equatorial belt (±20°N, 0–360°E; 76 sites). The results of the analysis of the diurnal and latitudinal dependencies of dI/I and dI/I distribution probability for 52 days with different levels of geomagnetic activity are presented. It was found that on average the relative amplitude of the TEC variations varies within the range 0–10% proportionally to the value of the K_p geomagnetic index. In quiet conditions the relative amplitude dI/I of the TEC variations at night significantly exceeds the daytime relative amplitude. At high levels of magnetic field disturbances, the geomagnetic control of the amplitude of TEC variations at high and middle latitudes is much more significant than the regular diurnal variations. At the equatorial belt, on average, the amplitude of TEC variations in quiet and disturbed periods almost does not differ. The obtained results may be useful for development of the theory of ionospheric irregularities.