廣西典型巖溶區(qū)流域降雨侵蝕力變化特征及影響因素研究

莫崇勛，藍(lán)華鯤，楊瑩，阮俞理，孫桂凱，黃亞

(1.廣西大學(xué) 土木建筑工程學(xué)院， 廣西 南寧 530004；2.工程防災(zāi)與結(jié)構(gòu)安全教育部重點(diǎn)實(shí)驗(yàn)室， 廣西 南寧 530004；3.廣西防災(zāi)減災(zāi)與工程安全重點(diǎn)實(shí)驗(yàn)室， 廣西 南寧 530004;4.北海藝術(shù)設(shè)計(jì)學(xué)院 建筑學(xué)院， 廣西 北海 536000)

0 引言

土壤侵蝕是現(xiàn)代農(nóng)業(yè)面臨的主要挑戰(zhàn)之一，也是一種重大的環(huán)境危害。降雨侵蝕力作為土壤侵蝕的重要參數(shù)之一，在估算和模擬土壤侵蝕及其對(duì)社會(huì)和環(huán)境的綜合影響方面具有重要作用。全球氣候變化正在急劇變化，特別是降水發(fā)生的變化，其對(duì)土壤侵蝕產(chǎn)生的影響亦趨復(fù)雜。氣候變化直接影響區(qū)域水文循環(huán)變化，水文循環(huán)氣象要素的年內(nèi)/年際間變化已日趨復(fù)雜[1-3]。變化環(huán)境下降雨量的研究為區(qū)域水資源的合理配置、開發(fā)、利用和災(zāi)害預(yù)警提供了一定的科學(xué)依據(jù)[4]，降雨是驅(qū)動(dòng)土壤侵蝕的重要因素[5]，降雨侵蝕力特征變化研究與區(qū)域降雨變化有著密切聯(lián)系。因此，研究變化環(huán)境下降雨量及降雨侵蝕力的變化特征具有重要意義。

當(dāng)前國(guó)內(nèi)外很多學(xué)者對(duì)降雨及降雨侵蝕力進(jìn)行了大量分析，并取得了一定的成果。提出預(yù)測(cè)土壤流失的最佳單變量是風(fēng)暴總降雨能量及其最大30強(qiáng)度的乘積，并且降雨侵蝕力是土壤方程USLE的主要影響因素之一[7]，RENARD等[8]研究發(fā)現(xiàn)降雨侵蝕力是土壤流失模型(the revised universal soil loss equation, RUSLE)的重要因子。降雨侵蝕力計(jì)算模型眾多，考慮到降雨侵蝕力的大小與降雨量、降雨強(qiáng)度、雨日數(shù)、暴雨數(shù)等降雨特征密切相關(guān)[9]，章文波等[10]提出了利用日降雨量估算降雨侵蝕力模型，與基于年、月尺度的降雨侵蝕力計(jì)算模型相比，日降雨資料易獲得，精度高，可以更好的反映降雨侵蝕力的特點(diǎn)[11]，因此得到廣泛應(yīng)用。謝坤堅(jiān)等[12]利用趨勢(shì)分析、氣候傾向率等方法研究了廣西降雨侵蝕力變化特征，發(fā)現(xiàn)大部分地區(qū)多年降雨侵蝕力日趨增加；林杰等[13]利用M-K檢驗(yàn)、重度極差法、變異系數(shù)等方法研究連云港市區(qū)域降雨侵蝕力的變化特征，結(jié)果發(fā)現(xiàn)該地區(qū)降雨侵蝕力分布集中在夏季；程金文等[14]研究了陜南地區(qū)近54 a降雨侵蝕力變化特征，結(jié)果表明陜南地區(qū)降雨侵蝕力變化規(guī)律為年內(nèi)分布差異較大且年際變化較大；周琳等[15]對(duì)四川省1955—2015年降雨侵蝕力的年內(nèi)時(shí)空變化趨勢(shì)特征進(jìn)行了分析，結(jié)果表明四川省年內(nèi)降雨侵蝕力與降雨量呈極顯著相關(guān)關(guān)系，降雨侵蝕力與降雨量的年內(nèi)分布規(guī)律基本一致。上述研究極少涉及巖溶區(qū)流域，而中國(guó)西南地區(qū)廣泛分布巖溶區(qū)，其土壤侵蝕特性與其他區(qū)域顯著不同。喀斯特地區(qū)生態(tài)系統(tǒng)敏感度高，災(zāi)變承受能力低，環(huán)境容量小，是典型的生態(tài)脆弱區(qū)[16]。澄碧河流域位于廣西壯族自治區(qū)百色市，是中國(guó)西南典型巖溶區(qū)流域。西南喀斯特地區(qū)的土壤侵蝕疊加了化學(xué)溶蝕、重力侵蝕和流水侵蝕的耦合作用，呈現(xiàn)地面流失和地下漏失的混合侵蝕機(jī)制[17-19]。鑒于此，論文基于典型巖溶區(qū)流域——澄碧河流域壩首站和平塘站1963-2017年逐日降雨量實(shí)測(cè)資料，對(duì)該流域降雨量和降雨侵蝕力時(shí)程變化、相關(guān)性及影響因素進(jìn)行分析，以期為巖溶地區(qū)有效防治土壤侵蝕和制定科學(xué)有效的水土保持方案提供理論依據(jù)。

