歐朝蓉,朱清科,孫永玉

(1.北京林業(yè)大學水土保持學院，北京 100083； 2.西南林業(yè)大學生態(tài)旅游學院，云南 昆明 650024；3.中國林業(yè)科學研究院資源昆蟲研究所，云南 昆明 650216)

元謀干熱河谷林地生態(tài)系統(tǒng)服務價值變化

歐朝蓉1,2,朱清科1,孫永玉3*

(1.北京林業(yè)大學水土保持學院，北京 100083； 2.西南林業(yè)大學生態(tài)旅游學院，云南 昆明 650024；3.中國林業(yè)科學研究院資源昆蟲研究所，云南 昆明 650216)

目的以元謀縣為研究區(qū)，基于生態(tài)學、RS和GIS理論和方法，研究干熱河谷林地利用及生態(tài)系統(tǒng)服務價值的變化特征，以期豐富干熱河谷林地生態(tài)研究，為區(qū)域生態(tài)環(huán)境調(diào)控提供理論依據(jù)。方法以2008年Landsat 7 ETM及2016年Landsat 8 OLI 2期遙感影像、國家基礎地理信息系統(tǒng)及ASTER GDEM數(shù)字高程模型為數(shù)據(jù)源，利用ENVI和ArcGIS 軟件對2期遙感影像進行解譯，獲取林地土地利用數(shù)據(jù)，分析林地的時空特征；調(diào)整林地生態(tài)系統(tǒng)服務當量因子系數(shù)，計算單位面積農(nóng)田食物生產(chǎn)功能價值及單位面積林地生態(tài)系統(tǒng)服務價值系數(shù)，分析林地生態(tài)系統(tǒng)靜態(tài)價值特征及其變化；以支付能力指數(shù)、支付意愿指數(shù)和環(huán)境能力指數(shù)為參數(shù)構(gòu)建林地生態(tài)系統(tǒng)服務動態(tài)價值調(diào)整系數(shù)，分析林地生態(tài)系統(tǒng)服務動態(tài)價值特征及其變化。結(jié)果林地主要分布在研究區(qū)南部和東部的中高山，在西部、西南部和北部及東部邊緣的中高山有少許分布，在河谷壩區(qū)有零星分布；2008—2016年林地面積增長了64.24 hm2，變化幅度達到11.42%；有90.53 km2的林地轉(zhuǎn)為非林地，有154.77 km2非林地轉(zhuǎn)為林地。林地向草地轉(zhuǎn)移了73.61 km2，草地向林地轉(zhuǎn)移123.11 km2；2016年研究區(qū)林地生態(tài)系統(tǒng)服務靜態(tài)價值為14.97億元，較2008年的13.44億元提高了1.53億元。2016年林地在研究區(qū)生態(tài)系統(tǒng)總靜態(tài)價值的比例為56.41%，較2008年的51.73%提高了4.68%；2008年和2016年林地生態(tài)系統(tǒng)服務動態(tài)價值分別為4.758億元和9.027億元，僅為同期林地生態(tài)系統(tǒng)服務靜態(tài)價值的35.40%和60.30%。結(jié)論研究區(qū)林地面積迅速增長，林地與非林地相互轉(zhuǎn)換頻繁，林地凈增長明顯；林地面積的增長提升了生態(tài)系統(tǒng)服務的靜態(tài)價值。林地的單項生態(tài)系統(tǒng)服務價值中，土壤形成與保護、氣體調(diào)節(jié)和生物多樣性保護的價值最大，食物生產(chǎn)的價值最小；林地生態(tài)系統(tǒng)服務動態(tài)價值雖然遠低于靜態(tài)價值，但其變化幅度明顯高于靜態(tài)價值。高生態(tài)系統(tǒng)服務價值的林地對社會和環(huán)境調(diào)節(jié)能力的增長更為敏感。

干熱河谷；林地；生態(tài)系統(tǒng)服務價值；當量因子；靜態(tài)價值；動態(tài)價值

Abstract:[Objective]Yuanmou county was selected to studied the change of forest land utilization and forest ecosystem service value in the dry-hot valley so as to enrich forest ecology research and provide theoretical basis for regional ecological environment regulation.[Method]Taking remote images of Landsat 7 ETM in 2008 and Landsat 8 OLI in 2016, national basic geographic information system and ASTER GDEM digital elevation model as data sources, the ENVI and ArcGIS software was used to interpret the remote sensing images, obtain the forest land-use data, and to analyze the temporal and spatial characteristics of forest land, and calculate the function value of farmland food production per unit area and coefficients per unit area of forest ecosystem service value. Then, the forest ecosystem service value and its changing characteristics were analyzed. The adjustment of coefficient of dynamic value of forest ecosystem service was constructed on the basis of payment ability, willingness to pay and land by using ENVI and ArcGIS software. Forest ecosystem equivalent coefficients were adjusted, as well as environmental capacity as main parameters in order to analyze the dynamic value characteristic and its change of forest ecosystem service.[Result]Forest land is mainly distributed in southern and eastern parts of the middle-high mountains in the study area. A small amount of forest scatters in western, southwest, north and eastern foothills of the middle-high mountains. Few parts are sporadic in the valley dam area. Forest area increased by 64.24 hm2, and variation range reached 11.42% from 2008 to 2016. The area of forest land converted into non-forestland reached 90.53 km2, while the forest area converted from other land-use types was 154.77 km2. There waere 73.61 km2of forest land converting to grassland and 123.11 km2converted from grassland. In 2016, the static value of forest ecosystem service in the study area was ￥1.497 billion, an increase of ￥153 million over ￥1.344 billion in 2008. In 2016, the proportion of forest land in the total static value of regional ecosystem service in the study area was 56.41%, which increased by 4.68% compared with 51.73% in 2008. The dynamic value of forest ecosystem service in 2008 and 2016 were respectively ￥ 475.8 million and ￥ 902.7 million, which only accounted for 35.40% and 60.30% of static value of forest ecosystem service in the same period. [Conclusion]The forest land area increased rapidly, the forest land and non-forest land converted to each other frequently, leading to an obvious net increase of forest area. The area increase of forest land resulted in static value growth of forest ecosystem service. Among individual value of forest ecosystem service, soil formation and protection, gas regulation and biodiversity conservation scored the highest, while the food production had the lowest score. Although the dynamic value of forest ecosystem service is far less than the static value, the range of dynamic value change is obviously higher than that of static value. The forest land with high ecosystem service value is more sensitive to the growth of social and environmental regulation.

