張瑜, 趙曉麗, 左麗君, 張增祥, 徐進(jìn)勇

(1.中國科學(xué)院遙感與數(shù)字地球研究所,北京 100101; 2.中國科學(xué)院大學(xué),北京 100049)

0 引言

黃土高原自然地理條件復(fù)雜，在自然因素和人類活動的綜合作用下，水資源短缺、水土流失及土地荒漠化等生態(tài)問題逐步加劇[1]。為保護(hù)黃土高原的生態(tài)環(huán)境，我國實施了一系列生態(tài)工程，改變了區(qū)域土地利用方式[2-3]，導(dǎo)致生態(tài)環(huán)境發(fā)生了巨大變化[4]。生態(tài)系統(tǒng)的變化反過來又會影響人類社會的發(fā)展，這種反饋在生態(tài)環(huán)境脆弱區(qū)和社會發(fā)展水平較高區(qū)十分顯著，如姚小英等[5]基于1971—2010年間甘肅玉米主產(chǎn)區(qū)降水及產(chǎn)量進(jìn)行研究，結(jié)果表明甘肅地區(qū)處于不同等級的干旱風(fēng)險狀態(tài)，隴東黃土高原地區(qū)玉米種植受干旱影響最大。不合理的土地利用方式造成的嚴(yán)重后果又將敦促人類對土地利用進(jìn)行管理，使其朝合理的方向改變[6]。在研究期內(nèi)，如果生態(tài)系統(tǒng)服務(wù)價值下降，表明研究區(qū)的土地利用方式不合理，應(yīng)對土地利用方式進(jìn)行調(diào)整以提升區(qū)域的生態(tài)環(huán)境質(zhì)量； 相反的，合理的土地利用方式會有利于區(qū)域生態(tài)環(huán)境質(zhì)量的提升，進(jìn)而促進(jìn)區(qū)域的可持續(xù)發(fā)展。生態(tài)系統(tǒng)服務(wù)功能和土地利用是一個矛盾體，彼此相互聯(lián)系、相互制約[7]。因此，研究生態(tài)系統(tǒng)服務(wù)價值的變化對于土地利用方式合理規(guī)劃與區(qū)域環(huán)境保護(hù)有著重要的意義。

眾多學(xué)者對生態(tài)系統(tǒng)服務(wù)功能已經(jīng)開展了大量的研究，價值評估方法日益完善。1997年，Costanza等[8]構(gòu)建了全球的生態(tài)系統(tǒng)服務(wù)價值評估法； 隨后謝高地等[9]據(jù)此構(gòu)建了中國大陸的生態(tài)系統(tǒng)服務(wù)價值當(dāng)量法。然而這些價值評估方法是針對全球和國家等大尺度的，在應(yīng)用于區(qū)域生態(tài)系統(tǒng)服務(wù)價值評估時，由于空間差異較大，評估結(jié)果的準(zhǔn)確性有待提高。此外，國內(nèi)大多數(shù)研究，在借鑒謝高地的價值當(dāng)量表時，對于不同時序往往采用相同的價值當(dāng)量[10-11]。這種靜態(tài)評價方法未能考慮由于時間變化導(dǎo)致服務(wù)價值發(fā)生的變化，因此有必要構(gòu)建基于時空尺度的生態(tài)系統(tǒng)服務(wù)價值動態(tài)評估模型，更客觀地闡明不同區(qū)域、不同時序下的生態(tài)系統(tǒng)服務(wù)價值。

本文綜合考慮生態(tài)系統(tǒng)服務(wù)價值的空間異質(zhì)性和時間演化性，選取黃土高原生態(tài)系統(tǒng)服務(wù)價值動態(tài)評估模型[12]，開展黃土高原2000年和2010年生態(tài)系統(tǒng)服務(wù)價值評估與時空特征分析。結(jié)合黃土高原2000—2010年間的土地利用變化，通過生態(tài)系統(tǒng)服務(wù)價值變率指數(shù)分析土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響。最后選取環(huán)境經(jīng)濟(jì)協(xié)調(diào)度指數(shù)(environmental economic coordination degree，CDEE)分析黃土高原不同區(qū)域的環(huán)境變化與社會經(jīng)濟(jì)發(fā)展之間的關(guān)系，評估黃土高原環(huán)境與經(jīng)濟(jì)協(xié)調(diào)發(fā)展?fàn)顩r，以期為黃土高原生態(tài)環(huán)境治理和保護(hù)成效評估提供科學(xué)依據(jù)。

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

