戴云哲, 李江風

(中國地質大學(武漢) 公共管理學院, 武漢 430074)

生態(tài)用地作為地球的“腎”和“肺”，具有多種生態(tài)功能[1]，同時為人類提供了供給服務、調節(jié)服務、支持服務和文化服務[2]。在《全國生態(tài)環(huán)境保護綱要》中，生態(tài)用地數(shù)量和質量的保護被重點強調。2010年頒布的《全國主體功能區(qū)規(guī)劃》單獨提出了生態(tài)功能區(qū)的概念并在地方主體功能區(qū)規(guī)劃中落實到地塊。隨著生產(chǎn)建設對土地需求的逐年增長，生態(tài)用地被擠占的速度也隨之加快，生態(tài)環(huán)境面臨挑戰(zhàn)[3]。將生態(tài)系統(tǒng)提供的服務量化為生態(tài)服務價值能夠清晰直觀地了解其時空變化的趨勢[4]，并進一步進行環(huán)境承載能力和生態(tài)風險[5]等評估。目前生態(tài)服務價值的研究在方法主要從生態(tài)系統(tǒng)產(chǎn)品價值[6-7]和土地利用覆被[8-9]兩個角度入手，研究主要集中于對區(qū)域尺度或特定地類生態(tài)服務價值總量的測算和變化特征的描述[10-12]。在地形對土地利用[13]的影響引起關注后，有學者對生態(tài)服務價值的垂直空間分布進行了探索，發(fā)現(xiàn)不同土地利用類型所提供的生態(tài)服務價值隨地形起伏呈現(xiàn)出明顯的垂直分布規(guī)律[14-15]，但土地利用結構和生態(tài)服務價值時間演化的地形梯度效應有待進一步研究。

洞庭湖區(qū)囊括了湖泊、平原、山地等多種地貌形態(tài)，生態(tài)用地數(shù)量隨地形變化顯著，根據(jù)洞庭湖區(qū)大型湖泊流域的特殊性，結合土地的空間形態(tài)[16]、生態(tài)功能[17]和主體功能[18]等不同視角將生態(tài)用地定義為：除生產(chǎn)生活用地外，以提供生態(tài)服務功能為主的用地類型，具體到洞庭湖區(qū)包括濕地、林地和草地。受長期高強度農(nóng)業(yè)生產(chǎn)與近年來快速城鎮(zhèn)化運動的影響，洞庭湖區(qū)存在湖泊萎縮、生態(tài)環(huán)境惡化等問題[19]。針對這些問題，退耕環(huán)湖、還林和濕地公園建設等生態(tài)修復與保護工程相繼實施，生態(tài)用地的變化也會更加活躍。通過不同時期生態(tài)用地及其生態(tài)服務價值垂直空間分布的變化，揭示洞庭湖區(qū)生態(tài)用地時空演化與生態(tài)服務價值的地形梯度效應，為合理配置生態(tài)用地的數(shù)量結構與空間布局和生態(tài)環(huán)境的保護與提升提供參考。

1　研究區(qū)概況

洞庭湖位于長江中游南岸，是我國第二大淡水湖，屬于過水性洪道型湖泊，綜合考慮了流域與行政區(qū)劃因素劃定的洞庭湖區(qū)涉及岳陽、益陽、常德市共計21個縣市區(qū)，以及湖北省的松滋、公安、石首等縣市，土地總面積為37 036.28 km2。洞庭湖區(qū)河網(wǎng)密布，生態(tài)要素高度聚集，對維護區(qū)域生態(tài)平衡有重要作用。20世紀大范圍的圍湖造田和造成了洞庭湖區(qū)濕地的劇烈萎縮，林地也因毀林開荒而明顯減少。1988年長江發(fā)生特大洪災后，為調節(jié)洞庭湖區(qū)水文調節(jié)和水源涵養(yǎng)功能，實施了退耕還湖、還林工程，洞庭湖區(qū)濕地和林地得到了有效保護和恢復。然而近年來隨著城鎮(zhèn)化進程的加速，建設用地擴張對耕地安全造成了威脅，在耕地紅線的壓力下，生態(tài)用地逐漸成為被占用的對象。

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

2.1　數(shù)據(jù)來源與處理

