周浩浩 徐清昊

(長(zhǎng)安大學(xué)土地工程學(xué)院，陜西 西安 710061)

土地利用是反映人類生活和生態(tài)系統(tǒng)之間耦合關(guān)系最根本的表現(xiàn)[1]。土地利用的變化會(huì)改變生態(tài)系統(tǒng)服務(wù)體系,進(jìn)而改變生態(tài)系統(tǒng)服務(wù)價(jià)值。生態(tài)系統(tǒng)服務(wù)價(jià)值(Ecosystem service value,ESV)首次由Costanza[2]提出,謝高地等[3]在此基礎(chǔ)上改進(jìn)了ESV系數(shù),估算了我國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值,建立了一套可廣泛應(yīng)用于評(píng)估我國(guó)生態(tài)系統(tǒng)服務(wù)的ESV系數(shù)；當(dāng)量因子法由于具有操作簡(jiǎn)單、綜合評(píng)價(jià)、適應(yīng)性強(qiáng)的優(yōu)點(diǎn)而被普遍采用[4]。

1 研究區(qū)概況

陜北地區(qū)位于黃土高原的核心地區(qū),在行政區(qū)劃上包含陜西省的延安和榆林兩市,北與內(nèi)蒙古的毛烏素沙地接壤,東部緊鄰山西省,西部接寧夏回族自治區(qū)和甘肅省。總面積92521.4km2,基本地貌類型是黃土塬、梁、峁、溝、塬,氣候干旱少雨,近年來(lái)經(jīng)濟(jì)發(fā)展迅速,多樣的地貌、流水侵蝕和各種程度的人類活動(dòng)作用下,生態(tài)系統(tǒng)比較脆弱,易受破壞。

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

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

本文的土地利用數(shù)據(jù)為自中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心獲取的2000年、2005年、2010年、2015年、2020年5期的土地利用遙感監(jiān)測(cè)數(shù)據(jù),社會(huì)經(jīng)濟(jì)數(shù)據(jù)來(lái)源為國(guó)家統(tǒng)計(jì)局官網(wǎng)的《陜西省統(tǒng)計(jì)年鑒》和《全國(guó)農(nóng)產(chǎn)品成本收益資料》等。

2.2 研究方法

2.2.1 土地利用轉(zhuǎn)移矩陣

也稱馬爾科夫矩陣,可全面反映區(qū)域內(nèi)各用地類型的轉(zhuǎn)移方向與轉(zhuǎn)移數(shù)量,被普遍應(yīng)用于土地利用變化研究,能很好地展示土地利用格局的時(shí)空演變過(guò)程[5],計(jì)算公式:

式中,Sij代表n×n矩陣；S代表面積；n代表土地類型數(shù)；i和j分別代表研究期初始和末尾的土地類型。

2.2.2 單一土地利用動(dòng)態(tài)度

單一土地利用動(dòng)態(tài)度可以直觀地反映某種土地類型的變化趨勢(shì)以及變化速度,計(jì)算公式:

式中,K表示某一類土地利用類型動(dòng)態(tài)度；Ua和Ub表示研究期初始和末尾某地類的面積；T表示研究時(shí)長(zhǎng),當(dāng)T的時(shí)段設(shè)定為年時(shí)；K為研究時(shí)段內(nèi)某一土地利用類型的年變化率。

2.2.3 綜合土地利用動(dòng)態(tài)度

綜合土地利用動(dòng)態(tài)度反映某種土地類型轉(zhuǎn)變?yōu)槠渌恋仡愋偷淖兓俾?計(jì)算公式:

式中,LC為研究期內(nèi)某區(qū)域的綜合土地利用動(dòng)態(tài)度；LUi為研究初期第i類土地利用類型面積；ΔLUi-j為第i類土地利用類型面積轉(zhuǎn)為非i類土地利用類型面積的絕對(duì)值；T為研究時(shí)長(zhǎng)。

2.2.4 生態(tài)服務(wù)價(jià)值計(jì)算

