利用TIE-GCM模式计算电离层电流及夜间电离层磁场

电离层电流产生的磁场是地磁场卫星测绘时需要剔除的干扰源.利用电离层热层模式TIE-GCM计算电离层中的中性风、重力驱动和压强梯度等形成的电离层电流的全球分布,分析电流在特定位置产生的磁场及磁场三分量随纬度的变化规律.结果表明,E层尤其是磁赤道和极区的电流密度较大,可达103nA·m-2量级,F层电流密度量级约为10nA·m-2.在磁静日（Kp≤ 1）夜间22:00LT-04:00LT,电离层电流在中低纬度（南北纬50°之间）产生的磁场量级为几个nT,且磁场的南北向分量和径向分量基本大于东西向分量.通过与CHAMP卫星磁测数据分析比较,发现TIE-GCM模式计算电离层干扰磁场在中低纬度可以取得较好的结果,但在高纬度地区的效果不理想,还需进一步改进模式以提高计算精度.      

CAS-DIS数字电离层测高仪系统研制

提出了将射频直接采样数字化的接收方式应用于数字电离层测高仪系统的设计方案, 并详细介绍了该方案在CAS-DIS (Chinese Academy of Sciences, Digital Ionosonde)电离层测高仪系统的应用, 给出了CAS-DIS电离层测高仪在中国武汉电离层探测标校试验场和武汉-北京斜向电离层探测系统中的实验结果. 观测结果表明, 采用本文提出的数字化方案设计的电离层测高仪系统性能优越, 且可以分别工作在垂直和斜向探测模式, 满足了电离层探测的多种应用需求.      

电离层数字测高仪被动接收观测模式研究

利用CADI(Canadian Advanced Digital Ionosonde)电离层数字测高仪平台,实现了新的电离层数字测高仪被动接收观测模式.利用新开发的观测模式,在观测台站开展了一系列实验观测研究,经过信号处理和信息提取,获得了电离层特征参量f0F2回归方程,高频信道背景噪声分布,电离层D层对电波的吸收等电离层探测信息.实验观测结果表明,所获取的f0F2与主动探测结果相关性在0.84以上,高频信道背景噪声分布以及电离层D层吸收状况与电离层实际分析结果相吻合.     

国际参考电离层模型输出参数峰值高度性能

电离层峰值高度H_mF₂是描述电离层形态的重要参数之一，国际参考电离层模型IRI-2016中融入了大量电离层测高仪和无线电掩星探测数据，用以提升H_mF₂的预测精度.本文利用太阳活动低年（2007—2010年）气象、电离层和气候卫星联合观测系统COSMIC探测数据描述全球范围内COSMIC H_mF₂的三维形态变化，对比分析了IRI-2016与IRI-2012模型的预测结果，同时分析了IRI-2016模型输出H_mF₂的性能.结果表明，IRI模型在中高纬度地区的输出结果高于COSMIC反演结果，而赤道及低纬地区则大都偏低.与IRI-2012模型相比，IRI-2016模型的输出结果在夜间至凌晨时段呈现较为明显的纬向梯度变化且大部分区域输出值偏高，但在白天时段赤道附近区域的输出值大都偏低.上述结果为电离层IRI模型的完善提供了一定参考.      

中纬电离层理论模式研究

本文利用谱方法建立起一维时变中纬度电离层理论模式．模式比较周密地考虑了中纬度地区电离层的主要动力学过程和光化学过程．模式的突出优点是计算耗时少．我们对日本Wakkanai站进行了模拟计算,模拟结果同实验结果符合得较好．     

1995年10月24日日食电离层效应的模式化研究

利用一个二维纬电离层理论模式模拟１９９５年１０月２４日日全食条件下东亚低纬电离层演变,探索低纬日食电离层效应。模式计算结果显示：日食散间日食区电离层电子与离子光致电离产生率减少使得Ｅ－Ｆ１区电离层光食效应显著,而较高部份电离层对日食影响的响应延迟。     

电离层短波射线追踪

本文以电离层等离子体参量的平均背景模式计算为基础, 给出一种电离层短波射线轨迹计算方法。供短波通信系统、工作于短波段的其它系统(如HF雷达的目标定位系统)及有关研究工作使用参考。      

