Evaluating Low-Cost GNSS Network Densification for Water-Vapor Tomography over an Urban Area: A Case Study over Lisbon

This study evaluates GNSS water-vapor tomography across the Lisbon metropolitan area and explores how increasing network density with low-cost receivers improves three-dimensional humidity fields for meteorological applications. Three configurations were tested for December 2022, a month characterized by several rainfall events, including a severe urban-impacting one: (i) a hybrid setup combining permanent and low-cost stations (TOMO_PL), (ii) a dense network of only low-cost stations (TOMO_L), (iii) a sparse arrangement using only permanent stations (TOMO_P). Tomographic water vapor density fields were compared with independent references from the Weather Research and Forecasting (WRF) model, ERA 5 reanalysis, and radiosonde data. All products show the expected exponential decline in water vapor with increasing altitude. Tomography consistently underestimates moisture in the lowest 2.0 to 2.5 km and tends to overestimate it at higher levels, with a weaker correlation above mid-tropospheric heights. Vertical RMSE remains below 2 g m−3 for all solutions, but TOMO_P performs the worst due to weak and uneven spatial geometry. Time–height analysis reveals that densified setups capture the changing moisture in the lower atmosphere, including increased near-surface humidity during December 11–13, when rainfall exceeded 120 mm in 24 h, although mid-level intrusions and dry layers observed by radiosondes are not captured. Mean PWV patterns show realistically low points over the Sintra mountain range and align best with TOMO_PL (spatial RMSE 0.6 g m−3, bias 0.4 g m−3, correlation 0.9), while TOMO_P creates artifacts that mimic mesoscale gradients. Categorized skill analysis shows the highest accuracy under high-moisture conditions and limited ability to detect dry conditions, with TOMO_PL showing the best overall performance against both ERA5 and WRF. Overall, low-cost densification significantly enhances boundary-layer humidity and PWV retrievals, supporting their use for urban heavy-rain monitoring and, with error-aware integration, for short-term forecasting.