喀斯特石漠化典型區(qū)藍(lán)水綠水轉(zhuǎn)化模擬

陳珂　楊勝天　江砥　管亞兵　黃勇奇　戴國富

摘要：由滲漏導(dǎo)致的降水大量流失是造成喀斯特石漠化地區(qū)地表水短缺的主要原因，將滲漏的難利用藍(lán)水轉(zhuǎn)化為可供植被生長利用的生態(tài)綠水是目前技術(shù)條件下提高水資源利用量的有益探索?；谒{(lán)水綠水理論和土壤-植物-大氣連續(xù)體（SPAC）原理，通過采用EcoHAT系統(tǒng)水循環(huán)模型模擬貴陽市非城鎮(zhèn)地區(qū)的藍(lán)水綠水轉(zhuǎn)化，發(fā)現(xiàn)：（1）研究區(qū)綠水占降水的份額不到50%，遠(yuǎn)低于全球65%的平均水平，具有很大的轉(zhuǎn)化利用潛力;（2）增加植被覆蓋度和增厚土層均能實現(xiàn)將滲漏藍(lán)水轉(zhuǎn)化為生態(tài)綠水，增厚土層的轉(zhuǎn)化效果更為顯著;（3）增加植被覆蓋度對減少滲漏的轉(zhuǎn)化作用存在一個明顯的閾值，當(dāng)植被覆蓋度增長到56%左右時，其減少滲漏的作用急遽減弱。鑒于研究區(qū)石漠化的典型性，本研究對類似地區(qū)提高水資源利用量具有一定參考意義。

關(guān)鍵詞：喀斯特;石漠化;藍(lán)水;綠水;生態(tài);EcoHAT;貴陽

中圖分類號：P426文獻(xiàn)標(biāo)志碼：A

Abstract：The leakage of large amount of precipitation? was the main cause of surface water shortage in karst rocky desertification areas.The conversion of the leaking hard-to-use blue water into ecological green water for vegetation growth and utilization was a useful exploration to improve water resource utilization under current technical level.Based on the Blue Water & Green Water Theory and the Soil-Plant-Atmosphere Continuum （SPAC） Principle，the conversion of blue water into green water in non-urban areas of Guiyang city was simulated with the Water Cycle Model of EcoHAT System in this study.It was found that：（1） The green water in the study area accounted for less than 50% of the precipitation，which was far below the global average of 65%，and with great potential for conversion and utilization.（2） Increasing vegetation coverage and thickening soil layer could enhance the process converting the leakage blue water into ecological green water，and the effects of thickened soil layer was more significant than those of increasing vegetation.（3） If the vegetation coverage was increased to 56% or above，which was an threshed，the effects on leakage was then limited.The research results of this paper would provide an add-value reference for improving water resources utilization in the similar rocky desertification areas.

Key words：karst;rock desertification;blue water;green water;ecology;EcoHAT;Guiyang

由歷史上植被破壞導(dǎo)致的土壤流失和巖石出露地表，是加劇喀斯特石漠化地區(qū)降水滲漏流失的直接原因[1-2]，致使當(dāng)?shù)厮Y源利用量偏低于同等降水條件的非喀斯特地區(qū)[3]。相當(dāng)一部分學(xué)者認(rèn)為，喀斯特石漠化地區(qū)屬于工程性缺水，可通過修建大中型蓄水工程來提高當(dāng)?shù)厮Y源的利用量[4-5]。實踐表明，在喀斯特地區(qū)修建大規(guī)模的蓄水工程，不僅成本巨大[6-7]，而且石灰?guī)r的滲透性對于水庫等蓄水工程也構(gòu)成潛在隱患[8-9]。為此，有必要探索適宜于喀斯特石漠化地區(qū)提高水資源利用的新途徑。

生態(tài)水文學(xué)將陸地表面的降水劃分為“綠水”和“藍(lán)水”兩部分，綠水指供給植被生長代謝用的氣態(tài)水或飽和土壤水，其中直接用于植物蒸散發(fā)的部分為綠水流，用于補(bǔ)充土壤水蓄變量的部分為綠水儲存;藍(lán)水指以液態(tài)形式供給水生生態(tài)系統(tǒng)或人類利用的另一部分降水[10-11]。研究發(fā)現(xiàn)，全球綠水約占降水總量的65%，藍(lán)水約占降水總量的35%[12]，越是濕潤的地區(qū)綠水所占降水比例越高[13]。在我國相對濕潤的西南喀斯特石漠化地區(qū)，由于植被破壞，導(dǎo)致原本應(yīng)參與植被生長代謝的大量綠水最終以藍(lán)水的形式滲漏流失。如果能夠采取適當(dāng)措施將這部分滲漏藍(lán)水轉(zhuǎn)化為可供植被生長利用的生態(tài)綠水，將有利于提高當(dāng)?shù)氐乃Y源利用量，促進(jìn)生態(tài)與經(jīng)濟(jì)可持續(xù)發(fā)展。

第17卷? 總第101期·南水北調(diào)與水利科技·2019年4月

陳珂，楊勝天，江砥，等·喀斯特石漠化典型區(qū)藍(lán)水綠水轉(zhuǎn)化模擬

