黃土高原生態系統服務價值動態評估與分析

張 瑜, 趙曉麗, 左麗君, 張增祥, 汪 瀟

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

黃土高原生態環境脆弱，區域自然、地理條件復雜，加之自然災害頻發及不合理的土地利用方式，使得生態環境不斷惡化，尤其是水資源短缺、水土流失及土地荒漠化等生態環境問題逐步加劇[1]。為保護黃土高原生態環境，實施了退耕還林還草[2]、三北防護林[3]等生態工程。已有研究表明，在一系列生態保護工程實施下，黃土高原的生態環境發生了變化[4]。生態環境的質量決定著生態系統服務價值[5]，因此，有必要深入研究黃土高原生態系統服務價值的時空演變過程，以此揭示該區域生態環境的變化過程。

隨著生態系統服務價值研究的不斷深入，價值評估方法日益完善。1997年價值當量法[6]首次被提出并被用于計算全球生態系統服務價值。隨后根據中國生態系統和社會經濟發展狀況，提出基于單位面積價值當量因子的生態系統服務價值評估體系[7]。但是該研究是針對全國尺度的，直接將其應用于計算區域生態系統服務價值時，評估結果的準確性有待提高。故學者們相繼提出基于專家知識當量修正[8]，生態區位修正[9]，社會生態經濟協調系數修正[10]及生態服務價值對土地利用變化的響應[11]等一系列研究。這些方法的改進使得生態系統服務價值計算更為準確，但只是單一方面的修正沒有綜合考慮同一生態系統位于不同地理位置的服務價值差異(即空間異質性)和社會經濟發展對生態系統服務價值的影響。

本文基于中國陸地生態系統單位面積生態服務價值當量表，綜合考慮生態系統服務價值的空間異質性和時間演化性，融合遙感與社會經濟多源數據，基于植被覆蓋度和NPP[12]開展空間異質修正[13]，基于社會發展系數和資源稀缺系數[14]開展時間修正，建立黃土高原生態系統服務價值動態評估模型；在此基礎上，開展黃土高原地區2000年、2005年、2010年間生態系統服務價值時空過程評估與動態特征分析，以期為黃土高原生態環境治理和保護成效評估提供科學依據。

1　研究區概況與數據來源

1.1　黃土高原概況

黃土高原地跨北緯33°41°—41°16′，東經100°52′—114°33′，東西長1 000余km，南北寬約750 km，處于我國第二級地形階梯上。黃土高原跨青海省、甘肅省、寧夏回族自治區、陜西省、內蒙古自治區、山西省及河南省7個省，總面積約62萬km2。黃土高原地勢西北高，東南低，黃土覆蓋厚度50～200 m[15]，為世界黃土面積覆蓋最大的高原。該區域植被稀疏，氣候較干旱，但夏季降水集中且雨量大，流水沖蝕作用強，因此在流水侵蝕作用下地表支離破碎，形成溝壑交錯其間，造成了嚴重的水土流失[16]。黃土高原水土流失面積達到42萬km2，年均輸入黃河泥沙16億t，是我國乃至世界上水土流失最嚴重、生態環境最脆弱的地區，對該區的治理乃是水土保持生態建設中最具挑戰性的問題[17]。

1.2　數據來源

所用數據包括遙感數據、專題數據和社會經濟數據三大類。其中，遙感數據有NPP和NDVI數據。NPP和NDVI數據分別來自美國NASA衛星中心發布的MOD17A3的凈初級生產力(NPP)和MOD13A3的歸一化植被指數(NDVI)，均為1 km的柵格數據。

土地覆蓋專題數據基于中國科學院遙感與數字地球研究所等單位利用遙感監測方法構建的中國1∶25萬比例尺長時間序列中國土地覆蓋數據庫，該數據庫是基于30 m分辨率的TM影像，采用人機交互目視解譯完成的全國多期土地覆蓋矢量數據，其土地覆蓋遙感監測分類系統包括6個一級分類，25個二級分類。一級分類為：森林、草地、農田、聚落、濕地水體、荒漠[18]。在此基礎上，利用GIS技術構建了2000年、2005年、2010年的黃土高原地區土地覆蓋一級分類的1 km柵格數據。

社會經濟數據包括恩格爾系數、GDP、城鎮化率、人口密度、“三廢”治理費用、三大糧食產量、種植面積比例及單位面積凈利潤。前五項來自中國統計年鑒和研究省份統計年鑒，后三項來自全國農產品成本收益資料匯編。

