Ionospheric response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

Mapping seasonal trends of electron temperature in the topside ionosphere based on DEMETER data

The diurnal, seasonal and latitudinal variations of the electron temperature in the Earth‘s topside ionosphere during relatively low solar activity period of 2005 – 2008 are investigated. In order to examine seasonal variations and morphology of the topside ionospheric plasma temperature, CNES micro-satellite DEMETER ISL data are used. Presented study is oriented on the dataset gathered in 2005 and 2008. Within conducted analysis, global maps of electron temperature for months of equinoxes and solstices have been developed. Furthermore, simultaneous studies on two-dimensional time series based on DEMETER measurements and predictions obtained with the IRI-2012 model supply examination of the topside ionosphere during recent deep solar minimum. Comparison with the IRI-2012 model reveals discrepancies between data and prediction, that are especially prominent during the periods of very low solar activity.     

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.     

A study on Nighttime Winter Anomaly (NWA) and other related Mid-latitude Summer Nighttime Anomaly (MSNA) in the light of International Reference Ionosphere (IRI) – Model

The ionospheric Nighttime Winter Anomaly (NWA) is a feature observed in the Northern Hemisphere at the American and in the Southern Hemisphere at the Asian longitude sector under low solar activity conditions. Jakowski et al. (2015) analyzed ground-based GPS derived TEC and peak electron density data from radio occultation measurements on Formosat-3/COSMIC satellites and confirmed the persistence of the phenomenon. Further, they assumed that Mid-latitude Summer Nighttime Anomaly (MSNA) and related special anomalies such as the Weddell Sea Anomaly (WSA) and the Okhotsk Sea Anomaly (OSA) are closely related to the NWA via enhanced wind-induced uplifting of the ionosphere. The aim of this paper is to study the factors causing these anomalies and also to investigate if these anomalies are re-produced by IRI. The results show that IRI model does include the NWA effect, though at a different longitude and could be improved for better predictions. The IRI-2016 model does show WSA in TEC but not in NmF2. Further, the IRI-2016 model could clearly predict the OSA both in NmF2 and TEC.     

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.     

国际参考电离层IRI－90与中国实测电离层月平均电子密度剖面形状的比较

将太阳活动峰年期间中国4个电离层站垂直探测得到的月平均电子浓度剖面与国际参考电离层IKI-90进行了系统的比较。结果表明国际参考电离层所计算的峰下电子含量(或峰下半厚)总体来说偏大。一天中白天符合较好,晚上较差,对中纬台站较好,对低纬台站较差。      

大耀斑期间向日面电离层总电子含量的响应个例分析

利用2001年4月15日1336UT耀斑爆发期间向日面GPS观测数据提取的总电子含量的时间变化曲线。分析了向日面电离层对这次耀斑的响应特点．结果表明，耀斑期间向日面电离层出现了总电子含量突增事件．最大总电子含量增加量约为2．6TECU，在0600LT和1800LT都观测到了总电子含量突增，世增加幅度仅为0．5-1TECU．在高纬地区，由于电离层闪烁，从TEC时间变化曲线提取不出来总电子含量增加值．从各卫星星下点处的TEC增加量和各星下点处的太阳天顶角的关系可以看到，TEC增加量与太阳天顶角有关，太阳天顶角越大，TEC增幅越小。另外，从总电子含量时间变化率曲线上还观测到了时间同步的小尺度扰动，通过与耀斑期间硬X射线辐射通量的比较，发现两者有明显的相关性，电离层中的这种扰动与耀斑期间的硬X射线或远紫外辐射有关．     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

基于双频GPS观测的电离层TEC与硬件延迟反演方法

利用全球定位系统（Global Positioning System,GPS）的双频观测数据反演得到电离层的总电子含量（Total Electron Content,TEC）,使得广域甚至全球范围高时空分辨率的电离层观测研究成为可能,但由于GPS卫星和接收机对信号的硬件延迟可导致TEC测量系统偏差,因此,需要探索反演TEC并估测GPS卫星与接收机硬件延迟的有效算法.本文根据电离层电波传播理论,阐述了基于双频GPS观测提取电离层TEC的方法,给出TEC与硬件延迟的基本关系.综合研究了TEC与硬件延迟的反演方法,进行分析与归纳分类,在此基础上提出了有待深入研究的问题.      