我國西南地區(qū)是世界上喀斯特地貌分布比較集中、石漠化危害較為嚴(yán)重的地區(qū)[14-15]，針對喀斯特地區(qū)石漠化治理及生態(tài)恢復(fù)，眾多學(xué)者展開了廣泛而深入的研究，其中通過調(diào)節(jié)綠水分配份額，從生態(tài)恢復(fù)角度實現(xiàn)對石漠化的治理，正逐漸被重視。

Yuan與萬軍等[16-17]認(rèn)為喀斯特生態(tài)環(huán)境相對脆弱、土壤貧瘠，植被生長過度依賴于生境條件;杜睿等[18]指出，是土壤含水量的變化導(dǎo)致了植被類型的演替，可見，植物生長對土壤水分依賴程度較大;王志強(qiáng)等[19]認(rèn)為，“天然植被經(jīng)過長期演替，具有自組織能力，能夠適應(yīng)一定的土壤水分狀況”，為綠水變化與植被演替之間的相關(guān)關(guān)系研究提供了寶貴的線索;王玉娟等[20]在對貴州喀斯特典型地貌區(qū)龍里實地實驗后認(rèn)為，應(yīng)當(dāng)通過增加植被蓋度和減少裸地面積，來加強(qiáng)對當(dāng)?shù)鼐G水資源消耗利用的調(diào)控;溫志群等[21]采用EcoHAT系統(tǒng)對植被蒸散量進(jìn)行估算后，認(rèn)為通過把更多的降水轉(zhuǎn)化為綠水，將更加有利于喀斯特石漠化地區(qū)的植被生態(tài)恢復(fù);楊勝天[22]基于EcoHAT系統(tǒng)實驗結(jié)果指出，貴州省喀斯特地區(qū)由于亂砍濫伐、陡坡墾殖等人類活動加劇了水土流失和石漠化，導(dǎo)致綠水存量逐漸減少，針對“藍(lán)-綠水”轉(zhuǎn)化與管理的研究亟待開展。

綜上所述，目前學(xué)術(shù)界對綠水與喀斯特地區(qū)植被生態(tài)系統(tǒng)相互作用的研究已經(jīng)達(dá)到一定深度，對綠水減少的機(jī)理也已有初步認(rèn)識，即由于土壤貧瘠、植被破壞、巖溶易滲透性等原因，導(dǎo)致盡管我國西南喀斯特地區(qū)降水豐富，卻大量滲漏流失。因此，當(dāng)前面臨的問題主要是如何將當(dāng)?shù)仉y以直接利用的滲漏藍(lán)水轉(zhuǎn)化為可供植被生長代謝用的生態(tài)用綠水、以及轉(zhuǎn)化多少。

基于此，本文擬開展以下方面的研究。

（1）通過對藍(lán)水綠水轉(zhuǎn)化進(jìn)行模擬并結(jié)合觀測數(shù)據(jù)，初步明晰研究區(qū)藍(lán)水與綠水的總量及空間分布。

（2）探索將喀斯特石漠化典型區(qū)滲漏藍(lán)水轉(zhuǎn)化為生態(tài)綠水的可行措施。

（3）基于模型計算，分析不同“藍(lán)-綠”水轉(zhuǎn)化措施的效果。

1 研究區(qū)概況

貴陽市是貴州省的省會，總面積約8 034 km2，地處貴州中部，位于云貴高原東北部長江與珠江的分水嶺區(qū)域，在東經(jīng)106°07′-107°17′、北緯26°11′-27°22′之間，喀斯特面積約為6 830.26 km2，超過全市總面積的80%[23-24]。由于歷史上大煉鋼鐵及開荒種地引發(fā)的亂砍濫伐，造成當(dāng)?shù)靥烊涣滞嘶瘒?yán)重，致使約85%的地區(qū)出現(xiàn)不同程度的石漠化[25]。貴陽市多年平均降水量為1 100 mm，2013年森林覆蓋率約為44.2%，鑒于當(dāng)?shù)厮亮魇?yán)重、地表石多土少，是典型的喀斯特石漠化地區(qū)[26-27]，本文選取貴陽市非城鎮(zhèn)地區(qū)（面積約7 495 km2）作為研究區(qū)。

2 研究方法

2.1 數(shù)據(jù)來源

本文對藍(lán)水綠水轉(zhuǎn)化過程進(jìn)行模擬的數(shù)據(jù)按屬性分為基礎(chǔ)地理空間數(shù)據(jù)和氣象水文數(shù)據(jù)，基礎(chǔ)地理空間數(shù)據(jù)主要通過MODIS遙感產(chǎn)品獲取，包括每日反照率、每日地表溫度、植被覆蓋度、土地利用、地表輻射、土壤質(zhì)地等數(shù)據(jù)。根據(jù)模型計算需要和遙感數(shù)據(jù)獲取情況，空間分辨率采用1 km精度;由于2003年之前研究區(qū)的部分MODIS產(chǎn)品數(shù)據(jù)缺失，而本研究于2015年啟動時，2014年之后的部分遙感數(shù)據(jù)尚不能下載，并且一部分作為驗證用的2014年之后的水文觀測數(shù)據(jù)，官方網(wǎng)站也尚未發(fā)布，因此選取2003年作為研究的初始年，2013年作為現(xiàn)狀年?？紤]到作為重要計算指標(biāo)的植被蓋度，在區(qū)域尺度上的年度變化并不明顯，為便于模型計算的比較分析，僅選取了2003年、2007年和2013年的實際數(shù)據(jù)參與計算。

