吳愛林, 陳 燕, 燕彩霞, 羅 嬋, 宋新山

(東華大學(xué) 環(huán)境學(xué)院, 上海 201620)

長江三角洲生態(tài)系統(tǒng)服務(wù)價值分析及趨勢預(yù)測

吳愛林, 陳 燕, 燕彩霞, 羅 嬋, 宋新山

(東華大學(xué)環(huán)境學(xué)院,上海201620)

[目的] 對長江三角洲生態(tài)系統(tǒng)服務(wù)價值進(jìn)行分析及趨勢預(yù)測，以期為該區(qū)未來城市發(fā)展以及制定生態(tài)環(huán)境建設(shè)補(bǔ)償政策提供支持。 [結(jié)果] 以長江三角洲16個地市2000—2014年的土地利用數(shù)據(jù)和相關(guān)統(tǒng)計資料為數(shù)據(jù)源，采用Costanza生態(tài)系統(tǒng)服務(wù)價值計算方法，參照修正后的生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價值系數(shù)，結(jié)合敏感性分析，探討生態(tài)系統(tǒng)服務(wù)價值的時空演化規(guī)律，并利用馬爾科夫模型對其未來3 a的發(fā)展趨勢進(jìn)行了預(yù)測分析。 [結(jié)果] 區(qū)域生態(tài)系統(tǒng)總服務(wù)價值從2000年的1 303.89億元減少到2014年的1 294.25億元，變化率為-0.74%，水源涵養(yǎng)、廢物處理和娛樂文化生態(tài)系統(tǒng)服務(wù)功能價值上升；在空間上，常州、南京、寧波、泰州、臺州、揚(yáng)州的生態(tài)系統(tǒng)服務(wù)價值不斷增加，而其他城市的則處于下降趨勢；對長江三角洲2015和2017年的土地利用構(gòu)成和生態(tài)系統(tǒng)服務(wù)價值進(jìn)行預(yù)測，呈現(xiàn)耕地、林地持續(xù)減少，建設(shè)用地持續(xù)增加的趨勢，生態(tài)系統(tǒng)總服務(wù)價值將減少276.65億元。 [結(jié)論] 在未來土地調(diào)整中應(yīng)更多地考慮到生態(tài)服務(wù)價值的重要意義，格外注意對耕地、水域的保護(hù)，控制建設(shè)用地占用耕地和林地，繼續(xù)加大植樹造林力度，提高土地利用集約度，維持生態(tài)系統(tǒng)的功能。

長江三角洲; 土地利用; 生態(tài)系統(tǒng)服務(wù)價值; 敏感性分析; 馬爾科夫模型

文獻(xiàn)參數(shù)： 吳愛林, 陳燕, 燕彩霞, 等.長江三角洲生態(tài)系統(tǒng)服務(wù)價值分析及趨勢預(yù)測[J].水土保持通報，2017,37(4)：254-259.DOI:10.13961/j.cnki.stbctb.2017.04.043； Wu Ailin, Chen Yan, Yancai Xia, et al. Ecosystem service value analysis and trend prediction in Yangtze River delta[J]. Bulletin of Soil and Water Conservation, 2017,37(4)：254-259.DOI:10.13961/j.cnki.stbctb.2017.04.043

