基于IRI-2012模型广州地区f_0F_2实测与预测的对比分析

利用广州台站(23.2°N,113.3°E)的实测数据和IRI-2012模型提供的预测数据,对比分析了2013年广州地区f_0F2的变化特征.结果表明,IRI-2012模型能够较好预测该地区f_0F_2的变化趋势,并且CCIR参数得到的预测值比URSI参数更接近实测值;预测值与实测值存在系统偏差,在11:00 LT-06:00 LT时段,观测值均比预测值大,其他时段则相反.在日落后至午夜前时段,预测值与实测值有较大差距.绝对偏差的极值点通常出现在20:00 LT左右,最大超过4 MHz.相对偏差变化比较明显的时段是午夜后至凌晨;在02:00 LT或04:00 LT及06:00 LT附近,可能会出现双误差峰值点,最大超过0.4;但在σ变化很大的20:00 LT附近,相对偏差却变化不大.夜间增强现象会使得偏差增大,导致预测值不能很好反映实测f_0F2的变化.     

Diurnal and seasonal variations of F2 layer characteristics over Irkutsk during the decrease in solar activity in 2003–2006: Observations and IRI-2001 model predictions

The diurnal and seasonal variations of F2 layer characteristics (critical frequency, peak height and bottomside thickness) over Irkutsk, Russia (52.3 N and 104.3 E) are studied by the method of running medians. The comparison with the IRI-2001 model during the decrease in solar activity in 2003–2006 revealed cases of both close agreement and systematic differences between predictions and observations. The systematic difference is not the only reason for disagreement between IRI and observations; there are also intrayear variations which are not associated with seasonal behavior. The period of observation was too short to make conclusions about solar activity dependence of the noon bottomside thickness and the modification of its diurnal behavior with decreasing solar activity.     

Comparison of observed ionospheric vertical TEC over the sea in Indian region with IRI-2016 model

The vertical ionospheric TEC values obtained from GAGAN grid based ionospheric delay correction values over the sea in the Indian equatorial region have been compared with the corresponding values derived from the International Reference Ionosphere model, IRI-2016. The objective of this work is to study the deviation of the vertical TEC derived from the IRI model from ground truths over the sea for different conditions. This will serve the basic intention of assessing the candidature of the IRI model as an alternative ionospheric correction model in navigation receivers in terms of accuracy. We have chosen different solar activity periods, seasons, geomagnetic conditions, locations etc. for our comparison and analysis. The TEC values by the IRI-2016 were compared with the actual measured values for the given conditions and errors were obtained. The measured vertical TEC values at the ionospheric grid points were derived from the GAGAN broadcast ionospheric delay data and used as reference. The IRI model with standard internal functions was used in estimating the TEC at the same ionospheric grid points. The errors in the model derived values are statistically analysed. Broadly, the results show that, for the Indian sector over the sea, the IRI model performs better on quiet days in off equatorial regions, particularly in the northern region. The overall performance degrades for other conditions with the model generally underestimating the true TEC values and most severely in the equatorial region. The performance is worst in this region for the disturbed days of the equinoctial period. The comparison study is also done with the TEC data measured directly by dual frequency GPS receivers. The results were found to be in general agreement with those obtained by comparing the model with GAGAN broadcast data as reference. This study will be useful in considering the IRI-2016 model for real time estimates of TEC as an alternative to the current parametric model in a satellite navigation receiver in absence of other options.     

Performance evaluation of IRI-2007 at equatorial latitudes and its Matlab version for GNSS applications

International Reference Ionosphere (IRI) model is the widely used empirical model for ionospheric predictions, especially TEC which is an important parameter for radio navigation and communication. The Fortran based IRI-2007 does not support real-time interactive visualization and debugging. Therefore, the source code is converted into Matlab and is validated for the purposes of this study. This facilitates easy representation of results and for near real-time implementation of IRI in the applications including spacecraft launching, now casting, pseudolite based navigation systems etc. In addition, the vertical delay results over the equatorial region derived from IRI and GPS data of three IGS stations namely Libreville (Garbon, Africa), Brasilia (Brazil, South America) and Hyderabad (India, Asia) are compared. As the IRI model does not account for plasmasphere TEC, the vertical delays are underestimated compared to vertical delays of GPS signals. Therefore, the model should be modified accordingly for precise TEC estimation.     

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

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

Comparison of GPS TEC variations with IRI-2007 TEC prediction at equatorial latitudes during a low solar activity (2009–2011) phase over the Kenyan region

This work presents an analysis of the Total Electron Content (TEC) derived from the International GNSS Service (IGS) receivers at Malindi (mal2: 2.9^oS, 40.1^oE, dip −26.813^o), Kasarani (rcmn: 36.89^oE, 1.2^oS, dip −23.970^o), Eldoret (moiu: 35.3^oE, 0.3^oN, dip −21.037^o) and GPS-SCINDA (36.8^oE, 1.3^oS, dip −24.117^o) receiver located in Nairobi for the period 2009–2011. The diurnal, monthly and seasonal variations of the GPS derived TEC (GPS-TEC) and effects of space weather on TEC are compared with TEC from the 2007 International Reference Ionosphere model (IRI-TEC) using the NeQuick option for the topside electron density. The diurnal peaks in GPS-TEC is maximum during equinoctial months (March, April, October) and in December and minimum in June solstice months (May, June, July). The variability in GPS-TEC is minimal in all seasons between 0:00 and 04:00 UT and maximum near noon between 10:00 and 14:00 UT. Significant variability in TEC at post sunset hours after 16:00 UT (19:00 LT) has been noted in all the seasons except in June solstice. The TEC variability of the post sunset hours is associated with the occurrence of the ionization anomaly crest which enhances nighttime TEC over this region. A comparison between the GPS-TEC and IRI-TEC indicates that both the model and observation depicts a similar trend in the monthly and seasonal variations. However seasonal averages show that IRI-TEC values are higher than the GPS-TEC. The IRI-TEC also depicts a double peak in diurnal values unlike the GPS-TEC. This overestimation which is primarily during daytime hours could be due to the model overestimation of the equatorial anomaly effect on levels of ionospheric ionization over the low latitude regions. The IRI-TEC also does not show any response to geomagnetic activity, despite the STORM option being selected in the model; the IRI model generally remains smooth and underestimates TEC during a storm. The GPS-TEC variability indicated by standard deviation seasonal averages has been presented as a basis for extending the IRI-model to accommodate TEC-variability.     

