蘇寧，范貴娟

（1.山東師范大學(xué)地理與環(huán)境學(xué)院，山東濟(jì)南250358；2.安徽農(nóng)業(yè)大學(xué)，安徽合肥230036）

基于土地利用變化的沂河流域生態(tài)服務(wù)價(jià)值研究

蘇寧1，范貴娟2

（1.山東師范大學(xué)地理與環(huán)境學(xué)院，山東濟(jì)南250358；2.安徽農(nóng)業(yè)大學(xué)，安徽合肥230036）

以沂河流域?yàn)槔?，在ARCGIS10.2中解譯沂河流域2005年和2015年遙感影像，得到兩期土地利用數(shù)據(jù)。并分析了2005—2015年沂河流域土地利用及其變化情況，以Costanza生態(tài)服務(wù)價(jià)值計(jì)算方法為參考，分析了沂河流域生態(tài)服務(wù)價(jià)值變化情況。結(jié)果表明：研究區(qū)耕地面積顯著減少，但生態(tài)服務(wù)總價(jià)值卻從2005年的2 338 964萬(wàn)元上升到2015年的2 427 039萬(wàn)元，上升了約3.77%，這主要是由于林地面積的增加帶來(lái)的較高的生態(tài)服務(wù)價(jià)值。

沂河流域；土地利用變化；生態(tài)系統(tǒng)服務(wù)價(jià)值

目前，對(duì)于生態(tài)服務(wù)價(jià)值的研究已成為評(píng)估可持續(xù)發(fā)展的關(guān)鍵環(huán)節(jié)[1]。生態(tài)服務(wù)這一概念指的是人類(lèi)從生態(tài)系統(tǒng)中直接或間接得到的利益[2]。關(guān)于生態(tài)系統(tǒng)的服務(wù)價(jià)值研究，有眾多學(xué)者對(duì)不同尺度上的生態(tài)系統(tǒng)服務(wù)價(jià)值定量評(píng)估并取得大量研究成果。Costanza等將全球各類(lèi)型土地利用基于對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行分類(lèi)及評(píng)估的基礎(chǔ)上進(jìn)行了生態(tài)服務(wù)價(jià)值計(jì)算[2]。歐陽(yáng)志云在論述生態(tài)系統(tǒng)服務(wù)的基礎(chǔ)上，結(jié)合我國(guó)實(shí)際，評(píng)估了中國(guó)陸地生態(tài)系統(tǒng)六種生態(tài)服務(wù)功能[3]。土地利用中生態(tài)服務(wù)功能表現(xiàn)為不同土地利用類(lèi)型的生態(tài)價(jià)值，具體來(lái)說(shuō)是人類(lèi)為滿(mǎn)足自身的生活及生產(chǎn)的需要而對(duì)土地利用的方式進(jìn)行改變，這種改變會(huì)影響到土地的生態(tài)服務(wù)價(jià)值[3-4]。謝高地等學(xué)者建立起中國(guó)不同陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值的計(jì)算方式，并計(jì)算出了我國(guó)單位面積土地上不同類(lèi)型土地的生態(tài)服務(wù)價(jià)值總和[5]。

作為自然界中受人類(lèi)干擾最大的部分，土地的利用及覆被的變化很大，進(jìn)而影響了生態(tài)系統(tǒng)的結(jié)構(gòu)與功能，具體而言表現(xiàn)為對(duì)區(qū)域水文、氣候、生物多樣性及生物地球化學(xué)循環(huán)等產(chǎn)生重大影響[6-9]。土地利用及覆被變化過(guò)程在維持生態(tài)系統(tǒng)服務(wù)功能上起是決定性的作用，因?yàn)槠鋵?duì)區(qū)域生態(tài)系統(tǒng)功能以及生態(tài)環(huán)境質(zhì)量影響深刻[10]。目前，鞏燦娟[11]研究了沂河流域中上游地區(qū)土地利用/覆被變化水文效應(yīng)的分析與模擬；陸建宇[12]等研究了沂河流域天然徑流變化規(guī)律；宋曉猛[13]等研究了沂河流域徑流演變特征及動(dòng)因分析；徐智廷[14]等研究了沂河流域地理特征及人文特點(diǎn)。

本研究通過(guò)沂河流域2005和2015年兩期遙感影像解譯出沂河流域2005—2015年的土地利用變化，通過(guò)研究方法計(jì)算出沂河流域研究期內(nèi)的生態(tài)服務(wù)價(jià)值變化，進(jìn)而分析出該流域土地利用中存在的問(wèn)題，并針對(duì)問(wèn)題提出沂河流域可持續(xù)發(fā)展的建議。

