皋 鵬 飛，李 玉 鳳，劉 紅 玉，季 香

(江蘇省地理環(huán)境演化國家重點實驗室培育建設(shè)點，江蘇 南京 210023；江蘇省地理信息資源開發(fā)與利用協(xié)同創(chuàng)新中心，江蘇 南京 210023；南京師范大學(xué)地理科學(xué)學(xué)院，江蘇 南京 210023)

水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力評估
——以風(fēng)嶺流域為例

皋 鵬 飛，李 玉 鳳*，劉 紅 玉，季 香

(江蘇省地理環(huán)境演化國家重點實驗室培育建設(shè)點，江蘇 南京 210023；江蘇省地理信息資源開發(fā)與利用協(xié)同創(chuàng)新中心，江蘇 南京 210023；南京師范大學(xué)地理科學(xué)學(xué)院，江蘇 南京 210023)

應(yīng)用InVEST模型和SCS模型并結(jié)合GIS技術(shù)，計算出儀征市風(fēng)嶺流域2015年的產(chǎn)水量、地表徑流量以及水塘系統(tǒng)的水源涵養(yǎng)量分別為92.07×106m3、20.11×106m3和36.47×106m3，在子流域上其總體變化均呈現(xiàn)為由流域下游向上游減少的趨勢。在此基礎(chǔ)上，研究評估了區(qū)域水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力，總體而言，流域下游的貢獻(xiàn)力均大于上游。在降水條件相同的情況下，建設(shè)用地分布較多的區(qū)域，產(chǎn)水量、地表徑流量以及水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力均大于林地覆蓋較多的區(qū)域；粘壤土分布區(qū)的地表徑流量大于壤土分布區(qū)；水塘系統(tǒng)水源涵養(yǎng)量隨地表徑流量的增加而增加。

InVEST模型；SCS模型；水塘系統(tǒng)；水源涵養(yǎng)；風(fēng)嶺流域

0 引言

水塘是指在平地鑿池、谷口或高地匯水處筑堤，就地蓄置雨水的小型水體，其面積在1～20 000 m2之間，具有自然、半自然和人工特征[1,2]。水塘系統(tǒng)由區(qū)域內(nèi)的水塘組成，是我國存在3 000多年的典型農(nóng)業(yè)小水利工程，其主要分布在南方丘陵地區(qū)。作為一種特殊的人工濕地，水塘系統(tǒng)具有調(diào)節(jié)氣候、水源涵養(yǎng)、水質(zhì)凈化、維持生物多樣性等重要的生態(tài)服務(wù)功能，在農(nóng)業(yè)發(fā)展過程中具有不可替代的作用和價值。其中，水塘系統(tǒng)的水源涵養(yǎng)功能在減緩水土流失、維持流域水量、彌補農(nóng)作物在需水時期的水資源短缺等方面起著非常重要的作用[3,4]，對于區(qū)域具有重要的水資源貢獻(xiàn)。因此定量評估區(qū)域內(nèi)水塘系統(tǒng)水源涵養(yǎng)的貢獻(xiàn)力就顯得尤為重要。

目前對于水源涵養(yǎng)功能的研究多集中在區(qū)域尺度上，如謝高地等對青藏高原水源涵養(yǎng)價值進(jìn)行評估[5]；陳龍等對瀾滄江流域的水源涵養(yǎng)功能進(jìn)行了研究[6,7]；傅斌等對都江堰市的水源涵養(yǎng)功能及其重要性進(jìn)行了研究[7]，Juszczak 等對波蘭西部的農(nóng)業(yè)森林流域水源涵養(yǎng)進(jìn)行了評估[8]，Pawattana 等對泰國池河流域的水源涵養(yǎng)能力進(jìn)行了研究[9]。但在這些研究中常常忽略水塘濕地的水源涵養(yǎng)功能及其對于整個區(qū)域的影響。因此圍繞水塘系統(tǒng)的水源涵養(yǎng)功能及其對于區(qū)域的貢獻(xiàn)程度進(jìn)行一個科學(xué)、定量的評估，是本文要解決的主要問題。

位于儀征市的風(fēng)嶺流域由于農(nóng)業(yè)灌溉的需要，區(qū)域內(nèi)分布了大量的水塘，形成了發(fā)達(dá)的水塘濕地系統(tǒng)，這為本文提供了理想的研究區(qū)域。本文利用相關(guān)模型和方法分析水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力及其空間分布特征，將有利于更好地認(rèn)識水塘濕地系統(tǒng)水資源利用過程及其生態(tài)效應(yīng)，更好地優(yōu)化區(qū)域內(nèi)水塘系統(tǒng)的布局，對于濕地研究、農(nóng)業(yè)開發(fā)和保護(hù)、區(qū)域水量平衡等具有重要意義。

