Comparison of electron density measurements from CSES and Swarm satellites with GNSS ionospheric tomography data

Electron density measurements obtained from China Seismo‐Electromagnetic Satellite (CSES) and Swarm-B can play an increasingly important role in the study of ionosphere above F2 peak height. This study presented a comprehensive comparison of electron density products obtained from Langmuir probe mounted on CSES and Swarm-B with ionospheric tomography for a whole year period of 2019. CSES was fully compared with Swarm-B on a global scale, including both absolute and relative differences, and a new index called NFI was developed to better quantify the similarity between two latitudinal profiles of electron density. CSES and Swarm-B were then compared with tomography respectively in four regions, roughly located in America, Europe, Australia and China. Results indicated that CSES data are consistent with Swarm-B, as NFI values exceed 0.6 for most of the analyses. Tomography and Swarm-B were found to have a good agreement as their biases are less than 0.2 × 10⁵ el/cm³ in general. For the comparison between CSES and tomography, the bias increased to around 0.6 × 10⁵ el/cm³ but the standard deviation changed slightly, validating the underestimation of electron density by CSES. The spatiotemporal comparisons of CSES and Swarm-B with tomography showed that: 1) the differences in electron density were relatively low in middle latitudes and increased rapidly in the regions of equatorial ionization anomaly; 2) Swarm-B has a better consistent with tomography than CSES, but both are capable of detecting ionosphere anomalies such as midlatitude arcs; and 3) CSES and Swarm-B both can capture the seasonal changes of electron density, while their values are basically smaller than those from tomography in Spring and Summer months.     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

电离层效应的测时误差及其修正

电磁波经过电离层传播时会受到电离层折射的影响而产生延迟, 星载接收机探测到的时间是信号延迟之后的到达时间. 某次实验数据显示, 一些波段的瞬态电测辐射信号的群时延之差可达10⁵ns数量级, 这在对信号源进行时差定位时是不能直接运用的. 为有效消除电离层延迟的影响, 将双频修正法应用于某项工程中, 利用接收到的实验数据求解出电离层的TEC(电离层总电子含量), 并在此基础上对信号到达星载接收机的时间进行修正. 最后,对修正结果进行了验证, 给出了误差来源.     

Tropospheric wet refractivity tomography using multiplicative algebraic reconstruction technique

Algebraic reconstruction techniques (ART) have been successfully used to reconstruct the total electron content (TEC) of the ionosphere and in recent years be tentatively used in tropospheric wet refractivity and water vapor tomography in the ground-based GNSS technology. The previous research on ART used in tropospheric water vapor tomography focused on the convergence and relaxation parameters for various algebraic reconstruction techniques and rarely discussed the impact of Gaussian constraints and initial field on the iteration results. The existing accuracy evaluation parameters calculated from slant wet delay can only evaluate the resultant precision of the voxels penetrated by slant paths and cannot evaluate that of the voxels not penetrated by any slant path. The paper proposes two new statistical parameters Bias and RMS, calculated from wet refractivity of the total voxels, to improve the deficiencies of existing evaluation parameters and then discusses the effect of the Gaussian constraints and initial field on the convergence and tomography results in multiplicative algebraic reconstruction technique (MART) to reconstruct the 4D tropospheric wet refractivity field using simulation method.     

GNSS radio tomography of the ionosphere: The problem with essentially incomplete data

We discuss the specific features of the spatiotemporal radio tomography of ionospheric electron density based on the data from high orbiting global navigational satellite system (GNSS). The main peculiarities are four-dimensionality of the problem and essential incompleteness of input data. The approach suggested for the solution of this problem is based on the search for the solution with minimal Sobolev’s norm (the most smooth solution), which provides, in particular, smooth extension of the solution to the area where no input data is available. Methods and algorithms for the solution of the problem are developed, and the resolution of the proposed technique is estimated. Examples of global and regional radio tomographic imaging based on GNSS data are presented. The obtained results are compared with the data of ionosondes and with the results of radio tomography based on the radio transmissions from low-orbiting satellites.     

An overview of radar soundings of the martian ionosphere from the Mars Express spacecraft

