电离层垂直TEC映射函数的实验观测与统计特性

利用GPS信标测量获得的电离层电子浓度总含量（TEC）是沿电波路径的斜向TEC.理论研究和实际应用中,常常需要通过映射函数将斜向TEC转换为垂直方向的TEC,这在当前主要采用对电子浓度分布模型的数值积分得到模型映射函数来实现.本文在考察现有不同模型映射函数的基础上,又提出了一种源于实际观测的实验映射函数的概念与估算方法.我们利用IGS的全球GPS观测站的斜向TEC和JPL提供的垂直TEC数据获得了2006年期间的实验映射函数,并对所得结果进行了初步统计分析.在卫星天顶角较小时,上述实验映射函数和模型映射函数之间相差甚微,均可很好描述垂直TEC与斜TEC之间关系;但卫星天顶角较大时,实验映射函数和常用的模型映射函数之间存在明显差异.本文认为,这种差异主要是因为现有模型映射函数中没有考虑到等离子体层的贡献.我们认为采用基于实验映射函数的模式,或者通过考虑等离子体层的贡献对现有模型映射函数进行改进,可以有效提高电离层TEC的估算精度.     

Reconstructing Regional Ionospheric Electron Density: A Combined Spherical Slepian Function and Empirical Orthogonal Function Approach

The computerized ionospheric tomography is a method for imaging the Earth’s ionosphere using a sounding technique and computing the slant total electron content (STEC) values from data of the global positioning system (GPS). The most common approach for ionospheric tomography is the voxel-based model, in which (1) the ionosphere is divided into voxels, (2) the STEC is then measured along (many) satellite signal paths, and finally (3) an inversion procedure is applied to reconstruct the electron density distribution of the ionosphere. In this study, a computationally efficient approach is introduced, which improves the inversion procedure of step 3. Our proposed method combines the empirical orthogonal function and the spherical Slepian base functions to describe the vertical and horizontal distribution of electron density, respectively. Thus, it can be applied on regional and global case studies. Numerical application is demonstrated using the ground-based GPS data over South America. Our results are validated against ionospheric tomography obtained from the constellation observing system for meteorology, ionosphere, and climate (COSMIC) observations and the global ionosphere map estimated by international centers, as well as by comparison with STEC derived from independent GPS stations. Using the proposed approach, we find that while using 30 GPS measurements in South America, one can achieve comparable accuracy with those from COSMIC data within the reported accuracy (1 × 10¹¹ el/cm³) of the product. Comparisons with real observations of two GPS stations indicate an absolute difference is less than 2 TECU (where 1 total electron content unit, TECU, is 10¹⁶ electrons/m²).     

Verification of the concept of seismoionospheric coupling under quiet heliogeomagnetic conditions,using the Wenchuan (China) earthquake of May 12, 2008, as an example

The variability of the ionosphere during April–May, 2008, has been analyzed in detail in order to reveal anomalous variations related to seismic activity, initiated by the strongest Wenchuan earthquake (M = 7.9) in the Sichuan province on May 12, 2008. Information about the total electron content (TEC) from the network of GPS receivers in the earthquake region, the global IONEX TEC maps, and the reconstructed vertical profiles of electron density according to the data of GPS receivers were used as a data source. The spatial and time localization of the observed anomalies, their morphological features, and the absence of geomagnetic disturbances during the observation period undoubtedly demonstrate that the observed variations were caused by seismic activity.     

Performance of different TEC models to provide GPS ionospheric corrections

The existence of a worldwide international GPS service (IGS) permanent network of dual-frequency receivers makes the computation of global ionospheric maps (GIMs) of total electron content (TEC) feasible. The GIMs computed by the IGS Associate Analysis Centers on a daily basis and by other kinds of forecast GIMs, which can be computed from, for instance, the international reference ionosphere (IRI) model, and the GPS broadcast models in the navigation message, can be applied to a broad diversity of fields, for instance as, navigation and time transfer.In this context, the performance of different kinds of models are presented in order to determine the accuracy of the different GIM. This is carried out by comparison with the TOPEX data that provides an independent and precise (at the level of few TECU) vertical TEC determination over the oceans and seas. Thus, the obtained accuracies, in terms of global relative error, ranging from 54% corresponding to the GPS broadcast model, to about 41% corresponding to IRI climatological model, and to less than 30% corresponding to GPS data driven models.     