1.1 研究區(qū)概況

黃土高原位于N33°41′～41°16′，E100°52′～114°33′之間，跨青海省、甘肅省、寧夏回族自治區(qū)、陜西省、內(nèi)蒙古自治區(qū)、山西省及河南省7個省區(qū)，總面積約62.4萬km2，四周高山環(huán)繞，南至秦嶺，北達(dá)陰山，西以日月山、賀蘭山為界，東到太行山。該區(qū)位于我國第二級地形階梯上，由東南向西北地勢逐漸升高，以六盤山和呂梁山為界將黃土高原劃分為東、中、西3個部分。東部海拔為500～1 000 m，中部海拔為1 000～2 000 m，西部海拔升至2 000～3 000 m。黃土高原地貌以黃土塬、梁、峁為主，地貌破碎且溝壑眾多[13]。全區(qū)年均氣溫在3.6～14.3 ℃之間，由東南向西北跨暖溫帶半濕潤區(qū)、暖溫帶半干旱區(qū)、中溫帶半干旱區(qū)[14]； 水資源匱乏且區(qū)域分布不均，年降水為150～750 mm，降水季節(jié)性差異大，7—8月多發(fā)暴雨，冬季枯水期長。土地覆蓋類型以草地和農(nóng)田為主，濕地水體面積小(圖1)。黃土高原生態(tài)環(huán)境問題嚴(yán)重，主要表現(xiàn)為自然災(zāi)害頻繁、水資源短缺、土地荒漠化及土壤侵蝕等[15]。

圖1 黃土高原生態(tài)系統(tǒng)空間分布Fig.1 Distribution of the ecosystems in Loess Plateau

1.2 數(shù)據(jù)源及其預(yù)處理

本文采用的數(shù)據(jù)包括土地利用、土地覆蓋專題數(shù)據(jù)和社會經(jīng)濟(jì)數(shù)據(jù)2大類。

土地利用、土地覆蓋專題數(shù)據(jù)是中國科學(xué)院遙感與數(shù)字地球研究所等單位利用遙感監(jiān)測方法分別構(gòu)建的1∶10萬比例尺長時間序列中國土地利用變化數(shù)據(jù)庫和1∶25萬比例尺長時間序列中國土地覆蓋數(shù)據(jù)[16-17]。在此基礎(chǔ)上，利用地理信息系統(tǒng)技術(shù)構(gòu)建了黃土高原地區(qū)2000—2010年間土地利用變化一級分類的100 m×100 m柵格數(shù)據(jù)、2000年和2010年土地覆蓋一級分類的1 km×1 km柵格數(shù)據(jù)以及2010年黃土高原生態(tài)系統(tǒng)空間分布，如圖1所示。土地覆蓋數(shù)據(jù)用于計算生態(tài)系統(tǒng)服務(wù)價值，土地利用變化數(shù)據(jù)用于分析黃土高原土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響。

社會經(jīng)濟(jì)數(shù)據(jù)包括國內(nèi)生產(chǎn)總值(gross domestic product，GDP)、三大糧食產(chǎn)量、種植面積比例及單位面積凈利潤等。GDP來自研究區(qū)內(nèi)各省區(qū)統(tǒng)計年鑒，用于環(huán)境經(jīng)濟(jì)協(xié)調(diào)度分析； 后3項來自全國農(nóng)產(chǎn)品成本收益資料匯編，用于生態(tài)系統(tǒng)服務(wù)價值計算。

2 研究方法

