全球电离层对2000年4月6-7日磁暴事件的响应

利用分布于全球的电离层台站的测高仪观测数据,对扰动期间,foF2值与其宁静期间参考值进行比较,研究了2000年4月6—7日磁暴期间全球不同区域电离层的响应形态,并通过对比磁扰期间NmF2的变化与由MSISR90经验模式估算的中性大气浓度比(no／nN2)的变化,探讨了本次事件期间的电离层暴扰动机制．结果表明,在磁暴主相和恢复相早期,出现了全球性的电离层F2层负相暴效应．最大负相暴效应出现在磁暴恢复相早期,即电离层暴恢复相开始时间滞后于磁暴恢复相开始时间．在磁暴恢复相后期,一些台站出现正相扰动．研究结果表明,本次事件期间的电离层暴主要是由磁暴活动而诱发的热层暴环流引起的．     

中低纬电离层不同季节地磁暴响应的事例分析

利用中国中低纬台站漠河（53.5°N，122.3°E）、北京（40.3°N，116.2°E）、武汉（30.5°N，114.2°E）和三亚（18.3°N，109.6°E）的电离层观测数据，对比分析了4个台站电离层参数在2015年不同季节4个地磁扰动事件期间的变化特征.结果表明，4个磁暴事件期间电离层的响应特征并不完全一致，有着明显的季节特征，春季、夏季和秋季电离层以负相扰动为主，冬季以正相扰动为主.分析发现，中性成分O/N₂的降低与电离层负相扰动有关，但三亚地区的负相扰动还与扰动发电机电场相关.正相扰动的机制在不同事件中并不相同，穿透电场可能是引起春季磁暴事件期间电离层短时正暴效应的原因，而冬季长时间的正暴效应则是扰动电场和中性风共同作用的结果.      

磁暴期间中纬度电离层剖面结构变化的数值模拟

利用电离层理论模型模拟了磁暴期间热层大气温度、成分、中性风和电场扰动对电离层电子密度剖面结构，特别是峰值密度和峰值高度变化的影响，结果表明，热层大气温度变化所引起光化反应系数的改变对电离层剖面结构影响不大；热层大气成分特别是Ｎ２／Ｏ的变化能有效地引起密度剖面变化，Ｎ２增加足以使峰值密度产生所观测到的负相暴；由中性风和电场引起等离子体漂移是峰值高度ｈｍＦ２变化的主要原因，但对电子密度的影响不足以抵消     

磁暴期间全球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变化的相关性较差.这可能是由于太阳高度角较低,光辐射通量较小,导致电子密度的增加不明显.      

中纬度地磁暴期间热层垂直风响应机制的模拟

基于TIMEGCM模型,研究了2005年9月10日中纬度地磁暴期间热层(100～650 km)水平风场变化对垂直风的影响.通过连续性方程诊断分析了暴时引起垂直风场变化的机制,结果表明:250 km以上的垂直风场取决于水平风场的变化,而250 km以下的垂直风场由较高高度的垂直风拉动;在地磁暴初相开始时,经向风场相比纬向...     

磁暴期间中国中低纬电离层不规则体与扰动分析

2017年9月8日发生了一次强磁暴，Kp指数最大值达到8.利用区域电离层格网模型（Regional Ionosphere Map，RIM）和区域ROTI（Rate of TEC Index）地图，分析了磁暴期间中国及其周边地区电离层TEC扰动特征和低纬地区电离层不规则体的产生与发展情况，同时利用不同纬度IGS（International GNSS Service）测站BJFS（39.6°N，115.9°E），JFNG（30.5°N，114.5°E）和HKWS（22.4°N，114.3°E）的GPS双频观测值，获取各测站的ROTI和DROT（Standard Deviation of Differential ROT）指数变化趋势.结果表明:此次磁暴发生期间电离层扰动先以正相扰动为主，主要发生在中低纬区域，dTEC（differential TEC）最大值达到14.9TECU，随后电离层正相扰动逐渐衰减，在低纬区域发生电离层负相扰动，dTEC最小值达到-7.2TECU；在12:30UT-13:30UT时段，中国南部低纬地区发生明显的电离层不规则体事件；相比BJFS和JFNG两个测站，位于低纬的HKWS测站的ROTI和DROT指数变化更为剧烈，这表明电离层不规则体结构存在纬度差异.      

