李佳林 廖家紅 王旭熙 劉守江 楊丹 鄢繼堯

摘要：根據(jù)嘉陵江流域2000、2005和2015年3期土地利用遙感影像數(shù)據(jù)，將地理數(shù)學方法與GIS相結(jié)合，研究了嘉陵江流域15年各種土地利用類型的面積變化、土地利用動態(tài)度等動態(tài)變化信息，并估算了各種土地利用類型的生態(tài)系統(tǒng)服務價值。結(jié)果表明，2000—2015年，嘉陵江流域耕地面積持續(xù)減少，而建設(shè)用地、草地、水域、林地和未利用地面積增加;嘉陵江流域生態(tài)系統(tǒng)服務價值由2000年的243 109.36百萬元減少到2015年的242 805.95百萬元;嘉陵江流域生態(tài)系統(tǒng)服務價值對價值系數(shù)缺乏彈性，反映耕地面積的變化對研究區(qū)生態(tài)系統(tǒng)服務價值的變化具有較大影響。

關(guān)鍵詞：土地利用;生態(tài)系統(tǒng)服務價值;嘉陵江流域;評估

中圖分類號：F301? ? ? ? ?文獻標識碼：A

文章編號：0439-8114（2019）02-0033-05

DOI：10.14088/j.cnki.issn0439-8114.2019.02.008? ? ? ? ? ?開放科學（資源服務）標識碼（OSID）：

生態(tài)系統(tǒng)服務功能是指生態(tài)系統(tǒng)服務與生態(tài)系統(tǒng)功能的綜合，指生態(tài)系統(tǒng)與生態(tài)過程所形成及維持的人類賴以生存的自然環(huán)境條件與效用[1]。當前有關(guān)生態(tài)系統(tǒng)服務價值（Ecosystem service value，ESV）的研究已成為生態(tài)經(jīng)濟學和環(huán)境經(jīng)濟學的研究熱點[2]。Costanza等[3]于1997年在《Nature》上發(fā)表了論文“全球生態(tài)系統(tǒng)服務價值和自然資本”，掀起了生態(tài)系統(tǒng)服務研究的熱潮，在全球、國家、市、縣、區(qū)域、流域等不同空間尺度上，針對森林、草地、農(nóng)田、濕地、海洋等不同生態(tài)系統(tǒng)類型的研究廣泛開展，把生態(tài)系統(tǒng)提供的商品和服務統(tǒng)稱為生態(tài)系統(tǒng)服務，同時對生態(tài)系統(tǒng)服務進行了分類，并對全球生態(tài)系統(tǒng)與自然資本的價值進行了估算。在中國實際情況下，謝高地等[4]綜合國內(nèi)數(shù)百位具有生態(tài)學背景的專業(yè)人員的調(diào)查問卷結(jié)果，制定了中國生態(tài)系統(tǒng)單位面積生態(tài)服務價值當量。

土地利用/覆蓋變化是人類活動對環(huán)境影響最直接的表現(xiàn)，土地利用現(xiàn)狀圖也是人地關(guān)系的一面鏡子[5]。土地利用/覆蓋變化通過改變地表覆蓋類型，從而影響生態(tài)系統(tǒng)的結(jié)構(gòu)、過程、功能，進而影響區(qū)域生態(tài)系統(tǒng)服務價值[6]。諸多學者也在不同空間尺度上就土地利用變化對生態(tài)系統(tǒng)服務價值的影響進行了相關(guān)研究[7，8]。對于快速城市化區(qū)域，由于人口迅速增長、城市迅速擴張、社會經(jīng)濟高速發(fā)展，生態(tài)系統(tǒng)服務價值對土地利用劇烈變化呈現(xiàn)出明顯和復雜的時空變化特征[9]。