1 研究區(qū)概況

高寶邵伯湖位于江蘇省中部，總面積876.76 km2，是江蘇省第三大湖，年平均水深1.46 m，是淮水入江的主要通道，且水質(zhì)良好(III-IV類水質(zhì))。研究區(qū)風(fēng)嶺流域位于江蘇省儀征市，在高寶邵伯湖的上游，是一個相對封閉的小流域，流域高程變化范圍為20～50 m，流域面積為17.31 km2。研究區(qū)土地利用方式主要包括林地、建設(shè)用地、水域濕地及耕地。氣候?qū)儆趤啛釒Ъ撅L(fēng)氣候，年均氣溫14 ℃，降雨豐沛，年降水量1 256.60 mm，降雨主要集中在6-9月，占全年降雨量的59.20%。土壤類型以黃棕壤和水稻土為主[10,11]。20世紀(jì)60-70年代，為了區(qū)域內(nèi)抗旱、防洪排澇，改良低洼易澇的鹽堿地，開始了大規(guī)模的水利設(shè)施建設(shè)。因此在風(fēng)嶺流域區(qū)域內(nèi)分布了大規(guī)模的水塘濕地系統(tǒng)，且至今持續(xù)使用。研究區(qū)具有水塘615個，平均面積0.28 hm2/個，總面積約175 hm2，占流域總面積10.11%。水塘系統(tǒng)為保證高寶邵伯湖的水質(zhì)發(fā)揮著重要作用，對于區(qū)域也具有較高的水源涵養(yǎng)貢獻(xiàn)力。

圖1 研究區(qū)地理位置

2 水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力評估方法

本文基于區(qū)域水量平衡法，對風(fēng)嶺流域內(nèi)水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力進(jìn)行評估。區(qū)域水量平衡法將研究對象看成一個“黑箱”，著眼于輸入水量與輸出水量，認(rèn)為降水量與蒸散量以及地表徑流量之差即為區(qū)域的水源涵養(yǎng)量[12]。本文應(yīng)用InVEST(the Integrate Valuation of Ecosystem Services and Tradeoffs Tool)模型計算研究區(qū)的產(chǎn)水量，即降水量減去蒸散量的部分；利用SCS模型計算研究區(qū)的地表徑流量。水塘系統(tǒng)的水源涵養(yǎng)量與區(qū)域水源涵養(yǎng)量的比值即為水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)率，其計算公式如下：

H=V/(W-E)

(1)

其中，H為水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)率(%)；V為水塘系統(tǒng)的水源涵養(yǎng)量(m3)，其利用GIS技術(shù)及相關(guān)數(shù)學(xué)方法等計算；W為區(qū)域產(chǎn)水量(m3)，包括地表徑流、土壤含水量、枯落物持水量和冠層截留量[13-15]；E為區(qū)域地表徑流量(m3)。

在計算出水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)率H的基礎(chǔ)上進(jìn)一步分級，分級標(biāo)準(zhǔn)如表1，從而評估出研究區(qū)水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力。

表1 水源涵養(yǎng)貢獻(xiàn)力分級

2.1 基于InVEST模型的產(chǎn)水量計算

InVEST模型是由美國斯坦福大學(xué)、世界自然基金會和大自然保護(hù)協(xié)會聯(lián)合開發(fā)的生態(tài)系統(tǒng)服務(wù)功能評估工具[16,17]，該模型包含水質(zhì)凈化、產(chǎn)水量、生物多樣性等多個模塊。本研究所用到的就是InVEST模型中的“產(chǎn)水量”模塊。

產(chǎn)水量模塊運行需要的參數(shù)如表2所示，其中，土壤有效含水量是指土壤中持有的并且可以被植物利用的那一部分水量[18]，取值范圍為0～1；參考作物蒸散量是指平坦地面被特定矮稈綠色植物(高0.12 m，地面發(fā)射率為0.23)全部遮蔽，同時土壤保持充分濕潤情況下的蒸散量[19]，它由Modified-Hargreaves法公式計算[20]：

ET0=0.0013×0.408×RA×(Tavg+17)×(TD-0.0123P)0.76

(2)

其中，ET0為參考作物蒸散量(mm/d)；RA為太陽大氣頂層輻射(MJ/(m2·d))，運用氣象站2015年太陽總輻射數(shù)據(jù)計算可得；Tavg是2015年日最高溫均值與日最低溫均值的均值(℃)，TD是日最高溫均值與日最低溫均值之差(℃)；P是降水量(mm)，降水量數(shù)據(jù)和氣溫數(shù)據(jù)收集于中國天氣網(wǎng)。

