Geostatistical Merging of a Single-Polarized X-Band Weather Radar and a Sparse Rain Gauge Network over an Urban Catchment

International audienceOptimal Quantitative Precipitation Estimation (QPE) of rainfall is crucial to the accuracy of hydrological models, especially over urban catchments. Small-to-medium size towns are often equipped with sparse rain gauge networks that struggle to capture the variability in rainfall over high spatiotemporal resolutions. X-band Local Area Weather Radars (LAWRs) provide a cost-effective solution to meet this challenge. The Clermont Auvergne metropolis monitors precipitation through a network of 13 rain gauges with a temporal resolution of 5 min. 5 additional rain gauges with a 6-minute temporal resolution are available in the region, and are operated by the national weather service Météo-France. The LaMP (Laboratoire de Météorologie Physique) laboratory's X-band single-polarized weather radar monitors precipitation as well in the region. In this study, three geostatistical interpolation techniques-Ordinary kriging (OK), which was applied to rain gauge data with a variogram inferred from radar data, conditional merging (CM), and kriging with an external drift (KED)-are evaluated and compared through cross-validation. The performance of the inverse distance weighting interpolation technique (IDW), which was applied to rain gauge data only, was investigated as well, in order to evaluate the effect of incorporating radar data on the QPE's quality. The dataset is comprised of rainfall events that occurred during the seasons of summer 2013 and winter 2015, and is exploited at three temporal resolutions: 5, 30, and 60 min. The investigation of the interpolation techniques performances is carried out for both seasons and for the three temporal resolutions using raw radar data, radar data corrected from attenuation, and the mean field bias, successively. The superiority of the geostatistical techniques compared to the inverse distance weighting method was verified with an average relative improvement of 54% and 31% in terms of bias reduction for kriging with an external drift and conditional merging, respectively (cross-validation). KED and OK performed similarly well, while CM lagged behind in terms of point measurement QPE accuracy, but was the best method in terms of preserving the observations' variance. The correction schemes had mixed effects on the multivariate geostatistical methods. Indeed, while the attenuation correction improved KED across the board, the mean field bias correction effects were marginal. Both radar data correction schemes resulted in a decrease of the ability of CM to preserve the observations variance, while slightly improving its point measurement QPE accuracy