Dependence of the ionospheric response on the solar flare parameters based on the theoretical modeling and GPS data

The measurements of an increase in the total electron content (TEC) of the ionosphere during solar flares, obtained based on the GPS data, indicated that up to 30% of TEC increments corresponded to the ionospheric regions above 300 km altitude in some cases, and TEC increased mainly below altitudes of 300 km in other cases. The theoretical model of the ionosphere and plasmasphere was used to study the obtained effects. The altitude-time variations in the charged particle density in the ionospheric region from 100 to 1000 km were used depending on the solar flare spectrum. An analysis of the modeling results indicated that an intensification of the flare UV emission in the 55–65 and 85–95 nm spectral ranges results in a pronounced increase in the electron density in the topside ionosphere (above 300 km). The experimental dependences of the ionospheric TEC response amplitude on the localization and peak power of flares on the Sun in the X-ray range, obtained based on the GPS data, are also presented in the work.     

顶部电离层和等离子体层电子密度分布——基于GRACE星载GPS信标测量的CT反演

本文利用两颗跟飞的GRACE卫星载GPS信标测量数据和基于差分相对TEC的层析算法,实现了全球范围的顶部电离层和等离子体层(450~5000 km) 层析成像.反演结果表明,利用低轨道卫星载GPS信标测量数据可以有效地重建顶部电离层和等离子体层的全球二维分布图像.对不同地磁活动条件下的天基层析反演结果表明,等离子体层电子密度随纬度的分布是不均匀的;在低纬赤道带,从顶部电离层向上延伸直到等离子体层,以及等离子体层中局地的电离增强云团,经常出现近似垂直于磁力线的电子密度柱状增强结构.     

Disturbances in the topside ionosphere during solar flares

The results of studying the ionospheric response to solar flares, obtained from the data of the GPS signal observations and incoherent scatter radars and as a result of the model calculations, are presented. It is shown that, according to the GPS data, a flare can cause a decrease in the electron content at altitudes of the topside ionosphere (h > 300 km). Similar effects of formation of a negative disturbance in the ionospheric F region were also observed during the solar flares of May 21 and 23, 1967, with the Arecibo incoherent scatter radar. The mechanism by which negative disturbances appear in the topside ionosphere during solar flares has been studied in this work based on the theoretical model of the ionosphere-plasmasphere coupling. It has been indicated that the formation of the electron density negative disturbance in the topside ionosphere is caused by an intense removal of O⁺ ions into the overlying plasmasphere under the action of an abrupt increase in the ion production rate and thermal expansion of the ionospheric plasma.     

Comparison of ionospheric equivalent slab thickness with bottomside digisonde profile over Wuhan

A comparison of the diurnal and seasonal variations in the ionospheric equivalent slab thickness (τ) and bottomside slab thickness (B0) is presented based on the observation during high solar activities at a mid-latitude station—Wuhan (114.4°E, 30.6°N). The investigated data include foF2, hmF2, B0, B1, and TEC, and are derived from the measured ionogram and GPS receiver over Wuhan from April 1999 to March 2000. The results show that τ and B0 are highly/weakly correlated during the day/night, respectively. Furthermore, a comprehensive discussion of the relation between τ, B0, and hmF2 for geomagnetic storm events is provided in this paper.     

Variations in ionospheric parameters during solar flares

Results of studying the ionospheric response to solar flares, obtained based on the incoherent scatter radar observations of the GPS signals and as a result of the model simulations, are presented. The method, based on the effect of partial “shadowing” of the atmosphere by the globe, has been used to analyze the GPS data. This method made it possible to estimate the value of a change in the electron content in the upper ionosphere during the solar flare of July 14, 2000. It has been shown that a flare can cause a decrease in the electron content at heights of the upper ionosphere (h > 300 km) according to the GPS data. Similar effects in the formation of a negative disturbance in the ionospheric F region were also observed during the solar flares of May 21 and 23, 1967, at the Arecibo incoherent scatter radar. The mechanism by which negative disturbances are formed in the upper ionosphere during solar flares has been studied based on the theoretical model of the ionosphere-plasmasphere coupling. It has been shown that an intense ejection of O⁺ ions into the above located plasmasphere under the action of a sharp increase in the ion production rate and the thermal expansion of the ionospheric plasma cause the formation of a negative disturbance in the electron concentration in the upper ionosphere.     

