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

F2层对地磁扰动的响应

利用37个电离层垂直探测站1974-1986年的数据,采用f0F2与地磁ap指数相关分析的方法,首次得到一个太阳活动周期各年东亚-澳大利亚扇区,欧洲-非洲扇区和美洲-东太平洋扇区F2层对地磁扰动响应随地磁纬度的分布.结果指出,地磁高纬和中纬地区为负响应,低纬和赤道地区为正响应,大约在±30°附近换向.最大正响应在磁赤道附近,最大负响应在地磁纬度±50°附近,最大负响应的幅度大于最大正响应的幅度.存在明显的经度差别和南北半球不对称性.     

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

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指数变化更为剧烈,这表明电离层不规则体结构存在纬度差异.      

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

利用中国中低纬台站漠河（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₂的降低与电离层负相扰动有关,但三亚地区的负相扰动还与扰动发电机电场相关.正相扰动的机制在不同事件中并不相同,穿透电场可能是引起春季磁暴事件期间电离层短时正暴效应的原因,而冬季长时间的正暴效应则是扰动电场和中性风共同作用的结果.      

2008—2009年电离层对重现型地磁活动的响应

利用2008—2009年的GPS TEC数据,分析了电离层对冕洞引起的重现型地磁活动的响应. 结果表明,在太阳活动低年,电离层TEC表现出与地磁 ap指数（采用全球3h等效幅度指数ap来表征）和太阳风速度相似的9天和13.5天短周期变化,表明TEC的这种短周期特性主要与重现型地磁活动相关. 地磁纬度和地方时分析表明,夜间高纬地区正负相扰动明显,中低纬地区则以正相扰动为主,较大的TEC变幅主要发生在南北半球高纬地区,夜间南半球高纬地区TEC变化相对ap指数变化有相位延迟. 白天中低纬地区正负相扰动明显,TEC短周期变化与ap指数变化相位基本一致. 2008年TEC的9天和13.5天周期变化幅度大于2009年.      

电离层等离子体泡的磁场扰动特征与识别

