日降水量空间插值方法研究
Spatial Interpolation Methods of Daily Precipitation
-
摘要: 采用反距离权重法和普通克里格方法对26°~34°N, 103°~115°E范围内2004年逐日降水量进行空间插值试验分析, 分辨率为1 km×1 km。采用交叉检验方法和准确率方法对两种方法插值的总体效果及不同等级的降水插值效果进行综合对比。结果表明:两种方法插值效果近似, 插值结果与实测值相关系数分别为0.83和0.82。但对日雨量较大的情况, 两种方法插值效果均有所降低, 相关系数为0.66和0.67。两种方法的实测值与插值结果的相关系数在不同季节非常接近, 并且均以春季最大, 其次为冬、秋季, 夏季相关系数最小; 通过采用平均误差、平均绝对误差和均方根3个指标衡量及不同等级雨量的插值准确率统计比较, 普通克里格方法插值效果略好于反距离权重法。Abstract:
High-resolution precipitation field is very useful in the study of precipitation-induced geological hazards. The quality of the precipitation field varies if different interpolation methods are used. In this study inverse distance weighting and ordinary Kriging methods are applied to interpolate daily precipitation of stations over 26°—34°N, 103°—115°E into 1 km×1 km grid.The interpolation results are then compared to indicate which method better represents the spatial pattern and intensity of precipitation. Results show that the two methods perform similarly well by their correlation coefficients of 0.83 and 0.82 between observed and interpolated precipitation in the cross validation test. When daily precipitation is equal to or greater than 10 mm, both methods are less efficient and the correlation is down to 0.66 and 0.67. The correlation coefficient, as an indicator of the interpolation quality, has a clear seasonal trend with a maximum in spring and a minimum in summer. From analysis to the accuracy of interpolation for heavier precipitation, both methods show unsatisfied results and Inverse distance weighting method with more error rate. The accuracy decrease from 80% for precipitation exceeding 10 mm to 65% for precipitation exceeding 25 mm, particular in rainstorm with above 50 mm precipitation, the accuracy is only 50%. The observed precipitation is systematically underestimated by both methods and the interpolation quality gradually decreases for heavy precipitation. Generally, ordinary Kriging is better than inverse distance weighting. But the running time is a disadvantage when the ordinary Kriging method is applied into operational system, inverse distance weighting method may be substituted for it.
-
表 1 两种插值方法的交叉检验结果 (日降水量≥10 mm)
Table 1 Cross validations of the two methods under daily precipitation≥10 mm
表 2 两种方法对不同降水级别以上的资料效果检验
Table 2 The verification of the two interpolation methods for different grades of precipitation
-
[1] David T Price, McKenney Daniel W, Nalder Ian A, et al. A comparison of two statistical methods for spatial interpolation of Canadian monthly mean climate data. Agricultural and Forest Meteorology.2000, 101 : 81-94. doi: 10.1016/S0168-1923(99)00169-0 [2] Samuel S P Shen, Dzikowski Peter, Li Guilong, et al. Interpolation of 1961-97 daily temperature and precipitation data onto alberta polygons of ecodistrict and soil landscapes of Canada. J Appl Meteorol. 2001.40:2162-2177. https://www.researchgate.net/publication/253964006_Interpolation_of_1961-97_Daily_Temperature_and_Precipitation_Data_onto_Alberta_Polygons_of_Ecodistrict_and_Soil_Landscapes_of_Canada [3] Muhammad Ashraf, loftis Jim C, Hubbard K G. Application of geostatistics to evaluate partial weather station networks. Agric For Meteor, 1997, 84 : 255-271. doi: 10.1016/S0168-1923(96)02358-1 [4] 尤卫红, 夏欣健, 赵宁坤.云南逐月雨量和气温的格点数据资料场建立.云南地理环境研究, 2004, 16(1):14-18. http://www.cnki.com.cn/Article/CJFDTOTAL-YNDL200401002.htm [5] Prudhomme C, Reed D W. Mapping extreme rainfall in a mountainous region using geostatistical techniques: A case study in Scotland. Int J Climatol, 1999, 19 : 1337-1356. doi: 10.1002/(ISSN)1097-0088 [6] 李新, 程国栋, 卢玲.空间内插方法比较.地球科学进展, 2000, 15(3):260-265. http://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ200003003.htm [7] Phillips D L, Dolph J, Marks D. A comparison of geostatistieal procedures for spatial analysis of precipitation in mountainous terrain. Agric For Meteor, 1992, 58: 119-141. doi: 10.1016/0168-1923(92)90114-J [8] Daly C. A statistical-topographic model for mapping climatological precipitation over mountainous terrain. J Appl Meteor, 1994, 33: 140-158. doi: 10.1175/1520-0450(1994)033<0140:ASTMFM>2.0.CO;2 [9] Creutin J D, Obled C. Objective analysis and mapping techniques for rainfall fields: An objective comparison. Water Resour Res, 1982, 18 (2): 413-431. doi: 10.1029/WR018i002p00413 [10] 庄立伟, 王石立.东北地区逐日气象要素的空间插值方法应用研究.应用气象学报, 2003, 14(5):605-615. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20030575&flag=1 [11] 魏凤英, 曹鸿兴.我国月降水和气温网格点资料的处理和分析.气象, 1994, 20(10):26-30. http://www.cnki.com.cn/Article/CJFDTOTAL-QXXX410.004.htm [12] Dirks K N, Hay J E, Stow C D, et al. High-resolution studies of rainfall on Norfolk Sland. Part Ⅱ: Interpolation of rainfall data. J Hydrol, 1998, 208: 187-193. doi: 10.1016/S0022-1694(98)00155-3 [13] 魏凤英, 曹鸿兴.地统计学分析技术及其在气象中的适用性.气象, 2002, 28(12):3-5. http://www.cnki.com.cn/Article/CJFDTOTAL-QXXX200212000.htm [14] 王政权.地统计学及在生态学中的应用.北京:科学出版社, 1999:1-197. [15] 林忠辉, 莫兴国, 李宏轩, 等.中国陆地区域气象要素的空间插值.地理学报, 2002, 57(1):47-56. http://www.cnki.com.cn/Article/CJFDTOTAL-DLXB200201005.htm