2.2 藍(lán)水綠水轉(zhuǎn)化模擬

根據(jù)“土壤-植物-大氣連續(xù)體（Soil-Plant-Atmosphere Continuum）”即SPAC原理[28-29]，綠水可以通過對區(qū)域蒸散發(fā)量的估算求得，而藍(lán)水亦可通過將降水量減去綠水量獲得。鑒于北京師范大學(xué)生態(tài)水文研究團(tuán)隊開發(fā)的EcoHAT系統(tǒng)，此前已多次應(yīng)用于貴州喀斯特石漠化地區(qū)的蒸散發(fā)量計算[21，30-31]，本文采用EcoHAT系統(tǒng)水循環(huán)模型，對研究區(qū)藍(lán)水與綠水轉(zhuǎn)化進(jìn)行模擬。該模型主要涉及蒸散發(fā)和土壤水運(yùn)移的運(yùn)算，分別介紹如下。

2.2.1 潛在蒸散發(fā)計算

EcoHAT系統(tǒng)水循環(huán)模型對潛在蒸散發(fā)量的計算采用Priestley-Taylor公式：

模型計算結(jié)果顯示，盡管初始年和現(xiàn)狀年降水量相差達(dá)40 mm，相差幅度為4.8%，但兩個年份綠水總量構(gòu)成之間基本沒有太大的差異，變化幅度為0.6%，體現(xiàn)了研究區(qū)藍(lán)綠水構(gòu)成的相對穩(wěn)定性。在植被蓋度增加22.5%的情況下，植被蒸騰量和截留量分別增長了11.69%和14.2%（參見表1），這與植被對降水截留與蒸騰等的作用特點相符;由于蒸騰和截留的增加導(dǎo)致土壤蒸發(fā)量、蓄變量、徑流量和滲漏量的相應(yīng)減少，也與水量平衡原理基本相符，這在一定程度反映了EcoHAT系統(tǒng)模型計算結(jié)果的合理性。初始年與現(xiàn)狀年藍(lán)綠水各分量值變化情況參見圖1。

從藍(lán)水與綠水的構(gòu)成來看，2003年藍(lán)水總量為506.87 mm，綠水總量為461.36 mm，藍(lán)水與綠水占當(dāng)年降水量的比例分別為52.35%和47.65%;2013年藍(lán)水總量為463.38 mm，綠水總量為458.68 mm，藍(lán)水與綠水占當(dāng)年降水量的比例分別為50.25%和49.75%?？梢姡跏寄旰同F(xiàn)狀年的綠水模擬值占降水量的比例均不到50%，遠(yuǎn)低于全球65%的平均水平;2003年和2013年的滲漏量占降水量的比例分別8.60%和7.61%，表明在將難利用的滲漏藍(lán)水轉(zhuǎn)化為對植被生長有用的生態(tài)綠水（蒸散發(fā)和土壤水蓄變量）方面，有較大可調(diào)控空間。

為便于直觀分析研究區(qū)藍(lán)水與綠水分布，基于模型計算結(jié)果，分別生成研究區(qū)初始年與現(xiàn)狀年滲漏量分布圖，以及相應(yīng)的植被蓋度圖和土壤質(zhì)地圖。如圖2所示，初始年與現(xiàn)狀年的滲漏地帶空間分布基本一致;滲漏地帶分布與植被覆蓋在空間上的分布并無明顯相關(guān)性，從側(cè)面反映出植被退化并不是導(dǎo)致滲漏的根本原因;嚴(yán)重滲漏區(qū)域（滲漏量大于180 mm）均位于土壤水傳導(dǎo)率相對高的壤土區(qū)內(nèi)，但壤土區(qū)并不都是滲漏嚴(yán)重的區(qū)域，表明滲漏是否嚴(yán)重應(yīng)與土壤質(zhì)地有一定關(guān)系，但也受到其他因素的影響。

3.2 模型計算結(jié)果驗證

對于模型各項計算結(jié)果數(shù)值的驗證，除藍(lán)水滲漏量的驗證，由于缺少觀測資料難于直接驗證外，其余各值均可通過查閱相關(guān)文獻(xiàn)資料獲取驗證值。因此，對滲漏量可采取間接驗證方式，通過驗證參與計算的各水文分量的值，再根據(jù)水量平衡原理，達(dá)到驗證滲漏量計算精度的目的。

徑流量驗證數(shù)據(jù)可通過貴陽市水務(wù)局發(fā)布的《水資源公報》獲得，植被降水截留量的驗證值來源于張喜等[45]在貴陽市開陽縣的實地觀測數(shù)據(jù);蒸騰量、土壤蒸發(fā)量和飽和水含量的驗證值根據(jù)張志才等[46]在同為喀斯特地區(qū)且與貴陽相鄰的普定縣陳旗小流域利用Perman-Monteith方法計算的結(jié)果，土壤水蓄變量驗證值根據(jù)溫志群等[21]在貴陽及周邊地區(qū)計算的土壤水蓄變率獲得。模型模擬值與驗證值參見表3。