1990年、1995年、2000年、2005年、2010年、2015年6期土地利用覆被數(shù)據(jù)為中國科學院資源環(huán)境科學數(shù)據(jù)中心制作的中國遙感影像解譯數(shù)據(jù)，該解譯數(shù)據(jù)對全國土地利用覆被劃分了6個一級地類和25個二級地類，柵格分辨率為30 m×30 m[20]。從中提取生態(tài)用地，根據(jù)柵格數(shù)量及分布變化可計算出生態(tài)用地面積變化和轉移去向與來源。計算地形梯度所需的高程和坡度數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺(http:∥www.gscloud.cn)。

2.2　土地利用轉移速率和動態(tài)度

土地利用轉移速率是某土地利用類型在一定時期內的轉移量。計算公式如下：

(1)

式中：Mi為第i種土地利用類型的轉移速率；Aib，Aia分別為研究基期、末期第i種土地利用類型的總面積；t為研究時間段。

土地利用動態(tài)度代表研究區(qū)土地利用變化的綜合活躍程度[21]。其定義為：

(2)

式中：D土地利用綜合動態(tài)度；Ai研究基期第i類土地利用類型的面積；ΔAi-j為基期到末期第i類土地利用類型轉化為其他類型的數(shù)量；t為時間段；n為土地利用類型總數(shù)。

2.3　地形位指數(shù)和地形分布指數(shù)

高程和坡度對地理要素的分布都具有顯著影響，參照以往研究成果[13,22]，綜合考慮高程和坡度的空間差異，引入地形位指數(shù)，計算公式為：

(3)

土地利用類型分布指數(shù)無量綱的屬性規(guī)避了不同土地利用類型的轉移在不同地形梯度上分布頻率的可比性受到各梯度總面積差異帶來的量綱影響[23]，能夠表征不同地形梯度上各土地利用類型分布的頻率高低，定義如下：

(4)

式中：Lie為土地轉移類型分布指數(shù)；A為土地總面積；i為土地利用類型；e為地形梯度。Lie>1的區(qū)間為第i類土地利用類型分布的優(yōu)勢地形區(qū)，Lie值越大，表明此土地利用類型在該地形梯度出現(xiàn)的頻率越高。土地轉移分布指數(shù)參考了土地利用類型分布指數(shù)，對不同地形梯度生態(tài)用地的轉移進行比較，反映生態(tài)用地的變化垂直分布特征，計算公式為：

(5)

式中：TsIke為土地轉移類型分布指數(shù)；TsA為土地轉移面積；k為土地轉移類型；e為地形梯度。TsAke>1的梯度為某土地轉移類型分布的優(yōu)勢地形區(qū)，且出現(xiàn)的頻率隨TsIke值的增加而升高。

2.4　生態(tài)服務價值及其地形分布指數(shù)

在歸納歷史生態(tài)系統(tǒng)公益價值評價研究的過程中，Costanza等人提出了生態(tài)系統(tǒng)的服務功能的貨幣化計算方法[24]：

(6)

式中：ESV為生態(tài)服務總價值；Ai為第i種土地利用類型的總面積；VCi為單位面積第i種土地利用類型的生態(tài)服務價值。

謝高地等在該生態(tài)系統(tǒng)服務功能分類的基礎上，結合中國實際情況，定義了農(nóng)田生態(tài)系統(tǒng)糧食生產(chǎn)的服務價值為1，其他生態(tài)系統(tǒng)服務價值與其比值為該生態(tài)系統(tǒng)當量因子的衡量標準[25]，并于2015年進一步改良和細化[26]。武愛彬等于2017年在壩上高原生態(tài)用地生態(tài)服務價值演變的研究中結合中國科學院資源環(huán)境科學數(shù)據(jù)中心制作的土地覆被數(shù)據(jù)中地類劃分的差異進行了相應調整[27]。綜合以上研究方法和研究成果，歸納得出洞庭湖區(qū)生態(tài)用地生態(tài)服務價值當量因子表(表1)，并以1 hm2農(nóng)田生態(tài)系統(tǒng)糧食生產(chǎn)的凈利潤作為糧食生產(chǎn)的服務價值。為使不同時期的生態(tài)服務價值更具可比性，以《全國農(nóng)產(chǎn)品成本收益資料匯編2015》中2015年洞庭湖區(qū)糧食收購價格均價為基礎(湖北省128.77元/50 kg，湖南省131.54元/50 kg)[28]，依據(jù)1990—2015年湖北省、湖南省統(tǒng)計年鑒中的糧食產(chǎn)量數(shù)據(jù)計算得出洞庭湖區(qū)在1990，1995，2000，2005，2010和2015年的農(nóng)田生態(tài)系統(tǒng)糧食生產(chǎn)的服務價值分別為1 336.54，2 142.04，1 970.92，2 148.40，2 236.79，2 150.02元/hm2。