The Mars Express spacecraft carries a low-frequency radar called MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) that is designed to study the subsurface and ionosphere of Mars. In this paper, we give an overview of the ionospheric sounding results after approximately one year of operation in orbit around Mars. Several types of ionospheric echoes are commonly observed. These include vertical echoes caused by specular reflection from the horizontally stratified ionosphere; echoes from a second layer in the topside ionosphere, possibly associated with O⁺ ions; oblique echoes from upward bulges in the ionosphere; and a variety of other echoes that are poorly understood. The vertical echoes provide electron density profiles that are in reasonable agreement with the Chapman photo-equilibrium model of planetary ionospheres. On the dayside of Mars the maximum electron density is approximately 2 × 10⁵ cm⁻³. On the nightside the echoes are often very diffuse and highly irregular, with maximum electron densities less than 10⁴ cm⁻³. Surface reflections are sometimes observed in the same frequency range as the diffuse echoes, suggesting that small isolated holes exist in the nightside ionosphere, possibly similar to those that occur on the nightside of Venus. The oblique echoes arise from upward bulges in the ionosphere in regions where the crustal magnetic field of Mars is strong and nearly vertical. The bulges tend to be elongated in the horizontal direction and located in regions between oppositely directed arch-like structures in the crustal magnetic field. The nearly vertical magnetic field lines in the region between the arches are thought to connect into the solar wind, thereby allowing solar wind electrons to heat the lower levels of the ionosphere, with an attendant increase in the scale height and electron density.     

The mean doppler velocity derived by an advanced ionospheric sounder — dynasonde

The lower ionosphere was examined by using a fixed sounding frequency 7.953 MHz of the EISCAT dynasonde. The phase parameter and Doppler spectrum techniques were applied on the echoes to obtain radial Doppler velocities. It is found that the two derived velocities are very different. A procedure is developed and introduced to apply on the Doppler velocities obtained from the phase parameter technique. Result shows that mean Doppler velocities derived by the weighted mean procedure are close to those of the typical spectrum technique.     

Fading in the HF ionospheric channel and the role of irregularities

It is well known that the ionosphere affects radio wave propagation especially in the high frequency (HF) range. HF radio waves reflected by the ionosphere can reach considerable distances, often with changes in amplitude, phase, and frequency. The ionosphere is a dispersive in frequency and time, bi-refractive, absorbing medium, in which multipath propagation due to traveling irregularities is very frequent. The traveling irregularities undulate the reflecting ionospheric layer, introducing variations in signal amplitude (fading). In this multipath time variant channel fading is mainly considered, even though it is not the sole effect. Echo signals from a single reflection, as in ionospheric vertical sounding (VIS) techniques, are affected by a certain degree of variability even in quiet ionospheric conditions. In this work the behavior of the ionospheric channel is studied and characterized by observing the power variation of received echoes using the VIS technique. Multipath fading was analyzed quantifying the power variation of the signal echo due to irregularities on a temporal scale from 0.5 to 256 s. An experimental set-up derived from an ionosonde was implemented and the analysis was performed employing a special numerical algorithm operating off-line on the acquired time sequence of the signal. The gain-loss of the irregularity shapes are determined in some special cases.     

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.     

Model based Computerized Ionospheric Tomography in space and time

Reconstruction of the ionospheric electron density distribution in space and time not only provide basis for better understanding the physical nature of the ionosphere, but also provide improvements in various applications including HF communication. Recently developed IONOLAB-CIT technique provides physically admissible 3D model of the ionosphere by using both Slant Total Electron Content (STEC) measurements obtained from a GPS satellite - receiver network and IRI-Plas model. IONOLAB-CIT technique optimizes IRI-Plas model parameters in the region of interest such that the synthetic STEC computations obtained from the IRI-Plas model are in accordance with the actual STEC measurements. In this work, the IONOLAB-CIT technique is extended to provide reconstructions both in space and time. This extension exploits the temporal continuity of the ionosphere to provide more reliable reconstructions with a reduced computational load. The proposed 4D-IONOLAB-CIT technique is validated on real measurement data obtained from TNPGN-Active GPS receiver network in Turkey.     

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.     

利用COSMIC掩星数据监测电离层的异常变化