采用的投影為雙標準緯線等面積割圓錐投影，用全國統一的中央經線和雙標準緯線，中央經線為東經105°，雙標準緯線為北緯25°和北緯47°，所采用的橢球體是KRASOVSKY橢球體。

2　研究方法及改進

本文基于中國陸地生態系統單位面積生態服務價值當量表，考慮到該表是針對全國平均水平，故應用于建立黃土高原生態系統服務價值動態評估模型時，需進行以下4個方面的修正。首先，為了體現同類生態系統在不同區域服務價值的差異，對價值當量表進行空間異質修正，實現從全國到區域再到單位面積的尺度轉換。其次，不同的社會經濟發展水平下，人類對經濟發展與生態環境保護關系的認識深淺不一，導致對生態服務價值的評估有所變化，故需對價值當量表進行社會發展系數修正。再者，由于生態資源供給關系變化時，價值也會隨之改變，故也需對價值當量進行資源稀缺度修正。以上通過對價值當量表的修正實現對生態系統服務功能量的確定，同時也需考慮生態系統服務功能單位面積價值，本文計算不同時期黃土高原總服務價值時均采用2005年不變價格作為依據，從而減弱市場經濟的影響。經過以上步驟，建立黃土高原生態系統服務價值動態評估模型。

圖1研究方法改進技術路線

2.1　空間異質修正

同一類型生態系統在不同地區的內部結構和外部形態差異明顯，導致其生態系統服務價值也有所不同。研究發現植被覆蓋度的不同會導致生態系統服務功能的不同[19]；生物量在其形成和累積的過程中，對生態系統的服務功能會產生重要影響[6]，該要素在黃土高原尤為重要。黃土高原東西部水熱條件不同，導致不同位置的同一類生態系統植被狀況差異大，因此，本文選取植被凈初級生產力(NPP)和植被覆蓋度[9]來進行空間異質修正。

生態系統服務價值基準價格是全國平均水平狀況下的農田每年自然糧食產量的經濟價值，這里黃土高原選取三大糧食主產物水稻、小麥和玉米(由統計年鑒數據可知，三大糧食產物種植面積約占農作物種植總面積70%)，因此首先需要用研究區的糧食產量對價值當量進行粗修正，以實現價值當量表從全國到黃土高原的尺度轉換，公式如下：

(1)

式中：Q1為研究區的糧食修正系數；r和R，w和W以及c和C分別為研究區和全國的水稻、小麥、玉米每公頃產量。

一定區域的生態系統服務價值總量可表述為：

(2)

式中：i=1,2，…，6，表示生態系統的類型；ESVi表示第i類生態系統服務價值。

(3)

式中：j=1,2,…，9，表示第i類生態系統的第j種生態服務功能；t=1,2，…，n,表示一定區域內第i類生態系統在空間上分布的像元數；Qjt為研究區任一像元NPP和植被覆蓋度(植被覆蓋度由NDVI計算所得)的修正系數，NPPt和ft為研究區t像元的凈初級生產力和植被覆蓋度；NPPmean和fmean為研究區第i類生態系統植被凈初級生產力的均值和植被覆蓋度的均值。

通過對黃土高原進行NPP和植被覆蓋度修正，實現價值當量表從黃土高原到像元尺度的轉換。因此，研究區的空間異質修正系數Q可表述為：

Q=Q1×Qjt

(4)

2.2　社會發展系數修正

人類對生態系統服務功能的認識是隨著時間不斷變化的，不同的社會經濟發展狀況下，人類對經濟發展與生態環境保護的認識深淺不一。黃土高原跨7個省域，不同省域社會經濟發展差異明顯，通過查閱大量相關文獻，這里選取黃土高原7省的恩格爾系數[14]、GDP[20]和城鎮化率[21]3個指標進行社會發展系數修正。恩格爾系數是食品支出總額占個人消費支出總額的比重。GDP是核算國民經濟的核心指標，也是衡量一個國家或地區總體經濟狀況的重要指標。城鎮化是人口從農村向城市遷移聚集的過程，同時也帶來地域景觀、產業結構及生活方式的變化。采用經濟學上的皮爾生長曲線模型來表示社會發展階段系數。

(5)

