利用精密单点定位求解电离层延迟

近年来,高时空分辨率的全球导航卫星系统(GNSS)观测信号已成为电离层研究的重要资源.利用GNSS研究电离层,需首先将观测资料转换成包含电离层信息的可观测量(Ionospheric 0bservables,称之为"电离层观测值").目前,最常用的电离层观测值一般采用联合无几何影响组合的码和相位观测,利用相位平滑伪距方法...     

A neural network approach for regional vertical total electron content modelling

A Neural Network model has been developed for estimating the total electron content (TEC) of the ionosphere. TEC is proportional to the delay suffered by electromagnetic signals crossing the ionosphere and is among the errors that impact GNSS (Global Navigation Satellite Systems) observations. Ionospheric delay is particularly a problem for single frequency receivers, which cannot eliminate the (first-order) ionospheric delay by combining observations at two frequencies. Single frequency users rely on applying corrections based on prediction models or on regional models formed based on actual data collected by a network of receivers. A regional model based on a neural network has been designed and tested using data sets collected by the Brazilian GPS Network (RMBC) covering periods of low and high solar activity. Analysis of the results indicates that the model is capable of recovering, on average, 85% of TEC values.     

Ionospheric space weather effects monitored by simultaneous ground and space based GNSS signals

Ionospheric space weather effects can degrade the performance of global navigation satellite systems (GNSS), i.e. their accuracy, reliability and availability. However, well established ground based and innovative space based GNSS measurements offer the unique chance for a permanent monitoring of the electron density structure of the global ionosphere–plasmasphere system. In this paper we review various types of perturbations in the ionospheric plasma density and distribution. In order to analyze these space weather effects we use 30 s sampled measurements provided by the global GPS ground tracking network of the IGS. Furthermore, to get a more comprehensive view on the perturbations analyzed also are simultaneously obtained GPS measurements onboard the LEO satellite CHAMP (challenging minisatellite payload). Whereas the ground based measurements show strong horizontal redistribution of plasma during ionospheric storms, the space-borne measurements indicate a severe vertical redistribution of the ionospheric plasma during the selected events. The role of the various dynamical forces such as meridional winds and electric fields is also discussed.     

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.     

地震电离层前兆综合监测技术

地震前后对电离层进行观测的方法正在逐渐发展成为一种新的地震前兆观测方法。应用多手段地基无线电遥感探测技术,通过软件无线电形式,实现一机多能,对电离层状态的综合探测,获取地震前电离层变态的多参数数据,结合对地震电离层前兆物理机制的研究,将能大大提高地震电离层前兆预报分析的可信度、可靠性,有助于对地震电磁现象机制的识别和前兆信息的提取和强震预报预测。     

Using GPS–GLONASS–GALILEO data and IRI modeling for ionospheric calibration of radio telescopes and radio interferometers

VHF and HF radio signals are widely used to observe the Sun and pulsars. Nowadays, large low-frequency radio astronomical arrays (LOFAR, 30–240 MHz; MIRA, 80–300 MHz) are being constructed to record radiation of pulsars at the maximum distance. registration of the solar radio emission intensity at fixed frequencies and in the spectral VHF band is very important along with other methods of monitoring of coronal mass ejections. Interpreting radio astronomical data is known to be necessary to take into account possible distortions of these signals in the Earth ionosphere. However, in contrast to modern navigation systems (Global Position System (GPS), GLObal NAvigation Satellite System (GLONASS), GALILEO), in which a very accurate reconstruction of ionosphere parameters is a built-in function, in present-day radio astronomy a retrieve of ionosphere transfer characteristics has not been appropriately worked out. This collides with increasing requirements to accuracy of the analysis of a radio emission amplitude profile and to the angular and polarizing resolution of radio telescopes of new generation (LOFAR, SKA, etc.). We have developed a method and software to calculate the ionosphere rotation measure (RM) and dispersion measure (DM). We used the ionosphere model IRI-2001, magnetic field model IGRF-10, and the ionosphere total electron content values obtained from GPS measurements. The obtained values of DM and RM were recalculated into characteristics of the phase delay, Faraday amplitude modulation, and polarization changes. We calculated ones for different levels of geomagnetic activity as well as different angular positions of radio sources. Our main idea is to use a signal of navigation satellites (GPS, GLONASS, GALILEO) as a testing signal from a “reference” source located at minimal angle distance from a source studied. Our project allows development of methods and systems of ADAPTIVE RADIO ASTRONOMY, adaptive to the non-uniform and non-stationary ionosphere, by analogy with known systems of adaptive optics intended to adapt optical telescopes to varying conditions of the optically non-uniform and non-stationary troposphere.     

