Spatial interpolation of precipitation from multiple rain gauge networks and weather radar data for operational applications in Alpine catchments

Increasing meteorological data availability and quality implies an adaptation of the interpolation methods for data combination. In this paper, we propose a new method to efficiently combine weather radar data with data from two heated rain gauge networks of different quality. The two networks being non-collocated (no common location between the two networks), pseudo cross-variograms are used to compute the linear model of coregionalization for kriging computation. This allows considering the two networks independently in a co-kriging approach. The methodology is applied to the Upper Rhône River basin, an Alpine catchment in Switzerland with a complex topography and an area of about 5300 km2. The analysis explores the newly proposed Regression cokriging approach, in which two independent rain gauge networks are considered as primary and secondary kriging variables. Regression co-kriging is compared to four other methods, including the commonly applied Inverse distance weighting method used as baseline scenario. Incorporation of additional networks located within and around the target region in the interpolation computation is also explored. The results firstly demonstrate the added value of the radar information as compared to using only ground stations. As compared to Regression kriging using only the network of highest quality, the Regression co-kriging method using both networks slightly increases the performance. A key outcome of the study is that Regression co-kriging performs better than Inverse distance weighting even for the data availability scenario when the radar network was providing lower quality radar data over the studied basin. The results and discussion underline that combining meteorological information from different rain gauge networks with different equipments remains challenging for operational purposes. Future research in this field should in particular focus on additional pre-processing of the radar data to account for example for areas of low visibility of the weather radars due to the topography.