不同强度磁暴期间热层大气密度响应特征分析

利用CHAMP卫星数据,对2002-2008年12个不同强度磁暴事件期间的热层大气密度变化特征进行分析,并研究对应磁暴期间大气模式NRLMSISE-00分布特征.结果表明,大磁暴期间日侧大气密度峰值从高纬到低纬的时间延迟为2h,中小磁暴期间的延迟时间为3～4h;春秋季暴时大气密度分布基本呈南北对称分布,而夏冬季大气密度的分布是夏半球大于冬半球,春秋季暴时大气密度大于夏冬季;NRLMSISE-00大气模式得到的热层大气密度很好的体现了半球分布以及季节分布的特征,但模式模拟结果偏小;Dst指数峰值比ap指数峰值更能反应大气密度的变化情况.      

Study of TEC changes during geomagnetic storms occurred near the crest of the equatorial ionospheric ionization anomaly in the Indian sector

GPS satellites data obtained at Bhopal (23.16° N, 77.36° E, geomagnetic latitude 14.23° N) India were analyzed to study the TEC changes during several geomagnetic storms (−300 nT < Dst < −50 nT) occurred in 2005–2007. We had segregated the storms according to the Dst value, i.e. moderate storms (−100 nT < Dst ? −50 nT), strong storms (−150 nT < Dst < −100 nT), and severe storms (Dst less than −150 nT). Total of 21 geomagnetic storms (10 moderate, 9 strong, 2 severe) are considered for the present study. Deviation in vertical total electron content (VTEC) during the main phase of the storm was found to be associated with the prompt penetration of electric field originated due to the under-shielding and over-shielding conditions for almost all geomagnetic storms discussed in this paper. For most of the storms VTEC shows the positive percentage deviation during the main phase while it shows positive as well as the negative deviation during the recovery phase of the storms. The −80% deviation in VTEC was found for geomagnetic storm occurred on July 17, 2005 and the negative trend continued for recovery phase of the storm. This was mainly due to the thermospheric composition changes by Joule heating effect at auroral latitudes that generate electric field disturbance at low latitudes. Traveling ionospheric disturbances (TIDs) were responsible for the formation of wave like nature in VTEC for the storms occurred on May 15, 2005, whereas it was not observed for storm occurred on August 24, 2005.     

Investigation of ionospheric response to two moderate geomagnetic storms using GPS–TEC measurements in the South American and African sectors during the ascending phase of solar cycle 24

The responses of the ionospheric F region using GPS–TEC measurements during two moderate geomagnetic storms at equatorial, low-, and mid-latitude regions over the South American and African sectors in May 2010, during the ascending phase of solar cycle 24, are investigated. The first moderate geomagnetic storm studied reached a minimum Dst value of −64 nT at 1500 UT on 02 May 2010 and the second moderate geomagnetic storm reached a minimum Dst value of −85 nT at 1400 UT on 29 May 2010. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from Global Positioning System (GPS) observations from the equatorial to mid-latitude regions in the South American and African sectors. Our results obtained during these two moderate geomagnetic storms from both sectors show significant positive ionospheric storms during daytime hours at the equatorial, low-, and mid-latitude regions during the main and recovery phases of the storms. The thermospheric wind circulation change towards the equator is a strong indicator that suggests an important mechanism is responsible for these positive phases at these regions. A pre-storm event that was observed in the African sector from low- to the mid-latitude regions on 01 May 2010 was absent in the South American sector. This study also showed that there was no generation or suppression of ionospheric irregularities by storm events. Therefore, knowledge about the suppression and generation of ionospheric irregularities during moderate geomagnetic storms is still unclear.     

Comparison of the observed and modeled low- to mid-latitude thermosphere response to magnetic activity: Effects of solar cycle and disturbance time delay

The performance of JB2008 and NRLMSISE-00 models, in describing the response of the thermosphere to magnetic activity are evaluated against total mass density retrieved from accelerometer measurements made onboard CHAMP satellite during 5 years. We show that the global low- to mid-latitude disturbance amplitude is correctly described by the JB2008 model for low solar activity conditions and by both the JB2008 and the NRLMSISE-00 models for high solar activity conditions. For low solar activity conditions, statistics based on almost 3 years of data confirm the large underestimation by the NRLMSISE-00 model quantified by Lathuillère et al. (2008) for the year 2004. We also found that the time delay between low- to mid-latitude global thermosphere disturbance and magnetic activity is statistically well estimated by the NRLMSISE-00 and JB2008 models for disturbed conditions. For moderately disturbed conditions however, the time delay estimated by the JB2008 model is too large by about 3 h. For very disturbed conditions, we found different time delays during day-time and night-time, using new geomagnetic proxies with a 30-min time resolution.     

An empirical relation to correct storm-time thermospheric mass density modeled by NRLMSISE-00 with CHAMP satellite air drag data