目前對于嘉陵江流域土地利用和生態(tài)系統(tǒng)服務價值的研究并不豐富，主要集中于應用GIS空間分析方法分析嘉陵江流域土地利用/覆被變化、影響因素及驅(qū)動力[10]，而對嘉陵江流域整體生態(tài)系統(tǒng)服務價值變化并沒有詳細分析。本研究以嘉陵江流域2000、2005和2015年土地利用變化為基礎(chǔ)，借鑒Costanza等[3]、謝高地等[4]的研究方法，評價了嘉陵江流域2000—2015年生態(tài)系統(tǒng)服務價值對土地利用變化的響應，并分析了該地區(qū)生態(tài)系統(tǒng)服務價值的變化及其產(chǎn)生原因，旨在為嘉陵江流域土地資源可持續(xù)利用和生態(tài)環(huán)境保護提供科學依據(jù)。

1? 研究區(qū)概況

嘉陵江（29°20′N-34°30′N，102°30′E-109°51′E）發(fā)源于陜西省秦嶺南麓，由北向南流經(jīng)陜西、甘肅、四川和重慶三省一市，其干流分東西兩源，東源起自陜西省鳳縣以北的秦嶺鎮(zhèn)，西源起自甘肅省天水平南川，習慣上東源為正源，西源稱為西漢水。東源向南流經(jīng)甘肅省徽縣，于陜西略陽與西源西漢水匯合，經(jīng)陽平關(guān)入川，然后向南流經(jīng)廣元、昭化、蒼溪、閬中、南部、蓬安、南充、武勝等兩市六縣，在合川入渝，在重慶匯入長江，全長1 119 km。嘉陵江流域?qū)賮啛釒Ъ撅L氣候，流域多年平均降水量為935.2 mm，由于地形復雜，降水時間分布很不均勻，降水量有明顯的旱季和雨季，6—9月降水量占年降水量的66%，12月至次年2月為全年降水量最少的季節(jié)，占年降水量的3.2%。整個流域多為紫紅色沙泥巖，質(zhì)地松脆，水土流失嚴重，成為流域內(nèi)泥沙產(chǎn)生的主要來源。嘉陵江流域自然條件復雜，自然資源豐富，人口眾多，是西部地區(qū)重要的糧、油、果、豬、蠶生產(chǎn)基地[11]。嘉陵江是長江上游重要的支流，聯(lián)系了中國西部的甘、陜、川、渝4個省市，該流域是聯(lián)系“陸上絲綢之路經(jīng)濟帶”和“長江經(jīng)濟帶”的重要紐帶[12]。

鑒于嘉陵江流經(jīng)范圍廣，且歷史流數(shù)據(jù)不易獲取，以縣（市）為基本地域單元，研究2000—2015年嘉陵江流域生態(tài)系統(tǒng)服務價值。基于數(shù)據(jù)收集的可得性，以嘉陵江流域2008年平均糧食單產(chǎn)及其相對應的2008年各作物全國的平均收購價格為基準，由當量因子價格和當量因子得出流域不同生態(tài)系統(tǒng)類型的生態(tài)價值系數(shù)。

2? 數(shù)據(jù)來源與研究方法

2.1? 土地利用動態(tài)數(shù)據(jù)獲取

根據(jù)研究區(qū)2000、2005和2015年3期的土地利用數(shù)據(jù)。該數(shù)據(jù)由陸地衛(wèi)星Land Sat同一時相遙感影像作為基礎(chǔ)數(shù)據(jù)，在ENVI 5.2遙感軟件支持下配準，進行遙感影像預處理、輻射糾正和集合糾正、圖像配準[13];在GIS的支持下，以柵格影像為背景，利用人機交互提取專題數(shù)據(jù)，并參考全國土地利用分類方法，將嘉陵江流域用地類型分為6類，分別是耕地、林地、草地、水域、建設(shè)用地和未利用地。應用ArcGIS 10.2軟件屬性表的空間計算功能，統(tǒng)計嘉陵江流域6類土地利用類型的面積，得到嘉陵江流域2000、2005、2015年土地利用分類。

2.2? 土地利用研究方法

土地利用類型動態(tài)度（K）指某研究區(qū)在一定時間范圍內(nèi)某種土地利用類型的數(shù)量變化情況[14，15]。為了反映嘉陵江流域各類土地的動態(tài)變化情況，采用土地利用類型動態(tài)度模型具體分析該流域土地利用變化情況。計算公式[16，17]如下。

