馬偉龍++任平++陶曉明

摘要：為更好地發(fā)揮耕地生態(tài)系統(tǒng)調(diào)節(jié)水源的作用，科學(xué)量化耕地涵養(yǎng)水源價值就顯得尤為重要。在收集整理相關(guān)數(shù)據(jù)的基礎(chǔ)上，借助GIS技術(shù)分析了四川省耕地涵養(yǎng)水源總量及空間特征，并且運用影子工程法、替代市場法估算了耕地涵養(yǎng)水源功能價值。結(jié)果表明，2011年四川省耕地生態(tài)系統(tǒng)涵養(yǎng)水源總量為309.70×108 m3，涵養(yǎng)水源的總價值為596.17億元，其中水田涵養(yǎng)水源總量為159.76×108 m3，價值為307.54億元；旱地涵養(yǎng)水源總量為149.94×108 m3，價值為288.63億元；農(nóng)作物的枝葉層、耕地土壤層是耕地涵養(yǎng)水源的主要作用層，其涵養(yǎng)水源的總量分別為1.48×108、308.22×108 m3，價值分別為2.84億、593.33億元。四川省耕地生態(tài)系統(tǒng)涵養(yǎng)水源價值較大，為當(dāng)年農(nóng)業(yè)總產(chǎn)值的24.29%；水田、旱地涵養(yǎng)水源能力分別占總價值的51.59%、48.41%，土層蓄水能力遠(yuǎn)高于農(nóng)作物枝葉涵養(yǎng)水源能力；耕地涵養(yǎng)水源空間分布特征為川東部盆地、平原地區(qū)涵養(yǎng)水源能力高于川西南山地、川西北高原。四川省耕地涵養(yǎng)水源功能及生態(tài)經(jīng)濟效益不容忽視，研究結(jié)果對四川省耕地資源保護有著重要的理論參考價值。

關(guān)鍵詞：涵養(yǎng)水源量；涵養(yǎng)水源價值；耕地生態(tài)系統(tǒng)；空間分析；四川省

中圖分類號： S181文獻標(biāo)志碼： A文章編號：1002-1302（2015）09-0399-04

耕地生態(tài)系統(tǒng)是受人類干預(yù)最為明顯的經(jīng)過人類改造的自然-人工復(fù)合生態(tài)系統(tǒng)[1]，也是地球上最重要的陸地生態(tài)系統(tǒng)之一，在為人類提供糧食的同時，也為人類提供涵養(yǎng)水源等生態(tài)服務(wù)功能[2]。涵養(yǎng)水源功能是耕地生態(tài)系統(tǒng)生態(tài)效益的重要組成部分[3]，被經(jīng)濟學(xué)家稱為耕地環(huán)境資源的非消耗性利用價值[4]。城市化進程加快導(dǎo)致城市人口快速增加、城鎮(zhèn)建設(shè)用地面積迅速擴張[5]、經(jīng)濟快速發(fā)展，對土地資源和水資源的需求不斷上升，進而出現(xiàn)土地資源尤其是耕地資源不斷減少[6]、水資源短缺[7-8]。不合理的耕地利用方式也導(dǎo)致耕地生態(tài)系統(tǒng)以前所未有的速度退化，耕地生態(tài)價值損失空前。因此，開展耕地生態(tài)系統(tǒng)涵養(yǎng)水源服務(wù)價值研究不僅可以認(rèn)識耕地生態(tài)系統(tǒng)的重要性[3]，豐富生態(tài)系統(tǒng)服務(wù)功能理論體系，而且可以在一定程度上對防治水土流失、保護耕地具有一定的現(xiàn)實意義。

