Dear Syed: This is a good question. In our experience there are two possible answers: 1) There is not an unique interpolation method for all the possible life situation when runoff mapping is required, so we need to try several techniques and every time, with a Jackknife assessment, select the method with the best performance for the studied situation. 2) The interpolation method should follow the nature of the studied process, in such case, being runoff a stochastic process then optimal interpolation methods should be used. In this line you need to follow methods like Kolmogorov optimal interpolation, Gandin interpolation and Kriging. Being the last the most well known due to broadly distributed tools where this method is implemented, nevertheless, the method that generalize all stochastic ones is the Kolmogorov optimal interpolation. That is from the theoretical point of view, from the practice, from our experience, there is no method able to provide best performance than 40% Mean Absolute Relative Error due to the low density of hydrological monitoring networks that are not enough to support a better result. Also, the majority of discharge observations are heavily influenced by human presure introducing additional variability that can not be handled by the interpolation methods. My sugestion: Try to find an interpolation code for the Kolmogorov optimal interpolation, if there is not any then apply Kriging.
For more research see the following links. Best regards,
Sir, As you know main sources of surface runoff are Precipitation in terms of Snowfall and rainfall. Satellite data interpretation gives only idea of Terrain, Landuse/Landcover and soil type. so if you want to calculate surface runoff you need precipitation information as well as water discharge in the river outlet.
Dear people, I would recommend the use of a simplified hydrological model capable of deriving runoff from precipitation and evapotranspiration or whatever hydrological process you would like to consider as well as the time step and spatial extent. I guess the interpolation of precipitation and potential evapotranspiration is easier being two disperse variables.
such a simple conceptual model is for instance GR2M, which proved very robust in many different kinds of climates.
I put an example of a paper on this topic, one about GR2M, and a second one about the Yates model, which needs no calibration and gave consistent results, compared to the Yates model.
In order to have better results, you will need some discharge measurements. Otherwise, everything else will be in the sphere of stochasticity or approach with great uncertainty...