Large-scale disturbances of auroral origin during strong magnetic storms of October 29–31, 2003, and November 7–11, 2004, according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered according to measurements of the total electron content (TEC) during the magnetic storms of October 29–31, 2003, and November 7–11, 2004, has been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region (foF2) were used for this purpose. The variations in TEC and foF2 were similar for all events mentioned above. The previous assumption that the region of thickness 150–200 km in the vicinity of the ionospheric F region mainly contributes to TEC modulation was confirmed for the cases when the electron density disturbance at an F region maximum was not more than 50%. However, this region probably becomes more extensive in vertical when the electron density disturbance in the vicinity of the ionospheric F region is about 85%.     

NeQuick bottomside analysis at low latitudes

NeQuick ionospheric electron density model produces the full electron density profile in the ionosphere using the F2 layer peak values (foF2 and hmF2) as anchor points. Each part of the profile is modeled using Epstein layer formalism. Simple empirical relations are used to compute the thicknesses of each semi-Epstein layer. It has been observed that when NeQuick model is used to estimate total electron content at low latitudes the modeled values tend to underestimate the observed ones. Beside the F2 peak values, the most important profile parameter is the thickness of the F2 layer bottomside. The present study focuses on NeQuick model behavior at low latitudes comparing modeled profiles parameters with the ones extracted from experimental data mostly from African and Indian sector at different levels of solar activity and different time of the day. Possible model improvements are discussed.     

Asymmetric response of the topside ionosphere to large-scale IGW generated during the November 30, 1979, substorm

We used bottomside ground observations and topside sounding data from the Intercosmos-19 satellite to study a Travelling Ionospheric Disturbance (TID) that occurred in response to Large-Scale Internal Gravity Wave (LSIGW) propagation during a substorm on November 30, 1979. We built a global scheme for the wavelike ionospheric variations during this medium substorm (AE_max ～800 nT). The area where the TID was observed looks like a wedge since it covers the nighttime hours at subauroral latitudes but contracts to a ～02 h local sector at low latitudes. The ionospheric response is strongly asymmetric because the wedge area and the TID amplitude are larger in the winter hemisphere than in the summer hemisphere. Clear evidence was obtained indicating that the more powerful TID from the Northern (winter) hemisphere propagated across the equator into the low latitude Southern (summer) hemisphere. Intercosmos-19 observations show that the disturbance covers the entire thickness of the topside ionosphere, from h_mF2 up to at least the 1000 km satellite altitude at post-midnight local times. F-layer lifting reached ～200 km, N_e increases in the topside ionosphere by up to a factor of ～1.9 and variations in N_mF2 of both signs were observed. Assumptions are made concerning the reason for the IGW effect at high altitudes in the topside ionosphere. The relationship between TID parameters and source characteristics determined from a global network of magnetometers are studied. The role of the dayside cusp in the generation of the TID in the daytime ionosphere is discussed. The magnetospheric electric field effects are distinguished from IGW effects.     

应用经验正交函数估算顶部电离层电子密度剖面

本文基于IRI模型、地面数字测高仪和GNSS TEC数据,提出了一种利用经验正交函数(Empirical Orthogonal Function,简称EOF)估算顶部电离层电子密度剖面的方法,并将其应用于美国Millstone Hill测高仪和GNSS数据以估算顶部电离层电子密度剖面.通过将估算的临界频率、峰值高度、400km以上电子密度分别与测高仪实测临界频率、测高仪实测峰值高度以及非相干散射雷达实测400km以上电子密度作对比以对方法的有效性进行验证.统计结果显示估算临界频率、峰值高度与测高仪实测数据基本一致,400km以上估算电子密度相较于非相干散射雷达实测的绝对误差平均值仅是测高仪推算400km以上电子密度绝对误差平均值的一半左右.所以本文提出的方法可以更加精确地估算顶部电离层电子密度.     

Modeling the responses of the middle latitude ionosphere to solar flares

