Development of a GNSS water vapour tomography system using algebraic reconstruction techniques

A GNSS water vapour tomography system developed to reconstruct spatially resolved humidity fields in the troposphere is described. The tomography system was designed to process the slant path delays of about 270 German GNSS stations in near real-time with a temporal resolution of 30 min, a horizontal resolution of 40 km and a vertical resolution of 500 m or better. After a short introduction to the GPS slant delay processing the framework of the GNSS tomography is described in detail. Different implementations of the iterative algebraic reconstruction techniques (ART) used to invert the linear inverse problem are discussed. It was found that the multiplicative techniques (MART) provide the best results with least processing time, i.e., a tomographic reconstruction of about 26,000 slant delays on a 8280 cell grid can be obtained in less than 10 min. Different iterative reconstruction techniques are compared with respect to their convergence behaviour and some numerical parameters. The inversion can be considerably stabilized by using additional non-GNSS observations and implementing various constraints. Different strategies for initialising the tomography and utilizing extra information are discussed. At last an example of a reconstructed field of the wet refractivity is presented and compared to the corresponding distribution of the integrated water vapour, an analysis of a numerical weather model (COSMO-DE) and some radiosonde profiles.     

GNSS water vapour tomography – Expected improvements by combining GPS,GLONASS and Galileo observations

The German Research Centre for Geosciences (GFZ) operates a GNSS water vapour tomography system using about 350 German GNSS stations. The GNSS data processing at the GFZ works in near real-time and provides zenith total delays, integrated water vapour and slant delay data operationally. This large data set of more than 50,000 slant delays per hour is used to reconstruct spatially resolved humidity fields by means of tomographic techniques. It can be expected that additional observations from the future Galileo system provide more information with improved quality. A simulation study covering 12 h at 14 July 2009 was therefore started to estimate the impact of GPS, Galileo and GLONASS data on the GNSS tomography. It is shown that the spatial coverage of the atmosphere with slant paths is highly improved by combining observations from two or three satellite systems. Equally important for a reliable tomographic reconstruction is the distribution of slant path intersections as they are required to locate the integrated delay information. The number of intersection points can be increased by a factor of 4 or 8 if two or three systems are combined and their distribution will cover larger regions of the atmosphere. The combined data sets can be used to increase the spatiotemporal resolution of the reconstructed humidity fields up to 30 km horizontally, 300 m vertically and 15 min. The reconstruction quality could not be improved considerably using the currently available techniques.     

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.     

Characterisation of residual ionospheric errors in bending angles using GNSS RO end-to-end simulations

Global Navigation Satellite System (GNSS) radio occultation (RO) is an innovative meteorological remote sensing technique for measuring atmospheric parameters such as refractivity, temperature, water vapour and pressure for the improvement of numerical weather prediction (NWP) and global climate monitoring (GCM). GNSS RO has many unique characteristics including global coverage, long-term stability of observations, as well as high accuracy and high vertical resolution of the derived atmospheric profiles. One of the main error sources in GNSS RO observations that significantly affect the accuracy of the derived atmospheric parameters in the stratosphere is the ionospheric error. In order to mitigate the effect of this error, the linear ionospheric correction approach for dual-frequency GNSS RO observations is commonly used. However, the residual ionospheric errors (RIEs) can be still significant, especially when large ionospheric disturbances occur and prevail such as during the periods of active space weather. In this study, the RIEs were investigated under different local time, propagation direction and solar activity conditions and their effects on RO bending angles are characterised using end-to-end simulations. A three-step simulation study was designed to investigate the characteristics of the RIEs through comparing the bending angles with and without the effects of the RIEs. This research forms an important step forward in improving the accuracy of the atmospheric profiles derived from the GNSS RO technique.     

Improving the estimation of zenith dry tropospheric delays using regional surface meteorological data