Zhang系數(shù)是表征降水特征的常數(shù)，取值范圍為1～10[21]。一般冬季降水集中的區(qū)域Zhang系數(shù)取10，終年濕潤或夏季降水集中的區(qū)域Zhang系數(shù)取1。研究區(qū)屬于夏季降水集中性的區(qū)域，因此Zhang系數(shù)取1[22,23]。生物物理系數(shù)表中需要確定的參數(shù)是植被蒸散系數(shù)K，它受到地表覆被的影響，因此，本文應(yīng)用FAO提出的適合于自然植被非完全覆蓋條件下的植被系數(shù)的計算方法估算[24]；子流域的范圍的劃分主要是根據(jù)風(fēng)嶺流域水文運移路徑(出水口位置、水流方向)確定子流域的邊界[25-27]。具體過程是在確保流域水塘濕地系統(tǒng)完整的情況下，以DEM為基礎(chǔ)，利用ArcGIS10.2的水文分析模塊，得到水流方向、匯流累積量等水文因子，并以此生成流域的河網(wǎng)。通過對大量實驗同現(xiàn)有地形進(jìn)行對比，研究得出當(dāng)河網(wǎng)的閾值為1 500時，河網(wǎng)的分布與實際最為接近，因此確定以1 500為最佳閾值劃分出35個子流域。

表2 產(chǎn)水量模型運行所需參數(shù)

2.2 基于SCS模型的地表徑流量計算

SCS(Soil Conservation Service)模型是美國農(nóng)業(yè)部水土保持局于1954年開發(fā)的流域水文模型，它能夠客觀反映土壤質(zhì)地、土地利用方式以及土壤含水量對降雨徑流的影響，是目前應(yīng)用最廣泛的流域水文模型之一[28]。模型計算公式：

E=R×K/1 000

(3)

R=(P-0.2L)2/(P+0.8L)，P≥0.2L

(4)

L=25 400/CN-254

(5)

其中，E是地表徑流量(m3)；R是地表產(chǎn)流深度(mm)；K是區(qū)域面積(m2)；P是降水量(mm);L是最大可能滯留深度(mm)，取決于土壤入滲速率和土壤儲水容量[29]；CN(Curve Number)是徑流曲線系數(shù)，是與土壤質(zhì)地和土地利用方式密切相關(guān)的一個徑流系數(shù)，CN值越大，越容易產(chǎn)生地表徑流[30]。

因此，模型所需要的參數(shù)分別是降水量P以及徑流曲線系數(shù)CN。徑流曲線系數(shù)CN值的大小主要受降水前濕潤程度、土壤質(zhì)地和土地利用方式3個因素的影響。降水前濕潤程度通過已經(jīng)獲得的降水量，經(jīng)過計算確定研究區(qū)域降水前土壤濕潤程度的等級為Ⅰ類(表3)；土壤質(zhì)地劃分通過1∶100萬土壤空間屬性數(shù)據(jù)查找獲得，結(jié)合研究區(qū)土壤質(zhì)地狀況，將研究區(qū)的土壤質(zhì)地劃分為粘壤土和壤土(圖2)；土地利用方式通過ArcGIS10.2對研究區(qū)的土地利用方式進(jìn)行解譯，并按照SCS模型特定的分類標(biāo)準(zhǔn)合并(圖3)。在此3個變量的基礎(chǔ)上，通過查找美國國家工程手冊第四章的CN值表得到研究區(qū)的CN值[31](表4)。在確定參數(shù)之后，按照式(3)-式(5)求出研究區(qū)的年地表徑流量。

表3 降水前濕潤程度等級劃分

圖2 研究區(qū)土壤質(zhì)地

圖3 研究區(qū)土地利用分類

表4 不同土地濕潤程度下各土地利用方式和土壤類型的CN值

Table 4 Runoff curve number of different type of land use and soil texture in different ARCs

土壤質(zhì)地濕潤程度森林農(nóng)田濕地/城鎮(zhèn)壤土粘壤土Ⅰ5175100Ⅱ7088100Ⅲ8595100Ⅰ5980100Ⅱ7791100Ⅲ8997100

2.3 水塘系統(tǒng)的水源涵養(yǎng)量計算

水塘作為特殊的人工濕地，在水源涵養(yǎng)方面具有重要作用[32]。水塘的年水源涵養(yǎng)量包括兩部分：一是水塘的基礎(chǔ)庫容量；二是水塘對村莊、農(nóng)田排水和降雨的截留增加量[33]，即：

V=V1+V2=(h1+Δh2)×S

(6)