磁暴期间电离层电磁离子回旋波的地方时分布差异

利用Swarm卫星的高精度（50 Hz）磁场观测数据，对2015年3月16—25日磁暴期间中纬度电离层电磁离子回旋（EMIC）波时空分布特征进行了研究.结果表明:晨侧EMIC波事件数与昏侧大致相当，午前时段明显多于子夜前时段.昏侧EMIC波高发生率与等离子体羽状结构有关，晨侧EMIC波高发生率与太阳风动压增强及稠密冷等离子体有关.晨侧-正午前EMIC波频率高于昏侧-子夜前，表明源区位置以及离子成分占比存在地方时差异.昏侧事件大多发生在早期恢复相，晨侧事件大多发生在晚期恢复相，晨-昏两侧的时间差异源于磁暴期间高能离子西向漂移所需时间及等离子体层顶位置的地方时差异.磁暴期间，EMIC波以H+波和He+为主，其中H+波主要分布在06:00 MLT—10:00 MLT（磁地方时）扇区，He+波主要分布在18:00 MLT—22:00 MLT扇区.在磁暴主相期间没有出现H+带波，但是出现He+-O+双波段EMIC波，表明磁暴主相期间环电流高浓度氧离子对H+带EMIC波具有抑制作用.      

第23太阳活动周太阳极紫外辐射与磁暴相关性分析

统计第23个太阳活动周内中等及以上强度（Dst_min<-50nT）的磁暴事件，线性拟合分析磁暴主相DDst_min和达到DDst_min前一个表征太阳极紫外辐射强度的F_10.7之间的相关性.结果表明:随着太阳极紫外辐射增强，DDst_min<-50nT的磁暴出现的总数增多，在弱、中等和强太阳极紫外辐射条件下，其数量分别为56，84和85；随着太阳极紫外辐射增强，强磁暴（-200nT ≤ Dst_min<-100nT）和大磁暴（Dst_min<-200nT）发生的数量和相对发生率呈增长趋势，尤其是大磁暴数目（1，4，12）和相对发生率（1.79%，4.76%，14.12%）明显呈增长趋势；大磁暴（|Dst_min|）与太阳极紫外辐射（F_10.7）之间存在中度正相关关系，其相关系数为0.532，并且主要体现在大磁暴（|Dst_min|）与强太阳极紫外辐射（F_10.7）之间的中度正相关性，其相关系数为0.582.大磁暴与强太阳极紫外辐射之间的相关性可为空间天气预报提供参考依据.      

地磁暴期间中国中低纬电离层的CIT研究

利用电离层层析成像技术(Computerized Ionospheric Tomography, CIT)处理115°E子午圈附近6个台站的GPS观测数据, 分析了2004年11月地磁暴期间中国中低纬电离层的响应情况. 结果表明, 电离层呈正相扰动, 且不同高度上的响应不同, 800 km以下电子密度有不同程度的增加, 且在峰值高度附近增幅最大, 800 km以上地磁暴的影响并不显著; 伴随地磁能量的注入, 赤道异常峰极向扩展; 随磁扰强度的降低, 电子密度也逐渐恢复至平静水平. 这些结果与以往的理论和观测结果一致, 初步估计扰动是由热层暴环流引起的, 并受到赤道异常峰移动的影响.      

平流层爆发性增温事件中大气准16日行星波