Keywords: dry-hot valley; forest land; ecosystem service value; equivalent factor; static value; dynamic value

生態(tài)系統(tǒng)服務價值是生態(tài)系統(tǒng)服務功能的貨幣化形式。區(qū)域生態(tài)系統(tǒng)服務價值越高在一定程度上表明區(qū)域生態(tài)服務功能越好，人們的生態(tài)意識越強，越有利于區(qū)域生態(tài)安全。20世紀90年代以來，隨著人類對生態(tài)系統(tǒng)功能認識的加強，生態(tài)系統(tǒng)服務價值研究受到廣泛關(guān)注，最具代表性的是Costanza等[1]提出的“生態(tài)系統(tǒng)服務價值與自然資本”理論，然而其研究中的“全球靜態(tài)部分平衡模型”是基于全球生態(tài)系統(tǒng)服務價值均值化的處理，在地域性生態(tài)系統(tǒng)服務價值評估中會產(chǎn)生較大誤差。謝高地等[2]在Costanza等[1]研究基礎上結(jié)合中國的自然特征制定了青藏高原單位面積生態(tài)系統(tǒng)服務價值當量因子,成為國內(nèi)生態(tài)系統(tǒng)服務價值研究的理論基礎和主要方法參考。近30年來，國內(nèi)外學者從景觀尺度、單一生態(tài)系統(tǒng)、物種和生物多樣性等方面開展了生態(tài)系統(tǒng)服務價值多角度研究[3-8],不僅夯實了生態(tài)系統(tǒng)功能及價值評估的理論基礎，同時也較有力的推動了該領(lǐng)域研究方法的創(chuàng)新[9-15]。森林是陸地生態(tài)系統(tǒng)的重要組成部分，在維護區(qū)域生態(tài)平衡中具有關(guān)鍵作用。國外研究側(cè)重于林地的單一生態(tài)服務功能的測評及森林生態(tài)系統(tǒng)服務與局地生態(tài)環(huán)境的關(guān)系分析等方面[16-20]，國內(nèi)研究則側(cè)重于從土地利用角度評估區(qū)域生態(tài)系統(tǒng)服務價值，王麗等[21]從涵養(yǎng)水源價值、土壤保持價值、固定CO2價值、凈化空氣價值、社會與文化價值5個方面評估了青城生態(tài)系統(tǒng)服務價值，鄭江坤等[22]利用馬爾科夫鏈模型評估了潮白河流域森林生態(tài)服務價值的動態(tài)變化，馬國軍等[23]和汪有奎等[24]對干旱和半干旱區(qū)的森林生態(tài)系統(tǒng)服務價值進行了評價。

干熱河谷是中國西南地區(qū)一種局地特殊的地理景觀和氣候類型，氣候炎熱、干燥，水分蒸發(fā)比極為失衡，土壤侵蝕和水土流失嚴重，植被類型為“河谷型薩瓦納植被”，特殊的自然生態(tài)環(huán)境使干熱河谷成為西南地區(qū)典型的生態(tài)脆弱區(qū)。林地是干熱河谷山地生態(tài)系統(tǒng)的重要組成部分，具有水源涵養(yǎng)和氣候調(diào)節(jié)等多種功能，是維護干熱河谷區(qū)生態(tài)系統(tǒng)功能的主要地類，林地的數(shù)量及時空分布對干熱河谷區(qū)域生態(tài)安全的影響顯著。現(xiàn)有的干熱河谷林地研究側(cè)重于林木育苗[25-26]、種群結(jié)構(gòu)[27]和植被恢復[28-29]等方面，這些研究雖從單一角度部分探討了林地的生態(tài)功能，但是無法有效系統(tǒng)量化林地生態(tài)系統(tǒng)服務價值。盡管有學者曾研究了林地生態(tài)系統(tǒng)服務靜態(tài)價值[30]，但沒能深入評估林地生態(tài)系統(tǒng)動態(tài)價值的變化，不能揭示社會、經(jīng)濟和環(huán)境因素對干熱河谷林地生態(tài)系統(tǒng)服務價值的影響。鑒于此，本文以具有干熱河谷典型環(huán)境特征的元謀縣為案例，研究干熱河谷林地的時空格局及生態(tài)系統(tǒng)服務靜態(tài)價值和動態(tài)價值的變化特征，辨明影響林地生態(tài)系統(tǒng)服務價值變化的主要因素，以期為干熱河谷林地可持續(xù)發(fā)展及區(qū)域生態(tài)環(huán)境保護提供科學依據(jù)。

1 研究區(qū)概況

圖1 研究區(qū)區(qū)位示意及地形圖Fig.1 Sketch map of geographical location and topography of the study area

2 研究方法

2.1數(shù)據(jù)源及數(shù)據(jù)預處理