1 研究區(qū)域概況

澄碧河發(fā)源于廣西壯族自治區(qū)凌云縣境內(nèi)，由縣城北水源洞流出，向南至彩架入洞，輾轉(zhuǎn)伏流約12 km，最終匯入澄碧河水庫(kù)，屬于珠江流域西江水系右江干流左岸一級(jí)支流。澄碧河水庫(kù)流域位于106°21′E～106°48′E，23°50′N～24°45′N，流域總面積2 087 km2，其中巖溶區(qū)面積400 km2，占流域面積19%，是中國(guó)西南典型巖溶區(qū)流域。澄碧河流域地處亞熱帶季風(fēng)氣候區(qū)，夏熱冬溫，四季分明，季風(fēng)發(fā)達(dá)。流域多年平均降雨量為1 560.0 mm，24 h最大降水量461.3 mm，降水主要集中在汛期。流域內(nèi)共設(shè)有10個(gè)氣象站和兩個(gè)水文站，其中壩首和平塘水文站記錄有1963—2017年55逐日降雨資料，10個(gè)氣象站點(diǎn)逐日降雨資料存在缺失情況。因此，壩首站和平塘站逐日降雨數(shù)據(jù)對(duì)澄碧河流域壩首區(qū)域降雨侵蝕力的研究具有良好的代表性和合理性。澄碧河流域DEM及水文站分布如圖1所示，其中框選區(qū)域?yàn)槌伪毯恿饔驂问讌^(qū)域。

圖1 澄碧河流域DEM及水文站分布

2 研究方法和數(shù)據(jù)來(lái)源

本文以澄碧河流域55逐日降雨資料為數(shù)據(jù)基礎(chǔ)，結(jié)合文獻(xiàn)[10]提出的降雨侵蝕力計(jì)算模型分析典型西南巖溶區(qū)降雨侵蝕力的變化特征及影響因素，主要從以下3個(gè)方面展開研究：①采用相對(duì)變率分析法[20]分析降雨量和降雨侵蝕力的年內(nèi)穩(wěn)定性特征；①采用趨勢(shì)分析方法和相關(guān)分析法分析[21]降雨量和降雨侵蝕力的年際變化特征；③最后利用交叉小波分析法研究降雨侵蝕力的影響因素。本章節(jié)僅介紹降雨侵蝕力模型計(jì)算方法及交叉小波分析方法。

2.1 降雨侵蝕力計(jì)算

本文降雨侵蝕力計(jì)算模型采用的是文獻(xiàn)[10]提出的日雨量計(jì)算降雨侵蝕力的簡(jiǎn)易算法模型，按照月尺度統(tǒng)計(jì)得到R值。計(jì)算公式為

(1)

(2)

α=21.586β-7.189 1，

(3)

(4)

2.2 交叉小波分析