為了更準確地定位生態(tài)服務價值的優(yōu)勢分布區(qū)間，參考土地利用類型分布指數(shù)的計算方法，獲得生態(tài)服務價值分布指數(shù)：

(7)

式中：LESVie為生態(tài)服務價值分布指數(shù)；ESV，A，i，e與土地利用類型分布指數(shù)的計算公式含義相同。若LESVie>1，則為優(yōu)勢分布區(qū)，且LESVie值越大，在此梯度該類生態(tài)系統(tǒng)的服務價值越高。

表1　洞庭湖區(qū)生態(tài)用地單位面積生態(tài)服務價值當量因子

3　結果與分析

3.1　生態(tài)用地演化基本特征

3.1.1生態(tài)用地數(shù)量變化特征洞庭湖區(qū)1990—2015 年間生態(tài)用地由1990年的17 215.77 km2、占總用地的46.48%，增加到2015年的17 258.28 km2、占土地總量的46.59%。其中濕地數(shù)量增加，由6 192.45 km2增加至6 354.20 km2，所占土地總量的比例分別由16.72%，增加至17.16%；林地和草地數(shù)量下降，分別由10 588.69 km2、434.63 km2減少到10 514.01 km2、390.07 km2，所占土地總量的比例分別由28.59%，1.17%，減少到28.39%，1.05%。從生態(tài)用地數(shù)量的變化歷程來看，可歸納為3個時期：1990—2000年為穩(wěn)定期，各類生態(tài)用地數(shù)量波動較小，總量小幅增加；2000—2010年為增長期，各類生態(tài)用地數(shù)量大幅增加，占土地總量的比例一度達到48.41%；2010—2015年為萎縮期，生態(tài)用地數(shù)量劇烈下降，抵消了前20 a的大部分增量(圖1)。

3.1.1生態(tài)用地結構變化的地形梯度效應按以上公式計算出地形位指數(shù)T(0～3.250 5)(圖2)，整體呈現(xiàn)出中間低，四周高的空間分布規(guī)律。將地形位指數(shù)等分為50個梯度，并劃分為低梯度(1～10)、中低梯度(10～20)、中高梯度(20～30)、高梯度(30～50)。由于水受重力影響向地勢低洼處匯集的自然特征，濕地主要分布在低梯度。林地分布指數(shù)從中梯度開始明顯增加，從中高梯度開始呈現(xiàn)為優(yōu)勢分布區(qū)，且優(yōu)勢程度隨地形梯度遞增。草地主要分布在高梯度，且比例隨地形梯度升高遞增，在梯度50出現(xiàn)遠高于其他地形梯度的峰值，中梯度是各類生態(tài)用地分布都較少的區(qū)域，也是人類生產(chǎn)生活的聚集區(qū)。各時期不同地類的分布規(guī)律相近，故僅列出1990年、2015年為例。

圖1洞庭湖區(qū)1990-2015年生態(tài)用地數(shù)量變化

圖2洞庭湖區(qū)地形位指數(shù)空間分布

對比1990年、2015年各地形梯度上的生態(tài)用地分布指數(shù)的變化可以發(fā)現(xiàn)：中低(10～20)、中高(20～30)梯度濕地分布指數(shù)減小，說明中低、中高梯度是人類活動的主要區(qū)域，耕地和建設用地的擴張占用了濕地；其中梯度(4～9)濕地分布指數(shù)增加明顯，是濕地面積增長的主要區(qū)域，也是20世紀圍湖造田集中發(fā)生的區(qū)間。林地分布指數(shù)在中低、中高梯度(10～30)減小，在低梯度(1～10)和高梯度(30～50)增大，這是由中低、中高梯度的非生態(tài)用地擴張與低梯度的生態(tài)林建設共同導致的現(xiàn)象。草地總量較小，因此波動幅度明顯大于濕地和林地，其分布指數(shù)在低梯度降低而在高梯度上升，一方面中梯度的人類活動導致草地面積減少，另一方面是由于洞庭湖水位的上升，湖畔大量草地轉化為沼澤(圖3)。