分析了COSMIC掩星数据反演电子密度的方法,利用实例研究反演方法的特点,并采用ISR非相干散射雷达获取的电子密度数据进行验证,进而反演了长三角区域SHAO（IGS）站上空在日全食和太阳风暴期间的电子密度廓线图. 通过与平静日期间电离层电子密度进行比较,发现日全食及太阳风暴导致电离层发生的异常变化,从而提出COSMIC掩星数据反演电子密度在监测电离层变化时所具有的优势.     

南京地区斜向加热电离层效应数值模拟

与垂直加热相比,斜向加热电离层具有灵活性高、影响范围广和便于实际操作等优势.在非偏区考虑欧姆吸收,基于电子能量方程和连续性方程构建电波斜向加热低电离层的物理模型,并利用IRI-2007和NRLMSISE-00经验模型提供的背景参数对南京地区斜向加热低电离层进行数值模拟,对比不同加热条件下电子温度和电子密度的扰动情况.研究结果表明：电波加热效果随入射仰角和有效辐射功率的增大而增大;电子温度和电子密度增幅随电波频率增大而减小;X波模比O波模造成的电子温度扰动幅度和电子密度扰动幅度更大,同时X波模比O波模能更快地使电子温度和电子密度达到稳定状态;一定范围内较小仰角、较低频率、较大有效辐射功率的电波能使电子密度更快达到稳定,后两者还能加快电子温度达到稳定的过程;电子温度达到稳定所需时间随入射电波仰角呈单峰变化,仰角为62°时达到最大.

     

Tomographic reconstruction of the ionospheric electron density as a function of space and time

Electron density distribution is the major determining parameter of the ionosphere. Computerized Ionospheric Tomography (CIT) is a method to reconstruct ionospheric electron density image by computing Total Electron Content (TEC) values from the recorded Global Positioning Satellite System (GPS) signals. Due to the multi-scale variability of the ionosphere and inherent biases and errors in the computation of TEC, CIT constitutes an underdetermined ill-posed inverse problem. In this study, a novel Singular Value Decomposition (SVD) based CIT reconstruction technique is proposed for the imaging of electron density in both space (latitude, longitude, altitude) and time. The underlying model is obtained from International Reference Ionosphere (IRI) and the necessary measurements are obtained from earth based and satellite based GPS recordings. Based on the IRI-2007 model, a basis is formed by SVD for the required location and the time of interest. Selecting the first few basis vectors corresponding to the most significant singular values, the 3-D CIT is formulated as a weighted least squares estimation problem of the basis coefficients. By providing significant regularization to the tomographic inversion problem with limited projections, the proposed technique provides robust and reliable 3-D reconstructions of ionospheric electron density.     

Improved predictions of equatorial TEC for real time applications

It is important to use models developed specifically for the equatorial ionospheric estimation for real-time applications, particularly in Satellite Navigation. This work demonstrates a methodology for improved predictions of VTEC in real time using the model developed for the equatorial ionosphere by the authors. This work has been done using TEC data of the low solar activity period of 2005 obtained using dual frequency GPS receivers installed under the GAGAN project of ISRO. For the purpose, the model is first used in conjunction with Kriging technique. Improvement in accuracy is observed when compared with the estimations from the model alone using the measurements as true reference. Further improvement is obtained by Bayesian combination of these estimates with independent Neural Network based predictions. Statistical performance of improvement is provided. An improvement of ∼1 m in confidence level of estimation of VTEC is obtained.     

星载SAR电离层探测及误差补偿

工作于低频波段的星载合成孔径雷达(SAR)信号会受到电离层的显著影响,因此在系统设计时必须考虑相应的补偿方法。基于此,首先开展了适用于低频大宽带SAR模式的电离层影响评估研究,建立了基于勒让德展开的五阶误差分析模型,可有效解决传统模型各阶次耦合问题。其次,针对背景电离层色散问题,开展了基于双频自聚焦算法的补偿研究,利用ALOS PALSAR数据进行了半物理仿真验证,电离层反演精度优于0.4TECU,可有效提升图像聚焦质量;针对法拉第旋转角误差,利用ALOS PALSAR回波散射矩阵信息开展了补偿研究,结果显示,相比官网提供的补偿参数,误差可进一步降低27%。上述基于回波数据本身的补偿研究,可避免第三方数据精度差、分辨率低、需地面接收机、传播路径不一致等问题,同时降低了载荷成本。最后,结合SAR高分辨率特性,基于回波的电离层反演可有效提高现有电离层探测能力,开展了PALSAR-垂测仪联合反演电子密度研究,其结果比仅垂测仪数据精度普遍提高了30%以上。研究成果可为未来低频星载SAR的系统设计提供技术支撑。     