地基GNSS VTEC约束的大气掩星电离层误差修正方法

介绍了一种考虑电离层沿GNSS掩星射线路径分布不对称且兼顾一阶和二阶项的大气掩星电离层误差修正新方法,它综合地基GNSS VTEC的水平变化信息和电离层模式的垂直变化信息,在沿"入射线"与"出射线"双边局部球对称假设下,估算GNSS掩星射线路径上的电子密度;进而计算一阶和二阶项弯曲角电离层误差廓线.采用太阳活动低年的2008年7月15日和太阳活动较高年的2013年7月15日两天的MetOp-A和GRACE掩星观测资料和IGS的GNSS VTEC数据,计算了GNSS大气掩星弯曲角一阶和二阶项电离层误差廓线.对比分析表明：新方法经验模型和理论模型的二阶项弯曲角电离层残差,以及经验模型与实测数据的一阶和二阶项弯曲角电离层误差均具有良好的一致性,因此该方法可用于L2信号数据质量较差或L2信号中断的掩星事件,同时修正大气掩星一阶和二阶电离层误差,从而提高弯曲角精度和掩星观测资料的利用率.     

Improving capabilities of broadband differential satellite navigation systems via radio occultation technology

The existent satellite system for radio occultation monitoring the Earth’s neutral atmosphere and ionosphere (COSMIC) provides data to consumers in the regions with limited possibilities of constructing dense measurement networks (e.g., in the World Ocean area). A forthcoming increase of LEO small spacecrafts and the deployment of new satellite radio navigation systems will result in a pronounced increase in the efficiency of radio occultation method and its space resolution. As a result, the Space-Based Augmentation Systems (SBAS) broadband differential system will become global, or the quality of corrections delivered to single-frequency consumers of individual systems, e.g., the Augmentation and Monitoring System, will be improved. Therefore, the methods for processing and analyzing obtained radio occultation data should be improved. A simple method to reconstruct the electron density profile at radio occultation points, based on the total electron content measurement on the satellite–satellite path and the IRI-type ionospheric model has been proposed. The method needs initial information, it does not require refraction measurements, and it is free of the assumption that the ionosphere is spherically stratified in the occultation region. Verification of the proposed method based on data for 121 radio occultation cases across Europe in May 2013 demonstrated good agreement with the vertical sounding data.     

Effects on the Ionosphere Due to Phenomena Occurring Below it

The terrestrial thermosphere and ionosphere form the most variable part of theEarth's atmosphere. Because our society depends on technological systems thatcan be affected by thermospheric and ionospheric phenomena, understanding,monitoring and ultimately forecasting the changes of the thermosphere–ionosphere system are of crucial importance to communications, navigation and the exploration of near-Earth space. The reason for the extreme variability of the thermosphere–ionosphere system isits rapid response to external forcing from various sources, i.e., thesolar ionizing flux, energetic charged particles and electric fields imposed via the interaction between the solar wind, magnetosphere and ionosphere, as well as coupling from below (meteorological influences) by the upward propagating, broad spectrum,internal atmospheric waves (planetary waves, tides, gravity waves) generated in thestratosphere and troposphere. Thunderstorms, typhoons, hurricanes, tornadoes andeven seismological events may also have observable consequences in the ionosphere.The release of trace gases due to human activity have the potential to cause changes inthe lower and the upper atmosphere.A brief overview is presented concerning the discoveries and experimentalresults that have confirmed that the ionosphere is subject to meteorologicalcontrol (especially for geomagnetic quiet conditions and for middle latitudes).D-region aeronomy, the winter anomaly of radiowave absorption, wave-liketravelling ionospheric disturbances, the non-zonality and regional peculiaritiesof lower thermospheric winds, sporadic-E occurrence and structure, spread-Fevents, the variability of ionospheric electron density profiles and Total ElectronContent, the variability of foF2, etc., should all be considered in connection withtropospheric and stratospheric processes. Ionospheric weather, as a part of spaceweather, (i.e., hour-to-hour and day-to-day variability of the ionospheric parameters)awaits explanation and prediction within the framework of the climatological, seasonal,and solar-cycle variations.     

