太阳活动对电离层TEC变化影响分析ormalsize

为研究太阳活动对电离层TEC变化的影响,从整体到局部分析了2000—2016年的太阳黑子数、太阳射电流量F_10.7指数日均值与电离层TEC的关系,并重点分析了2017年9月6日太阳爆发X9.3级特大耀斑前后15天太阳活动与电离层TEC变化的相关性.结果表明:由2000—2016年的数据整体看来,太阳黑子数、太阳F_10.7指数、TEC两两之间具有很强的整体相关性,但局部相关性强弱不均;此次耀斑爆发前后太阳黑子数、太阳F_10.7指数和TEC具有很强的正相关特性,太阳活动对TEC的影响时延约为2天;太阳活动对全球电离层TEC的影响不同步,从高纬至低纬约有1天的延迟,且对低纬度的影响远大于中高纬度.太阳活动是影响电离层TEC变化的主要原因,但局部也可能存在其他重要影响因素.      

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

The solar eclipse and its associated ionospheric TEC variations over Indian stations on January 15, 2010

An annular solar eclipse occurred over the Indian subcontinent during the afternoon hours of January 15, 2010. This event was unique in the sense that solar activity was minimum and the eclipse period coincides with the peak ionization time at the Indian equatorial and low latitudes. The number of GPS receivers situated along the path of solar eclipse were used to investigate the response of total electron content (TEC) under the influence of this solar eclipse. These GPS receivers are part of the Indian Satellite Based Augmentation System (SBAS) named as ‘GAGAN’ (GPS Aided Geo Augmented Navigation) program. The eight GPS stations located over the wide range of longitudes allows us to differentiate between the various factors induced due to solar eclipse over the equatorial and low latitude ionosphere. The effect of the eclipse was detected in diurnal variations of TEC at all the stations along the eclipse path. The solar eclipse has altered the ionospheric behavior along its path by inducing atmospheric gravity waves, localized counter-electrojet and attenuation of solar radiation intensity. These three factors primarily control the production, loss and transport of plasma over the equatorial and low latitudes. The localized counter-electrojet had inhibited the equatorial ionization anomaly (EIA) in the longitude belt of 72°E–85°E. Thus, there was a negative deviation of the order of 20–40% at the equatorial anomaly stations lying in this ‘inhibited EIA region’. The negative deviation of only 10–20% is observed for the stations lying outside the ‘inhibited EIA region’. The pre-eclipse effect in the form of early morning enhancement of TEC associated with atmospheric gravity waves was also observed during this solar eclipse. More clear and distinctive spatial and temporal variations of TEC were detected along the individual satellite passes. It is also observed that TEC starts responding to the eclipse after 30 min from start of eclipse and the delay of the maximum TEC deviation from normal trend with respect to the maximum phase of the eclipse was close to one hour in the solar eclipse path.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

基于IGS电离层TEC格网的扰动特征统计分析

电离层总电子含量（TEC）是研究空间天气特性的重要参量,通过分析电离层TEC,可以了解空间环境的变化特征.利用IGS提供的1999—2016年全球电离层TEC格网数据,按照地磁纬度将全球划分为高、中、中低、低磁纬四个区域,计算不同区域的电离层扰动;利用大量统计数据选取电离层扰动事件的判定阈值,分析电离层扰动与太阳活动、时空之间的关系;计算电离层扰动指数与地磁活动之间的相关系数.结果显示:电离层扰动与太阳活动变化具有较强的正相关特性.在太阳活动低年,电离层扰动事件发生的概率约为1.79%,在太阳活动高年发生扰动的概率约为10.18%.在空间分布上,无论是太阳活动高年还是低年,高磁纬地区发生扰动事件的概率均大于其他磁纬出现扰动事件的概率.计算得到的中磁纬和中低磁纬地区电离层扰动指数与全球地磁指数Ap的相关系数分别为0.57和0.56,说明电离层扰动指数与Ap具有较好的相关关系;高磁纬电离层扰动指数与Ap的相关系数为0.44;低磁纬扰动指数与Ap的相关系数为0.39.以上结果表明,不同区域电离层扰动与全球地磁指数Ap的相关性不同,测定区域地磁指数可能会提高与电离层扰动的相关性.      