謝高地等學(xué)者的當(dāng)量因子法首次測(cè)算了全國(guó)的生態(tài)服務(wù)價(jià)值,為后續(xù)的研究提供了思路,并且具有較強(qiáng)的可操作性,被廣泛采用。劉浩等[6]指出,傳統(tǒng)當(dāng)量因子法計(jì)算出的結(jié)果存在誤差。本文結(jié)合研究區(qū)的實(shí)際情況,選取小麥和玉米2種主要糧食產(chǎn)物為代表[7],根據(jù)農(nóng)產(chǎn)品的單位面積的產(chǎn)值進(jìn)行修正,根據(jù)修正結(jié)果,1個(gè)標(biāo)準(zhǔn)生態(tài)系統(tǒng)服務(wù)經(jīng)濟(jì)價(jià)值等效系數(shù)為1039.25元·hm-2,進(jìn)一步計(jì)算陜北地區(qū)的生態(tài)服務(wù)價(jià)值總量。計(jì)算公式:

式中,ESV是研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值；VCi是第i種地類的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值；Ai為第i種地類的面積；i為土地利用類型；ESVf是第f項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值；VCfi是第i種地類的單位面積的第f項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值；f為生態(tài)系統(tǒng)服務(wù)價(jià)值的項(xiàng)數(shù)。

表1 陜北生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量

3 結(jié)果與分析

3.1 土地利用時(shí)空變化

近20a內(nèi),陜北地區(qū)的土地利用類型在空間的分布上,主要以林地、耕地和草地為主,三者之和達(dá)研究區(qū)面積的93%左右,林地主要分布于研究區(qū)的南部地區(qū),且集中分布于黃龍山和洛河河網(wǎng)地區(qū),耕地和草地在研究區(qū)廣泛分散分布,且在研究區(qū)的東北部地區(qū)有集中趨勢(shì),未利用地主要為與毛烏素沙地連接的沙漠化土地,集中分布于榆林北部地區(qū),建設(shè)用地呈中心發(fā)散狀,集中分布于延安和榆林市的中心城鎮(zhèn)地區(qū)。分布位置見(jiàn)圖1。

圖1 陜北土地利用類型

陜北地區(qū)土地利用在時(shí)間上發(fā)生的變化,從總體上來(lái)看,近20a間的變化主要在林草地與耕地之間發(fā)生,表現(xiàn)為耕地的集聚減少與林草地的快速增加同步進(jìn)行。從具體地類上來(lái)看,耕地面積從2000年的28453.44km2減少至2020年的25189.90km2,共計(jì)減少了3263.54km2,占研究區(qū)總面積的4%；與之相對(duì)應(yīng)的是林地,由最初的11052.39km2增加至12289.15 km2,凈增加1236.76km2,增幅約10%；草地由2000年的34753.43km2增 加 至36083.63km2,凈 增 加1330.20km2,增加面積占總面積的2%；建設(shè)用地呈現(xiàn)急劇擴(kuò)張態(tài)勢(shì),由初始的297.90km2增加至1375.49km2,增幅達(dá)450%；以沙地為代表的未利用地,減少了約378km2,水域面積基本上保持不變。具體變化見(jiàn)表2。

表2 陜北土地利用類型面積

3.2 土地利用結(jié)構(gòu)變化

20a來(lái),除去水域面積波動(dòng)不大外,其他5種地類均有明顯的地類轉(zhuǎn)移變化。2000—2020年,耕地呈大幅減少態(tài)勢(shì),主要流出為林地、草地和建設(shè)用地,對(duì)其進(jìn)行補(bǔ)充的地類主要為草地、未利用地和林地；林地轉(zhuǎn)出面積較小,主要為耕地和草地對(duì)其流入；草地轉(zhuǎn)出方向主要是耕地、林地和建設(shè)用地,轉(zhuǎn)入的主要是耕地；水域面積有小幅度波動(dòng),基本保持不變；建設(shè)用地呈現(xiàn)快速增加,各個(gè)地類均有轉(zhuǎn)入；未利用地流出方向主要為耕地、草地和建設(shè)用地,轉(zhuǎn)入的主要地類為草地。究其原因,各個(gè)地類之間的轉(zhuǎn)變主要受政策影響較大,耕地的大幅度轉(zhuǎn)出,與黃土高原的退耕還林還草政策緊密相關(guān),建設(shè)用地轉(zhuǎn)入也有充分利用沙漠化荒漠化的劣地。具體的結(jié)構(gòu)變化見(jiàn)表3。