2.1 生態(tài)系統(tǒng)服務(wù)價值變率

在社會經(jīng)濟(jì)發(fā)展過程中，人們頻繁改變土地利用方式，使得區(qū)域生態(tài)系統(tǒng)內(nèi)部結(jié)構(gòu)和外部形態(tài)發(fā)生變化，進(jìn)而影響生態(tài)系統(tǒng)服務(wù)價值。通過構(gòu)建生態(tài)系統(tǒng)服務(wù)價值變率(ecosystem service value change rate，CR)[18-19]，研究土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響，計算公式為

(1)

式中：ESVm和ESVn分別為研究區(qū)某一類土地利用調(diào)整前后的生態(tài)系統(tǒng)服務(wù)價值；Am和An分別為研究區(qū)某一類土地利用調(diào)整前后的面積，m=1,2,…,6，n=1,2,…,6。

CR能夠表示某一土地利用類型的面積變化對其生態(tài)系統(tǒng)服務(wù)價值的重要性。CR>1，表示該類土地利用方式的變化對生態(tài)系統(tǒng)服務(wù)價值的影響明顯，值越大，影響越大；CR<1，表示該類土地利用方式的變化對生態(tài)系統(tǒng)服務(wù)價值的影響較弱[20-21]。

2.2 環(huán)境經(jīng)濟(jì)協(xié)調(diào)分析

環(huán)境與經(jīng)濟(jì)的協(xié)調(diào)發(fā)展是指環(huán)境子系統(tǒng)與經(jīng)濟(jì)子系統(tǒng)各自組成要素之間相互配合，形成互利共生的關(guān)系，二者協(xié)調(diào)發(fā)展，從而促進(jìn)區(qū)域的整體利益不斷得到提升[22]。協(xié)調(diào)度是衡量系統(tǒng)或要素之間發(fā)展?fàn)顩r的量化指標(biāo)。根據(jù)生態(tài)環(huán)境與經(jīng)濟(jì)協(xié)調(diào)發(fā)展的關(guān)系，構(gòu)建CDEE[23-25]。CDEE指的是在研究階段，單位面積生態(tài)系統(tǒng)服務(wù)價值變化率ESVpr與單位面積GDP變化率GDPpr的比值，即

(2)

(3)

(4)

式中：ESVpm和ESVpn分別為研究階段始、末年份的單位面積生態(tài)系統(tǒng)服務(wù)價值；GDPpm和GDPpn分別為研究階段始、末年份的單位面積GDP。CDEE值越大，表明研究區(qū)生態(tài)系統(tǒng)服務(wù)功能變化與經(jīng)濟(jì)發(fā)展協(xié)調(diào)度越高，具體可分為以下4種情況： ①CDEE≥1，表示該區(qū)域生態(tài)經(jīng)濟(jì)發(fā)展水平為高度協(xié)調(diào)，即生態(tài)系統(tǒng)服務(wù)價值增速不低于經(jīng)濟(jì)發(fā)展增速； ②0

黃土高原整體的CDEE計算方法為： 首先，以各省區(qū)占黃土高原面積比例為權(quán)重，計算黃土高原單位面積GDP變化率； 然后，結(jié)合生態(tài)系統(tǒng)服務(wù)價值變化率，計算黃土高原整體的CDEE。

2.3 生態(tài)系統(tǒng)服務(wù)價值估算技術(shù)路線

本文技術(shù)路線為： 首先，采用黃土高原生態(tài)系統(tǒng)服務(wù)價值評估模型，計算2000年和2010年黃土高原生態(tài)系統(tǒng)服務(wù)價值，并分析其隨時間和空間的變化特征； 然后，分析2000—2010年間黃土高原土地利用變化的時空特征，在此基礎(chǔ)上評估土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響； 最后，采用CDEE分析2000—2010間年黃土高原生態(tài)系統(tǒng)服務(wù)價值變化與經(jīng)濟(jì)發(fā)展之間的協(xié)調(diào)關(guān)系。

3 結(jié)果與分析