式中，Ua、Ub分別為研究初期及研究末期某一種土地利用類型的數(shù)量;T為時段，當T被設(shè)定為年時，K就是該研究區(qū)某種土地利用類型的年變化率。

2.3? 生態(tài)系統(tǒng)服務價值研究方法

根據(jù)研究區(qū)土地資源特征與景觀變化差異，參照中國《土地利用現(xiàn)狀調(diào)查技術(shù)規(guī)程》和前人采用的土地分類系統(tǒng)[3，17，18]，將研究區(qū)土地利用類型劃分為6個類別，即耕地（水田、水澆地和旱地等）、林地（有林地、灌木林、疏林地和其他林地）、草地、建設(shè)用地（城鄉(xiāng)居民用地與工礦交通用地等）、水域（河湖溝渠、水庫坑塘、灘涂灘地等）、未利用地（鹽堿地、沙地、戈壁、裸巖石質(zhì)地和沼澤地等）。對建設(shè)用地按照Costanza等[3]的方法，不估算其生態(tài)系統(tǒng)服務功能經(jīng)濟價值。

2.3.1? 生態(tài)系統(tǒng)服務價值計算方法? 結(jié)合 Costanza等[3]提出的生態(tài)系統(tǒng)服務價值（ESV）計算公式和謝高地等[4]2007年基于500份調(diào)查問卷修訂后的中國生態(tài)系統(tǒng)服務價值當量因子，估算流域的生態(tài)系統(tǒng)服務價值。根據(jù)嘉陵江流域統(tǒng)計年鑒，為消除農(nóng)作物價格波動對總價值量的影響，以嘉陵江流域2008年平均糧食單產(chǎn)5 431.14 kg/hm2，及其相對應的2008年各作物全國的平均收購價格1.69元/kg為基準，考慮到?jīng)]有人力投入的自然生態(tài)系統(tǒng)提供的經(jīng)濟價值是現(xiàn)有單位面積農(nóng)田提供的食物生產(chǎn)服務經(jīng)濟價值的1/7[19，4]，計算出嘉陵江流域一個當量因子的價格約為1 311.23元/hm2。由當量因子價格和當量因子，得出流域不同生態(tài)系統(tǒng)類型的生態(tài)價值系數(shù)（表1）。主要計算公式[17]如下。

式中，ESV為研究區(qū)生態(tài)系統(tǒng)服務總價值（元）;VCK為第k類土地利用類型單位面積的生態(tài)功能服務價值系數(shù)（元/hm2）;Ak為研究區(qū)內(nèi)第k類土地利用類型的面積（hm2）;n為土地利用類型數(shù)。

2.4? 敏感性分析方法

借用敏感性指數(shù)CS（Sensitivity coefficient）驗證生態(tài)系統(tǒng)類型對各種土地類型的代表性和價值的準確性，同時檢驗謝高地等[4]的生態(tài)系統(tǒng)服務價值系數(shù)對本研究區(qū)的適合程度[20，21]，進而確定生態(tài)系統(tǒng)服務價值隨時間的變化對于生態(tài)服務功能價值系數(shù)（VC）的依賴程度。敏感性指數(shù)CS的含義是指VC變動引起的ESV的變化情況。如果CS>1，表明ESV相對于VC是富有彈性的;如果CS<1，ESV則被認為是缺乏彈性的。敏感性系數(shù)越大，表明生態(tài)系統(tǒng)服務價值系數(shù)的準確性越關(guān)鍵，過高或過低賦值的生態(tài)系統(tǒng)服務價值系數(shù)也很可能在很大程度上影響生態(tài)系統(tǒng)服務價值隨時間變化的真實性。敏感性指數(shù)的計算公式如下。

式中，ESV為總生態(tài)系統(tǒng)服務價值;VC為生態(tài)系統(tǒng)價值系數(shù);i和j分別代表初始價值和生態(tài)系統(tǒng)價值系數(shù)調(diào)整后價值;k為各土地利用類型。

