張海波，張明陽，王克林，秦建新 ，符 靜

隨著水資源需求量的不斷增加以及水環(huán)境的急劇惡化，生態(tài)系統(tǒng)水源涵養(yǎng)作用已引起高度重視[1]。雖然一些人為工程在一定程度上能緩解水資源短缺的危機，但隨之而來的卻是各種環(huán)境問題。自然生態(tài)系統(tǒng)中森林、灌叢和草地等有著巨大的水源涵養(yǎng)功能，具有“天然綠色水庫”的美譽，并且能通過植被冠層、枯枝落葉層和土壤層的綜合作用來緩和地表徑流、補充地下水和調(diào)節(jié)和流量[2]，而濕地更是區(qū)域生態(tài)環(huán)境的穩(wěn)定器和調(diào)節(jié)器，對地區(qū)的經(jīng)濟和社會發(fā)展有著重要意義[3]。因此，生態(tài)系統(tǒng)水源涵養(yǎng)功能的分析對了解地區(qū)水資源利用效益及生態(tài)環(huán)境具有重大的意義。目前評價水源涵養(yǎng)功能的方法主要有區(qū)域水量平衡法[4]、土壤蓄水法[5]、地下徑流增長法[6]和降雨儲存量法[7]，同時一些學(xué)者采用InVEST模型進行定量評估[8-10]，也有一些學(xué)者利用轉(zhuǎn)移矩陣，結(jié)合景觀格局變化來分析生態(tài)系統(tǒng)服務(wù)功能的變化特征[11]。然而，當(dāng)前相關(guān)的研究多注重單一時相定量評估，對多時相動態(tài)變化分析略顯不足。

南方丘陵山地帶作為我國主體生態(tài)功能區(qū)劃中“兩屏三帶”國家生態(tài)安全格局的重要組成部分，位于長江流域與珠江流域的分水嶺及源頭區(qū)，主要為加強植被修復(fù)和水土流失防治，從而發(fā)揮華南和西南地區(qū)生態(tài)安全屏障作用，但目前涉及本區(qū)域的相關(guān)研究還十分匱乏[12]。因此，研究其十年（2000-2010年）間生態(tài)系統(tǒng)水源涵養(yǎng)功能及其變化特征，對了解其生態(tài)安全屏障的作用發(fā)揮與否以及生態(tài)環(huán)境工程的效益具有重要現(xiàn)實意義。本文在遙感（Remote Sensor,RS）與地理信息系統(tǒng)（Geographic Information System,GIS）的支持下，以2000、2005和2010年三期生態(tài)系統(tǒng)類型數(shù)據(jù)為基礎(chǔ)，采用降水貯存量法，替代工程法及轉(zhuǎn)移矩陣的方法，對該區(qū)域生態(tài)系統(tǒng)水源涵養(yǎng)功能及其變化特征進行分析，獲得水源涵養(yǎng)功能變化的基本特征及變化的方向，以期為該區(qū)將生態(tài)系統(tǒng)管理和生態(tài)環(huán)境綜合評估提供科學(xué)參考。

1 研究方法

1.1 研究區(qū)概況

研究區(qū)位于 102.5°-116.9°E，22.4°-26.7°N（圖 1），涉及云南、貴州、廣西、廣東、湖南、江西六省共114個縣/區(qū)，總面積28.85萬km2。地貌類型復(fù)雜多樣，平原臺地、丘陵、低山、中山和高山兼有，其中以丘陵為主（13.43萬km2，面積比例45.26%）。坡地以緩坡和斜坡為主（面積比例分別為36.10%和29.65%）。地勢西高東低，北高南低，海拔為-25 m-3040 m。河流眾多，河網(wǎng)密度大，屬于亞熱帶季風(fēng)氣候區(qū)，年平均氣溫在14-22℃之間，年降水量為1300-1610 mm，降雨分布不均勻，呈由東部、東南部向西部、西北部減少的趨勢。土壤類型以黃壤、紅壤、赤紅壤和磚紅壤為主，另有少量石灰土、紫色土水稻土和潮土[13]。

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

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

以2000年、2005年、2010年的遙感影像為基礎(chǔ)數(shù)據(jù)（來源于中國科學(xué)院遙感與數(shù)字地球研究所），經(jīng)幾何校正、影像融合、遙感解譯等數(shù)據(jù)處理，最終將地表分為森林、農(nóng)田、草地、灌叢、濕地、人工表面和稀疏地7種生態(tài)系統(tǒng)類型，其中人工表面和稀疏地不參與計算；年平均降雨數(shù)據(jù)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc.cma.gov.cn/home.do)。