With the help of STAR (Spatial Triaxial Accelerometer for Research) accelerometer measurements on board CHAMP (Challenging Minisatellite Payload), the global distributions of total mass density changes at about 400 km height during major magnetic storms are studied, aiming to improve the capability of current thermospheric model for predicting the storm-time mass density distribution. The density calculated by the NRLMSISE-00 model without using the geomagnetic active index as input is taken as a reference on top of which the storm-time changes are added. In total 19 storm events during 2001–2004 are used to perform a comprehensive statistical analysis. A relative calibration of drag coefficient along with accelerometer calibration parameters is made by fitting the CHAMP observed initial mass densities in with the NRLMSISE-00 model on quiet days before each storm. The dependences of the storm-time changes in mass density on both the total global Joule heating power, ∑_Qj$\sum Q_{\text{j}}$ and the high-resolution ring current index, Sym-H, are investigated. The lag times of mass density changes with respect to the Joule heating and Sym-H variation are obtained as a function of latitude and sunlight. By using a multiple linear regression analysis with proper time shift, an empirical relation connecting storm-time changes in mass density for 400 km height with the two parameters, ∑_Qj$\sum Q_{\text{j}}$ and Sym-H, has been worked out for different latitude and sunlight conditions (day-side or night-side). Adding a correction calculated from the empirical relation to the NRLMSISE-00 model reference leads to a better prediction of storm-time thermospheric mass density distribution.     

Storm-time atmospheric density modeling using neural networks and its application in orbit propagation

Upper atmospheric densities during geomagnetic storms are usually poorly estimated due to a lack of clear understanding of coupling mechanisms between the thermosphere and magnetosphere. Consequently, the orbit determination and propagation for low-Earth-orbit objects during geomagnetic storms have large uncertainties. Artificial neural networks are often used to identify nonlinear systems in the absence of rigorous theory. In the present study, an attempt has been made to model the storm-time atmospheric density using neural networks. Considering the debate over the representative of geomagnetic storm effect, i.e. the geomagnetic indices ap and Dst, three neural network models (NNM) are developed with ap, Dst and a combination of ap and Dst respectively. The density data used for training the NNMs are derived from the measurements of the satellites CHAMP and GRACE. The NNMs are evaluated by looking at: (a) the mean residuals and the standard deviations with respect to the density data that are not used in training process, and (b) the accuracy of reconstructing the orbits of selected objects during storms employing each model. This empirical modeling technique and the comparisons with the models NRLMSIS-00 and Jacchia-Bowman 2008 reveal (1) the capability of neural networks to model the relationship between solar and geomagnetic activities, and density variations; and (2) the merits and demerits of ap and Dst when it comes to characterizing density variations during storms.     

磁暴期间全球TEC扰动特性分析

磁暴期间白天电离层总电子含量（TEC）大幅度扰动.TEC扰动与磁暴发生时的世界时（UT）有关.利用7年的数据对TEC对磁暴的响应进行统计研究.结果显示,磁暴期间白天TEC增大明显,且在午后TEC的增大比例有一个高峰.在18:00UT-04:00UT,南美地区与其他地区相比TEC增长较大,这可能与白天的光照有关.为了研究TEC变化与磁暴的关系,结合同样时间段的Dst指数,把TEC数据分为磁暴日（Dst<-100nT）和平静日（Dst>-50nT）.研究发现,将TEC前移2h,低纬日侧地区TEC增大值随着世界时的变化与Dst变化的负相关性较好,相关系数为-0.75.在中纬度地区,将TEC扰动前移1h,相关系数为-0.61.这可能是行进式大气扰动携带着赤道向的子午风,由极区向低纬传播引起.可以认为,TEC的变化可能是由磁暴引起的.在高纬地区,TEC增大值随着世界时的变化与Dst变化的相关性较差.这可能是由于太阳高度角较低,光辐射通量较小,导致电子密度的增加不明显.      

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.     

Ionospheric response of equatorial and low latitude F-region during the intense geomagnetic storm on 24–25 August 2005