3.3 土地利用動(dòng)態(tài)度

從單一土地利用動(dòng)態(tài)度方面來(lái)講,2000—2020年,耕地、水域和未利用地的單一動(dòng)態(tài)度為負(fù)值,建設(shè)用地、林地和草地的動(dòng)態(tài)度為正值,且建設(shè)用地的動(dòng)態(tài)度最大為18.09%,說(shuō)明耕地、未利用地和水域呈現(xiàn)減少趨勢(shì),建設(shè)用地、草地和林地呈增加趨勢(shì),且建設(shè)用地?cái)U(kuò)張最為迅速。分時(shí)段來(lái)看,2000—2010年,此前10a耕地單一土地利用動(dòng)態(tài)度為較小的負(fù)值,表明前10a耕地急劇減少；與之相對(duì)應(yīng)的是2000—2005年和2005—2010年,林地和草地的動(dòng)態(tài)度為各自的峰值,說(shuō)明在這2個(gè)時(shí)段分別為退耕還林和退耕還草最為迅速的時(shí)段；水域在整體上呈現(xiàn)波動(dòng)姿態(tài),在2005—2010年急劇減少,在2010—2015年快速增加；建設(shè)用地在2000—2005年為較小動(dòng)態(tài)度,從2005年開(kāi)始均以較大動(dòng)態(tài)變化,說(shuō)明陜北地區(qū)建設(shè)用地急速擴(kuò)張是從2005年左右開(kāi)始；未利用地一直以負(fù)值變化,說(shuō)明陜北地區(qū)治沙行動(dòng)從未停止,且在2005—2010年5a間成效最為顯著。從綜合土地利用總態(tài)度方面分析,2000—2010年10a間,綜合土地利用動(dòng)態(tài)度較大,且之后緩慢減少,顯而易見(jiàn),陜北地區(qū)在這段時(shí)間的地類轉(zhuǎn)移變化最為頻繁。具體動(dòng)態(tài)度數(shù)據(jù)見(jiàn)表4。

表4 陜北土地利用動(dòng)態(tài)度

3.4 生態(tài)系統(tǒng)服務(wù)價(jià)值分析

陜北地區(qū)生態(tài)服務(wù)價(jià)值在這近20a間呈先增加后趨于穩(wěn)定的趨勢(shì),從2000年的888.53億元增加至2020年的914.86億元,在整體上凈增加了26.33億元,說(shuō)明陜北地區(qū)整體的生態(tài)得到了改善。

從不同地類的生態(tài)系統(tǒng)服務(wù)價(jià)值來(lái)言,林草地為生態(tài)服務(wù)價(jià)值貢獻(xiàn)的主體,兩者之和占總體服務(wù)價(jià)值的65%以上,且林草的生態(tài)服務(wù)價(jià)值的變化趨勢(shì)與陜北地區(qū)整體的生態(tài)服務(wù)價(jià)值變化趨勢(shì)保持一致。耕地的服務(wù)價(jià)值呈現(xiàn)連續(xù)下降趨勢(shì),未利用地和水域的生態(tài)服務(wù)價(jià)值在較小幅度保持減少,變化不大,未利用地的生態(tài)服務(wù)占比最少。具體各地類的生態(tài)服務(wù)價(jià)值如表5所示。

表5 陜北各地類生態(tài)服務(wù)價(jià)值

就不同的生態(tài)服務(wù)功能方面來(lái)看,一級(jí)服務(wù)功能中調(diào)節(jié)服務(wù)占比最大,達(dá)總量的50%以上,接下來(lái)為支持、供給和文化服務(wù)。除供給服務(wù)價(jià)值呈現(xiàn)降低趨勢(shì)外,其他3種均呈現(xiàn)緩慢增加趨勢(shì)；在二級(jí)服務(wù)功能中氣候調(diào)節(jié)和水文調(diào)節(jié)服務(wù)功能對(duì)生態(tài)系統(tǒng)總體的服務(wù)功能貢獻(xiàn)較大,分別為217.23億元和204.25億元,并且在變化量上,氣候調(diào)節(jié)生態(tài)服務(wù)價(jià)值的變化量最大,為11.58億元,且變化率最高,說(shuō)明在整體的生態(tài)服務(wù)價(jià)值變化中,二級(jí)服務(wù)價(jià)值氣候調(diào)節(jié)功能的改善做出巨大的貢獻(xiàn)。在所有的二級(jí)服務(wù)功能變化中,食物生產(chǎn)為唯一負(fù)值,這與耕地的減少密不可分,其他生態(tài)服務(wù)價(jià)值均有增加,說(shuō)明退耕在陜北地區(qū)對(duì)絕大部分的生態(tài)服務(wù)功能都有積極影響。具體生態(tài)服務(wù)價(jià)值變化見(jiàn)表6。