Effects of X2-class solar flare events on ionospheric GPS-TEC and radio waves over Brazilian sector

In this investigation, we present and discuss the effects of 6 X2-class solar flare events in the ionospheric F region over Brazilian sector that occurred during 2013 to 2015. For this investigation, we present the vertical total electron content (VTEC) observations from nearly 120 Global Positioning System (GPS) receivers all over the Brazilian sector for each event. Also, ionospheric sounding observations obtained in São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S; hereafter referred to as SJC), under the southern crest of the equatorial ionospheric anomaly (EIA), Brazil, are presented. The observations show that the greatest TEC impact occurs with the EUV fluxes increases lasting for more than one hour and when the solar active region is located close to the solar disc center. We present a detailed study of the efficiency of the EUV flux with wavelengths ranging from 0.1 to 190?nm for the F region ionization. The largest increase of ΔTEC occurs below the magnetic equator line, covering mainly the central, northeast, southeast and south regions, which includes the equatorial ionospheric anomaly (EIA) region. The ionograms show partial or total fade out in the echoes traces observed causing blackouts of radio signals of up to 60?min, which can have serious consequences to technological systems of public and private agencies around Brazilian sector. This study can help to better understand the effects of solar flares in the ionospheric F region.     

Variation of ionospheric total electron content in Indian low latitude region of the equatorial anomaly during May 2007–April 2008

The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.     

TEC variability near northern EIA crest and comparison with IRI model

Monthly median values of hourly total electron content (TEC) is obtained with GPS at a station near northern anomaly crest, Rajkot (geog. 22.29°N, 70.74°E; geomag. 14.21°N, 144.9°E) to study the variability of low latitude ionospheric behavior during low solar activity period (April 2005 to March 2006). The TEC exhibit characteristic features like day-to-day variability, semiannual anomaly and noon bite out. The observed TEC is compared with latest International Reference Ionosphere (IRI) – 2007 model using options of topside electron density, NeQuick, IRI01-corr and IRI-2001 by using both URSI and CCIR coefficients. A good agreement of observed and predicted TEC is found during the daytime with underestimation at other times. The predicted TEC by NeQuick and IRI01-corr is closer to the observed TEC during the daytime whereas during nighttime and morning hours, IRI-2001 shows lesser discrepancy in all seasons by both URSI and CCIR coefficients.     

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

On the uncertainties in the measurement of absolute (true) TEC over Indian equatorial and low latitude sectors

The Indian sector encompasses the equatorial and low latitude regions where the ionosphere is highly dynamic and is characterized by the equatorial ionization anomaly (EIA) resulting in large latitudinal electron density gradients causing errors and uncertainties in the estimation of range delays in satellite based navigation systems. The diurnal and seasonal variations of standard deviations in the TEC data measured during the low sunspot period 2004–2005 at 10 different Indian stations located from equator to the anomaly crest region and beyond are examined and presented. The day-to-day variability in TEC is found to be lowest at the equatorial station and increases with latitude up to the crest region of EIA and decreases beyond.     

Statistical study of the time delay of ionospheric TEC storms to geomagnetic storms in Taoyuan,Taiwan

We have studied the time delay of ionospheric storms to geomagnetic storms at a low latitude station Taoyuan (25.02°N, 121.21°E), Taiwan using the Dst and TEC data during 126 geomagnetic storms from the year 2002 to 2014. In addition to the known local time dependence of the time delay, the statistics show that the time delay has significant seasonal characteristics, which can be explained within the framework of the seasonal characteristics of the ionospheric TEC. The data also show that there is no correlation between the time delay and the intensity of magnetic storms. As for the solar activity dependence of the time delay, the results show that there is no relationship between the time delay of positive storms and the solar activity, whereas the time delay of negative storms has weakly negative dependence on the solar activity, with correlation coefficient −0.41. Especially, there are two kinds of extreme events: pre-storm response events and long-time delay events. All of the pre-storm response events occurred during 15–20 LT, manifesting the Equator Ionospheric Anomaly (EIA) feature at Taoyuan. Moreover, the common features of the pre-storm response events suggest the storm sudden commencement (SSC) and weak geomagnetic disturbance before the main phase onset (MPO) of magnetic storms are two main possible causes of the pre-storm response events. By analyzing the geomagnetic indices during the events with long-time delay, we infer that this kind of events may not be caused by magnetic storms, and they might belong to ionospheric Q-disturbances.     