In this investigation, we present and discuss the response of the ionospheric F-region in the South American and East Asian sectors during an intense geomagnetic storm in August 2005. The geomagnetic storm studied reached a minimum Dst of −216 nT at 12:00 UT on 24 August. In this work ionospheric sounding data obtained of 24, 25, and 26 August 2005 at Palmas (PAL; 10.2° S, 48.2° W; dip latitude 6.6° S), São José dos Campos (SJC, 23.2° S, 45.9° W; dip latitude 17.6° S), Brazil, Ho Chi Minh City, (HCM; 10.5° N, 106.3° E; dip latitude 2.9° N), Vietnam, Okinawa (OKI; 26.3° N, 127.8° E; dip latitude 21.2° N), Japan, are presented. Also, the GPS observations obtained at different stations in the equatorial and low-latitude regions in the Brazilian sector are presented. On the night of 24–25 August 2005, the h′F variations show traveling ionospheric disturbances associated with Joule heating in the auroral zone from SJC to PAL. The foF2 variations show a positive storm phase on the night of 24–25 August at PAL and SJC during the recovery phase. Also, the GPS-VTEC observations at several stations in the Brazilian sector show a fairly similar positive storm phase on 24 August. During the fast decrease of Dst (between 10:00 and 11:00 UT) on 24 August, there is a prompt penetration of electric field of magnetospheric origin that result in abrupt increase (∼12:00 UT) in foF2 at PAL, SJC (Brazil) and OKI (Japan) and in VTEC at IMPZ, BOMJ, PARA and SMAR (Brazil). OKI showed strong oscillations of the F-region on the night 24 August resulted to the propagation of traveling atmospheric disturbances (TADs) by Joule heating in the auroral region. These effects result a strong positive observed at OKI station. During the daytime on 25 August, in the recovery phase, the foF2 observations showed positive ionospheric storm at HCM station. Some differences in the latitudinal response of the F-region is also observed in the South American and East Asian sectors.     

Morphology of OI 8446 Å dayglow emission

A comprehensive model is developed to study the atomic oxygen OI 8446 Å dayglow emission. The emission rate profiles and intensities are obtained using the recently developed Solar2000 EUV (Extreme Ultra Violet) flux model. These emission profiles are used to construct the morphology of the 8446 Å emission between equator and 45°N in the northern hemisphere. The longitudinal variation of 8446 Å dayglow emission is found about 5% and is not included in the presentation of morphology. A span of five years is chosen to study the effect of varying solar activity on the morphology of the OI 8446 Å dayglow emission. The morphology is studied on April 3 which lies under the equinox conditions. In year 2001 the solar F10.7 index on the chosen date was as high as 223.1 which is the case of solar maximum. From the present calculations it is found that the intensity does not vary linearly with the F10.7 solar index. The morphology shows that the region of maximum emission rate expands towards the higher latitudes as F10.7 solar index increases. The similar effects have also been found in the morphology of 7320 Å dayglow emission (Sunil Krishna and Singh, 2009). The similarities in the morphology of 7320 Å and 8446 Å dayglow emissions further suggest that the photoelectron flux has strong bearing on the production of these emissions.     

Variations in the thermosphere and ionosphere response to the 17–20 April 2002 geomagnetic storms

The responses of the thermospheric density and ionospheric foF2 to the intense magnetic storms event on 17–20 April were analyzed by using data from CHAMP/STAR and ionosonde stations respectively, and NRLMSISE-00 and IRI-2007 models were used to simulate. The models can capture the tendency of changes, especially under quiet or moderate geomagnetic conditions, but are less accurate under geomagnetic storms. The thermospheric density is sensitive to the EUV emission and geomagnetic activity, and double-peak structure appeared in the dayside. On 19 April dayside, TADs traveled toward the equator with phase speeds of the order of 300–750 m/s, interfered near the equator to produce a total density perturbation of 25%, and then passed through each other and into the opposite hemisphere. For ionospheric foF2, there are non-symmetric hemispheres’ features during the intense geomagnetic activities. In details, middle latitudes in the north and high latitudes in both hemispheres are negative ionospheric storms, and the maximum amplitudes of δfoF2$\delta \mathit{foF}2$ is about 60%, but the amplitudes decrease from the higher to lower latitudes in the Southern Hemisphere. Meanwhile, the equatorial station shows positive phase, and the maximum value is about 100%. Finally, the mechanisms for these features will be discussed in this study.     

Initial phase of mid-latitude aurora during strong geomagnetic storms

Analysing the initial mid-latitude aurora phase during strong geomagnetic storms we found that the initial phase of the mid-latitude aurorae observed at 630 nm emission during the strong geomagnetic storms on March 24, 1991, April 6, 2000, October 30 and November 20, 2003 is characterized by a short (∼1 h) wave-like disturbance. This disturbance corresponds to the beginning of main phase of the magnetic storms. The marked effect of the mid-latitude aurorae is analyzed using data on magnetosphere and ionosphere conditions in observation periods. The features of the dynamics of the 630 nm emission intensity and its connection with the dynamics of magnetospheric–ionospheric structures are considered. Possible excitation mechanisms of the atomic oxygen emission (630 nm) during these disturbances are discussed.     

海南地区强磁暴期间强区域扩展F事件研究

