Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model

Current theories of F-layer storms are discussed using numerical simulations with the Upper Atmosphere Model, a global self-consistent, time dependent numerical model of the thermosphere-ionosphere-plasmasphere-magnetosphere system including electrodynamical coupling effects. A case study of a moderate geomagnetic storm at low solar activity during the northern winter solstice exemplifies the complex storm phenomena. The study focuses on positive ionospheric storm effects in relation to thermospheric disturbances in general and thermospheric composition changes in particular. It investigates the dynamical effects of both neutral meridional winds and electric fields caused by the disturbance dynamo effect. The penetration of short-time electric fields of magnetospheric origin during storm intensification phases is shown for the first time in this model study. Comparisons of the calculated thermospheric composition changes with satellite observations of AE-C and ESRO-4 during storm time show a good agreement. The empirical MSISE90 model, however, is less consistent with the simulations. It does not show the equatorward propagation of the disturbances and predicts that they have a gentler latitudinal gradient. Both theoretical and experimental data reveal that although the ratio of [O]/[N₂] at high latitudes decreases significantly during the magnetic storm compared with the quiet time level, at mid to low latitudes it does not increase (at fixed altitudes) above the quiet reference level. Meanwhile, the ionospheric storm is positive there. We conclude that the positive phase of the ionospheric storm is mainly due to uplifting of ionospheric F2-region plasma at mid latitudes and its equatorward movement at low latitudes along geomagnetic field lines caused by large-scale neutral wind circulation and the passage of travelling atmospheric disturbances (TADs). The calculated zonal electric field disturbances also help to create the positive ionospheric disturbances both at middle and low latitudes. Minor contributions arise from the general density enhancement of all constituents during geomagnetic storms, which favours ion production processes above ion losses at fixed height under day-light conditions.     

Thermosphere–ionosphere response at middle and high latitudes during perturbed conditions: A case study

Neutral gas composition and ionospheric measurements taken by the Dynamic Explorer 2 satellite at F2-region heights (280–300 km) during an intense geomagnetic storm (peak Dst=−187 nT) were used to analyze the role of some possible physical mechanisms responsible for the changes of electron density at high and middle latitudes. The storm considered in this study occurred on 26 September 1982. The main features observed were increases of electron density during the initial stages of the storm at middle latitudes; followed by decreases of electron density at high and mid-high latitudes during the main phase of the storm and the first phase of the recovery. Delayed increases of electron density during the recovery phase have also been observed at mid-high latitudes (50–60°). Several mechanisms were discussed in explaining the features observed for the electron density variations.     

磁暴期间热层大气密度变化

基于CHAMP卫星资料,分析了2002—2008年267个磁暴期间400 km高度大气密度变化对季节、地方时与区域的依赖以及时延的统计学特征,得到暴时大气密度变化的一些新特点,主要结论如下: 1)两半球大气密度绝对变化(δρa)结果在不同强度磁暴、不同地方时不同.受较强的焦耳加热和背景中性风共同作用,在北半球夏季,中等磁暴过程中夜侧和大磁暴中,夏半球的δρa强于冬半球;由于夏季半球盛行风环流造成的扰动传播速度快,北半球夏季日侧30°附近大气,北(夏)半球到达峰值的时间早于南(冬)半球.而可能受半球不对称背景磁场强度所导致的热层能量输送率影响,北半球夏季强磁暴和中磁暴个例的日侧,南半球δρa强于北半球;春秋季个例中日侧30°附近大气,北半球先于南半球1~2 h达到峰值. 2)受叠加在背景环流上的暴时经向环流影响,春秋季暴时赤道大气密度达到峰值的时间最短,日/夜侧大气分别在Dstmin后1 h和2 h达到峰值.至点附近夜侧赤道大气达到峰值时间一致,为Dstmin后3 h;不同季节日侧结果不同,在北半球冬季时赤道地区经过更长的时间达到峰值. 3)日侧赤道峰值时间距离高纬度峰值时间不受季节影响,为3 h左右.在春秋季和北半球冬季夜侧,赤道大气密度先于高纬度达到峰值,且不同纬度大气密度的峰值几乎无差别,表明此时低纬度存在其他加热源起着重要作用.     