生態(tài)系統(tǒng)涵養(yǎng)水源服務(wù)價值研究始于1997年Costanza等的報道[9]，目前生態(tài)系統(tǒng)涵養(yǎng)水源服務(wù)價值的研究對象多為森林生態(tài)系統(tǒng)[10]、濕地生態(tài)系統(tǒng)、草地生態(tài)系統(tǒng)[11]，以耕地生態(tài)系統(tǒng)為對象研究其涵養(yǎng)水源價值的案例較少?，F(xiàn)有研究主要采用區(qū)域水量平衡法[12-13]、土壤蓄水估算法[14]、多因子回歸法[15]、地下徑流增長法[16]計算耕地涵養(yǎng)水源物質(zhì)量，然后運用影子工程法[17-18]、替代市場法[1]等方法來評價生態(tài)系統(tǒng)涵養(yǎng)水源的效益。但耕地生態(tài)系統(tǒng)與人類的關(guān)系最為密切且受人類活動的影響較為深刻，至今尚未形成準(zhǔn)確、統(tǒng)一的評價模型[3]。本研究在總結(jié)前人有關(guān)生態(tài)系統(tǒng)涵養(yǎng)水源研究成果基礎(chǔ)上，進一步考慮降水強度因素、農(nóng)閑田土季節(jié)性差異，即冬季休耕、春夏秋播種對耕地生態(tài)系統(tǒng)涵養(yǎng)水源量的影響，同時以研究區(qū)內(nèi)每個研究單元平均蓄水庫容成本的加權(quán)平均值修正單位蓄水成本；采用修正后的土壤蓄水估算法測算耕地涵養(yǎng)水源物質(zhì)量[16，19]，用影子工程法、替代市場法測算耕地涵養(yǎng)水源價值；并在地理信息系統(tǒng)（GIS）技術(shù)支持下，對耕地生態(tài)系統(tǒng)涵養(yǎng)水源價值空間特征進行分析，以期為四川省耕地資源保護提供參考。

1研究區(qū)概況與數(shù)據(jù)來源

四川省位于我國西南腹地，地理位置為97°21′～108°31′E，26°03′～34°19′N，東西長1 075 km，南北寬 921 km，土地面積4 850萬 hm2。2011年末實有耕地面積398.34萬 hm2，其中水田面積206.93萬 hm2，旱地面積191.41萬 hm2，分別占耕地面積的51.95%、48.05%，年內(nèi)共減少2.73萬 hm2。該區(qū)域地形復(fù)雜，氣候類型多樣，以亞熱帶季風(fēng)氣候為主，降水集中在夏季，年均降水量900～1 200 mm。

數(shù)據(jù)主要來源于《四川省統(tǒng)計年鑒：2012》[20]、《四川省土地利用總規(guī)劃》，結(jié)合四川省地形圖、行政區(qū)劃圖、耕地類型圖、氣候圖、四川省多個氣象站觀測資料和氣候公報等資料，在GIS軟件平臺下，對所得數(shù)據(jù)進行數(shù)字化處理，建立空間數(shù)據(jù)庫、屬性數(shù)據(jù)庫。

2研究方法

生態(tài)系統(tǒng)涵養(yǎng)水源功能主要涉及林冠層、枯枝落葉層、土壤層，由于耕地生態(tài)系統(tǒng)中基本沒有枯枝落葉，所以本研究重點考慮農(nóng)作物枝葉層、土壤層的水文特征。耕地主要由水田、旱地組成，水田以種植水稻為主，旱地以種植玉米為主，因此，計算時水田以水稻為標(biāo)準(zhǔn)，旱地以玉米為標(biāo)準(zhǔn)，其他農(nóng)作物換算成這2種農(nóng)作物來測算。計算枝葉截留降水量時以農(nóng)作物播種面積為準(zhǔn)，計算土壤持水能力時以旱地水田實有面積為準(zhǔn)，使測算結(jié)果更加符合實際。

2.1耕地涵養(yǎng)水源量的估算方法

2.1.1農(nóng)作物截留降水量測算農(nóng)作物枝葉截留降水量與降水量強度、植被覆蓋度、枝葉葉面積指數(shù)密切相關(guān)。已有眾多學(xué)者提出了植被截留降水量估算方法[21-22]，但何東進等認(rèn)為，這些方法存在唯象因素使得結(jié)果存在不確定性，并提出不含唯象性質(zhì)的植被截留降水量測算方法[23-24]。不唯象研究方法也存在1個問題，即對降水量因素考慮不夠深入。本研究以最大截留降水量為臨界點分段計算截留降水量，替代平均降水量測算截留降水量，修正后的測算模型如下：