3? 結(jié)果與分析

3.1? 土地利用變化分析

2000—2015年嘉陵江流域土地利用情況發(fā)生了較明顯變化，2000—2015年嘉陵江流域耕地集中分布在南部，區(qū)域集約化程度高，北部東部有少許耕地，主要原因是山區(qū)地形條件限制;水域和草地在這15年間呈增加趨勢;林地和未利用地和水域變化不大;建設(shè)用地集中在中部、南部，2015年中部、東部建設(shè)用地增加，區(qū)域集約化程度提高。

由表2可知，2000—2015年嘉陵江流域土地利用類型發(fā)生了較大變化，其總體趨勢為耕地減少，林地、草地、建設(shè)用地、水域和未利用地面積增加。2000—2015年耕地面積大幅度持續(xù)減少，減少量達1 620 km2，變化率為-2.06%，且2005—2015年的變化率大于2000-2005年的變化率;林地面積在2000—2015年先增加后持平，整體增加了88 km2;2000—2015年草地面積增加了190 km2，變化率為0.36%;2000—2005年建設(shè)用地面積增加，變化量達227 km2，變化率為15.77%，2005—2015年持續(xù)增加，增加了956 km2，變化率為 57.38%，因此，2000—2015年建設(shè)用地面積整體呈上升趨勢，增加了1 183 km2，變化率為82.21%;2000年水域面積僅占總用地面積的0.91%，2000—2005年水域面積增加了19 km2，2005—2015年增加了77 km2，2000—2015年水域面積整體呈增加趨勢;2000—2015年未利用地變化也較為明顯，總體上呈增加趨勢，增加了63 km2，變化率為2.48%。

嘉陵江流域在2000—2015年建設(shè)用地變化量最大，變化速度最快，動態(tài)度為5.48%;其次是水域、未利用地、耕地，它們的動態(tài)度分別為0.36%、0.17%、-0.14%;林地和草地的總體變化較平緩，動態(tài)度分別為0.01%、0.02%（表3）。

建設(shè)用地面積的大幅度增加主要與流域內(nèi)城市規(guī)劃逐漸成熟以及后期旅游業(yè)發(fā)展迅速有關(guān)[22]，城市建設(shè)速度加快，交通用地和商業(yè)用地不斷擴大，在此進程中，同時增加了對耕地的需求與開發(fā)，導致其面積明顯減少。

3.2? 土地利用變化對生態(tài)系統(tǒng)服務價值的影響

根據(jù)嘉陵江流域不同陸地生態(tài)系統(tǒng)服務價值（表1）和嘉陵江流域各土地利用類型面積（表2），計算得到嘉陵江流域生態(tài)系統(tǒng)服務價值的變化情況（表4）。從表4可以看出，嘉陵江流域2000—2015年生態(tài)系統(tǒng)服務總價值呈現(xiàn)減少的變化趨勢，2000、2005和2015年的生態(tài)系統(tǒng)服務總價值分別為243 109.36百萬、243 236.49百萬和242 805.95百萬元，2000—2005年增加了127.13百萬元。在研究期間，耕地大面積減少，2000—2005年和2005—2015年其生態(tài)系統(tǒng)服務價值損失量分別達到了365.60百萬和693.92百萬元;草地生態(tài)系統(tǒng)服務價值先增加后減少，2005—2015年面積減少，其生態(tài)系統(tǒng)服務價值損失量達到了32.85百萬元，而2000—2015年整體增加了183.57百萬元;未利用地面積變化不明顯，其生態(tài)系統(tǒng)服務價值增加了7.25百萬元;2000—2015年水域面積擴大，水域生態(tài)系統(tǒng)服務價值增加了360.43百萬元;林地價值先增加后無變化，整體增加了204.86百萬元。嘉陵江流域在2000—2015年耕地向水域、草地、林地、未利用地和建設(shè)用地轉(zhuǎn)化，雖然水域、草地、林地所提供服務功能的經(jīng)濟價值高于耕地，但因轉(zhuǎn)化為水域、草地、林地的面積遠小于耕地的減少面積，大部分耕地轉(zhuǎn)化為建設(shè)用地，因而整體生態(tài)系統(tǒng)服務價值減少。