其中，V為水塘的水源涵養(yǎng)量(m3)；V1為水塘的基礎(chǔ)庫容(m3)；V2為水塘一年內(nèi)對于水分的截留增加量(m3)；h1為水塘的基礎(chǔ)水深(m),Δh2為一年內(nèi)水塘的水位上升量(m)；S為水塘的表面積(m2)。

水塘的基礎(chǔ)水深測量方法是:在水塘兩岸選點設(shè)腳架，懸掛鋼尺進(jìn)入水中，通過鋼尺拉動測量基礎(chǔ)水深，并利用DEM進(jìn)行校驗獲得。DEM數(shù)據(jù)來源于江蘇省測繪局，其精度為5 m×5 m；水塘的水位上升量是指水塘一年內(nèi)水位上升的總和，其來源于2015年一年的野外實地調(diào)查，調(diào)查方法是基于60個水塘樣點，每個樣點確定一個基準(zhǔn)點，量取每周的水位變化量，統(tǒng)計獲得一年內(nèi)水塘水位的上升量，該數(shù)據(jù)能夠反映出水塘蓄水供水的情況。再根據(jù)該數(shù)據(jù)得到水塘的水源涵養(yǎng)總深度，并建立水塘系統(tǒng)的水源涵養(yǎng)深度與水塘表面積的線性模型：h1+Δh2=0.177×10-6×S+3.354，以此求出其他水塘的水源涵養(yǎng)量。

3 結(jié)果分析

3.1 子流域劃分

根據(jù)模型的要求，計算出產(chǎn)水量柵格結(jié)果后，要在子流域上進(jìn)行統(tǒng)計，因此，要對研究區(qū)風(fēng)嶺流域內(nèi)各個子流域的范圍進(jìn)行劃分。利用ArcGIS10.2水文分析模塊并根據(jù)實際地形，經(jīng)過大量的實驗對比劃分出35個子流域(圖4)，子流域的面積范圍在0.15～1.02 km2之間，坡度范圍在0°～35.12°之間。

圖4 研究區(qū)子流域分布

3.2 不同子流域的產(chǎn)水量

輸入模型所需要的參數(shù)后，運行InVEST模型，得到2015年研究區(qū)的產(chǎn)水量，研究區(qū)總體的產(chǎn)水量為92.07×106m3，產(chǎn)水量較大。為使產(chǎn)水量的結(jié)果進(jìn)行子流域之間的比較而不受子流域面積的影響，本文選用單位時間單位面積上的產(chǎn)水量(即產(chǎn)水模數(shù))來表達(dá)(圖5)。由圖5可知，產(chǎn)水量較大的子流域有13個，主要分布在流域的下游，包括子流域6～8、18～23、34～35以及1和4，產(chǎn)水量在605.41～1 066.64 m3/(km2·a)之間，子流域4的產(chǎn)水量最大。產(chǎn)水量一般的子流域有12個，主要分布在流域的下游，包括子流域2～3、11～15以及5、17、24、30、32，產(chǎn)水量在409.12～583.93 m3/(km2·a)之間。產(chǎn)水量較小的子流域有10個，主要分布在流域的上游，包括子流域9～10、25～29以及16、31、33，產(chǎn)水量在293.07～389.14 m3/(km2·a)之間，子流域28的產(chǎn)水量最小。

圖5 研究區(qū)各子流域產(chǎn)水量分布

因此可以得出結(jié)論，產(chǎn)水量在流域內(nèi)表現(xiàn)出明顯的空間差異，總體變化呈現(xiàn)由流域的下游向上游方向減少的趨勢。產(chǎn)水量的大小與降水量、土地利用方式、土壤有效含水量、土壤深度等有關(guān)，研究區(qū)的降水量總體一致，因此可以結(jié)合研究區(qū)土地利用方式等因素分析其原因：首先，流域的下游相對于上游分布了更加密集的建設(shè)用地，如城鎮(zhèn)、工廠及村莊等，建設(shè)用地多為不透水地面，降水難以下滲，相對于上游較容易產(chǎn)流；其次，流域下游的農(nóng)田作物根系深度較小，造成土壤有效含水量偏低，地表蒸發(fā)量較小，而上游的土地利用方式除了農(nóng)田以外，植被類型還包含了較多的林地，土壤的深度和植被的根系深度均較大，提高了土壤的有效含水量和地表蒸發(fā)量。因此總體上，流域下游的產(chǎn)水量比上游偏大。

3.3 不同子流域的地表徑流量

