基于小波分解与重构方法研究电离层偶发E层

电离层不规则体对卫星导航、通信、雷达系统等有重要影响.通过数值模拟及与实测数据的对比,论证基于小波分解与重构方法实现利用掩星数据反演电离层不规则体的可行性.以电离层偶发E层为例,利用国际参考电离层(IRI)模拟背景电离层电子密度分布,利用掩星探测的水平电子密度总含量δht反演不规则体信息,并与模拟数据进行比较.对200...     

电离层掩星反演的TEC修正反演法

电离层掩星数据反演的传统方法是采用改正TEC的Abel 变换反演法, 实际电离层的非球对称性会给电子密度的反演结果带来误差. 文中研究了利用TEC修正方法结合背景场来剔除TEC 受电子密度水平变化的影响, 改善球对称假设适用性, 提高反演精度, 并应用此方法于模拟掩星事件的反演. 结果表明, 与传统的Abel 变换反演相比, TEC 修正反演法能够减小反演误差. 用TEC 修正反演法对不同方法获取的背景场的反演结果比较表明, 背景场与实际场吻合的程度越高, 反演效果越好.      

TIEGCM集合卡尔曼滤波同化模型设计及初步试验

选择参数化的电离层热层理论模型TIEGCM作为背景模型,基于COSMIC掩星观测的电子密度廓线数据,应用集合卡尔曼滤波方法建立全球电离层电子密度同化模型,实现了全球电离层的电子密度同化.同化结果表明,该同化模型能将观测资料有效同化到背景模式中,获得全球三维电离层电子密度.与背景模式相比,同化得到的电子密度相对于观测值的偏差显著下降.对于有同化和无同化参与的试验,N_mF₂的标准偏差分别降低约60%和20%.此外,分组同化与同时同化的结果对比显示,平均偏差改善基本一致,同时同化后的标准偏差在峰值高度以上略有减小.      

利用COSMIC数据构建的全球垂直TEC图的验证研究

在日固坐标系(地磁纬度和地方时)下, 累积地方时过去24h的COSMIC(Constellation Observing System for Meteorology Ionosphere andClimate)观测资料, 通过对110$sim$750km高度范围内的电子密度进行数值积分得到各掩星点的垂直TEC值, 进而利用Kriging方法插值产生近实时的全球地方时MAGLat2.5°×2h的COSMIC TEC图. 利用2008年1月1日至2010年6月30日共30个月的COSMIC数据, 逐日构建COSMICTEC图, 将其与全球导航卫星系统服务组织(International GNSS Service,IGS)发布的全球电离层TEC图(Global Ionospheric Maps, GIMs)以及OSTM/JASON-2卫星高度计观测值分别进行比对,证明利用COSMIC掩星资料构建全球电离层垂直TEC图是可行的.      

风云三号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掩星探测资料的综合利用具有重大意义.      

基于GPS掩星平均弯曲角的折射率反演研究

利用Abel积分变换通过掩星弯曲角计算折射率需要对高层弯曲角进行统计优化.目前由于所使用的背景场资料和具体反演方法不同,导致所发布的掩星大气数据气候统计值存在一定差异.本文使用2008年1,4,7,10月共4个月的COSMIC大气掩星附加相位数据,从纬圈平均弯曲角廓线反演相应月平均折射率,对反演结果进行比较分析.研究表明,利用掩星折射数据进行气候研究时无需逐一对掩星探测廓线进行统计优化,在40 km以下高度基于平均弯曲角的反演方法与传统统计相比能够获得几乎一致的月平均折射率,在50 km以上高度基于平均弯曲角的反演结果更加接近ECMWF资料统计.      

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

电离层峰值高度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模型的完善提供了一定参考.      

基于电子密度同化的热层大气密度预报方法

采用热层电离层耦合模式TIEGCM和集合卡尔曼滤波同化方法,利用同化COSMIC电离层掩星电子密度数据优化热层电离层参量,并将模式预报的大气密度与CHAMP卫星大气密度数据进行对比,分别开展模拟和实测数据的同化预报实验.在模拟数据同化实验中,状态向量包含温度、风场和离子成分的实验结果表明,仅优化温度即可达到最优的热层大气密度预报效果.在实测数据同化实验中,将温度作为状态向量参数,优化结果表明,循环同化过程中模式预报的大气密度相对偏差的均方根误差在48h内从38%减小到27%,同化稳定时间至少需要30h.预报过程中大气密度预报效果的改善持续时间为34h.这表明电子密度同化能够改善热层大气密度的预报精度,设计的实验方案合理可行,可获得较长的预报时效.      

基于STK的GPS/LEO无线电掩星技术仿真研究