式中：l為與現實支付意愿有關的社會發展階段系數；L為極富階段的支付意愿，取值為1；e為自然對數；En為恩格爾系數。

由于各地區的福利待遇、價格體系等差異較大，只用恩格爾系數來反應社會發展階段系數可能導致誤差，因此，引入GDP和城鎮化率來進一步修正社會發展階段系數，如下：

(6)

式中：D為研究區的社會發展系數，l，U和GDP為研究區各省的社會發展階段系數、城鎮化水平和GDP產值；lmean，Umean和GDPmean為全國平均水平的社會發展階段系數、城鎮化水平和GDP均值。

2.3　資源稀缺系數修正

資源稀缺度是反映某區域在經濟、社會和環境的可持續發展過程中生態資源需求量與供給量之間的關系[22]。生態系統資源稀缺度變化時，依據經濟學上的供求關系，其價值也會發生變化。黃土高原人口密度差異大，生態資源需求量與人口總數[12]呈線性相關性，因此選用人口密度來進行資源稀缺系數修正，如下：

(7)

式中：S是資源稀缺系數；p和P是研究區人口密度和全國平均人口密度。

2.4　黃土高原生態系統服務功能單位面積價值的確定

1個標準單位生態系統生態服務價值當量因子是指1 hm2全國平均產量的農田每年自然糧食產量的經濟價值。以三大糧食主產物稻谷、小麥和玉米單位面積農田生態系統糧食生產的凈利潤按種植面積比例加權分配，作為1個標準生態服務價值當量因子的價值量。2000年、2005年及2010年糧食的凈利潤受市場經濟的影響很大，所以如果用當年凈利潤計算生態系統服務價值的話，無法區分生態系統服務價值的變化是糧食凈利潤的變化還是生態系統服務功能的變化導致的。因此這里采用2005年糧食的凈利潤作為不變價格來計算黃土高原生態服務價值當量因子的價值量[14]。公式如下：

(8)

式中：P是黃土高原1個標準單位當量因子的價值量；m=1,2,3分別指稻谷、小麥及玉米；N是1 kg糧食主產物的凈利潤；Y是單位面積糧食產量；am是黃土高原稻谷、小麥及玉米的種植面積；A是黃土高原三大糧食主產物種植總面積。

經過計算，黃土高原在2005年1個標準單位當量因子的價值量是1 843.1元/hm2，并根據公式10計算得到的黃土高原2000年、2005年及2010年當量的修正系數，結合中國生態系統單位面積生態服務價值當量表，可得到森林、草地、農田、濕地水體以及荒漠的生態系統服務功能單位面積價值。

由于中國陸地生態系統服務價值當量表中沒有考慮聚落的服務價值，因此單獨考慮聚落的生態系統服務功能單位面積價值。其中，聚落在氣體調節價值、廢物處理價值和水文調節價值采用防治成本法計算，通過各省的廢氣、廢水、廢物治理投入金額，以研究區與各省聚落面積比例作為權重，通過面積加權計算黃土高原聚落的生態系統服務價值；聚落提供美學景觀價值，同樣的按照面積加權方式取各省旅游收入的1/3[23]。

2.5　黃土高原生態系統服務價值評估模型建立

在進行了空間異質修正、社會發展系數修正、資源稀缺系數修正及確定生態系統服務功能單位面積價值后，得到研究區生態系統服務功能總價值為：

(9)

Rjt=Q×D×S

(10)

式中：i,j,t如上所述；ESVij為第i類生態系統的第j種生態服務功能的單位面積價值；Rjt表示每個像元的修正系數；Ajt表示各個像元的面積，如果為等面積投影，則為給定的常數，ESVs為生態系統服務總價值。

在ArcGIS平臺支持下，對空間異質修正結果按省邊界進行分割，與計算所得的各省社會發展系數、資源稀缺系數進行疊加，得到黃土高原價值當量修正結果。然后按6類土地覆蓋類型進行分割，與中國陸地生態系統服務價值當量表及生態系統服務功能單位面積價值進行疊加，計算可得黃土高原2000年、2005年及2010年的生態系統服務價值。

3　結果與分析

3.1　生態系統服務價值當量修正結果

