Abstract
The signal of the GNSS satellites can be used to estimate the amount of water vapor in the atmosphere. For this reason, GNSS observations are nowadays routinely used by several meteorological institutes (e.g., MetOffice, Meteo France) to monitor weather events and to improve their weather forecasts quality. The analysis of a whole network of GNSS stations to estimate a full three-dimensional model of the water vapor content is a challenging and computationally demanding task. For this purpose, a tomographic system SEGAL GNSS Water Vapour Reconstruction Image Software (SWART) was developed and tested. The new method makes use of parallelized algebraic reconstruction techniques (ARTs) and supersedes other implementations in terms of speed by at least 50% for small networks. For SWART, the computation time grows linearly with the number of observations. As a result, the new method makes possible to estimate the water vapor for larger GNSS networks and can be used for near-real-time weather predictions. To show its potential, data from 26 stations in Poland were analyzed using data from a period of 56 days. Good agreement in the estimated water vapor between SWART and radiosondes solutions was obtained, with a mean RMS of 1.5 g/m3 for the lower layers and an overall improvement of 5% until the layer 6750 m when compared with the atmospheric model (WRF). Furthermore, rapid and strong variations observed by radiosondes were not modeled by the WRF but were detected by GPS tomography.
Similar content being viewed by others
Data availability
The data are available upon reasonable request to the corresponding author.
Notes
Noise is seen on the image and also known as impulse noise. It can be caused by sharp and sudden disturbance in the image signal. It presents itself as a sparsely occurring white-and-black pixel.
Hamming window is one of many window functions for smoothing values; it is recommended for smoothing and truncate the auto-covariance function in the time domain.
References
Ahmed F, Václavovic P, Teferle FN, Douša J, Bingley R, Laurichesse D (2016) Comparative analysis of real-time precise point positioning zenith total delay estimates. GPS Solutions 20:187. https://doi.org/10.1007/s10291-014-0427-z
Andersen A, Kak A (1982) Digital ray tracing in two-dimensional refractive fields. J Acoust Sot Am 72:1593–1606
Andersen A, Kak A (1984) Simultaneous algebraic reconstruction technique (SART): a superior implementation of the art algorithm. Ultrason Imaging 6(1):81–94
Andersson E (2018) Statement of guidance for global Numerical Weather Prediction (NWP). World Meteorological Organisation
Bender M, Dick G, Ge M, Deng Z, Wickert J, Kahle HG, Raabe A, Tetzlaff G (2011a) Development of a GNSS water vapour tomography system using algebraic reconstruction techniques. Adv Space Res 47(10):1704–1720
Bender M, Stosius R, Zus F, Dick G, Wickert J, Raabe A (2011b) GNSS water vapour tomography—expected improvements by combining GPS, GLONASS and Galileo observations. Adv Space Res 47(5):886–897
Bento F, Sá A, Crocker P, Fernandes RMS (2013) Analysis of algebraic reconstruction algorithms performance in CPU and GPU. ICEUBI 2013, Engineering for Economic Development, Covilhã, Portugal
Bevis M, Businger S, Herring TA, Rocken C, Anthes RA, Ware RH (1992) GPS meteorology: remote sensing of atmospheric water vapor using the Global Positioning System. J Geophys Res Atmos 97(D14):15787–15801
Bi Y, Mao J, Li C (2006) Preliminary results of 4-D water vapor tomography in the troposphere using GPS. Adv Atmospheric Sci 23:551–560
Boehm J, Schuh H (2013) Atmospheric effects in space geodesy. Springer, Heidelberg
Boehm J, Werl B, Schuh H (2006) Troposphere mapping functions for gps and very long baseline interferometry from european centre for medium-range weather forecasts operational analysis data. J Geophys Res Solid Earth 111(B2): https://doi.org/10.1029/2005JB00362978
Brenot H, Walpersdorf A, Reverdy M, Van Baelen J, Ducrocq V, Champollion C, Masson F, Doerflinger E, Collard P, Giroux P (2014) A GPS network for tropospheric tomography in the framework of the Mediterranean hydrometeorological observatory Cévennes-Vivarais (southeastern France). Atmos Meas Techn 7(2):553–578. https://doi.org/10.5194/amt-7-553-2014
Brenot H, Rohm W, Kacmarick M, Moller G, Sá A, Tondas D, Rapant L, Biondi R, Manning T, Champollion C (2020) Cross-comparison and methodological improvement in GPS tomography. Remote Sensing. https://doi.org/10.3390/rs12010030
Champollion C, Masson F, Bouin MN, Walpersdorf A, Doerflinger E, Bock O, Van Baelen J (2005) GPS water vapour tomography: preliminary results from the ESCOMPTE field experiment. Atmos Res 74(1):253–274
Ding N, Zhang S, Wu S, Wang X, Kealy A, Zhang K (2018) A new approach for GNSS tomography from few GNSS stations. Atmos Meas Tech 11:3511–3522. https://doi.org/10.5194/amt-11-3511-2018
Douša J, Dick G, Kačmařík M, Bro Řková R, Zus F, Brenot H, Stoycheva A, Möller G, Kaplon J (2016) Benchmark campaign and case study episode in central Europe for development and assessment of advanced GNSS tropospheric models and products. Atmos Meas Tech 9:2989–3008
Douša J, Václavovic P, Zhao L, Kačmařík M (2018) New adaptable all-in-one strategy for estimating advanced tropospheric parameters and using real-time orbits and clocks. Remote Sensing 10:232. https://doi.org/10.3390/rs10020232
Ducrocq V, Ricard D, Lafore JP, Orain F (2002) Storm-scale numerical rainfall prediction for five precipitations events over France: on the importance of the initial humidity field. Weather Forecast 17:1236–1256
Elgered G, Davis JL, Herring TA, Shapiro II (1991) Geodesy by radio interferometry: water vapour radiometry for estimation of the wet delay. J Geophys Res 96:6541–6555. https://doi.org/10.1029/90JB00834
Flores A, Ruffini G, Rius A (2000) 4D tropospheric tomography using GPS slant wet delays. In: Annales geophysicae, vol 18(2). Springer, Berlin, pp 223–234
Gradinarsky LP, Jarlemark P (2004) Ground-based GPS tomography of water vapor: analysis of simulated and real data. J Meteorol Soc Japan 82(1B):551–560
Guerova G, Jones J et al (2016) Review of the state of the art and future perspectives of the ground based GNSS meteorology in Europe. Atmos Meas Techn Discuss. https://doi.org/10.5194/amt-2016-125
Hadas T, Hobiger T, Hordyniec P (2020) Considering different recent advancements in GNSS on real-time zenith troposphere estimates. GPS Solut 24:99. https://doi.org/10.1007/s10291-020-01014-w
Haltiner G, Martin F (1957) Dynamical and physical meteorology. McGraw-Hill Book Company, New York
Hansen P, Saxild-Hansen M (2012) AIR tools—a MATLAB package of algebraic iterative reconstruction methods. J Comput Appl Math 236(8):2167–2178. https://doi.org/10.1016/j.cam.2011.09.039
Hordyniec P (2018) Simulation of liquid water and ice contributions to bending angle profiles in the radio occultation technique. Adv Space Res 62(5):1075–1089
Hordyniec P, Bosy J, Rohm W (2015) Assessment of errors in precipitable water data derived from global navigation satellite system observations. J Atmos Solar Terr Phys 129:69–77
Immler FJ, Dykema J, Gardiner T, Whiteman DN, Thorne P, Vomel H (2010) Reference quality upper-air measurements: guidance for developing GRUAN data products. Atmos Meas Tech 3(5):1217–1231
Innes PM, Dorling S (2012) Operational weather forecasting. Wiley, Hoboken
Jones J, Guerova G, Douša J, Dick G, de Haan S, Pottiaux E, Bock O, Pacione R, van Malderen R (eds) (2020) Advanced GNSS tropospheric products for monitoring severe weather events and climate. COST action ES1206 final action dissemination report. Springer, Berlin. https://doi.org/10.1007/978-3-030-13901-8
Kak A, Slaney M (1988) Principles of computerized tomographic imaging. IEEE Press, New York, pp 5–8
Kačmařík M, Dousa J, Dick G, Zus F, Brenot H, Möller G, Pottiaux E, Kapaon J, Hordyniec P, Václavovic P, Morel L (2017) Inter-technique validation of tropospheric slant total delays. Atmos Meas Tech 10(6):2183–2208
Kryza M, Werner M, Wałszek K, Dore AJ (2013) Application and evaluation of the WRF model for high-resolution forecasting of rainfall—a case study of SW Poland. Meteorol Z 22(5):595–601
Li X, Zus F, Lu C, Dick G, Ning T, Ge M, Wickert J, Schuh H (2015) Retrieving of atmospheric parameters from multi-GNSS in real time: validation with water vapor radiometer and numerical weather model. Geophys Res Lett Atmos. https://doi.org/10.1002/2015JD023454
Manning T, Zhang K, Rohm W, Choy S, Hurter F (2012) Detecting severe weather in Australia using GPS tomography. J Global Position Syst 11(1):58–70
McCarthy DD (ed) (1996) IERS Conventions
Miranda PMA, Alves JMR, Serra N (2013) Climate change and upwelling: response of Iberian upwelling to atmospheric forcing in a regional climate scenario. Clim Dyn 40(11–12):2813–2824
Möller G, Landskron D (2019) Atmospheric bending effects in GNSS tomography. Atmos Meas Tech 12:23–34. https://doi.org/10.5194/amt-12-23-2019
Nirvikar (2012) Algebraic reconstruction technique for computerized tomography with parallel Beam projection data. Doctoral dissertation, Shobhit University, Meerut
Notarpietro R, Cucca M, Gabella M, Venuti G, Perona G (2011) Tomographic reconstruction of wet and total refractivity fields from GNSS receiver networks. Adv Space Res 47(5):898–912
Peña A, Hasager CB, Gryning SE, Courtney M, Antoniou I, Mikkelsen T (2009) Offshore wind profiling using light detection and ranging measurements. Wind Energy Int Jo Prog Appl Wind Power Convers Technol 12(2):105–124
Reinders J (2007) Intel threading building blocks - outfitting C++ for multi-core processor parallelism. O’Reilly, 2007. [Online]. Available: http://www.oreilly.com/catalog/9780596514808/index.html
Reverdy M (2008) Estimation des paramètres atmosphériques par GPS: analyse de la variabilité spatio-temporelle de la vapeur d'eau. Doctoral dissertation, Université Blaise Pascal-Clermont-Ferrand II
Rius A, Ruffini G, Cucurull L (1997) Improving the vertical resolution of ionospheric tomography with GPS occultations. Geophys Res Lett 24(18):2291–2294
Rohm W (2013) The ground GNSS tomography-unconstrained approach. Adv Space Res 51(3):501–513
Rohm W, Bosy J (2009) Local tomography troposphere model over mountains area. Atmos Res 93(4):777–783
Rohm W, Zhang K, Bosy J (2014) Limited constraint, robust Kalman filtering for GNSS troposphere tomography. Atmos Meas Tech 7:1475–1486. https://doi.org/10.5194/amt-7-1475-2014
Saastamoinen J (1972) Atmospheric correction for the troposphere and stratosphere in radio ranging satellites. In: The use of artificial satellites for geodesy, pp 247–251
Schenewerk M (2003) A brief review of basic gps orbit interpolation strategies. GPS Solutions 6(4):265–267. https://doi.org/10.1007/s10291-002-0036-0
Shepp L, Logan B (1974) The Fourier reconstruction of a head section. IEEE Trans Nuclear Sci 21(3):21–43
Skamarock W, Klemp J, Dudhia J, Gill D, Barker D, Duda M, Huang X, Wang W, Powers J (2018) A description of the advanced research WRF version 3: NCAR technical note TN-475+str. National Center for Atmospheric Research
Sonntag D (1999) Important new values of the physical constants of 1986, vapour pressure formulations based on the ITS-90, and psychrometer formulae. Z Meterol 70(340):1990
Sá A (2019) PhD thesis “Tomographic determination of the spatial distribution of Water Vapor using GNSS observations for real-time applications”. https://www.researchgate.net/profile/Andre_Sa5/publications
Tregoning P, Boers R, O'Brien D, Hendy M (1998) Accuracy of absolute precipitable water vapor estimates from GPS observations. J Geophys Res Atmos 103(D22):28701–28710
Troller M, Geiger A, Brockmann E, Bettems JM, Bürki B, Kahle HG (2006a) Tomographic determination of the spatial distribution of water vapor using GPS observations. Adv Space Res 37(12):2211–2217
Troller M, Geiger A, Brockmann E, Kahle G (2006b) Determination of the spatial and temporal variation of tropospheric water vapour using CGPS networks. Geophys J Intern 167(2):509–520
Van Baelen J, Reverdy M, Tridon F, Labbouz L, Dick G, Bender M, Hagen M (2011) On the relationship between water vapour field evolution and the life cycle of precipitation systems. Q J R MeteorolSoc 137(S1):204–223
Wang P (2008) Compare windows* threads, openmp*, intel® threading building blocks for parallel programming. http://software.intel.com/en-us/blogs/2008/12/16/compare-windows-threads-openmp-intel-threading-building-blocks-for-parallel-programming/23. Accessed 15 Dec 2012
Wang Z, Li Z, Wang L, Wang X, Yuan H (2018) Assessment of multiple GNSS real-time SSR products from different analysis centers ISPRS. Int J Geo-Inf 7:85
Wilgan K, Rohm W, Bosy J (2015) Multi-observation meteorological and GNSS data comparison with numerical weather prediction model. Atmos Res 156:29–42
Wu L (2003) A parameter choice method for Tikhonov regularization ENTA Kent State University. Electron Trans Numer Anal 16:107–128
Yeh T, Hong J, Wang C, Hsiao T, Fong C (2014) Applying the water vapor radiometer to verify the precipitable water vapor measured by GPS. Terr Atmos Ocean Sci 25(2):189–201. https://doi.org/10.3319/TAO.2013.10.14.01(A)
Author information
Authors and Affiliations
Contributions
W.R., R.F. and A.S. designed the research; A.S. and E.T. prepared the data for the study; A.S., F.B. and M.B. performed the development of the software and parallelization implementation; A.S., E.T. and W.R. performed the GNSS tomography and tests for refinement; A.S., W.R., R.F. and M.B. analyzed and discussed the output results. A.S. and W.R. wrote the manuscript; all the authors read, revised and approved the final manuscript.
Corresponding author
Appendices
Appendix 1
1.1 Algebraic reconstruction techniques
Algebraic reconstruction technique is an iterative approach for imaging reconstruction using data obtained from a series of projections such as those obtained from electron microscopy, X-ray photography and in medical imaging like in computed axial tomography—CAT scans (Kak and Slaney 1988). The projection data are usually obtained by a small-angle rotation around the object. It is not possible to measure a large number of projections in every situation, and it may also happen that they are not uniformly distributed on complete angular distance, i.e., on 180° or 360° (used in 3D), which is the case of GNSS tomography (Nirvikar 2012).
There are two different techniques for algebraic reconstruction, namely algebraic reconstruction techniques (ARTs) and simultaneous iterative reconstruction techniques (SIRTs).
1.2 ART methods
The algebraic reconstruction techniques (ARTs) are row action methods that treat the equations one at the time. This means that in each iteration, each equation is solved individually. The updates in each ith iteration are made using Eq. 14 plus a relaxation parameter \({\uplambda }_{{\text{k}}}\) described by:
What distinguishes the various methods is the order in each row is processed (Hansen and Saxild-Hansen 2012). The ART reconstructions usually suffer from salt-and-pepper noise,Footnote 1 which is caused by the inconsistencies introduced in the set of equations by the \({\text{w}}_{{{\text{ij}}}}\) approximations as they are usually measured by experiments (Kak and Slaney 1988).
1.3 SIRT methods
The simultaneous iterative reconstruction techniques (SIRTs) are named “simultaneous” as all equations are solved at the same time in one iteration, based on matrix multiplications (Hansen and Saxild-Hansen 2012). Usually, these methods converge slower to the solution than the ART methods; however, they result in better looking images (Kak and Slaney 1988). The SIRT methods work by first solving all the equations, without changing the cell values, and then only at the end of each iteration, the cell values are changed/updated, being the average value of all computed changes for that cell.
The general form of these methods is as follows:
where A is the matrix with the various projections, \({\uplambda }_{{\text{k}}}\) is a relaxation parameter and the matrices M and T are diagonal and symmetric positive definite. The various SIRT methods depend on these matrices. \({\uplambda }_{{\text{k}}}\) is defined as \(\frac{2}{{{\uprho }\left( {{\text{TA}}^{{\text{T}}} {\text{MA}}} \right)}} - {\upvarepsilon }\), where \({{ \rho }}\) is the spectral radius and \({\upvarepsilon }\) the machine´s epsilon.
1.4 Simultaneous algebraic reconstruction technique (SART)
SART is a kind of variation on the algebraic approaches already discussed (ART and SIRT) that combine the best of both. This technique was first reported in Andersen and Kak (1982). The main features can be described as: First, to reduce errors in the approximation of ray integrals of a smooth image by finite sums, the traditional pixel basis is abandoned in favor of bilinear elements. To further reduce the noise resulting from the unavoidable, but now presumably considerably smaller inconsistencies with real data projection, the correction terms are simultaneously applied for all the rays in one projection; this contrasts with the ray-by-ray updates in ART. In addition, a heuristic procedure is used to improve the quality of reconstructions: A longitudinal Hamming windowFootnote 2 is used to emphasize the corrections applied near the middle of a ray relative to those applied near its ends.
Appendix 2
Statistics for all levels and both 1-h and 24-h solutions with respective Student’s p values. The p values were calculated as a two-sided test results for mean (bias) of residuals between SWART and RS as well as WRF and RS.
Appendix 3
Number of examples selected from all 56 days of Case 2, containing inversions in vertical profiles.
See Figs.
8,
9,
10,
11,
12,
13,
14 and
15.
Rights and permissions
About this article
Cite this article
Sá, A., Rohm, W., Fernandes, R.M. et al. Approach to leveraging real-time GNSS tomography usage. J Geod 95, 8 (2021). https://doi.org/10.1007/s00190-020-01464-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00190-020-01464-7