Anomalous modification of the ionospheric total electron content prior to the 26 September 2005 Peru earthquake

This paper investigates the features of pre-earthquake ionospheric anomalies in the total electron content (TEC) data obtained on the basis of regular GPS observations from the International GNSS Service (IGS) network. For the analysis of the ionospheric effects of the 26 September 2005 Peru earthquake, Global Ionospheric Maps (GIMs) of TEC were used. The possible influence of the earthquake preparation processes on the main low-latitude ionosphere peculiarity—the equatorial anomaly—is discussed. Analysis of the TEC maps has shown that modification of the equatorial anomaly occurred a few days before the earthquake. In previous days, during the evening and night hours (local time—LT), a specific transformation of the TEC distribution had taken place. This modification took the shape of a double-crest structure with a trough near the epicenter, though usually in this time the restored normal latitudinal distribution with a maximum near the magnetic equator is observed. Additional measurements (CHAMP satellite) have also confirmed the presence of this structure. To compare the vertical TEC measurements obtained with GPS satellite signals (GPS TEC), the International Reference Ionosphere, IRI-2001, was used for calculating the IRI TEC.     

Large-scale ionospheric gradients over Europe observed in October 2003

It is well known that ionospheric perturbations are characterised by strong horizontal gradients and rapid changes of the ionisation. Thus, space weather induced severe ionosphere perturbations can cause serious technological problems in Global Navigation Satellite Systems (GNSS) such as GPS. During the severe ionosphere storm period of 29–31 October 2003, reported were several significant malfunctions due to the adverse effects of the ionosphere perturbations such as interruption of the WAAS service and degradation of mid-latitudes GPS reference services. To properly warn service users of such effects, a quick evaluation of the current signal propagation conditions expressed in a suitable ionospheric perturbation index would be of great benefit. Preliminary results of a comparative study of ionospheric gradients including vertical sounding and Total Electron Content (TEC) data are presented. Strong enhancements of latitudinal gradients and temporal changes of the ionisation are observed over Europe during the 29–30 October storm period. The potential use of spatial gradients and rate of change of foF2 and TEC characterising the actual perturbation degree of the ionosphere is discussed. It has been found that perturbation induced spatial gradients of TEC and foF2 strongly enhance during the ionospheric storm on 29 October over the Central European region in particular in North–South direction exceeding the gradients in East–West direction by a factor of 2.     

大气掩星反演误差特性初步分析

GPS大气掩星探测技术可以获得全球大气折射率、气压、密度、温度和水汽压等气象参数,该技术基本原理是基于几何光学近似的Abel积分反演.地球扁率、电离层传播时间延迟、大气大尺度水平梯度、多路径传播现象等因素在某些高度范围影响大气反演的精度.本文采用模拟的方法,分析其中地球扁率及电离层对反演结果的影响,并讨论局部圆弧修正及电离层修正的效果.利用CHAMP掩星实测轨道数据和有关电离层和大气经验模式、采用三维射线追踪方法模拟计算几种情形下的GPS掩星观测附加相位数据,对模拟数据进行反演,将反演气象参量剖面与模拟时给定模式剖面进行比较,得到了0～60 km高度范围内的反演误差.误差统计分析结果表明,局部圆弧中心的修正以及电离层修正,对于高精度的GPS掩星反演是非常重要的;电离层修正残差仍是制约30～60 km高度范围内反演精度的重要因素,进一步完善和优化大气掩星反演需要发展新的电离层修正算法.     

O<Superscript>+</Superscript>-Reduced Models of the High Latitude Ionosphere and the Ionospheric Alfvén Resonator in Broadband Pc1 Events

