榆中縣土地利用轉(zhuǎn)型與生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空演變研究

屈雯, 陳英, 謝保鵬, 裴婷婷, 任璽錦, 賀瑩

榆中縣土地利用轉(zhuǎn)型與生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空演變研究

屈雯, 陳英, 謝保鵬*, 裴婷婷, 任璽錦, 賀瑩

甘肅農(nóng)業(yè)大學(xué)管理學(xué)院, 蘭州 730070

生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)是區(qū)域生態(tài)安全的重要表征, 識(shí)別土地利用轉(zhuǎn)型引起的生態(tài)系統(tǒng)服務(wù)價(jià)值變化, 對(duì)區(qū)域生態(tài)環(huán)境保護(hù)具有重大意義。基于2008年和2018年榆中縣土地利用數(shù)據(jù), 借助ArcGIS 軟件, 運(yùn)用價(jià)值當(dāng)量因子法、生態(tài)系統(tǒng)服務(wù)價(jià)值變化指數(shù)、空間統(tǒng)計(jì)分析等方法, 分析了榆中縣土地利用轉(zhuǎn)型對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值變化的影響。結(jié)果表明: (1)2008—2018年, 榆中縣耕地、草地面積減少, 建設(shè)用地面積增加, 耕地—建設(shè)用地、草地—耕地、草地—建設(shè)用地是榆中縣土地轉(zhuǎn)型的主要模式。(2)研究期間, 榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值總體呈下降趨勢(shì), 其中草地ESV占比最大, 占比達(dá)54%, 單項(xiàng)ESV以水文調(diào)節(jié)和氣候調(diào)節(jié)為主, 二者累計(jì)占比接近ESV總量的50%。(3)榆中縣ESV空間分布具有明顯的空間正相關(guān)性, ESV的熱點(diǎn)區(qū)位于榆中縣南部, 冷點(diǎn)區(qū)位于榆中縣中部, ESV空間分布呈現(xiàn)南部高, 中部低, 北部較高的特征, 且ESV低值區(qū)在榆中縣中部進(jìn)一步擴(kuò)散。研究結(jié)果可以為區(qū)域生態(tài)環(huán)境保護(hù)和土地管理提供科學(xué)參考。

土地轉(zhuǎn)型; 生態(tài)系統(tǒng)服務(wù)價(jià)值; 空間自相關(guān); 熱點(diǎn)分析; 榆中縣

0 前言

土地利用轉(zhuǎn)型是指區(qū)域土地利用形態(tài)在時(shí)序上發(fā)生轉(zhuǎn)變的過(guò)程, 其通常會(huì)受社會(huì)經(jīng)濟(jì)的驅(qū)動(dòng)而發(fā)生改變[1]。作為全球變化研究的重要內(nèi)容, 土地利用轉(zhuǎn)型不僅影響景觀格局中土壤、水體、大氣、生物等環(huán)境要素的生態(tài)過(guò)程使地表自然景觀發(fā)生改變[2-3], 同時(shí)也造成水土流失、土壤污染、生態(tài)退化等一系列問(wèn)題[4], 進(jìn)而導(dǎo)致生態(tài)系統(tǒng)服務(wù)功能發(fā)生顯著的變化。生態(tài)系統(tǒng)服務(wù)是指直接或間接獲得的人類生活必需品以及人類生活質(zhì)量保證的生命支持產(chǎn)品及服務(wù)[5]。對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行定量評(píng)估是衡量生態(tài)系統(tǒng)服務(wù)功能好壞的關(guān)鍵[6-7], 是研究人類活動(dòng)對(duì)生態(tài)系統(tǒng)影響程度的重要手段[4]。近年來(lái), 隨著我國(guó)經(jīng)濟(jì)的快速發(fā)展和城市化進(jìn)程的不斷推進(jìn), 土地利用轉(zhuǎn)型過(guò)程發(fā)生了顯著的變化, 使得土地利用結(jié)構(gòu)和功能快速轉(zhuǎn)變[8-10], 進(jìn)而影響各類生態(tài)系統(tǒng)的類型、面積以及空間分布[11], 是生態(tài)系統(tǒng)服務(wù)價(jià)值變化的重要驅(qū)動(dòng)力。因此, 加強(qiáng)區(qū)域土地利用轉(zhuǎn)型和生態(tài)系統(tǒng)服務(wù)價(jià)值變化的研究, 對(duì)促進(jìn)區(qū)域的協(xié)調(diào)發(fā)展有重大的意義。

Cosatanza[5]等學(xué)者, 最早建立了定量評(píng)估生態(tài)系統(tǒng)服務(wù)價(jià)值的方法。我國(guó)學(xué)者歐陽(yáng)志云[12]等最先引進(jìn)生態(tài)系統(tǒng)服務(wù)的概念。謝高地[13]等基于對(duì)Cosatanza評(píng)價(jià)模型的改進(jìn)制定了中國(guó)陸地生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值表, 并在2015年對(duì)該表進(jìn)行了進(jìn)一步的修正[14]。劉永強(qiáng)[15]等基于土地利用轉(zhuǎn)型對(duì)湖南省的生態(tài)系統(tǒng)服務(wù)價(jià)值效應(yīng)進(jìn)行了分析, 并提出了土地利用優(yōu)化的配置建議。王權(quán)[8]等分析了土地利用轉(zhuǎn)型對(duì)喀斯特槽谷區(qū)ESV的影響。戴文遠(yuǎn)[16]等從“三生空間”的角度出發(fā), 分析了土地利用轉(zhuǎn)型對(duì)福州新區(qū)ESV的影響。從已有的研究來(lái)看, 大多數(shù)研究分析了土地利用變化引起的生態(tài)系統(tǒng)服務(wù)價(jià)值轉(zhuǎn)變[17-21], 以土地利用轉(zhuǎn)型為視角分析生態(tài)系統(tǒng)服務(wù)價(jià)值的研究較少, 研究尺度主要為省、市、城市群、流域, 對(duì)微觀尺度網(wǎng)格單元的關(guān)注較少, 且研究側(cè)重于分析ESV數(shù)量上的變化, 對(duì)ESV空間格局上的研究較少。因此, 本研究嘗試以網(wǎng)格單元為研究視角, 從數(shù)量和空間兩方面對(duì)ESV進(jìn)行分析。