根據(jù)SCS模型計算得到研究區(qū)2015年地表徑流量為20.11×106m3，地表徑流量較大。為了使地表徑流量的結(jié)果進(jìn)行子流域之間的比較而不受子流域面積的影響，本文選用單位時間單位面積上的地表徑流量(即地表徑流模數(shù))來表達(dá)(圖6)。由圖6可知，地表徑流量較大的子流域有12個，其分布多集中在流域的下游，包括子流域14～15、18～23、34～35以及5、32，地表徑流量在97.70～99.60 m3/(km2·a)之間，子流域35的地表徑流量最大。地表徑流量一般的子流域有13個，主要分布在流域的下游，包括子流域3～4、6～8、11～13、29～30以及1、17、24，地表徑流量在96.19～97.69 m3/(km2·a)之間。地表徑流量較小的子流域有10個，集中分布在流域的上游，包括子流域9～10、25～28以及2、16、31、33，地表徑流量在92.01～95.18 m3/(km2·a)之間，子流域25的地表徑流量最小。

圖6 研究區(qū)各子流域地表徑流量分布

因此可以得出結(jié)論，地表徑流量在流域內(nèi)表現(xiàn)出明顯的空間差異，總體變化呈現(xiàn)由流域的下游向上游方向減少的趨勢。分析出現(xiàn)這種情況的原因：在降水相同的條件下，地表徑流量的大小與土壤質(zhì)地、土地利用方式密切相關(guān)。首先，流域下游的土壤質(zhì)地主要為粘壤土，質(zhì)地黏重，降水不易下滲，容易形成地表徑流而流失；而上游的土壤質(zhì)地為壤土，下滲能力較粘壤土要大，因此上游相對于下游形成的地表徑流量要小。其次，研究區(qū)的建設(shè)用地(如城鎮(zhèn)、工廠及村莊等)在流域下游的分布相對于上游更加密集，且下游分布了大量農(nóng)田，在這樣的土地利用方式下，降水更容易形成地表徑流；而流域上游的林地相對較多，對于降水的截留能力較強，地表徑流量較小。因此總體上，流域下游的地表徑流量比上游偏大。

3.4 不同子流域內(nèi)水塘系統(tǒng)的水源涵養(yǎng)量

通過GIS技術(shù)以及相關(guān)數(shù)學(xué)模型的計算，獲得2015年水塘系統(tǒng)的水源涵養(yǎng)量。研究區(qū)總體的水塘系統(tǒng)水源涵養(yǎng)量為36.47×106m3，水源涵養(yǎng)量較大。為了使水塘系統(tǒng)的水源涵養(yǎng)量的結(jié)果進(jìn)行子流域之間的比較而不受子流域面積的影響，本文選用單位時間單位面積上的水源涵養(yǎng)量(即水源涵養(yǎng)模數(shù))來表達(dá)(圖7)。由圖7可知，水塘系統(tǒng)的水源涵養(yǎng)量較大的子流域有12個，其分布集中在流域下游，包括子流域1～2、5～8、11～13、23～24以及18，水源涵養(yǎng)量在264.57～402.60 m3/(km2·a)之間，子流域1的地表徑流量最大。水源涵養(yǎng)量一般的子流域有13個，流域的上游和下游均有分布，包括子流域14～17、19～21、31～32以及3、9、27、29，水源涵養(yǎng)量在180.47～245.29 m3/(km2·a)之間。水源涵養(yǎng)量較小的子流域有10個，主要分布在流域的上游，包括子流域25～26、33～35以及4、10、22、28、30，水源涵養(yǎng)量在22.92～166.56 m3/(km2·a)之間，子流域22的地表徑流量最小。

因此可以得出結(jié)論，水塘的水源涵養(yǎng)量在流域內(nèi)表現(xiàn)出明顯的空間差異，總體變化呈現(xiàn)由流域的下游向上游方向減少的趨勢。分析出現(xiàn)這種情況的原因：水塘系統(tǒng)的水源涵養(yǎng)量主要來源于直接降水、農(nóng)田以及村鎮(zhèn)的排水徑流等，而流域的降水量是一致的，所以水塘系統(tǒng)的水源涵養(yǎng)量的空間差異主要受到地表徑流量的影響。經(jīng)過前面的分析，由于研究區(qū)的土壤質(zhì)地、土地利用方式等原因，導(dǎo)致流域產(chǎn)生的地表徑流量不同，且由流域的下游向上游減少。而流域內(nèi)的水塘系統(tǒng)能夠攔截地表徑流(如農(nóng)田排水)，并進(jìn)行存儲，使得徑流不至于全部流入河流，進(jìn)而實現(xiàn)水源涵養(yǎng)，地表徑流量越大，水塘系統(tǒng)的水源涵養(yǎng)量越大。因此總體上，流域下游水塘系統(tǒng)的水源涵養(yǎng)量大于上游。

圖7 各子流域水塘系統(tǒng)的水源涵養(yǎng)量分布

3.5 不同子流域內(nèi)水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力