圖31990年、2015年不同地形梯度上生態(tài)用地分布指數(shù)

3.2　生態(tài)用地轉移特征

3.2.1生態(tài)用地轉移基本特征洞庭湖區(qū)1990—2015年期間生態(tài)用地綜合動態(tài)度為1.94%，草地轉移速率最快，為每年—0.69%；濕地和林地則以年均0.18%，-0.05%的速率轉移。1990—2000年期間生態(tài)用地綜合動態(tài)度為2.01%，生態(tài)用地的數(shù)量變化并不明顯，轉入和轉出處于穩(wěn)定平衡的狀態(tài)，濕地和林地的轉入量略高于轉出量，只有草地的轉出規(guī)模較大，且去向以濕地為主。2000—2010年期間生態(tài)用地綜合動態(tài)度為3.24%，生態(tài)用地轉移速率明顯加快，總量以轉入為主，濕地數(shù)量的增長最為顯著，林地次之、草地增長面積較少，但增長比例較高。2010—2015年期間生態(tài)用地綜合動態(tài)度為3.57%，生態(tài)用地轉移速率進一步加快，綜合活躍程度更高，但生態(tài)用地的轉移方向變?yōu)橐赞D出為主，非生態(tài)用地擴張迅猛。

表2　洞庭湖區(qū)1990-2015生態(tài)用地轉移速率及綜合動態(tài)度　%

3.2.2生態(tài)用地及轉移的地形梯度分布特征由圖4可見，1990—2000年期間生態(tài)用地轉出和轉入基本平衡，轉出主要發(fā)生在中低和中高梯度(10～30)，峰值位于梯度18，是轉出比例最高的區(qū)域。生態(tài)用地轉入的優(yōu)勢分布區(qū)則主要位于中高梯度(20～30)，然而峰值卻出現(xiàn)在梯度3，濕地的大量轉入是造成這一現(xiàn)象的直接原因。這段時期雖然生態(tài)用地總量波動較小，但在垂直分布上向更高地形梯度轉移。高梯度的生態(tài)用地的轉出和轉入在高梯度的活躍程度逐級遞減。

2000—2010年期間生態(tài)用地轉移以轉入為主，轉出在中梯度(10～20)的分布指數(shù)較上一期明顯下降，于梯度3出現(xiàn)峰值，而優(yōu)勢分布區(qū)向中高梯度(20～30)集中，生態(tài)用地轉入的優(yōu)勢分布區(qū)則明顯偏向中低梯度(10～20)。隨著退田還湖工程的推進，中低梯度的生態(tài)用地數(shù)量明顯增加。2010—2015年期間生態(tài)用地轉出和轉入優(yōu)勢區(qū)的同步性較強，主要都活躍在中高梯度(10～30)，生態(tài)用地的保護與蠶食展開拉鋸，但從總量變化上看轉出明顯占優(yōu)，生態(tài)用地數(shù)量大幅下降，以中梯度損失最為嚴重，是人類活動較為集中的區(qū)域。1990—2015年期間生態(tài)用地轉出主要發(fā)生在中低和中高低梯度(10～30)，也是城鎮(zhèn)和農(nóng)業(yè)空間擴張的主要梯度，生態(tài)用地轉入主要集中于低梯度(3～10)，是退耕還湖的重點區(qū)域。梯度3在各時段生態(tài)用地的轉出和轉入都出現(xiàn)波峰，表明以濕地為主的梯度3是濕地和非生態(tài)用地相互轉移拉鋸最重要的梯度。

圖4　1990-2015年不同地形位上生態(tài)用地轉移分布指數(shù)

3.3　生態(tài)用地生態(tài)服務價值演化特征

3.3.1生態(tài)服務價值演化的基本特征洞庭湖區(qū)生態(tài)用地生態(tài)服務價值在1990—2015年期間由856.72億元增加到1 393.28億元，漲幅為62.63%。其中1990—2000年期間年均增長幅度為4.7%，受1988年、1998年兩次洪澇災害的沖擊，糧食減產(chǎn)嚴重，生態(tài)用地生態(tài)服務價值明顯低于其他年份。對比1990年和2000年的生態(tài)用地與生態(tài)服務價值變化，在洪水流量更大的情況下[29]，1998年受到的影響要明顯小于1988年，洞庭湖水文調節(jié)功能提高是重要因素。2000—2010年期間增長速率明顯放緩，但生態(tài)用地在各種生態(tài)建設工程的實施下大面積恢復，生態(tài)服務功能有明顯的實質性提升。2010—2015年期間生態(tài)用地受到快速城鎮(zhèn)化的影響快速減少，生態(tài)服務價值也有所下降。濕地和林地提供的生態(tài)服務價值較高，草地由于總量較小，提供的生態(tài)服務價值也較低。