1 研究區(qū)概況

沂河流域?yàn)橐屎恿饔蚺R沂站以上河道流域，大致位于117°25′E～119°49′E，33°30′N(xiāo)～36°20′N(xiāo)之間，流域面積10 025 km2，流域涵蓋沂源縣大部分、沂水縣面積的60%、沂南縣大部分、蒙陰縣全部、平邑縣大部分、費(fèi)縣大部分及臨沂市面積的30%。在這一流域中，河道長(zhǎng)約223 km，集水面積達(dá)到10 315 km2，大約占沂河流域總面積的87.3%。就地形分布而言，沂河流域中山丘區(qū)面積約為70%，平原區(qū)面積約占30%。山丘區(qū)主要分布在流域上游，土質(zhì)較為疏松，風(fēng)化層淺而薄，植被狀況相對(duì)較差，不利于涵養(yǎng)水分，平原區(qū)主要分布在流域中下游，土質(zhì)多為壤土，植被狀況較好。由于流域峰時(shí)漫溢、枯時(shí)斷流的特性，沂河流域成為旱澇災(zāi)害多發(fā)區(qū)。目前，沂河流域經(jīng)濟(jì)水平相對(duì)落后，在東部沿海地帶屬于經(jīng)濟(jì)低谷區(qū)[15]。

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

2.1 數(shù)據(jù)來(lái)源

以研究區(qū)2005年和2015年LandsatTM遙感影像為基礎(chǔ)數(shù)據(jù)，以沂河流域地形圖、土地利用現(xiàn)狀圖及其他輔助圖件為參考，結(jié)合流域氣候、土壤、植被等自然地理資料，將遙感影像進(jìn)行人機(jī)交互目視解譯，在遙感影像解譯過(guò)程中，對(duì)于判讀不準(zhǔn)的圖斑進(jìn)行野外實(shí)地考察。最后得到沂河流域2005年和2015年矢量格式的全區(qū)土地利用數(shù)據(jù)。

參照全國(guó)土地利用分類(lèi)方法，將研究區(qū)土地利用類(lèi)型分為耕地、林地、草地、水域、建設(shè)用地和未利用土地六大類(lèi)。從研究區(qū)的統(tǒng)計(jì)年鑒中獲得其他所需要的自然和社會(huì)經(jīng)濟(jì)數(shù)據(jù)。

圖1 沂河流域地理位置

2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)價(jià)方法

1997年，Costanza等最先提出了生態(tài)服務(wù)價(jià)值評(píng)估的方法和原理，但其研究高估了濕地生態(tài)價(jià)值而低估了耕地價(jià)值，引起國(guó)內(nèi)外學(xué)者廣泛爭(zhēng)議[2]。2003年，謝高地等參考Costanza的部分研究成果，結(jié)合中國(guó)實(shí)際，對(duì)我國(guó)200多位專(zhuān)業(yè)學(xué)者進(jìn)行了問(wèn)卷調(diào)查，制定出中國(guó)陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量表[5]。2007年，謝高地等在有關(guān)問(wèn)卷調(diào)查基礎(chǔ)上，又建立了新的生態(tài)系統(tǒng)服務(wù)功能評(píng)估體系，并且這次調(diào)查結(jié)果更符合我國(guó)生態(tài)環(huán)境的實(shí)際。謝高地2007年生態(tài)系統(tǒng)單位面積價(jià)值當(dāng)量表只適合大尺度區(qū)域的研究，但此研究的沂河流域只是個(gè)中小尺度流域[16]。因此本研究根據(jù)謝高地等提出的生態(tài)系統(tǒng)服務(wù)價(jià)值的區(qū)域修正系數(shù)表，獲得山東省農(nóng)田生態(tài)系統(tǒng)的生物量因子為1.38，對(duì)謝高地等在2007年制定的生態(tài)服務(wù)價(jià)值表進(jìn)行修正，最終制定了適合沂河流域的生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量如表1所示[17]。

表1 沂河流域生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量

生態(tài)系統(tǒng)服務(wù)價(jià)值的計(jì)算方法[17]：