榆中縣屬典型的西北干旱區(qū), 位于三大自然帶的過(guò)渡區(qū)域[22], 其土地利用具有明顯的地域性特征, 縣域內(nèi)具有興隆山國(guó)家級(jí)自然保護(hù)區(qū), 是蘭州市乃至黃河上游區(qū)域重要的生態(tài)屏障區(qū)。上位規(guī)劃將榆中縣作為蘭州市發(fā)展的承接區(qū), 縣域人為活動(dòng)不斷強(qiáng)化, 如何協(xié)調(diào)土地利用與生態(tài)保護(hù)并推進(jìn)建設(shè)發(fā)展是該區(qū)面臨的重大問(wèn)題。以榆中縣為研究區(qū), 探討快速城鎮(zhèn)化背景下的土地利用轉(zhuǎn)型與生態(tài)系統(tǒng)服務(wù)價(jià)值的變化, 以期為榆中縣的可持續(xù)發(fā)展和生態(tài)文明建設(shè)提供參考。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

榆中縣(103°49′15″—104°34′40″E, 35°34′20″—36°26′30″N)位于甘肅省中部地區(qū), 隸屬蘭州市管轄, 地處黃土高原西部, 同時(shí)也是西北干旱區(qū)、青藏高寒區(qū)和東部季風(fēng)區(qū)的交匯地帶[22], 地勢(shì)南高北低, 中部平坦, 呈馬鞍形分布, 南北長(zhǎng)97 km, 東西寬66 km, 現(xiàn)轄11鎮(zhèn)9鄉(xiāng)。榆中縣地形復(fù)雜, 土地利用方式多種多樣, 中部為川塬河谷平原區(qū), 地勢(shì)平坦, 土壤肥沃, 人民出行方便, 是榆中縣主要的農(nóng)業(yè)生產(chǎn)區(qū), 北部系干旱山區(qū), 氣候干旱, 降水稀少, 植被稀疏, 是榆中縣水土流失山最為嚴(yán)重的區(qū)域, 南部為二陰黑土層和高寒陰濕區(qū), 馬銜山和興隆山從西北向東南橫列整個(gè)縣域的南部, 森林資源豐富, 是榆中縣主要的生態(tài)保障區(qū)。

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

本研究包括的基本數(shù)據(jù)有: (1)榆中縣糧食平均單產(chǎn)數(shù)據(jù)來(lái)源于2008—2018年《蘭州市統(tǒng)計(jì)年鑒》和2008—2018年榆中縣國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào); (2)中國(guó)1個(gè)生態(tài)系統(tǒng)價(jià)值當(dāng)量的經(jīng)濟(jì)價(jià)值數(shù)據(jù)來(lái)源于前人研究; (3)榆中縣2008年和2018年土地利用數(shù)據(jù)來(lái)源于榆中縣土地利用變更數(shù)據(jù)庫(kù)。參照全國(guó)土地利用覆蓋分類體系和謝高地學(xué)者的研究, 并結(jié)合研究區(qū)的土地利用特征, 利用ArcGIS 10.2軟件平臺(tái)對(duì)變更數(shù)據(jù)中的用地類型進(jìn)行進(jìn)一步劃分, 為耕地、林地、草地、水域、濕地、建設(shè)用地和未利用地7大類。

1.3 研究方法

1.3.1 土地利用結(jié)構(gòu)變化分析方法

土地利用結(jié)構(gòu)變化選用土地利用動(dòng)態(tài)度和土地利用變化重要性指數(shù)進(jìn)行分析。其中土地利用動(dòng)態(tài)度可以直觀的反映某一種土地利用類型面積的增加或減少[23], 計(jì)算方法如下:

式中:為研究期內(nèi)第類土地利用類型的動(dòng)態(tài)度,為某一種土地利用類型的變化年份,為研究初期第種土地利用類型的面積,為研究末期第種土地利用類型的面積。

土地利用變化重要性指數(shù)用來(lái)反映研究區(qū)主要的土地利用變化類型[24], 計(jì)算方法如下:

式中:為第種土地利用類型的土地利用變化重要性指數(shù),為第類土地利用類型變化了的面積,為研究區(qū)各土地利用類型變化了的面積之和。

1.3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值測(cè)算方法

根據(jù)謝高地學(xué)者2015年修正的價(jià)值當(dāng)量表, 結(jié)合研究區(qū)實(shí)際情況, 選取榆中縣2000—2018年主要糧食作物的單位面積產(chǎn)量與同期全國(guó)單位面積糧食產(chǎn)量的比值作為修正系數(shù), 確定榆中縣的生態(tài)系統(tǒng)服務(wù)當(dāng)量?jī)r(jià)值的修正系數(shù)為0.66, 得到榆中縣1個(gè)標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量為2260.63(元·hm-2)。考慮到榆中縣水澆地面積較大, 結(jié)合趙鴻雁[25]對(duì)甘肅省水澆地生態(tài)系統(tǒng)服務(wù)價(jià)值的研究, 對(duì)榆中縣耕地的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行修正, 其他土地利用類型對(duì)應(yīng)的生態(tài)系統(tǒng)服務(wù)價(jià)值的確定參照謝高地2015年修訂的當(dāng)量因子表。在確定單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量的基礎(chǔ)上, 結(jié)合各土地利用類型的面積得到榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值, 計(jì)算公式為:

式中:為研究區(qū)的生態(tài)系統(tǒng)服務(wù)總價(jià)值,為第類生態(tài)系統(tǒng)的面積,VC為第類生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù),為研究區(qū)第項(xiàng)生態(tài)系統(tǒng)服務(wù)功能的價(jià)值,為第類第項(xiàng)生態(tài)服務(wù)價(jià)值功能系數(shù)。

生態(tài)系統(tǒng)服務(wù)變化指數(shù)(ESCI), 可以用來(lái)表征每項(xiàng)生態(tài)系統(tǒng)服務(wù)的相對(duì)增益或損失[26]。計(jì)算方法如下:

值為0時(shí), 說(shuō)明生態(tài)系統(tǒng)服務(wù)價(jià)值沒有變化[27], 即指生態(tài)系統(tǒng)服務(wù)價(jià)值既沒有增益也沒有減損; 當(dāng)?shù)闹敌∮?時(shí), 說(shuō)明生態(tài)系統(tǒng)服務(wù)價(jià)值有減損; 當(dāng)ESCI的值大于0時(shí), 說(shuō)明生態(tài)系統(tǒng)服務(wù)價(jià)值有增益。

式中:為第項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值的變化指數(shù),CURx為最終狀態(tài)下的生態(tài)服務(wù)價(jià)值,HISx最初狀態(tài)下的生態(tài)服務(wù)價(jià)值。

表1 榆中縣單位面積生態(tài)服務(wù)價(jià)值當(dāng)量

1.3.3 空間統(tǒng)計(jì)分析