基于研究區(qū)的年產(chǎn)水量、年地表徑流量以及水塘系統(tǒng)的年水源涵養(yǎng)量，得出水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)率(圖8a、表5)，并進(jìn)一步進(jìn)行分級，評估出水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力，評估結(jié)果如圖8b所示。

圖8 不同子流域水塘系統(tǒng)水源涵養(yǎng)貢獻(xiàn)率與貢獻(xiàn)力分布

由圖8b可知，水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力極大的子流域有9個，包括子流域1～2、12～13、27～29以及5、16，其分布在流域的下游，水源涵養(yǎng)貢獻(xiàn)率大于70%。水源涵養(yǎng)貢獻(xiàn)力較大的子流域有13個，在流域的上游和下游均有分布，包括子流域6～7、9～11、24～25、31～33以及3、14、17，水源涵養(yǎng)貢獻(xiàn)率在50%～70%之間。水源涵養(yǎng)貢獻(xiàn)力一般的子流域有8個，主要分布在流域的下游，包括子流域18～21以及8、15、23、26，水源涵養(yǎng)貢獻(xiàn)率在30%～50%之間。水源涵養(yǎng)貢獻(xiàn)力較小的子流域有5個，其主要分布在流域的上游，包括子流域34～35以及4、22、30，水源涵養(yǎng)貢獻(xiàn)率小于30%。

因此可以得出結(jié)論，水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力在流域內(nèi)表現(xiàn)出明顯的空間差異，總體來說，流域下游的水塘系統(tǒng)水源涵養(yǎng)貢獻(xiàn)力大于上游。分析出現(xiàn)這種情況的原因：在降水條件相同的情況下，水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力與土地利用方式密切相關(guān)。流域的下游分布著較密集的建設(shè)用地，其對降水量截留的較少，從而水塘系統(tǒng)截留的降水量較多，水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力較大[34];而流域的上游分布著較多的林地，其對降水量截留的較多，水塘系統(tǒng)截留的較少，因此流域上游的水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力一般甚至較小。

表5 研究區(qū)水塘系統(tǒng)水源涵養(yǎng)貢獻(xiàn)率百分比

4 結(jié)論與討論

本文以水量平衡法為基礎(chǔ)，應(yīng)用InVEST、SCS模型并結(jié)合GIS 技術(shù)，對研究區(qū)水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力進(jìn)行了科學(xué)、定量的評估分析，得出如下結(jié)論：1)風(fēng)嶺流域具有大規(guī)模的水塘系統(tǒng)，區(qū)域內(nèi)具有較大的年產(chǎn)水量(92.07×106m3)和年地表徑流量(20.11×106m3)，其總體變化呈現(xiàn)由流域的下游向上游方向減少的趨勢。在降水條件相同的情況下，建設(shè)用地覆蓋較多區(qū)域的產(chǎn)水量、地表徑流量大于林地覆蓋較多的區(qū)域；粘壤土分布區(qū)域的地表徑流量大于壤土分布的區(qū)域。2)區(qū)域內(nèi)水塘系統(tǒng)具有較大的水源涵養(yǎng)量(36.47×106m3)，其總體變化呈現(xiàn)由流域的下游向上游方向減少的趨勢。在降水條件相同的情況下，地表徑流量較大的區(qū)域，水塘系統(tǒng)的水源涵養(yǎng)量較大，反之水塘系統(tǒng)的水源涵養(yǎng)量較小。3)不同子流域內(nèi)水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力具有明顯的空間差異，總體而言，流域下游的水塘系統(tǒng)水源涵養(yǎng)貢獻(xiàn)力大于上游。在降水條件相同的情況下，建設(shè)用地分布較多的區(qū)域水塘系統(tǒng)水源涵養(yǎng)貢獻(xiàn)力大于林地覆蓋較多的區(qū)域。

本文對風(fēng)嶺流域內(nèi)水塘系統(tǒng)的水源涵養(yǎng)貢獻(xiàn)力進(jìn)行了評估，補充了以往在水塘濕地系統(tǒng)水源涵養(yǎng)功能評估方面的缺失。通過分析水源涵養(yǎng)貢獻(xiàn)力，能夠更好地優(yōu)化區(qū)域內(nèi)的水塘濕地系統(tǒng)，為水資源的存儲、水源供給利用以及農(nóng)業(yè)開發(fā)保護(hù)提供了理論基礎(chǔ)。結(jié)合InVEST模型和SCS模型能夠?qū)δ骋粎^(qū)域的產(chǎn)水量、地表徑流量及其差異進(jìn)行科學(xué)、定量的分析，兩種模型均以水量平衡法為基礎(chǔ)，具有較高的可行性，并且能夠反映出土地利用方式、土壤質(zhì)地等因素對于計算結(jié)果的影響。其中SCS模型雖然能夠計算出區(qū)域地表徑流量的差異，但是模型中的各變量關(guān)系是經(jīng)驗性的，各參數(shù)之間的關(guān)系以及CN值的確定需要進(jìn)一步討論。