生態(tài)系統(tǒng)服務(wù)功能，是生態(tài)系統(tǒng)、生態(tài)過程形成的自然環(huán)境條件及其效用[1]，是維持人類生存的根本。近幾年，生態(tài)服務(wù)價值已經(jīng)成為可持續(xù)發(fā)展研究的熱點，普遍研究者已將生態(tài)服務(wù)價值作為評價生態(tài)環(huán)境變化的主要指標(biāo)之一[2-3]，生態(tài)系統(tǒng)服務(wù)價值(ecosystem service value， ESV)的定量評估已成為生態(tài)經(jīng)濟(jì)學(xué)和環(huán)境經(jīng)濟(jì)學(xué)的研究重點[4-5]。對生態(tài)系統(tǒng)服務(wù)價值及其變化進(jìn)行評估，不僅反映了人類對自然資源周期性利用情況，更重要是突出在經(jīng)濟(jì)發(fā)展過程中人類對生態(tài)系統(tǒng)平衡的重要性的認(rèn)知[6]。生態(tài)系統(tǒng)服務(wù)價值是指通過生態(tài)系統(tǒng)的結(jié)構(gòu)、過程和功能直接或間接提供的生命支持產(chǎn)品(如原材料和食物等)和服務(wù)(如棲息地提供等)，包括人類生活所必需的生態(tài)產(chǎn)品和保證人類生活質(zhì)量的生態(tài)功能兩部分[7]。國內(nèi)外學(xué)者已經(jīng)分別從全球、國家、區(qū)域等不同尺度以及單個生態(tài)系統(tǒng)服務(wù)價值的評估進(jìn)行了研究，其中Costanza等發(fā)表的“全球生態(tài)系統(tǒng)服價值和自然資本”對之后的相關(guān)研究產(chǎn)生了深遠(yuǎn)的影響，使生態(tài)系統(tǒng)服務(wù)價值評估的原理與方法從科學(xué)意義上得以明確。本文基于ENVI軟件平臺，采用轉(zhuǎn)移矩陣、Costanza等[1]提出的生態(tài)系統(tǒng)服務(wù)價值計算公式和馬爾科夫預(yù)測方法，結(jié)合敏感度分析，探討了2000—2014年15 a間不同空間尺度上的生態(tài)系統(tǒng)服務(wù)價值變化，并對2015—2025年的生態(tài)系統(tǒng)服務(wù)價值進(jìn)行預(yù)測，以期為長江三角洲未來城市發(fā)展以及制定生態(tài)環(huán)境建設(shè)補(bǔ)償政策提供支持。

1 研究區(qū)概況

長江三角洲地區(qū)包括杭州、蘇州、臺州、紹興、寧波、上海、南通、揚(yáng)州、湖州、南京、無錫、泰州、常州、鎮(zhèn)江、嘉興和舟山兩省一市共16個城市。是中國綜合實力最強(qiáng)的經(jīng)濟(jì)中心。截至2014年末，長江三角洲國內(nèi)生產(chǎn)總值是2000年的6.62倍，達(dá)到了128 802.8億元，占全國的20.25%。經(jīng)濟(jì)的快速發(fā)展，城鎮(zhèn)化的快速推進(jìn)，城市用地類型及空間結(jié)構(gòu)也隨之發(fā)生變化。研究期間，耕地、林地、草地面積減少，其余類型面積呈增加趨勢。其中，耕地的變化幅度最大，面積減少了8 129.45 km2，其次為城鎮(zhèn)用地，城鎮(zhèn)用地凈增加了5 162.18 km2，從而嚴(yán)重影響了區(qū)域內(nèi)的生態(tài)系統(tǒng)服務(wù)價值。

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

2.1 研究方法

2.1.1 生態(tài)系統(tǒng)服務(wù)價值評價方法 本研究以謝高地等[8]提出的“中國生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價值當(dāng)量”為基礎(chǔ)，結(jié)合長三角的實際情況進(jìn)行修訂，根據(jù)長三角2000—2014年統(tǒng)計年鑒的相關(guān)數(shù)據(jù)，計算得到長三角2000—2014年的平均糧食產(chǎn)量為2 087.42 kg/hm2，2014年間的平均糧食價格為2.604 6元/kg，參考謝高地[8]關(guān)于單位生態(tài)服務(wù)價值當(dāng)量因子的經(jīng)濟(jì)價值量等于當(dāng)年全國平均糧食產(chǎn)量市場價值的1/7的方法，計算出長三角一個當(dāng)量因子的經(jīng)濟(jì)價值約為776.82元/hm2。由此計算出長三角不同生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價值(表1)。結(jié)合相關(guān)土地利用類型面積，根據(jù)表1對長三角生態(tài)系統(tǒng)服務(wù)價值進(jìn)行評估。其中耕地和農(nóng)田對應(yīng)，林地和森林對應(yīng)，未利用土地和荒漠對應(yīng)，城鎮(zhèn)用地、農(nóng)村居民點和其他建設(shè)用地賦值為零。