Global Navigation Satellite Systems (GNSS) are emerging as possible tools for remote sensing high-resolution atmospheric water vapour that improves weather forecasting through numerical weather prediction models. Nowadays, the GNSS-derived tropospheric zenith total delay (ZTD), comprising zenith dry delay (ZDD) and zenith wet delay (ZWD), is achievable with sub-centimetre accuracy. However, if no representative near-site meteorological information is available, the quality of the ZDD derived from tropospheric models is degraded, leading to inaccurate estimation of the water vapour component ZWD as difference between ZTD and ZDD. On the basis of freely accessible regional surface meteorological data, this paper proposes a height-dependent linear correction model for a priori ZDD. By applying the ordinary least-squares estimation (OLSE), bootstrapping (BOOT), and leave-one-out cross-validation (CROS) methods, the model parameters are estimated and analysed with respect to outlier detection. The model validation is carried out using GNSS stations with near-site meteorological measurements. The results verify the efficiency of the proposed ZDD correction model, showing a significant reduction in the mean bias from several centimetres to about 5 mm. The OLSE method enables a fast computation, while the CROS procedure allows for outlier detection. All the three methods produce consistent results after outlier elimination, which improves the regression quality by about 20% and the model accuracy by up to 30%.     

High temporal resolution global PWV dataset of 2005–2016 by using a neural network approach to determine the mean temperature of the atmosphere

Atmospheric water vapour plays an important role in phenomena related to the global hydrologic cycle and climate change. However, the rapid temporal–spatial variation in global tropospheric water vapour has not been well investigated due to a lack of long-term, high-temporal-resolution precipitable water vapour (PWV). Accordingly, this study generates an hourly PWV dataset for 272 ground-based International Global Navigation Satellite System (GNSS) Service (IGS) stations over the period of 2005–2016 using the zenith troposphere delay (ZTD) derived from global-scale GNSS observation. The root mean square (RMS) of the hourly ZTD obtained from the IGS tropospheric product is approximately 4 mm. A fifth-generation reanalysis dataset of the European Centre for Medium-range Weather Forecasting (ECMWF ERA5) is used to obtain hourly surface temperature (T) and pressure (P), which are first validated with GNSS synoptic station data and radiosonde data, respectively. Then, T and P are used to calculate the water vapour-weighted atmospheric mean temperature (Tm) and zenith hydrostatic delay (ZHD), respectively. T and P at the GNSS stations are obtained via an interpolation in the horizontal and vertical directions using the grid-based ERA5 reanalysis dataset. Here, Tm is calculated using a neural network model, whereas ZHD is obtained using an empirical Saastamoinen model. The RMS values of T and P at the collocated 693 radiosonde stations are 1.6 K and 3.1 hPa, respectively. Therefore, the theoretical error of PWV caused by the errors in ZTD, T and P is on the order of approximately 2.1 mm. A practical comparison experiment is performed using 97 collocated radiosonde stations and 23 GNSS stations equipped with meteorological sensors. The RMS and bias of the hourly PWV dataset are 2.87/?0.16 and 2.45/0.55 mm, respectively, when compared with radiosonde and GNSS stations equipped with meteorological sensors. Additionally, preliminary analysis of the hourly PWV dataset during the EI Niño event of 2014–2016 further indicates the capability of monitoring the daily changes in atmospheric water vapour. This finding is interesting and significant for further climate research.     

Tropospheric delays from GNSS for application in coastal altimetry

In the scope of the development of an improved methodology for the computation of the wet tropospheric correction for coastal altimetry, based on the use of tropospheric delays derived from GNSS (Global Navigation Satellite Systems), various studies have been conducted aiming to improve the estimation, at global scale, of GNSS-derived tropospheric delays.     

The verification of GNSS tropospheric tomography model in a mountainous area