1.3 計算與分析方法

1.3.1 水源涵養(yǎng)量 水源涵養(yǎng)功能物質(zhì)量估算采用降水貯存量法[14],即用生態(tài)系統(tǒng)的截留降雨和蓄水效應(yīng)來衡量其涵養(yǎng)水分的功能，生態(tài)系統(tǒng)的水源涵養(yǎng)功能可理解為，相對于裸地降水過程中產(chǎn)流的減少，即僅在降雨強度大于產(chǎn)流降雨時，其水源涵養(yǎng)功能才得以發(fā)揮[15]。計算公式如下：

式中：Q為與裸地相比較，森林、草地、濕地、耕地、灌叢等生態(tài)系統(tǒng)涵養(yǎng)水分的增加量(mm/hm2·a)；為生態(tài)系統(tǒng)面積(hm2)（表 1）；J為區(qū)內(nèi)多年均產(chǎn)流降雨量（P>20 mm）(mm)；Jo為區(qū)內(nèi)多年均降雨總量 (1601.74 mm)；K為區(qū)內(nèi)產(chǎn)流降雨量占降雨總量的比例(秦嶺-淮河以北取0.4，以南取0.6)；R為生態(tài)系統(tǒng)減少徑流的效益系數(shù)；Ro為產(chǎn)流降雨條件下裸地降雨徑流率；Rg為產(chǎn)流降雨條件下生態(tài)系統(tǒng)降雨徑流率。

1.3.2 水源涵養(yǎng)經(jīng)濟價值 采用替代工程法對研究區(qū)生態(tài)系統(tǒng)水源涵養(yǎng)量價值化，公式為[16]：

表1 研究區(qū)各生態(tài)系統(tǒng)面積

式中：Vw為區(qū)域生態(tài)系統(tǒng)涵養(yǎng)水源的價值，Vw(x)為象元處涵養(yǎng)水源的價值，Q(x)為象元x處水源涵養(yǎng)量，Pw(x)為建成單位庫容的成本，我國1 m3庫容的成本為0.67元。

1.3.3 生態(tài)系統(tǒng)水源涵養(yǎng)功能分級 將生態(tài)系統(tǒng)水源涵養(yǎng)量結(jié)果進行標(biāo)準(zhǔn)化，將標(biāo)準(zhǔn)化后的生態(tài)系統(tǒng)水源涵養(yǎng)功能評估單元劃分為高（0.8-1.0）、較高（0.6-0.8）、中（0.4-0.6）、較低（0.2-0.4）、低（0-0.2）五個等級。

式中：SSC表示標(biāo)準(zhǔn)化之后的生態(tài)系統(tǒng)水源涵養(yǎng)功能值；SCx表示各評價單元（此處為柵格）生態(tài)系統(tǒng)水源涵養(yǎng)量；SCmax和SCmin表示生態(tài)系統(tǒng)水源涵養(yǎng)量的最大值和最小值。

1.3.4 生態(tài)系統(tǒng)水源涵養(yǎng)功能變化 本文借助轉(zhuǎn)移矩陣分析區(qū)域生態(tài)系統(tǒng)水源涵養(yǎng)功能變化的結(jié)構(gòu)特征以及各級別水源涵養(yǎng)功能變化的方向，根據(jù)地圖代數(shù)原理，對任意兩期水源涵養(yǎng)功能等級數(shù)據(jù)的轉(zhuǎn)移矩陣計算公式如下：

該公式使用范圍類型小于10種，式中Ci×j為k時期到k+1時期水源涵養(yǎng)功能等級變化數(shù)據(jù)為 k時期水源涵養(yǎng)功能等級數(shù)據(jù)，為k+1時期水源涵養(yǎng)功能等級數(shù)據(jù)，i和j分別表示不同的水源涵養(yǎng)等級[17-19]。

2 結(jié)果與分析

2.1 生態(tài)系統(tǒng)水源涵養(yǎng)功能