(1)

(2)

式中：ESV， ESVf——生態(tài)系統(tǒng)服務(wù)價值和第f項服務(wù)價值；Ai——土地利用類型i的面積； VCfi——單位面積土地利用類型i的第f項服務(wù)功能價值系數(shù)； VCi——生態(tài)系統(tǒng)價值系數(shù)。

2.1.2 敏感性分析方法 生態(tài)系統(tǒng)價值敏感性指數(shù)(coefficient of sensitivity， CS)是為了驗證生態(tài)系統(tǒng)類型對于土地利用類型的代表性以及生態(tài)價值系數(shù)的準(zhǔn)確性而引入的[9-10]。該指數(shù)主要是用來確定ESV對VC的依賴程度，當(dāng)CS<1，表明ESV對于VC是缺乏彈性的，當(dāng)CS>1，表明ESV對于VC是富有彈性的。CS越大，表明VC的準(zhǔn)確性越關(guān)鍵。本文通過對各土地利用類型的VC進(jìn)行上下調(diào)整50%來說明ESV的變化情況。其計算公式為：

(3)

式中：b,a——初始價值和生態(tài)價值系數(shù)調(diào)整后的價值。

表1 長江三角洲地區(qū)生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價值

2.1.3 馬爾科夫預(yù)測方法 馬爾科夫預(yù)測過程是一種具有“無后效應(yīng)”的特殊隨機(jī)過程[11]，它表示隨機(jī)過程X(t)時刻(n+(1)的狀態(tài)概率分布只與時刻n的狀態(tài)概率分布有關(guān)，而與時刻n之前的狀態(tài)概率分布無關(guān)[12]?？捎棉D(zhuǎn)移概率矩陣P和初始分布S(0)并利用如下公式對土地利用狀態(tài)進(jìn)行預(yù)測[13-16]：

(4)

式中：Pij——i類土地利用類型轉(zhuǎn)化為j類土地利用類型的概率，Pij滿足2個條件：

平穩(wěn)馬爾柯夫鏈在時刻T+1的狀態(tài)采用下式計算

S(T+(1)=S(T)P

(5)

2.2 數(shù)據(jù)來源

以長江三角洲為研究區(qū)域，以2000和2014年Landsat TM衛(wèi)星影像為主要數(shù)據(jù)源，數(shù)據(jù)預(yù)處理包括數(shù)據(jù)讀取、輻射定標(biāo)、大氣校正(FLAASH模型)、波段合成、影像鑲嵌、研究區(qū)裁剪等。將長江三角洲的土地利用劃分為6大類17小類，分別是耕地(包括水田和旱地)；林地(包括有林地、灌木林、疏林地和其他林地)；草地(包括高覆蓋度草地、中覆蓋度草地和地覆蓋度草地)、水域(包括河渠、湖泊、水庫坑塘、灘涂和灘地)；城鄉(xiāng)、工礦、居民用地(包括城鎮(zhèn)用地、農(nóng)村居民點和其他建設(shè)用地)；未利用土地。解譯精度約為89.4%，并以當(dāng)年的植被圖、地形圖、地貌圖、等專題圖為依據(jù)進(jìn)行數(shù)據(jù)校驗。其他數(shù)據(jù)包括長江三角行政邊界、社會經(jīng)濟(jì)數(shù)據(jù)等。圖像處理及數(shù)據(jù)分析軟件分采用ArcGIS 10.2，ENVI5.1，Matlab等。

3 結(jié)果與分析

3.1 生態(tài)系統(tǒng)服務(wù)價值計算

根據(jù)公式(1)，可計算得到2000和2014年長三角生態(tài)系統(tǒng)服務(wù)總價值量(表2)。結(jié)果表明：長三角15 a生態(tài)系統(tǒng)服務(wù)總價值減少了9.65億元，下降了約0.74%，2000年的生態(tài)系統(tǒng)服務(wù)總價值為1 303.89億元，2014年為1 294.25億元。其中，林地對生態(tài)系統(tǒng)服務(wù)價值總量的貢獻(xiàn)最大，分別為41.15%和41.01%，而未利用地的貢獻(xiàn)最小。從各類用地生態(tài)系統(tǒng)服務(wù)價值的總體變化情況來看，水域和未利用土地的生態(tài)系統(tǒng)服務(wù)價值呈增長趨勢，水域的增加量明顯，占增加總量的93.95%；耕地、林地和草地生態(tài)系統(tǒng)服務(wù)價值呈降低趨勢，以耕地的降低量最大，占三者降低總量的86.56%，其次是林地，草地最低。研究期內(nèi)，水域生態(tài)系統(tǒng)服務(wù)價值占長江三角生態(tài)系統(tǒng)服務(wù)總價值的比例增加了3.4%，其生態(tài)系統(tǒng)服務(wù)價值所占比例均達(dá)30%以上，但其面積比例卻遠(yuǎn)小于此，2000，2014年長江三角洲水域面積比例分別為11.29%和12.18%。表明生態(tài)價值系數(shù)較高的水域面積增加已經(jīng)彌補(bǔ)不了由于耕地和林地減少引起的生態(tài)系統(tǒng)服務(wù)總價值的下降，使得長江三角洲生態(tài)系統(tǒng)服務(wù)總價值近年來有所下降。長江三角洲在今后土地利用中應(yīng)注重對林地和耕地生態(tài)系統(tǒng)的保護(hù)，繼續(xù)加強(qiáng)生態(tài)環(huán)境建設(shè)。

表2 長江三角洲地區(qū)生態(tài)系統(tǒng)服務(wù)價值(ESV)變化

3.2 單項生態(tài)系統(tǒng)服務(wù)總價值計算

根據(jù)公式(2)，可計算得到2000和2014年長三角單項生態(tài)系統(tǒng)服務(wù)價值量(表3)?？梢钥闯鲅芯科趦?nèi)生態(tài)系統(tǒng)各項服務(wù)價值的結(jié)構(gòu)順序為:水源涵養(yǎng)>廢物處理>土壤形成與保護(hù)>生物多樣性保護(hù)>氣候調(diào)節(jié)>氣體調(diào)節(jié)>娛樂文化>原材料>食物生產(chǎn)。前4項生態(tài)系統(tǒng)服務(wù)功能都占據(jù)各時期內(nèi)所有功能的10%以上，并且總和也超過了生態(tài)系統(tǒng)服務(wù)功能的60%，說明水源涵養(yǎng)、廢物處理、土壤形成與保護(hù)、生物多樣性保護(hù)是長三角地最主要的生態(tài)系統(tǒng)服務(wù)功能。水源涵養(yǎng)、廢物處理和娛樂文化呈不斷增加的趨勢，尤其是水源涵養(yǎng)功能增加程度最顯著，達(dá)到13.3億元，占增加量的60.56%。這與長三角水系發(fā)達(dá)、水量充沛密切相關(guān)。土壤形成與保護(hù)價值損失最大(10.5億元)，氣候調(diào)節(jié)次之(6.05億元)，其損失均主要由林地和高覆蓋地草地退化所致。

表3 長江三角洲地區(qū)生態(tài)系統(tǒng)單項服務(wù)功能價值(ESV)變化

3.3 生態(tài)系統(tǒng)服務(wù)價值空間分異變化

以長三角16地級市為單位，對生態(tài)系統(tǒng)服務(wù)價值變化量和變化率進(jìn)行計算和比較，結(jié)果詳見表4?？梢钥闯觯?000—2014年14 a間，杭州、湖州等10個市的生態(tài)服務(wù)價值處于下降趨勢，其中杭州市減少最多，達(dá)到5.57億元，其次是上海、紹興、無錫等市。其他地市生態(tài)系統(tǒng)服務(wù)價值均有不同程度的增加。其中，寧波市的增加量最大，為9.92億元，占全省生態(tài)系統(tǒng)服務(wù)價值總增加量的53.94%，年均增加0.71億元，其次是臺州市，增加了4.97億元，占全省生態(tài)系統(tǒng)服務(wù)價值總增加量的27.03%，年均增加0.36億元，然后依次為泰州市、南京市、揚(yáng)州市，常州市的增加量最少，為0.02億元從變化幅度來看，生態(tài)系統(tǒng)服務(wù)價值變化率下降最快的是嘉興市，達(dá)到了-9.07%，其次是上海市、無錫市，分別達(dá)到了-5.46%，-5.17%；其余地級市的變化幅度均低于±5%。

3.4 敏感性分析

根據(jù)式(3)，通過對各土地利用類型的VC進(jìn)行上下調(diào)整50%的生態(tài)系統(tǒng)價值系數(shù)來計算敏感性指數(shù)，結(jié)果詳見表5。由表5可以看出：ESV對VC的敏感性指數(shù)都小于1，由高到低依次為林地、水域、耕地、草地、未利用土地，其中，草地敏感性指數(shù)約為0.01，說明當(dāng)草地的VC增加1%時，ESV將增加0.01%。這表明，研究區(qū)域ESV對VC缺乏彈性，研究結(jié)果是可信的，可為長三角土地利用規(guī)劃提供借鑒和參考。

3.5 生態(tài)系統(tǒng)服務(wù)價值的馬爾科夫預(yù)測

運(yùn)用馬爾科夫過程進(jìn)行預(yù)測的關(guān)鍵在于轉(zhuǎn)移概率的確定。本研究先用ENVI 5.1軟件中的Class Statistics工具統(tǒng)計長三角2014年土地利用類型面積，形成初始狀態(tài)矩陣，再用ENVI5.1軟件對長三角2000和2014年的遙感影像分類結(jié)果圖進(jìn)行轉(zhuǎn)移矩陣計算，獲得了各種土地利用類型之間的轉(zhuǎn)化關(guān)系(表6)。由表6可以看出，2014年長江三角洲城鎮(zhèn)用地面積共增加了5 345.60 km2。有4 166.05 km2來自耕地，92.52 km2來自林地，11.26 km2來自草地，134.35 km2來自水域，826.60 km2來自農(nóng)村居民地，還有112.76，2.06 km2分別來自其他建設(shè)用地和未利用土地。同樣了解到林地、草地、水域等土地利用類型的轉(zhuǎn)變及其構(gòu)成。其中2014年耕地面積減少最多，主要轉(zhuǎn)化為城鎮(zhèn)用地和農(nóng)村居民用地。利用已求出的2000—2014年內(nèi)的土地利用各類型面積轉(zhuǎn)移矩陣，求出該時間段內(nèi)各個土地利用類型的年平均轉(zhuǎn)移概率，構(gòu)成1個年平均轉(zhuǎn)移概率矩陣，該矩陣即為初始狀態(tài)轉(zhuǎn)移矩陣。利用P(n)=Pn(n=1,2,3)求得最終轉(zhuǎn)移概率矩陣。最后用初始狀態(tài)矩陣與轉(zhuǎn)移概率矩陣相乘(運(yùn)算過程采用Matlab軟件進(jìn)行)，即可得最終狀態(tài)矩陣，即2015，2020，2025年預(yù)測的長三角各土地類型面積，并求得未來10 a生態(tài)系統(tǒng)總服務(wù)價值。

表4 長江三角洲地區(qū)各市生態(tài)系統(tǒng)服務(wù)價值空間分異

表5 調(diào)整生態(tài)價值系數(shù)后生態(tài)系統(tǒng)服務(wù)總價值(ESV)的變化情況(VC)及敏感性指數(shù)(CS)

表6 長江三角洲地區(qū)2000-2014年土地利用各類型轉(zhuǎn)化情況 km2

根據(jù)表7對長三角2015-2017年各類土地利用面積變化進(jìn)行預(yù)測，可以看出城鎮(zhèn)用地、農(nóng)村居民點、未利用土地面積出現(xiàn)不同程度增加，變化率分別達(dá)到42.87%，15.20%，21.02%。這表明隨著經(jīng)濟(jì)高速發(fā)展、城市化進(jìn)程加快，長江三角洲建設(shè)用地占用土地的強(qiáng)度還將不斷加大。

由預(yù)測結(jié)果可以看出，長江三角洲2015—2017年總的生態(tài)系統(tǒng)服務(wù)價值呈下降的趨勢，從1 223.5億元減少至946.85億元，共減少276.65億元。長江三角洲不同土地利用類型的生態(tài)系統(tǒng)服務(wù)價值變化規(guī)律略有不同:2015—2017年內(nèi)，僅有未利用土地的生態(tài)系統(tǒng)服務(wù)價值不變，其他的均出現(xiàn)下降。其中，耕地的生態(tài)系統(tǒng)服務(wù)價值量減少最多，達(dá)到46.20億元，占減少總量的71.74%；其次是林地，減少量達(dá)到11.95億元，占減少總量的18.56%；水域、草地的生態(tài)系統(tǒng)服務(wù)價值分別減少了5.03,1.22億元，所占土地生態(tài)服務(wù)價值減少總量的比例分別為7.81%和1.89%。

表7 研究區(qū)2015,2020，2025年各類土地利用面積及生態(tài)系統(tǒng)服務(wù)價值預(yù)測

4 結(jié) 論

2000—2014年長三角土地生態(tài)系統(tǒng)服務(wù)價值總體呈不斷下降的時序變化趨勢，年下降率為0.74%。這與長三角耕地、林地、草地面積的大幅度減少，引起土壤形成與保護(hù)、食物生產(chǎn)、氣候調(diào)節(jié)等價值快速下降密切相關(guān)。即使高生態(tài)價值的水域面積增加也彌補(bǔ)不了由于耕地、林地、草地減少引起的生態(tài)系統(tǒng)服務(wù)總價值的下降。在空間上，杭州市、蘇州市、臺州市的生態(tài)系統(tǒng)服務(wù)價值在各地級市排序中排前3位，鎮(zhèn)江市、嘉興市、舟山市的生態(tài)系統(tǒng)服務(wù)價值最低，尤其是嘉興市在研究期間生態(tài)系統(tǒng)服務(wù)價值的下降率最高，迫切需要保護(hù)耕地、水域，繼續(xù)加大植樹造林力度，維持生態(tài)系統(tǒng)的功能。

根據(jù)馬爾柯夫過程對2015—2017年長三角的土地利用結(jié)構(gòu)和生態(tài)系統(tǒng)服務(wù)價值進(jìn)行了預(yù)測。未來土地利用生態(tài)系統(tǒng)服務(wù)價值變化趨勢與研究區(qū)主要構(gòu)成者耕地、林地、水域的變化趨勢保持一致，研究區(qū)域生態(tài)系統(tǒng)總服務(wù)價值將減少276.65億元。因此，在未來土地調(diào)整中應(yīng)更多地考慮到生態(tài)服務(wù)價值的重要意義，格外注意對森林生態(tài)系統(tǒng)的保護(hù)，控制建設(shè)用地占用耕地和林地，提高土地利用集約度。

[1] Daily G C, Myers J P, Reichert J. Nature’s service:societal dependence on natural ecosystems[M].Washington: Island Press, 1997.

[2] Seppelt R, Dormann C F, Eppink F V, et al. A quantitative review of ecosystem service studies: Approaches,shortcomings and the road ahead [J].Journal of Applied Ecology, 2011,48(3):630-636.

[3] Fisher B, Turner R K, Morling P. Defining and classifying ecosystem services for decision making [J].Ecological Economics, 2009,68(3):643-653.

[4] Chen Zhongxin, Zhang Xinshi. Value of ecosystem services in China[J]. Chinese Science Bulletin, 2000,45(10):870-876.

[5] Pearce D. Auditing the Earth: The value of the world’s ecosystem services and natural capital[J]. Environment, 1998,40(2):23-28.

[6] 黃湘,陳亞寧,馬建新.西北干旱區(qū)典型流域生態(tài)系統(tǒng)服務(wù)價值變化[J].自然資源學(xué)報,2011,26(8):1364-1376.

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

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

[9] Li Renqiang, Dong Ming, Cui Jianyong, et al. Quantification of the impact of land-use changes on ecosystem services: A case study in Pingbian County, China[J]. Environmental Monitoring & Assessment, 2007, 128(4):503-510.

[10] 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.

[11] Hulst R. On the dynamics of vefetation: Markovchains as models of succession[J]. Vegetation, 1979,40(1):3-14.

[12] 胡汝曉,趙松義,譚周進(jìn),等.煙草連作對稻田土壤微生物及酶的影響[J].核農(nóng)學(xué)報,2007,21(5):494-497.

[13] 王曉峰,任志遠(yuǎn),黃青.農(nóng)牧交錯區(qū)縣城土地利用變化及驅(qū)動力分析:以陜北神木縣為例[J].干旱區(qū)地理,2003,26(4):402-407.

[14] 郜魯豪,徐旌.基于Markov模型的安寧市土地利用預(yù)測[J].云南地理環(huán)境研究,2010,22(2):87-89.

[15] 藍(lán)永超,丁永建,康爾酒,等.黑河流域水資源動態(tài)變化及其趨勢的灰色Markov鏈預(yù)測[J].中國沙漠,2003,24(3):435-440.

[16] 李黔湘,王華斌.基于馬爾科夫模型的漲渡湖流域土地利用變化預(yù)測[J].資源科學(xué),2008,30(10):1542-1543.

Ecosystem Service Value Analysis and Trend Prediction in Yangtze River Delta

WU Ailin, CHEN Yan, Yancai Xia, LUO Chan, SONG Xinshan

(College of Environment, Donghua University, Shanghai 201620, China)

[Objective] This paper examined the spatio-temporal variations in the value of ecosystem services and then predicted it to provide support for formulating compensation policy in ecological environment construction and future urban development. [Methods] Based on land use data and relevant statistical data in the Yangtze River delta from 2000 to 2014, with the method proposed by Costanza wherein the values of ecosystem services are calculated with regard to the value coefficients of ecosystem services per unit area in the Chinese terrestrial ecosystem and the sensitivity analysis, this paper examined the spatio-temporal variations in the value of ecosystem services and then predicted it by Markov model. [Results] The total value of ecosystem services decreased from 130 389 million yuan in 2000 to 129 425 million yuan, reduced by 0.74%. In addition, the ecosystem service value of water supply, recreation and culture function and waste treatment increased while other ecosystem functions declined. From the perspective of spatial analysis, the ecosystem service values of Changzhou, Nanjing, Ninbo, Taizhou, Taizhou, Yangzhou City increased, while values of other cities decreased. The land use structures of the Yangtze River delta in 2015 and 2017 were predicted using the Markov model. The trends of woodland and cultivated land will maintain a continual increase while the construction land will be in a continual decrease in the coming three years. The total ecosystem services value will have a decrement of 276.65 million yuan. [Conclusion] Therefore, the importance of ecological service value should be given more consideration in future land adjustment, especially for the protection of cultivated land and water and improvement of land utilization intensity. it is very important to control the occupancies of cultivated and forest land by construction, to continually increase the effort of afforestation and maintain the function of ecosystem.

YangtzeRiverDelta;landuse;ecosystemservicevalue;sensitivityanalysis;Markovmodel

A

： 1000-288X(2017)04-0254-06

： F301.24， Q149

2016-10-15

：2016-12-22

國家自然科學(xué)基金項目“長江三角洲區(qū)域生態(tài)安全時空演變及調(diào)控圖譜研究”(41471089); 中央高?；究蒲袠I(yè)務(wù)費(fèi)項目

吳愛林(1990—),女(漢族),江西省玉山縣人,碩士研究生,研究方向為長江三角洲的土地生態(tài)敏感性評價。E-mail:270134083@qq.com。

陳燕(1983—),女(漢族),上海市人，博士，副教授，主要從事生態(tài)安全、納米技術(shù)(材料)在水處理中的應(yīng)用以及環(huán)境遙感與環(huán)境信息系統(tǒng)。E-mail：chenyanhjxy@dhu.edu.cn。