Upper-atmospheric effects of magnetic storms: a brief tutorial

The physical processes underlying several phenomena of upper-atmospheric storms are described: magnetospherically driven ion convection and Joule heating and their impact on the high-latitude thermosphere and ionosphere; global changes in thermospheric circulation and composition; traveling atmospheric disturbances; and effects of electric-field penetration to middle and low latitudes. Examples from the 1997 January 10–11 storm are used to illustrate some of these features. It is pointed out that not only the magnitude, but also the sign of many storm-time changes at any given location depend sensitively on the temporal and spatial variations of auroral particle precipitation and high-latitude electric fields. In order for simulation models to be able to predict upper-atmospheric storm effects accurately, improved determination of the high-latitude inputs will be required.     

基于经验正交分析法的暴时热层大气密度时空分布规律

本文选取2002—2006年期间的36个强磁暴为研究对象,对CHAMP卫星加速度仪反演的实测大气密度进行经验正交分解,研究暴时热层大气密度的纬度分布特征,以及大气密度与ap指数、Dst指数的关系.结果表明,大气密度的纬度分布与季节相关,夏季半球的密度大于冬季半球,春秋季节南北半球的大气密度几乎对称分布;春秋季节白天大气密度在低纬地区呈现出赤道密度异常结构,在中高纬地区密度随纬度增加而减小,夜间则呈现抛物线的形状,赤道附近密度值最小.大气密度的纬度分布特征在若干天内具有良好的稳定性,发生时间相近的磁暴事件,纬度分布曲线非常相似,并且暴前与暴时的纬度分布变化不大.相关性分析表明,大气密度滞后ap指数2~6 h,相对Dst指数平均提前0~1 h,对磁暴的响应速度在日照区比在阴影区快,大气密度与ap指数、Dst指数具有较好的相关性.     

2006年4月13～17日西太平洋地区电离层暴时特性研究

本文利用西太平洋地区的15个电离层台站的测高仪数据和国际GPS服务中心IGS 36个站提供的TEC数据,以及由美国喷气推进动力学实验室提供的Jason-1 TEC数据对2006年4月13～17日间一次由冕洞高速流所引发的磁暴所造成的电离层效应进行了分析.分析结果表明这次电离层暴呈现出显著的纬度效应,foF2和TEC等参量显示在磁暴主相期间对称分布的强正暴效应中心在磁纬±30°～±40°,且持续时间超过12 h.负暴效应被限制在中高纬地区,在磁暴进入恢复相时,开始向低纬渗透,且具有明显的地方时效应.TIMIED卫星测量的Σ［O/N2］显示磁暴发生后,暴时环流使得中低纬地区的Σ［O/N2］有大幅增加,而中高纬地区则显著下降.通过对hmF2的分析发现磁暴主相期间,有磁层电场向中低纬地区穿透,且持续时间较长为1～3 h.因此这次强正暴效应可能是由风场、电场和化学成分这三个因素的共同作用造成的.这次磁暴造成的电离层暴响应非常复杂,对造成各种正负暴的物理和化学机制还需要进一步的研究.     

2003年11月超强磁暴热层大气密度扰动及其与焦耳加热和环电流指数的关系——CHAMP卫星观测

本文利用CHAMP卫星加速度仪测量数据,计算和分析2003年11月20～21日大磁暴期间大气质量密度扰动的全球分布特征;研究暴时变化与极区大尺度对流引起的全球焦耳加热总功率及环电流指数SYM-H之间的关系.结果表明,磁暴期间400 km高度上热层大气质量密度大幅度上升, NRLMSISE-00模式预测值与此相比有很大差别;暴时大气密度的增大存在昼/夜半球不对称性：白天强于夜晚,且白天随纬度的分布呈现出比较复杂的图像,在赤道附近和南半球中低纬区（10°N ～50°S）大气密度增大较强,并呈双峰分布,两个峰分别位于0°和45°S,另外在极区也出现大气密度扰动的局部极大,而在夜晚,大气密度变化南北半球比较对称,在赤道低纬区大气密度增大较强;互相关分析表明,中低纬区大气密度变化滞后于全球焦耳加热总功率3～7 h,滞后于环电流指数（SYM-H）0～3 h,与二者存在很强的相关,表明极区焦耳加热和赤道环电流过程对暴时热层大气密度扰动有重要影响.     