2014年1月上中旬高纬平流层发生弱增温事件,增温幅度约25K,纬向西风减弱并于2月初转向.行星波在平流层爆发性增温（SSW）事件产生中具有重要作用.利用北半球近东经120°链上中低纬5个流星雷达探测的风场数据,研究了此SSW事件发生前和发生期间中间层和低热层区（MLT）大气风场的行星波状况.结果显示,极区平流层增温前MLT区大气呈现出明显增强的准16日波动,增温达到最大时,16日波也最强,表明中低纬MLT区的行星波变化与SSW事件存在耦合关系.进一步利用欧洲中心平流层再分析资料数据,分析SSW期间北半球平流层的波动和零风线状况,发现平流层准16日波和零风线随时间由低纬向高纬移动,反映出16日波与SSW之间存在某种动力学联系.      

Ionospheric parameters in the European sector during the magnetic storm of August 25–26, 2018

Variations of ionospheric parameters Total Electron Content (TEC) by GNSS, critical frequency (foF2) by vertical sounding and electron density (Ne) by low-altitude satellite were studied at high, mid and low latitudes of the European sector during the magnetic storm of August 25–26, 2018. During the main phase of the storm the ionospheric F2-layer was under the positive disturbance at mid and low latitudes. Then the transition from the positive to negative ΔfoF2 values occurred at all latitudes. The recovery phase was characterized by negative ionospheric disturbance at all latitudes. This is due to the decrease of thermospheric O/N2 ratio during the recovery phase of the storm. The intense Es layers screened the reflections from the F2-layer on August 26th at high and at low latitudes but at different times. Some blackouts occurred due to the high absorption level at high latitudes. In general, foF2 and TEC data were highly correlated. The major Ne changes were at the low latitudes. In general, Ne data confirmed the ionospheric dynamics revealed with foF2 and TEC.     

Aspects related to variability of radiative cooling by NO in lower thermosphere,TEC and O/N2 correlation,and diffusion of NO into mesosphere during the Halloween storms

Nitric Oxide is a very important trace species which plays a significant role acting as a natural thermostat in Earth’s thermosphere during strong geomagnetic activity. In this paper, we present various aspects related to the variation in the NO Infrared radiative flux (IRF) exiting the thermosphere by utilizing the TIMED/SABER (Thermosphere Ionosphere Mesosphere Energetics and Dynamics/ Sounding of the Atmosphere using Broadband Emission Radiometry) observational data during the Halloween storm which occurred in late October 2003. The Halloween storm comprised of three intense-geomagnetic storms. The variability of NO infrared flux during these storm events and its connection to the strength of the geomagnetic storms were found to be different in contrast to similar super storms. The connection between the quantum of energy outflux from the upper atmosphere into space in terms of NO IRF and the duration of storms is established. The NO radiative cooling, and the closely correlated depletion in O/N₂ ratio are controlled by the Joule heating intensity (proxied by AE-index). The collisional excitation rate of NO, calculated using the modelled datasets of WACCM-X (Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), correlates well with the observed pattern of radiative emission by NO. Observational datasets from TIMED/GUVI (Global Ultra-Violet Imager) and MIT Haystack observatory madrigal GNSS (Global navigation satellite system) total electron content (TEC) database shows that the TEC and O/N₂ enhancement in low-mid northern hemispheric latitudes are mainly controlled by the z-component of Interplanetary magnetic field (IMF-B_z). The penetration of eastward electric field during the storm events is found to be responsible for the overall enhancement of TEC. The contribution of enhanced day-side TEC in observed variation of O/N₂ ratio by GUVI is also reported. It is also seen that during substorms related events the night-time polar region experiences more cooling due to NO than the daytime polar region. The connections between the mid- and low-latitude enhancement in NO IRF with the propagation of LSTIDs (Large-scale traveling ionospheric disturbances) in combination with the O/N₂ variability, and the altitudinal variation in NO flux with the progression of the storm is also investigated. This study presents the evidence on the role of diffusion processes in the large scale enhancement of NO in the mesospheric altitudes.     

Some peculiarities of longitudinal distribution of proton fluxes at high latitudes