2004年11月一次磁暴期间全球电离层TEC扰动分析

利用全球分布的GPS原始观测数据提取的电离层总电子含量(TEC)分析了2004年11月6日至12日期间全球电离层暴的形态特点与发展过程.结果表明,11月8日磁暴主相期间电离层暴以大范围的强烈正暴为主,在11月10日的恢复相,Dst又一次降到最低值前后期间,电离层再次受到很强的扰动,大范围的正暴和负暴交替出现.这次磁暴期间夏季半球的负暴更加强烈,反映出负暴偏向于在夏季半球发生的季节变化特点.另外,磁暴期间,夜晚TEC值普遍比磁暴前的平静期要低,具体是什么机制导致还需要进一步收集数据和分析.      

Empirical model of the TEC response to the geomagnetic activity over the North American region

The paper presents an empirical model of the total electron content (TEC) response to the geomagnetic activity described by the K_p-index. The model is built on the basis of TEC measurements covering the region of North America (50°W–150°W, 10°N–60°N) for the period of time between October 2004 and December 2009. By using a 2D (latitude-time) cross-correlation analysis it is found that the ionospheric response to the geomagnetic activity over the considered geographic region and at low solar activity revealed both positive and negative phases of response. The both phases of the ionospheric response have different duration and time delay with respect to the geomagnetic storm. It was found that these two parameters of the ionospheric response depend on the season and geographical latitude. The presence of two phases, positive and negative, of the ionospheric response imposed the implementation of two different time delay constants in order to properly describe the two different delayed reactions. The seasonal dependence of the TEC response to geomagnetic storms is characterized by predominantly positive response in winter with a short (usually ∼5–6 h) time delay as well as mainly negative response in summer with a long (larger than 15 h) time delay. While the TEC response in March and October is more close to the winter one the response in April and September is similar to the summer one.     

Comparison of observed ionospheric vertical TEC over the sea in Indian region with IRI-2016 model

The vertical ionospheric TEC values obtained from GAGAN grid based ionospheric delay correction values over the sea in the Indian equatorial region have been compared with the corresponding values derived from the International Reference Ionosphere model, IRI-2016. The objective of this work is to study the deviation of the vertical TEC derived from the IRI model from ground truths over the sea for different conditions. This will serve the basic intention of assessing the candidature of the IRI model as an alternative ionospheric correction model in navigation receivers in terms of accuracy. We have chosen different solar activity periods, seasons, geomagnetic conditions, locations etc. for our comparison and analysis. The TEC values by the IRI-2016 were compared with the actual measured values for the given conditions and errors were obtained. The measured vertical TEC values at the ionospheric grid points were derived from the GAGAN broadcast ionospheric delay data and used as reference. The IRI model with standard internal functions was used in estimating the TEC at the same ionospheric grid points. The errors in the model derived values are statistically analysed. Broadly, the results show that, for the Indian sector over the sea, the IRI model performs better on quiet days in off equatorial regions, particularly in the northern region. The overall performance degrades for other conditions with the model generally underestimating the true TEC values and most severely in the equatorial region. The performance is worst in this region for the disturbed days of the equinoctial period. The comparison study is also done with the TEC data measured directly by dual frequency GPS receivers. The results were found to be in general agreement with those obtained by comparing the model with GAGAN broadcast data as reference. This study will be useful in considering the IRI-2016 model for real time estimates of TEC as an alternative to the current parametric model in a satellite navigation receiver in absence of other options.     

基于因果卷积与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预报模型的实际应用提供了参考。      

2009年7月22日日全食期间电离层总电子含量变化研究