式中：Q1為單位面積農(nóng)作物截留降水量，m3；VEG為單位面積農(nóng)作物覆蓋度，%；LAI為單位面積農(nóng)作物葉面積指數(shù)；J為1個降水過程水量，mm；α為農(nóng)作物類型；Q1為單位面積農(nóng)作物最大截留降水量，mm；A為葉面上平均最大持水深度，mm，變動范圍約為0.1～0.3 mm[25]。endprint

2.1.2土壤蓄水量測算土壤降水儲存量與土壤類型及面積、土壤地下水深度、土壤粗孔隙率成正比[26]，按土壤類型測算后再求和。測算土壤蓄水量時以土壤的飽和含水量為臨界點分段計算土壤蓄水量，彌補過去研究中不考慮降水強度因素的不足，修正后模型如下：

式中：Q2為單位面積土壤降水儲存量，m3；J為降水過程的單位面積水雨量，mm；Hi為第i類耕地土層深度，m；ρi為第i類耕地的土壤容重，kg/m3；λi為第i類耕地的田間持水量，m3；Q2*為單位面積土壤最大蓄水量，m3。

2.2耕地涵養(yǎng)水源價值評估

2.2.1耕地涵養(yǎng)水源總量測算耕地生態(tài)系統(tǒng)涵養(yǎng)水源總量關(guān)系式為：

QW=∑i（QSi×mi） 。 （5）

式中：QW為四川省耕地生態(tài)系統(tǒng)涵養(yǎng)水源總量，m3；QSi為第i類耕地單位面積土壤涵養(yǎng)水源量，m3/hm2；mi為第i類耕地面積，hm2；i為耕地類型。

2.2.2耕地水源涵養(yǎng)量經(jīng)濟價值測算本研究在考慮人們對生態(tài)價值的支付意愿、替代工程易于計價的基礎(chǔ)上，運用替代市場法計算耕地生態(tài)系統(tǒng)涵養(yǎng)水源價值總量，測算時用經(jīng)過生態(tài)價值發(fā)展階段系數(shù)修正后的每個研究單元平均蓄水庫容成本的加權(quán)平均值，替代研究區(qū)單位蓄水量庫容成本，公式如下：

V= QW×P 。（6）

式中：V為耕地涵養(yǎng)水源價值，元；QW為四川省耕地生態(tài)系統(tǒng)涵養(yǎng)水源總量，m3；P為研究單元平均蓄水庫容成本的加權(quán)平均值修正單位蓄水成本，為1.925元/m3。

3結(jié)果與分析

3.1四川省耕地生態(tài)系統(tǒng)水源涵養(yǎng)量及其價值分析

運用式（1）、式（2）、式（3）、式（4）分別計算Q1、Q2值，根據(jù)式（5）、式（6）以及水田、旱地的面積計算四川省耕地水源涵養(yǎng)能力（表1）。由表1可知，四川省2011年耕地生態(tài)系統(tǒng)水源涵養(yǎng)總量為309.70×108 m3，總價值為596.17億元。其中農(nóng)作物截留水量1.48×108 m3，價值為2.84億元；土層蓄水量308.22×108 m3，價值為593.33億元。表12011年四川省耕地水源涵養(yǎng)總量及價值

耕地類別2011年耕地面積

（萬hm2）農(nóng)作物截留降水量土層蓄水量合計總量（×108 m3）價值（億元）總量（×108 m3）價值（億元）總量（×108 m3）價值（億元）水田206.930.841.62158.92305.92159.76307.54旱地191.410.641.22149.30287.41149.94288.63合計398.341.482.84308.22593.33309.70596.17

由表1還可以看出，四川省耕地分為水田、旱地，其面積分別占總耕地面積的51.95%、48.05%，水田、旱地涵養(yǎng)水源總量分別為159.76×108、149.94×108 m3；價值分別為30754億、288.63億元，分別占總價值的51.59%、48.41%。水田涵養(yǎng)水源價值大于旱地涵養(yǎng)水源的價值（表1、圖1）。2種耕地生態(tài)系統(tǒng)均以土層涵養(yǎng)水源價值為主，農(nóng)作物截留降水量價值相對土層蓄水價值較小，但也不能忽視農(nóng)作物截留降水總價值。

