Original Articles

Study on Spatial Interpolation of Precipitation with Large Scale Samples: A Case Study on 2009’s Precipitation of China

  • 1. Institute of Geographical Sciences and Natural Resources Research, CAS, Beijing 100101, China;
    2. Graduate University of Chinese Academy of Sciences , Beijing 100049, China;
    3. School of Geography, Beijing Normal University, Beijing 100875, China

Received date: 2010-08-01

  Revised date: 2010-12-01

  Online published: 2011-07-25


This paper studies spatial interpolation of precipitation at a national scale, the precipitation data comes from 2203 meteorological stations of China in 2009. The research content includes three parts as follows. At first, this paper divided the stations into different groups by stochastic methods, including 90%, 80%, 70%, 60%, 50%, 40%, 30%, 20%, 10%, 5% and 4% sampling rate groups, and analyzed the impact of sampling data on interpolation result by inverse distances weighting methods. Secondly, the stations with 80%, 50% and 20% of sampling rates were treated by the same sampling interval method, which extracted sampling data with gaps equal 1 or 4 according to ID serial number. By comparing interpolation result with the result by stochastic method in the same sampling rate, this paper analyzed the relationship between sampling methods of data and the result of interpolation. Finally, the paper analyzed the differences of interpolation result by IDW, Kriging and Co-Kriging in the same sampling rate groups. We can draw some conclusions as follows: (1) Using IDW methods, MAE and RMSE decreased gradually as the amount of sampling data increased, while the correlation coefficient decreased at the same time. The increase from 50% to 90% was slow with slight fluctuations, and that from 20% to 50% became obvious. Especially, when the sampling fraction <20%, MAE and RMSE were increased significantly, and the correlation coefficient was significantly reduced. (2) We found that the relationship between uniform of stations distribution and precipitation interpolation results was complex after cross-validation, and sometimes it did not have better interpolation result under more uniform distribution of stations. (3) With not only the random sampling data, but also the same interval sampling data, MAE and RMSE using IDW methods were large than Kriging interpolation method, while R2 was smaller. It is suggested that Kriging interpolation was better than IDW method in this paper. Taking 50% and 20% sampling rate groups as an example, using IDW and Kriging spatial interpolation methods, we obtained the precipitation spatial distribution of China, and the interpolation results were consistent with the actual situation.

Cite this article

ZENG Hongwei, LI Lijuan, ZHANG Yongxuan, LIU Yumei . Study on Spatial Interpolation of Precipitation with Large Scale Samples: A Case Study on 2009’s Precipitation of China[J]. PROGRESS IN GEOGRAPHY, 2011 , 30(7) : 811 -818 . DOI: 10.11820/dlkxjz.2011.07.005


[1] 余钟波. 流域分布式水文学原理及应用. 北京: 科学出版社, 2008: 28-33.

[2] 朱会义, 贾邵凤. 降雨信息空间插值的不确定性分析. 地理科学进展, 2004, 23(2): 34-42.

[3] Xie P, Arkin P A. Global precipitation: A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs. Bulletin of the American Meteorological Society, 1997, 78(11): 2539-2558.

[4] Arpe K. The hydrological cycle in the ECMWF short range forecasts. Dynamics of Atmospheres and Oceans, 1991, 16(1-2): 33-59.

[5] New M, Todd M. Precipitation measurements and trends in the twentieth century. International Journal of Climatology, 2001, 21(15): 1889-1922.

[6] 高歌, 龚乐冰, 赵珊珊, 等. 日降水量空间插值方法研究. 应用气象学报, 2007, 18(5): 732-736.

[7] Dirks K N, Hay J E, Stow C D, et al. High-resolution studies of rainfall on Norfolk Sland. Part Ⅱ: Interpolation of rainfall data. Hydro, 1998, 208(3-4): 187-193.

[8] LAM N. Spatial interpolation methods: A review. The American Cartographer, 1983, 10(2): 129-149.

[9] 石朋, 芮孝芳. 降雨空间插值方法的比较与改进. 河海大学学报: 自然科学版, 2005, 33(4): 361-365.

[10] Dubrule O. Two methods with different objectives: Spline and Kriging. Mathematical Geology, 1983, 15(2): 245-257.

[11] Puente C E, Bras R L. Disjunctive Kriging, universal kriging, or no kriging: Small sample results with simulated fields. Mathematical Geology, 1986, 18(3): 287-205.

[12] 封志明, 杨艳昭, 丁晓强, 等. 气象要素插值方法优化. 地理研究, 2004, 23(3): 357-364.

[13] 刘登伟, 封志明, 杨艳昭. 海河流域降水空间插值方法的选取. 地球信息科学, 2006, 8(4): 75-83.

[14] Pardo-Igúzquiza E. Optimal selection of number and location of rainfall gauges for areal rainfall estimation using geostatistics and simulated annealing. Journal of Hydrology, 1998, 210(1-4): 206-220.

[15] Grimes D I F, Pardo-Igúzquiza E, Bonifacio R. Optimal areal rainfall estimation using raingauges and satellite data . Journal of Hydrology, 1999, 222(1-4): 93-108.

[16] 李海滨, 林忠辉, 刘苏峡. Kriging 方法在区域土壤水分估值中的应用. 地理研究, 2001, 20(4): 446-452.

[17] 汤国安, 杨昕. ArcGIS 地理信息系统空间分析实验教程. 北京: 科学出版社, 2006: 363-422.

[18] 刘志红, Li Lingtao, McVicar T R, 等. 专用气候数据空间插值软件ANUSPLIN 及其应用. 气象, 2008, 34(2): 92-100.

[19] 姜晓剑, 刘小军, 黄芬, 等. 逐日气象要素空间插值方法的比较. 应用生态学报, 2010, 21(3): 624-630.

[20] 李新, 程国栋, 卢玲. 空间内插方法比较. 地球科学进展, 2000, 15(3): 260-265.

[21] 林忠辉, 莫兴国, 李宏轩, 等. 中国陆地区域气象要素的空间插值. 地理学报, 2002, 57(1): 47-56.

[22] Holdaway M R. Spatial modeling and interpolation of monthly temperature using Kriging. Climate Research, 1996, 6(3): 215-225.

[23] 林云萍, 赵春生. 中国地区不同强度降水的变化趋势. 北京大学学报: 自然科学版, 2009, 45(6): 995-1002.