利用中国区域内五个GPS台站(一个台站处于日全食区域、四个台站处于日偏食区域)观测数据, 研究2009年7月22日日全食期间电离层总电子含量(TEC)的变化, 结果表明, 日全食期间, 电离层TEC值经历了下降和恢复的过程, 最小TEC相对于最大食偏的时间延迟约为1～10min; 台站测得最小TEC的星下点(IPP)越靠近日全食带TEC下降量越大, 在日食期间武汉站(114.35°E, 30.53°N) TEC相对于各参考日期的TEC, 其平均下降量最大, 达到4.58TECU.      

2000年7月和2003年10月大磁暴期间东亚地区中低纬电离层的GPS TEC的响应研究

利用武汉电离层观象台研制的GPS TEC的现报方法及现报系统,对东亚地区GPS台网的观测数据进行处理分析,特别对2000年7月14-18日和2003年10月28日至11月1日两次特大磁暴期间的数据进行了对比考察,文中分析了两次磁暴间的电离层响应,得到对应不同磁暴时段电离层TEC的不同变化情况,着重揭示了TEC赤道异常峰的压缩和移动以及赤道异常随时间的压缩—反弹—恢复的过程,并结合高纬电离层的部分响应机制进行了说明,结果显示,两次磁暴期的电离层响应表现出了各自不同的特点,从而反映出因季节变化引起的高纬电离层暴时能量注入的不同而造成的全球性电离层扰动的不同形态,由此看出,磁暴期间电离层TEC的变化直接与太阳扰动发生的时间及其对高纬电离层的耦合有关,若短时期内连续发生多次磁暴,则电离层反应更加复杂,不能简单地当做单一磁暴叠加处理。     

Analysis of ionospheric TEC response to solar and geomagnetic activities at different solar activity stages

Studying the relationship of total electron content (TEC) to solar or geomagnetic activities at different solar activity stages can provide a reference for ionospheric modeling and prediction. On the basis of solar activity indices, geomagnetic activity parameters, and ionospheric TEC data at different solar activity stages, this study analyzes the overall variation relationships of solar and geomagnetic activities with ionospheric TEC, the characteristics of the quasi-27-day periodic oscillations of the three variables at different stages, and the delayed TEC response of solar activity by conducting correlation analysis, Butterworth band-pass filtering, Fourier transform, and time lag analysis. The following results are obtained. (1) TEC exhibits a significant linear relationship with solar activity at different solar activity stages. The correlation coefficients |R| are arranged as follows: |R|_EUV > |R|_F10.7 > |R|sunspot number. No significant linear relationship exists between TEC and geomagnetic activity parameters (|R| < 0.35). (2) TEC, solar activity indices, and geomagnetic activity parameters have a period of 10.5 years. The maximum amplitudes of the Fourier spectrum for TEC and solar activity indices are nearly 27 days and those of geomagnetic activity parameters are nearly 27 and 13.5 days. (3) The deviations of the quasi-27-day significant periodic oscillation of TEC and solar activity indices are consistent. (4) No evident relationship exists between the quasi-27-day periodic oscillation of TEC and geomagnetic activity parameters. (5) The delay time of TEC for the 10.7 cm solar radio flux and extreme ultraviolet is always consistent, whereas that for sunspot number varies at each stage.     

Validate the IRI2007 model by the COSMIC slant TEC data during the extremely solar minimum of 2008

During 2008, the solar activity is extremely low. The satellite observations show that the ionospheric height and electron density is much lower than the predictions by the international reference ionosphere (IRI) model. In this paper, we compared the slant total electron content (TEC) observed by the COSMIC satellites during 2008 with the IRI model results. It is found that the IRI model with IRI2001 and IRI2001 Cor. topside options will always overestimate the electron density in both lower and higher altitudes. But the rest two topside options (NeQuick, and TTS) tend to overestimate the electron density in the F layer and underestimate it in the topside altitudes. The switch altitude between overestimation and underestimation and the latitude-local time distribution of the model deviation depend on the topside option. The current investigation might be useful for the model improvement as well as data assimilation work based on the IRI model and the LEO TEC data.