3.3.2生態(tài)用地生態(tài)服務價值及變化的地形梯度分布特征由于濕地主要分布于低梯度，林地數(shù)量則隨梯度升高遞增，洞庭湖區(qū)生態(tài)用地生態(tài)服務價值的優(yōu)勢分布區(qū)集中在低梯度(1～10)和高梯度(30～50)。而低梯度生態(tài)服務價值明顯高于高梯度則與濕地的單位面積當量因子較高有明顯的關聯(lián)，是典型的大型湖泊地區(qū)的共有特征。相比供給服務和支持服務，文化服務和調節(jié)服務更集中在低梯度。生態(tài)服務價值在地形梯度8，13分別出現(xiàn)兩次波峰，與濕地分布指數(shù)在梯度8，13出現(xiàn)波峰響應，說明低梯度和中梯度的生態(tài)服務價值由濕地主導。在梯度(20～50)，生態(tài)服務價值的分布趨勢與林地極為接近，表明林地貢獻了高梯度絕大部分的生態(tài)服務價值。草地由于總量很小，在各地形梯度提供的生態(tài)服務價值也較低，對生態(tài)用地生態(tài)服務價值的整體分布沒有顯著影響。

表3　1990-2015年生態(tài)用地生態(tài)服務價值變化

1990—2000年期間生態(tài)服務價值總量大幅增加，生態(tài)用地由中低梯度(10～20)和中高梯度(20～30)向高梯度(30～50)轉移的趨勢明顯，生態(tài)服務價值的流動方向與之基本同步，并在梯度3與生態(tài)用地轉入產(chǎn)生波峰響應，總量增加的同時在垂直空間上完成了重新配置(圖5)。2000—2010年期間中低梯度(10～20)生態(tài)用地大量恢復，生態(tài)服務價值分布指數(shù)也出現(xiàn)明顯增加，并在在梯度10形成波峰，生態(tài)用地轉入明顯高于轉出，是生態(tài)用地恢復成效最明顯的區(qū)域。梯度8，13的濕地數(shù)量較大，但變化并不明顯，生態(tài)服務價值分布指數(shù)變化出現(xiàn)波谷，增長并不如附近其他梯度明顯。梯度3為波谷，與生態(tài)用地轉出的波峰響應。2010—2015生態(tài)服務價值總量下降明顯，各梯度生態(tài)用地都有不同程度的減少，但主要集中在中梯度，表現(xiàn)為中低梯度(10～20)是生態(tài)服務價值下降最明顯的區(qū)域，生態(tài)服務價值的大幅下降在梯度10形成波谷，與2000—2010年期間分布指數(shù)在梯度10的波峰形成鮮明對比，非生態(tài)用地在梯度10附近的分布優(yōu)勢極為顯著，這說明梯度10附近不僅是人類活動最為密集的地區(qū)，生態(tài)用地的流轉頻率也最高，穩(wěn)定性最差。在中低梯度生態(tài)服務價值整體下滑的情況下，梯度8，13仍然與濕地分布產(chǎn)生波峰響應，生態(tài)服務價值下降較少，證明濕地對洞庭湖區(qū)生態(tài)服務價值的分布具有決定性作用，也是其穩(wěn)定性的重要保障。

圖6　1990-2015年生態(tài)用地生態(tài)服務價值地形分布指數(shù)變化

1990—2015年期間，生態(tài)服務價值的增加主要集中在低梯度(1～10)，低梯度濕地的大規(guī)?；謴?，減輕了核心湖區(qū)的生態(tài)服務壓力，高梯度林地的良好維護與保養(yǎng)使洞庭湖區(qū)各地形梯度的生態(tài)服務價值分布更加平衡。然而這樣的變化主要來源于2000—2010年期間退湖還田和生態(tài)林地建設，2010年后生態(tài)用地銳減的趨勢導致洞庭湖區(qū)的生態(tài)安全形勢依舊嚴峻。

4　結 論