研究區(qū)2000、2005和2010年水源涵養(yǎng)量空間分布大致呈現(xiàn)東高西低、東南高西北低的空間分布特征（圖2）。十年間水源涵養(yǎng)總量上升了1.61×108m3，經(jīng)濟價值上升了1.08×108元，三個年份水源涵養(yǎng)總量分別為856.40×107m3、864.96×107m3和872.47×107m3，水源涵養(yǎng)功能的經(jīng)濟價值分別為57.38×108元、57.95×108元和 58.46×108元（見表2）。不同生態(tài)系統(tǒng)水源涵養(yǎng)總量排序都是森林>農(nóng)田>灌叢>草地>濕地，其中森林水源涵養(yǎng)量比重大（三個年份分別69.17%、69.51%、69.91%），濕地生態(tài)系統(tǒng)由于面積比例?。?.19%、1.20%和1.24%）而水源涵養(yǎng)總量少（僅占3.84%、3.71%和3.66%）。不同生態(tài)系統(tǒng)的水源涵養(yǎng)總量變化趨勢不盡相同，其中森林和草地生態(tài)系統(tǒng)呈上升趨勢，濕地、農(nóng)田和灌叢生態(tài)系統(tǒng)不同程度下降。

不同生態(tài)系統(tǒng)的水源涵養(yǎng)能力除灌叢無明顯變化外，其余都發(fā)生了不同程度的變化，其中森林和草地生態(tài)系統(tǒng)水源涵養(yǎng)能力持續(xù)上升，濕地和農(nóng)田生態(tài)系統(tǒng)水源涵養(yǎng)能力呈下降趨勢，十年間，各生態(tài)系統(tǒng)水源涵養(yǎng)能力大小排序為：濕地>森林>草地>灌叢>農(nóng)田（表2）。

圖2 研究區(qū)水源涵養(yǎng)量空間分布圖

表2 研究區(qū)不同生態(tài)系統(tǒng)水源涵養(yǎng)量

2.2 生態(tài)系統(tǒng)水源涵養(yǎng)功能分級

研究區(qū)生態(tài)類型水源涵養(yǎng)功能高級別以濕地生態(tài)系統(tǒng)為主（表3），占全區(qū)總面積的比例較小，三個年份分別為1.18%、1.19%和1.23%，這主要是由于濕地生態(tài)系統(tǒng)所占的面積最?。?1.3%）導(dǎo)致；水源涵養(yǎng)功能在較高級別的生態(tài)系統(tǒng)以森林生態(tài)系統(tǒng)中的常綠闊葉林及濕地生態(tài)系統(tǒng)中的灌叢沼澤和草本沼澤為主，三者之和面積比例在三個年份分別為24.86%、24.92%和24.95%（表4）；水源涵養(yǎng)功能在中級別的生態(tài)系統(tǒng)主要有森林生態(tài)系統(tǒng)中的落葉闊葉林、常綠針葉林和針闊混交林以及草地生態(tài)系統(tǒng)，其區(qū)內(nèi)分布范圍廣泛；水源涵養(yǎng)功能在較低級別的生態(tài)系統(tǒng)以農(nóng)田生態(tài)系統(tǒng)為主，三個年份占全區(qū)總面積的比例分別為22.17%、21.89%和21.76%；低級別的生態(tài)系統(tǒng)分布較零散，其面積占全區(qū)總面積的比例較小，主要包括居住地、工業(yè)用地、交通用地、裸土等硬質(zhì)地面或植被稀少的地貌。

2.3 生態(tài)系統(tǒng)水源涵養(yǎng)功能變化

2000-2005年和2005-2010年這兩個時段中，水源涵養(yǎng)功能高等級轉(zhuǎn)換成較低等級的面積分別為87.22 km2和202.49 km2，轉(zhuǎn)換比率分別為2.6%和5.8%，轉(zhuǎn)換成低等級的面積分別為10.09 km2和25.30 km2，轉(zhuǎn)換比率分別為0.3%和0.7%。表明，受人類活動的影響濕地生態(tài)系統(tǒng)有部分面積被轉(zhuǎn)換成了耕地和人工表面，且這種現(xiàn)象有持續(xù)發(fā)展的趨勢。

表3 研究區(qū)生態(tài)類型水源涵養(yǎng)量標(biāo)準(zhǔn)化值

表4 研究區(qū)水源涵養(yǎng)功能各等級面積及其比例

兩個時段中低等級的水源涵養(yǎng)功能大量轉(zhuǎn)換成較低和較高等級（占其自身轉(zhuǎn)出率的63.13%和59.18%），其中轉(zhuǎn)換成較低等級的面積分別為79.60 km2和181.26 km2，轉(zhuǎn)換成較高等級的面積分別為106.66 km2和201.96 km2。有可能是部分裸土地被開發(fā)成耕地用于農(nóng)業(yè)生產(chǎn)，而稀疏地的植被逐漸成長，使得原本水源涵養(yǎng)功能低的稀疏地轉(zhuǎn)換成涵養(yǎng)功能較高的林地。受國家退耕還林，還草的政策影響，研究區(qū)部分農(nóng)田轉(zhuǎn)換成林地或草地。較低等級轉(zhuǎn)換成中等級的面積分別為2738.46 km2和1974.65 km2，轉(zhuǎn)換成較高等級的面積分別為1063.73km2和698.44km2，占自身轉(zhuǎn)出率的84.68%和79.82%。