從生態(tài)系統(tǒng)服務價值構(gòu)成（表5）來看，研究時段內(nèi)林地所產(chǎn)生的生態(tài)服務價值所占比重最大，2000、2005和2015年分別為54.98%、55.03%、55.12%，占整個生態(tài)系統(tǒng)的主導地位;耕地和草地的生態(tài)系統(tǒng)服務價值也較高，所占比例分別為21.20%、21.04%、20.79%和20.98%、21.05%、21.08%，二者生態(tài)系統(tǒng)服務價值變化小;水域的生態(tài)系統(tǒng)服務價值相對較小，在嘉陵江流域生態(tài)系統(tǒng)服務總價值中處于第四位;未利用地生態(tài)系統(tǒng)服務價值的貢獻率非常低，在研究時段貢獻率僅0.12%。

3.3? 生態(tài)服務價值敏感度分析

根據(jù)敏感性指數(shù)的計算公式，把VC分別上下調(diào)整50%，計算出研究區(qū)2000、2005和2015年的敏感性。結(jié)果（表6）表明，流域所有土地利用類型生態(tài)系統(tǒng)服務價值系數(shù)的敏感性指數(shù)（CS）均小于1，說明生態(tài)系統(tǒng)服務價值對生態(tài)價值系數(shù)缺乏彈性，研究結(jié)果可信。未利用地的CS最小，為0.001 2，即當未利用地的生態(tài)價值系數(shù)增加1%時，流域總的生態(tài)系統(tǒng)服務價值僅增加0.001 2%。林地因面積較大，其在流域總生態(tài)系統(tǒng)服務價值中的貢獻最大，敏感性指數(shù)CS也最大，為0.549 7～0.551 2;草地和耕地的價值系數(shù)較大，CS次之;未利用地面積小于林地，其價值系數(shù)較低，CS較小。研究期間各類型土地的CS變化趨勢與各自面積的變化趨勢基本一致，年際變化不明顯。

4? 小結(jié)與討論

本研究將地理數(shù)學方法與GIS相結(jié)合，對嘉陵江流域2000—2015年的土地利用變化情況和生態(tài)系統(tǒng)服務價值進行分析，為嘉陵江流域土地資源可持續(xù)利用和生態(tài)環(huán)境保護提供科學依據(jù)。通過研究主要得出以下結(jié)論。

1）2000—2015年嘉陵江流域土地利用類型發(fā)生了較大變化，耕地、林地和草地是嘉陵江流域的主要土地類型，研究期內(nèi)各種土地利用類型變化較明顯。其總體趨勢為耕地減少，減少量達1 620 km2，變化率為-2.06%，而林地、草地、建設(shè)用地、水域和未利用地面積增加。

2）嘉陵江流域生態(tài)系統(tǒng)服務價值由2000年的243 109.36百萬元減少到2015年的242 805.95百萬元，說明嘉陵江流域生態(tài)系統(tǒng)服務價值呈下降趨勢，耕地向水域、草地、林地、未利用地和建設(shè)用地轉(zhuǎn)化，雖然水域、草地、林地所提供服務功能的經(jīng)濟價值高于耕地，但因轉(zhuǎn)化為水域、草地、林地的面積遠小于耕地的減少面積，大部分耕地轉(zhuǎn)化為建設(shè)用地，因而整體生態(tài)系統(tǒng)服務價值減少，生態(tài)系統(tǒng)服務價值對土地利用變化表現(xiàn)出了較強的敏感性。

3）嘉陵江流域生態(tài)系統(tǒng)服務價值相對于價值系數(shù)的敏感性都小于1，說明前者對后者是缺乏彈性的，本研究所選的生態(tài)價值系數(shù)是較合理的，結(jié)果可信。