In this paper, we investigate the solar flare effects of the ionosphere at middle latitude with a one-dimensional ionosphere theoretical model. The measurements of solar irradiance from the SOHO/Solar EUV Monitor (SEM) and GOES satellites have been used to construct a simple time-dependent solar flare spectrum model, which serves as the irradiance spectrum during solar flares. The model calculations show that the ionospheric responses to solar flares are largely related to the solar zenith angle. During the daytime most of the relative increases in electron density occur at an altitude lower than 300 km, with a peak at about 115 km, whereas around sunrise and sunset the strongest ionospheric responses occur at much higher altitudes (e.g. 210 km for a summer flare). The ionospheric responses to flares in equinox and winter show an obvious asymmetry to local midday with a relative increase in total electron content (TEC) in the morning larger than that in the afternoon. The flare-induced TEC enhancement increases slowly around sunrise and reaches a peak at about 60 min after the flare onset.     

基于非相干散射雷达和GPS观测研究Millstone Hill地区等离子体层电子含量

本文尝试结合非相干散射雷达和GPS TEC观测数据提取等离子体层总电子含量(PTEC).我们首先描述所用的技术方法,然后具体利用了Millstone Hill台站的观测数据研究该地区上空等离子体层总电子含量(PTEC)的变化情况.我们采用变化标高的Chapman函数对非相干散射雷达测得的电子浓度剖面数据进行拟合,然后通过对剖面积分得到100 km到1000 km高度范围的电离层总电子含量.GPS提供的TEC数据为高度达20200 km的总电子含量,两者之差可近似看成等离子体层的电子含量.本文分别选取太阳活动高年(2000, 2002年)和太阳活动低年(2005,2008年)Millstone Hill台站的静日数据进行研究.结果表明,等离子体层电子含量及其所占GPS TEC的比例具有明显的周日变化.PTEC含量在白天高于夜间,而所占GPS TEC的百分比,夜间明显高于白天.太阳活动高年所选月份等离子体层电子含量在4~14 TECU (1TECU=10¹⁶el/m²) 范围内变化,夜间所占比例可达60%左右.太阳活动低年所选月份等离子体层电子含量在3~7 TECU范围内变化,所占比例夜间最高可达80%左右.我们所得到的结果与前人基于其它观测手段所得结果在变化趋势上一致,在量级上也大致相当.因此,这从一个侧面证明了我们所用方法的可靠性.非相干散射雷达能够探测包括F2层峰值以下及以上高度的电子浓度,利用这一设备所观测得到的资料来推算电离层电子含量将比前人基于电离层垂测仪观测资料进行的推算更具真实性,由此得到的等离子体层电子含量也将更为接近真实情况.     

Response of the ionosphere to the Baltic Sea earthquake of 21 September 2004

GPS observations of the European permanent network were used to identify seismo-ionospheric precursors of Baltic Sea earthquake of 21 September 2004. It is a very rare event for this region of Europe (magnitude of about 5.0). This value is the threshold for the occurrence of seismic effects in the ionosphere. In total electron content (TEC) data over the region of the earthquake, a specific ionospheric anomaly appeared one day before the earthquake was detected. The ionospheric variability had a positive sign with an enhancement of about 4–5 TECU (1 TECU = 10¹⁶ electrons/m²) relative to the non-disturbed state of the ionosphere. The anomaly had a duration of 4–5 hours in the day time. The special size of this anomaly was about 1000 km. The characteristic parameters of the anomaly show that it can be associated with ionospheric precursors of an earthquake.     

Prompt derivation of TEC from GEONET data for space weather monitoring in Japan

Continuous monitoring of ionospheric conditions is essential to monitoring and forecasting space weather. The worldwide use of global navigation satellite systems like the Gobal Positioning System (GPS) makes it possible to continuously monitor the total electron content (TEC) of the ionosphere and plasmasphere up to a height of about 20,000 km. We have developed a system for deriving the TEC from GEONET data rapidly and we use the TEC distribution over Japan in the daily operations of the Space Weather Forecast Center at NICT (RWC Tokyo of ISES). Using instrumental biases from a few days before enables us to drastically shorten the processing time for deriving TEC. The latest TEC values (with a delay of about 1 h) are obtained every 3 h, and most of the values are within 2 TEC units of the actual TEC. We have found our system for deriving TEC rapidly to be useful for continuously monitoring the progress of ionospheric storms under any ionospheric conditions, even those under which the usual ionosonde observations are unable to obtain F-region profiles.     