Dynamical features of proton fluxes at high and middle latitudes were studied based on data measured by Sun-synchronous low-altitude (800 km height) polar-orbiting three NOAA series satellites: POES 15, 16, and 17 during the geomagnetic storm on January, 21–22, 2005. Data from three satellites that passed the Northern hemisphere along different MLTs allow reconstructing the longitudinal distribution of the proton fluxes. Measurements of protons with energies of 30–80 keV and 80–240 keV (the ring current energy range) by 0- and 90-detectors were used to evaluate and compare the longitudinal asymmetry of proton flux distribution measured in the regions equatorward and poleward of the isotropic boundary. It was found that during all the phases of the geomagnetic storm distribution of the maximum flux of precipitating protons (0-detector data) is sufficiently asymmetric. The maximal flux position along MLT is moving from pre-midnight sector in quiet time to post-midnight one before and during SSC and moving back during recovery phase. The longitudinal distribution of precipitation maxima demonstrates the local increase in afternoon sector (approximately at 13:30 MLT) and decrease in the dusk one during SSC. These features are evident consequence of the magnetosphere compression. To identify the origin of the particles, the locations of maximum fluxes have been projected to the magnetosphere. It was determined that during geomagnetic storm main and recovery phases maximum fluxes were measured at latitudes poleward of the isotropic boundary. To evaluate the trapped particle flux asymmetry, the particle fluences (90-detector data) were calculated along the satellite orbit from L = 2 to the isotropic boundary. The total fluences of trapped particles calculated along the satellite orbit show regular asymmetry between dusk and dawn during main and recovery phases. The maximal intensity of proton fluxes of both investigated populations located poleward and equatorward of the isotropic boundary is achieved during SSC. The total flux measured during crossing the anisotropic region can be considered as a proxy for ring current injection rate.     

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.     

Seasonal variations of GPS derived TEC at three different latitudes of the southern hemisphere during geomagnetic storms

Data from the archive of the International GNSS Services (IGS) were used to study the seasonal variations of Total Electron Content (TEC) over three stations located at different latitudes in the southern hemisphere during the geomagnetic storms of 11 January, 6 April, 8 June, and 13 October 2000, representing storms that occurred in summer, autumn equinox, winter and spring equinox, respectively. The percentage TEC deviation with respect to reference values differs substantially from season to season. A strong seasonal anomaly and clear equinoctial asymmetry in TEC response to the storms were observed. Weak and short-lived positive TEC deviations as well as strong and long-lasting negative trends were observed in summer storm during the main and recovery phases respectively over the high and low latitudes whereas in winter storm, the highest positive TEC deviations was recorded during the main phase over the entire latitudes. TEC enhancement dominated all the stations during the autumn (March) equinox storm while TEC depletion was majorly observed during the spring (September) equinox. All these variations find their explanations in the thermospheric composition change and circulation. Future work with direct or modeled measurement of atomic Oxygen to molecular Nitrogen ratio (O/N₂), large number of storms and other possible factors such as variations in storm’s intensity and local time dependence of the storm onset is expected to validate the observations in this study.     

Response of equatorial,low- and mid-latitude F-region in the American sector during the intense geomagnetic storm on 24–25 October 2011