表3中數(shù)值顯示，徑流量的模擬值和驗證值之間相差較大，原因是驗證值417.9 mm來自2013年水資源公報數(shù)據(jù)，而同年貴陽市統(tǒng)計年鑒顯示，徑流量數(shù)值為50.21億m3，換算為徑流量應(yīng)為624.9 mm，兩者相差高達(dá)200 mm。由于統(tǒng)計年鑒數(shù)值與水資源公報數(shù)值同樣來源于水行政主管部門，這種差異應(yīng)該產(chǎn)生于統(tǒng)計上報過程中，可歸類為觀測誤差，不影響對模型計算結(jié)果的可信度評價。其次，土壤水蓄變量的驗證值比模擬值稍偏大，根據(jù)溫志群等的解釋，主要原因為其觀測時段從9月至次年3月，該時段的平均次降雨雨強(qiáng)較小，降雨歷時相對長、氣溫低，植被蒸散發(fā)較弱，處于植被生長緩慢期，因此降水大多用于補(bǔ)給綠水儲存部分，導(dǎo)致土壤含水量增長較快，土壤水蓄變量也較其他時段要高;而模擬值計算的是全年的平均數(shù)，自然偏低于驗證值[20]。另外，飽和含水量的驗證值也偏高于模擬值，原因應(yīng)與驗證值來自對小流域的觀測，而模擬值由較大區(qū)域數(shù)值平均后得到有關(guān)。其他變量的模擬值和驗證值之間的差異不是很明顯，顯示由EcoHAT模型計算的各水文分量值與驗證值之間具有良好的相似度。模型模擬值與驗證值的相似度參見圖3。

圖3中模擬值的分布均接近于45°線，顯示由EcoHAT系統(tǒng)計算而得的蒸散發(fā)量及土壤水運(yùn)移量與驗證值之間具有良好相似度，根據(jù)水量平衡原理，滲漏量模擬值也應(yīng)該具有良好的可信度。

3.3 藍(lán)水綠水轉(zhuǎn)化模擬

根據(jù)SPAC原理，在水循環(huán)的三大環(huán)節(jié)中，以現(xiàn)有技術(shù)而言，人類難以改變大氣運(yùn)動或降水過程，卻可以通過改變植被覆蓋度和土壤厚度等地表因素來影響降水的分配，從而實現(xiàn)將喀斯特石漠化地區(qū)滲漏的難利用藍(lán)水轉(zhuǎn)化為可供植被代謝利用的生態(tài)綠水。

3.3.1 增加植被對“藍(lán)-綠”水轉(zhuǎn)化的模擬分析

鑒于在土壤條件不變的情況下，能夠?qū)﹃懨嫠h(huán)產(chǎn)生影響的另外兩大因素分別為大氣降水和植被變化?？杉俣ń邓坎蛔儯M植被變化下的滲漏量變化。利用MODIS數(shù)據(jù)反演植被蓋度發(fā)現(xiàn)，2003年至2013年間，研究區(qū)植被蓋度增長幅度接近10%;為便于比較，在進(jìn)行模擬時，按植被蓋度每增長5%作為一個時間間隔，由于2007年研究區(qū)植被蓋度剛好在2003年的基礎(chǔ)上增長約5%，將2007年數(shù)據(jù)納入模型計算。同時，為分析植被蓋度小幅變化對滲漏量的影響，在2013年基礎(chǔ)上，增加模擬植被蓋度增長1%和2%的影響，模擬結(jié)果參見表4。

通過對表4的分析發(fā)現(xiàn)：

（1）在保持降水量不變的條件下，隨著植被蓋度的增加，徑流量不斷減少，綠水總量保持相同數(shù)值的增長，表明在植被增長模式下，綠水的增長來源于藍(lán)水徑流量的轉(zhuǎn)化;

（2）隨著植被蓋度的增加，研究區(qū)滲漏量呈不斷減少的趨勢，盡管減少的絕對值與綠水總量的增長相差懸殊，但仍表明植被增長能夠在一定程度上將滲漏藍(lán)水轉(zhuǎn)化為生態(tài)綠水;

（3）當(dāng)植被蓋度由50%增長到56%左右時，滲漏量減少的幅度最大（1.19%），約為1 mm/km2;

（4）根據(jù)表3中數(shù)據(jù)，按植被蓋度每增長5%計算所對應(yīng)滲漏減少的幅度，匯總?cè)绫?所示。發(fā)現(xiàn)當(dāng)植被蓋度由50%增長至56%時，滲漏量減幅較大，此后滲漏量減幅顯著下降，且下降幅度越來越小。

通過對植被蓋度增加與滲漏量減少幅度變化的相關(guān)性分析發(fā)現(xiàn)，植被蓋度增長與滲漏量變化之間呈明顯負(fù)相關(guān)的冪指數(shù)關(guān)系，且植被蓋度接近56%時，滲漏量減少的拐點作用表現(xiàn)明顯（參見圖4）。表明降水不變時，植被蓋度增加對減少滲漏的作用可能存在一個閾值，植被蓋度增加越接近這一閾值，滲漏量減少幅度變得越小，因此需要借助其他更有效的途徑來減少滲漏量。

3.3.2 增加土壤厚度對藍(lán)-綠水轉(zhuǎn)化的模擬分析