隨著城市化進程加快，城市建設(shè)占用大量耕地以及人類不合理利用耕地導(dǎo)致四川省耕地面積逐漸減少。截至2011年年底，四川省年末實有耕地面積398.34萬hm2，年內(nèi)減少了27 300 hm2（其中，水田減少26 500 hm2，旱地減少800 hm2），導(dǎo)致耕地涵養(yǎng)水源總量減少了2.11×108 m3，總價值減少了4.06億元。從圖1 還可以看出，在氣候、地形等因素的影響下，四川省的攀枝花市、綿陽市、廣元市、遂寧市、南充市、雅安市、資陽市和3大自治州耕地涵養(yǎng)水源價值以旱地為主，均達(dá)到50%以上，其中3大自治州旱地涵養(yǎng)水源能力更為突出，占其耕地涵養(yǎng)水源價值的80%以上；其他市（州）則以水田涵養(yǎng)水源能力為主。

3.2四川省耕地生態(tài)系統(tǒng)涵養(yǎng)水源功能空間特征分析

根據(jù)表2計算結(jié)果，在Arcgis10.1軟件中建立空間數(shù)據(jù)庫，并用自然斷點分級法對耕地涵養(yǎng)水源價值進行分級，得到四川省耕地涵養(yǎng)水源價值空間分布情況（圖2），將單位面積涵養(yǎng)水源價值與行政區(qū)劃圖疊加得到耕地單位面積涵養(yǎng)水源價值分布情況（圖3）。

在耕地涵養(yǎng)水源價值中，涼山州涵養(yǎng)水源價值最高，為52.97億元，占四川省耕地涵養(yǎng)水源總價值的8.89%；其次為成都市，價值為48.56億元，占8.15%；攀枝花市涵養(yǎng)水源價值最低，為6.10億元，占1.02%；雅安市涵養(yǎng)水源價值比攀枝花市多2.26億元，名列倒數(shù)第2（表2）。耕地生態(tài)系統(tǒng)涵養(yǎng)水源能力受多方面因素影響，不同研究單元耕地涵養(yǎng)水源能力差異主要體現(xiàn)在耕地面積、降水量等方面。涼山州耕地涵養(yǎng)水源價值最大，究其原因是其耕地面積最大，且降水豐富；成都市價值僅次于涼山州是因為成都市地處成都平原，耕地面積廣，尤其是水田面積占比大；攀枝花市價值最低有2個方面原因，即耕地面積有限、年降水量相對較少；雅安市雖然降水量豐富，但退耕還林導(dǎo)致耕地面積大為減少，所以耕地涵養(yǎng)水源價值也不高。

從單位面積涵養(yǎng)水源價值看，四川省耕地涵養(yǎng)水源生態(tài)服務(wù)功能單位面積價值為14 966元/hm2，受耕地面積數(shù)量的影響，單位面積價值高于15 030元/hm2的有阿壩州、甘孜州；單位面積價值低于14 931元/hm2有成都市、瀘州市、廣安市、眉山市；其他介于14 931～15 030元/hm2（圖3）。

從整體上看，四川省耕地涵養(yǎng)水源功能價值空間分布特征為：平原、盆地>川西南山地>川西北高原。在盆地內(nèi)部尤其是成都平原因地勢平坦且耕地面積廣闊，利于耕地涵養(yǎng)水源，耕地涵養(yǎng)水源總價值相對較高；而川西南山地和川西北高原地區(qū)受地形、氣候等因素影響，耕地涵養(yǎng)水源能力不如盆地底部及邊緣地區(qū)。endprint

4結(jié)論與討論

4.1結(jié)論

從耕地生態(tài)系統(tǒng)作為人類干預(yù)最為深刻的自然-人工復(fù)合生態(tài)系統(tǒng)角度出發(fā)，綜合考慮不同農(nóng)作物、不同耕地類型、耕地復(fù)種指數(shù)、天氣狀況、人類經(jīng)濟活動等因素，使修正后的指標(biāo)更符合耕地涵養(yǎng)水源價值指標(biāo)體系，最大限度反映耕地特有屬性，具有一定的可行性和科學(xué)性。

