As an important GIS function, spatial interpolation is one of the most often used geographic techniques for spatial query, spatial data visualization, and spatial decision-making processes in GIS and environmental science. However, less attention has been paid on the comparisons of available spatial interpolation methods, although a number of GIS models including inverse distance weighting, spline, radial basis functions, and the typical geostatistical models (i.e. ordinary kriging, universal kriging, and cokriging) are already incorporated in GIS software packages. In this research, the conceptual and methodological aspects of regression kriging and GIS built-in interpolation models and their interpolation performance are compared and evaluated. Regression kriging is the combination of multivariate regression and kriging. It takes into consideration the spatial autocorrelation of the variable of interest, the correlation between the variable of interest and auxiliary variables (e.g., remotely sensed images are often relatively easy to obtain as auxiliary variables), and the unbiased spatial estimation with minimized variance. To assess the efficiency of regression kriging and the difference between stochastic and deterministic interpolation methods, three case studies with strong, medium, and weak correlation between the response and auxiliary variables are compared to assess interpolation performances. Results indicate that regression kriging has the potential to significantly improve spatial prediction accuracy even when using a weakly correlated auxiliary variable.

The dynamics of surface water play a crucial role in the hydrological cycle and are sensitive to climate change and anthropogenic activities, especially for the agricultural zone. As one of the most populous areas in China’s river basins, the surface water in the Huai River Basin has significant impacts on agricultural plants, ecological balance, and socioeconomic development. However, it is unclear how water areas responded to climate change and anthropogenic water exploitation in the past decades. To understand the changes in water surface areas in the Huai River Basin, this study used the available 16,760 scenes Landsat TM, ETM+, and OLI images in this region from 1989 to 2017 and processed the data on the Google Earth Engine (GEE) platform. The vegetation index and water index were used to quantify the spatiotemporal variability of the surface water area changes over the years. The major results include: (1) The maximum area, the average area, and the seasonal variation of surface water in the Huai River Basin showed a downward trend in the past 29 years, and the year-long surface water areas showed a slight upward trend; (2) the surface water area was positively correlated with precipitation (p < 0.05), but was negatively correlated with the temperature and evapotranspiration; (3) the changes of the total area of water bodies were mainly determined by the 216 larger water bodies (>10 km²). Understanding the variations in water body areas and the controlling factors could support the designation and implementation of sustainable water management practices in agricultural, industrial, and domestic usages.