根據(jù)SPAC原理，假定降水和植被蓋度不變，選取便于模型進(jìn)行計算與分析的土壤層厚度數(shù)值，在原有40 cm平均土層厚度基礎(chǔ)上，分別增加10 cm和20 cm，對藍(lán)水綠水變化情況進(jìn)行模擬，經(jīng)計算得到各主要指標(biāo)值見表6。

通過分析表6的模擬結(jié)果發(fā)現(xiàn)：

（1）隨著土層的增厚，土壤蒸發(fā)量不斷升高，土層增厚20 cm所增加的土壤蒸發(fā)量幾乎相當(dāng)于土層增厚10 cm導(dǎo)致增加的蒸發(fā)量的2倍，表明土壤增加的厚度與所增加的蒸發(fā)量的增量之間呈較好的正相關(guān)關(guān)系;

（2）當(dāng)土壤增厚10 cm時，土壤水蓄變量所增加的45.81 mm水量剛好與滲漏減少的量相等，從水量平衡的角度來看，表明增厚土層所減少的滲漏藍(lán)水已全部轉(zhuǎn)化為綠水儲存量，此時徑流量所減少的53.22 mm水量中，45.81 mm來自于滲漏量的減少，而另外7.41 mm用于土壤水蒸發(fā)，所減少的藍(lán)水總量全部轉(zhuǎn)化為綠水;

（3）當(dāng)土壤增厚20 cm后，導(dǎo)致減少的114.82 mm徑流總量中，有100.11 mm轉(zhuǎn)化[HJ2.1mm]為土壤水蓄變量，另外14.71 mm轉(zhuǎn)化為土壤水蒸發(fā)，土壤水蓄變量所增加的水量除一部分來自于減少的64.47 mm滲漏量外，剩余部分應(yīng)來自于地表徑流量的減少;

（4）當(dāng)土壤在原有基礎(chǔ)上增厚20 cm后，導(dǎo)致除滲漏藍(lán)水外，高達(dá)114.82 mm的徑流量減少。盡管根據(jù)Tennant的研究，徑流量減少在40%以內(nèi)，不會對水生生態(tài)系統(tǒng)產(chǎn)生明顯不利影響[47]，但考慮到土壤增厚20 cm后，綠水占降水的份額已達(dá)61.57%，接近全球65%的平均水平，并且總徑流量減少幅度已高達(dá)23.58%，可能會對生態(tài)環(huán)境產(chǎn)生一定負(fù)面作用，在對進(jìn)行生態(tài)評估之前，不建議再增厚土層。

藍(lán)水綠水變化對增厚土層的響應(yīng)關(guān)系參見圖5。

圖5顯示，滲漏量與土壤厚度變化呈較顯著的負(fù)指數(shù)相關(guān)關(guān)系，同時兩者之間變化也呈稍弱的線性負(fù)相關(guān)關(guān)系，原因應(yīng)與作為滲漏藍(lán)水轉(zhuǎn)化對象的土壤水蓄變量同時受到蓄變能力增強(qiáng)和徑流量減少雙重影響有關(guān)。

3.4 藍(lán)水綠水轉(zhuǎn)化效果對比分析

此前的模型計算結(jié)果顯示，增加植被蓋度和增厚土層均能實現(xiàn)將喀斯特石漠化地區(qū)滲漏的難利用藍(lán)水轉(zhuǎn)化為生態(tài)綠水，但哪一種方法轉(zhuǎn)化效果更明顯，尚需進(jìn)一步比較分析。為此，假定降水量不變，將兩種方法對轉(zhuǎn)化綠水的貢獻(xiàn)數(shù)值列表見表7，為便于比較，對增厚土層的模擬，增加了45 cm和55 cm兩種土壤厚度。

通過表7的數(shù)值對比發(fā)現(xiàn)，隨著植被蓋度和土壤厚度的不斷增加，綠水總量均呈快速增長的趨勢，增厚土層所轉(zhuǎn)化的綠水量明顯多于增加植被蓋度的轉(zhuǎn)化量。

4 結(jié)論

本文基于藍(lán)水綠水理論，采用EcoHAT系統(tǒng)水循環(huán)模型對喀斯特石漠化典型地區(qū)—貴陽市非城鎮(zhèn)區(qū)的藍(lán)水和綠水變化進(jìn)行了模擬，分析了增加植被和增厚土層兩種方案對將研究區(qū)滲漏的難利用藍(lán)水轉(zhuǎn)化為可供植被生長利用的生態(tài)綠水的不同效果，研究結(jié)論如下。

（1）研究區(qū)在初始年和現(xiàn)狀年綠水占降水的份額分別為47.65%和49.75%，均遠(yuǎn)低于全球65%的平均水平，且2003年和2013年的滲漏量占降水量的比例分別8.60%和7.61%，具有將滲漏藍(lán)水轉(zhuǎn)化為生態(tài)綠水的巨大潛力。

（2）增加植被覆蓋度和增厚土層是實現(xiàn)將滲漏藍(lán)水轉(zhuǎn)化為生態(tài)綠水的有益途徑，增厚土層的轉(zhuǎn)化效果較增加植被更為顯著。

（3）增加植被覆蓋度對減少滲漏的轉(zhuǎn)化作用存在一個明顯的閾值，即在植被覆蓋度增長到56%左右后，增加植被對減少滲漏的作用急遽減弱。