3.1.1空間異質修正結果黃土高原3 a空間異質系數大于1的區域主要分布在黃土高原東部、南部和西南部(圖2)，表明這些區域的植被覆蓋度和凈初級生產力高于黃土高原平均水平。空間異質系數低于0.5的區域集中在黃土高原西北、西部，主要是這些區域的自然條件比黃土高原平均水平差；結合土地覆蓋類型圖發現，這些區域與荒漠分布大致吻合，荒漠的水熱條件相對差，導致植被覆蓋度和凈初級生產力低于黃土高原平均水平。

比較黃土高原不同時期的空間異質修正系數圖，可以看到修正系數為1～2的區域逐漸由東南向西北擴張，修正系數為0～0.5的區域表現為先增加后減小。受地理環境因素的影響，不同地區同種土地覆蓋類型的植被覆蓋度和凈初級生產力差異明顯，進而影響生態系統服務價值，因此對黃土高原進行空間異質修正是很有必要的。

3.1.2社會發展系數修正結果除內蒙古自治區2005年和2010年的社會發展系數大于1外，黃土高原其余6省在2000年、2005年及2010年社會發展系數均低于1(圖3A)，說明黃土高原社會經濟發展整體水平低于全國平均水平。這主要是因為黃土高原地理環境和氣候條件較差，經濟發展速度低于全國平均水平，而內蒙古自治區能源礦產豐富、畜牧業發達，所以社會發展系數高。從社會發展系數變化角度分析，2000—2010年黃土高原各省的社會發展系數不同程度的增加。

3.1.3資源稀缺系數修正結果黃土高原7省在2000—2010年各省資源稀缺系數變化不大，主要是各省人口總數變化不明顯(圖3B)。比較黃土高原7省的資源稀缺系數，差異明顯。資源稀缺系數最大的是河南省，其次是山西、陜西，并且這3個省的資源稀缺系數大于1，即相比于全國平均水平，這3個省由于人口總數大導致資源更緊缺；寧夏回族自治區和甘肅省資源稀缺系數介于0.8～1；內蒙古自治區資源稀缺系數約為0.6，低于全國平均水平，是由于內蒙古自治區自然資源儲量大且人口密度小；青海省資源稀缺系數最低，主要是由于人口密度小。

圖2　黃土高原空間異質修正系數

圖3　黃土高原社會發展系數修正(A)及黃土高原資源稀缺系數修正(B)

3.2　黃土高原生態系統服務價值時空特征分析

采用黃土高原生態系統服務價值評估模型，計算黃土高原2000年、2005年及2010年的生態系統服務價值(圖4)。

圖4　黃土高原生態系統服務價值

3.2.1黃土高原生態系統服務價值時間變化分析黃土高原在2000年、2005年及2010年的生態系統服務總價值分別為6 333.91億元、7 679.95億元及9 390.46億元，整個時段增幅達48%。2010年黃土高原9項生態系統服務功能價值由大到小依次為：保持土壤>維持生物多樣性>水文調節>氣候調節>廢物處理 >食物生產>提供美學景觀>原材料生產>氣體調節(表1)。黃土高原生態系統服務功能中保持土壤價值居首位，約占黃土高原生態系統服務功能總價值的1/4；其中，森林和草地提供的保持土壤價值約占80%，這說明黃土高原治理水土流失問題的關鍵在于保持和提高森林和草地的質量。維持生物多樣性價值居第二位，約占總服務價值的1/5。黃土高原生態系統服務價值中，氣體調節價值較低，只占總服務價值2%左右，主要是因為黃土高原濕地和森林少，而氣體調節主要依賴于濕地和森林生態系統[5]。

2000—2005年，黃土高原森林生態系統和草地生態系統服務價值增量約占服務價值總增量的70%，表明黃土高原生態環境的改善主要得益于森林和草地面積的增加。其次，對黃河河道清淤、上游生態建設減少水土流失等措施收效顯著，故濕地水體生態服務價值增加，且增量約占總增量的18%。農田生態系統服務價值也有較大增加，約占服務價值總增量的11%，主要源于農藥化肥的使用、農田水利灌溉設施的建設及機械化耕作等方面。國家對荒漠的治理使得荒漠生態環境逐漸改善，其服務價值增加近2.8億元。

表1　2010年黃土高原各土地覆蓋類型生態系統服務價值　億元