本研究的數(shù)據(jù)源為：(1)遙感數(shù)據(jù)，選取2008年Landsat7ETM及2016年Landsat8OLI衛(wèi)星遙感影像，過境日期分別為2008年1月28日和2016年1月4日，分辨率為30m。兩景遙感影像時段均處于研究區(qū)旱季。數(shù)據(jù)由中國科學院地理空間數(shù)據(jù)云平臺提供。(2)地理信息系統(tǒng)數(shù)據(jù)，選取研究區(qū)1∶100000的國家地理基礎信息系統(tǒng)數(shù)據(jù)，包括水系、道路和居民點等基本信息。數(shù)據(jù)由國家基礎地理信息中心提供。(3)DEM數(shù)據(jù)，選取ASTER GDEM V2產(chǎn)品，分辨率為30m。研究區(qū)共涉及到4個圖幅(ASTGTM2_N26E101、ASTGTM2_N26E102、ASTGTM2_N25E101、ASTGTM2_N25E102)。數(shù)據(jù)由中國科學院地理空間數(shù)據(jù)云平臺提供。

在ArcGIS中利用Mosaic工具對4個圖幅的DEM進行拼接，拼接圖幅包含了研究區(qū)區(qū)域范圍。利用水系和高程點為控制點完成拼接的DEM圖幅與1∶100000基礎地理信息系統(tǒng)數(shù)據(jù)的空間配準，以研究區(qū)范圍的矢量數(shù)據(jù)裁剪拼接的DEM圖幅獲得研究區(qū)范圍DEM，利用ArcGIS柵格計算功能提取研究區(qū)4個垂直自然帶的空間范圍，利用ENVI軟件完成兩景遙感影像的去云、去條帶、標準假彩色合成及幾何精校正等圖像處理工作。完成2個時段遙感影像與1∶100000地形圖的空間配準，并以研究區(qū)矢量邊界為掩膜對裁切遙感影像獲得了研究區(qū)范圍內(nèi)的遙感影像。

2.2土地利用分類體系

根據(jù)《國家森林資源連續(xù)清查技術(shù)規(guī)定》地類劃分標準和研究區(qū)遙感影像特征，結(jié)合前人研究成果[30-33]和實際研究需要，將研究區(qū)分為林地和非林地2個大類別，非林地包括草地、耕地、建設用地、未利用地和水域5種一級土地利用類型。

2.3土地利用數(shù)據(jù)的提取

綜合運用非監(jiān)督分類、監(jiān)督分類及目視解譯等方法進行遙感圖像解譯。以土地利用類型歷史圖件與GPS野外采點調(diào)研驗證了土地利用類型的解譯精度。2期遙感圖像的解譯精度均在88%以上，符合解譯精度要求(表1)。在ENVI軟件中以像元二分法提取研究區(qū)2008年和2016年的植被覆蓋度。利用ArcGIS軟件提取主要土地利用類型的基本統(tǒng)計數(shù)據(jù)，計算林地土地利用指數(shù)和土地狀態(tài)指數(shù)[34]，利用Dissolve和Intersect工具對2008年和2016年2個時段土地利用類型進行疊置分析，得到研究區(qū)土地利用轉(zhuǎn)移矩陣。

表1 土地利用類型解譯精度

2.4林地生態(tài)系統(tǒng)服務靜態(tài)價值系數(shù)調(diào)整

2.4.1林地生態(tài)系統(tǒng)服務當量因子系數(shù)調(diào)整 謝高地等[2]的研究以一級土地利用類型為研究基礎，將不同的林地類型歸為一種，不能有效揭示植被種類之間的差異性導致的生態(tài)系統(tǒng)服務價值的差別。李曉賽等[35]在方精云等[36]對植被生物量差異性研究的基礎上，基于材積源生物量法提出了有林地、灌木林和其它林地的生態(tài)系統(tǒng)服務當量因子系數(shù)的調(diào)整方案，取3種林地的功能性調(diào)整系數(shù)平均值為林地生態(tài)系統(tǒng)服務靜態(tài)價值當量因子調(diào)整系數(shù)，即1.032(表2)。研究區(qū)的林地主要類型為有林地、灌木林和其它林地。為了更好地反映干熱河谷不同林地的生態(tài)系統(tǒng)服務價值的差異性，本研究采用李曉賽等[35]的研究方法設定研究區(qū)林地生態(tài)系統(tǒng)服務當量因子的調(diào)整系數(shù)為1.032。

表2 林地生態(tài)系統(tǒng)服務當量因子系數(shù)調(diào)整

注：本表來源于李曉賽等[35]和方精云[36]等的研究成果。

2.4.2 單位面積食物生產(chǎn)功能價值當量因子 單位面積食物生產(chǎn)功能價值當量因子是指1 hm2農(nóng)田每年自然糧食的平均產(chǎn)量的產(chǎn)值, 林地生態(tài)系統(tǒng)服務價值當量因子則以林地生態(tài)系統(tǒng)產(chǎn)生該生態(tài)服務相對于農(nóng)田食物生產(chǎn)服務貢獻的大?。?/p>

Ea=1/7×Ti/Mi

(1)

式中：Ea為單位面積農(nóng)田生態(tài)系統(tǒng)提供食物生產(chǎn)服務功能的經(jīng)濟價值(元·hm-2)；Ti是研究區(qū)每年糧食作物的經(jīng)濟產(chǎn)值(元·a-1)；Mi是研究區(qū)糧食作物的總面積。1/7是指在沒有人工干預的自然生態(tài)系統(tǒng)提供的經(jīng)濟價值與現(xiàn)有單位面積農(nóng)田提供的食物生產(chǎn)服務經(jīng)濟價值的比值。以2008年和2016年研究區(qū)農(nóng)田糧食作物(水稻、玉米和雜糧)的價格為生態(tài)系統(tǒng)服務價值系數(shù)計算的基本依據(jù)，依據(jù)統(tǒng)計資料計算得到2個年份的單價分別為979.36、1 142.55元·hm-2，取二者的平均值1 060.96元·hm-2作為單位面積農(nóng)田食物生產(chǎn)功能價值。結(jié)合調(diào)整后的林地生態(tài)系統(tǒng)服務當量因子，通過公式(1)計算研究區(qū)林地及其它地類單位面積生態(tài)系統(tǒng)服務價值。