等离子体泡是赤道区电离层常见的密度耗空结构,同时伴随着背景磁场的增强和垂直背景磁场方向上的扰动,这些扰动具有阿尔芬波的特征.本文利用Swarm卫星的标量、矢量磁场以及等离子体观测数据,考察了等离子体泡产生的磁场扰动信号特征,发现磁场强度和平均场向分量的增强均与等离子体泡密度下降有较强的负相关性,说明磁场标量和矢量(平均...     

地球磁场对太阳风扰动响应的研究

分别对行星际激波、太阳风动压增大事件和减小事件的地球磁场响应进行了比较. 分析结果表明, 同步轨道磁场对太阳风扰动在向阳面产生较强的正响应, 在背阳面 响应较弱且有时会出现负响应, 地磁指数SYM-H对太阳风扰动的响应为正响应. 同时还得出, 向阳侧同步轨道磁场响应幅度d Bz与地磁指数响应幅度d SYM-H、上下游动压均方差均具有较好的相关性. 地磁指数响应幅度与同步轨道磁场响应幅度相关关系在激波和动压增大事件中具有一致性, 动压减小事件出 现明显差异, 这说明激波和动压增大事件在影响地球磁场方面具有某种共性.      

1994年2月21日行星际激波引起的磁暴

利用Imp-8,Geotail和Goes-6等卫星资料,研究了1994年2月21日0900UT到达地球磁层的行星际激波引起的磁暴期间,从太阳风向磁层传输能量的有关问题.结果指出:（1）南向行星际磁场（IMF）的长持续时间不是太阳风向磁层输能的必要条件.南北振荡的,较强IMF也能产生显著的能量传输;（2）行星际扰动磁场通过弓激波和磁层顶后扰动磁能增加,增幅将近5倍;（3）在磁层内扰动磁场的Bz分量在1×10-4Hz附近显著被吸收.这一低频扰动磁场可能是磁暴期间导致氧离子和质子等环电流粒子向内扩散并被加速的原因之一.     

风云三号C星GNOS北斗掩星电离层探测初步结果

利用风云三号卫星C星GNOS掩星探测仪电离层数据,分析了2013年10月FY-3C GNOS探测的北斗掩星电离层廓线分布,将2013年10月1日至2015年10月10日期间FY-3C GNOS观测的F₂层峰值电子密度（N_mF₂）与地面电离层测高仪观测结果进行对比,验证了FY-3C GNOS北斗电离层掩星的探测精度.结果表明,FY3-C GNOS北斗电离层掩星与电离层测高仪探测的N_mF₂数据相关系数为0.96,平均偏差为10.21%,标准差为19.61%.在不同情况下其数据精度有如下特征:白天精度高于夜晚;夏季精度高于分季,分季精度高于冬季;中纬地区精度高于低纬地区,低纬地区精度高于高纬地区; BDS倾斜同步轨道（IGSO）卫星精度高于同步轨道（GEO）卫星和中轨道（MEO）卫星.FY-3C GNOS北斗电离层掩星与国际上其他掩星电离层数据精度的一致性对GNSS掩星探测资料的综合利用具有重大意义.      

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

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

亚暴起始位置和膨胀相持续时间对行星际磁场Bz分量的响应

磁层亚暴是太阳风–磁层–电离层耦合过程中的重要爆发性事件,其特性受太阳风参数的影响很大。本文利用对IMAGE卫星在2000 － 2005年观测到的4193个亚暴起始事件,统计研究了在不同的行星际磁场（IMF）B_z 条件下亚暴起始位置和膨胀相持续时间。结果表明,南向IMF发生的亚暴比北向IMF下发生的亚暴要多。南向IMF条件下亚暴AE指数最大值的平均值基本上>600 nT,并有随南向IMF持续时间增大而增大的趋势。北向IMF条件下亚暴AE指数最大值的平均值基本上<500 nT,并有随北向IMF持续时间增大而减小的趋势。亚暴的起始磁纬度基本上位于65° － 70°之间。当南向IMF或北向IMF的持续时间增大,超过80 min时,北半球的亚暴起始磁纬度会降低。亚暴起始磁地方时大部分位于22:15 － 23:15 MLT之间。但整体分布比较分散,显示不出特别清晰的随IMF B_z持续时间变化的趋势。相比于南向的IMF,北向IMF期间发生亚暴的平均膨胀相持续时间增大了将近10 min,表明南向IMF期间,亚暴强度虽然较大,但其膨胀相持续时间较短,亚暴能量释放和耗散的速度更快。      

Statistical average estimates of high latitude field-aligned currents from the STARE and SABRE coherent VHF radar systems

Two bistatic VHF radar systems, STARE and SABRE, have been employed to estimate ionospheric electric fields in the geomagnetic latitude range 61.1 – 69.3° (geographic latitude range 63.8 – 72.6°) over northern Scandinavia. 173 days of good backscatter from all four radars have been analysed during the period 1982 to 1986, from which the average ionospheric divergence electric field versus latitude and time is calculated. The average magnetic field-aligned currents are computed using an AE-dependent empirical model of the ionospheric conductance. Statistical Birkeland current estimates are presented for high and low values of the Kp and AE indices as well as positive and negative orientations of the IMF B_z component. The results compare very favourably to other ground-based and satellite measurements.     

The solar wind control of the equatorial ionosphere dynamics during geomagnetic storms

The interplanetary magnetic field, geomagnetic variations, virtual ionosphere height h′F, and the critical frequency foF2 data during the geomagnetic storms are studied to demonstrate relationships between these phenomena. We study 5-min ionospheric variations using the first Western Pacific Ionosphere Campaign (1998–1999) observations, 5-min interplanetary magnetic field (IMF) and 5-min auroral electrojets data during a moderate geomagnetic storm. These data allowed us to demonstrate that the auroral and the equatorial ionospheric phenomena are developed practically simultaneously. Hourly average of the ionospheric foF2 and h′F variations at near equatorial stations during a similar storm show the same behavior. We suppose this is due to interaction between electric fields of the auroral and the equatorial ionosphere during geomagnetic storms. It is shown that the low-latitude ionosphere dynamics during these moderate storms was defined by the southward direction of the B_z-component of the interplanetary magnetic field. A southward IMF produces the Region I and Region II field-aligned currents (FAC) and polar electrojet current systems. We assume that the short-term ionospheric variations during geomagnetic storms can be explained mainly by the electric field of the FAC. The electric fields of the field-aligned currents can penetrate throughout the mid-latitude ionosphere to the equator and may serve as a coupling agent between the auroral and the equatorial ionosphere.     

TC-1和Cluster对向阳侧磁层顶通量传输事件的联合观测研究

2004年2至4月期间,探测一号(TC-1)卫星和Cluster卫星有25次同时处在向阳侧磁层顶附近的磁鞘内,TC-1卫星在低纬区,Cluster卫星在中高纬区.利用这一期间两卫星探测到的27个通量传输事件(FTE),分析行星际磁场(IMF)横向分量BT={By,Bz}对磁层顶重联发生位置的影响,以及分量重联的观测事实,得到如下主要结果.(1)当IMF南向分量Bz占优势(|Bz|＞|By|)时,FTE大多(约占87.5%)能在低纬观测到,而当IMF By分量占优势(|Bz|＜|By|)时,则FTE大部分能在中高纬观测到(占84.2%);(2)很少观测到相关联的事件(关联事件指在低纬生成的FTE,向高纬运动中先后被TC-1卫星和Cluster卫星探测到的事件),表明在低纬形成的FTE可能大多沿磁层顶两侧滑向磁尾,只有少数可能运动到高纬地区;(3)中纬地区探测到的FTE大多是以分量重联方式产生于该区,而非来自磁赤道附近成对形成的FTE.      

Electrodynamic Coupling in the Solar Wind-Magnetosphere-Ionosphere System

This paper presents a brief summary of our recent work based on global MHD simulations of the Solar wind-Magnetosphere-Ionosphere (SMI) system with emphasis on the electrodynamic coupling in the system. The main conclusions obtained are summarized as follows. (1) As a main dynamo of the SMI system, the bow shock contributes to both region 1 Field-Aligned Current (FAC) and cross-tail current. Under strong interplanetary driving conditions and moderate Alfven Mach numbers, the bow shock's contribution may exceed more than fifty percent of the total of either region 1 or cross-tail currents. (2) In terms of more than 100 simulation runs with due southward Interplanetary Magnetic Field (IMF), we have found a combined parameter f = E_swP_swM_A-1/2 (E_sw, P_sw, and M_A are the solar wind electric field, ram pressure, and Alfven Mach number, respectively): both the ionospheric transpolar potential and the magnetopause reconnection voltage vary linearly with f for small f, but saturate for large f. (3) The reconnection voltage is approximately fitted by sin3/2θ_IMF/2, where θ_IMF is the IMF clock angle. The ionospheric transpolar potential, the voltage along the polar cap boundary, and the electric fields along the merging line however defined they may be, respond differently to θ_IMF, so it is not justified to take them as substitutes for the reconnection voltage.      

Magnetic field disturbances in the mid-altitude cusp

Magnetic field disturbances and hot particles in the cusp as seen mainly by both the IMAP-3 magnetometer and PROMICS -3 spectrometer aboard the INTERBALL-AURORAL PROBE (AU) (perigee 4 Re, inclination 65 deg) are discussed. Orbits from March, 1997 are analysed accordingly. The INTERBALL-AU transits cross the middle of the cusp from low to high latitude from pre-noon to post-noon hours 11–13 MLT. Although the crossings are not exactly meridional, they reach 80–81 ILAT and permitting detection of the cusp, its subregions and their equatorward and poleward boundaries. Magnetic field disturbances are analysed for field-aligned currents (FACs) during different IMF conditions. Regions with structured forms are emphasised, in particular where a more intense current is concentrated. The following cases are compared: (a) the most intense current in the boundary cusp, adjacent to the ring current (March 7, southward IMF); (b) FACs are concentrated completely inside the cusp (March 13, southward -downward IMF) and (c) FACs with pronounced intensities equatorward from the cusp proper (March 11, duskward IMF).     

On the association of quiet-time Pi2 pulsations with IMF variations

The association of quiet-time Pi2 pulsations with the variations of the interplanetary magnetic field (IMF) has been examined by using three reported events, occurring during extremely quiet intervals, of which the first was on 10 March 1997, the second 27 December 1997, and the third 11 May 1999. For the first event, the onset time of ground Pi2s maps to the IMF structure bearing a variation cycle of north-to-south and north again as seen by Wind in the upstream region and Geotail in the magnetosheath. Likewise, the second and the third events have respectively, four and three recurrent turnings propagating to the Earth sensed by multiple satellites. The comparison of geomagnetic perturbations, auroral brightenings, and energetic particle data in the magnetotail with the IMF observations shows successive substorm-like activations accompanied by ground Pi2 onsets. For a clear variation cycle, the first Pi2 burst appears 36 ± 8 min after southward turning of the IMF and the second one follows14 ± 4 min after a northward turning. Moreover, ground Pi2 onsets recur under low IMF clock angle conditions. These observational results can be interpreted with the prevailing models of externally triggered substorm. But the solar wind coupling to the magnetosphere under quiet conditions proceeds in a less efficient way than under substorm time conditions. Consequently, we suggest that recurrent quiet-time Pi2s can be associated with IMF variations and their cause can be the same as those for substorm times.     

Evidence for ion energy dispersion in the polar cusp related to a northward-directed IMF

Data from the particle experiment aboard the AUREOL-3 polar satellite show that about 30% of the summer cusp crossings are characterised by a clear latitudinal energy dispersion of the solar wind ions. This energy-latitude correlation is observed at very high latitudes, 80° – 85°, near the polar boundary of the cusp, as an increase of the ion average energy with latitude. These structures have a typical latitude extent of 1° – 2° at ionospheric heights and correspond to a northward-directed IMF. These observations are consistent with a sunward convection of the foot of the magnetic flux tubes recently merged with a northward directed interplanetary magnetic field.     

A comparative analysis of the role of interplanetary magnetic field (IMF) and sudden impulse (SI) in triggering a substorm

During a typical Akasofu-type of substorm, the southward component of IMF Bz is necessary prior to the onset. However, a sudden compression of solar wind, if intense enough, can also sometimes trigger a substorm, and is independent of the IMF orientation. The Akasofu-type substorm and the Impulse-induced substorm may differ in their occurrence mechanism and ground-based observations. This is shown using the initial four substorm events discussed in this paper having distinctly different IMF and sudden impulse conditions. A question then arises is how will these signatures vary when both sudden impulse and a southward component of IMF Bz are present prior to the onset. To account for the same, we analyze two substorm events of 05th April 2010 and 22nd June 2015. The substorm onsets on these days not just coincided with the sudden impulse but also a southward component of IMF Bz was present prior to the onsets. The present study accounts for the similarities and differences among isolated IMF induced substorms, isolated impulse-induced substorms and when both sudden impulse and a southward component of IMF Bz are present. We examined the relative dominance between the two factors in triggering a substorm using ground-based and satellite-based observations. If IMF Bz is near zero, a strong pressure pulse and/or large IMF By can lead to particle precipitation away from the usual midnight. To further ensure whether a pressure pulse or IMF By predominantly influences the substorm onset location, a statistical analysis of isolated substorms will be needed.     

Interplanetary magnetic field control of plasma distribution in the polar ionosphere

We used the TEC (Total electron content) data of 5 min resolution obtained from the Madrigal database during solar-maximum winter (Nov. 6, 2000–Feb. 4, 2001) to study statistically the polar ionospheric plasma distribution response to different intensity and orientation of IMF By/Bz components. The sunlit high-density plasma extension from dayside to nightside is favored in negative IMF By and Bz conditions. With the magnitude of the negative Bz increasing, the time range corresponding to the distinct high-density extension feature expands, and the plasma density along the extension path enhances, which can be attributed to the interaction between dayside solar-produced ionization whose poleward limit is decided by terminator and convection extent mainly modulated by IMF Bz component. As for IMF By component influence on the sunlit plasma extension, the combination effect of convection and corotation electric fields is necessary to be considered.