從2000-2010總體上看，十年間，研究區(qū)水源涵養(yǎng)功能的變化方向主要是向中和較高等級轉(zhuǎn)移，且各生態(tài)系統(tǒng)以向森林生態(tài)系統(tǒng)轉(zhuǎn)換為主。除自身轉(zhuǎn)移的面積（分別為123272.87 km2和69715.72 km2）外，轉(zhuǎn)入的面積分別為8356.23 km2和 5829.38 km2。

3 討論與結(jié)論

3.1 討論

南方丘陵山地帶是“兩屏三帶”中的重要組成部分，其生態(tài)環(huán)境的好壞，直接影響著華南和西南地區(qū)的生態(tài)安全。本文分析結(jié)果一定程度表明，國家實施退耕還林、還草等生態(tài)工程取得了重大的效益。區(qū)內(nèi)生態(tài)系統(tǒng)水源涵養(yǎng)能力整體得到了提升，生態(tài)環(huán)境得到了進一步改善。但須指出的是區(qū)內(nèi)水源涵養(yǎng)功能變化依然存在由高往低轉(zhuǎn)換的現(xiàn)象，有些生態(tài)系統(tǒng)如濕地生態(tài)系統(tǒng)仍然沒有得到有力的保護。

有效分析區(qū)內(nèi)生態(tài)系統(tǒng)水源涵養(yǎng)功能是建立長效生態(tài)補償機制的基礎(chǔ)，可為進一步生態(tài)建設(shè)提供保障。上世紀(jì)九十年代起，國家在云南、貴州、廣西和廣東開展了長江、珠江流域防護林體系，1999年，區(qū)內(nèi)廣西、云南、貴州作為西部省份納入退耕還林工程建設(shè)范圍，2001年湖南和江西作為中部省份納入退耕還林工程范圍[20]。這一系列措施直接影響土地利用方式及土地覆被結(jié)構(gòu)變化。本文對水源涵養(yǎng)功能的變化方向進行分析，一方面反映了水源涵養(yǎng)功能變化的趨勢，另一方面也初步反映了土地覆被類型間的轉(zhuǎn)換對水源涵養(yǎng)功能變化的影響。但文中涉及的生態(tài)系統(tǒng)類型較多，只能粗略對其說明，在以后的研究中，可進一步結(jié)合土地覆被類型的轉(zhuǎn)換分析其對水源涵養(yǎng)能力的影響強度，為優(yōu)化土地利用結(jié)構(gòu)提供參考。

表5 不同級別水源涵養(yǎng)功能生態(tài)系統(tǒng)轉(zhuǎn)移矩陣（單位：km2）

不同生態(tài)系統(tǒng)水源涵養(yǎng)能力與總量是不相同的，從而它們對整個區(qū)域水源涵養(yǎng)功能的貢獻率也有所區(qū)別。本文研究表明，研究區(qū)不同生態(tài)系統(tǒng)水源涵養(yǎng)總量大小排序為森林>農(nóng)田>灌叢>草地>濕地，生態(tài)系統(tǒng)水源涵養(yǎng)能力大小排序為：濕地>森林>草地>灌叢>農(nóng)田。雖然生態(tài)系統(tǒng)水源涵養(yǎng)總量大小與其涵養(yǎng)能力大小并無緊密關(guān)系，但森林生態(tài)系統(tǒng)依然是水源涵養(yǎng)功能的主要貢獻者（69.17%、69.51%、69.91%），因此森林生態(tài)系統(tǒng)仍是重點保護對象。

水源涵養(yǎng)功能物質(zhì)量估算的方法不盡相同，用不同方法估算的結(jié)果存在一定差異。本研究采用降水貯存量法對區(qū)域生態(tài)系統(tǒng)水源涵養(yǎng)功能物質(zhì)量進行估算，并在此基礎(chǔ)上對水源涵養(yǎng)功能變化特征進行分析。該方法所需數(shù)據(jù)容易獲取，計算過程也較簡單，從分析結(jié)果來看，能有效的反映區(qū)域水源涵養(yǎng)功能變化特征，但由于沒有考慮到不同生態(tài)系統(tǒng)對此方法的限制性因素，一定程度上影響了物質(zhì)量估算結(jié)果的精度。因此今后可針對不同生態(tài)系統(tǒng)采用多種方法進行計算，以便更好的研究水源涵養(yǎng)功能及其經(jīng)濟價值。