The GNSS (Global Navigation Satellite System) has not been developed as a meteorological data source provider, but with a careful and sophisticated processing strategy it might be used as one. The term GNSS tomography refers to the usage of the ray traced GNSS signal as scanning rays in the tomographic model input. The model is divided into a number of voxels. The system is inverted and value of refractivity is obtained. Typically, as in the most of the inverse processing, there is a problem of the undetermined system and as a consequence the cofactor matrix is close to singular. To avoid singularity additional conditions or constrains should be added to the system. Here, additional parameters are derived with the help of the air flow analysis in the Sudety mountains (south-west region of Poland), and special Slant Wet Delay (SWD) trimming procedure. The flow’s synthetic parameters like the Bruint-Väisälä frequency and the Froude number are determined. This way the type of the flow is recognized and the analysis of the impact of orographic barrier has been quantified. The SWDs from the GNSS observations were tested against, SWD from raytracing through the COAMPS model field. The modified GNSS tomography model was tested for the real GNSS observations delivered from the GNSS network Karkonosze located in the Sudety mountains and compared with the COAMPS model. The solution shows a considerable improvement in comparison with plain tomographic model results.     

Near real-time estimation of tropospheric water vapour content from ground based GNSS data and its potential contribution to weather now-casting in Austria

The importance of high resolution meteorological analysis of the atmosphere increased over the past years. A detailed analysis of the humidity field is an important precondition for a better monitoring of local and regional extreme precipitation events and for forecasts with improved spatial resolution. For this reason, the Austrian Meteorological Agency (ZAMG) is operating the spatial and temporal high resolution INCA system (Integrated Now-casting through Comprehensive Analysis) since begin of 2005. Errors in this analysis occur mainly in the areas of rapidly changing and hard to predict weather conditions or rugged topography with extreme differences in height such as the alpine area of Austria. The aim of this work is to provide GNSS based measurements of the tropospheric water vapour content with a temporal resolution of 1 h and a temporal delay of less than 1 h to assimilate these estimates into the INCA system. Additional requirement is an accuracy of better than 1 mm of the precipitable water (PW) estimates.     

An estimation of tropospheric corrections using GPS and synoptic data: Improving Urmia Lake water level time series from Jason-2 and SARAL/AltiKa satellite altimetry

Tropospheric correction is one of the most important corrections in satellite altimetry measurements. Tropospheric wet and dry path delays have strong dependence on temperature, pressure and humidity. Tropospheric layer has particularly high variability over coastal regions due to humidity, wind and temperature gradients. Depending on the extent of water body and wind conditions over an inland water, Wet Tropospheric Correction (WTC) is within the ranges from a few centimeters to tens of centimeters. Therefore, an extra care is needed to estimate tropospheric corrections on the altimetric measurements over inland waters. This study assesses the role of tropospheric correction on the altimetric measurements over the Urmia Lake in Iran. For this purpose, four types of tropospheric corrections have been used: (i) microwave radiometer (MWR) observations, (ii) tropospheric corrections computed from meteorological models, (iii) GPS observations and (iv) synoptic station data. They have been applied to Jason-2 track no. 133 and SARAL/AltiKa track no. 741 and 356 corresponding to 117–153 and the 23–34 cycles, respectively. In addition, the corresponding measurements of PISTACH and PEACHI, include new retracking method and an innovative wet tropospheric correction, have also been used. Our results show that GPS observation leads to the most accurate tropospheric correction. The results obtained from the PISTACH and PEACHI projects confirm those obtained with the standard SGDR, i.e., the role of GPS in improving the tropospheric corrections. It is inferred that the MWR data from Jason-2 mission is appropriate for the tropospheric corrections, however the SARAL/AltiKa one is not proper because Jason-2 possesses an enhanced WTC near the coast. Furthermore, virtual stations are defined for assessment of the results in terms of time series of Water Level Height (WLH). The results show that GPS tropospheric corrections lead to the most accurate WLH estimation for the selected virtual stations, which improves the accuracy of the obtained WLH time series by about 5%.     

Analysis of precipitable water vapour characteristics from GNSS measurements during the snow season in Liaoning Province,China

