Evaluation of IWV from the numerical weather prediction WRFmodel with PPP GNSS processing for Bulgaria

Abstract. Global Navigation Satellite Systems (GNSS) meteorology is an established operational service providing hourly updated GNSS tropospheric products to the National Meteorologic Services (NMS) in Europe. In the last decade through the ground-based GNSS network densification and new processing strategies like Precise Point Positioning (PPP) it has become possible to obtain sub-hourly tropospheric products for monitoring severe weather events. In this work one year (January–December 2013) of sub-hourly GNSS tropospheric products (Zenith Total Delay) are computed using the PPP strategy for seven stations in Bulgaria. In order to take advantage of the sub-hourly GNSS data to derive Integrated Water Vapour (IWV) surface pressure and temperature with similar temporal resolution is required. As the surface observations are on 3 hourly basis the first step is to compare the surface pressure and temperature from numerical weather prediction model Weather Forecasting and Research (WRF) with observations at three synoptic stations in Bulgaria. The mean difference between the two data-sets for 1) surface pressure is less than 0.5 hPa and the correlation is over 0.989, 2) temperature the largest mean difference is 1.1 °C and the correlation coefficient is over 0.957 and 3) IWV mean difference is in range 0.1–1.1 mm. The evaluation of WRF on annual bases shows IWV underestimation between 0.5 and 1.5 mm at five stations and overestimation at Varna and Rozhen. Varna and Rozhen have also much smaller correlation 0.9 and 0.76. The study of the monthly IWV variation shows that at those locations the GNSS IWV has unexpected drop in April and March respectively. The reason for this drop is likely problems with station raw data. At the remaining 5 stations a very good agreement between GNSS and WRF is observed with high correlation during the cold part of 2013 i.e. March, October and December (0.95) and low correlation during the warm part of 2013 i.e. April to August (below 0.9). The diurnal cycle of the WRF model shows a dry bias in the range of 0.5-1.5 mm. Between 00 and 01 UTC the GNSS IWV tends to be underestimate IWV which is likely due to the processing window used. The precipitation efficiency from GNSS and WRF show very good agreement on monthly bases with a maximum in May-June and minimum in August–September. The annual precipitation efficiency in 2013 at Lovech and Burgas is about 6 %.

[1]  I. Shapiro,et al.  Geodesy by radio interferometry: Effects of atmospheric modeling errors on estimates of baseline length , 1985 .

[2]  Jan Douša,et al.  Real-time zenith tropospheric delays in support of numerical weather prediction applications , 2014 .

[3]  Witold Rohm,et al.  Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data , 2017, Journal of Geodesy.

[4]  R. Van Malderen,et al.  A multi-site intercomparison of integrated water vapour observations for climate change analysis , 2014 .

[5]  Eric Pottiaux,et al.  Review of the state of the art and future prospects of the ground-based GNSS meteorology in Europe , 2016 .

[6]  Guergana Guerova,et al.  Validation of NWP Mesoscale Models with Swiss GPS Network AGNES , 2003 .

[7]  Roscoe R. Braham,et al.  THE WATER AND ENERGY BUDGETS OF THE THUNDERSTORM AND THEIR RELATION TO THUNDERSTORM DEVELOPMENT , 1952 .

[8]  G. D. Thayer,et al.  An improved equation for the radio refractive index of air , 1974 .

[9]  L. S. Pereira,et al.  Ground-Based GPS Measurements of Precipitable Water Vapor and Their Usefulness for Hydrological Applications , 2014, Water Resources Management.

[10]  Chris Rizos,et al.  The International GNSS Service in a changing landscape of Global Navigation Satellite Systems , 2009 .

[11]  Christian Melsheimer,et al.  A multi-instrument comparison of integrated water vapour measurements at a high latitude site , 2012 .

[12]  K. Rannat,et al.  Column water vapour: an intertechnique comparison of estimation methods in Estonia , 2014 .

[13]  Song‐You Hong,et al.  The WRF Single-Moment 6-Class Microphysics Scheme (WSM6) , 2006 .

[14]  Urs Hugentobler,et al.  The International GNSS Real-Time Service , 2012 .

[15]  M. Tomassini,et al.  Monitoring IWV from GPS and limited - area forecast model , 2003 .

[16]  D. Laurichesse,et al.  Comparative analysis of real-time precise point positioning zenith total delay estimates , 2016, GPS Solutions.

[17]  J. Dudhia,et al.  A New Vertical Diffusion Package with an Explicit Treatment of Entrainment Processes , 2006 .

[18]  E. Mlawer,et al.  Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave , 1997 .

[19]  H. Schuh,et al.  Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data , 2006 .

[20]  Jordan G. Powers,et al.  A Description of the Advanced Research WRF Version 2 , 2005 .

[21]  Witold Rohm,et al.  Real‐time retrieval of precipitable water vapor from GPS precise point positioning , 2014 .

[22]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[23]  Y. Xue,et al.  Modeling of land surface evaporation by four schemes and comparison with FIFE observations , 1996 .

[24]  Galina Dick,et al.  Monitoring of Integrated Water Vapour from ground-based GPS observations and their assimilation in a limited-area NWP model , 2002 .

[25]  Qile Zhao,et al.  Assessment of precipitable water vapor derived from ground-based BeiDou observations with Precise Point Positioning approach , 2015 .

[26]  Gunnar Elgered,et al.  Multi-technique comparisons of 10 years of wet delay estimates on the west coast of Sweden , 2012, Journal of Geodesy.