2.5林地生態(tài)系統(tǒng)動態(tài)服務價值系數(shù)調(diào)整

經(jīng)濟學研究表明, 對物品的支付水平與區(qū)域內(nèi)某時段人們的支付能力和基于物品的需求程度的支付意愿有關(guān)[35]，同時由于研究區(qū)所提供生態(tài)資源并非一成不變，而是受生態(tài)環(huán)境變化的影響，因此，本研究將支付能力、支付意愿[37]和環(huán)境能力等指標納入生態(tài)系統(tǒng)服務動態(tài)價值評估體系，對動態(tài)價值進行系數(shù)調(diào)整，獲得研究區(qū)生態(tài)系統(tǒng)服務動態(tài)價值系數(shù)。

2.5.1支付能力指數(shù)(Pi) 支付能力是個人經(jīng)濟能力的直接表現(xiàn)。個人支付能力越強，越有可能承受對生態(tài)系統(tǒng)服務功能的支付。本研究將研究區(qū)人均國內(nèi)生產(chǎn)總值與當年中國人均國內(nèi)生產(chǎn)總值的比值作為支付能力指數(shù)，其計算公式如下:

Pi=GDPi/GDPimean(i=1,2, …,10)

(2)

式中：GDPi表示第i年元謀縣人均國內(nèi)生產(chǎn)總值；GDPimean表示第i年中國人均國內(nèi)生產(chǎn)總值。利用統(tǒng)計數(shù)據(jù)通過公式(2)計算研究區(qū)2008年和2016年生態(tài)系統(tǒng)服務價值支付能力指數(shù)。

2.5.2支付意愿指數(shù)(Ni) 人們對生態(tài)系統(tǒng)價值的認識和需求通常與社會發(fā)展階段相關(guān)[35]。社會發(fā)展階段越高，人們的生態(tài)環(huán)境意識及對生態(tài)服務功能的需求越多(直至飽和狀態(tài))，支付意愿也越強。因此，本研究對支付意愿指數(shù)(Ni)進行調(diào)整以反映社會發(fā)展階段對生態(tài)系統(tǒng)服務價值的動態(tài)影響。支付意愿通常用Logistics生長曲線模型來刻畫，其表達式為：

(3)

式中：Ni為代表生長特性的參數(shù)，在此表示第i年社會發(fā)展階段系數(shù)；t在此表示社會經(jīng)濟發(fā)展階段；a、b為常數(shù)取值為1，e為自然對數(shù)。當t值很小時，即社會發(fā)展水平很低，Ni值趨于0；當t值很大時，社會發(fā)展水平很高，Ni值趨于飽和值[35]。社會發(fā)展階段的計算通常與反映生活支出結(jié)構(gòu)的恩格爾系數(shù)(En)的倒數(shù)對應起來，表達式如下：

t=1/En-2.5

(4)

利用統(tǒng)計數(shù)據(jù)由公式(4)計算得到研究區(qū)2008年和2016年的t值分別為-0.46、-0.32，以公式(3) 計算研究區(qū)2008年和2016年生態(tài)系統(tǒng)服務價值支付意愿指數(shù)。

2.5.3環(huán)境能力指數(shù)(Si) 植被的生物量不僅與植被的類型有關(guān)，同時受當年水分和熱量等環(huán)境條件的影響。植被覆蓋度是反映植被生長狀況的主要指標，能較客觀的反映植被的生產(chǎn)力和生物量的變化。本研究以植被覆蓋度作為基本參數(shù)構(gòu)建環(huán)境能力指數(shù)(Si)反映研究區(qū)生態(tài)系統(tǒng)服務價值因環(huán)境條件而產(chǎn)生的動態(tài)變化，其公式為：

(5)

式中：Sd表示該年度研究區(qū)植被覆蓋度，Smean表示研究區(qū)多年植被覆蓋度均值。根據(jù)遙感圖像運算結(jié)果，以公式(5)計算研究區(qū)2008年和2016年的生態(tài)系統(tǒng)服務價值環(huán)境能力指數(shù)。

基于上述研究成果，同時參考相關(guān)[37-40]研究，提出生態(tài)系統(tǒng)服務動態(tài)價值調(diào)整系數(shù)(Mij)，在此基礎上構(gòu)建研究區(qū)生態(tài)系統(tǒng)服務動態(tài)價值評估模型(Qij)，其具體模型如下：

Mij=Pij×Nij×Sij

(6)

(7)

式中：i為第i個年份，j為第j種土地利用類型；Pij表示居民在第i年第j種土地利用類型生態(tài)系統(tǒng)服務價值的支付能力指數(shù)；Nij表示居民在第i年第j種土地利用類型生態(tài)系統(tǒng)服務價值的支付意愿指數(shù)；Sij表示第i年第j種土地利用類型生態(tài)系統(tǒng)服務價值的環(huán)境能力指數(shù)；Qij為研究區(qū)第i年第j種土地利用類型生態(tài)系統(tǒng)服務動態(tài)價值；Cij為第i年第j種土地利用類型的面積；Eij是第i年第j種土地利用類型單位面積生態(tài)系統(tǒng)服務靜態(tài)價值(元·hm-2)；通過公式(6)計算2008年和2016年研究區(qū)生態(tài)系統(tǒng)服務價值動態(tài)調(diào)整系數(shù)值，林地和其它地類被賦予相同的生態(tài)系統(tǒng)服務動態(tài)價值調(diào)整系數(shù)值。利用公式(7)分別計算研究區(qū)2008和2016年及林地和其它地類生態(tài)系統(tǒng)服務動態(tài)價值及區(qū)域生態(tài)系統(tǒng)服務總動態(tài)價值。

3 結(jié)果與分析

3.1林地利用變化

圖2 研究區(qū)遙感分類結(jié)果Fig.2 Remote sensing classification result of the study area