探索將喀斯特石漠化地區(qū)滲漏的難利用藍(lán)水轉(zhuǎn)化為可供植被生長利用的生態(tài)綠水，對類似地區(qū)提高水資源利用量的研究，將具有一定的參考意義。

參考文獻(xiàn)（References）：

[1] GUO F，JIANG G，YUAN D et al.，Evolution of major environmental geological problems in karst areas of Southwestern China[J].Environmental Earth Sciences，2013，69（7）：2427-2435.DOI：10.1007/s12665-012-2070-8.

[2] BALLESTEROS D，MALARD A，JEANNIN P et al.KARSYS hydrogeological 3D modeling of alpine karst aquifers developed in geologically complex areas：Picos de Europa National Park （Spain）[J].Environmental Earth Sciences，2015，74（12）：7699-7714.DOI：10.1007/s12665-015-4712-0.

[3] 張軍以，王臘春，馬小雪等.西南巖溶地區(qū)地下水污染及防治途徑[J].水土保持通報，2014（2）：245-249.（ZHANG J Y，WANG L C，MA X X，et al.Groundwater pollution and controlling measures in karst mountainous areas of Southwestern China[J].Bulletin of Soil and Water Conservation，2014（2）：245-249.（in Chinese）） DOI：10.13961/j.cnki.stbctb.2014.02.051.

[4] QIN X，JIANG Z.Situation and comprehensive treatment strategy of drought in karst mountain areas of southwest China[Z].Springer，2011：383-389.DOI：10.1007/978-3-642-19902-8_45.

[5] PARISE M，CLOSSON D，GUTIERREZ F et al.，Anticipating and managing engineering problems in the complex karst environment[J].Environmental Earth Sciences，2015，74（12）：7823-7835.DOI：10.1007/s12665-015-4647-5.

[6] QIN L，BAI X，WANG S，et al.Major problems and solutions on surface water resource utilisation in karst mountainous areas[J].Agricultural Water Management，2015，159：55-65.DOI：10.1016/j.agwat.2015.05.024.

[7] 朱生亮，張建利，吳克華，等.巖溶工程性缺水區(qū)農(nóng)村飲用儲存水凈化方法[J].長江科學(xué)院院報，2013（11）：20-23.（ZHU S L ，ZHANG J L，WU K H，et al.，Device of purifying storage drinking water in rural karst engineering water shortage area[J].Journal of Yangtze River Scientific Research Institute，2013（11）：20-23.（in Chinese）） DOI：10.3969/j.issn.1001-5485.2013.11.005.

[8] LU Y，LIU Q，ZHANG F.Environmental characteristics of karst in China and their effect on engineering[J].Carbonates and Evaporites，2013，28（1-2SI）：251-258.DOI：10.1007/s13146-013-0158-1.

[9] 朱文孝，李坡，賀衛(wèi)，等.貴州喀斯特山區(qū)工程性缺水解決的出路與關(guān)鍵科技問題[J].貴州科學(xué)，2006（1）：1-7.（ZHU W X，LI P，HE W，et al.Solutions and major scientific problems on engineering water shortage in karst mountain area of Guizhou[J].Guizhou Science，2006（1）：1-7.（in Chinese）） DOI：10.3969/j.issn.1003-6563.2006.01.001.