电离层薄层高度对电离层模型化的影响

利用IRI2012模型分析了电离层薄层高度的时空变化规律,提出了基于应用中STEC的电离层改正误差分析理论,分析了电离层薄层高度变化的相关影响.结果表明,电离层薄层高度变化对电离层穿刺点位置、投影映射函数值、电离层建模结果、电离层模型精化和电离层模型精度评估结果的影响较大.高度截止角为10°时,电离层薄层高度变化导致电离层穿刺点的经纬度差异最大可达3.2°,投影映射函数最高可引入约15.46%的误差,电离层建模结果差异和建模实用误差最高分别达9.71%,3.64%,采用不同薄层高度数据的电离层模型参数拟合和模型精化结果最大可引入约9.26%的误差,采用不同电离层薄层高度数据进行模型精度评定时最大可引入约9.62%的误差.根据这些研究结果可知:在实际应用中应采用电离层薄层高度模型,并选取较大的卫星高度截止角来减小薄层高度变化引入的误差;采用固定高度时,区域电离层建模采用与实际电离层薄层一致的固定高度;进行精度评估时,参考数据的电离层薄层高度与需要精度评估的电离层模型薄层高度相等.      

用经验模式参数研究电离层发电机效应

本文利用MSIS-86和IRI-86模式的基本参数,求得热层风系和电导率的三维分布;继而从发电机理论得到了电离层电位及层电流密度在北半球的分布和变化。本研究将大量数据统计平均的大气模式参数与热层、电离层理论研究联系了起来。     

化学物质释放人工改变电离层

考虑中性气体在电离层高度的扩散过程和相应的电离层离子化学过程，研究了利用主动化学物质释放来改变电离层的方法，理论计算了H2O和SF6两种气体释放后电离层随时间的响应过程，结果表明，在电离层高度上气体的扩散过程非常迅速，电离层F区的电子密度有很大程度的减少，而扩散慢且化学反应快的气体对电离层的影响更大，就更加有利于电离层洞的形成。     

国际参考电离层IRI－90与中国实测电离层月平均电子密度剖面形状的比较

将太阳活动峰年期间中国4个电离层站垂直探测得到的月平均电子浓度剖面与国际参考电离层IKI-90进行了系统的比较。结果表明国际参考电离层所计算的峰下电子含量(或峰下半厚)总体来说偏大。一天中白天符合较好,晚上较差,对中纬台站较好,对低纬台站较差。      

Evaluation of the improvement of IRI-2016 over IRI-2012 at the India low-latitude region during the ascending phase of cycle 24

This paper investigated the performance of the latest International Reference Ionosphere model (IRI-2016) over that of IRI-2012 in predicting total electron content (TEC) at three different stations in the Indian region. The data used were Global Positioning System (GPS) data collected during the ascending phase of solar cycle 24 over three low-latitude stations in India namely; Bangalore (13.02°N Geographic latitude, 77.57°E Geographic longitude), Hyderabad (17.25°N Geographic latitude, 78.30°E Geographic longitude) and Surat (21.16°N Geographic latitude, 72.78°E Geographic longitude). Monthly, the seasonal and annual variability of GPS-TEC have been compared with those derived from International Reference Ionosphere IRI-2016 and IRI-2012 with two different options of topside electron density: NeQuick and IRI01-corr. It is observed that both versions of IRI (i.e., IRI-2012 and IRI-2016) predict the GPS-TEC with some deviations, the latest version of IRI (IRI-2016) predicted the TEC similar to those predicted by IRI-2012 for all the seasons at all stations except for morning hours (0500 LT to 1000?LT). This shows that the effect of the updated version is seen only during morning hours and also that there is no change in TEC values by IRI-2016 from those predicted by IRI-2012 for the rest of the time of the day in the Indian low latitude region. The semiannual variations in the daytime maximum values of TEC are clearly observed from both GPS and model-derived TEC values with two peaks around March-April and September-October months of each year. Further, the Correlation of TEC derived by IRI-2016 and IRI-2012 with EUV and F10.7 shows similar results. This shows that the solar input to the IRI-2016 is similar to IRI 2012. There is no significant difference observed in TEC, bottom-side thickness (B0) and shape (B1) parameter predictions by both the versions of the IRI model. However, a clear improvement is visible in hmF2 and NmF2 predictions by IRI-2016 to that by IRI-2012. The SHU-2015 option of the IRI-2016 gives a better prediction of NmF2 for all the months at low latitude station Ahmedabad compare to AMTB 2013.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