式中：Vij——第i種生態(tài)系統(tǒng)的第j種生態(tài)服務(wù)功能的經(jīng)濟(jì)價(jià)值；Si——流域第i種生態(tài)系統(tǒng)土地利用類(lèi)型的面積；Eij——第i種生態(tài)系統(tǒng)的第j種生態(tài)服務(wù)功能單位面積的經(jīng)濟(jì)價(jià)值；ESV——流域生態(tài)系統(tǒng)服務(wù)總價(jià)值。

3 結(jié)果分析

3.1 土地利用變化分析

通過(guò)兩期遙感影像解譯結(jié)果，統(tǒng)計(jì)出2005—2015年沂河流域土地利用現(xiàn)狀，見(jiàn)表2。從表2可以看出，耕地占據(jù)流域面積比例最大，占到60%以上，其次是林地，約占12%以上。草地、水域和建設(shè)用地所占比例較少。從2005—2015年，沂河流域耕地減少56 161.93 hm2，草地減少3 255.47 hm2，林地增加17 924.47 hm2，建設(shè)用地面積增加了38 086.62 hm2，表明流域內(nèi)主要為耕地、林地和建設(shè)用地之間的相互轉(zhuǎn)化。在ARCGIS10.2中將兩期遙感影像疊合得到土地利用轉(zhuǎn)移矩陣，根據(jù)土地利用轉(zhuǎn)移矩陣可以看出耕地主要向建設(shè)用地和林地轉(zhuǎn)化，草地主要轉(zhuǎn)化為林地。

為了更全面地分析沂河流域的土地利用狀況，本研究對(duì)土地利用程度進(jìn)一步加以分析。土地利用程度不僅體現(xiàn)土地本身的自然屬性，同時(shí)也反映土地利用過(guò)程中人為因素與自然環(huán)境因素的綜合效應(yīng)[18]。根據(jù)劉紀(jì)遠(yuǎn)等[19]提出的數(shù)量化綜合分析方法，將沂河流域土地利用程度分為四個(gè)等級(jí)，如表3[20]。

表22005 年與2015沂河流域土地利用狀況比較

表3 沂河流域土地利用程度分級(jí)指數(shù)

根據(jù)土地利用程度綜合指數(shù)模型和土地利用程度變化模型，得到沂河流域2005—2015年的土地利用程度綜合指數(shù)及土地利用程度變化量，見(jiàn)表4。從表4可以看出，2005—2015年，沂河流域土地利用綜合指數(shù)總體呈現(xiàn)增加趨勢(shì)，從283.41上升到285.37，△L的變化量為1.96，表明該流域土地利用程度不斷加強(qiáng)，土地利用正處于開(kāi)發(fā)時(shí)期。主要原因可能是區(qū)域城市化水平不斷提高，人民生活水平的提高和收入的增多改變了人們的消費(fèi)結(jié)構(gòu)和消費(fèi)觀念，間接推動(dòng)了土地利用結(jié)構(gòu)的改變。

表4 沂河流域各時(shí)期土地利用程度綜合指數(shù)

3.2 生態(tài)服務(wù)價(jià)值變化分析

文獻(xiàn)資料和調(diào)查結(jié)果顯示，沂河流域主要的糧食作物包括小麥、谷子、玉米、高粱、大豆、綠豆、稻谷和薯類(lèi)。根據(jù)2005—2015年沂河流域主要糧食作物的平均播種面積、面積比重、平均單產(chǎn)量以及單價(jià)（其中各糧食作物的單價(jià)參考2010年臨沂市糧食價(jià)格），利用單位面積耕地食物生產(chǎn)服務(wù)價(jià)值計(jì)算公式計(jì)算出沂河流域單位面積耕地的自然糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值為1 876.84元/（hm2·a)。結(jié)合其它生態(tài)系統(tǒng)服務(wù)功能單位面積價(jià)值計(jì)算公式，最終可獲得沂河流域生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值系數(shù)表，進(jìn)而得到2005年和2015年沂河流域生態(tài)服務(wù)價(jià)值表。

從表6和表7可以看出2005—2015年間，沂河流域生態(tài)服務(wù)總價(jià)值顯著上升，從2 338 964萬(wàn)元上升到2 427 039萬(wàn)元，增加了88 075萬(wàn)元，這期間，草地生態(tài)服務(wù)價(jià)值基本不變，耕地生態(tài)服務(wù)價(jià)值顯著下降，但林地、水域和建設(shè)用地生態(tài)服務(wù)價(jià)值顯著增加，增加的生態(tài)服務(wù)價(jià)值要遠(yuǎn)遠(yuǎn)多于減少的，所以整體是增加的。耕地、林地的面積變化在沂河流域生態(tài)服務(wù)價(jià)值的研究中起到了決定性的作用。