3.2 結(jié)論

本文通過分析得到以下主要結(jié)論：2000-2010年間，研究區(qū)水源涵養(yǎng)量空間分布呈現(xiàn)東高西低、東南高西北低的特征，與降雨變化趨勢基本一致。其水源涵養(yǎng)總量呈上升趨勢，十年間總量上升了1.61×108m3，其上升的經(jīng)濟價值為1.08×108元。水源涵養(yǎng)功能主要集中在較高、中和較低等級，三者總面積比例三個年份分別為97.19%、97.00%和96.87%，水源涵養(yǎng)功能的變化方向主要是向中和較高等級轉(zhuǎn)移，整體呈現(xiàn)由低往高的變化特征。

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

[2]秦嘉勵,楊萬勤,張健.岷江上游典型生態(tài)系統(tǒng)水源涵養(yǎng)量及價值評估[J].應(yīng)用與環(huán)境生物學(xué)報,2009,15(4):453-458.

[3]姜昊,程磊磊,尹昌斌.洪澤湖濕地生態(tài)系統(tǒng)服務(wù)功能價值評估研究[J].農(nóng)業(yè)現(xiàn)代化研究,2008,29(3):331-334.

[4]肖寒,歐陽志云,趙景柱,等.森林生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟價值評估初探—以海南島尖峰嶺熱帶森林為例[J].應(yīng)用生態(tài)學(xué)報,2000,11(4):481-484.

[5]馬雪華.森林水文學(xué)[M].北京:中國林業(yè)出版社,1993.

[6]候元兆,張穎,曹克瑜.森林資源核算(上卷)[M].北京:中國科學(xué)技術(shù)出版社,2005.

[7]張三煥,趙國柱,田允哲,等.長白山琿春林區(qū)森林資源資產(chǎn)生態(tài)環(huán)境價值的評估研究[J].延邊大學(xué)學(xué)報 (自然科學(xué)版),2001,27(2):126-134.

[8]Erik Nelson,GuillermoMendoza,JamesRegetz,etal.Modelingmultiple ecosystem services,biodiversity conservation,commodity production,and tradeoffs at landscape scales[J].Frontiers in Ecology and the Environment,2009,7(1):4-11.

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

[10]Peter Kareiva,Heather Tallis,Taylor HRicketts,etal.Natural Capital:Theory and Practice of Mapping Ecosystem Services[M].London:Oxford University Press,2011.

[11]楊愛霞.甘家湖濕地邊緣帶景觀格局變化及生態(tài)功能價值研究[D].烏魯木齊:新疆師范大學(xué),2012.

[12]歐陽貝思,張明陽,王克林,等.2000-2010年南方丘陵山地帶土地覆被及景觀格局變化特征[J].農(nóng)業(yè)現(xiàn)代化研究,2013,34(4):467-471.

[13]陳勇,王亞麗.我國東南丘陵區(qū)土壤侵蝕現(xiàn)狀及防治對策[J].熱帶農(nóng)業(yè)科學(xué),2003,23(1):38-42.

[14]李金昌.生態(tài)價值論[M].重慶:重慶大學(xué)出版社.1999

[15]陽柏蘇.景區(qū)土地利用格局及生態(tài)系統(tǒng)服務(wù)功能研究—以張家界國家森林公園為例（1990-2000）[D].長沙:中南林業(yè)科技大學(xué),2005.

[16]張明陽.基于景觀格局分析的喀斯特區(qū)生態(tài)服務(wù)功能研究-以桂西北為例[D].北京:中國科學(xué)院,2009.

[17]李雁,趙坤,方精云,等.城市湖泊地區(qū)的土地利用變化—以武漢東湖為例[J].長江流域資源與環(huán)境,2004,13(3):229-233.

[18]Shi Yusheng,Xiao Jieying,Shen Yanjun,etal.Quantifyingthespatial differences of landscape change in the Hai River Basin,China,in the1990s[J].International Jouranl of Remote Sensing.2012,33(14):4482-4501.

[19]黃浦江,劉艷芳,劉暢,等.基于RS與GIS的武漢城市湖泊演化研究[J].生態(tài)環(huán)境學(xué)報,2012,21(9):1588-1593.

[20]歐陽貝思.南方丘陵山地帶十年（2000-2010）土地覆被格局變化及驅(qū)動機制研究[D].長沙:湖南農(nóng)業(yè)大學(xué),2013.