3.1 黃土高原生態(tài)系統(tǒng)服務(wù)價值評估

2000年黃土高原的生態(tài)系統(tǒng)服務(wù)價值如表1所示。各生態(tài)系統(tǒng)服務(wù)價值總共為6 333.91億元。生態(tài)系統(tǒng)服務(wù)價值由大到小依次為草地、森林、農(nóng)田、濕地水體、荒漠和聚落。其中，草地生態(tài)系統(tǒng)和森林生態(tài)系統(tǒng)為黃土高原提供了約71%的生態(tài)系統(tǒng)服務(wù)價值，荒漠僅提供約0.1%的服務(wù)價值，聚落提供的生態(tài)系統(tǒng)服務(wù)價值為負(fù)值。從提高黃土高原的生態(tài)系統(tǒng)服務(wù)功能角度，保護(hù)和提升森林和草地的服務(wù)功能、荒漠治理、控制城鎮(zhèn)化過程中聚落的擴(kuò)張速度是有效的環(huán)境保護(hù)措施。按照服務(wù)功能的價值由大到小依次為保持土壤、維持生物多樣性、水文調(diào)節(jié)、氣候調(diào)節(jié)、廢物處理、食物生產(chǎn)、提供美學(xué)景觀、原材料生產(chǎn)和氣體調(diào)節(jié)。保持土壤服務(wù)價值居首位，其中森林生態(tài)系統(tǒng)和草地生態(tài)系統(tǒng)提供了約79%的保持土壤服務(wù)價值，因此，促進(jìn)荒漠向草地、林地的轉(zhuǎn)化，利用植被的固土能力減少荒漠的土壤流失是黃土高原水土流失治理的關(guān)鍵； 氣體調(diào)節(jié)服務(wù)價值相對最低，只占總服務(wù)價值的2%左右，主要是由于氣體調(diào)節(jié)需依賴于濕地，而黃土高原的濕地面積僅約52 km2，不足總面積的1%。

表1 2000年黃土高原生態(tài)系統(tǒng)服務(wù)價值Tab.1 Ecosystem service value of Loess Plateau in 2000 (億元)

從時間上來看，2000—2010年間，黃土高原生態(tài)系統(tǒng)服務(wù)價值增加了約3 056.56億元，增幅高達(dá)48%(表2)。其中，增量最大的是草地生態(tài)系統(tǒng)，增加了1 158.71億元，主要源于黃土高原草地面積最大且草地質(zhì)量有所提升，草地生態(tài)系統(tǒng)服務(wù)價值增加明顯； 其次是森林生態(tài)系統(tǒng)，服務(wù)價值的大量增加一方面得益于森林單位面積的服務(wù)價值高，另一方面是林地面積有較大的增加； 濕地水體位居第三，濕地水體的單位面積服務(wù)價值最大，加上研究階段濕地水體的面積有所增加，因此服務(wù)價值相應(yīng)增加； 農(nóng)田生態(tài)系統(tǒng)位居第四，農(nóng)田面積大、機(jī)械化耕種、農(nóng)藥化肥的使用是農(nóng)田生態(tài)系統(tǒng)服務(wù)價值增加的主要原因； 荒漠的服務(wù)價值略有增加，主要是社會發(fā)展系數(shù)增加使得服務(wù)價值計算的基數(shù)變大； 聚落的服務(wù)價值明顯下降，原因是聚落的服務(wù)價值為負(fù)值且黃土高原聚落的面積明顯增加。

表2 2000—2010年間黃土高原生態(tài)系統(tǒng)服務(wù)價值變化Tab.2 Ecosystem service value change in Loess Plateau from 2000 to 2010 (億元)