[10] [ZK（#]FALKENMARK M.Coping with water scarcity under rapid population growth[C].Proceedings of Conference of SADC Ministers，Pretoria，South Africa，1995.

[11] FALKENMARK M，ROCKSTROM J.The new blue and green water paradigm：Breaking new ground for water resources planning and management[J].Journal of Water Resources Planning and Management-Asce，2006，132（3）：129-132.DOI：10.1061/（ASCE）0733-9496（2006）132：3（129）.

[12] LATHUILLIERE M J，COE M T，JOHNSON M S.A review of green- and blue-water resources and their trade-offs for future agricultural production in the Amazon Basin：what could irrigated agriculture mean for Amazonia[J].Hydrology and Earth System Sciences，2016，20（6）：2179-2194.DOI：10.5194/hess-20-2179-2016.

[13] STEWART B A，PETERSON G A.Managing green water in dryland agriculture[J].Agronomy Journal，2015，107（4）：1544-1553.DOI：10.5194/hess-2016-71.

[14] HUANG Q，CAI Y，XING X.Rocky Desertification，antidesertification，and sustainable development in the karst mountain region of Southwest China[J].Ambio：A Journal of the Human Environment，2008，37（5）：390-392.DOI：10.1579/08-s-493.1.

[15] 楊勝天，宋敘言，徐明峰，等.扶貧攻堅V4.0版：貴州省臺江縣脫貧致富的智力“流”模式[J].貴州師范大學(xué)學(xué)報（自然科學(xué)版），2018，36（4）：7-12.（YANG S T，SONG S Y，XU M F，et al.，4.0 Version of poverty alleviation：Study on intelligent“flow”of poverty alleviation in Taijiang county[J].Journal of Guizhou Normal University （ Natural Sciences），2018，36（4）：7-12.（in Chinese）） DOI：10.16614 /j.gznuj.zrb.2018.04.002.

[16] YUAN D.On the karst ecosystem[J].Acta Geologica Sinica （English Edition），2001，75（3）：336-338.DOI：10.1111/j.1755-6724.2001.tb00541.x.

[17] 萬軍，蔡運(yùn)龍.喀斯特生態(tài)脆弱區(qū)的土地退化及生態(tài)重建——以貴州省關(guān)嶺縣為例[J].中國人口.資源與環(huán)境，2003（2）：55-59.（WAN J，CAI Y L.Land degradation and eco-reconstruction in fragile karst ecosystem：The Caseof Guanling County，Guizhou[J].China Population，Resources and Environment，2003（2）：55-59.（in Chinese）） DOI：10.3969/j.issn.1002-2104.2003.02.011.

[18] 杜睿，周宇光，王庚辰，等.土壤水分對溫帶典型草地N_2O排放過程的影響[J].自然科學(xué)進(jìn)展，2003（09）：45-51.（DU R，ZHOU Y G，WANG G C，et al.，Effects of soil moisture on N_2O emission from temperate typical grassland[J].Natural science progress，2003（9）：45-51.（in Chinese））DOI：10.3321/j.issn：1002-008X.2003.09.008.

[19] 王志強(qiáng)，劉寶元，王曉蘭.黃土高原半干旱區(qū)天然錦雞兒灌叢對土壤水分的影響[J].地理研究，2005（01）：113-120.（WANG Z Q，LIU B Y，WANG X L.Effects of natural caragana shrubs on soil moisture in semi-arid region of the Loess Plateau[J].Geographic Research，2005（1）：113-120.（in Chinese））DOI：10.3321/j.issn：1000-0585.2005.01.013.

[20] 王玉娟，杜迪，楊勝天，等.貴州龍里典型喀斯特地區(qū)綠水資源耗用研究[J].中國巖溶，2008（4）：340-346.（WANG Y J，DU D，YANG S T，et al.，Study on the consumption of green water resources in the typical karst area of Longli，Guizhou[J].Chinese Karst，2008（4）：340-346.（in Chinese））DOI：10.3969/j.issn.1001-4810.2008.04.008.

[21] 溫志群，楊勝天，宋文龍，等.典型喀斯特植被類型條件下綠水循環(huán)過程數(shù)值模擬[J].地理研究，2010（10）：1841-1852.（WEN Z Q，YANG S T，SONG W L，et al.The numerical simulation on green water cycle of typical vegetation types in karst area[J].Geographical Research，2010（10）：1841-1852.（in Chinese）） DOI：10.11821/yj2010100012.

[22] 楊勝天.喀斯特地區(qū)綠水利用研究[M].北京：科學(xué)出版社，2014，100-106.（YANG S T.Study on the utilization of green water in karst area[M].Beijing：Science Press，2014，100-106.（in Chinese））

[23] GAO X，PENG S，WANG W，et al.，Spatial and temporal distribution characteristics of reference evapotranspiration trends in Karst area：A case study in Guizhou Province，China[J].Meteorology and Atmospheric Physics，2016，128（5）：677-688.DOI：10.1007/s00703-016-0442-4.

[24] 王桂萍，韓堂松，朱華等.貴陽地區(qū)石漠化植被恢復(fù)造林樹種的篩選[J].林業(yè)科技，2012（4）：37-41.（WANG G P，HAN T S，ZHU H，et al.，Study on afforestation tree species selection for rocky desertification vegetation restoration in Guiyang[J].Forestry Science & Technology，2012（4）：37-41.（in Chinese）） DOI：10.3969/j.issn.1001-9499.2012.04.014.

[25] LIU Y，HUANG X，YANG H，et al.，Environmental effects of land-use/cover change caused by urbanization and policies in Southwest China Karst area–A case study of Guiyang[J].Habitat International，2014，44：339-348.DOI：10.1016/j.habitatint.2014.07.009.

[26] WANG S J，LIU Q M，ZHANG D F.Karst rocky desertification in southwestern China：Geomorphology，landuse，impact and rehabilitation[J].Land Degradation & Development，2004，15（2）：115-121.DOI：10.1002/ldr.592.

[27] YANG P，TANG Y，ZHOU N，et al.，Characteristics of red clay creep in karst caves and loss leakage of soil in the karst rocky desertification area of Puding County，Guizhou，China[J].Environmental Earth Sciences，2011，63（3）：543-549.DOI：10.1007/s12665-010-0721-1.

[28] KANG S Z，ZHANG F C，ZHANG J H.A simulation model of water dynamics in winter wheat field and its application in a semiarid region[J].Agricultural Water Management，2001，49（2）：115-129.DOI：10.1016/S0378-3774（00）00137-2.

[29] KONRAD W，ROTH-NEBELSICK A.Integrating Plant Gas Exchange，Soil，and Hydrological Parameters in an Analytical Model：Potential Use and Limitations[J].Vadose Zone Journal，2011，10（4）：1196-1204.DOI：10.2136/vzj2010.0157.

[30] 王玉娟，杜迪，楊勝天，等.貴州龍里典型喀斯特地區(qū)綠水資源耗用研究[J].中國巖溶，2008（4）：340-346.（WANG Y J，DU D，YANG S T，et al.，Consumption of green water resources in typical karst area-A case study in Longli，Guizhou[J].CARSOLOGICA SINICA，2008（04）：340-346.（in Chinese）） DOI：10.3969/j.issn.1001-4810.2008.04.008.

[31] 楊勝天.遙感水文數(shù)字實驗 ：EcoHAT使用手冊[M].北京：科學(xué)出版社，2015.（YANG S T ，Digital experiment for remeto sensing application in hydrology-EcoHAT Handbook[M].Beijing：Science Press，2015 （in Chinese））

[32] BARELLA-ORTIZ A，POLCHER J，TUZET A，et al.，Potential evaporation estimation through an unstressed surface-energy balance and its sensitivity to climate change[J].Hydrology and Earth System Sciences，2013，17（11）：4625-4639.DOI：10.5194/hess-17-4625-2013.

[33] CLULOW A D，EVERSON C S，MENGISTU M G，et al.，Extending periodic eddy covariance latent heat fluxes through tree sap-flow measurements to estimate long-term total evaporation in a peat swamp forest[J].Hydrology and Earth System Sciences，2015，19（5）：2513-2534.DOI：10.5194/hess-19-2513-2015.

[34] MA N，ZHANG Y，SZILAGYI J，et al.，Evaluating the complementary relationship of evapotranspiration in the alpine steppe of the Tibetan Plateau[J].Water Resources Research，2015，51（2）：1069-1083.DOI：10.1002/2014WR015493.

[35] PRIESTLEY C，TAYLOR R J.Assessment of surface heat-flux and evaporation using large-scale parameters[J].Monthly Weather Review，1972，100（2）：81.DOI：10.1175/1520-0493（1972）100%3C0081：OTAOSH%3E2.3.CO;2.

[36] 趙玲玲，王中根，夏軍，等.Priestley-Taylor公式的改進(jìn)及其在互補(bǔ)蒸散模型中的應(yīng)用[J].地理科學(xué)進(jìn)展，2011（7）：805-810.（ZHAO L L，WANG Z G，et al.，Improved priestley-taylor method and its application in complementary relationship evapotranspiration model[J].Progress in Geography，2011（7）：805-810.（in Chinese）） DOI：10.11820/dlkxjz.2011.07.004.

[37] RITCHIE J T，HANKS J.Modeling plant and soil systems[J].Modeling Plant and Soil Systems，1991.DOI：10.2134/agronmonogr31.

[38] ASTON A R.Rainfall interception by eight small trees[J].Journal of Hydrology，1979，42（3）：383-396.DOI：10.1016/0022-1694（79）90057-X.

[39] KERKIDES P，KARGAS G，ARGYROKASTRITIS I.The effect of different methods used for hysteretic K（H） determination on the infiltration simulations[J].Irrigation and Drainage，2006，55（4）：403-418.DOI：10.1002/ird.266.

[40] 李保國.農(nóng)田土壤水的動態(tài)模型及應(yīng)用[M].北京：科學(xué)出版社，2000.（LI B G.Dynamic model and application of farmland soil water[J].Beijing：Science Press，2015 （in Chinese））

[41] VAN GENUCHTEN M T.A closed-form equation for predicting the hydraulic conductivity of unsaturated[J].Soil Science Society of America Journal，1980，44（5）：892-898.DOI：10.2136/sssaj1980.03615995004400050002x.

[42] 吳擎龍，雷志棟，楊詩秀.壓力入滲儀測定導(dǎo)水率的理論及其應(yīng)用[J].水利學(xué)報，1996（2）：56-62.（WU Q L，LEI Z D，YANG S X，et al.，Theory and application of pressure infiltration instrument for measuring hydraulic conductivity[J].Journal of Hydraulic Engineering，1996（2）：56-62.（in Chinese）） DOI：10.13243/j.cnki.slxb.1996.02.009.

[43] VEREECKEN H，MAES J，F(xiàn)EYEN J et al.，Estimating the soil moisture retention characteristic from texture，Bulk Density，and Carbon Content[J].Soil Science，1989，148（6）：389-403.DOI：10.1097/00010694-198912000-00001.

[44] 吳擎龍.田間騰發(fā)條件下水熱遷移數(shù)值模擬的研究[D].北京：清華大學(xué)，1993.（WU Q L.Numerical simulation of hydrothermal migration under field transpiration and evaporation conditions[D].Beijing：Tsinghua University，1993（in Chinese））

[45] 張喜，薛建輝，生原喜久雄，等.黔中山地喀斯特森林的水文學(xué)過程和養(yǎng)分動態(tài)[J].植物生態(tài)學(xué)報，2007（5）：757-768.（ZHANG X，XUE J H，HIROHIRO INOGUCHI，et al.，Distributed hydrological model and eco-hydrological effect of vegetation in Karst watershed[J].Advances in Water Science，2007（5）：757-768.（in Chinese）） DOI：10.17521/cjpe.2007.0097.

[46] 張志才，陳喜，石朋，等.喀斯特流域分布式水文模型及植被生態(tài)水文效應(yīng)[J].水科學(xué)進(jìn)展，2009（6）：806-811.（ZHANG Z C，CHEN X，SHI P，et al.，Distributed hydrological model and eco-hydrological effect of vegetation in Karst watershed[J].Advances in Water Science，2009（6）：806-811.（in Chinese）） DOI：10.14042/j.cnki.32.1309.2009.06.003.

[47] TENNANT D L.Instream flow regimens for fish，wildlife，recreation and related environmental resources[J].Fisheries，1976，1（4）：6-10.DOI：10.1577/1548-8446（1976）001%3C0006：IFRFFW%3E2.0.CO;2.