Cross-hemisphere comparison of mid-latitude ionospheric variability during 1996–2009: Juliusruh vs. Hobart

Analysis of a long-time series of hourly median characteristics of the ionospheric plasma at two mid-latitude locations in the Northern and Southern hemisphere, Juliusruh (54.6N; 13.4E) and Hobart (42.9S; 147.3E), reveals patterns of their synchronous and independent variability. We studied timelines of GPS vTEC, ionogram-derived F2-layer peak electron density NmF2, ionospheric equivalent slab thickness τ, and their ratios at two locations during the complete 23rd solar cycle and its following period of the extremely low solar activity in 2008–2009. This study has also involved the comparative analysis of the observed data versus the model predictions by IRI-2012. During the high solar activity in 2000–2002, seasonal variations show a complicated cross-hemisphere behavior influenced by the winter and semi-annual anomalies, with the largest noon-time values of TEC and NmF2 observed around equinoxes. Strength of the winter anomaly in NmF2 was significantly greater at Juliusruh in comparison with Hobart. The winter anomaly in GPS vTEC values was much weaker than in NmF2 for the Northern hemisphere mid-latitudes and was entirely absent at the Southern hemisphere. Cross-hemisphere analysis of the equivalent slab thickness shows its clear seasonal dependence for all levels of solar activity: the day-time maximum τ_max is observed during local summer, whereas the day-time minimum τ_min is observed during local winter. The night-time values of τ were higher compared to the day-time values during the winter and equinox seasons. Comparative model-data study shows rather good IRI performance of the day-time NmF2 for mid-latitudes of both hemispheres and rather noticeable overestimations for the mid-night NmF2 values during high solar activity. Analysis of IRI vTEC demonstrates the model limitations, related with the absence of the plasmaspheric part, and actual demand in a reliable and standard ionosphere–plasmasphere model for analysis of GPS vTEC.     

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.     

Prediction of total electron content using the International Reference Ionosphere

The International Reference Ionosphere (IRI) is a model of the ionosphere, based on experimental data, which has been proposed as a standard ionospheric model. As such, it should be tested extensively to determine its range of validity. One of the ways in which the electron denisty profile given by the IRI, especially above the peak of the F layer, can be tested is to compare calculated and observed values of total electron content (TEC). We have therefore studied the discrepancies between calculated and observed values of TEC recorded at 15 stations covering a wide range of longitudes and latitudes, mainly in the northern hemisphere, and mainly for high levels of solar activity. W have found that the IRI produces reasonably accurate values of TEC at mid and high latitudes, but that it greatly underestimates the daytime values of TEC at low latitudes. We conclude therefore that the daytime electron density profile given by the IRI is reasonably accurate at mid and high latitudes, at least above the peak of the F2 layer. The situation at low latitudes clearly requires more work, and we have suggested two possible lines of study. The generally low discrepancies at night indicate that the night-time electron density profiles given by the IRI correspond fairly closely to the actual profiles.     

Ionospheric effects of the solar flares as deduced from global GPS network data

Results derived from analysing the ionospheric response to faint and bright solar flares are presented. The analysis used novel technology of a global detection of ionospheric effects from solar flares as developed by the authors (Afraimovich, 2000a; Afraimovich, 2000b), on the basis of phase measurements of the total electron content (TEC) in the ionosphere using an international GPS network. The essence of the method is that use is made of appropriate filtering and a coherent processing of variations in the TEC which is determined from GPS data, simultaneously for the entire set of visible GPS satellites at all stations used in the analysis. This technique is useful for identifying the ionospheric response to faint solar flares (of X-ray class C) when the variation amplitude of the TEC response to separate line-on-sight (LOS) is comparable to the level of background fluctuations. The dependence of the TEC variation response amplitude on the flare location on the Sun is investigated.     

Evaluating the IRI topside model for the South African region: An overview of the modelling techniques

The representation of the topside ionosphere (the region above the F2 peak) is critical because of the limited experimental data available. Over the years, a wide range of models have been developed in an effort to represent the behaviour and the shape of the electron density (Ne) profile of the topside ionosphere. Various studies have been centred around calculating the vertical scale height (VSH) and have included (a) obtaining VSH from Global Positioning System (GPS) derived total electron content (TEC), (b) calculating the VSH from ground-based ionosonde measurements, (c) using topside sounder vertical Ne profiles to obtain the VSH. One or a combination of the topside profilers (Chapman function, exponential function, sech-squared (Epstein) function, and/or parabolic function) is then used to reconstruct the topside Ne profile. The different approaches and the modelling techniques are discussed with a view to identifying the most adequate approach to apply to the South African region’s topside modelling efforts. The IRI-2001 topside model is evaluated based on how well it reproduces measured topside profiles over the South African region. This study is a first step in the process of developing a South African topside ionosphere model.     