GPS/LEO无线电掩星技术反演地球大气参数剖面已经具有较高的精度. 国外开展了多个GPS/LEO掩星项目, 但中国还尚未深入进行相关的实验, 这制约了中国掩星技术的发展. 本文提出基于STK进行GPS/LEO掩星技术研究的方法; 根据GPS/LEO掩星的原理, 推导出掩星事件发生的条件和掩星切点的计算公式; 利用STK对掩星过程进行模拟, 得到掩星数据. 在大气球对称假设和大气模型已知的条件下, 反演得到中性大气折射指数. 通过比较模型和反演数据, 表明反演数据精度较高, 验证了利用STK模拟GPS/LEO掩星实验方法的可行性.      

ROPP反演软件算法及其精度分析

介绍了ROPP反演软件中无线电掩星反演的算法与精度分析. 采用COSMIC卫星2008年1 月1日全天的附加相位数据, 反演得到折射率、温度、压强与湿度等参数, 并与CDAAC 相应结果进行对比. 实验结果表明, 在30km高度以下, 折射率、压强和湿度的相对 误差在2%以内, 温度误差不超过2K.      

Ionospheric responses to the 21 August 2017 great American solar eclipse – A multi-instrument study

Herein, we report on the ionospheric responses to a total solar eclipse that occurred on 21 August 2017 over the US region. Ground-based GPS total electron content (TEC) data along with ground-based measurements (Millstone Hill Observatory (MHO) and digital ionosondes) and space-based measurements (COSMIC radio occultation (RO) technique) allowed us to identify eclipse-associated ionospheric responses. TEC data at ~20°, ~30°, and ~40°N latitudes from the west to east longitudes show not only considerable depression but also wave-like characteristics in TEC both in the path of totality and away from it, exclusively on the day of eclipse. Interestingly, the observed depressions are associated with lesser (higher) magnitudes at stations over which the solar obscuration percentage was meager (significant), a clear indication of bow-wave-like features. The MHO observes a 30% reduction in F₂-layer electron densities between 180 and 220 km on eclipse day. Ionosonde-scaled parameters over Boulder (40.4°N, 100°E) and Austin (30.4°N, 94.4°E) show a significant decrease in critical frequencies while an altitude elevation is seen in the virtual heights of the F-layer only during the eclipse day and that decreases are associated with wave-like signatures, which could be attributed to eclipse-generated waves. The estimated vertical electron density profile from the COSMIC RO-based technique shows a maximum depletion of 40%. Relatively intense and moderate depths of TEC depression, considerable reductions in the F₂-layer electron densities measured by the MHO and COSMIC RO-measured densities at the F₂-layer peak, and elevations in virtual heights and reduction in the critical frequencies measured by ionosondes during the eclipse day could be due to the eclipse-induced dynamical effects such as gravity waves (GWs) and their associated electro-dynamical effects (modification of ionospheric electric fields due to GWs).     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

Validation of FORMOSAT-3/COSMIC radio occultation electron density profiles by incoherent scatter radar data

Low Earth Orbiting satellites carrying a dual frequency GPS receiver onboard offer a unique opportunity to remote sensing of the global ionosphere on a continuous basis. No other profiling technique unifies profiling through the entire F2-layer with global coverage. The FORMOSAT-3/COSMIC data can make a positive impact on the global ionosphere study providing essential information about the height electron density distribution and particularly over regions that are not accessible with ground-based measuring instruments such as ionosondes and GPS dual frequency receivers. Therefore, it is important to verify occultation profiles with other techniques and to obtain experience in the reliability of their derivation. In the given study we present results of comparison of the electron density profiles derived from radio occultation measurements on-board FS-3/COSMIC and from the Kharkov incoherent scatter radar sounding.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Electron density retrieval from occulting GNSS signals using a gradient-aided inversion technique

In the coming years, opportunities for remote sensing of electron density in the Earth’s ionosphere will expand with the advent of Galileo, which will become part of the global navigation satellite system (GNSS). Methods for accurate electron density retrieval from radio occultation data continue to improve. We describe a new method of electron density retrieval using total electron content measurements obtained in low Earth orbit. This method can be applied to data from dual-frequency receivers tracking the GPS or Galileo transmitters. This simulation study demonstrates that the method significantly improves retrieval accuracy compared to the standard Abel inversion approach that assumes a spherically symmetric ionosphere. Our method incorporates horizontal gradient information available from global maps of Total Electron Content (TEC), which are available from the International GNSS Service (IGS) on a routine basis. The combination of ground and space measurements allows us to improve the accuracy of electron density profiles near the occultation tangent point in the E and F regions of the ionosphere.     

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).     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

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.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.