豪雨誘發山崩的研究中，雨量為重要的促崩因子。由於雨量站的空間密度不足，導致雨量的內插結果無法反映真實情況。因此，推求準確的雨量空間分布是相當重要的課題。本研究以2004年艾利颱風事件為例進行石門水庫集水區雨量空間分布特性之分析，分別以單變量的距離平方反比法與克利金法，以及結合雷達回波使用聯合克利金法、Kriging with varying local means及Kriging with varying local means（trend）等多變量地質統計方法進行雨量空間內插，以改善雨量站密度不足的問題。
　　針對篩選過的艾利颱風雨量資料及雷達回波資料進行半變異圖分析，顯示雨量的影響範圍約43公里，雷達回波資料的影響範圍約35公里；對石門水庫集水區及其鄰近地區地表雨量站總雨量值與該點上方2000公尺高度的總雷達回波值進行迴歸分析，顯示兩者相關性甚高，R2可達0.7以上。
　　本研究並針對五種內插方法以研究區域內10個雨量站作總雨量的驗證，發現對全區而言Kriging with varying local means（trend）方法所得的整體誤差最小，距離平方反比法的誤差最大；聯合克利金法與Kriging with varying local means內插方法推估的雨量結果，能表現較細微的雨量空間分布。除了研究區南區缺乏雷達資料的區域外，對其餘地區而言使用聯合克利金法估計誤差明顯小於克利金法的估計誤差，顯示使用雷達回波資料做為輔助變數確實能改善空間內插之結果。

　　How to estimate the spatial distribution of precipitation is a critical subject in the study of storm-induced landslides. Because the density of rain stations is not adequate, the estimated result can not reflect the real feature. This study presents characteristic of the rainfall spatial distribution of Typhoon Aere in the Shihmen Reservoir watershed by using two univariate techniques (inverse square distance method and kriging) and three multivariate geostatistical algorithms (cokriging, kriging with varying local means and kriging with varying local means(trend)). The multivariate geostatistics incorporate rainfall data and radar reflectivity data in the interpolation.
　　The spatial variability of rainfall and radar data, were processed by semivariogram analysis. It shows that the effective range of rainfall is 43 km and the effective range of radar reflectivity is 35 km. Rainfall and radar reflectivity data at the elevation 2000 m shows high correlation in the study area (R2>0.7).
　　The performances of the five interpolators were validated by comparing with data of each of the ten rain station in the Shihmen Reservoir watershed, when the target station is not included in the interpolation. For the entire watershed, Kriging with varying local means(trend) method provides the smallest mean absolute error, whereas the inverse square distance method provides the largest one. The estimated precipitation map by using cokriging and krigng with varying local means show more detailed spatial variation of precipitation. For the area, besides the southern area where lack of radar reflectivity data, the estimated error is smaller by using cokriging method than by kriging apparently. It appears that radar reflectivity data could offer more detailed rainfall variability, and should be a good auxiliary variable for rainfall interpolation.