[1] BIGGS J,WALKER D,WHITFIELD M.Pond action:Promoting the conservation of ponds in Britain[J].Freshwater Forum,1991,1(2):114-118.

[2] GAO J,WANG R,HUANG J.Ecological engineering for traditional Chinese agriculture—A case study of Beitang[J].Ecological Engineering,2015,76(1):7-13.

[3] 尹云鶴，吳紹洪，李華友，等.SRES情景下青藏高原生態(tài)功能保護(hù)區(qū)水源涵養(yǎng)功能的變化研究[J].資源科學(xué),2013,35(10):2003-2010.

[4] 劉賢德，李效雄，張學(xué)龍，等.干旱半干旱區(qū)山地森林類型的土壤水文特征[J].干旱區(qū)地理,2009,32(5):691-697.

[5] 謝高地，魯春霞，冷允法，等.青藏高原生態(tài)資產(chǎn)的價值評估[J].自然資源學(xué)報,2003,18(2):189-196.

[6] 陳龍，謝高地，張昌順，等.瀾滄江流域生態(tài)系統(tǒng)水源涵養(yǎng)功能研究(英文)[J].Journal of Resources and Ecology,2011,2(4):322-327.

[7] 傅斌，徐佩，王玉寬，等.都江堰市水源涵養(yǎng)功能空間格局[J].生態(tài)學(xué)報,2013,33(3):789-797.

[8] JUSZCZAK R,KE｀DZIORA A,OLEJNIK J.Assessment of water retention capacity of small ponds in Wyskoc＇ agricultural-forest catchment in western Poland[J].Polish Journal of Environmental Studies,2007,16(5):685-695.

[9] PAWATTANA C,TRIPATHI N K,HTWE S L.Development of potential floodwater retention zones using AHP and GIS:A case study in the Chi River Basin,Thailand[J].International Journal of Geoinformatics,2009,5(4):17-24.

[10] 李玉鳳，劉紅玉，朱麗娟，等.農(nóng)村多水塘系統(tǒng)中污染物來源評價——以陳集鎮(zhèn)為例[J].南京師大學(xué)報(自然科學(xué)版),2009,32(1):136-140.

[11] 李玉鳳，劉紅玉，郝敬鋒，等.農(nóng)村小流域景觀結(jié)構(gòu)與水質(zhì)耦合關(guān)系分析——以儀征市陳集鎮(zhèn)為例[J].生態(tài)與農(nóng)村環(huán)境學(xué)報,2010,26(1):15-19.

[12] 李盈盈.陜西省北洛河流域水源涵養(yǎng)生態(tài)服務(wù)功能及其價值估算[D].西安：西北大學(xué),2015.

[13] 馬良，金陶陶，文一惠，等.InVEST模型研究進(jìn)展[J].生態(tài)經(jīng)濟,2015,31(10):126-131.

[14] 唐堯，祝煒平，張慧，等.InVEST模型原理及其應(yīng)用研究進(jìn)展[J].生態(tài)科學(xué),2015,34(3):204-208.

[15] 潘韜，吳紹洪，戴爾阜，等.基于InVEST模型的三江源區(qū)生態(tài)系統(tǒng)水源供給服務(wù)時空變化[J].應(yīng)用生態(tài)學(xué)報,2013,24(1):183-189.

[16] 余新曉，周彬，呂錫芝，等.基于InVEST模型的北京山區(qū)森林水源涵養(yǎng)功能評估[J].林業(yè)科學(xué),2012,48(10):1-5.

[17] 陳駿宇，劉鋼，白楊.基于InVEST模型的太湖流域水源涵養(yǎng)服務(wù)價值評估[J].水利經(jīng)濟,2016,34(2):25-29.

[18] 李盈盈，劉康，胡勝，等.陜西省子午嶺生態(tài)功能區(qū)水源涵養(yǎng)能力研究[J].干旱區(qū)地理,2015,38(3):636-642.

[19] 劉鈺，蔡林根.參照騰發(fā)量的新定義及計算方法對比[J].水利學(xué)報,1997,4(6):28-34.

[20] 賈芳芳.基于InVEST模型的贛江流域生態(tài)系統(tǒng)服務(wù)功能評估[D].北京：中國地質(zhì)大學(xué)(北京),2014.