從空間上來看，黃土高原生態(tài)系統(tǒng)服務(wù)價值總體趨勢東南高、西北低(圖2)，原因是東南部的水熱條件更適宜植被生長。由于濕地水體的單位面積生態(tài)系統(tǒng)服務(wù)價值最大，生態(tài)系統(tǒng)服務(wù)價值最高的區(qū)域分布在河流水系地區(qū)，主要包括黃河河段、汾河、涇河及渭河； 其次是黃土高原東南部，因為對應(yīng)區(qū)域土地覆蓋類型基本為森林，森林的生態(tài)系統(tǒng)服務(wù)功能強(qiáng)，因此東南部的生態(tài)系統(tǒng)服務(wù)價值相應(yīng)較高； 生態(tài)系統(tǒng)服務(wù)價值較低的區(qū)域分布在黃土高原西部、西北部和黃土高原的甘肅地區(qū)，由于這些區(qū)域荒漠面積較大，荒漠的生態(tài)系統(tǒng)服務(wù)價值低； 生態(tài)系統(tǒng)服務(wù)價值最低的區(qū)域與聚落分布基本一致，聚落的生態(tài)系統(tǒng)服務(wù)價值為負(fù)值，故而對應(yīng)區(qū)域的生態(tài)系統(tǒng)服務(wù)價值最低。

(a) 2000年 (b) 2010年

圖2 2000年和2010年黃土高原生態(tài)系統(tǒng)服務(wù)價值

Fig.2EcosystemservicevalueofLoessPlateauin2000and2010

3.2 生態(tài)系統(tǒng)服務(wù)價值對土地利用變化的響應(yīng)分析

3.2.1 黃土高原區(qū)土地利用變化的時空特征

2000年，黃土高原土地覆蓋類型構(gòu)成中，草地比例最大，占41.9%； 農(nóng)田次之，占33.1%； 森林比例居中，占14.7%； 再次是荒漠，占6.4%； 聚落和濕地水體比例最小，分別只有2.4%和1.5%。可見，草地和農(nóng)田是黃土高原土地覆蓋類型的主體。表3為2000—2010年間黃土高原土地利用變化轉(zhuǎn)移矩陣。由表3可知，2000—2010年間，黃土高原土地利用變化面積約為15 376.51 km2，約占黃土高原總面積的2.5%。耕地的轉(zhuǎn)入和轉(zhuǎn)出面積分別為2 222.26 km2和7 358.54 km2，其中流失耕地主要轉(zhuǎn)出為草地、城鄉(xiāng)工礦/居民用地以及林地，新增耕地主要由草地轉(zhuǎn)入； 林地凈增加面積為2 708.95 km2，主要來源是耕地和草地，表明黃土高原退耕還林政策落實良好，部分草地逐漸向森林過渡； 草地面積凈減少1 360.92 km2，草地主要轉(zhuǎn)出為未利用土地、耕地和林地； 水域面積凈增加306.73 km2，得益于河道的治理疏通； 城鄉(xiāng)工礦/居民用地面積大幅增加，增加2 784.79 km2，主要是占用耕地和草地； 未利用土地面積凈增加696.73 km2，主要是沙地增加，表明黃土高原土地沙漠化在研究階段尚十分嚴(yán)重。

表3 2000—2010年間黃土高原土地利用變化轉(zhuǎn)移矩陣Tab.3 land use change transfer matrix in Loess Plateau from 2000 to 2010 (km2)

黃土高原土地利用變化的主要類型有8種： 耕地轉(zhuǎn)草地、耕地轉(zhuǎn)城鄉(xiāng)工礦/居民用地、耕地轉(zhuǎn)林地、草地轉(zhuǎn)未利用土地、草地轉(zhuǎn)林地、草地轉(zhuǎn)耕地、未利用土地轉(zhuǎn)草地以及草地轉(zhuǎn)城鄉(xiāng)工礦/居民用地，占總變化面積的80%。

3.2.2 黃土高原土地利用變化對生態(tài)系統(tǒng)服務(wù)價值變化的影響

