Blending multi-satellite, atmospheric reanalysis and gauge precipitation products to facilitate hydrological modelling

Abstract Satellite-retrieved and atmospheric reanalysis precipitation can bridge the spatiotemporal gaps of in-situ gauging networks, but estimation biases can limit their reliable applications in hydrological monitoring and modelling. To correct precipitation occurrence and intensity simultaneously, this study develops a three-stage blending approach to integrate three multi-satellite precipitation datasets (IMERG Final, TMPA 3B42V7 and PERSIANN-CDR), the ERA5 atmospheric reanalysis product and a gauge dataset within a dynamic framework. Firstly, the systematic biases of the four members were individually corrected by combining the local intensity scaling and ratio bias correction methods. Then, the “state weights” used for determining wet/dry events were optimized by evaluating a score function of the four bias-corrected members. Thirdly, the “intensity weights” were optimized using the cuckoo search (CS) algorithm and the Bayesian Model Averaging (BMA) method, respectively. The three-stage blending approach produced dynamic weights varying both spatially and temporally, and the performance was thoroughly evaluated over mainland China. Results show that the three-stage dynamic scheme performs better than individual datasets and two-stage blending methods in terms of all eight statistical metrics, and the CS algorithm outperforms the BMA method in the third stage. By randomly sampling validation sites using K-fold experiments, the developed algorithm also demonstrates a superior performance in ungauged regions. After interpolating and normalizing blending parameters of all gauges to entire domain using ordinary kriging, a new blended precipitation dataset with a daily 0.25° scale was produced. Four hydrological models are forced by blended and primary precipitations in 238 catchments over China, further confirming that the developed approach can facilitate hydrological modelling demonstrated by improving the Kling-Gupta efficiency of simulated streamflow by 12-35%.