Accurate estimation of relative atmospheric density error on the example of uncertain geomagnetic activity information

Errors in neutral atmospheric density are the dominant contributor to unrealistic orbital state-vector covariances in low Earth orbits (LEO). Density uncertainty is caused by model-uncertainty at spatial scales below and within the model resolution, as well as input-uncertainty of the environmental parameters supplied to the semi-empirical atmospheric model.The paper at hand provides multiple contributions. First, analytic equations are derived to estimate the relative density error due to an input parameter uncertainty in any of the environmental parameters supplied to the model. Second, it is shown on the example of uncertain geomagnetic activity information, how to compute the required inputs to facilitate the accurate estimation of the relative density error.A clamped cubic splining approach for the conversion from geomagnetic amplitude (ap) to the kp index is postulated to perform this uncertainty propagation, as other algorithms were found unsuitable for this task. Results of numerical simulations with three popular semi-empirical models are provided to validate the set of derived equations. It is found that geomagnetic input uncertainty is especially important to consider in case of low global geomagnetic activity. The findings seamlessly integrate with prior work by the authors to perform density-uncertainty considering orbit estimation.     

On the variabilities of the Total Electron Content (TEC) over the Indian low latitude sector

The Total Electron Content (TEC) from four locations in the Indian sector namely, Trivandrum (8.47°N, 76.91°E, Geomag.0.63°S, 0.3° dip), Waltair (17.7° N, 83.3°E, Geomag. 6.4°N, 20° dip), Bhopal (23.28°N, 77.34°E, Geomag.14.26°N, 33.2° dip), and Delhi (28.58°N, 77.21°E, Geomag.19.2°N, 43.4° dip) during a low sunspot year of 2004 are used to study the variabilities of the TEC. The day time TEC values are higher over Waltair and Bhopal compared to those at Trivandrum and Delhi. Considerable day-to-day variations in the diurnal values of TEC are observed at the anomaly crest locations. The observed GPS-TEC has been compared with the IRI-2007 model derived TEC considering three different options (IRI-2001, IRI-2001 corrected and Ne-Quick) available in the model for the topside electron density. The TEC derived with Ne-Quick and IRI-01 corrected options show better agreement with GPS-TEC while the TEC from IRI-01 method shows larger deviations. From the correlation analysis carried out between TEC value at 1300 h LT and solar indices parameters namely sunspot number (SSN), F10.7 and EUV, it is observed that the correlation is more during equinoctial months and less during summer months. The correlation coefficients observed over the anomaly locations, Bhopal and Delhi are lower compared to those at Trivandrum and Waltair.     

The subauroral electron density trough: Comparison between satellite observations and IRI-2007 model estimates

We compare electron density predictions of the International Reference Ionosphere (IRI-2007) model with in situ measurements of the satellites CHAMP and GRACE for the years from 2005 to 2010 over the subauroral regions. The electron density between 58° and 68° Mlat are considered. The trough region Ne peaks during local summers and attain the valley during local winter. Around −100°E and 60°E, two larger electron density sectors features can be seen in both hemispheres during all three seasons, which attributed to the electron extending from middle latitude to trough region. From 2005 to the beginning of 2010, the model overestimates the trough region Ne by 20% on average and the decrease of Ne in this region can also be seen during the last solar minimum. In the southern hemisphere, the model prediction shows quite well consistence with the observation during all three seasons while the huge difference between observations and model estimation implies that the IRI-2007 model needs significant improvement to predict better the trough region in northern hemisphere.     

A topside investigation over a mid-latitude digisonde station in Cyprus

We examine the systematic differences between topside electron density measurements recorded by different techniques over the low-middle latitude operating European station in Nicosia, Cyprus (geographical coordinates: 35.14^oN, 33.2^oE), (magnetic coordinates 31.86^oN, 111.83 ^oE). These techniques include space-based in-situ data by Langmuir probes on board.European Space Agency (ESA) Swarm satellites, radio occultation measurements on board low Earth orbit (LEO) satellites from the COSMIC/FORMOSAT-3 mission and ground-based extrapolated topside electron density profiles from manually scaled ionograms. The measurements are also compared with International Reference Ionosphere Model (IRI-2016) topside estimations and IRI-corrected NeQuick topside formulation (method proposed by Pezzopane and Pignalberi (2019)). The comparison of Swarm and COSMIC observations with digisonde and IRI estimations verifies that in the majority of cases digisonde underestimates while IRI overestimates Swarm observations but in general, IRI provides a better topside representation than the digisonde. For COSMIC and digisonde profiles matched at the F layer peak the digisonde systematically underestimates topside COSMIC electron density values and the relative difference between COSMIC and digisonde increases with altitude (above hmF2), while IRI overestimates the topside COSMIC electron density but after a certain altitude (~150 km above hmF2) this overestimation starts to decrease with altitude. The IRI-corrected NeQuick underestimates the majority of topside COSMIC electron density profiles and relative difference is lower up to approximately 100 km (above the hmF2) and then it increases. The overall performance of IRI-corrected NeQuick improves with respect to IRI and digisonde.