(1) 洞庭湖區(qū)生態(tài)用地數(shù)量在研究期間呈先增加后減少的變化趨勢，其中濕地數(shù)量明顯增加，林地、草地數(shù)量減少。轉移速率最高的為草地，濕地次之，林地最低。洞庭湖區(qū)在2000—2010年期間實施的退耕還林、還草、還湖等生態(tài)工程使各類生態(tài)用地數(shù)量都有顯著增加，綜合動態(tài)度達到3.24%，活躍程度遠高于1990—2000年期間的2.01%，而2010—2015年城鎮(zhèn)化建設的快速推進，生態(tài)用地大幅減少，綜合動態(tài)度達到3.57%。生態(tài)用地轉移主要表現(xiàn)為草地向非生態(tài)用地、林地轉移，林地向草地轉移，濕地向非生態(tài)用地轉移。非生態(tài)用地主要轉移為草地和林地。

(2) 生態(tài)用地生態(tài)服務價值在1990—2015年期間從856.72億元增加到1 393.28億元，漲幅為62.63%，農(nóng)田生態(tài)系統(tǒng)糧食生產(chǎn)的服務價值從1990年的1 336.54元/hm2提高到2015 年的2 150.02元/hm2是決定性因素，但在生態(tài)用地總量減少的情況下，濕地恢復也起到了至關重要的作用。

(3) 不同生態(tài)用地在地形梯度上表現(xiàn)出顯著差異。垂直分布特征大致保持三段模式：低梯度(1～10)為濕地的優(yōu)勢分布區(qū)；高梯度(30～50)是林地和草地的優(yōu)勢分布區(qū)。生態(tài)用地轉出和轉入在不同時期和地形梯度上交替進行。1990—2000年期間主要變現(xiàn)為生態(tài)用地總量變化不大，空間上向較高梯度轉移。2000—2010年各梯度生態(tài)用地數(shù)量都有明顯增加，以中低梯度為主，濕地恢復成效可觀。而2010—2015年期間非生態(tài)用地快速擴張，在中低梯度占據(jù)了數(shù)量龐大的生態(tài)用地。

(4) 1990—2015年期間低梯度濕地大面積增加，生態(tài)服務價值分布指數(shù)增長明顯，中梯度由于非生態(tài)用地的擴張導致分布指數(shù)的下降幅度大于高梯度區(qū)域。其中1990—2000年。2000—2010年期間梯度3濕地的大規(guī)?；謴蜆O大地緩解了梯度1，2湖泊核心區(qū)的生態(tài)服務壓力，濕地生態(tài)服務功能得到顯著增強。但2010年后，中梯度濕地和林地大量減少，生態(tài)服務價值也隨之降低，與2000—2010年以生態(tài)用地恢復為主形成強烈反差，之前恢復的生態(tài)用地遭到快速蠶食，相比低梯度濕地生態(tài)系統(tǒng)有較高的政策保障和生態(tài)建設投入，中高梯度的林地、草地的保護存在較大隱患。

洞庭湖區(qū)是長江流域中段最重要的水文調節(jié)緩沖區(qū)，同時也是兩湖地區(qū)重要的生物多樣性寶庫，對維持兩湖地區(qū)的生態(tài)平衡有著不可或缺的作用。通過對最近25 a生態(tài)用地格局演變與生態(tài)服務價值變化在不同地形梯度上差異化分布的研究，揭示其垂直空間分布規(guī)律，有助于環(huán)保機構和部門進行更精準的生態(tài)管控和建設，在生態(tài)環(huán)境遭受快速城鎮(zhèn)化高強度沖擊的新形勢下采取切實有效的措施積極應對。

參考文獻:

[1]喻鋒,李曉波,張麗君,等.中國生態(tài)用地研究:內涵、分類與時空格局[J].生態(tài)學報,2015,35(14):4931-4943.

[2]Leemans H B J, Groot R S D. Millennium ecosystem assessment:ecosystems and human well-being:A framework for assessment[J]. Physics Teacher, 2003,34(9):534-534.

[3]周朕,蒙吉軍,齊楊,等.中國生態(tài)用地重要性及其格局優(yōu)化研究進展[J].生態(tài)學雜志,2016,35(1):218-225.

[4]陳芳淼,田亦陳,袁超,等.基于供給生態(tài)服務價值的云南土地資源承載力評估方法研究[J].中國生態(tài)農(nóng)業(yè)學報,2015,23(12):1605-1613.