表6 陜北生態(tài)服務(wù)功能價(jià)值變化

4 結(jié)論

土地利用變化是生態(tài)系統(tǒng)服務(wù)功能轉(zhuǎn)變的重要驅(qū)動(dòng)因素[8],會(huì)改變生態(tài)系統(tǒng)服務(wù)的價(jià)值。根據(jù)土地利用進(jìn)行生態(tài)服務(wù)價(jià)值的研究切實(shí)可行的。通過(guò)本文的研究,得出以下結(jié)論。

近20a內(nèi)土地利用變化主要在耕地、林草地和建設(shè)用地之間發(fā)生較為劇烈的變化,具體表現(xiàn)為耕地的減少,林草的增加和建設(shè)用地以城市中心為基點(diǎn)的呈發(fā)散狀向外部急劇擴(kuò)張,并且在北部的城市擴(kuò)張同時(shí)利用了沙地,有較高程度的節(jié)約集約利用土地。在土地利用結(jié)構(gòu)轉(zhuǎn)變方面,耕地減少,并且主要轉(zhuǎn)變?yōu)椴莸亍⒘值睾臀蠢玫?較少量流向建設(shè)用地,林地增加,絕大部分是耕地、草地轉(zhuǎn)變?yōu)榱值?草地也增加,來(lái)源主要為耕地和未利用地,建設(shè)用地劇烈增加,來(lái)源主要是草地、耕地和未利用地中的沙地。在變化速率方面,只有建設(shè)用地最為突出,且從2005年后增幅巨大。深入探究原因可知,陜北地區(qū)的退耕還林還草政策,是耕地以及未利用地變化的主要原因,建設(shè)用地的變化與陜北地區(qū)的經(jīng)濟(jì)發(fā)展密不可分,并且未利用地的范圍線有明顯的后退,且在原來(lái)其所屬范圍內(nèi),有建設(shè)用地的明顯擴(kuò)張,呈現(xiàn)人進(jìn)沙退的風(fēng)貌。綜合來(lái)看陜北地區(qū)的風(fēng)沙治理,生態(tài)改善政策成效顯著,且經(jīng)濟(jì)活動(dòng)比較過(guò)熱,有較高的節(jié)約集約利用土地的意識(shí)。

近20a來(lái),根據(jù)量表測(cè)算的結(jié)果顯示,陜北地區(qū)的生態(tài)服務(wù)價(jià)值有較大的增加,由最初的888.53億元增加至914.86億元,凈增約26億元生態(tài)服務(wù)價(jià)值。從地類來(lái)看,林地和草地這2種地類對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的貢獻(xiàn)率最高,變化趨勢(shì)控制了整體服務(wù)價(jià)值的走勢(shì),且在前10a增幅明顯,后10a內(nèi)趨于穩(wěn)定。從生態(tài)服務(wù)功能方面來(lái)看,調(diào)節(jié)服務(wù)功能占比最大,并且逐年增加,在二級(jí)服務(wù)功能中氣候調(diào)節(jié)和水文調(diào)節(jié)占比較大,食物生產(chǎn)服務(wù)功能是所有二級(jí)服務(wù)功能中唯一減小的,其他功能均有所增加,追根溯源,與耕地的大幅減少和林草的增加密不可分?？偠灾?陜北地區(qū)在整體上生態(tài)環(huán)境改善,生態(tài)系統(tǒng)服務(wù)價(jià)值有所增加,除食物生產(chǎn)功能有所下降外,其他生態(tài)服務(wù)功能均得到改善提升。