Steep latitudinal gradients of thermospheric composition during magnetic storms: a possible formation mechanism

Mass spectrometer satellite observations show that a narrow region with steep latitudinal gradients of neutral composition is formed in the subauroral winter thermosphere during magnetic storms. In order to analyze the relative importance of individual terms in the continuity equation for atomic oxygen, a two-dimensional model was used to simulate the thermospheric disturbance formation in response to intense Joule heating imposed in the auroral oval. Such an approach allowed three characteristic zones to be distinguished in the high-latitude thermosphere at heights of about 250 km. It was shown that vertical transport is of greatest importance within the local heating region. Horizontal transport dominates at subauroral latitudes near the mid-night edge of the auroral oval. Propagation of the disturbances to middle latitudes is prohibited near the noon edge of the oval by a strong counteraction of a poleward meridional wind. Here is a “relaxation zone” defined as the region which is spread to the equator from the boundary between the local heating area and the subauroral zone in the noon sector LT. It is at this boundary that composition distributions with steep latitudinal gradient are formed within the first few hours of Joule heating source action. Perturbations transported to middle latitudes during the periods when the meridional wind is directed equatorward begin to relax in this zone with a characteristic time scale of about 7 h, independent of season. However, in winter, composition at subauroral latitudes recovers to unperturbed N₂/O values before the wind again turns equatorward, giving rise to a distribution with steep latitudinal gradient recovering. In summer, a complete relaxation cannot be reached due to a shorter time interval with poleward wind and a larger disturbance amplitude. These two factors result in an effective smoothing of the initial steep gradient and a more regular latitudinal distribution of composition is observed in the summer thermosphere.     

Daytime F2-layer positive storm effect at middle and lower latitudes

Daytime F2-layer positive storm effects at middle and lower latitudes in the winter thermosphere are analyzed using AE-C, ESRO-4 neutral gas composition data, ground-based ionosonde observations and model calculations. Different longitudinal sectors marked by the storm onset as ‘night-time’ and ‘daytime’ demonstrate different F2-layer positive storm mechanisms. Neutral composition changes in the ‘night-time’ sector with increased [O] and [N₂] absolute concentrations, while (N₂/O)_storm/(N₂/O)_quiet\approx1 at F2-layer heights, are shown to contribute largely to the background N_mF2 increase at lower latitudes lasting during daytime hours. Storm-induced surges of the equatorward wind give rise to an additional N_mF2 increase above this background level. The mid-latitude F2-layer positive storm effect in the ‘daytime’ sector is due to the vertical plasma drift increase, resulting from the interaction of background (poleward) and storm-induced (equatorward) thermospheric winds, but not to changes of [O] and [N₂] concentrations.     

A three-dimensional model of ionospheric F-region

The three-dimensional stationary in solar-terrestrial system numerical model of the global ionosphere at F-region altitudes is developed. The input parameters are: the structural parameters of empirical thermospheric models (temperature, composition), electric fields, solar UV-radiation spectrum, corpuscular flows at high latitudes. The model includes the calculations of thermospheric circulation, electron density, electron and ion temperatures. The model reproduces the main morphological peculiarities of the distribution of thermospheric circulation and ionospheric plasma parameters. A comparison analysis of the results of the model calculation corresponding to different thermospheric models (DTM, MSIS, Jacchia-77, MSIS-83) is carried out. It is shown, that thermospheric circulation systems are the effective indicator of faithfulness of thermospheric models.     

Disturbances at F2-region heights of equatorial anomaly during geomagnetic storms