表5 沂河流域生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值系數(shù)元/（hm2·a）

表62005 年沂河流域生態(tài)系統(tǒng)服務(wù)價(jià)值萬(wàn)元

表72015 年沂河流域生態(tài)系統(tǒng)服務(wù)價(jià)值萬(wàn)元

4 結(jié)論與討論

4.1 結(jié)論

4.1.1 土地利用變化特征

耕地在該流域所占面積比例最大，在整個(gè)研究時(shí)期都占據(jù)了流域面積的60%以上，其次是林地，約占總比例的12%，然后是草地，占總比例的10%，水域和未利用地所占面積比例較少，水域約占3.5%，未利用地占0.13%。近10年來(lái)沂河流域耕地顯著減少；建設(shè)用地在整個(gè)研究時(shí)期面積不斷增加，隨著城市化水平的提高，建設(shè)用地急劇擴(kuò)張；林地面積也不斷增加，這主要是國(guó)家退耕還林政策的實(shí)施和流域果園林的增多造成的。在研究期間，草地、水域和未利用地變化不大。從土地利用動(dòng)態(tài)度和土地利用程度可以看出，該流域土地利用變化屬于緩慢變化型，正處于蓬勃發(fā)展階段。

4.1.2 流域生態(tài)服務(wù)價(jià)值變化特征

2005—2015年，流域生態(tài)服務(wù)價(jià)值顯著升高，從2 338 964萬(wàn)元上升到2 427 039萬(wàn)元，升高了88 075萬(wàn)元，在這段時(shí)期耕地的面積顯著減少，林地和建設(shè)用地的面積顯著增加，流域土地利用處于快速發(fā)展階段?？赡苁怯捎诔鞘谢降牟粩嗵岣撸嗣裆钏降奶岣吆褪杖氲脑龆喔淖兞巳藗兊南M(fèi)結(jié)構(gòu)和消費(fèi)觀念，間接推動(dòng)了土地利用結(jié)構(gòu)的改變。合理的規(guī)劃布局和生態(tài)保護(hù)帶來(lái)了生態(tài)服務(wù)價(jià)值的顯著提升。

通過(guò)本研究可以看出，沂河流域的耕地面積正在加速減少，流域的土地資源不斷被侵蝕，長(zhǎng)此以往，流域的耕地資源將難以保障當(dāng)?shù)鼐用竦男枰瑹o(wú)法實(shí)現(xiàn)區(qū)域的可持續(xù)發(fā)展。所以為了區(qū)域的可持續(xù)發(fā)展，沂河流域應(yīng)合理的分配各種土地資源，保護(hù)耕地、水源、林地這些生態(tài)服務(wù)價(jià)值較高的土地利用類(lèi)型，發(fā)揮資源潛力。政府應(yīng)制定土地保護(hù)政策，限制建設(shè)用地的過(guò)度擴(kuò)張，使該區(qū)域處在一個(gè)健康的發(fā)展模式下。

4.2 討論

本研究通過(guò)謝高地等在2007年制定的生態(tài)服務(wù)價(jià)值表對(duì)沂河流域進(jìn)行研究，雖然針對(duì)山東地區(qū)進(jìn)行了修正，但由于生態(tài)環(huán)境的復(fù)雜多變，區(qū)域之間氣候條件各有不同，因此在運(yùn)用生態(tài)服務(wù)價(jià)值當(dāng)量表計(jì)算生態(tài)服務(wù)價(jià)值時(shí)會(huì)存在一定的誤差。本研究數(shù)據(jù)量過(guò)大，有些數(shù)據(jù)在獲取時(shí)存在困難，雖然參考了大量資料，但數(shù)據(jù)來(lái)源的精確性依然存在不足。遙感影像在解譯時(shí)有時(shí)候草地、林地的界限難以明確區(qū)分，會(huì)給土地利用數(shù)據(jù)帶來(lái)一定誤差。

[1]蔣小榮，李丁，李智勇.基于土地利用的石羊河流域生態(tài)服務(wù)價(jià)值[J].中國(guó)人口·資源與環(huán)境，2010，20（6）：68-73.

[2]COSTANZA R D，ARGE R，GROOT R D，et al.The value of the world's ecosystem services and natural capital[J].Nature，1997，387：253-260．

[3]歐陽(yáng)志云，王如松，趙景柱.生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟(jì)價(jià)值評(píng)價(jià)[J].應(yīng)用生態(tài)學(xué)報(bào)，1999，10（5）：636-640.

