基于小波分析與地統(tǒng)計學的盤山縣降水量時空特征分析

曹君宇

(遼寧省水資源管理集團有限責任公司，沈陽 110003)

降水是生態(tài)環(huán)境演替的基本要素，也是人類社會發(fā)展必不可少的資源支撐。在全球變化背景下，降水時空分布模式正在潛移默化地改變[1]。作為水資源重要補給，降水的時空變化對水資源的數(shù)量、質(zhì)量及其供需平衡具有重要意義，并深刻影響作物生產(chǎn)、居民生活、水資源保護。掌握區(qū)域降水特征是水生態(tài)文明建設(shè)的基礎(chǔ)工作[2]。本研究擬從時間、空間維度對盤山縣降水量特征進行量化分析，以期為區(qū)域水資源優(yōu)化管理提供參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

盤山地處東北遼河下游，地理位置121°34′-122°29′E， 40°50′-41°27′N，面積1735km2。受西北太平洋暖濕氣流與西伯利亞高壓控制，形成溫帶大陸新季風氣候，季候分明、雨熱同期，年平均氣溫8.8℃，降水量605mm，無霜期170d。區(qū)域地形平坦開闊，平均海拔4m，發(fā)育有雙臺子河、外遼河等13條河流，另有沼澤、灘涂廣布。該縣蘊藏豐富的水資源，年平均徑流量達2.58億m3，河川徑流總量72.04億m3。

1.2 小波方差

將降水時間序列經(jīng)距平處理后，從不同尺度a對所有小波系數(shù)的平方進行積分，得到小波方差V(a):

(1)

式中：Wf為小波基函數(shù)；b為小波中心位置。

1.3 地統(tǒng)計學

地統(tǒng)計學是面向空間分布變量的一種空間統(tǒng)計方法。其中，半方差函數(shù)(Semivariogram)是地統(tǒng)計學的理論基礎(chǔ)，其認為空間上臨近為主上變量屬性呈一定相關(guān)性，并隨距離發(fā)生變化。

(2)

式中：h為空間臨近點之間距離；N(h)∈1,2……N，表示距離為h的所有臨近點數(shù)量；Z(xi)和Z(xi+h)分別為變量在位置xi、xi+h處的值，Z為空間擬合函數(shù)。

1.4 GIS空間插值法

基于點位獲取的降雨觀測資料并不能直接描述區(qū)域上降雨分布具體信息?；贕IS的空間插值法成為降雨信息空間可視化的普遍技術(shù)方法。其在ArcGIS10.5環(huán)境下的操作過程如下:Geostatistical AnalysisGeostatistical Wizard(Input data, Method )General propertiesFinish。

1.4 數(shù)據(jù)來源與處理