將2000—2010年黃土高原生態(tài)系統(tǒng)服務(wù)價值變化圖(圖3(a))與2000—2010年黃土高原土地利用變化主要類型圖(圖3(b))疊加分析，發(fā)現(xiàn)耕地轉(zhuǎn)城鄉(xiāng)工礦/居民用地、草地轉(zhuǎn)未利用土地、草地轉(zhuǎn)耕地、以及草地轉(zhuǎn)城鄉(xiāng)工礦/居民用地區(qū)域，使對應(yīng)的生態(tài)系統(tǒng)服務(wù)價值呈現(xiàn)不同程度的下降，其中草地轉(zhuǎn)城鄉(xiāng)工礦/居民用地生態(tài)系統(tǒng)服務(wù)價值下降程度最大； 耕地轉(zhuǎn)草地、耕地轉(zhuǎn)林地、草地轉(zhuǎn)林地以及未利用土地轉(zhuǎn)草地區(qū)域的生態(tài)系統(tǒng)服務(wù)價值表現(xiàn)為不同程度的增加，其中耕地轉(zhuǎn)林地生態(tài)系統(tǒng)服務(wù)價值上升最顯著。

(a) 生態(tài)系統(tǒng)服務(wù)價值變化 (b) 土地利用變化主要類型

圖3 2000—2010年黃土高原生態(tài)系統(tǒng)服務(wù)價值變化及土地利用變化主要類型

Fig.3EcosystemservicevaluechangeandmajortypesoflandusechangeinLoessPlateaufrom2000to2010

2000—2010年間，黃土高原6類土地利用類型的CR均大于1(圖4)，表明黃土高原各土地利用類型變化均會對生態(tài)系統(tǒng)服務(wù)價值產(chǎn)生重要影響。

圖4 黃土高原CR空間分布Fig.4 Spatial distribution of CR in Loess Plateau

由圖4，CR最高的是草地，達(dá)13.45，草地面積變化率為3.6%，其生態(tài)系統(tǒng)服務(wù)價值變化率為48.2%。由于黃土高原獨特的地理條件，草地作為其最主要的土地覆蓋類型，在維持黃土高原生態(tài)系統(tǒng)方面有著十分重要的作用。林地變化導(dǎo)致其CR位居第二，林地生態(tài)系統(tǒng)的服務(wù)價值高，在水土保持、風(fēng)沙防治、維持生物多樣性等多個方面發(fā)揮著重要作用。城鄉(xiāng)工礦/居民用地和未利用土地對應(yīng)的CR并列第三，這2類生態(tài)系統(tǒng)的服務(wù)功能差，因此在發(fā)展社會經(jīng)濟(jì)的同時應(yīng)注意生態(tài)環(huán)境的保護(hù)，同時加大土地荒漠化治理力度。耕地單位面積的生態(tài)系統(tǒng)服務(wù)價值較低，其CR位居第五，黃土高原耕地面積大，減少的耕地主要轉(zhuǎn)為單位面積服務(wù)價值較高的林地和草地。水域的CR最低，但仍對黃土高原生態(tài)系統(tǒng)服務(wù)價值有重要影響。

3.3 環(huán)境經(jīng)濟(jì)協(xié)調(diào)分析

經(jīng)計算，2000—2010年間，黃土高原的CDEE為0.115 3，環(huán)境經(jīng)濟(jì)協(xié)調(diào)水平為低度協(xié)調(diào)，表明黃土高原的生態(tài)系統(tǒng)服務(wù)價值以及GDP均表現(xiàn)為增加，但是生態(tài)系統(tǒng)服務(wù)價值的增速低于GDP的增速。黃土高原各省區(qū)的環(huán)境狀況以及經(jīng)濟(jì)發(fā)展水平存在差異，各省區(qū)的環(huán)境經(jīng)濟(jì)協(xié)調(diào)度略有區(qū)別(表4)。

表4 2000—2010年間黃土高原CDEETab.4 Environmental economic coordination index in Loess Plateau from 2000 to 2010