交叉小波分析可以直觀分析氣象和水文要素之間的相干性[23]，且已有研究表明我國(guó)氣候受ENSO、太陽(yáng)活動(dòng)等大氣環(huán)流異常因子影響[24-25]。鑒于此，本文嘗試采用交叉小波分析探究多元ENSO指數(shù)(multivariate ENSO index, MEI)、海洋尼諾指數(shù)(ocean nino index, ONI)和太陽(yáng)黑子三個(gè)大氣環(huán)流異常指標(biāo)對(duì)澄碧河流域降雨侵蝕力變化的影響，文中分析所使用的交叉變換工具包括小波功率譜(cross wavelet, XWT)和小波凝聚譜(wavelet transform coherence, WTC)。

兩個(gè)時(shí)間序列x(t)和y(t)間的XWT定義為[26-28]

(5)

WTC的定義為

(6)

2.3 數(shù)據(jù)來(lái)源

本論文采用1963—2017年澄碧河流域壩首和平塘水文站逐日降雨資料；降雨侵蝕力影響分析中采用的氣候指標(biāo)太陽(yáng)黑子、MEI和ONI數(shù)據(jù)分別來(lái)源于美國(guó)航空航天局和美國(guó)大氣局物理科學(xué)實(shí)驗(yàn)室。

3 結(jié)果與分析

3.1 年內(nèi)分布特征

壩首站和平塘站多年各月降雨量與侵蝕性降雨量2，由圖2知，1963—2017年各月的降雨量和侵蝕性降雨量時(shí)間分布趨勢(shì)基本一致，降雨量和侵蝕性降雨量主要集中在5～8月，其中壩首站侵蝕性降雨量占全年侵蝕性降雨量的72.48%，平塘站侵蝕性降雨量占全年侵蝕性降雨量的75.07%。壩首站的多年平均月侵蝕性降雨量的最大值出現(xiàn)在6月，172.25 mm，平塘站的多年平均月侵蝕性降雨量的最大值亦出現(xiàn)在6月，值為234.89 mm。綜上所述，侵蝕性降雨的年內(nèi)變化十分顯著，且年內(nèi)分布不均勻，呈單峰型變化，侵蝕性降雨集中在汛期。此外，降雨量和降雨侵蝕力的相對(duì)變率可以表示降雨量和降雨侵蝕力的年內(nèi)穩(wěn)定性特征，對(duì)壩首站和平塘站月降雨和月降雨侵蝕力進(jìn)行相對(duì)變率計(jì)算分析，降雨量和降雨侵蝕力相對(duì)變率如圖3所示。從圖3(a)可以看出壩首站降雨量相對(duì)變率最大值出現(xiàn)在12月份，達(dá)84.76%，可知壩首站12月份降雨量最不穩(wěn)定，8月份相對(duì)變率最小，值為37.38%，壩首站8月份的降雨量最穩(wěn)定；從圖3(b)可以看出平塘站降雨量相對(duì)變率最大值亦是出現(xiàn)在12月份，達(dá)86.68%，可知平塘站12月份降雨量最不穩(wěn)定，6月份相對(duì)變率最小，值為35.80%，平塘站6月份的降雨量最穩(wěn)定；從圖3(c)可以看出壩首站降雨侵蝕力相對(duì)變率最大值出現(xiàn)在12月份，達(dá)149.09%，可知壩首站12月份降雨侵蝕力最不穩(wěn)定，8月份相對(duì)變率最小，值為53.25%，壩首站8月份的降雨侵蝕力最穩(wěn)定；從圖3(d)可以看出平塘站降雨侵蝕力相對(duì)變率最大值亦是出現(xiàn)在2月份，達(dá)145.45%，可知平塘站12月份降雨侵蝕力最不穩(wěn)定，7月份相對(duì)變率最小，值為50.22%，平塘站7月份的降雨侵蝕力最穩(wěn)定。

(a) 壩首站多年各月降雨量與侵蝕性降雨量

(a) 壩首站降雨量相對(duì)變率

3.2 年際變化特征