Global Navigation Satellite System (GNSS) remote sensing precipitable water vapour (PWV) data from November 2015 to March 2019 were combined with snowfall observation data and used to analyse PWV characteristics in Liaoning Province during the snow season (from November to March the following year) and their relationship with snowfall. The potential of using GNSS for PWV measurements was demonstrated using sounding data with a correlation coefficient higher than 0.9 and a mean bias error lower than 0.5 mm. According to the GNSS PWV data gathered at 30-min intervals from 68 GNSS stations in Liaoning during the snow season, the monthly PWV average was highest in November and lowest in January. Negative correlations were found between PWV and altitude. Most of the water vapour was concentrated in the low layer of the atmosphere, and the contribution of this vapour to the PWV was higher during the snow season than in summer. A total of 43 snow cases were identified using the snowfall records from 53 GNSS stations, and the characteristics of PWV during these snowfalls were analysed. An increase in PWV was observed before snowfall events. Moreover, the influence of synoptic systems and air mass origins on PWV was analysed based on National Centers for Environmental Prediction (NCEP) reanalysis data and the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. The results show that the water vapour condition was better when the synoptic systems or air masses came from areas south of Liaoning.     

卫星精密测定轨中的对流层延迟估算系统

基于V isual Basic编程,设计了卫星精密测定轨中的实时气象采集与对流层延迟估算系统。利用实时采集ZQZ-A型自动气象站的温度、相对湿度、大气压等气象要素,采用Saastamoinen模型进行估算天顶方向对流层静力延迟和湿延迟,同时采用地基双波段微波辐射计测量大气水汽总量和云液态水总量,进而估算更高精度的湿延迟,在其有效时与模型计算的静力延迟相加,得到天顶方向总延迟;无效时采用模型计算的总延迟。最后将其投影到卫星方向得到总延迟。实际应用表明,设计的对流层延迟估算系统人机界面友好,具有易用性,能实时、稳定地为卫星精密测定轨提供服务。     

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.     

Real-time zenith tropospheric delays in support of numerical weather prediction applications

The Geodetic Observatory Pecný (GOP) routinely estimates near real-time zenith total delays (ZTD) from GPS permanent stations for assimilation in numerical weather prediction (NWP) models more than 12 years. Besides European regional, global and GPS and GLONASS solutions, we have recently developed real-time estimates aimed at supporting NWP nowcasting or severe weather event monitoring. While all previous solutions are based on data batch processing in a network mode, the real-time solution exploits real-time global orbits and clocks from the International GNSS Service (IGS) and Precise Point Positioning (PPP) processing strategy. New application G-Nut/Tefnut has been developed and real-time ZTDs have been continuously processed in the nine-month demonstration campaign (February–October, 2013) for selected 36 European and global stations. Resulting ZTDs can be characterized by mean standard deviations of 6–10 mm, but still remaining large biases up to 20 mm due to missing precise models in the software. These results fulfilled threshold requirements for the operational NWP nowcasting (i.e. 30 mm in ZTD). Since remaining ZTD biases can be effectively eliminated using the bias-reduction procedure prior to the assimilation, results are approaching the target requirements in terms of relative accuracy (i.e. 6 mm in ZTD). Real-time strategy and software are under the development and we foresee further improvements in reducing biases and in optimizing the accuracy within required timeliness. The real-time products from the International GNSS Service were found accurate and stable for supporting PPP-based tropospheric estimates for the NWP nowcasting.     

The impact of microwave absorber and radome geometries on GNSS measurements of station coordinates and atmospheric water vapour

