Spatial interpolation of precipitation data is of great importance for hydrological modelling. Geostatistical methods (kriging) are widely applied in spatial interpolation from point measurement to continuous surfaces. The first step in kriging computation is the semi-variogram modelling which usually used only one variogram model for all-moment data. The objective of this paper was to develop different algorithms of spatial interpolation for daily rainfall on 1 km2 regular grids in the catchment area and to compare the results of geostatistical and deterministic approaches. This study leaned on 30-yr daily rainfall data of 70 raingages in the hilly landscape of the Ourthe and Ambleve catchments in Belgium (2908 km2). This area lies between 35 and 693 m in elevation and consists of river networks, which are tributaries of the Meuse River. For geostatistical algorithms, seven semi-variogram models (logarithmic, power, exponential, Gaussian, rational quadratic, spherical and penta-spherical) were fitted to daily sample semi-variogram on a daily basis. These seven variogram models were also adopted to avoid negative interpolated rainfall. The elevation, extracted from a digital elevation model, was incorporated into multivariate geostatistics. Seven validation raingages and cross validation were used to compare the interpolation performance of these algorithms applied to different densities of raingages. We found that between the seven variogram models used, the Gaussian model was the most frequently best fit. Using seven variogram models can avoid negative daily rainfall in ordinary kriging. The negative estimates of kriging were observed for convective more than stratiform rain. The performance of the different methods varied slightly according to the density of raingages, particularly between 8 and 70 raingages but it was much different for interpolation using 4 raingages. Spatial interpolation with the geostatistical and Inverse Distance Weighting (IDW) algorithms outperformed considerably the interpolation with the Thiessen polygon, commonly used in various hydrological models. Integrating elevation into Kriging with an External Drift (KED) and Ordinary Cokriging (OCK) did not improve the interpolation accuracy for daily rainfall. Ordinary Kriging (ORK) and IDW were considered to be the best methods, as they provided smallest RMSE value for nearly all cases. Care should be taken in applying UNK and KED when interpolating daily rainfall with very few neighbourhood sample points. These recommendations complement the results reported in the literature. ORK, UNK and KED using only spherical model offered a slightly better result whereas OCK using seven variogram models achieved better result.

An integrated hydrological-hydraulic model employing the 2-D local inertial equation as the core is established for effective numerical simulation of surface water flows in a great lake and its floodplain. The model is a cascade of validated hydrological and hydraulic sub-models. The model was applied to simulating the surface water flows of the Tonle Sap Lake and its floodplain in Cambodia using the roughness coefficient value calibrated comparing with a remote-sensing data set. The resulting model reasonably handles backwater flows during the rainy season and simulates the propagations of wet and dry interfaces without numerical instability, owing to a proper setting of time step supported by a novel numerical stability analysis. Sensitivity analysis of the surface water dynamics focusing on the setting of roughness coefficient and the backwater effect was also carried out. Overall, utilizing the 2-D local inertial equation in the assessment of lake water dynamics is a new modelling approach, which turns out to be an efficient simulation tool.

The applicability of multivariate interpolation and information entropy to optimize the raingauge network in the Mekong River Basin (MRB) is investigated. Three different spatial interpolation methods are tested: inverse distance squared (IDS), ordinary kriging (OK) and gradient plus inverse distance squared (GIDS). The validated results confirm that the GIDS method outperformed IDS and OK. The application of information entropy together with GIDS on a network of 57 gauges provided the same information content (7.34 nat) as could be obtained using all 6788 gauges in the MRB. Combining this result with meteorological and hydrological indicators revealed that the number of gauges for the optimum raingauge network could be reduced to 40. The results imply good applicability of the proposed method, which may be used to help prioritize efforts and funds to maintain the raingauge network in a given river basin.