The state of the art in D-region IRI modelling

While describing the formulation of IRI-79, in so far as it is applicable to the sub E-peak electron density profile, an attempt is made to compare its predictions with the available experimental evidence. The improvements needed for D-region electron density modelling are summarized. The current D-region IRI modelling effort is illustrated with a typical example.     

Comparison between d- and lower e-region electron density profiles and IRI-79

The D and E-region electron density profiles obtained by different techniques are compared with the IRI-79 model to see how they fit. The rocket data showed good agreement. However discrepancies between the observed and model values were found especially for solar zenith angle greater than 50 degrees.     

Improvement of the solar-cycle variation of IRI lower ionosphere models by means of radio wave propagation data

Data on day-time and night-time radio wave absorption in the frequency range 50 to 2614 kHz, obtained in long-term observational programmes in Central Europe, are compared with corresponding absorption values calculated from electron density profiles of the International Reference Ionosphere (IRI-1979) using the full-wave method.Discrepancies between calculated and observed absorption values were found for the diurnal and the solar-cycle variation, the amplitudes of the solar-cycle variation of absorption being considerably larger than the observed variation.A modification of the solar-activity dependence of the D-region electron density parameters is derived, which provides an improvement of the solar-cycle variation as well as the diurnal variation of the IRI electron density profiles.     

Comparison of observed TEC values with IRI-2007 TEC and IRI-2007 TEC with optional foF2 measurements predictions at an equatorial region,Chumphon, Thailand

In this research, as part of working towards improving the IRI over equatorial region, the total electron content (TEC) derived from GPS measurements and IRI-2007 TEC predictions at Chumphon station (10.72°N, 99.37°E), Thailand, during 2004–2006 is analyzed. The seasonal variation of the IRI-2007 TEC predictions is compared with the TEC from the IRI-2007 TEC model with the option of the actual F2 plasma frequency (foF2) measurements as well as the TEC from the GPS and International GNSS service (IGS). The Chumphon station is located at the equatorial region and the low latitude of 3.22°N. For a declining phase of the solar cycle (2004–2006), the study shows that the IRI-2007 TEC underestimates the IRI-2007 TEC with the foF2 observation at the nighttime by about 5 TECU. The maximum differences are about 15 TECU during daytime and 5 TECU during nighttime. The overestimation is more evident at daytime than at nighttime. When compared in terms of the root-mean square error (RMSE), we find that the highest RMSE between GPS TEC and IRI 2007 TEC is 14.840 TECU at 1230 LT in 2004 and the lowest average between them is 1.318 TECU at 0630 LT in 2006. The noon bite-out phenomena are clearly seen in the IRI-2007 TEC with and without optional foF2 measurements, but not on the GPS TEC and IGS TEC. The IRI TEC with optional foF2 measurements gives the lowest RMSE values between IRI TEC predicted and TEC measurement. However, the TEC measurements (GPS TEC and IGS TEC) are more correct to use at Chumphon station.     

Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

Comparison of peak characteristics of F2 ionospheric layer over Tehran region at a low solar activity period with IRI-2001 and IRI-2007 models predictions

In the present work values of peak electron density (NmF2) and height of F2 ionospheric layer (hmF2) over Tehran region at a low solar activity period are compared with the predictions of the International Reference Ionosphere models (IRI-2001 and IRI-2007). Data measured by a digital ionosonde at the ionospheric station of the Institute of Geophysics, University of Tehran from July 2006 to June 2007 are used to perform the calculations. Formulations proposed by  and  are utilized to calculate the hmF2. The International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) options are employed to run the IRI-2001 and IRI-2007 models. Results show that both IRI-2007 and IRI-2001 can successfully predict the NmF2 and hmF2 over Tehran region. In addition, the study shows that predictions of IRI-2007 model with CCIR coefficient has closer values to the observations. Furthermore, it is found that the monthly average of the percentage deviation between the IRI models predictions and the values of hmF2 and NmF2 parameters are less than 10% and 21%, respectively.     

Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.