總體來看，2011年四川省耕地涵養(yǎng)水源總量為30970×108 m3，總價值為596.17億元，占當(dāng)年農(nóng)業(yè)總產(chǎn)值 2 454.26 億元的24.29%，可見耕地生態(tài)系統(tǒng)涵養(yǎng)水源價值巨大，而當(dāng)年由于耕地面積減少導(dǎo)致的耕地涵養(yǎng)水源價值損失為2.59億元，損失也相對較大。

從耕地類型來看，四川省水田涵養(yǎng)水源價值為30754億元，旱地涵養(yǎng)水源價值為288.63億元，分別占耕地涵養(yǎng)水源總價值的51.59%、48.41%，說明水田和旱地涵養(yǎng)水源能力相當(dāng)。另外，土層蓄水量價值為593.33億元，占涵養(yǎng)水源總價值的99.52%，其涵養(yǎng)水源能力遠(yuǎn)高于農(nóng)作物枝葉涵養(yǎng)水源能力。

從空間上來看，四川省耕地涵養(yǎng)水源價值空間分布特征表2四川省耕地水源涵養(yǎng)價值（V）

地區(qū)面積

表現(xiàn)為川東部盆地、平原地區(qū)高于川西南山地和川西北高原，涼山州、成都市、達(dá)州市、南充市屬于總價值最高的區(qū)域，而甘孜州、阿壩州、雅安市、攀枝花市屬于總價值最低的區(qū)域，說明耕地涵養(yǎng)水源價值深受耕地面積、地形、氣候等因素的影響。

4.2討論

耕地生態(tài)系統(tǒng)受較多外來不確定因素的影響，目前尚未形成公認(rèn)統(tǒng)一的測算耕地涵養(yǎng)水源模型。本研究在參考前人有關(guān)涵養(yǎng)水源測算模型的基礎(chǔ)上，對測算模型的相關(guān)參數(shù)進行改進，即以最大截留降水量為臨界點分段計算截留降水量，替代平均降水量測算截留降水量；以土壤的飽和含水量為臨界點分段計算土壤蓄水量，彌補了已有研究中不考慮降水強度因素的不足；用經(jīng)過生態(tài)價值發(fā)展階段系數(shù)修正后的每個研究單元平均蓄水庫容成本的加權(quán)平均值，替代單位蓄水量庫容成本，形成新的耕地涵養(yǎng)水源測評模型，更能準(zhǔn)確反映四川省耕地生態(tài)系統(tǒng)涵養(yǎng)水源價值及空間分布特征，為以后相關(guān)研究提供了一種新思路。值得說明的是，耕地實際水源涵養(yǎng)價值如同價值規(guī)律一般圍繞理論水源涵養(yǎng)價值上下波動，在后續(xù)的研究中應(yīng)圍繞理論值與實際值的差距或者實際值圍繞理論值波動的規(guī)律及其深層原因進行討論，將評價單元放在具體的環(huán)境中研究，考慮功能的區(qū)域差異，引入耕地水循環(huán)模型，同時考慮各分量間以及耕地生態(tài)系統(tǒng)與人類活動的相互關(guān)系，才能更準(zhǔn)確地度量耕地生態(tài)系統(tǒng)涵養(yǎng)水源的價值，這是今后進一步研究的方向。

參考文獻：

[1]王勇，駱世明，黃門福，等. 農(nóng)林生態(tài)系統(tǒng)截流蓄水的功能及其核算方法[J]. 廣東農(nóng)業(yè)科學(xué)，2007，34（5）：51-54.

[2]張穎聰. 四川省耕地生態(tài)服務(wù)價值時空變化研究[D]. 雅安：四川農(nóng)業(yè)大學(xué)，2012.