The significance of the O⁺-ion density altitude profile of the outer ionosphere for determination of the Ionospheric Alfvén Resonator (IAR) lower harmonic structure has been demonstrated. The O⁺-reduced and exponentially extrapolated ionosphere models at high altitudes are generally acceptable for the IAR interpretation of subauroral broadband Pc1 events. Instantaneous ionospheric plasma data based on simultaneous EISCAT (CP-1, CP-7) measurements should be most suitable for the interpretation of different pulsation events. The limited applicability of the averaged International Reference Ionosphere (IRI) models has also been demonstrated.     

Ionosphere Disturbances Preceding Earthquakes according to the Data of Ground Based Station of the Vertical Ionospheric Sounding Wakkanai

The paper analyzes the data of manual ionograms processing of hourly measurements of the critical frequency foF2 of the F2 ionospheric layer at the Wakkanai ionospheric vertical sounding station (Japan) in a geomagnetically quiet environment before a series of earthquakes with magnitude M > 6.0, for which the station entered the earthquake preparation zone, in order to detect possible Ionospheric Disturbances Preceding Earthquakes (IDPE), and to determine their quantitative characteristics. Detected IDPE, in the opinion of the authors, can be related to the processes of preparation of the corresponding earthquakes, i.e., to be Ionospheric Precursors of Earthquakes (IPE).     

Kalman filter-based algorithms for monitoring the ionosphere and plasmasphere with GPS in near-real time