Neutral gas composition and ionospheric measurements taken by the Dynamic Explorer 2 satellite at F2-region heights during two geomagnetic storms are used to analyze the role of some possible physical mechanisms responsible for the changes of electron density at equatorial and low geomagnetic latitudes. The storms considered occurred on October 2, 1981 (storm 1) and July 13, 1982 (storm 2). During storm 1 (weak), vertical plasma drifts and equatorward storm-time winds operated increasing of the electron density at the trough of equatorial anomaly and the decreases at the crest region. During storm 2 (intense) changes of composition (increase of molecular nitrogen and atomic oxygen) played a fundamental role for the changes of electron density observed at low latitudes in summer hemisphere. It is concluded that different physical processes seem to have varying degrees of importance depending on the intensity of the storm.     

Annual and semiannual variations in the ionospheric F2-layer: II. Physical discussion

The companion paper by Zou et al. shows that the annual and semiannual variations in the peak F2-layer electron density (NmF2) at midlatitudes can be reproduced by a coupled thermosphere-ionosphere computational model (CTIP), without recourse to external influences such as the solar wind, or waves and tides originating in the lower atmosphere. The present work discusses the physics in greater detail. It shows that noon NmF2 is closely related to the ambient atomic/molecular concentration ratio, and suggests that the variations of NmF2 with geographic and magnetic longitude are largely due to the geometry of the auroral ovals. It also concludes that electric fields play no important part in the dynamics of the midlatitude thermosphere. Our modelling leads to the following picture of the global three-dimensional thermospheric circulation which, as envisaged by Duncan, is the key to explaining the F2-layer variations. At solstice, the almost continuous solar input at high summer latitudes drives a prevailing summer-to-winter wind, with upwelling at low latitudes and throughout most of the summer hemisphere, and a zone of downwelling in the winter hemisphere, just equatorward of the auroral oval. These motions affect thermospheric composition more than do the alternating day/night (up-and-down) motions at equinox. As a result, the thermosphere as a whole is more molecular at solstice than at equinox. Taken in conjunction with the well-known relation of F2-layer electron density to the atomic/molecular ratio in the neutral air, this explains the F2-layer semiannual effect in NmF2 that prevails at low and middle latitudes. At higher midlatitudes, the seasonal behaviour depends on the geographic latitude of the winter downwelling zone, though the effect of the composition changes is modified by the large solar zenith angle at midwinter. The zenith angle effect is especially important in longitudes far from the magnetic poles. Here, the downwelling occurs at high geographic latitudes, where the zenith angle effect becomes overwhelming and causes a midwinter depression of electron density, despite the enhanced atomic/molecular ratio. This leads to a semiannual variation of NmF2. A different situation exists in winter at longitudes near the magnetic poles, where the downwelling occurs at relatively low geographic latitudes so that solar radiation is strong enough to produce large values of NmF2. This circulation-driven mechanism provides a reasonably complete explanation of the observed pattern of F2 layer annual and semiannual quiet-day variations.     

1989年3月特强磁暴期间的电离层暴

用欧亚大陆地面电离层垂测站资料考察1989年3月12～16日磁暴期间的电离层暴形态及其发展变化. 特强磁暴引发的电离层暴是全球性的,但自磁层沉降的高能粒子对热层低部的加热程度及区域分布不同,因而各经度链区域内电离层暴的特征也有所差异. 本文研究表明,与理论推断对照,欧洲地区内F2层最大电子密度NmaxF2(或f0F2)并不出现正暴现象,而负暴自高纬向低纬的发展则与典型的热层环流结果相符. 此外,此磁暴过程期间在中低纬区存在明显的波动过程. 在亚洲高纬地区,磁暴初期13日有约10 h的正暴,而负暴过程则与欧洲地区类似,但不太清晰;且无波动现象. 磁暴期间,同一经度链的中低纬地区,夜间常发生多站同时的h′F突增. 本文再次证实这是一般磁暴期间常出现的普遍现象.     

大磁暴期间TIEGCM模式和CHAMP卫星热层大气密度扰动特性的统计研究

本文统计分析了2001—2005年的39次大磁暴事件(Dst-100nT)期间TIEGCM模式和CHAMP卫星大气密度数据.研究结果表明,模式结果与实测数据具有较好的一致性,但仍存在一定的偏差.大气密度及增量与SYM-H指数相关性较好,并且随纬度、光照条件和地磁活动水平变化.模式低估了磁暴期间大气密度的增幅,特别是在地磁活动水平较强时模式与实测的偏差较大.模式的偏差在高纬地区高于低纬地区,日侧高于夜侧.Dst指数越低,偏差越大,而当Dst指数低于-150nT以后,绝对偏差和相对偏差变化不明显.     