Ionospheric slab thickness and its seasonal variations observed by GPS

The ionospheric slab thickness, the ratio of the total electron content (TEC) to the F2-layer peak electron density (NmF2), is closely related to the shape of the ionospheric electron density profile Ne (h) and the TEC. Therefore, the ionospheric slab thickness is a significant parameter representative of the ionosphere. In this paper, the continuous GPS observations in South Korea are firstly used to study the equivalent slab thickness (EST) and its seasonal variability. The averaged diurnal medians of December–January–February (DJF), March–April–May (MAM), June–July–August (JJA) and September–October–November (SON) in 2003 have been considered to represent the winter, spring, summer and autumn seasons, respectively. The results show that the systematic diurnal changes of TEC, NmF2 and EST significantly appeared in each season and the higher values of TEC and NmF2 are observed during the equinoxes (semiannual anomaly) as well as in the mid-daytime of each season. The EST is significantly smaller in winter than in summer, but with a consistent variation pattern. During 14–16 LT in daytime, the larger EST values are observed in spring and autumn, while the smaller ones are in summer and winter. The peaks of EST diurnal variation are around 10–18 LT which are probably caused by the action of the thermospheric wind and the plasmapheric flow into the F2-region.     

Seasonal variation of the total electron content, maximum electron density and equivalent slab thickness at a South-American station

Total electron content (TEC) and foF2 ionosonde data obtained at Tucumán (26.9°S; 65.4°W) from April 1982 to March 1983 (high solar activity period) are analyzed to show the seasonal variation of TEC, NmF2 (proportional to square of foF2) and the equivalent slab thickness EST. Bimonthly averages of the monthly median for January–February, April–May, July–August and October–November have been considered to represent summer, autumn, winter and spring seasons, respectively. The results show that the higher values of TEC and maximum electron density of F2-layer NmF2 are observed during the equinoxes (semiannual anomaly). During daytime, both in TEC and in NmF2 the seasonal or winter anomaly can be seen. At nighttime, this effect is not observed. Also, the observed NmF2 values are used to check the validity of International Reference Ionosphere (IRI) to predict the seasonal variability of this parameter. In general, it is found that averaged monthly medians (obtained with the IRI model) overestimate averaged monthly median data for some hours of the day and underestimate for the other hours.     

Comparison of the measured and modelled electron densities and temperatures in the ionosphere and plasmasphere during 20/30 January, 1993

We present a comparison of the electron density and temperature behaviour in the ionosphere and plasmasphere measured by the Millstone Hill incoherent-scatter radar and the instruments on board of the EXOS-D satellite with numerical model calculations from a time-dependent mathematical model of the Earths ionosphere and plasmasphere during the geomagnetically quiet and storm period on 20/30 January, 1993. We have evaluated the value of the additional heating rate that should be added to the normal photoelectron heating in the electron energy equation in the daytime plasmasphere region above 5000 km along the magnetic field line to explain the high electron temperature measured by the instruments on board of the EXOS-D satellite within the Millstone Hill magnetic field flux tube in the Northern Hemisphere. The additional heating brings the measured and modelled electron temperatures into agreement in the plasmasphere and into very large disagreement in the ionosphere if the classical electron heat flux along magnetic field line is used in the model. A new approach, based on a new effective electron thermal conductivity coefficient along the magnetic field line, is presented to model the electron temperature in the ionosphere and plasmasphere. This new approach leads to a heat flux which is less than that given by the classical Spitzer-Harm theory. The evaluated additional heating of electrons in the plasmasphere and the decrease of the thermal conductivity in the topside ionosphere and the greater part of the plasmasphere found for the first time here allow the model to accurately reproduce the electron temperatures observed by the instruments on board the EXOS-D satellite in the plasmasphere and the Millstone Hill incoherent-scatter radar in the ionosphere. The effects of the daytime additional plasmaspheric heating of electrons on the electron temperature and density are small at the F-region altitudes if the modified electron heat flux is used. The deviations from the Boltzmann distribution for the first five vibrational levels of N₂(v) and O₂(v) were calculated. The present study suggests that these deviations are not significant at the first vibrational levels of N₂ and O₂ and the second level of O₂, and the calculated distributions of N₂(v) and O₂(v) are highly non-Boltzmann at vibrational levels v > 2. The resulting effect of N₂(v > 0) and O₂(v > 0) on NmF2 is the decrease of the calculated daytime NmF2 up to a factor of 1.5. The modelled electron temperature is very sensitive to the electron density, and this decrease in electron density results in the increase of the calculated daytime electron temperature up to about 580 K at the F2 peak altitude giving closer agreement between the measured and modelled electron temperatures. Both the daytime and night-time densities are not reproduced by the model without N₂(v > 0) and O₂(v > 0), and inclusion of vibrationally excited N₂ and O₂ brings the model and data into better agreement.     