Automated daily process for global ionospheric total electron content maps and satellite ocean altimeter ionospheric calibration based on Global Positioning System data

The accuracy of single-frequency ocean altimeters benefits from calibration of the total electron content (TEC) of the ionosphere below the satellite. Data from a global network of Global Positioning System (GPS) receivers provides timely, continuous, and globally well-distributed measurements of ionospheric electron content. For several months we have been running a daily automatic Global Ionospheric Map process which inputs global GPS data and climatological ionosphere data into a Kalman filter, and produces global ionospheric TEC maps and ocean altimeter calibration data within 24 h of the end-of-day. Other groups have successfully applied this output to altimeter data from the GFO satellite and in orbit determination for the TOPEX/Poseidon satellite. Daily comparison of the global TEC maps with independent TEC data from the TOPEX altimeter is performed as a check on the calibration whenever the TOPEX data are available. Comparisons of the global TEC maps against TOPEX data will be discussed. Accuracy is best at mid-to-high absolute latitudes (∣latitude∣>30°) due to the better geographic distribution of GPS receivers and the relative simplicity of the ionosphere. Our highly data-driven technique is relatively less accurate at low latitudes and especially during ionospheric storm periods, due to the relative scarcity of GPS receivers and the structure and volatility of the ionosphere. However, it is still significantly more accurate than climatological models.     

Global plasmaspheric TEC and its relative contribution to GPS TEC

The plasmaspheric electron content is directly estimated from the global positioning system (GPS) data onboard JASON-1 Satellite for the first time. Similarly, the ground-based GPS total electron content (TEC) is estimated using about 1000 GPS receivers distributed around the globe. The relative contribution of the plasmaspheric electron content to the ground-based GPS TEC is then estimated globally using these two independent simultaneous measurements; namely ground-based GPS TEC and JASON-1 GPS TEC. Results presented here include data from 3 months of different solar cycle conditions (October 2003, May 2005, and December 2006). The global comparison between the two independent measurements was performed by dividing the data into three different regions; equatorial, mid- and high-latitude regions. This division is essential as the GPS raypaths traverse different distances through the plasmasphere at different latitudes. The raypath length through the plasmasphere decreases as latitude increases. The relative contribution of the plasmaspheric electron content exhibits a diurnal variation that depends on latitude with minimum contribution (10%) during daytime and maximum (up to 60%) at night. The contribution is also maximum at the equatorial region where the GPS raypath traverses a long distance through the plasmasphere compared to its length in mid- and high-latitude regions. Finally, the solar cycle variation of plasmaspheric contribution is also reported globally.     

基于等离子体GCPM模型对电离层薄壳模型高度的仿真研究

在基于GPS数据提取电离层总电子含量(TEC)的过程中,电离层薄壳高度的选择对解算电离层垂直TEC的精度有很大的影响.但由于不可能获得一个真实的从电离层D层到GPS卫星高度的电子密度剖面,关于电离层薄壳高度的选择一直是基于GPS数据解算电离层TEC方法中关注的一个问题.本文利用等离子体GCPM模型,对太阳活动高年(2002)和太阳活动低年(2008)情况下电离层有效薄壳高度的选择进行了仿真计算.结果表明,最佳的薄壳高度在2002年为560 km,而在2008年为695 km.通过对全球八个具有代表性地点的仿真计算,揭示了有效薄壳高度更复杂的变化特点.在白天,最佳薄壳的高度变化不大(500 km至750 km);但在夜晚,最佳薄壳高度变化范围很大,甚至可以超过2000 km.此外,本文还对不同卫星仰角的情况下斜向TEC转换为垂直TEC的误差进行了分析,结果表明:随着卫星仰角的增加,薄壳模型带来的转换误差基本上是单调减少的.因而,在实际应用中,尽可能地采用大仰角的卫星数据有助于提高解算的电离层垂直TEC的精度.最后,对全球不同地点的电离层TEC的仿真研究表明,在电子密度水平梯度较大的地区,应用电离层薄壳模型时会导致电子密度较高处的TEC被高估,而电子密度较低处的TEC被低估,在分析基于GPS数据提取的电离层TEC空间变化时要认识到这一点.     

