Based on the geo-statistical theory and ArcGIS geo-statistical module, datas of 30 groundwater level observation wells were used to estimate the decline of groundwater level in Beijing piedmont. Seven different interpolation methods (inverse distance weighted interpolation, global polynomial interpolation, local polynomial interpolation, tension spline interpolation, ordinary Kriging interpolation, simple Kriging interpolation and universal Kriging interpolation) were used for interpolating groundwater level between 2001 and 2013. Cross-validation, absolute error and coefficient of determination (R2) was applied to evaluate the accuracy of different methods. The result shows that simple Kriging method gave the best fit. The analysis of spatial and temporal variability suggest that the nugget effects from 2001 to 2013 were increasing, which means the spatial correlation weakened gradually under the influence of human activities. The spatial variability in the middle areas of the alluvial–proluvial fan is relatively higher than area in top and bottom. Since the changes of the land use, groundwater level also has a temporal variation, the average decline rate of groundwater level between 2007 and 2013 increases compared with 2001–2006. Urban development and population growth cause over-exploitation of residential and industrial areas. The decline rate of the groundwater level in residential, industrial and river areas is relatively high, while the decreasing of farmland area and development of water-saving irrigation reduce the quantity of water using by agriculture and decline rate of groundwater level in agricultural area is not significant.

In the karst area of southern China, karst water is important for supporting the sustainable production and home living for the local residents. Consequently, it is of significance to fully understand the water cycle, so as to make full use of water resources. In karst areas, epikarst and conduits are developed, participating in the hydrological cycle actively. For conventional lumped hydrologic models, it is difficult to simulate the hydrological cycle accurately. These models neglect to consider the variation of underlying surface and weather change. Meanwhile, for the original distributed hydrological model, the existence of epikarst and underground conduits as well as inadequate data information also make it difficult to achieve accurate simulation. To this end, the framework combining the advantages of lumped model–reservoir model and distributed hydrologic model–Soil and Water Assessment Tool (SWAT) model is established to simulate the water cycle efficiently in a karst area. Xianghualing karst watershed in southern China was selected as the study area and the improved SWAT model was used to simulate the water cycle. Results show that the indicators of E_NS and R² in the calibration and verification periods are both above 0.8, which is evidently improved in comparison with the original model. The improved SWAT model is verified to have better efficiency in describing the hydrological cycle in a typical karst area.