[3]李晶，張微微. 關(guān)中-天水經(jīng)濟區(qū)農(nóng)田生態(tài)系統(tǒng)涵養(yǎng)水源價值量時空變化[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報，2014，35（3）：52-57.

[4]李晶，任志遠(yuǎn). 陜北黃土高原生態(tài)系統(tǒng)涵養(yǎng)水源價值的時空變化[J]. 生態(tài)學(xué)雜志，2008，27（2）：240-244.

[5]曾杰，李江風(fēng)，姚小薇. 武漢城市圈生態(tài)系統(tǒng)服務(wù)價值時空變化特征[J]. 應(yīng)用生態(tài)學(xué)報，2014，25（3）：883-891.

[6]封志明，劉寶勤，楊艷昭. 中國耕地資源數(shù)量變化的趨勢分析與數(shù)據(jù)重建：1949～2003[J]. 自然資源學(xué)報，2005，20（1）：35-43.

[7]鄧偉，何巖. 水資源：21世紀(jì)全球更加關(guān)注的重大資源問題之一[J]. 地理科學(xué)，1999，19（2）：97-101.

[8]馮起，程國棟，三上正男，等. 中國水資源的問題和對策[J]. AMBIO-人類環(huán)境雜志，1999，28（2）：202-203.

[9]Costanza R，dArge R，de Groot R，et al.The value of the worlds ecosystem services and natural capital[J]. Nature，1997，387：253-260.

[10]張彪，李文華，謝高地，等. 森林生態(tài)系統(tǒng)的水源涵養(yǎng)功能及其計量方法[J]. 生態(tài)學(xué)雜志，2009，28（3）：529-534.

[11]于格，魯春霞，謝高地，等. 基于RS和GIS的青藏高原草地生態(tài)系統(tǒng)土壤水分保持功能及其經(jīng)濟價值評估——以生長季為例[J]. 山地學(xué)報，2006，24（4）：498-503.

[12]肖玉，謝高地，安凱，等. 華北平原小麥-玉米農(nóng)田生態(tài)系統(tǒng)服務(wù)評價[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報，2011，19（2）：429-435.

[13]鄒君，劉蘭芳，謝小立. 湖南亞熱帶丘崗區(qū)稻田生態(tài)系統(tǒng)水分平衡研究[J]. 熱帶地理，2002，22（3）：270-274.

[14]白楊，歐陽志云，鄭華，等. 海河流域農(nóng)田生態(tài)系統(tǒng)環(huán)境損益分析[J]. 應(yīng)用生態(tài)學(xué)報，2010，21（11）：2938-2945.

[15]王禮先，王斌瑞，朱金兆，等. 森林生態(tài)工程學(xué)[M]. 北京：中國林業(yè)出版社出版，2000.

[16]張燦強，李文華，張彪，等. 基于土壤動態(tài)蓄水的森林水源涵養(yǎng)能力計量及其空間差異[J]. 自然資源學(xué)報，2012，27（4）：697-704.

[17]莊樹宏，初洋，卞福花，等. 廟島群島海岸帶農(nóng)田生態(tài)系統(tǒng)服務(wù)價值的評估[J]. 煙臺大學(xué)學(xué)報：自然科學(xué)與工程版，2008，21（4）：273-280.

[18]王冬銀，楊慶媛，何濤. 重慶市耕地資源非市場價值估算[J]. 中國土地科學(xué)，2013，27（10）：76-82.

[19]宋吉紅，王百田，林富榮. 黃土高原旱地果園土壤蓄水保墑技術(shù)定量研究[J]. 水土保持學(xué)報，2000，14（4）：95-98.

[20]四川省統(tǒng)計局.四川統(tǒng)計年鑒：2012[M]. 北京：中國統(tǒng)計出版社，2012.

[21]Street-Perrott A，Beran M，Ratcliffe R. Variations in the global water budget[M]. Boston：Reidel Publishing Company，1983：147-171.

[22]Mintz Y，Walker G K. Global fields of soil moisture and land surface evapotranspiration derived from observed precipitation and surface air temperature[J]. J Appl Meteor，1993，32（8）：1305-1334.

[23]何東進，洪偉. 植被截留降水量公式的改進[J]. 農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，1999，15（3）：200-202.

[24]儀垂祥，劉開瑜，周濤. 植被截留降水量公式的建立[J]. 土壤侵蝕與水土保持學(xué)報，1996，2（2）：47-49.

[25]Manzi A O，Planto S.Implentation of the ISBA parameterization scheme for land surface processes in a GCM-an annual cycle experiment[J]. Journal of Hydrology，1994，155（3/4）：353-387.

[26]任志遠(yuǎn)，李晶. 陜南秦巴山區(qū)植被生態(tài)功能的價值測評[J]. 地理學(xué)報，2003，58（4）：503-511.endprint