Multi-GNSS contributions to differential code biases determination and regional ionospheric modeling in China

Due to the limited number and uneven distribution globally of Beidou Satellite System (BDS) stations, the contributions of BDS to global ionosphere modeling is still not significant. In order to give a more realistic evaluation of the ability for BDS in ionosphere monitoring and multi-GNSS contributions to the performance of Differential Code Biases (DCBs) determination and ionosphere modeling, we select 22 stations from Crustal Movement Observation Network of China (CMONOC) to assess the result of regional ionospheric model and DCBs estimates over China where the visible satellites and monitoring stations for BDS are comparable to those of GPS/GLONASS. Note that all the 22 stations can track the dual- and triple-frequency GPS, GLONASS, and BDS observations. In this study, seven solutions, i.e., GPS-only (G), GLONASS-only (R), BDS-only (C), GPS + BDS (GC), GPS + GLONASS (GR), GLONASS + BDS (RC), GPS + GLONASS + BDS (GRC), are used to test the regional ionosphere modeling over the experimental area. Moreover, the performances of them using single-frequency precise point positioning (SF-PPP) method are presented. The experimental results indicate that BDS has the same ionospheric monitoring capability as GPS and GLONASS. Meanwhile, multi-GNSS observations can significantly improve the accuracy of the regional ionospheric models compared with that of GPS-only or GLONASS-only or BDS-only, especially over the edge of the tested region which the accuracy of the model is improved by reducing the RMS of the maximum differences from 5–15 to 2–3 TECu. For satellite DCBs estimates of different systems, the accuracy of them can be improved significantly after combining different system observations, which is improved by reducing the STD of GPS satellite DCB from 0.243 to 0.213, 0.172, and 0.165 ns after adding R, C, and RC observations respectively, with an increment of about 12.3%, 29.4%, and 32.2%. The STD of GLONASS satellite DCB improved from 0.353 to 0.304, 0.271, and 0.243 ns after adding G, C, and GC observations, respectively. The STD of BDS satellite DCB reduced from 0.265 to 0.237, 0.237 and 0.229 ns with the addition of G, R and GR systems respectively, and increased by 10.6%, 10.4%, and 13.6%. From the experimental positioning result, it can be seen that the regional ionospheric models with multi-GNSS observations are better than that with a single satellite system model.     

Precision analysis of troposphere sensing using GPS single-frequency signals

Various studies have been performed to investigate the accuracy of troposphere zenith wet delays (ZWDs) determined from GPS. Most of these studies use dual-frequency GPS data of large-scale networks with long baselines to determine the absolute ZWDs. For small-scale networks the estimability of the absolute ZWDs deteriorates due to high correlation between the solutions of the ZWDs and satellite-specific parameters as satellite clocks. However, as relative ZWDs (rZWDs) can always be estimated, irrespective of the size of the network, it is of interest to understand how the large-scale network rZWD-performance of dual-frequency GPS using an ionosphere-float model compares to the small-scale network rZWD-performance of single-frequency GPS using an ionosphere-weighted model. In this contribution such an analysis is performed using undifferenced and uncombined network parametrization modelling. In this context we demonstrate the ionosphere weighted constraints, which allows the determination of the rZWDs independent from signals on the second frequency. Based on an analysis of both simulated and real data, it is found that under quiet ionosphere conditions, the accuracy of the single-frequency determined rZWDs in the ionosphere-weighted network is comparable to that of the large-scale dual-frequency network without ionospheric constraints. Making use of the real data from two baselines of 15?days, it was found that the absolute differences of the rZWDs applying the two strategies are within 1?cm in over 90% and 95% of the time for ambiguity-float and -fixed cases, respectively.     

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.