3.1.2林地的時間變化 通過本研究的數(shù)據(jù)源獲得的研究區(qū)范圍面積為2032.84km2。2008—2016年，研究區(qū)林地的面積從2008年的562.08km2增長至2016年的626.32km2，增長了64.24km2，變化幅度達到11.42%，年變化率達1.42%。林地占研究區(qū)總面積比例由2008年的27.65%增至30.81%，增長了3.16%，增幅較為明顯。研究區(qū)林地的土地利用綜合指數(shù)由2008年的55.30增至60.20，增長了4.90，表明林地土地利用指數(shù)有所增長且增幅較大。

從表3中可以看出：林地與非林地之間相互轉(zhuǎn)化較為頻繁。林地向非林地的轉(zhuǎn)移面積為90.53km2，由非林地轉(zhuǎn)入為林地的面積為154.77km2。林地主要與非林地中的草地、耕地和未利用地和建設用地相互轉(zhuǎn)化。林地轉(zhuǎn)為草地的面積最多，達到73.61km2；其次是林地向耕地和建設用地的轉(zhuǎn)移，分別達到14.02、2.9km2。非林地向林地的轉(zhuǎn)移中，草地的轉(zhuǎn)移面積最大，為123.11km2；其次是耕地和未利用地，分別向林地轉(zhuǎn)移27.22、4.44km2。林地雖與草地相互轉(zhuǎn)化量最大，但林地以轉(zhuǎn)入為主，有49.5km2的凈增長，主要原因在于研究區(qū)林地以灌木林為主，有林地少。灌木林和灌草叢的生境相似，易受水分條件的波動及人為干擾的影響而相互轉(zhuǎn)換，發(fā)展人工林使更多的草地轉(zhuǎn)為林地；林地與耕地的相互轉(zhuǎn)化中，林地以轉(zhuǎn)出為主，凈減少13.2km2，主要受耕地空間分布變化的影響。退耕還林和天然林工程及水土保持工作的開展使山區(qū)坡度≥25°的耕地轉(zhuǎn)為林地，而壩周低山區(qū)和河谷壩區(qū)地勢較平，坡度小，土地開發(fā)難度相對較小，農(nóng)業(yè)經(jīng)濟發(fā)展對耕地的強烈需求使更多的林地開發(fā)成耕地，因而林地在與耕地的相互轉(zhuǎn)化中表現(xiàn)為凈減少；林地與未利用地的相互轉(zhuǎn)化中，林地以轉(zhuǎn)入為主，凈增長4.44km2，主要是由于經(jīng)濟利益驅(qū)使的土地開發(fā)和生態(tài)恢復工程的實施使未利用地向林地轉(zhuǎn)化，因而林地獲得一定凈增長；林地與建設用地的相互轉(zhuǎn)化中，林地以轉(zhuǎn)出為主，凈減少2.9km2，這主要是由于交通運輸?shù)陌l(fā)展使部分林地轉(zhuǎn)化為道路，居民用地增長使部分林地轉(zhuǎn)為城鄉(xiāng)居民用地，因而林地表現(xiàn)為一定的凈減少。林地的土地利用狀態(tài)指數(shù)為0.26，指數(shù)為正且數(shù)值較小，表明林地以轉(zhuǎn)入為主，且轉(zhuǎn)入轉(zhuǎn)出較為頻繁，這與林地轉(zhuǎn)移狀況的研究結(jié)果一致。

表3 土地利用轉(zhuǎn)移矩陣

表4 林地生態(tài)系統(tǒng)服務當量因子及林地單位面積生態(tài)系統(tǒng)服務價值

3.2林地生態(tài)系統(tǒng)服務靜態(tài)價值變化

通過對林地生態(tài)系統(tǒng)服務當量因子系數(shù)的調(diào)節(jié)，結(jié)合謝高地等[2]的生態(tài)系統(tǒng)服務當量因子基礎值確定研究區(qū)林地生態(tài)系統(tǒng)服務當量因子及林地單位面積生態(tài)系統(tǒng)服務價值(表4)。

綜合研究區(qū)土地利用變化(表3)及林地單位面積生態(tài)系統(tǒng)服務價值研究(表4)，分析2008年和2016年研究區(qū)生態(tài)系統(tǒng)服務靜態(tài)價值的變化狀況(表5)。2008年和2016年研究區(qū)生態(tài)系統(tǒng)服務總靜態(tài)價值分別為25.98億元和26.54億元，2016年較2008年增長了0.56億元，增長了2.16%，平均增幅為0.27%，總價值增長幅度不大。2016年林地生態(tài)系統(tǒng)服務靜態(tài)價值為14.97億元，較2008年的13.44億元增長了1.53億元。2016年林地在研究區(qū)生態(tài)系統(tǒng)服務靜態(tài)服務總價值比例為56.41%，較2008年的51.73%提高了4.68%。2個年份林地在研究區(qū)生態(tài)系統(tǒng)服務靜態(tài)價值的比重都在50%以上，遠高于草地及其它地類(表5)。林地在區(qū)域生態(tài)系統(tǒng)服務靜態(tài)價值的比重大且增幅明顯，表明林地是干熱河谷生態(tài)安全服務價值的關(guān)鍵類型，且其作用力不斷增強，主要原因在于林地自身的單位面積生態(tài)系統(tǒng)服務靜態(tài)價值高，林地面積的增長使其生態(tài)系統(tǒng)服務靜態(tài)價值增長更為明顯，在區(qū)域生態(tài)系統(tǒng)服務價值的比重也提升較快。自20世紀50年代以來，人類大幅砍伐森林用于薪材及其它生活和生產(chǎn)活動，導致20世紀80年代研究區(qū)的森林覆蓋率不足10%，森林的生態(tài)功能遭到毀滅性的破壞，極大的危害了長江流域的水土安全，引起了國家和各級政府的高度關(guān)注。作為長江上游水土保持重點區(qū)域之一，元謀強化了天然林工程、退耕還林還草工程及其它植被恢復措施的實施，林地比例的不斷增加促進了元謀干熱河谷生態(tài)環(huán)境的恢復。