2005—2010年，黃土高原生態服務價值持續增加，其中約40%的服務價值增量由草地貢獻，退耕還草生態工程效益明顯，加上合理放牧使得草地生態恢復。森林生態服務價值增速較前一時段有所降低，下降8%左右，原因是前期造林多、收效快，后一階段造林速度減緩。濕地水體服務價值增加約318億元。農田服務價值增量約占價值總增量的13%，增速較2000—2005年提高2%，主要是農民耕作水平的進一步提高和機械化的推廣。聚落的生態系統服務價值減少速度變緩，原因是對于工業三廢的治理力度加強及居民生活垃圾集中處理，降低對環境的污染。

3.2.2黃土高原生態系統服務價值空間差異分析從空間上來看，生態系統服務價值高的區域分布在山西、河南及陜西中部和南部，生態系統服務價值增速最快的是內蒙古自治區。對比不同時期黃土高原生態系統服務價值圖發現，黃土高原東南部單位面積生態系統服務價值為200～600的區域面積在增大，并且單位面積價值也有所增加(圖4)。這是由于東南部黃土厚度高于西北部，并且東部海拔比西部低，更適宜植被生長；此外，由于三北防護林[9]的生態工程建設，東南部森林的面積增加，也使得該區域生態系統服務價值增加。黃土高原西北部的生態系統服務價值也表現為增加，結合土地覆蓋數據分析，主要是由于退耕還草的政策[3]落實良好，使得草地面積增加，生態環境質量改善，該區生態系統服務價值增加。2000—2005年，土地覆蓋類型變化較大，主要發生在內蒙古、寧夏北部及甘肅東部、陜西北部，主要方式是農田轉換為草地和林地。2005—2010年土地覆蓋類型變化相對減弱，比較明顯的是內蒙古自治區的黃河河段附近區域植被覆蓋度增加。

結合土地覆蓋變化分析各省生態系統服務價值的空間特點，發現內蒙古自治區生態系統服務價值增加顯著，這與該區域社會經濟發展快、草場面積增加等關系密切。陜西省中部和南部森林和草地分布廣且質量改善，因此生態服務價值高、增加明顯；陜西北部草地覆蓋度從2000—2010年呈增加趨勢，所以生態服務價值相應增加。山西省耕地以水田和水澆地為主，水田和水澆地相比于旱地植被長勢更好，且山西省森林面積大，故生態系統服務價值最高。寧夏回族自治區土地覆蓋類型以草地和農田為主，全省生態系統服務價值增加明顯，農田耕種水平提高、社會發展系數提高是其增加的主要原因。甘肅省主要土地覆蓋類型是草地和農田，但由于社會發展系數低、海拔高及水熱條件不適宜，導致生態系統服務功能較低且增加不明顯。青海省土地覆蓋類型以草地、農田和森林為主，該區自然條件較差，導致糧食產量低，并且資源稀缺度低，因此青海省生態系統服務價值低。河南省只有小部分區域屬于黃土高原，主要土地覆蓋類型為農田和森林，該區域水源充足，氣候適宜，植被生長狀況好，且資源稀缺度高，因此生態系統服務價值高。

4　結論與討論

(1) 建立了黃土高原生態系統服務價值動態評估模型。本文結合黃土高原地域特色，對生態系統服務價值當量進行修正，提高了原有生態系統服務價值評估模型的時空適用性。空間尺度上利用植被覆蓋度和NPP數據進行修正，能夠直觀的顯示不同區域的生態服務價值的差異；同時，考慮社會經濟發展和資源稀缺度，探索不同時期生態系統服務價值當量隨時間的變化，使得生態服務價值評估模型能夠反映社會發展的映射作用。

(2) 黃土高原生態系統服務價值模型計算2000年、2005年及2010年的生態系統服務價值分別為6 333.91億元、7 679.95億元及9 390.46億元，從2000—2010年，增幅達48%。2000—2005年，黃土高原森林生態系統和草地生態系統服務價值增量約占服務價值總增量的70%。2005—2010年，生態系統服務價值總增量約有40%由草地貢獻。在黃土高原各項生態系統服務功能中，保持土壤價值居首位，并且森林和草地提供了近80%的服務價值，因此黃土高原治理水土流失問題的關鍵在于保持和提高森林和草地的質量。