采用線性趨勢(shì)分析法對(duì)澄碧河流域壩首站和平塘站年降雨量和年降雨侵蝕力進(jìn)行趨勢(shì)分析，壩首站和平塘站降雨量降雨侵蝕力年變化如圖4所示。從圖4(a)可以看出壩首站多年平均降雨量為1 099.12 mm，年最大降雨量為1 636.50 mm，年最小降雨量為721.00 mm，壩首站降雨量年際變化總體上呈現(xiàn)下降趨勢(shì)；從圖4(b)可以看出平塘站多年平均降雨量為1 320.67 mm，年最大降雨量為2 293.90 mm，年最小侵蝕性降雨量為816.00 mm，平塘站降雨量年際變化總體上也是呈現(xiàn)下降趨勢(shì)。通過(guò)上述數(shù)據(jù)可以看出，澄碧河流域降雨侵蝕力R值和降雨量年際變化總體呈現(xiàn)下降趨勢(shì)，且降雨侵蝕力R值的年際變化更強(qiáng)烈，變化幅度更大?？傮w說(shuō)來(lái)，澄碧河流域壩首區(qū)域降雨侵蝕力和降雨量的年際波動(dòng)明顯，逐年的變化幅度較大且不均勻。從圖4(c)可以看出澄碧河流域壩首站多年平均R值為6 879.13 (MJ·mm)/(hm2·h·a)，壩首站年最大R值為14 720.19 (MJ·mm)/(hm2·h)，年最小R值為3 443.76 (MJ·mm)/(hm2·h)，壩首站降雨侵蝕力總體呈現(xiàn)下降趨勢(shì)；從圖4(d)可以看出平塘站多年平均R值為9 507.43 (MJ·mm)/(hm2·h·a)，年最大R值達(dá)到21 023.35 (MJ·mm)/(hm2·h)，年最小R值為4 796.31 (MJ·mm)/(hm2·h)，平塘站降雨侵蝕力年際變化總體上呈現(xiàn)下降趨勢(shì)。

(a) 壩首站降雨年變化

3.3 相關(guān)性

1963—2017年間，澄碧河流域壩首站多年平均降雨量為1 099.12 mm，壩首站降雨侵蝕力變化范圍為3 434.76～14 720.19 (MJ·mm)/(hm2·h)；平塘站多年平均降雨量為1 320.67 mm，平塘站降雨侵蝕力變化范圍為4 796.31～21 023.35 (MJ·mm)/(hm2·h)。圖5所示為壩首站、平塘站降雨與降雨侵蝕力相關(guān)關(guān)系。從圖5中可以看出壩首站和平塘站的多年降雨量與降雨侵蝕力之間呈顯著正相關(guān)關(guān)系，R2值分別為0.89和0.85，二者擬合效果較好。分別計(jì)算壩首站、平塘站降雨和降雨侵蝕力的皮爾遜相關(guān)系數(shù)，其中壩首站為0.923，平塘站為0.941，同樣證明澄碧河流域的降雨與降雨侵蝕力兩者極其顯著相關(guān)。

(a) 壩首站降雨與降雨侵蝕力相關(guān)關(guān)系

3.4 影響因素

3.4.1 MEI對(duì)降雨侵蝕力的影響

壩首站、平塘站降雨侵蝕力與MEI的交叉小波能量譜如圖6所示，從圖6(a)可以看出，70年代初在小波功率譜高能量區(qū)存在3.1～3.9的周期，降雨侵蝕力和MEI相關(guān)顯著，存在負(fù)相位的相關(guān)關(guān)系，1984—1990年和1992—2001年在小波功率譜高能量區(qū)存在3.5～4.3和4～5的周期，降雨侵蝕力和MEI相關(guān)顯著，且存在正相位的相關(guān)關(guān)系；從圖6(b)可以看出，壩首站1994—2006年在低能量區(qū)顯著周期為6.1～6.5，存在正相位的相關(guān)關(guān)系。從圖6(c)可以看出，70年代初在小波功率譜高能量區(qū)存在3.1～3.9的周期，降雨侵蝕力和MEI相關(guān)顯著，存在負(fù)相位的相關(guān)關(guān)系，1983—1990年和1995—2001年在小波功率譜高能量區(qū)存在3.6～4.2和3.9～5.1的周期，降雨侵蝕力和MEI相關(guān)顯著，且存在正相位的相關(guān)關(guān)系；從圖6(d)可以看出，平塘站1981—2015年在低能量區(qū)顯著周期為3.0～7.1，并存在正相位的相關(guān)關(guān)系。

(a) 壩首站降雨侵蝕力和MEI交叉小波功率譜

3.4.2 ONI對(duì)降雨侵蝕力的影響