本研究中空間數(shù)據(jù)為地理國情監(jiān)測云平臺(http://www.dsac.cn)提供的降水量空間柵格文件。文件格式為tiff，空間分辨率達500mm，示意圖如圖1。鑒于該數(shù)據(jù)過于粗糙，為得到縣域降水量空間分布詳細信息，先將柵格像元提取至點，再運用地統(tǒng)計學方法進行空間插值?；诿绹h(huán)境系統(tǒng)研究所公司(Environmental Systems Research Institute, ESRI)開發(fā)的ArcGIS10.5軟件對柵格文件進行處理，空間分析過程參照《AicGIS空間分析教程》。

圖1 降水要素柵格文件平面信息示意圖

降水量時間數(shù)據(jù)來自國家氣象中心基礎(chǔ)數(shù)據(jù)云平臺(http://data.cma.cn/user/toLogin.html)。由于權(quán)限限制，僅獲取了1950-2017年值序列資料。通過Excel 2016軟件進行經(jīng)典統(tǒng)計分析，運用Matlab 2018的小波工具箱進行小波方差分析。

2 結(jié)果與分析

2.1 盤山縣降水量時間序列變化特征

圖2為1950-2017年盤山縣降水量變化。該時域內(nèi)降水量呈波動特征，于1964年、1985年、2010年和2012年出現(xiàn)了明顯的峰值，年降水量達965.94mm、907.18mm、892.2mm、864.2mm，高出平均水平的50%-75%；在1958年、1978年、1992年、2014年出現(xiàn)極低值，為390.5mm、335.6mm、302.8mm、353.7mm，低于平均值的20%-30%。近67a內(nèi)年降水量平均值為623.2mm，低于平均值的有32a，高于平均值的年份達35次，總體來看，盤山縣降水量異常分布相對均衡，極高極低值呈對稱分布。統(tǒng)計顯示，年降雨量標準差為139.6mm，變異系數(shù)達24.63%，表明區(qū)域降水量年際呈中等程度變異性。線性回歸宏觀反映了降水量時間變化趨勢，其傾向斜率為-1.1914 mm/a，說明近半世紀以來該地降水量呈現(xiàn)微弱減少特征，但并未達到5%水平的顯著性檢驗。

圖2 盤山縣降水量時間序列特征

2.2 盤山縣降水量時間周期性特征

降水是氣候系統(tǒng)的要素之一，其年際波動性通常是由氣象過程演變與外部干擾引起的。時間序列研究經(jīng)驗表明，降水波動性變化具有一定規(guī)律。小波分析法能削弱波動的降雨信號噪聲，從中挖掘隱含信息熵(圖3)。小波方差表明，當震蕩周期為5a時，小波方差達到極大值，為1.94，說明5a是區(qū)域降水量時間變化的主周期。此外，當周期為3a、9a、22a時，小波方差也分別達到極值1.23、1.16、0.65，這些周期是降水變化的次周期。以上表明，盤山縣降水具有長短周期疊合特征。

圖3 盤山縣降水量小波方差

2.3 盤山縣降水量空間模型

一定尺度空間上降水量具有穩(wěn)定的分布規(guī)律，即空間結(jié)構(gòu)性。針對空間變量，地統(tǒng)計學中提供了高斯、線性、球面和指數(shù)4種模型定量捕捉其結(jié)構(gòu)性規(guī)律。以柵格像元點值為數(shù)據(jù)輸入進行空間擬合，得到結(jié)果如圖4所示。不同模型的擬合的散點分布與空間結(jié)構(gòu)存在一定差異。其中高斯模型的R2為0.762，RSS為0.124，精度高于其他3種模型，表明該模型是盤山縣降水量空間擬合最佳模型。依圖可知，在高斯模型中隨著步長的增加，半方差值呈現(xiàn)增加后平穩(wěn)趨勢；同時像元空間離散度增大，導致鄰域之間的自相關(guān)性減弱。步長達到2.5km時，半方差值趨于飽和。由此表明，2.5km是縣域降水空間自相關(guān)的變程。

2.4 盤山縣降水量空間分布特征

圖5(a)為區(qū)域降水像元點值分布總體趨勢，具有明顯的空間地帶性規(guī)律。圖5(b)為區(qū)域降水量空間分布圖，圖中顯示區(qū)域降水量介于695-604mm之間，空間平均值為621.8mm，標準差為198.6mm，空間變異性達31.93%，呈中度變異性。其中降水高值區(qū)集中于東南部，低值區(qū)位于西北部，表現(xiàn)出條帶狀分布格局。這種分布規(guī)律與海陸位置存在密切耦合關(guān)系密。盤山縣屬于前海低地向內(nèi)陸延伸地帶，東南部瀕臨海洋，因而受海洋水汽活動影響較強烈；而西北部距離海洋較遠產(chǎn)生水汽運移損失，導致降雨量較小。雖然縣域空間尺度較小，但依然反映了降水在小尺度上的空間異質(zhì)性。

圖4 盤山縣降水量半方差函數(shù)特征

(a)

(b)圖5 盤山縣降水量柵格像元與空間分布

結(jié)合圖4和圖5可知，盤山縣降水量空間分布具有一定的結(jié)構(gòu)性特征。

3 小 結(jié)

基于數(shù)理統(tǒng)計與空間分析方法研究了區(qū)域降水量時空分布規(guī)律。結(jié)果顯示，盤山1950-2017呈波動下降趨勢，變化率為-0.19mm/a，但并不顯著(P>0.05)，時間序列變化存在多個震蕩周期，長周期為22a、短周期為3a、5a、9a，以短周期為主。地統(tǒng)計學方法提供了多種空間分布模型，盤山縣降雨量最佳空間擬合模型為高斯模型，具有良好的空間自相關(guān)性，空間變程為2.5km；其空間分布呈現(xiàn)自東南向西北減少的地帶性規(guī)律。該時空規(guī)律性分析能幫助理解區(qū)域降水變化特征，為水資源優(yōu)化配置提供科學基礎(chǔ)。