Ionospheric behavior over Europe during the solar eclipse of 3 October 2005

An annular eclipse occurred over Europe in the morning hours of 3 October 2005. The well-defined obscuration function of the solar radiation during the eclipse provided a good opportunity to study the ionospheric/thermospheric response to solar radiation changes. Since the peak electron density behavior of the ionospheric F2 layer follows the local balance of plasma production, loss and transport, the ionospheric plasma redistribution processes significantly affect the shape of the electron density profile. These processes are discussed here based on a comparison of vertical incidence sounding (VS) and vertical total electron content (TEC) data above-selected ionosonde stations in Europe. The equivalent slab thickness, derived with a time resolution of 10 min, provides relatively good information on the variation of the electron density profile during the eclipse. The computations reveal an increased width of the ionosphere around the maximum phase. As indicated by the available measurements over Spain, the photo production is significantly reduced during the event leading to a slower increase of the total ionization in comparison with the neighboring days. The supersonic motion of the Moon's cool shadow through the atmosphere may generate atmospheric gravity waves that propagate upward and are detectable as traveling ionospheric disturbances at ionospheric heights. High-frequency (HF) Doppler shift spectrograms were recorded during the eclipse showing a distinct disturbance along the eclipse path. Whereas the ionosonde measurements at the Ebro station/Spain in the vicinity of the eclipse path reveal the origin of the wave activity in the lower thermosphere below about 180 km altitude, the similar observations at Pruhonice/Czech Republic provide arguments to localize the origin of the abnormal waves in the middle atmosphere well below the ionospheric heights. Although ionosonde and HF Doppler measurements show enhanced wave activity, the TEC data analysis does not, which is an indication that the wave amplitudes are too small for detecting them via this interpolation method. The total ionization reduces up to about 30% during the event. A comparison with similar observations from the solar eclipse of 11 August 1999 revealed a quite different ionospheric behavior at different latitudes, a fact that needs further investigation.     

Altitude variation of the plasmapause signature in the main ionospheric trough

The projection of the plasmapause magnetic-field lines to low altitudes, where the light-ion chemistry is dominated by O⁺, tends to occur near the minimum electron density in the main (midlatitude) electron density trough at night. With increasing altitude in the trough, where H⁺ emerges as the dominant ion on the low-latitude boundary, we have found cases where the plasmapause field lines are located on the sharp low-latitude side of the trough as expected if this topside ionosphere H⁺ distribution varies in step with the plasmapause gradient in the distant plasmasphere. These conclusions are based on near-equatorial crossings of the plasmapause (corresponding to the steep gradient in the dominant species H⁺) by the Explorer-45 satellite as determined from electric-field measurements by Maynard and Cauffman in the early 1970s and ISIS-2 ionospheric topside-sounder measurements. The former data have now been converted to digital form and made available at http://nssdcftp.gsfc.nasa.gov. The latter provide samples of nearly coincident observations of ionospheric main trough crossings near the same magnetic-field lines of the Explorer 45-determined equatorial plasmapause. The ISIS-2 vertical electron density profiles are used to infer where the F-region transitions from an O⁺ to a H⁺ dominated plasma through the main trough boundaries.     

Solar and latitude dependence of TEC and SLAB thickness

Several years of the vertical total electron content (TEC) and of the equivalent slab thickness data measured at different European locations are analysed by using the linear regression technique to demonstrate the annual response of these two ionospheric parameters to changes in solar activity and to their dependence on latitude. It is shown that TEC is highly correlated with the latitude and solar activity, whereas the slab thickness is independent of both solar activity and latitude. Their intercepts relative to solar activity show clearly a seasonal dependence.