[5]康鵬,陳衛(wèi)平,王美娥.基于生態(tài)系統(tǒng)服務的生態(tài)風險評價研究進展[J].生態(tài)學報,2016,36(5):1192-1203.

[6]徐偉平,康文星,何介南.洞庭湖區(qū)生態(tài)系統(tǒng)服務功能價值分析[J].草業(yè)學報,2016,25(1):217-229.

[7]陳陽,張建軍,杜國明,等.三江平原北部生態(tài)系統(tǒng)服務價值的時空演變[J].生態(tài)學報,2015,35(18):6157-6164.

[8]姚小薇,曾杰,李旺君.武漢城市圈城鎮(zhèn)化與土地生態(tài)系統(tǒng)服務價值空間相關特征[J].農(nóng)業(yè)工程學報,2015,31(9):249-256.

[9]姜翠紅,李廣泳,程滔,等.青海湖流域生態(tài)服務價值時空格局變化及其影響因子研究[J].資源科學,2016,38(8):1572-1584.

[10]石培禮,李文華,何維明,等.川西天然林生態(tài)服務功能的經(jīng)濟價值[J].山地學報,2002,20(1):75-79.

[11]劉興元,龍瑞軍,尚占環(huán).草地生態(tài)系統(tǒng)服務功能及其價值評估方法研究[J].草業(yè)學報,2011,20(1):167-174.

[12]孫新章,謝高地,成升魁,等.中國農(nóng)田生產(chǎn)系統(tǒng)土壤保持功能及其經(jīng)濟價值[J].水土保持學報,2005,19(4):156-159.

[13]張靜,任志遠.秦巴山區(qū)土地利用時空格局及地形梯度效應[J].農(nóng)業(yè)工程學報,2016,32(14):250-257.

[14]劉欣,趙艷霞,武愛彬,等.河北省淺山丘陵區(qū)土地利用變化及人文驅動分析[J].農(nóng)業(yè)現(xiàn)代化研究,2015,36(1):126-131.

[15]郜紅娟,韓會慶,俞洪燕,等.烏江流域重要生態(tài)系統(tǒng)服務地形梯度分布特征分析[J].生態(tài)科學,2016,35(5):154-159.

[16]岳健,張雪梅.關于我國土地利用分類問題的討論[J].干旱區(qū)地理,2003,26(1):78-88.

[17]龍花樓,劉永強,李婷婷,等.生態(tài)用地分類初步研究[J].生態(tài)環(huán)境學報,2015,24(1):1-7.

[18]鄧紅兵,陳春娣,劉昕,等.區(qū)域生態(tài)用地的概念及分類[J].生態(tài)學報,2009,29(3):1519-1524.

[19]王琦,湯放華.洞庭湖區(qū)生態(tài)—經(jīng)濟—社會系統(tǒng)耦合協(xié)調發(fā)展的時空分異[J].經(jīng)濟地理,2015,35(12):161-167.

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

[21]劉紀遠,匡文慧,張增祥,等.20世紀80年代末以來中國土地利用變化的基本特征與空間格局[J].地理學報,2014,69(1):3-14.

[22]孫丕苓,許月卿,王數(shù).環(huán)京津貧困帶土地利用變化的地形梯度效應分析[J].農(nóng)業(yè)工程學報,2014,30(14):277-288.

[23]喻紅,曾輝,江子灜.快速城市化地區(qū)景觀組分在地形梯度上分布特征研究[J].地理科學,2001,21(1):64-68.

[24]Costanza R, Arge R, Groot R, et al. The value of the world′s ecosystem services and nature capital[J]. Nature,1997,387(15):253-260.

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

[26]謝高地,張彩霞,張雷明,等.基于單位數(shù)量價值當量因子的生態(tài)系統(tǒng)服務價值化方法改進[J].自然資源學報,2015,30(8):1243-1253.

[27]武愛彬,趙艷霞.壩上高原生態(tài)用地時空格局演變與生態(tài)系統(tǒng)服務價值分析[J].農(nóng)業(yè)工程學報,2017,33(2):283-290.

[28]國家發(fā)展和改革委員會價格司.全國農(nóng)產(chǎn)品成本收益匯編(2015)[M].北京:中國統(tǒng)計出版社,2015.

[29]高俊峰,張琛,姜加虎,等.洞庭湖的沖淤變化和空間分布[J].地理學報,2001,56(3):269-277.