In this paper, we present and discuss the response of the ionospheric F-region in the American sector during the intense geomagnetic storm which occurred on 24–25 October 2011. In this investigation ionospheric sounding data obtained of 23, 24, 25, and 26 October 2011 at Puerto Rico (United States), Jicamarca (Peru), Palmas, São José dos Campos (Brazil), and Port Stanley, are presented. Also, the GPS observations obtained at 12 stations in the equatorial, low-, mid- and high-mid-latitude regions in the American sector are presented. During the fast decrease of Dst (about ∼54 nT/h between 23:00 and 01:00 UT) on the night of 24–25 October (main phase), there is a prompt penetration of electric field of magnetospheric origin resulting an unusual uplifting of the F region at equatorial stations. On the night of 24–25 October 2011 (recovery phase) equatorial, low- and mid-latitude stations show h′F variations much larger than the average variations possibly associated with traveling ionospheric disturbances (TIDs) caused by Joule heating at high latitudes. The foF2 variations at mid-latitude stations and the GPS-VTEC observations at mid- and low-latitude stations show a positive ionospheric storm on the night of 24–25 October, possibly due to changes in the large-scale wind circulation. The foF2 observations at mid-latitude station and the GPS-VTEC observations at mid- and high-mid-latitude stations show a negative ionospheric storm on the night of 24–25 October, probably associated with an increase in the density of molecular nitrogen. During the daytime on 25 October, the variations in foF2 at mid-latitude stations show large negative ionospheric storm, possibly due to changes in the O/N2 ratio. On the night of 24–25, ionospheric plasma bubbles (equatorial irregularities that extended to the low- and mid-latitude regions) are observed at equatorial, low- and mid-latitude stations. Also, on the night of 25–26, ionospheric plasma bubbles are observed at equatorial and low-latitude regions.     

Response of the lower atmosphere to intense geomagnetic storms

In this short paper we examine the possible connection between atmospheric parameters measured at low and middle altitudes and geomagnetic storms occurred in 2000 and 2003. For that, from a chain of stations located near the meridian 60°W we compare the storm time values of temperature and wind speed with their standard deviation 2σ obtained from quiet time values. We observed statistically significant variations at several altitudes during the storm recovery phase and after it, both in neutral wind speed and temperature. The results obtained suggest that atmospheric parameters could be affected by geomagnetic storms.     

Ionospheric response to the 17 March and 22 June 2015 geomagnetic storms over Wuhan region using GNSS-based tomographic technique

By using the data of GNSS (Global Navigation Satellite System) observation from Crustal Movement Observation Network of China (CMONOC), ionospheric electron density (IED) distributions reconstructed by using computerized ionospheric tomography (CIT) technique are used to investigate the ionospheric storm effects over Wuhan region during 17 March and 22 June 2015 geomagnetic storm periods. F-region critical frequency (foF2) at Wuhan ionosonde station shows an obvious decrease during recovery phase of the St. Patrick’s Day geomagnetic storm. Moreover, tomographic results present that the decrease in electron density begins at 12:00 UT on 17 March during the storm main phase. Also, foF2 shows a long-lasting negative storm effect during the recovery phase of the 22 June 2015 geomagnetic storm. Electron density chromatography presents the evident decrease during the storm day in accordance with the ionosonde observation. These ionospheric negative storm effects are probably associated with changes of chemical composition, PPEF and DDEF from high latitudes.     

平流层臭氧和辐射场的季节分布特征

利用美国NCAR化学气候耦合模式WACCM3对平流层温度场、风场、臭氧及辐射场进行了模拟.结果表明,在适宜飞艇长期驻留的准零风层高度20～22km(对应大气压强范围为50～30hPa,以下均采用气压值表征对应大气高度),7-8月风速小于5m·s-1的风带可长期稳定在40°N以北.臭氧空间分布显示,在30hPa气压高度处中国地区臭氧浓度出现了带状分布,30hPa高度以下低纬度地区臭氧浓度低于中纬度地区.平流层太阳加热率的时空变化表明,在平流层上层,太阳加热率可达100×10-6K·s-1,而在平流层下层,只有10×10-6K·s-1.6-8月中国区域的太阳加热率大于9月;在100～30hPa高度内,中纬度地区太阳加热率高于低纬度地区,在30hPa高度以上,低纬度地区太阳加热率高于中纬度地区;8-9月30～40hPa高度处,太阳加热率的空间变化较小.在30hPa高度上,太阳加热率在40°N昼夜变化最大;50hPa高度处,太阳加热率的昼夜变化小于30hPa高度处,而且白天太阳加热率出现极大值的纬度明显靠北.平流层低纬度地区的长波加热率小于中纬度地区.青藏高原由于地形特殊,其6-7月的臭氧浓度、太阳加热率和长波加热率均小于同纬度其他地区.