Imaging of structures in the high-latitude ionosphere: model comparisons

The tomographic reconstruction technique generates a two-dimensional latitude versus height electron density distribution from sets of slant total electron content measurements (TEC) along ray paths between beacon satellites and ground-based radio receivers. In this note, the technique is applied to TEC values obtained from data simulated by the Sheffield/UCL/SEL Coupled Thermosphere/Ionosphere/Model (CTIM). A comparison of the resulting reconstructed image with the input modelled data allows for verification of the reconstruction technique. All the features of the high-latitude ionosphere in the model data are reproduced well in the tomographic image. Reconstructed vertical TEC values follow closely the modelled values, with the F-layer maximum density (NmF2) agreeing generally within about 10%. The method has also been able successfully to reproduce underlying auroral-E ionisation over a restricted latitudinal range in part of the image. The height of the F2 peak is generally in agreement to within about the vertical image resolution (25 km).     

Test of GPS for permanent ionospheric TEC monitoring at high latitudes

The Global Positioning System (GPS) observables are affected by the ionosphere. The dispersive nature of this effect and the use of two frequencies in the GPS observations make possible to measure the ionospheric total electron content (TEC) from dual frequency GPS data. In this work we test the concept of permanent monitoring of TEC using a network of GPS receivers at high latitudes. We have used GPS data from five permanent receivers in Scandinavia, from 1–30 January 1994, with geographic latitudes ranging from 57.4^°N to 78.9^°N. The results show the capability of the method to monitor the evolution of TEC as a function of time and geographical location. We have detected night-time enhancements almost every night for some of the stations, and we have also been able to produce maps of the instantaneous TEC as a function of both latitude and longitude around the GPS network. We also present some of the current limitations in the use of GPS for estimating TEC at high latitudes such as the difficulties in solving for cycle-slips, and the necessity of reliable values for the receiver and satellite differential instrumental biases.     

GPS-based ionospheric corrections for single frequency radar altimetry

The global ionospheric total electron content maps (GIMs) provide integrated electron densities between the ground and the GPS satellite altitude (20,200 km). Satellite altimeter ionospheric delay corrections require integrated electron densities between the ground and altimeter satellite altitude. In the case of the Geosat Follow-On (GFO) spacecraft, flying at 800 km, we estimated that using GIM TEC data alone, up to a 2 cm path delay can be introduced into the GFO measurements for high solar activity period by not taking into account the electron content above this altitude. Furthermore, the GIMs can have errors of 20–30 TECU in low latitudes for high solar activity in areas where there is little GPS data (such as over the oceans). In this paper, we describe the results of ingesting GIM TEC data into the International Reference Ionosphere model (IRI-95) to mitigate these two effects.     

Large-scale disturbances of auroral origin during strong magnetic storms of October 29–31, 2003, and November 7–11, 2004, according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered according to measurements of the total electron content (TEC) during the magnetic storms of October 29–31, 2003, and November 7–11, 2004, has been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region (foF2) were used for this purpose. The variations in TEC and foF2 were similar for all events mentioned above. The previous assumption that the region of thickness 150–200 km in the vicinity of the ionospheric F region mainly contributes to TEC modulation was confirmed for the cases when the electron density disturbance at an F region maximum was not more than 50%. However, this region probably becomes more extensive in vertical when the electron density disturbance in the vicinity of the ionospheric F region is about 85%.     

顶部电离层和等离子体层电子密度分布——基于GRACE星载GPS信标测量的CT反演

本文利用两颗跟飞的GRACE卫星载GPS信标测量数据和基于差分相对TEC的层析算法,实现了全球范围的顶部电离层和等离子体层(450~5000 km) 层析成像.反演结果表明,利用低轨道卫星载GPS信标测量数据可以有效地重建顶部电离层和等离子体层的全球二维分布图像.对不同地磁活动条件下的天基层析反演结果表明,等离子体层电子密度随纬度的分布是不均匀的;在低纬赤道带,从顶部电离层向上延伸直到等离子体层,以及等离子体层中局地的电离增强云团,经常出现近似垂直于磁力线的电子密度柱状增强结构.     

Tomographic imaging of the ionosphere using the GPS/MET and NNSS data