表5 不同年份土地利用類型的生態(tài)系統(tǒng)服務靜態(tài)價值

林地9個單項生態(tài)系統(tǒng)服務功能按價值量的大小排序依次為:土壤形成與保護>氣體調(diào)節(jié)>生物多樣性保護>水源涵養(yǎng)>氣候調(diào)節(jié)>原材料>廢物處理>娛樂文化>食物生產(chǎn)(表6)。林地在土壤形成與保護、氣體調(diào)節(jié)和生物多樣性保護的作用最大，食物生產(chǎn)價值最弱。由于研究區(qū)林地以灌木林為主，灌木林是研究區(qū)的景觀基質(zhì)，較好的適應了當?shù)氐母蔁釟夂颉Ａ帜靖翟谕寥乐械难由鞛橥寥赖男纬珊捅Ｗo提供了豐富的養(yǎng)分和水分來源，是固土和水土保持的關(guān)鍵因素，因而林地的土壤形成與保護單向價值高。林木的光合作用在固定CO2和釋放O2方面具有重要作用，故氣體調(diào)節(jié)功能價值高。林地植被覆蓋度高，生態(tài)環(huán)境質(zhì)量好，為研究區(qū)動植物生存提供了良好的生境，因而生物多樣性保護價值高。灌木林與有林地相比生物量偏低，食物產(chǎn)出少，以灌木林為主的林地結(jié)構(gòu)導致林地食物生產(chǎn)價值弱。各單項生態(tài)系統(tǒng)服務價值雖有一定差異，但林地總體具有較高的生態(tài)系統(tǒng)服務價值。

表6 不同年份林地的單項生態(tài)系統(tǒng)服務靜態(tài)價值

3.3林地生態(tài)系統(tǒng)服務動態(tài)價值變化

2008年和2016年生態(tài)系統(tǒng)服務動態(tài)價值調(diào)整系數(shù)值分別為0.354和0.603(表7)。2008年和2016年研究區(qū)林地生態(tài)系統(tǒng)服務動態(tài)價值分別為4.758億元和9.027億元，僅為同期林地生態(tài)系統(tǒng)服務靜態(tài)價值的35.40%和60.30%，這主要是因為研究區(qū)經(jīng)濟發(fā)展水平較落后，人均GDP與全國平均水平的差異較大，恩格爾系數(shù)高，因此人們對林地生態(tài)系統(tǒng)服務價值的支付能力和支付意愿較低，導致林地生態(tài)系統(tǒng)服務動態(tài)價值總體偏低。2008年和2016年研究區(qū)生態(tài)系統(tǒng)服務總動態(tài)價值分別為9.197億元和16.004億元，林地生態(tài)系統(tǒng)服務動態(tài)價值占研究區(qū)生態(tài)系統(tǒng)服務總動態(tài)價值的比例分別為51.73%和56.41%，表明林地是區(qū)域生態(tài)系統(tǒng)服務價值的關(guān)鍵類型。

表7 林地生態(tài)系統(tǒng)服務動態(tài)價值調(diào)整系數(shù)

從林地生態(tài)系統(tǒng)服務動態(tài)價值變化狀況看，8年間其值增長了4.269億元,增幅達89.72%，遠超過同期林地生態(tài)系統(tǒng)靜態(tài)服務價值增長率(11.38%)。這是由于較2008年相比，2016年研究區(qū)人均GDP增長非常迅速,與全國平均GDP水平的差異縮小了25.6%，支付能力提高明顯。社會和經(jīng)濟的迅速發(fā)展較大的提升了社會發(fā)展階段，恩格爾系數(shù)降低，逐漸接近0.40的拐點，支付意愿得到提升，因而生態(tài)系統(tǒng)服務動態(tài)價值顯著提升。隨著支付能力和支付意愿的迅速提升，研究區(qū)林地生態(tài)系統(tǒng)服務動態(tài)價值將會顯著增加。林地生態(tài)系統(tǒng)服務動態(tài)價值增長速度(89.72%)高于區(qū)域生態(tài)系統(tǒng)服務動態(tài)價值的增長均值(74.01%)，表明高生態(tài)價值的土地利用類型對社會和環(huán)境調(diào)節(jié)能力增長更為敏感。

4 討論