[21] ZHANG L,DAWES W R,WALKER G R.Response of mean annual evapotranspiration to vegetation changes at catchment scale[J].Water Resources Research,2001,37(3):701-708.

[22] ZHANG C Q,LI W H,ZHANG B,et al.Water yield of Xitiaoxi River basin based on InVEST modeling[J].Journal of Resources and Ecology,2012,3(1):50-54.

[23] 張燦強，李文華，張彪，等.基于InVEST模型的西苕溪流域產(chǎn)水量分析(英文)[J].Journal of Resources and Ecology,2012,3(1):50-54.

[24] ZHOU W,LIU G,PAN J,et al.Distribution of available soil water capacity in China[J].Journal of Geographical Sciences,2005,15(1):3-12.

[25] 程崢，李永勝，高微微.基于ArcGIS的DEM流域劃分[J].地下水,2011,33(6):128-130.

[26] 王俊啟.基于1:1萬地形圖的數(shù)字流域劃分研究[D].杭州：浙江大學(xué),2013.

[27] 宋曉猛，張建云，占車生，等.基于DEM的數(shù)字流域特征提取研究進(jìn)展[J].地理科學(xué)進(jìn)展， 2013,32(1):31-40.

[28] 劉家福，蔣衛(wèi)國，占文鳳，等.SCS模型及其研究進(jìn)展[J].水土保持研究,2010,17(2):120-124.

[29] 鄭暢，倪九派，魏朝富.基于DEM和SCS模型的四川盆地丘陵區(qū)局地徑流研究[J].水土保持學(xué)報,2008,22(5):73-77.

[30] 符素華，王向亮，王紅葉，等.SCS-CN徑流模型中CN值確定方法研究[J].干旱區(qū)地理,2012,35(3):415-421.

[31] 董文濤，程先富，張群，等.SCS-CN模型在巢湖流域地表產(chǎn)流估算中的應(yīng)用[J].水土保持通報,2012,32(3):174-177.

[32] 熊遠(yuǎn)清，吳鵬飛，張洪芝，等.若爾蓋濕地退化過程中土壤水源涵養(yǎng)功能[J].生態(tài)學(xué)報,2011,31(19):5780-5788.

[33] 毛戰(zhàn)坡，尹澄清，單寶慶，等.水塘系統(tǒng)對農(nóng)業(yè)流域水資源調(diào)控的定量化研究[J].水利學(xué)報,2003,1(12):76-83.

[34] WANG C,ZHAO C Y,XU Z L,et al.Effect of vegetation on soil water retention and storage in a semi-arid alpine forest catchment[J].Journal of Arid Land,2013,5(2):207-219.

Evaluation on the Contribution Capacity of Water Conservation of Pond System:A Case Study of Fengling Watershed

GAO Peng-fei,LI Yu-feng,LIU Hong-yu,JI Xiang

(StateKeyLaboratoryCultivationBaseofGeographicalEnvironmentEvolutionofJiangsuProvince,Nanjing210023;JiangsuCenterforCollaborativeInnovationinGeographicalInformationResourceDevelopmentandApplication,Nanjing210023;CollegeofGeographicalScience,NanjingNormalUniversity,Nanjing210023,China)

In this paper,taking Fengling watershed in Yizheng as an example,water yield and surface runoff in 2015 were obtained by applying the InVEST model and SCS model.Based on the data above,water conservation of pond system was calculated.The results were as follows:water yield,surface runoff and water conservation of pond system were 92.07×106m3,20.11×106m3and 36.47×106m3respectively.Generally speaking,the value of water yield,surface runoff and water conservation tended to decrease from downstream to upstream.Based on these,the contribution capacity of water conservation of pond system was evaluated.On the whole,the contribution capacity of water conservation in the downstream was greater than that in the upstream.Under the same condition of rainfall,water yield,surface runoff and the contribution capacity of water conservation of the areas covered by more construction was larger than that of the areas covered by more forests;surface runoff of the clay loam area was larger than that of the loam area;water conservation of pond system increased with surface runoff.

InVEST model;SCS model;pond system;water conservation;Fengling watershed

2016-06-18；

2016-08-18

國家自然科學(xué)基金項目(41401205、31570459)；江蘇省自然科學(xué)基金項目(BK20140921)；江蘇省高校自然科學(xué)研究重點項目(15KJA170002)；江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目(164320H116)；江蘇省地理信息資源開發(fā)與利用協(xié)同創(chuàng)新中心資助

皋鵬飛(1992-)，男，碩士研究生，主要從事濕地景觀與GIS應(yīng)用研究。*通訊作者E-mail：pandalee_0826@163.com

10.3969/j.issn.1672-0504.2016.06.016

P333.1

A

1672-0504(2016)06-0094-07