The earlier experiments of ionospheric tomography were conducted by receiving satellite signals from ground-based stations and then reconstructing electron density distribution from measures of the total electron content (TEC). In June 1994, National Central University built up the low-latitude ionospheric tomography network (LITN) including six ground stations spanning a range of 16.7° (from 14.6°N to 31.3°N) in latitude within 1° of 121°E longitude to receive the naval navigation satellite system (NNSS) signals (150 and 400 MHz). In the study of tomographic imaging of the ionosphere, TEC data from a network of ground-based stations can provide detailed information on the horizontal structure, but are of restricted utility in sensing vertical structure. However, an occultation observation mission termed the global positioning system/meteorology (GPS/MET) program used a low Earth orbiting (LEO) satellite (the MicroLab-1) to receive multi-channel GPS carrier phase signals (1.5 and 1.2 GHz) and demonstrate active limb sounding of the Earth's atmosphere and ionosphere. In this paper, we have implemented the multiplicative algebraic reconstruction technique (MART) to reconstruct and compare two-dimensional ionospheric structures from measured TECs through the receptions of the GPS signals, the NNSS signals, and/or both of the systems. We have also concluded the profiles retrieved from tomographic reconstruction showing much reasonable electron density results than the original vertical profiles retrieved by the Abel transformation and being in more agreement in peak electron density to nearby ionosonde measurements.     

基于非相干散射雷达和GPS观测研究Millstone Hill地区等离子体层电子含量

本文尝试结合非相干散射雷达和GPS TEC观测数据提取等离子体层总电子含量(PTEC).我们首先描述所用的技术方法,然后具体利用了Millstone Hill台站的观测数据研究该地区上空等离子体层总电子含量(PTEC)的变化情况.我们采用变化标高的Chapman函数对非相干散射雷达测得的电子浓度剖面数据进行拟合,然后通过对剖面积分得到100 km到1000 km高度范围的电离层总电子含量.GPS提供的TEC数据为高度达20200 km的总电子含量,两者之差可近似看成等离子体层的电子含量.本文分别选取太阳活动高年(2000, 2002年)和太阳活动低年(2005,2008年)Millstone Hill台站的静日数据进行研究.结果表明,等离子体层电子含量及其所占GPS TEC的比例具有明显的周日变化.PTEC含量在白天高于夜间,而所占GPS TEC的百分比,夜间明显高于白天.太阳活动高年所选月份等离子体层电子含量在4~14 TECU (1TECU=10¹⁶el/m²) 范围内变化,夜间所占比例可达60%左右.太阳活动低年所选月份等离子体层电子含量在3~7 TECU范围内变化,所占比例夜间最高可达80%左右.我们所得到的结果与前人基于其它观测手段所得结果在变化趋势上一致,在量级上也大致相当.因此,这从一个侧面证明了我们所用方法的可靠性.非相干散射雷达能够探测包括F2层峰值以下及以上高度的电子浓度,利用这一设备所观测得到的资料来推算电离层电子含量将比前人基于电离层垂测仪观测资料进行的推算更具真实性,由此得到的等离子体层电子含量也将更为接近真实情况.     

Ionospheric imaging at mid-latitudes using both GPS and ionosondes

Measurements from the Global Positioning System (GPS) satellites provide a valuable source of information about the ionosphere in the form of ray-path integrations of electron density. Total electron content (TEC) through the ionosphere can be estimated for specific satellite-to-ground paths using the two GPS frequencies and knowledge of the dispersive properties of the ionosphere. One approach used is the ionospheric imaging tool Multi Instrument Data Analysis System (MIDAS), which uses differential phase data from a number of GPS satellites and receivers to create an ionospheric movie of electron density. This paper addresses the accuracy with which MIDAS images the electron density at the F-layer peak. Firstly, the image accuracy is tested using a simulation of the imaging technique, representative of 1 year of data. Experimental GPS phase data are then used to image the electron density during a period of disturbed geomagnetic activity during April 2002. The images are compared to independent measurements from three ionosondes located across Europe and confirm the underestimate in peak electron density that was found in the simulation. Regardless of the peak density errors the vertical TEC in the images remains accurate. The accuracy of the imaged peak electron density is shown to improve across the image when measurements from ionosondes are included in the inversion process.     

GPS L-band scintillations and ionospheric irregularity zonal drifts inferred at equatorial and low-latitude regions

