 This study aimed to develop algorithms for spatial interpolation of daily rainfall data on 1 km and LTE, SUP and GT, 2 and LTE slash SUP and GT, regular grids in the Arda and Amblev catchments in Belgium. The authors compared the results of geostatistical and deterministic approaches using seven semivariogram models, including Gaussian, logarithmic, power, exponential, rational quadratic, spherical, and pentaspherical. They found that the Gaussian model was the most frequently best fit for daily rainfall data. The negative estimates of cryging were observed for convective more than stratiform rain. Spatial interpolation with geostatistical and inverse distance weighting, IDW, algorithms outperformed the Thiessen polygon method, while integrating elevation into cryging with an external drift, KED, and ordinary co-croaging, OCK, did not improve the interpolation accuracy for daily rainfall. The authors recommended using ordinary cryging, ORK, UNK, and KED when interpolating daily rainfall with very few neighboring sample points. This article was authored by S. Ly, C. Charles, and A. Degra. We are article.tv, links in the description below.