空間自相關(guān)分析(Moran's)一般用來(lái)表示某一地理事物或現(xiàn)象在某一區(qū)域的不同位置上的空間相關(guān)程度[28], 可以衡量ESV的分布是否具有集聚性, 冷熱點(diǎn)分析Hot Spot Analysis(Getis-Ord Gi*)用以衡量ESV空間變化的聚集與分異特征, 探究ESV的空間變化是否具有高值集聚(熱點(diǎn))和低值集聚(冷點(diǎn))的現(xiàn)象[29], 通過(guò)熱點(diǎn)分析, 可以確定ESV高值區(qū)或低值區(qū)在空間上發(fā)生聚類的位置[30], 計(jì)算公式為:

式中:為研究區(qū)的空間網(wǎng)格單元數(shù)量,和分別為空間單元和空間單元的觀測(cè)值, (x-`)為第個(gè)空間單元上的觀測(cè)值與平均值的偏差,為空間單元和的權(quán)重矩陣。

2 結(jié)果與分析

2.1 土地利用轉(zhuǎn)型分析

2.1.1 土地利用轉(zhuǎn)型時(shí)空特征

榆中縣土地利用以草地和耕地為主(圖1), 耕地分布在縣域中部平原區(qū)和北部山區(qū), 草地分布在縣域中部。研究期內(nèi), 榆中縣土地利用類型變化比較顯著(表2, 圖2)。從土地利用變化數(shù)量上看, 草地、耕地、林地、濕地、未利用地面積有所減少, 其中, 耕地面積減少最多, 減少1798.00 hm2, 占比從38.69%減少到38.14%, 其次是草地, 面積減少621.93 hm2, 占比從41.04%減少到40.85%, 林地、濕地、未利用地面積有小幅度減少, 分別減少88.46 hm2、60.48 hm2、4.62 hm2; 建設(shè)用地面積增加2566.90 hm2, 占比從4.43%增加到5.21%。從土地利用動(dòng)態(tài)度來(lái)看, 耕地和濕地動(dòng)態(tài)度分別為-0.141, -0.265, 說(shuō)明耕地和濕地減少速度最快, 水域和建設(shè)用地的動(dòng)態(tài)度分別為0.121、1.756, 說(shuō)明建設(shè)用地和水域面積有所增加, 其中建設(shè)用地的動(dòng)態(tài)度最大, 增速最快。總體來(lái)看, 榆中縣耕地、草地面積減少最多, 建設(shè)用地面積增加最多。從空間變化上看, 建設(shè)用地變化最為明顯, 主要是由于蘭州市城市“東擴(kuò)”, 縣域人口持續(xù)增加, 建設(shè)用地需求不斷上漲, 從而導(dǎo)致境內(nèi)建設(shè)用地?cái)?shù)量快速增加。

2.1.2 土地利用轉(zhuǎn)型模式

從表3可以看出, 2008—2018年, 榆中縣約有0.8799%的土地利用類型發(fā)生了轉(zhuǎn)變。其中, 耕地轉(zhuǎn)出面積最多, 共轉(zhuǎn)出2120.68 hm2, 主要轉(zhuǎn)為建設(shè)用地, 占耕地轉(zhuǎn)出總面積的99.80%, 耕地流失主要集中在榆中縣中部的和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、連搭鎮(zhèn)、城關(guān)鎮(zhèn)、金崖鎮(zhèn)、夏官營(yíng)鎮(zhèn)、小康營(yíng)鄉(xiāng)、清水驛鄉(xiāng)、甘草店鎮(zhèn)和高崖鎮(zhèn)(圖2)。其次是草地, 主要轉(zhuǎn)為耕地和建設(shè)用地, 分別占草地轉(zhuǎn)出面積的50.06%、49.79%, 草地轉(zhuǎn)為耕地主要集中在青城鎮(zhèn)、和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、金崖鎮(zhèn)、高崖鎮(zhèn), 在甘草店鎮(zhèn)、小康營(yíng)鄉(xiāng)、馬坡鄉(xiāng)等其他鄉(xiāng)鎮(zhèn)有零星的分布, 草地轉(zhuǎn)為建設(shè)用地分布在榆中縣的中部地區(qū), 和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、連搭鎮(zhèn)、城關(guān)鎮(zhèn)、夏官營(yíng)鎮(zhèn)較為集中。林地、建設(shè)用地、濕地分別轉(zhuǎn)出87.55 hm2、5.48 hm2、57.96 hm2, 林地主要轉(zhuǎn)為建設(shè)用地, 占林地轉(zhuǎn)出面積的97.62%, 主要發(fā)生在和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、連搭鎮(zhèn)、城關(guān)鎮(zhèn)。建設(shè)用地轉(zhuǎn)為耕地, 主要發(fā)生在榆中縣中部和西北部(青城鎮(zhèn))的少數(shù)區(qū)域, 濕地主要轉(zhuǎn)為建設(shè)用地, 占濕地轉(zhuǎn)出面積的94.42%, 集中在榆中縣中部地區(qū), 未利用地轉(zhuǎn)為建設(shè)用地, 轉(zhuǎn)出面積為4.62 hm2, 分布比較零散, 水域主要轉(zhuǎn)為建設(shè)用地, 占水域轉(zhuǎn)出面積的94.22%。總體來(lái)看, 整個(gè)研究期內(nèi), 榆中縣土地利用類型均發(fā)生了變化, 耕地、草地、林地、建設(shè)用地和濕地的變化較為明顯。耕地—建設(shè)用地、草地—耕地、草地—建設(shè)用地是榆中縣土地轉(zhuǎn)型的主要類型, 重要性指數(shù)累計(jì)達(dá)94.43%(表4)。近年來(lái)榆中縣城鎮(zhèn)化速度不斷加快, 尤其是中部川地, 鄰近蘭州市主城區(qū), 人口不斷向該區(qū)集聚, 導(dǎo)致建設(shè)用地需求增加, 建設(shè)用地從城市近郊區(qū)不斷向外擴(kuò)散, 占用了城市周圍的耕地和草地, 從而使得城區(qū)外圍的耕地和草地減少。此外, 由于中部地區(qū)地形平坦, 土壤肥沃, 所以部分耕地占補(bǔ)平衡的項(xiàng)目在該區(qū)進(jìn)行了用地結(jié)構(gòu)調(diào)整, 導(dǎo)致部分草地轉(zhuǎn)為耕地。

表2 2008—2018年榆中縣土地利用及動(dòng)態(tài)度變化

2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)變化特征分析

2.2.1 不同土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值變化

2008—2018年榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值基本上保持平穩(wěn), 但有小幅度減少(表5), 生態(tài)系統(tǒng)服務(wù)價(jià)值由2008年的68.34億元減少到2018年的67.93億元, 減少0.41億元。從土地利用方式來(lái)看, 林地、草地構(gòu)成了榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值的主體, 2018年的貢獻(xiàn)率累計(jì)達(dá)到77.82%。研究期間草地和林地面積減少, 導(dǎo)致草地ESV減少0.17億元, 林地ESV減少0.04億元。耕地ESV減少0.15億元, 主要受到耕地面積變化的影響。濕地ESV減少0.07億元, 未利用地ESV減少0.00002億元, 減少幅度最小。縣域水庫(kù)修建使得水域面積增加, 從而導(dǎo)致水域ESV增加了0.02億元, 由此可見, 耕地和草地面積變化是造成榆中縣土地利用生態(tài)系統(tǒng)服務(wù)價(jià)值變化的主要原因。

圖1 榆中縣土地利用類型圖

Figure 1 Land use types in Yuzhong County

表3 2008—2018年榆中縣土地利用變化轉(zhuǎn)移矩陣 (hm2)

表4 2008—2018年榆中縣土地利用變化重要性指數(shù)

圖2 2008—2018年榆中縣土地利用變化空間分布圖

Figure 2 Spatial distribution of land use change in Yuzhong County from 2008 to 2018

表5 2008—2018年榆中縣各類土地生態(tài)系統(tǒng)服務(wù)價(jià)值變化

2.2.2 單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化

2008—2018年榆中縣單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化如表6所示。從單項(xiàng)生態(tài)系統(tǒng)一級(jí)服務(wù)功能來(lái)看, 對(duì)榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值貢獻(xiàn)率最大的是調(diào)節(jié)服務(wù), 2018年的貢獻(xiàn)率達(dá)68.51%, 其次是支持服務(wù), 2018年貢獻(xiàn)率為22.34%, 整個(gè)研究期內(nèi)榆中縣生態(tài)系統(tǒng)一級(jí)服務(wù)有所下降。從單項(xiàng)生態(tài)系統(tǒng)二級(jí)服務(wù)功能來(lái)看, 單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值以氣候調(diào)節(jié)和水文調(diào)節(jié)為主, 二者的累計(jì)貢獻(xiàn)率達(dá)到49%以上, 占比接近榆中縣總ESV的一半。其次是氣體調(diào)節(jié)和土壤保持, 2008年占ESV總量的比例分別為11.61%、11.05%, 2018年占ESV總量的比例分別為11.60%、11.05%。整個(gè)研究期內(nèi), 單項(xiàng)二級(jí)服務(wù)功能中, 水資源供給的ESV出現(xiàn)了增加, 主要是受到水域面積增加的影響, 其他二級(jí)服務(wù)功能均出現(xiàn)降低的趨勢(shì)。ESV減少最多的是水文調(diào)節(jié)和氣候調(diào)節(jié), 分別由2008年的17.60億元、16.01億元減少到2018年的17.48億元、15.93億元, 其次是氣體調(diào)節(jié)和土壤保持, ESV分別減少了0.06億元和0.04億元, 食物生產(chǎn), 原料生產(chǎn)和維持養(yǎng)分循環(huán)的變化最小。

2.2.3 生態(tài)系統(tǒng)服務(wù)變化指數(shù)空間分布

為查明榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布特征, 研究借助ArcGIS 10.2軟件中的漁網(wǎng)工具, 對(duì)研究區(qū)進(jìn)行了整體劃分, 在查閱相關(guān)文獻(xiàn)的基礎(chǔ)上[7, 31], 嘗試構(gòu)建了250 m×250 m、500 m×500 m、800 m×800 m、1000 m×1000 m的網(wǎng)格, 但考慮到研究區(qū)的自然地理特征以及數(shù)據(jù)處理效率等因素, 最終將確定為500 m×500 m的正方形網(wǎng)格, 并分別計(jì)算了每個(gè)網(wǎng)格單元中土地利用類型的面積及其所對(duì)應(yīng)的生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)和生態(tài)系統(tǒng)服務(wù)變化指數(shù)(ESCI), 應(yīng)用自然斷點(diǎn)法對(duì)ESV和ESCI分別制圖, 如圖3所示。

從ESV的空間分布來(lái)看(圖3), 榆中縣ESV的空間差異明顯, 從南到北總體上呈現(xiàn)出“高—低—高”的格局, 榆中縣南部區(qū)域?yàn)榕d隆山, 林地和草地是該區(qū)主要的土地利用類型, ESV值較高; 中部區(qū)域, 土地利用以耕地和建設(shè)用地為主, ESV值偏低;北部區(qū)域, 氣候干燥, 地形破碎, 土地利用以草地為主, ESV值較高。從2008年和2018年ESV的損益變化(ESCI)來(lái)看(圖4), 榆中縣ESCI的變化主要集中在中部地區(qū), 和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、連搭鎮(zhèn)、金崖鎮(zhèn)、城關(guān)鎮(zhèn)和夏官營(yíng)鎮(zhèn)的ESV變化明顯, 損失較大,且ESV低值區(qū)域分布明顯擴(kuò)大。此外和平鎮(zhèn)、城關(guān)鎮(zhèn)、連搭鎮(zhèn)、清水驛鄉(xiāng)、夏官營(yíng)鎮(zhèn)、貢井鎮(zhèn)、新營(yíng)鎮(zhèn)等區(qū)域ESV有零星的增加, 其他區(qū)域ESV的變化總體上不是特別明顯。

2.3 生態(tài)系統(tǒng)服務(wù)價(jià)值空間統(tǒng)計(jì)分析

由表7可知, 2008年和2018年榆中縣ESV的全局Moran's值均大于0,值均小于0.001, 說(shuō)明榆中縣ESV的分布格局整體上呈較強(qiáng)的空間正向自相關(guān)關(guān)系, ESV網(wǎng)格之間具有顯著的聚集分布特征。2008年和2018年的全局Moran's值分別為0.5986、0.6046, 整個(gè)研究期間有了一定的增加, 說(shuō)明榆中縣ESV的空間自相關(guān)性在逐漸增強(qiáng)。總的來(lái)看(表8), 2008年林地ESV的全局Moran's值為0.8316, 其空間集聚性最強(qiáng), 耕地、草地、未利用地的全局Moran's值依次次之, 值為0.7618、0.7385、0.6908, 說(shuō)明耕地、草地、未利用地ESV的空間分布也具有集聚性, 濕地、水域的全局Moran's I值分別為0.3704、0.4167, 濕地和水域的分布未呈現(xiàn)出明顯的空間集聚性, 不同土地利用類型的全局值排序?yàn)? 林地>耕地>草地>未利用地>濕地>水域, 說(shuō)明整個(gè)研究期間, 榆中縣ESV的空間分布總體上具有集聚性, 但集聚性是小范圍的、比較分散的。

表6 2008—2018年榆中縣土地生態(tài)系統(tǒng)單項(xiàng)服務(wù)價(jià)值變化

圖3 2008—2018年榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值空間分布圖

Figure 3 Spatial distribution of ecosystem service value in Yuzhong County from 2008 to 2018

圖4 2008—2018年榆中縣ESCI的空間分布圖

Figure 4 Spatial distribution of ESCI in Yuzhong County from 2008 to 2018

在分析ESV的空間格局的基礎(chǔ)上, 借助熱點(diǎn)分析, 選取具有統(tǒng)計(jì)意義的冷點(diǎn)和熱點(diǎn), 得到2008年和2018年榆中縣生ESV變化的熱點(diǎn)圖, 從而進(jìn)一步分析ESV空間上的差異性與集聚性。從圖5可知, 榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值總體上呈現(xiàn)“南熱中冷, 北無(wú)顯著”的分布特點(diǎn), 其中增值熱點(diǎn)區(qū)域主要位于南部的興隆山自然保護(hù)區(qū), 該區(qū)林地和草地面積占比大, 生態(tài)系統(tǒng)服務(wù)價(jià)值高, 冷點(diǎn)主要位于中部的和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、連搭鎮(zhèn)、城關(guān)鎮(zhèn)、小康營(yíng)鄉(xiāng)等, 是主要的城市建設(shè)區(qū), 建設(shè)用地面積大, 生態(tài)系統(tǒng)服務(wù)價(jià)值普遍偏低, 北部大部分區(qū)域變化不是特別明顯, 除此之外, 和平鎮(zhèn)北部、金崖鎮(zhèn)和夏官營(yíng)鎮(zhèn)中部、貢井鎮(zhèn)南部、青城鎮(zhèn)中部和北部、園子岔鄉(xiāng)西部和北部也有小范圍的熱點(diǎn)、次熱點(diǎn)分布, 主要是由于草地、水域的影響。2008—2018年榆中縣ESV熱點(diǎn)分布格局基本上保持穩(wěn)定, 冷點(diǎn)在和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、城關(guān)鎮(zhèn)有所增加, 主要是建設(shè)用地增加的影響。

表7 2008—2018年榆中縣ESV及其變化的空間自相關(guān)性

3 討論

基于土地利用變更數(shù)據(jù), 運(yùn)用當(dāng)量因子法, 對(duì)榆中縣的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行了測(cè)算, 榆中縣2018年的ESV為67.93億元, 這與周文霞[32]等對(duì)蘭州市生態(tài)系統(tǒng)服務(wù)價(jià)值(值為80.6億元)的研究, 以及吳翠霞等對(duì)白銀市生態(tài)系統(tǒng)服務(wù)價(jià)值[6](值為56.75億元)的研究結(jié)果接近, 說(shuō)明本文的研究結(jié)果應(yīng)該比較符合榆中縣的實(shí)際情況, 應(yīng)該能反映近年來(lái)榆中縣土地利用轉(zhuǎn)型引起的生態(tài)系統(tǒng)服務(wù)功能的變化。生態(tài)系統(tǒng)具有復(fù)雜性, 生態(tài)系統(tǒng)服務(wù)功能會(huì)受到諸多因素的影響, 生態(tài)系統(tǒng)服務(wù)價(jià)值的定量化測(cè)算一直是學(xué)術(shù)界研究的熱點(diǎn)和難點(diǎn)問(wèn)題, 已有研究雖然在不斷創(chuàng)新研究方法, 但生態(tài)系統(tǒng)服務(wù)價(jià)值的測(cè)算仍然存在許多現(xiàn)階段無(wú)法避免的主觀因素[33], 不同評(píng)價(jià)方法的測(cè)算結(jié)果可能會(huì)存在較大的差異[34-35][35][35]。因此, 基于當(dāng)量法對(duì)榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值測(cè)算的結(jié)果和其他方法的測(cè)算結(jié)果可能會(huì)存在一定的誤差, 如潘翔[36]等應(yīng)用InVEST模型和影子工程法計(jì)算了秦巴山(甘肅部分)25個(gè)集中連片特困區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值, 值為1525.05億元。

表8 2008—2018年榆中縣各土地利用類型ESV的空間自相關(guān)性

圖5 2008—2018年榆中縣ESV變化冷點(diǎn)和熱點(diǎn)的空間分布

Figure 5 Spatial distribution of cold spots and hot spots of ESV changes in Yuzhong County from 2008 to 2018

土地利用轉(zhuǎn)型包括顯性和隱性土地利用形態(tài)的變化。顯性形態(tài)變化指土地利用結(jié)構(gòu)和功能的變化[8] [9][10], 隱性形態(tài)變化包括權(quán)屬、質(zhì)量、經(jīng)營(yíng)方式等引起的土地利用形態(tài)的轉(zhuǎn)變。本文僅對(duì)對(duì)顯性形態(tài)的土地利變化做了分析, 對(duì)顯性形態(tài)的質(zhì)量、等級(jí)、強(qiáng)度等其他變化未做分析, 但實(shí)際上這些變化必然會(huì)影響到各類生態(tài)系統(tǒng)的數(shù)量、空間分布、結(jié)構(gòu)以及功能的改變[11]。所以僅從不同土地利用類型的數(shù)量角度出發(fā)對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行測(cè)算, 可能會(huì)存在一定的誤差。本文在測(cè)算生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)沒有考慮土地質(zhì)量、土地等級(jí)等其他因素, 這在今后的研究中還需要進(jìn)一步改進(jìn), 并且Cosatanza[5]等學(xué)者認(rèn)為, 建設(shè)用地也具有大氣調(diào)節(jié)、水分調(diào)節(jié)等生態(tài)功能, 在研究時(shí)應(yīng)該考慮建設(shè)用地的生態(tài)系統(tǒng)服務(wù)價(jià)值, 本研究中未對(duì)建設(shè)用地的生態(tài)系統(tǒng)服務(wù)價(jià)值做相應(yīng)的分析, 后續(xù)研究中還要進(jìn)行修正。

4 結(jié)論

(1)草地和耕地是榆中縣最主要的生態(tài)系統(tǒng), 占比將近約80%。研究期間, 榆中縣耕地、草地、林地、濕地和未利用地面積均有所減少, 但總體而言, 土地利用轉(zhuǎn)型時(shí)空變化較為平穩(wěn), 僅有0.8834%的土地利用發(fā)生了變化。其中耕地面積損失最多, 主要轉(zhuǎn)換為建設(shè)用地, 其次是草地, 主要轉(zhuǎn)為耕地和建設(shè)用地。相反, 建設(shè)用地和水域面積有所增加, 建設(shè)用地面積增加都主要來(lái)源于耕地和草地, 由于榆中縣城鎮(zhèn)化的發(fā)展, 使得縣域建設(shè)用地需求快速增加; 水域面積增加主要來(lái)源于耕地, 主要受到縣域水庫(kù)建設(shè)增加的影響。總的來(lái)看, 榆中縣土地利用轉(zhuǎn)化最為活躍的區(qū)域?yàn)榭h域中部平原區(qū)域的和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、連搭鎮(zhèn)、城關(guān)鎮(zhèn)、金崖鎮(zhèn)等鄉(xiāng)鎮(zhèn), 該區(qū)域城鎮(zhèn)化進(jìn)程和經(jīng)濟(jì)發(fā)展速度較快, 是蘭州市城市“東擴(kuò)”的重要方向。

(2)2008—2018年榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值總體上保持穩(wěn)定的狀態(tài), 有小幅度的下降趨勢(shì), 生態(tài)系統(tǒng)服務(wù)價(jià)值凈減少了0.41億元。從土地利用類型來(lái)看, 耕地、林地、草地、濕地、未利用地的ESV均有所減少, 但水域的ESV有所增加。草地和耕地生態(tài)系統(tǒng)服務(wù)價(jià)值減少是造成榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值下降的主要原因。從生態(tài)系統(tǒng)服務(wù)功能來(lái)看, 水資源供給功能的服務(wù)價(jià)值有所增加, 除此之外, 其他各功能的生態(tài)系統(tǒng)服務(wù)價(jià)值均呈下降趨勢(shì), 單項(xiàng)ESV以水文調(diào)節(jié)和氣候調(diào)節(jié)為主。近10年來(lái), 榆中縣生態(tài)系統(tǒng)服務(wù)價(jià)值變化不大, 說(shuō)明區(qū)域生態(tài)保護(hù)總體而言比較好。

(3)榆中縣ESV的空間分布總體上呈現(xiàn)南部高, 中部低, 北部較高的特征, 從ESCI的分布來(lái)看, ESV低值區(qū)域在和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)、城關(guān)鎮(zhèn)有比較明顯的擴(kuò)大, 這主要是由于這些區(qū)域城鎮(zhèn)化的加快, 導(dǎo)致區(qū)域建設(shè)用地?cái)?shù)量的增加。2008年和2018年ESV的全局Moran's指均大于0, 說(shuō)明兩個(gè)時(shí)期ESV的空間分布具有明顯的空間正相關(guān)性, 熱點(diǎn)分析表明ESV的高值區(qū)主要位于縣域南部, 低值區(qū)主要位于縣域中部, 這與榆中縣中部側(cè)重生產(chǎn)、生活以及南部側(cè)重生態(tài)保護(hù)的現(xiàn)狀相吻合。

[1] 龍花樓. 論土地利用轉(zhuǎn)型與鄉(xiāng)村轉(zhuǎn)型發(fā)展[J]. 地理科學(xué)進(jìn)展, 2012, 31(2): 131–138.

[2] 于倩茹, 徐煖銀, 孫思琦, 等. 鄂西北地區(qū)土地利用格局及其生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空變化[J]. 生態(tài)科學(xué), 2019, 38(4): 119–128.

[3] 史慧慧, 程久苗, 費(fèi)羅成, 等. 1990—2015年長(zhǎng)三角城市群土地利用轉(zhuǎn)型與生態(tài)系統(tǒng)服務(wù)功能變化[J]. 水土保持研究, 2019, 26(1): 301–307.

[4] 路昌, 韓笑, 韓會(huì)芹. 山東省土地利用轉(zhuǎn)型與生態(tài)服務(wù)功能交叉敏感性研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào), 2020, 51(10): 223-231.

[5] COSTANZA R, D'ARGE R, DE GROOT R, et al. The value of the world's ecosystem services and natural capital[J]. Nature, 1997, 387 (6630): 253–260.

[6] 吳翠霞, 馮永忠, 陶卓琳, 等. 白銀市土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2020, 36(11): 74–81.

[7] 張學(xué)斌, 羅君, 石培基, 等. 格網(wǎng)尺度下張掖生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空演變及地形梯度分異特征[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2020, 31(2): 543–553.

[8] 王權(quán), 李陽(yáng)兵, 黃娟, 等. 喀斯特槽谷區(qū)土地利用轉(zhuǎn)型過(guò)程對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[J]. 水土保持研究, 2019, 26(3): 192–198.

[9] CHENG Jie, HAN Jichang, WANG Huanyuan, et al. Background of land development and opportunity of land use transition[J]. Asian Agricultural Research, 2015, 7(12): 45–48.

[10] LIU Zhaoxu, LIU Liming. Characteristics and driving factors of rural livelihood transition in the east coastal region of China: A case study of suburban Shanghai[J]. Journal of Rural Studies, 2016, 43: 145–158.

[11] 孫思琦, 郭濼, 薛達(dá)元. 重慶市巫山縣土地利用格局及其生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空變化[J]. 生態(tài)科學(xué), 2019, 38(1): 176–185.

[12] 歐陽(yáng)志云, 王效科, 苗鴻. 中國(guó)陸地生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟(jì)價(jià)值的初步研究[J]. 生態(tài)學(xué)報(bào), 1999, 19(5): 607–613.

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

[14] 謝高地, 張彩霞, 張昌順, 等. 中國(guó)生態(tài)系統(tǒng)服務(wù)的價(jià)值[J]. 資源科學(xué), 2015, 37(9): 1740–1746.

[15] 劉永強(qiáng), 廖柳文, 龍花樓, 等. 土地利用轉(zhuǎn)型的生態(tài)系統(tǒng)服務(wù)價(jià)值效應(yīng)分析——以湖南省為例[J]. 地理研究, 2015, 34(4): 691–700.

[16] 戴文遠(yuǎn), 江方奇, 黃萬(wàn)里, 等. 基于“三生空間”的土地利用功能轉(zhuǎn)型及生態(tài)服務(wù)價(jià)值研究——以福州新區(qū)為例[J]. 自然資源學(xué)報(bào), 2018, 33(12): 2098–2109.

[17] 朱治州, 鐘業(yè)喜. 長(zhǎng)江三角洲城市群土地利用及其生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空演變研究[J]. 長(zhǎng)江流域資源與環(huán)境, 2019, 28(7): 1520–1530.

[18] 薛嵩嵩, 高凡, 何兵, 等. 近30年烏倫古河流域土地利用與生態(tài)系統(tǒng)服務(wù)價(jià)值變化研究[J]. 水土保持通報(bào), 2019, 39(6): 223–229.

[19] 傅伯杰, 張立偉. 土地利用變化與生態(tài)系統(tǒng)服務(wù):概念、方法與進(jìn)展[J]. 地理科學(xué)進(jìn)展, 2014, 33(4): 441–446.

[20] 李濤, 甘德欣, 楊知建, 等. 土地利用變化影響下洞庭湖地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空演變[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(12): 3787–3796.

[21] 郭榮中, 楊敏華. 長(zhǎng)株潭地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值分析及趨勢(shì)預(yù)測(cè)[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(5): 238-246.

[22] 劉春芳, 王川, 劉立程. 三大自然區(qū)過(guò)渡帶生境質(zhì)量時(shí)空差異及形成機(jī)制——以榆中縣為例[J]. 地理研究, 2018, 37(2): 419–432.

[23] 王秀蘭, 包玉海. 土地利用動(dòng)態(tài)變化研究方法探討[J]. 地理科學(xué)進(jìn)展, 1999(1): 3–5.

[24] 羅婭, 楊勝天, 劉曉燕, 等. 黃河河口鎮(zhèn)—潼關(guān)區(qū)間1998-2010年土地利用變化特征[J]. 地理學(xué)報(bào), 2014, 69(1): 42–53.

[25] 趙鴻雁, 陳英, 楊潔, 等. 基于改進(jìn)當(dāng)量的甘肅省耕地生態(tài)系統(tǒng)服務(wù)價(jià)值及其與區(qū)域經(jīng)濟(jì)發(fā)展的空間關(guān)系研究[J]. 干旱區(qū)地理, 2018, 41(4): 851–858.

[26] 李晶, 李紅艷, 張良. 關(guān)中—天水經(jīng)濟(jì)區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系[J]. 生態(tài)學(xué)報(bào), 2016, 36(10): 3053–3062.

[27] 錢彩云, 鞏杰, 張金茜, 等. 甘肅白龍江流域生態(tài)系統(tǒng)服務(wù)變化及權(quán)衡與協(xié)同關(guān)系[J]. 地理學(xué)報(bào), 2018, 73(5): 868–879.

[28] 公云龍, 張紹良, 章蘭蘭. 城市地價(jià)空間自相關(guān)分析——以宿州市為例[J]. 經(jīng)濟(jì)地理, 2011, 31(11): 1906– 1911.

[29] Getis A, ORD J K. The analysis of spatial association by use of distance statistics[J]. Geographic Analysis, 1992, 24: 189–206.

[30] 王蓓, 趙軍, 胡秀芳. 基于InVEST模型的黑河流域生態(tài)系統(tǒng)服務(wù)空間格局分析[J]. 生態(tài)學(xué)雜志, 2016, 35(10): 2783–2792.

[31] 雷金睿, 陳宗鑄, 陳小花, 等.1980—2018年海南島土地利用與生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空動(dòng)態(tài)變化[J]. 生態(tài)學(xué)報(bào), 2020, 40(14): 4760–4773.

[32] 周文霞, 石培基, 王永男, 等. 河谷型城市生態(tài)系統(tǒng)服務(wù)價(jià)值效應(yīng)——以蘭州為例[J]. 干旱區(qū)研究, 2017, 34(1): 232–241.

[33] 徐麗芬, 許學(xué)工, 羅濤, 等. 基于土地利用的生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量修訂方法——以渤海灣沿岸為例[J]. 地理研究, 2012, 31(10): 1775–1784.

[34] 萬(wàn)利, 陳佑啟, 譚靖, 等. 土地利用變化對(duì)區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值的影響——以北京市為例[J]. 地域研究與開發(fā), 2009, 28(4): 94–99.

[35]鄧楚雄, 鐘小龍, 謝炳庚, 等. 洞庭湖區(qū)土地生態(tài)系統(tǒng)的服務(wù)價(jià)值時(shí)空變化[J]. 地理研究, 2019, 38(4): 844–855.

[36] 潘翔, 吳娜, 徐中民, 等. 精準(zhǔn)扶貧視角下生態(tài)系統(tǒng)服務(wù)與多維貧困耦合協(xié)調(diào)差異:以甘肅秦巴山貧困核心區(qū)為例[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2020, 36(7): 879–889.

Spatio-temporal dynamics of land use transition and ecosystem service value in Yuzhong County

Qu Wen, CHEN Ying,XIE Baopeng*, PEI Tingting, REN Xijin, He Ying

College of Management, Gansu Agricultural University, Lanzhou 730070, China

Ecosystem service value is an important representation of regional ecological security. It aimed to identify dynamics of ecosystem services value caused by land use transformation, which had great significance for regional ecological environmental protection. Based on land use data in Yuzhong County of 2008 and 2018, this study applied ArcGIS software, value equivalence factor method, ecosystem service value change index, and spatial statistical analysis method to research the impact of land use transition on ecosystem service value dynamics in Yuzhong County. The results showed that during the period from 2008 to 2018, the area of cropland and grassland decreased, but the area of build-up land increased; the main modes of land transformation in Yuzhong County were cropland and grassland to build-up land, grassland to cropland. In this period, overall ESV in Yuzhong County showed a downward trend. Among them the ESV of grassland accounted for the largest proportion with 54%. The single ESV was dominated by hydrological regulation and climate regulation, which accounted for nearly 50% of the total ESV. The spatial distribution of ESV in Yuzhong County had obvious positive spatial correlation. The hot spot area of ESV was located in the south of Yuzhong County, and the cold spot area was located in the middle of Yuzhong County. The spatial distribution of ESV was characterized by high in the south and north, low in the middle. The low-value area of ESV was further spread in the middle of Yuzhong County. At last, the research results would provide scientific reference for regional ecological environment protection and land management.

land use transition; ecosystem service value; spatial autocorrelation; hotspot analysis; Yuzhong County

屈雯, 陳英, 謝保鵬, 等. 榆中縣土地利用轉(zhuǎn)型與生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空演變研究[J]. 生態(tài)科學(xué), 2022, 41(5): 197–207.

Qu Wen, CHEN Ying,XIE Baopeng, et al. Spatio-temporal dynamics of land use transition and ecosystem service value in Yuzhong County[J]. Ecological Science, 2022, 41(5): 197–207.

10.14108/j.cnki.1008-8873.2022.05.024

S157.2

A

1008-8873(2022)05-197-11

2020-09-07;

2020-11-01

甘肅農(nóng)業(yè)大學(xué)公招博士科研啟動(dòng)基金(GAU- KYQD-2017-34); 甘肅省高等學(xué)校科研項(xiàng)目(2019A-46)

屈雯(1995—), 女, 甘肅天水人, 碩士研究生, 主要從事土地資源管理研究, E-mail: 1914687138@qq.com

謝保鵬, 男, 講師, 主要從事土地資源管理研究, E-mail: xiebp@gsau.edu.cn