Costanza等[1]和謝高地等[2]的研究為基于土地利用類型的生態(tài)服務價值評估奠定了基礎。在不同地域和類型的林地生態(tài)系統(tǒng)服務價值研究中，學者們多結(jié)合研究區(qū)的自然和植被特征對生態(tài)系統(tǒng)服務當量因子進行適當?shù)恼{(diào)整，李曉賽等[35]的生態(tài)系統(tǒng)服務當量因子調(diào)整方法較好的體現(xiàn)了林地類型間生物量差異對林地生態(tài)系統(tǒng)服務價值的影響，因此，本研究借鑒了李曉賽等[35]的研究方法，調(diào)整了研究區(qū)的林地生態(tài)系統(tǒng)服務當量因子系數(shù)。研究表明，林地對研究區(qū)水土、氣候和生物保護起關(guān)鍵作用，是提供區(qū)域生態(tài)系統(tǒng)服務功能的主要類型，生態(tài)系統(tǒng)服務價值高，這與干熱河谷生態(tài)系統(tǒng)實際狀況是相符的；然而，干熱河谷林地植被類型較復雜，且不同的林地植被演替階段的差異亦會導致生物量和生態(tài)系統(tǒng)服務價值的差別，但此方面的研究卻處于空白階段，因此，難以完全真實的反映干熱河谷林地生態(tài)系統(tǒng)的功能和價值。為提高干熱河谷林木生態(tài)系統(tǒng)服務價值研究的精確性，應結(jié)合實際調(diào)查更細致的測量不同植被類型和生長階段的林地生態(tài)系統(tǒng)單向服務價值。在生態(tài)系統(tǒng)動態(tài)價值系數(shù)研究中，恩格爾系數(shù)常作為反映支付意愿的指標。雖然恩格爾系數(shù)和支付意愿之間存在一定的相關(guān)關(guān)系，但二者之間是否存在明確的定量相關(guān)性卻沒有得到有效研究，而支付意愿同時也是居民的社會屬性和心理特征的反映，應開展社會發(fā)展階段、社會屬性特征及心理特征與支付意愿的量化關(guān)系研究以深入探討支付意愿對生態(tài)系統(tǒng)服務動態(tài)價值的影響。

5 結(jié)論

(1)中高山地帶是研究區(qū)林地的主要分布區(qū)，人為干擾使林地空間分布更破碎。林地主要分布在南部和東部中高山，在西部、西南部和北部、東部邊緣中高山地帶有少許分布，河谷壩區(qū)也有零星分布。退耕還林等生態(tài)保護工程促進了林地的生態(tài)恢復，而人為干擾的作用(主要是人工林的建設)加深了林地的破碎化程度。

(2)林地面積有明顯增長，林地與非林地之間相互轉(zhuǎn)化頻繁。8年來研究區(qū)的林地面積，增長了64.24 km2，變化幅度達到11.42%，年變化率達1.42%，總體幅度變化較大。林地以轉(zhuǎn)入為主，但轉(zhuǎn)入轉(zhuǎn)出較為頻繁，林地與草地的相互轉(zhuǎn)化尤為明顯，土地利用類型的變化對林地的干擾作用加強。

(3)林地是提供干熱河谷生態(tài)系統(tǒng)服務價值的關(guān)鍵類型。林地占生態(tài)系統(tǒng)服務的靜態(tài)價值，在區(qū)域生態(tài)系統(tǒng)服務總靜態(tài)價值中的比重在50%以上，且比重不斷提高，對區(qū)域生態(tài)系統(tǒng)服務的作用力不斷增強。林地的土壤形成與保護、氣體調(diào)節(jié)和生物多樣性保護單向生態(tài)系統(tǒng)服務價值最高，食物生產(chǎn)價值最弱。

(4)林地生態(tài)系統(tǒng)服務靜態(tài)價值和動態(tài)價值差異性大?，F(xiàn)階段人們對林地生態(tài)系統(tǒng)服務價值的支付能力和支付意愿較低，使林地生態(tài)系統(tǒng)服務動態(tài)價值總體偏低。研究區(qū)林地生態(tài)系統(tǒng)服務動態(tài)價值遠低于同期靜態(tài)價值，但其變化幅度遠高于靜態(tài)價值，反映人們支付能力和生態(tài)環(huán)境意識的提升。具有高生態(tài)系統(tǒng)服務價值的林地對社會和環(huán)境調(diào)節(jié)能力的增長更敏感。

[1] Costanza R, d’Arge R, de Groot R S,etal. The value of the world's ecosystem services and natural capital[J]. Nature, 1997, 387: 253-260.

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

[3] Mcneely J A, Miller K R, Reid W V,etal. Conserving the world’s biological diversity[M]. Gland：International Union for Conservation of Nature and Natural Resources, 1990．

[4] Pearce D W. Blueprint 4: Capturing global environmental value[M].London: Earthscan, 1995.

[5] Sutton P C, Costanza R. Global estimates of market and non-market values derived from night time satellite imagery, land cover, and ecosystem service valuation[J]. Ecological Economics, 2002, 41(3): 509-527.

[6] Potschin M B, Hainesyoung R H . Ecosystem services[J].Progress in Physical Geography, 2011, 35(5): 575-594.

[7] Lal P. Economic valuation of mangroves and decision-making in the in the Pacific[J].Ocean & Coastal Management, 2003,46(9-10): 823-844.

[8] Kroeger T, Casey F. An assessment of market-based approaches to providing ecosystem services on agricultural lands[J]．Ecological Economics, 2007, 64(2): 321-332.

[9] Jenkins W A, Murray B C,Kramer R A,etal. Valuing ecosystem services from wetlands restoration in the Mississippi Alluvial valley[J]. Ecological conomics, 2010, 69(5): 1051-1061.

[10] Chee Y. An ecological perspective on the valuation of ecosystem services[J]. Biological Conservation, 2004, 120: 549-565.

[11] Mendoza-González G, Martínez M L, Lithgow D,etal.Land use change and its effects on the value of ecosystem services along the coast of the Gulf of Mexico[J].Ecological Economics, 2012, 82(20): 23-32.

[12] Toman M. Why not to calculate the value of the world’s ecosystem services and natural capital[J]. Ecological Economics, 1998, 25(1): 57-60.

[13] Heal G. Valuing ecosystem services[J]. Ecological Economics, 2000, 35(3): 24-30.

[14] Wilson M A, Howarth R. Discourse based valuation of ecosystem services: establishing fair outcomes through group deliberation[J]. Ecological Economics, 2002, 41(3): 431-443.

[15] Mathias J D, Bonté B,Cordonnier T,etal. Using the viability theory to assess the flexibility of forest managers under ecological intensification[J].Environmental Management, 2015, 56(5):1170-1183.

[16] Vedel S E, Jacobsen J B,Bo J T .Forest owners’willingness to accept contracts for ecosystem service provision is sensitive to additionality. Ecological Economics, 2015, 113: 15-24.

[17] Dawson N, Martin A. Assessing the contribution of ecosystem services to human wellbeing: A disaggregated study in western Rwanda[J]. Ecological Economics, 2015, 117: 62-72.

[18] Foelker C J,Christopher Standley C R, Parry D,etal. Complex ecological relationships among an assemblage of indigenous hymenopterans parasitoids, the exotic European wood wasp (Sirex noctilio; Hymenoptera: Siricidae), and a native congener[J].Canadian Entomologist, 2016, 148(5): 532-542.

[19] Vibrans A C, Moser P, Oliveira L Z,etal.Generic and specific stem volume models for three subtropical forest types in southern Brazil[J].Annals of Forest Science, 2015, 72(6): 865-874.

[20] Matsushita K, Tomotsune M, Sakamaki Y,etal. Effects of management treatments on the carbon cycle of a cool-temperate broad-leaved deciduous forest and its potential as a bioenergy source[J].Ecological Research, 2015, 30(2): 293-302

[21] 王 麗.青城山高保護價值森林的生態(tài)服務功能評價研究[D].成都：四川農(nóng)業(yè)大學, 2012.

[22] 鄭江坤, 余新曉, 夏 兵, 等.潮白河流域林地轉(zhuǎn)化及森林生態(tài)服務價值動態(tài)分析[J].農(nóng)業(yè)工程學報, 2010, 26 (S1):308-314.

[23] 馬國軍, 林 棟. 石羊河流域生態(tài)系統(tǒng)服務功能經(jīng)濟價值評估[J].中國沙漠,2009, 29(6):1173 -1177.

[24] 汪有奎,郭生祥,汪 杰,等.甘肅祁連山國家級自然保護區(qū)森林生態(tài)系統(tǒng)服務價值評估[J].中國沙漠, 2013, 33(6): 1905-1911.

[25] 龍會英,沙毓滄,朱紅業(yè),等. 8份圭亞那柱花草在元謀干熱河谷的引種研究[J].西南農(nóng)業(yè)大學學報, 2007, 20(5): 1078-1083.

[26] 馬姜明,李 昆,張昌順. 元謀干熱河谷蘇門答臘金合歡、新銀合歡人工林天然更新初步研究[J].應用生態(tài)學報, 2006, 17(8): 1365-1369.

[27] 方海東, 段昌群, 紀中華,等. 金沙江干熱河谷自然恢復區(qū)植物種群生態(tài)位特征[J]. 武漢大學學報：理學版, 2008, 54(2): 177-182.

[28] 高文學,王志和,周慶宏，等.金沙江干熱河谷稀樹灌草叢植被恢復方式研究[J].林業(yè)調(diào)查規(guī)劃, 2005, 30(3): 87-91.

[29] 李 彬, 唐國勇,李 昆，等.元謀干熱河谷20年生人工恢復植被生物量分配與空間結(jié)構(gòu)特征[J]. 應用生態(tài)學報, 2013, 24(6): 1479-1486.

[30] 周紅藝,熊東紅,楊 忠.元謀干熱河谷土地利用變化對生態(tài)系統(tǒng)服務價值的影響[J].山地學報, 2008, 24(3): 135-138.

[31] 楊 忠，熊東紅，周紅藝，等．干熱河谷不同巖土組成坡地的降水入滲與林木生長[J].中國科學E輯，2003，33(增刊)：110-119.

[32] 歐曉昆. 元謀干熱河谷的自然生態(tài)特點及開發(fā)利用意見[J].西部林業(yè)科學, 1987(1): 17-19.

[33] 何錦峰, 蘇春江, 舒 蘭，等.基于3S技術(shù)的金沙江干熱河谷區(qū)LUCC研究——以云南省元謀縣為例[J].山地學報, 2009, 27(3): 341-348.

[34] 孫長安.香溪河流域土地利用與水土流失的關(guān)系研究[D].北京: 北京林業(yè)大學, 2008.

[35] 李曉賽, 朱永明, 趙 麗, 等.基于價值系數(shù)動態(tài)調(diào)整的青龍縣生態(tài)系統(tǒng)服務價值變化研究[J].中國生態(tài)農(nóng)業(yè)學報, 2015,23(3): 373-381.

[36] 方精云,柯金虎, 唐志堯,等.生物生產(chǎn)力的“4P”概念、估算及其相互關(guān)系[J].植物生態(tài)學報, 2001, 25(4): 414-419.

[37] 宗躍光,陳紅春,郭瑞華,等.地域生態(tài)系統(tǒng)服務功能的價值結(jié)構(gòu)分析——以寧夏靈武市為例[J]. 地理研究, 2000, 19(2): 148-155.

[38] 粟曉玲, 康邵忠, 佟 玲.內(nèi)陸河流域生態(tài)系統(tǒng)服務價值的動態(tài)估算方法[J].生態(tài)學報, 2006, 26(6): 2011-2019.

[39] 杜金龍. 土地利用變化及其對生態(tài)系統(tǒng)服務價值影響研究[D]. 武漢: 華中農(nóng)業(yè)大學, 2010.

[40] 鄧舒洪. 區(qū)域土地利用變化與生態(tài)系統(tǒng)服務價值動態(tài)變化研究[D]. 杭州:浙江大學, 2012.

(責任編輯：詹春梅)

ForestEcosystemServiceValueChangeinYuanmouDry-hotValley

OUZhao-rong1,2，ZHUQing-ke1,SUNYong-yu3

(1.School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China; 2.School of Ecotourism, Southwest Forestry University, Kunming 650224, Yunnan, China; 3.Research Institute of Resources Insects, Chinese Academy of Forestry, Kunming 650216, Yunnan, China)

S718.55

A

1001-1498(2017)05-0831-10

10.13275/j.cnki.lykxyj.2017.05.017

2017-05-25

國家重點研發(fā)計劃(2017YFC0505102)；中國林科院中央級公益性科研院所基本科研業(yè)務費專項資金資助項目(CAFYBB2017ZA002-4)

歐朝蓉(1978—)，女，湖南衡陽人，博士研究生，講師.主要研究方向：生態(tài)環(huán)境地理學.E-mail: flymoon97@163.com

* 通訊作者：孫永玉(1974—)，男，山東濟寧人，博士，副研究員.主要研究方向為干熱河谷恢復生態(tài)學.E-mail:cafsdren@163.com