壩首站、平塘站降雨侵蝕力與ONI的交叉小波能量譜如圖7所示，由圖7(a)可知，壩首站的降雨侵蝕力和ONI在1965—1973年存在2.2～3.8的顯著周期，存在負(fù)相位的相關(guān)關(guān)系，1984—1991年和1994—2001年分別存在3.6～4.2和5.5～6.1的顯著周期，此時(shí)間段降雨侵蝕力和ONI相關(guān)關(guān)系顯著，且為正相位相關(guān)關(guān)系；由圖7(b)所以看出，在低能量區(qū)，降雨侵蝕力和ONI在1993—2005年以6.3～7.2為顯著周期，同時(shí)1996—1998年疊加2.1～2.3，且正相位關(guān)系明顯；由圖7(c)可以看出,平塘站降雨侵蝕力和ONI具有的高能量區(qū)顯著周期2.5～3.7、3.8～4.2和6.0～7.0分布在1965—1971年、1983—1989年和1996—2002年；由圖7(d)可以看出，低能量區(qū)在1982—2015年顯著周期為3.1～7.5。

(a) 壩首站降雨侵蝕力和ONI交叉小波功率譜

3.4.3 太陽(yáng)黑子對(duì)降雨侵蝕力的影響

壩首站、平塘站降雨侵蝕力與太陽(yáng)黑子的交叉小波能量譜如圖8所示，由圖8(a)知，降雨侵蝕力和太陽(yáng)黑子在1963—2013年的顯著周期為8～10.2，且存在負(fù)相位的相關(guān)關(guān)系；在低能量區(qū)[圖8(b)]，降雨侵蝕力和太陽(yáng)黑子在1985—1997年以2.0～3.5為顯著周期，負(fù)相位關(guān)系明顯；由圖8(c)可以看出，平塘站降雨侵蝕力和太陽(yáng)黑子具有的高能量區(qū)的顯著周期為8.2～12.9分布在1963—2003年，其中1963—1985年存在周期為8.1～12.9顯著正相關(guān)關(guān)系，在1985—2003年存在周期為9.1～10.5顯著負(fù)相關(guān)關(guān)系；由圖8(d)可以看出低能量區(qū)在1994—1998年顯著周期為2.9～3.8，同時(shí)在2009—2017年疊加有3.8～4.3的顯著周期。

(a) 壩首站降雨侵蝕力和太陽(yáng)黑子交叉小波功率譜

由以上結(jié)果可以看出，降雨侵蝕力與大氣環(huán)流異常有著密切的聯(lián)系，且MEI、ONI、太陽(yáng)黑子對(duì)降雨侵蝕力的影響大致相似，其周期和相關(guān)性基本相同。對(duì)于降雨侵蝕力的影響，太陽(yáng)黑子最強(qiáng)，ONI次之，MEI最弱。

4 結(jié)論

① 壩首站和平塘站侵蝕性降雨與降雨侵蝕力的年內(nèi)變化十分顯著，且年內(nèi)分布不均勻，均為單峰型，主要集中在汛期，壩首站和平塘站的降雨侵蝕力傾向率均在6月達(dá)到最高。

② 壩首站和平塘站的降雨量和降雨侵蝕力年際變化極不均勻，存在較大的變幅，平塘站多年平均R值大于壩首站。

③ 澄碧河流域壩首站和平塘站降雨量和降雨侵蝕力之間存在顯著的正相關(guān)關(guān)系，R2值分別為0.89和0.85，皮爾遜相關(guān)系數(shù)分別為0.923和0.941。

④ MEI、ONI和太陽(yáng)黑子對(duì)降雨侵蝕力均具有一定的影響，其中太陽(yáng)黑子影響最大。在高能區(qū)，降雨侵蝕力和MEI及ONI存在正相位共振周期均為3～6，降雨侵蝕力與太陽(yáng)黑子的共振周期為8～13；在低能區(qū)，降雨侵蝕力和MEI及ONI均存在3～7.5的正相位共振周期，與太陽(yáng)黑子存在2.9～4.3的共振周期。

本文僅針對(duì)流域壩首區(qū)域展開研究，存在一定的局限性，建議進(jìn)一步研究澄碧河流域降雨侵蝕力空間分布特征，為全流域水土流失防治、水土保持建設(shè)和生態(tài)修復(fù)等工作提供科學(xué)依據(jù)。