Data collected from a GPS receiver located at low latitudes in the American sector are used to investigate the performance of the WinTEC algorithm [Anghel et al., 2008a, Kalman filter-based algorithm for near realtime monitoring of the ionosphere using dual frequency GPS data. GPS Solutions, accepted for publication; for different ionospheric modeling techniques: the single-shell linear, quadratic, and cubic approaches, and the multi-shell linear approach. Our results indicate that the quadratic and cubic approaches perform much better than the single-shell and multi-shell linear approaches in terms of post-fit residuals. The performance of the algorithm for the cubic approach is then further tested by comparing the vertical TEC predicted by WinTEC and USTEC [Spencer et al., 2004. Ionospheric data assimilation methods for geodetic applications. In: Proceedings of IEEE PLANS, Monterey, CA, 26–29 April, pp. 510–517] at five North American stations. In addition, since the GPS-derived total electron content (TEC) contains contributions from both ionospheric and plasmaspheric sections of the GPS ray paths, in an effort to improve the accuracy of the TEC retrievals, a new data assimilation module that uses background information from an empirical plasmaspheric model [Gallagher et al., 1988. An empirical model of the Earth's plasmasphere. Advances in Space Research 8, (8)15–(8)24] has been incorporated into the WinTEC algorithm. The new Kalman filter-based algorithm estimates both the ionospheric and plasmaspheric electron contents, the combined satellite and receiver biases, and the estimation error covariance matrix, in a single-site or network solution. To evaluate the effect of the plasmaspheric component on the estimated biases and total TEC and to assess the performance of the newly developed algorithm, we compare the WinTEC results, with and without the plasmaspheric term included, at three GPS receivers located at different latitudes in the American sector, during a solar minimum period characterized by quiet and moderate geomagnetic conditions. We also investigate the consistency of our plasmaspheric results by taking advantage of the specific donut-shaped geometry of the plasmasphere and applying the technique at 12 stations distributed roughly over four geomagnetic latitudes and three longitude sectors.     

Anamolous ionospheric TEC variations prior to the Indonesian earthquake (M 7.1) of November 15, 2014

This paper investigates preearthquake ionospheric variations with the use of TEC of Global Ionospheric Maps (GIMs) and regional maps based on Precise Point Positioning (PPP) during the 7.1-M Indonesian earthquake that occurred on November 15, 2014. TEC maps corresponding to 10 days before to 4 days after the event were examined. In addition, a time series of TEC values according to the PPP maps were also evaluated. In addition to GIMs, it was possible to detect TEC variations with PPP maps. The results showed that ionospheric TEC decreased strikingly 4 days prior to the earthquake. This TEC variation was highly likely related to seismic activity.     

Features of High-Latitude Ionospheric Irregularities Development as Revealed by Ground-Based GPS Observations,Satellite-Borne GPS Observations and Satellite In Situ Measurements over the Territory of Russia during the Geomagnetic Storm on March 17–18, 2015

The dynamic picture of the response of the high- and mid-latitude ionosphere to the strong geomagnetic disturbances on March 17–18, 2015, has been studied with ground-based and satellite observations, mainly, by transionospheric measurements of delays of GPS (Global Positioning System) signals. The advantages of the joint use of ground-based GPS measurements and GPS measurements on board of the Swarm Low-Earth-Orbit satellite mission for monitoring of the appearance of ionospheric irregularities over the territory of Russia are shown for the first time. The results of analysis of ground-based and space-borne GPS observations, as well as satellite, in situ measurements, revealed large-scale ionospheric plasma irregularities observed over the territory of Russia in the latitude range of 50°–85° N during the main phase of the geomagnetic storm. The most intense ionospheric irregularities were detected in the auroral zone and in the region of the main ionospheric trough (MIT). It has been found that sharp changes in the phase of the carrier frequency of the navigation signal from all tracked satellites were recorded at all GPS stations located to the North from 55° MLAT. The development of a deep MIT was related to dynamic processes in the subauroral ionosphere, in particular, with electric fields of the intense subauroral polarization stream. Analysis of the electron and ion density values obtained by instruments on board of the Swarm and DMSP satellites showed that the zone of highly structured auroral ionosphere extended at least to heights of 850–900 km.     

基于北斗、GLONASS和GPS系统的中低纬电离层特性联合探测

基于GNSS(Global Navigation Satellite Systems)的发展,我们利用具有北斗、GLONASS和GPS三系统信号接收功能的接收机观测的数据,结合电离层总电子含量(Total Electron Content, TEC)的反演算法,提取出GNSS三系统观测的电离层TEC;同时,将GNSS三系统获取的TEC应用到电离层TEC地图、行进式扰动、不规则体结构和电离层的太阳耀斑响应等方面的研究中,这也是首次使用三种GNSS系统数据对电离层进行联合探测研究.研究结果表明,增加了北斗系统的GNSS三系统在研究中国地区电离层TEC地图、周日变化、逐日变化,行进式扰动以及电离层的实时监测等方面较单系统的GPS具有明显的优势.     

Effect of horizontal gradients on ionospherically reflected or transionospheric paths using a precise homing-in method

A homing-in method is presented for determining ionospheric reflected or transionospheric paths between fixed transmitter and receiver locations in the presence of ionospheric gradients or ripples. Both initial elevation and azimuth are automatically adjusted to find the path that arrives exactly at the receiver. The method can be used for any 3D ionospheric model to find precise ray paths and phase and group delays for both magneto-ionic modes. The method takes full account of path location, geomagnetic field orientation and the bending of the ray path resulting from horizontal as well as vertical gradients of electron density. It can also find multiple paths e.g. low and high angle, 1- and 2-hops for both ordinary and extraordinary modes. Examples of its use are given for both terrestrial HF links and Earth to Satellite paths. For paths reflected from the ionosphere, the effect of gradients of both critical frequency and height of maximum electron density are determined and the comparative effect of gradients on high and low angle and 1- and 2-hops paths for both magneto-ionic modes investigated. Path variation with frequency for a fixed link is also studied and the bandwidth of the ionospheric background channel (dispersive bandwidth) and its reciprocal (the pulse rise time), important for wideband digital HF broadcasting or spread spectrum HF communications, is estimated for a range of frequencies, for high- and low-angle rays and 1- and 2-hop paths. For Earth–satellite paths, the effect of the ionosphere and horizontal ionospheric gradients is determined for a range of frequencies and elevation angles. It is shown that the method can also enable the determination of second-order errors in satellite navigation methods, such as GPS, due to ionospheric gradients and the effect of the geomagnetic field.     

Prospects for Scrutiny of Pulsars with Polish Part of LOFAR

LOFAR (LOw Frequency ARray) is a new generation digitally controlled radio telescope consisting of phased array antenna stations with sensitivity, bandwidth, range of frequency, and digital processing power that makes it an excellent tool for observations of pulsars. This interferometric instrument is able to work in a single-station mode as well as in group-of-selected-stations mode. This article discusses the great opportunity for conducting unique and independent research of pulsar sources with the three LOFAR stations located in Poland and maintained by the POLFAR consortium.