根據(jù)表4，CDEE最高的是寧夏回族自治區(qū)，一方面是由于黃河由南至北橫貫寧夏，合理灌溉使得農(nóng)田生態(tài)系統(tǒng)服務(wù)價值增加，加上草地質(zhì)量的改善，使其生態(tài)系統(tǒng)服務(wù)價值增速較大； 另一方面是由于該地區(qū)經(jīng)濟(jì)發(fā)展相對于其余6省區(qū)較緩，大幅增加的生態(tài)系統(tǒng)服務(wù)價值使得其CDEE相對較高。內(nèi)蒙古自治區(qū)的CDEE僅次于寧夏，內(nèi)蒙古自治區(qū)草地質(zhì)量改善，加上荒漠治理效益逐漸顯露，生態(tài)系統(tǒng)服務(wù)價值增加明顯； 另一方面該地區(qū)社會經(jīng)濟(jì)快速發(fā)展，GDP變化率高于生態(tài)系統(tǒng)服務(wù)價值變化率。

4 結(jié)論與討論

通過構(gòu)建時空修正的生態(tài)系統(tǒng)服務(wù)價值評估模型，對黃土高原生態(tài)系統(tǒng)服務(wù)價值進(jìn)行評估，進(jìn)而分析土地利用變化對生態(tài)系統(tǒng)服務(wù)價值影響以及環(huán)境與經(jīng)濟(jì)的協(xié)調(diào)度，研究得到以下結(jié)論：

1)2000年和2010年黃土高原生態(tài)系統(tǒng)服務(wù)價值分別為6 333.91億元和9 390.46億元，增幅達(dá)48%。草地生態(tài)系統(tǒng)服務(wù)價值增量位居榜首，增加了158.71億元，主要源于黃土高原草地面積最大且草地質(zhì)量有所提升，草地生態(tài)系統(tǒng)服務(wù)價值增加明顯。

2)2000—2010年間，黃土高原6類土地利用類型的CR均大于1，表明黃土高原各土地利用類型變化均會對生態(tài)系統(tǒng)服務(wù)價值產(chǎn)生重要影響。CR最高的是草地，為13.45。結(jié)合黃土高原2000—2010年間土地利用變化主要類型分析可知，導(dǎo)致草地面積減少的主要原因是草地轉(zhuǎn)為林地、未利用土地以及耕地。草地轉(zhuǎn)林地，區(qū)域的生態(tài)系統(tǒng)服務(wù)功能增強(qiáng)，有利于黃土高原生態(tài)環(huán)境建設(shè)； 部分草地退化轉(zhuǎn)為沙地，表明黃土高原土地荒漠化治理任務(wù)尚十分艱巨； 在草地轉(zhuǎn)為耕地時，應(yīng)注意維持區(qū)域生態(tài)系統(tǒng)的穩(wěn)定。

3)黃土高原在2000—2010年間的CDEE為0.115 3，處于低度協(xié)調(diào)水平，表明黃土高原生態(tài)系統(tǒng)服務(wù)價值的增速低于GDP的增速，但經(jīng)濟(jì)發(fā)展沒有使得生態(tài)系統(tǒng)服務(wù)價值降低。CDEE最高的是寧夏回族自治區(qū)，一方面是因為農(nóng)田生態(tài)效益增加明顯，加上草地質(zhì)量的改善，使其生態(tài)系統(tǒng)服務(wù)價值增速較大； 另一方面是因為該地區(qū)經(jīng)濟(jì)發(fā)展相對于其余6省較緩，大幅增加的生態(tài)系統(tǒng)服務(wù)價值使得其CDEE相對較高。

本文在研究黃土高原土地利用變化對生態(tài)系統(tǒng)服務(wù)價值影響時，并未進(jìn)行定量分析，在后續(xù)工作中，希望能構(gòu)建生態(tài)系統(tǒng)服務(wù)價值對土地利用變化的響應(yīng)模型，開展土地利用變化與生態(tài)系統(tǒng)服務(wù)價值變化之間的定量分析； 另一方面，隨著國家一系列重大工程的實施，黃土高原生態(tài)環(huán)境會受到重要影響，后續(xù)工作可以嘗試重點分析重大生態(tài)工程的環(huán)境效益。