The main scientific objective of this research is to study the spatial variability and dynamics of the F-region irregularities. To achieve this, amplitude scintillations at the L-band, total electron content (TEC) and irregularity drifts were measured, as part of the Conjugate Point Equatorial Experiment (COPEX) campaign, by a network of ground-based global positioning system (GPS) receivers. The observations reveal a strong variability in the evolution of the irregularities from the equator to low-latitudes, and their zonal velocities at conjugate sites present a decrease with local time, and also with latitude. Moreover, the scintillations appear to be correlated with strong TEC gradients in the equatorward edge of the enhanced equatorial anomaly peaks. Other relevant aspects of the observations are highlighted and discussed.     

New approaches in global ionospheric determination using ground GPS data

Since 1 June 1998, the group of Astronomy and Geomatics of the Polytechnic University of Catalonia (gAGE/UPC) is contributing to the international project of defining an ionospheric product (Total Electron Content, TEC) from the data gathered by the permanent ground GPS receivers of the International GPS Service (IGS) network. The strategy and algorithms related to such a preliminary product, its calibration with synthetic observations generated from the International Reference Ionosphere (IRI), and the comparison with TOPEX TEC data are presented. Finally, these methods are applied combining ionosonde with ground GPS data, in order to obtain the vertical structure of the free electron distribution.     

GPS地面台网和掩星观测结合的时变三维电离层层析

本文给出GPS地面台网和掩星观测结合的时变三维电离层层析的原理、算法和基于实测数据的反演结果.反演结果的比较表明,联合地基GPS与掩星观测数据进行重建,电子密度整体图像的重建质量特别是其垂直结构的重建质量得到了明显改善.在平静日和磁暴期间两种条件下利用实测数据的重建结果表明,GPS地面台网和掩星观测结合的电离层层析可以获得电离层电子密度在高度-纬度-经度-时间四维空间中的变化.重建结果清晰地显示了磁暴期间电离层负相暴效应,表明结合GPS地面台网和掩星观测的时变三维电离层层析可以有效地监测扰动条件下的大尺度电离层结构.     

Equatorial ionization anomaly of the total electron content and equatorial electrojet of ground-based geomagnetic field strength

Measurements from ground-based receiver chains of the global positioning system (GPS) and magnetometers of the Circum-pan Pacific Magnetometer Network (CPMN) in the west Pacific region during 1999–2003 are examined. The ionospheric total electron content (TEC) derived from the GPS receivers is used to observe the strength, location, and occurrence time of the equatorial ionization anomaly (EIA) crests, which resulted from the equatorial plasma E×B drift fountain. The magnetic field strength of CPMN is employed to monitor the equatorial electrojet (EEJ), and to further estimate the effectiveness of the E×B drift to the EIA crests. Results show that the strength and location of the EIA crests are proportional to the EEJ strength.     

The November 2004 superstorm: Comparison of low-latitude TEC observations with LLIONS model results

We investigate the effects of penetration electric fields, meridional thermospheric neutral winds, and composition perturbation zones (CPZs) on the distribution of low-latitude plasma during the 7–11 November 2004 geomagnetic superstorm. The impact on low-latitude plasma was assessed using total electron content (TEC) measurements from a latitudinally distributed array of ground-based GPS receivers in South America. Jicamarca Radio Observatory incoherent scatter radar measurements of vertical E×B drift are used in combination with the Low-Latitude IONospheric Sector (LLIONS) model to examine how penetration electric fields and meridional neutral winds shape low-latitude TEC. It is found that superfountain conditions pertain between ～1900 and 2100 UT on 9 November, creating enhanced equatorial ionization anomaly (EIA) crests at ±20° geomagnetic latitude. Large-amplitude and/or long-duration changes in the electric field were found to produce significant changes in EIA plasma density and latitudinal location, with a delay time of ～2–2.5 h. Superfountain drifts were primarily responsible for EIA TEC levels; meridional winds were needed only to create hemispherical crest TEC asymmetries. The [O/N₂] density ratio (derived from the GUVI instrument, flown on the TIMED satellite) and measurements of total atmospheric density (from the GRACE satellites), combined with TEC measurements, yield information regarding a likely CPZ that appeared on 10 November, suppressing TEC for over 16 h.