Application of neural networks to South African GPS TEC modelling

The propagation of radio signals in the Earth’s atmosphere is dominantly affected by the ionosphere due to its dispersive nature. Global Positioning System (GPS) data provides relevant information that leads to the derivation of total electron content (TEC) which can be considered as the ionosphere’s measure of ionisation. This paper presents part of a feasibility study for the development of a Neural Network (NN) based model for the prediction of South African GPS derived TEC. The South African GPS receiver network is operated and maintained by the Chief Directorate Surveys and Mapping (CDSM) in Cape Town, South Africa. Vertical total electron content (VTEC) was calculated for four GPS receiver stations using the Adjusted Spherical Harmonic (ASHA) model. Factors that influence TEC were then identified and used to derive input parameters for the NN. The well established factors used are seasonal variation, diurnal variation, solar activity and magnetic activity. Comparison of diurnal predicted TEC values from both the NN model and the International Reference Ionosphere (IRI-2001) with GPS TEC revealed that the IRI provides more accurate predictions than the NN model during the spring equinoxes. However, on average the NN model predicts GPS TEC more accurately than the IRI model over the GPS locations considered within South Africa.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

GPS scintillation and TEC gradients at equatorial latitudes in April 2006

We use observations of ionospheric scintillation at equatorial latitudes from two GPS receivers specially modified for recording, at a sampling rate of 50 Hz, the phase and the amplitude of the L1 signal and the Total Electron Content (TEC) from L1 and L2. The receivers, called GISTM (GPS Ionospheric Scintillation and TEC Monitor), are located in Vietnam (Hue, 16.4°N, 107.6°E; Hoc Mon, 10.9°N, 106.6°E). These experimental observations are analysed together with the tomographic reconstruction of the ionosphere produced by the Multi-Instrument Data Analysis System (MIDAS) for investigating the moderate geomagnetic storm which occurred on early April 2006, under low solar activity. The synergic adoption of the ionospheric imaging and of the GISTM measurements supports the identification of the scale-sizes of the ionospheric irregularities causing scintillations and helps the interpretation of the physical mechanisms generating or inhibiting the appearance of the equatorial F layer irregularities. In particular, our study attributes to the turning of the IMF (Interplanetary Magnetic Field) between northward and southward direction an important role in the inhibition of the generation of spread F irregularities resulting in a lack of scintillation enhancement in the post-sunset hours.     

On the performance of IRI-2016 to predict the North-South asymmetry of the equatorial ionization anomaly around 73°E longitude

The ionospheric total electron content (TEC) in both northern and southern Equatorial anomaly regions are examined by using the Global Positioning System (GPS) based TEC measurements around 73°E Longitude in the Asian sector. The TEC contour charts obtained at SURAT (21.16°N; 72.78°E; 12.9°N Geomagnetic Lat.) and DGAR (7.27°S; 72.37°E; 15.3°S Geomagnetic Lat.) over 73°E longitude during a very low solar activity phase (2009) and a moderate solar activity (2012) phase are used in this study. The results show the existence of hemispheric asymmetry and the effects of solar activity on the EIA crest in occurrence time, location and strength. The results are also compared with the TEC derived by IRI-2016 Model and it is found that the North-South asymmetry at the EIA region is clearly depicted by IRI-2016 with some discrepancies (up to 20% in the northern hemisphere at SURAT and up to 40% in the southern hemisphere at DGAR station for June Solstice and up to 10% both for SURAT and DGAR for December Solstice). This discrepancy in the IRI-2016 model is found larger during the year 2012 than that during the solar minimum year 2009 at both the hemispheres. Further, an asymmetry index, (A_i) is determined to illustrate the North-South asymmetry observed in TEC at EIA crest. The seasonal, annual and solar flux dependence of this index are investigated during both solstices and compared with the TEC derived by IRI.     

基于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的相关性不同,测定区域地磁指数可能会提高与电离层扰动的相关性.