李錫堤、潘國樑、林銘郎（2003）山崩調查與危險度評估-山崩潛感分析之研究（1/3），經濟部中央地質調查所報告，第92-11號，共154頁。
李錫堤、潘國樑、林銘郎（2004）山崩調查與危險度評估-山崩潛感分析之研究（2/3），經濟部中央地質調查所報告，第93-17號，共264頁。
李錫堤、潘國樑、林銘郎（2005）山崩調查與危險度評估-山崩潛感分析之研究（3/3），經濟部中央地質調查所報告，第94-18號，共268頁。
林彥享（2003）運用類神經網路進行地震誘發山崩之潛感分析，國立中央大學應用地質研究所碩士論文，共89頁。
林清豐（1993）雨量空間分佈特性之研究，國立台灣大學土木工程研究所碩士論文，共135頁。
林淑玲（2001）宜蘭地區颱風降雨與地形、空間分布關係之探討，國立中興大學水土保持學系碩士論文，共89頁。
孫恆強（1998）雷達與地面雨量計整合估計與應用－以中正機場雷達為例，國立台灣大學土木工程研究所碩士論文，共87頁。
張智昌（2005）整合氣象雷達與即時降雨資料於颱風降雨推估之研究，國立台灣大學地理環境資源學研究所碩士論文，共90頁。
許玉金（2005）台灣北部地區雨滴粒徑分佈特性與降雨估計之探討，國立中央大學大氣物理研究所碩士論文，共93頁。
黃威雄（2000）應用類神經網路於颱風期間雷達降雨模擬之研究，國立台灣大學土木工程學研究所碩士論文，共127頁。
黃琮逢（1995）模擬雷達與地面雨量資料整合方法之研究，國立台灣大學土木工程學研究所碩士論文，共97頁。
楊政潭（2003）雷達回波應用於颱風降雨空間分佈與總量之研究－以納莉颱風為例，國立中央大學水文科學研究所碩士論文，共105頁。
葉惠中（1999）區域化變數理論與隨機變域模擬在雨量站網設計之研究，國立台灣大學農業工程學研究所博士論文，共150頁。
鄭士仁、洪君伯、謝惠紅、王育民（2005）區域化變數理論應用於降雨之時間與空間特性之研究，農業工程學報，51（4），第69-83頁。
Bent, A.E. (1943) Radar echoes from atmospheric phenomena”, MIT radiation laboratory rep., 173, 10.
Buytaert, W., Celleri, R., Willems, P, Bievre, B. D., Wyseure, G.(2006) Spatial and temporal rainfall variability in mountainous areas: A case study from the south Ecuadorian Andes, Journal of Hydrology, 329, 413-421.
Creutin, J.D., Delrieu, G., Lebel, T. (1988) Rain measurement by　raingage-radar combination: a geostatistical approach. Journal of Atmosphere and Oceanic Technology, 5, 1, 102–115.
Daly, C., Neilson, R.P., Phillips, D.L. (1994) A statistical-topographic model for mapping climatological precipitation over mountainous terrain, Journal of Applied Meteorology, 33, 140-158.
Diodato N. and Ceccarelli M.(2005) Interpolation process using multivariate geostatistics for mapping of climatological precipitation mean in the Sannio Mountains(southern Italy), Earth Surface Processes and Landforms, 30, 259-268.
Goovaerts, P. (2000) Geostatistcal approaches for incorporating elevation　into the spatial interpolation of rainfall, Journal of Hydrology, 228, 113-129.
Haberlandt, U.(2007) Geostatistical interpolation of hourly precipitation from rain gauge and radar for a large-scale extreme rainfall event, Journal of Hydrology, 332, 144-157.
Hevesi, J.A., Flint, A.L., Istok, J.D.(1992a) Precipitation estimation in mountainous terrain using multivariate geostatistics. Part I: structural analysis. Journal of Applied Meteorology, 31, 661–676.
Hevesi, J.A., Flint, A.L., Istok, J.D.(1992b) Precipitation estimation in mountainous terrain using multivariate geostatistics. Part II: Isohyetal maps. Journal of Applied Meteorology, 31, 677–688.
Joss, J., J.C. Thams and A. Waldvogel(1968) The accuracy of daily rainfall measurement by radar. Preprints 13th Radar Meterology Conf., Montreal, and Amer. Meteor. Soc., 448-451.
Journel, A. G., C. J. Huijbregts (1978) Mining geostatistics: Academic Press, New York.
Krajewski, W.F. (1987) Cokriging radar-rainfall and rain gage data, Journal of Geophysical Research, 92, D8, 9571-9580.
Marshall, J.M., and Palmer, W. McK. (1948) The distribution of raindrops with size, J. Metero., 5, 165-166.
Marshall, J.S., Langille, R.C. and Palmer, W. Mck. (1947) Measurement of rainfall by radar, J. Metero., 4, 186-192.
Martinez-Cob, A.(1996): Multivariate geostatistical analysis of evapotranspiration and precipitation in mountainous terrain. Journal of Hydrology. 174,19-35.
Teng. J.H., C.-S. Chen, T.-C. Chen Weang, and Y.-L. Chen (2000) Orographic effects on a aquall-line system over Taiwan. Mon. Wea. Rev., 128, 1123-1138.
Twomey, S.(1953): On the measurement of precipitation intensity by radar, J. Meteor. 10: 66-67.
Wilson, J. P. and Gallant, J. C. (2000) Digital terrain analysis, In: Wilson, J. P. and Gallant, J. C. (eds.) Terrain Analysis- Principles and Applications, New York, John Wiley & Sons.