[4]王廣成，李中才.基于時(shí)空尺度及利益關(guān)系的生態(tài)服務(wù)功能[J].生態(tài)學(xué)報(bào)，2007，27（11）：4758-4765.

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

[6]Ayres R U.The price value Paradox[J].Ecological Economics，1998，25（1）：17-19.

[7]SERAFY S.Pricing the invaluable：the value of the world's ecosystem service and natural capital[J].Ecological Economics，1998，25（1）：25-27.

[8]吳鋼，肖寒，趙景柱，等.長(zhǎng)白山森林生態(tài)系統(tǒng)服務(wù)功能[J].中國(guó)科學(xué)（C輯：生命科學(xué)），2001，31（5）：471-480.

[9]TURNER B LⅡ，SKOLE D，F(xiàn)ISCHER G，et al，ed.Landuse landcover change：science/change plan[C]//IGBP Report No.35 and HDP Report No.7.Stockholm and Geneva，1995：14-22.

[10]陳姝，居為民，李顯風(fēng).常熟市土地利用變化對(duì)生態(tài)服務(wù)價(jià)值的影響[J].水土保持研究，2009，16（5）：93-97.

[11]鞏燦娟.沂河流域中上游地區(qū)土地利用/覆被變化水文效應(yīng)的分析與模擬[D].山東：山東師范大學(xué)，2014.

[12]陸建宇，陸寶宏，朱從飛，等.沂河流域天然徑流變化規(guī)律分析[J].中國(guó)農(nóng)村水利水電，2014（7）：67-71.

[13]宋曉猛，孔凡哲.沂河流域徑流演變特征及動(dòng)因分析[J].安徽農(nóng)業(yè)科學(xué)，2010（16）：8563-8565，8572.

[14]徐智廷，趙光武，孫廷璽，等.沂河流域地理特性及人文特點(diǎn)分析[J].中國(guó)水利，2006，38（18）：50-51.

[15]薛麗芳，譚海樵.沂河流域水文特征變化及其驅(qū)動(dòng)因素[J].地理科學(xué)進(jìn)展，2011，30（11）：1354-1360.

[16]謝高地，甄霖，魯春霞，等.一個(gè)基于專(zhuān)家知識(shí)的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法[J].自然資源學(xué)報(bào)，2008（5）：911-919.

[17]錢(qián)翠.基于土地利用變化的祊河流域生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空動(dòng)態(tài)研究[D].山東：山東師范大學(xué)，2015.

[18]王思遠(yuǎn)，劉紀(jì)遠(yuǎn)，張?jiān)鱿?，?中國(guó)土地利用時(shí)空特征分析[J].地理學(xué)報(bào)，2001，56（6）：631-639.

[19]劉紀(jì)遠(yuǎn).中國(guó)資源環(huán)境遙感宏觀調(diào)查與動(dòng)態(tài)研究[M].北京：中國(guó)科學(xué)技術(shù)出版社，1996：262-274.

[20]莊大方，劉紀(jì)遠(yuǎn).中國(guó)土地利用程度的區(qū)域分異模型研究[J].自然資源學(xué)報(bào)，1997，12（2）：105-111.

（編輯：周利海）

A Study based on Ecological Service Value of Land Use Change of Yi River Basin

Su Ning1,Fan Guijuan2
（Shandong Normal University,Jinan Shandong 250358,China）

Taking Yi River Basin as an example,interpreted the 2005 and 2015 remote sensing image of Yi River Basin in ARCGIS10.2,finally got two issues of land use data.Analyzed the 2005—2015 Yi River Basin land use change,and referring to the calculation method of ecological service value of Costanza,analyzed the change of ecosystem service value of Yi River Basin 2005—2015.The results showed that the cultivated land in the study area was significantly reduced,but the total value of ecosystem services has increased from 23 389.64 million yuan in 2005 to 24 270.39 million yuan in 2015,mainly due to the ecological service value brought by the increase of ecological service forest area.

Yi River Basin,land use change,ecological service value

X196

A

1008-813X（2017）04-0032-05

10.13358 /j.issn.1008-813x.2017.04.09

2017-06-14

蘇寧（1991-），男，山東濟(jì)寧人，山東師范大學(xué)自然地理學(xué)專(zhuān)業(yè)（土地利用及覆被變化）碩士研究生在讀，主要從事土地利用及覆被變化等方面的研究。