(3) 黃土高原各省的生態服務價值存在差異，社會經濟和自然地理環境共同影響著生態系統服務價值。生態系統服務價值高的區域分布在山西、河南及陜西中部和南部，生態系統服務價值增速最快的是內蒙古自治區。寧夏回族自治區生態服務價值增加明顯，青海、甘肅省服務價值低且無明顯增加。因此，為了提高黃土高原的生態系統服務價值，應該綜合各省的社會經濟和自然地理環境，因地制宜的改善黃土高原的生態系統狀況。

對生態系統服務價值評估方法進行修正時，考慮了空間異質性、社會發展狀況及資源稀缺度，在選取指標時具有一定的主觀性，缺乏對眾多因素的比較評價，在接下來的工作中，希望在這方面做出進一步的研究。另一方面為了分析黃土高原生態系統服務價值與社會經濟之間的相互影響機制，嘗試利用敏感性分析[21]、社會經濟協調度分析[9]等方法。

參考文獻：

[1]孟慶香.基于遙感、GIS和模型的黃土高原生態環境質量綜合評價[D].陜西楊凌:西北農林科技大學,2006.

[2]許炯心.黃土高原生態環境建設的若干問題與研究需求[J].水土保持研究,2000,7(2):10-13.

[3]Duan H, Yan C, Tsunekawa A, et al. Assessing vegetation dynamics in the Three-North Shelter Forest region of China using AVHRR NDVI data[J]. Environmental Earth Sciences, 2011,64(4):1011-1020.

[4]Liu J, Li S, Ouyang Z, et al. Ecological and socioeconomic effects of China′s policies for ecosystem services[J]. Proceedings of the National Academy of Sciences, 2008,105(28):9477-9482.

[5]Du X, Huang Z. Ecological and environmental effects of land use change in rapid urbanization:The case of hangzhou, China[J]. Ecological Indicators, 2017,81:243-251.

[6]Costanza R, d′Arge R, Groot R D, et al. The value of the world′s ecosystem services and natural capital[J]. Nature, 1998,25(1):3-15.

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

[8]謝高地,甄霖,魯春霞,等.一個基于專家知識的生態系統服務價值化方法[J].自然資源學報,2008，23(5):911-919.

[9]唐秀美,陳百明,路慶斌,等.生態系統服務價值的生態區位修正方法[J].生態學報,2010,30(13):3526-3535.

[10]宋佳楠,梅建屏,金曉斌,等.基于協調系數修正的區域生態系統服務價值測算研究[J].地理與地理信息科學,2010，26(1):86-89.

[11]曹銀貴,白中科,景明,等.土地利用變化對區域生態服務價值的影響研究[J].水土保持研究,2013,20(6):256-261.

[12]潘耀忠,史培軍,朱文泉,等.中國陸地生態系統生態資產遙感定量測量[J].地球科學,2004,34(4):375-384.

[13]張曉楠,宋宏利,李振杰.基于地統計學的區域生態服務價值空間分異規律研究[J].水土保持研究,2012,19(6):168-171.

[14]胡喜生,洪偉,吳承禎.土地生態系統服務價值動態估算模型的改進與應用:以福州市為例[J].資源科學,2013,35(1):30-41.

[15]高海東,李占斌,李鵬,等.基于土壤侵蝕控制度的黃土高原水土流失治理潛力研究[J].地理學報,2015,70(9):1503-1515.

[16]Fu B, Liu Y, Lü Y, et al. Assessing the soil erosion control service of ecosystems change in the Loess Plateau of China[J]. Ecological Complexity, 2011,8(4):284-293.

[17]姚文藝,李勉.黃土高原土壤侵蝕及綜合治理研究評述[J].中國水土保持,2005(4):15-17.

[18]張,増祥.中國土地覆蓋遙感監測[M].北京:星球地圖出版社,2010.

[19]畢曉麗,葛劍平.基于IGBP土地覆蓋類型的中國陸地生態系統服務功能價值評估[J].山地學報,2004,22(1):48-53.

[20]蘇飛,張平宇.基于生態系統服務價值變化的環境與經濟協調發展評價:以大慶市為例[J].地理科學進展,2009,28(3):471-477.

[21]姚成勝,劉耀彬.福建省生態系統服務價值變化對土地利用變化驅動因子的敏感性分析[J].農業系統科學與綜合研究,2010,26(1):80-85.

[22]徐中民,龍愛華.中國社會化水資源稀缺評價[J].地理學報,2004,59(6):982-988.

[23]武文婷.杭州市城市綠地生態服務功能價值評估研究[D].南京:南京林業大學,2011.