4）隨著嘉陵江流域經(jīng)濟快速發(fā)展，城市化建設(shè)將不可避免，但流域生態(tài)系統(tǒng)服務價值下降。因此，協(xié)調(diào)社會經(jīng)濟發(fā)展與生態(tài)環(huán)境保護之間的矛盾，制定合理的土地利用規(guī)劃和決策，是恢復和保育嘉陵江流域生態(tài)系統(tǒng)服務的有效途徑。

參考文獻：

[1] DAILY G. Natures Services：Societal Dependence on Natural Ecosystems[M].Washington D C：Island Press，1997.

[2] CHEN Z X，ZHANG X S. Value of ecosystem services in China[J].Chinese Science Bulletin，2000，45（10）：870-876.

[3] COSTANZA R，DARGE R，GROOT R DE，et al. The value of the worlds ecosystem services and natural capital[J].Nature，1997，387：253-260.

[4] 謝高地，甄? 霖，魯春霞，等.一個基于專家知識的生態(tài)系統(tǒng)服務價值化方法[J].自然資源學報，2008，23（5）：911-919.

[5] TUAN Y F. Geography，phenomenology，and the study of human nature[J].Canadian Geographer，2010，15（3）：181-192.

[6] 劉紀遠，劉明亮，莊大方，等.中國近期土地利用變化的空間格局分析[J].中國科學：地球科學（D輯），2002，32（12）：1031-1040.

[7] YUE S P，ZHANG S W，YAN Y C. Impacts of land use change on ecosystem services value in the northeast China transect （NECT）[J].Acta Geographica Sinica，2007，62（8）：879-886.

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

[9] ZHANG Z，WU C F，TAN R. Application of ecosystem service value in land use change research：Bottlenecks and prospects[J].Chinese Journal of Applied Ecology，2013，24（2）：556-562.

[10] 齊師杰，張行南，夏達忠，等.嘉陵江流域土地利用/覆被變化特征及其驅(qū)動力分析[J].長江科學院院報，2013，30（1）：1-7.

[11] 翟有龍.嘉陵江流域農(nóng)業(yè)產(chǎn)業(yè)化經(jīng)營模式的選擇[J].經(jīng)濟地理，2004，24（3）：387-390.

[12] 何博汶，楊顯明.嘉陵江流域（川渝地區(qū)）城鎮(zhèn)經(jīng)濟聯(lián)系空間格局演化研究[J].四川文理學院學報，2017，27（4）：113-118.

[13] 唐古拉，周萬村，馮文蘭，等.基于遙感影像多尺度分析技術(shù)的LUCC研究[J].遙感技術(shù)與應用，2005，20（6）：581-585.

[14] 陳百明，劉新衛(wèi)，楊? 紅.LUCC研究的最新進展評述[J].地理科學進展，2003，22（1）：22-29.

[15] 王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討[J].地理科學進展，1999，18（1）：81-87.

[16] 張燦榮，伍素貞.利用遙感進行土地利用動態(tài)變化監(jiān)測——以東莞市長安鎮(zhèn)為例[J].地礦測繪，2007，23（3）：17-19.

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

[18] 劉紀遠.中國資源環(huán)境遙感宏觀調(diào)查與動態(tài)研究[M].北京：中國科學技術(shù)出版社，1996.

[19] 肖? 玉，謝高地，安? 凱.莽措湖流域生態(tài)系統(tǒng)服務功能經(jīng)濟價值變化研究[J].應用生態(tài)學報，2003，14（5）：676-680.

[20] KREUTER U P，HARRIS H G，MATLOCK M D，et al. Change in ecosystem service values in the San Antonio area，Texas[J].Ecological Economics，2001，39（3）：333-346.

[21] 王宗明，張? 柏，張樹清.吉林省生態(tài)系統(tǒng)服務價值變化研究[J].自然資源學報，2004，19（1）：55-61.

[22] 徐榮悅.嘉陵江經(jīng)濟走廊全域旅游發(fā)展研究[D].四川南充：西華師范大學，2017.