地磁扰动的时空变化规律和产生原因：统计和模拟研究

地磁感应电流（Geomagnetically Induced Currents,GICs）是发生在地球表面的一种空间天气现象,对石油管道、电缆等长距离地面基础设施会产生不利影响.研究GICs的产生原因及其与太阳风驱动条件的关系对于预报灾害性空间天气具有重要意义.由于GICs的产生与地磁场扰动紧密相关,本文利用北半球100多个地磁台站的长期监测数据,统计分析了地磁扰动及其时间变化率（磁扰率）与各种太阳风参数/地磁指数之间的相关性,并以2013年3月17日的磁暴事件为例,采用全球空间天气模型框架（space weather modeling framework,SWMF）模拟了磁暴期间地磁场扰动在北半球的分布以及日地空间中各电流体系对地磁扰动的作用.研究结果表明：（1）中低纬度的地磁扰动北向分量B_x与表征环电流强度的SYM-H指数呈现较高的正相关性（相关系数CC=0.75）,高于它与亚暴AE指数以及其他太阳风参数的相关性,说明磁暴期间环电流是导致中低纬度北向磁场减弱的主因,而在平静期间东向的磁层顶电流是中低纬度北向磁场增强的原因;（2）地磁扰动率与太阳风动压、行星际磁场IMF B_z、亚暴AE指数或者磁暴SYM-H指数均没有强相关性;（3）高纬地区的磁扰率通常大于低纬地区,而较强的磁扰率倾向于出现在中等磁暴或者中等-强亚暴期间的中高纬地区;（4）SWMF模型能较好地反演地磁平静时的北向地磁扰动和磁暴时的东向地磁扰动;（5）磁暴期间,磁层电流对中低纬度北向地磁扰动的贡献最大,而电离层霍尔电流对高纬地区的北向地磁扰动有着很强的支配地位;另外,高纬地区的东向地磁扰动主要由霍尔电流控制,而中低纬地区则受制于场向电流.

     

Ionospheric response to the geomagnetic storm on 13–17 April 2006 in the West Pacific region

This paper presents an investigation of geomagnetic storm effects in the equatorial and middle-low latitude F-region in the West Pacific sector during the intense geomagnetic storm on 13–17 April, 2006. The event, preceded by a minor storm, started at 2130 UT on April 13 while interplanetary magnetic field (IMF) B_z component was ready to turn southward. From 14–17 the ionosphere was characterized by a large scale enhancement in critical frequency, foF2 (4～6 MHz) and total electron content (TEC) (～30TECU, 1TECU=1×10¹⁶el/m²) followed by a long-duration negative phase observed through the simultaneous ionospheric sounding measurements from 14 stations and GPS network along the meridian 120°E. A periodic wave structure, known as traveling ionospheric disturbances (TIDs) was observed in the morning sector during the initial phase of the storm which should be associated with the impulsive magnetospheric energy injection to the auroral. In the afternoon and nighttime, the positive phase should be caused by the combination of equatorward winds and disturbed electric fields verified through the equatorial F-layer peak height variation and modeled upward drift of Fejer and Scherliess [1997. Empirical models of storm time equatorial electric fields. Journal of Geophysical Research 102, 24,047–24,056]. It is shown that the large positive storm effect was more pronounced in the Southern Hemisphere during the morning-noon sector on April 15 and negative phase reached to lower magnetic latitudes in the Northern Hemisphere which may be related to the asymmetry of the thermospheric condition during the storm.     

Ionospheric effects of magnetospheric and thermospheric disturbances on March 17–19, 2015