利用2004年海南DPS-4数字测高仪观测到的强区域扩展F(SSF)数据,分析研究了5个强磁暴(Dst＜-100 nT)事件期间海南SSF的响应特征.结果发现,在海南地区,5个强磁暴事件中有3个磁暴Dst最小值位于2200-0200 LT之间,在磁暴主相及恢复相初期均出现了SSF现象,这种触发作用可能源于磁层直接渗透电场的作用,而另两个磁暴Dst最小值均发生在白天,一个SSF现象出现在磁暴的恢复相晚间,另一个SSF现象出现在超强磁暴的初相晚间,后者可能由该超强磁暴的急始造成的直接渗透电场所触发;5个强磁暴期间发生的SSF现象或者仅出现在午夜前,或者先出现在午夜前并持续到午夜后;同时,还就这些观测结果与Dabas等人有关磁暴对ESF影响的结论进行对比和讨论.      

Feasibility of developing an ionospheric E-region electron density storm model using TIMED/SABER measurements

We present a new technique for improving ionospheric models of nighttime E-region electron densities under geomagnetic storm conditions using TIMED/SABER measurements of broadband 4.3 μm limb radiance. The response of E-region electron densities to geomagnetic activity is characterized by SABER-derived NO⁺(v) 4.3 μm Volume Emission Rates (VER). A storm-time E-region electron density correction factor is defined as the ratio of storm-enhanced NO⁺(v) VER to a quiet-time climatological average NO⁺(v) VER, which will be fit to a geomagnetic activity index in a future work. The purpose of this paper is to demonstrate the feasibility of our technique in two ways. One, we compare storm-to-quiet ratios of SABER-derived NO⁺(v) VER with storm-to-quiet ratios of electron densities measured by Incoherent Scatter Radar. Two, we demonstrate that NO⁺(v) VER can be parameterized by widely available geomagnetic activity indices. The storm-time correction derived from NO⁺(v) VER is applicable at high-latitudes.     

Moderate geomagnetic storms of January 22–25, 2012 and their influences on the wave components in ionosphere and upper stratosphere-mesosphere regions

Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF _Bz$B_z$ are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.     

Midlatitude ionospheric F2-layer response to eruptive solar events-caused geomagnetic disturbances over Hungary during the maximum of the solar cycle 24: A case study

In our study we analyze and compare the response and behavior of the ionospheric F2 and of the sporadic E-layer during three strong (i.e., Dst?<??100nT) individual geomagnetic storms from years 2012, 2013 and 2015, winter time period. The data was provided by the state-of the art digital ionosonde of the Széchenyi István Geophysical Observatory located at midlatitude, Nagycenk, Hungary (IAGA code: NCK, geomagnetic latitude: 46.17° geomagnetic longitude: 98.85°). The local time of the sudden commencement (SC) was used to characterize the type of the ionospheric storm (after Mendillo and Narvaez, 2010). This way two regular positive phase (RPP) ionospheric storms and one no-positive phase (NPP) storm have been analyzed. In all three cases a significant increase in electron density of the foF2 layer can be observed at dawn/early morning (around 6:00 UT, 07:00 LT). Also we can observe the fade-out of the ionospheric layers at night during the geomagnetically disturbed time periods. Our results suggest that the fade-out effect is not connected to the occurrence of the sporadic E-layers.     

利用人工神经网络预报大磁暴

本文采用阈值预报的策略和人工神经网络BP模型,以13个太阳风参量和地磁AE,Dst指数作为输入,以0或1作为输出,提前4h预报大磁暴主相发生的时刻．结果表明,采用神经网络方法的阈值预报可以对灾害性磁暴的发生提前数小时做出比较准确的预报．     

Evidence for the lowering of the centroid of daytime thermospheric O(D) 630.0 nm emission over the magnetic equator: First results

It is shown in this paper for the first time that the intensity of the daytime thermospheric O(¹D) 630.0 nm airglow as measured by the ground-based dayglow photometer over Trivandrum (8.5°N; 77°E; dip lat. 0.5°N), a geomagnetic dip equatorial station, exhibit a direct correlation with the electron density at 180 km. This altitude is about ∼40 km lower than the believed centroid of the O(¹D) 630.0 nm dayglow emission i.e. 220 km. This observation is contrary to the understanding of the behavior of O(¹D) 630.0 nm dayglow over equatorial/low latitudes. Over these latitudes, the variations of the measured intensity of O(¹D) 630.0 nm dayglow are known to be associated with the changes in the electron density at altitudes around 220 km, the centroid of this emission. In this context, the present results indicating the lowering of the peak altitude of O(¹D) 630.0 nm emission from ∼220 to ∼180 km over the dip equator is new. Recent results on solar XUV flux indicate that this could be an important parameter that controls the O(¹D) 630.0 nm dayglow excitation rates through modulations in the neutral and ionic composition in lower thermosphere-ionosphere region. However, the lowering of the centroid of O(¹D) 630.0 nm emission, as shown in this study, has been ascribed primarily to the fountain effect associated with the equatorial ionization anomaly.