We have used microwave absorbing material in different geometries around ground-based Global Navigation Satellite System (GNSS) antennas in order to mitigate multipath effects on the estimates of station coordinates and atmospheric water vapour. The influence of a hemispheric radome – of the same type as in the Swedish GPS network SWEPOS – was also investigated. Two GNSS stations at the Onsala Space Observatory were used forming a 12 m baseline. GPS data from October 2008 to November 2009 were analyzed by the GIPSY/OASIS II software using the Precise Point Positioning (PPP) processing strategy for five different elevation cutoff angles from 5° to 25°. We found that the use of the absorbing material decreases the offset in the estimated vertical component of the baseline from ∼27 mm to ∼4 mm when the elevation cutoff angle varies from 5° to 20°. The horizontal components are much less affected. The corresponding offset in the estimates of the atmospheric Integrated Water Vapour (IWV) decreases from ∼1.6 kg/m² to ∼0.3 kg/m². Changes less than 5 mm in the offsets in the vertical component of the baseline are seen for all five elevation cutoff angle solutions when the antenna was covered by a hemispheric radome. Using the radome affects the IWV estimates less than 0.4 kg/m² for all different solutions. IWV comparisons between a Water Vapour Radiometer (WVR) and the GPS data give consistent results.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

GNSS reflectometry and remote sensing: New objectives and results

The Global Navigation Satellite System (GNSS) has been a very powerful and important contributor to all scientific questions related to precise positioning on Earth’s surface, particularly as a mature technique in geodesy and geosciences. With the development of GNSS as a satellite microwave (L-band) technique, more and wider applications and new potentials are explored and utilized. The versatile and available GNSS signals can image the Earth’s surface environments as a new, highly precise, continuous, all-weather and near-real-time remote sensing tool. The refracted signals from GNSS radio occultation satellites together with ground GNSS observations can provide the high-resolution tropospheric water vapor, temperature and pressure, tropopause parameters and ionospheric total electron content (TEC) and electron density profile as well. The GNSS reflected signals from the ocean and land surface could determine the ocean height, wind speed and wind direction of ocean surface, soil moisture, ice and snow thickness. In this paper, GNSS remote sensing applications in the atmosphere, oceans, land and hydrology are presented as well as new objectives and results discussed.     

Effect of the troposphere on surface neutron counter measurements

Surface neutron counter data are often used as a proxy for atmospheric ionisation from cosmic rays in studies of extraterrestrial effects on climate. Neutron counter instrumentation was developed in the 1950s and relationships between neutron counts, ionisation and meteorological conditions were investigated thoroughly using the techniques available at the time; the analysis can now be extended using modern data. Whilst surface neutron counts are shown to be a good proxy for ionisation rate, the usual meteorological correction applied to surface neutron measurements, using surface atmospheric pressure, does not completely compensate for tropospheric effects on neutron data. Residual correlations remain between neutron counts, atmospheric pressure and geopotential height, obtained from meteorological reanalysis data. These correlations may be caused by variations in the height and temperature of the atmospheric layer at ∼100 hPa. This is where the primary cosmic rays interact with atmospheric air, producing a cascade of secondary ionising particles.     

Tomographic determination of the spatial distribution of water vapor using GPS observations

With the advent of the GPS navigation system, a promising ground based technique has been introduced which makes it possible to estimate the amount of water vapor in the troposphere from operational GPS networks at relatively low additional costs. While the estimation of the integrated amount is currently well established, the determination of the spatial water vapor distribution and its temporal variation are still a major challenge. To account for the vertical resolution, several tomographic approaches were pursued. We developed the software package AWATOS (atmospheric water vapor tomography software) which is based on the assimilation of double differenced GPS observations. Applying a least-squares inversion, the inhomogeneous spatial distribution of water vapor is determined. An extensive investigation has been carried out in Switzerland. GPS measurements are performed by the dense permanent Swiss national GPS network AGNES of the Swiss Federal Office of Topography (swisstopo). A total of 40 equally distributed water vapor profiles have been estimated on an hourly basis. For the purpose of validation, 22 radiosonde profiles were used at the GPS and meteorological station Payerne. Furthermore, data of the numerical weather model aLMo (alpine model in Switzerland, MeteoSwiss) were compared with the tomographic results. An overall good agreement of the three methods with an rms of better than 1.6 g/m³ absolute humidity was achieved. The results show that AGNES can be used as a dedicated network for the purpose of GPS-tomography, using a horizontal resolution of approximately 50 km and height layers of 300–500 m thickness in the lower troposphere.     

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.