2009年1月平流层爆发性增温期间全球电离层响应的研究

2009年1月平流层爆发性增温(Stratospheric Sudden Warming, SSW)事件是有记录以来最强、持续时间最长的一次主增温事件(Major Warming Event, MWE), 期间太阳活动和地磁活动均处于较低的水平, 因此非常有利于研究电离层对平流层增温事件的响应情况. 本文利用COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate)系统提供的掩星数据, 使用Kriging方法分别构建了此次SSW期间及平静期的电离层N_mF₂, h_mF₂和110～750km高度范围的垂直积分TEC (简称VTEC)地图, 同时从全球定位导航卫星系统服务组织(International GNSS Service, IGS)发布的全球电离层TEC地图(Global Ionospheric Maps, GIMs)中提取了日固坐标系(Sun-fixed, 地磁纬度和地方时)下全球TEC地图. 通过对比发现, SSW期间与平静期相比, 地磁纬度中低纬电离层参数存在早晨上升, 下午和夜间下降的现象. 利用OSTM/JASON-2卫星高度计观测值进行验证后的结果显示, 此前研究均未有提及的夜间时段电离层参数N_mF₂, h_mF₂和TEC (VTEC和IGS TEC)的下降现象的确存在.      

A Prediction Method for△f0F2min of a Single Station From Interplanetary Parameters Based on Statistical Study

Using 86 CME-interplanetary shock events,the correlation between the peak values of （a） the solar wind parameters(B_z,E_y,P_dyn) and the geomagnetic indices（SYM-H,ASY-H,Kp）, （b） the coupling functions（Borovsky,Akasofu,Newell） and the geomagnetic indices,（c） the solar wind parameters/coupling functions/geomagnetic indices and the ionospheric parameter(Δf₀F_2min), are investigated.The statistical results show that in group（a）,B_{z min} and SYM-H_min have the best correlation,that in group（b）,the best correlation is between the peak values of Akasofu function (A_min) and SYM-H_min,and that in group（c）,the best correlation is between K_pmax andΔf₀F_2min. Based on the statistical results,a method for predicting f₀F₂ of a single station is attempted to be set up.The input is modified B_{z min} and the outputs are SYM-H_min andΔf₀F_2min.Then 25 CME-IPS events that caused geomagnetic storms in 1998 and 2009 are used to check the prediction method. The results show that our method can be used to predict SYM-H_min andΔf₀F_2min.     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

Magnetospheric disturbances associated with the 13 December 2006 solar flare and their ionospheric effects over North-East Asia

We present an observational study of magnetospheric and ionospheric disturbances during the December 2006 intense magnetic storm associated with the 4В/Х3.4 class solar flare. To perform the study we utilize the ground data from North–East Asian ionospheric and magnetic observatories (60–72°N, 88–152°E) and in situ measurements from LANL, GOES, Geotail and ACE satellites. The comparative analysis of ionospheric, magnetospheric and heliospheric disturbances shows that the interaction of the magnetosphere with heavily compressed solar wind and interplanetary magnetic field caused the initial phase of the magnetic storm. It was accompanied by the intense sporadic E and F2 layers and the total black-out in the nocturnal subauroral ionosphere. During the storm main phase, LANL-97A, LANL 1994_084, LANL 1989-046 and GOES_11 satellites registered a compression of the dayside magnetosphere up to their orbits. In the morning–noon sector the compression was accompanied by an absence of reflections from ionosphere over subauroral ionospheric station Zhigansk (66.8°N, 123.3°E), and a drastic decrease in the F2 layer critical frequency (foF2) up to 54% of the quite one over subauroral Yakutsk station (62°N, 129.7°E). At the end of the main phase, these stations registered a sharp foF2 increase in the afternoon sector. At Yakutsk the peak foF2 was 1.9 time higher than the undisturbed one. The mentioned ionospheric disturbances occurred simultaneously with changes in the temperature, density and temperature anisotropy of particles at geosynchronous orbit, registered by the LANL-97A satellite nearby the meridian of ionospheric and magnetic measurements. The whole complex of disturbances may be caused by radial displacement of the main magnetospheric domains (magnetopause, cusp/cleft, plasma sheet) with respect to the observation points, caused by changes in the solar wind dynamic pressure, the field of magnetospheric convection, and rotation of the Earth.     

Low-latitude ionospheric scintillations and total electron content obtained with the CITRIS instrument on STPSat1 using radio transmissions from DORIS ground beacons