Using vertical and oblique radio-sounding data, we analyze the ionospheric and thermospheric disturbances during the magnetic storm that occurred in northeastern Russia on March 17–19, 2015. We consider the heliospheric sources that induced the magnetic storm. During the main and early recovery phases, the midlatitude stations are characterized by extremely low values of electron density at the F2 layer maximum. Using oblique sounding data, we recorded signals that propagated outside the great circle arc. In evening and night hours, no radio signals were found to pass along the Norilsk–Irkutsk and Magadan–Irkutsk paths. The observed ionospheric effects are shown to be caused by a sharp shift of the boundaries of the main ionospheric trough to the invariant latitude 46° N during the main phase of the magnetic storm. The negative ionospheric disturbance during the recovery phase of the storm, which was associated with significant variations in the composition of the neutral atmosphere, led to a change in the mode composition of received radio signals and a decline in observed maximal frequencies in daytime hours of March 18, 2015 by more than 2 times.     

The recovery phase of the superstrong magnetic storm of July 15–17, 2000: Substorms and ULF pulsations

The geomagnetic observations, performed at the global network of ground-based observatories during the recovery phase of the superstrong magnetic storm of July 15–17, 2000 (Bastille Day Event, Dst = ?301 nT), have been analyzed. It has been indicated that magnetic activity did not cease at the beginning of the storm recovery phase but abruptly shifted to polar latitudes. Polar cap substorms were accompanied by the development of intense geomagnetic pulsations in the morning sector of auroral latitudes. In this case oscillations at frequencies of 1–2 and 3–4 mHz were observed at geomagnetic latitudes higher and lower than ～62°, respectively. It has been detected that the spectra of variations in the solar wind dynamic pressure and the amplitude spectra of geomagnetic pulsations on the Earth’s surface were similar. Wave activity unexpectedly appeared in the evening sector of auroral latitudes after the development of near-midnight polar substorms. It has been established that the generation of Pc5 pulsations (in this case at frequencies of 3–4 mHz) was spatially asymmetric about noon during the late stage of the recovery phase of the discussed storm as took place during the recovery phase of the superstrong storms of October and November 2003. Intense oscillations were generated in the morning sector at the auroral latitudes and in the postnoon sector at the subauroral and middle latitudes. The cause of such an asymmetry, typical of the recovery phase of superstrong magnetic storms, remains unknown.     

Investigation of Ionospheric Response to Geomagnetic Storms over a Low Latitude Station,Ile-Ife,Nigeria

Due to several complexities associated with the equatorial ionosphere, and the significant role which the total electron content (TEC) variability plays in GPS signal transmission, there is the need to monitor irregularities in TEC during storm events. The GPS SCINDA receiver data at Ile-Ife, Nigeria, was analysed with a view to characterizing the ionospheric response to geomagnetic storms on 9 March and 1 October 2012. Presently, positive storm effects, peaks in TEC which were associated with prompt penetration of electric fields and changes in neutral gas composition were observed for the storms. The maximum percentage deviation in TEC of about 120 and 45% were observed for 9 March and 1 October 2012, respectively. An obvious negative percentage TEC deviation subsequent to sudden storm commencement (SSC) was observed and besides a geomagnetic storm does not necessarily suggest a high scintillation intensity (S₄) index. The present results show that magnetic storm events at low latitude regions may have an adverse effect on navigation and communication systems.     

Interplanetary origins of November 2004 superstorms

The sun was very active in the declining phase of solar cycle 23. Large sunspot active regions gave origin to multiple flare and coronal mass ejection (CME) activity in the interval 2003–2005. On November 2004, the active region AR 10696 was the origin of dozens of flares and many CMEs. Some events of this solar activity region resulted in two large geomagnetic storms, or superstorms (Dst??250 nT) on November 8, peak Dst=?373 nT, and on November 10, peak Dst=?289 nT. It is the purpose of this article to identify the interplanetary origins of these two superstorms. The southward-directed interplanetary magnetic fields (IMF Bs) that caused the two superstorms were related to a magnetic cloud (MC) field for the first superstorm, and a combination of sheath and MC fields for the second superstorm. However, this simple, classic picture is complicated by the presence of multiple shocks and waves. Six fast-forward shocks and, at least, two reverse waves were observed in the period of the two superstorms. A detailed analysis of these complex interplanetary features is performed in this work.