The primary objective of the Scintillation and Tomography Receiver in Space (CITRIS) is to detect ionospheric irregularities from space at low latitude. For this purpose, the satellite receiver uses the UHF and S-Band transmissions of the ground network of Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) beacons. CITRIS, developed at the Naval Research Laboratory, differs from the normal DORIS receiver by being able to capture and store the complex amplitude of the 401.25 and 2036.25 MHz transmissions at 200 Hz sample rate. Ground processing of the CITRIS data yields total electron content (TEC) and both phase and amplitude scintillations. With CITRIS flying on the US Space Test Program (STP) satellite STPSat1, 2 years of data were collected and processed to determine the fluctuations in ionospheric TEC and radio scintillations associated with equatorial irregularities. CITRIS flights over DORIS transmitters yield direct measurements of the horizontal plasma density fluctuations associated with equatorial plasma bubbles. Future flights of CITRIS can provide valuable complements to other satellite instruments such as GPS occultation receivers used to estimate vertical electron density profiles in the ionosphere.     

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

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

基于空间变异特性的格网电离层延迟估计方法

广域增强系统采用格网模型校正用户电离层延迟误差,为保证理想的校正效果,必须精确估计格网点电离层延迟。电离层延迟的空间相关性是完成这一估计的理论依据。针对传统反距离加权法因不能准确描述这种空间相关性而使得估计效果受限的问题,依据空间统计学原理建立了电离层延迟空间变异函数模型,提出一种新的权值计算方法,并通过实际观测数据进行了仿真验证。结果表明,相比于传统方法,这一方法能明显改善格网点电离层延迟估计效果,尤其在穿透点观测数据稀疏的区域,其改善效果更为明显。     

Ionospheric electron density over Resolute Bay according to E-CHAIM model and RISR radar measurements

In this study, predictions of the E-CHAIM ionospheric model are compared with measurements by the incoherent scatter radars RISR at Resolute Bay, Canada, in the northern polar cap. Reasonable coverage was available for all seasons except winter for which no conclusions were drawn. It is shown that ratios of the model-to measured electron densities are close to unity in the central part of the F layer, around its peak. This is particularly evident for summer daytime. Distributions of the ratios are wider for other seasons indicating larger number of cases when the model underestimates or overestimates. E-CHAIM underestimates the electron density at ionospheric topside and bottomside by ～ 10–20 %. At the bottomside, the underestimations are strongest in summer and equinoctial nighttime. At the topside, the underestimations are strongest in autumn nighttime. Model overestimations are noticeable in the middle part of the F layer during dawn hours in autumn. Overall, the model tends to not predict highest-observed peak electron densities and the largest-observed heights of the peak.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

Symmetry and asymmetry of ionospheric weather at magnetic conjugate points for two midlatitude observatories

Variations of the ionospheric weather W-index for two midlatitude observatories, namely, Grahamstown and Hermanus, and their conjugate counterpart locations in Africa are studied for a period from October 2010 to December 2011. The observatories are located in the longitude sector, which has consistent magnetic equator and geographic equator so that geomagnetic latitudes of the line of force are very close to the corresponding geographic latitudes providing opportunity to ignore the impact of the difference of the gravitational field and the geomagnetic field at the conjugate points on the ionosphere structure and dynamics. The ionosondes of Grahamstown and Hermanus provide data of the critical frequency (foF2), and Global Ionospheric Maps (GIM) provide the total electron content (TECgps) along the magnetic field line up to the conjugate point in the opposite hemisphere. The global model of the ionosphere, International Reference Ionosphere, extended to the plasmasphere altitude of 20,200 km (IRI-Plas) is used to deliver the F2 layer peak parameters from TECgps at the magnetic conjugate area. The evidence is obtained that the electron gas heated by day and cooled by night at the summer hemisphere as compared with the opposite features in the conjugate winter hemisphere testifies on a reversal of plasma fluxes along the magnetic field line by the solar terminator. The ionospheric weather W-index is derived from NmF2 (related with foF2) and TECgps data. It is found that symmetry of W-index behavior in the magnetic conjugate hemispheres is dominant for the equinoxes when plasma movement along the magnetic line of force is imposed on symmetrical background electron density and electron content. Asymmetry of the ionospheric storm effects is observed for solstices when the plasma diffuse down more slowly into the colder winter hemisphere than into the warmer summer hemisphere inducing either plasma increase (positive phase) or decrease (negative phase of W-index) in the ionospheric and plasmaspheric plasma density.