近20年南充市土地利用與生態(tài)系統(tǒng)服務(wù)價值時空變化

摘 要：在城鎮(zhèn)化推進(jìn)和退耕還林還草等政策影響下，南充市土地利用結(jié)構(gòu)發(fā)生顯著改變，基于2000—2019年3個時期土地利用數(shù)據(jù)，運(yùn)用價值當(dāng)量因子和空間統(tǒng)計分析等方法，探索南充市土地利用與生態(tài)系統(tǒng)服務(wù)價值（ESV）的時空變化。結(jié)果表明：（1）耕地是南充市主要土地利用類型，占82%以上，2000—2019年耕地和灌木呈減少趨勢，而林地、水域、草地和建設(shè)用地呈增長趨勢，建設(shè)用地（109.44%）和水域（40.62%）增長最為顯著。（2）2000—2019年南充市ESV總量先急速增長后小幅回落，凈增加73.8億元，增長率為22.89%，耕地ESV占38%以上；單項(xiàng)ESV以水文調(diào)節(jié)和土壤保持為主，兩者占總ESV的44%以上；各縣（市、區(qū)）ESV總量和單位面積ESV呈北高南低的特點(diǎn)。（3）2000—2019年3個時期ESV及其變化的Moran’s I指數(shù)均大于0，具有顯著自相關(guān)性；ESV高值區(qū)和低值區(qū)不斷擴(kuò)展，空間聚集性逐步增強(qiáng)；ESV減損區(qū)和變化冷點(diǎn)區(qū)集中于城鎮(zhèn)建設(shè)區(qū)域，增益區(qū)和變化熱點(diǎn)區(qū)集中于水域和林地覆蓋區(qū)域。研究結(jié)果可為南充市國土規(guī)劃和生態(tài)文明建設(shè)提供參考。

關(guān)鍵詞：生態(tài)系統(tǒng)服務(wù)價值；時空變化；政策驅(qū)動；價值當(dāng)量因子法；空間統(tǒng)計分析

中圖分類號：X171.1;F301.2"" 文獻(xiàn)標(biāo)志碼：A""" 文章編號：1673-5072（2023）02-0121-10

土地利用/覆被變化（Land-Use and Land-Cover Change，LUCC）是人類最基本的實(shí)踐活動［1］，對維持生態(tài)系統(tǒng)服務(wù)功能起決定作用［2］，是研究生態(tài)系統(tǒng)服務(wù)的關(guān)鍵切入點(diǎn)［3］。生態(tài)系統(tǒng)服務(wù)是人類賴以生存與發(fā)展的自然環(huán)境條件［4］，與人類福祉密切相關(guān)［5］，其可持續(xù)供給是社會可持續(xù)發(fā)展的基礎(chǔ)［6］。正確認(rèn)識生態(tài)系統(tǒng)服務(wù)之間關(guān)系，是開展多方位生態(tài)系統(tǒng)服務(wù)可持續(xù)管理決策的前提［7-8］。隨著1997年Costanza等［9］提出生態(tài)系統(tǒng)服務(wù)價值（Ecosystem service value，ESV）的估算方法及謝高地等［10］結(jié)合中國實(shí)際提出并修改中國單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量表以來，ESV已成為國內(nèi)外土地利用變化與生態(tài)系統(tǒng)間相互影響的研究熱點(diǎn)。目前眾多研究［11-12］均采用謝高地等［10］提出的中國生態(tài)系統(tǒng)服務(wù)價值當(dāng)量法進(jìn)行ESV估算；部分學(xué)者則充分考慮生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性，采用InVEST模型對各項(xiàng)生態(tài)系統(tǒng)服務(wù)估算［13-14］。價值當(dāng)量法因具有數(shù)據(jù)需求量少［7］、操作簡單、結(jié)果易于比較［15］等特點(diǎn)，被廣大學(xué)者運(yùn)用于流域［1］、島嶼［7］、城市群［16］、?。?7］等各級區(qū)域尺度的研究，成果豐碩，但這些研究往往忽視建設(shè)用地的作用，缺乏對建設(shè)用地ESV的評估。

南充是國家重要的商品糧和農(nóng)副產(chǎn)品生產(chǎn)基地，也是成渝經(jīng)濟(jì)圈中心城市之一，在經(jīng)濟(jì)建設(shè)和城鎮(zhèn)化的快速推進(jìn)過程中，生態(tài)系統(tǒng)的健康發(fā)展尤為重要。本文以南充市為研究對象，基于單位面積價值當(dāng)量因子法和空間統(tǒng)計分析，利用2000—2019年土地利用數(shù)據(jù)，從時空角度探索南充市土地利用和生態(tài)系統(tǒng)服務(wù)價值的動態(tài)變化特征，尤其是建設(shè)用地對生態(tài)系統(tǒng)服務(wù)價值的影響，以期為南充市生態(tài)文明建設(shè)和區(qū)域協(xié)同發(fā)展提供科學(xué)依據(jù)。

1 研究地區(qū)與研究方法

1.1 研究區(qū)概況

南充市位于四川省東北部，嘉陵江中游，地理位置介于30°27′—31°51′N，105°26′—106°57′E，總面積為12 477.20 km2，占全省面積的2.57%，下轄9個縣級行政區(qū)（3市轄區(qū)、5縣、1縣級市）。南充屬中亞熱帶濕潤季風(fēng)氣候區(qū)，四季分明，雨熱同期，降水集中于夏秋季節(jié)，年均溫為15.5～17.4 ℃，年均降水量為980～1 150 mm，地勢北高南低，北部為低山，南部為丘陵。南充2019年末常住人口為643.5萬人，地區(qū)生產(chǎn)總值（GDP）2 322.22億元，分別占四川省的7.68%和4.98%。21世紀(jì)以來，南充社會經(jīng)濟(jì)增長迅速，城鎮(zhèn)面積快速擴(kuò)張，土地利用變化加快，深刻影響著地區(qū)生態(tài)系統(tǒng)的結(jié)構(gòu)和功能。

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

文中2000年、2010年和2019年土地利用數(shù)據(jù)來源于武漢大學(xué)遙感與信息工程學(xué)院黃昕教授團(tuán)隊［18］基于Landsat衛(wèi)星衍生的中國年度土地覆蓋產(chǎn)品（https：//doi.org/10.5281/zenodo.4417810），空間分辨率為30 m，按照該產(chǎn)品分類體系，南充市土地利用類型劃分為耕地、林地、灌木、草地、水域、裸地和建設(shè)用地。2010年南充市糧食產(chǎn)量和播種面積數(shù)據(jù)來源于四川省統(tǒng)計局公布的《2011年四川省統(tǒng)計年鑒》（https：//tjj.sc.gov.cn/scstjj/c105855/nj.shtml）。

1.3 土地利用動態(tài)度分析

土地利用動態(tài)度可反映一定時間內(nèi)土地利用變化的區(qū)域差異［19］，采取土地利用動態(tài)度分析2000—2019年南充市各類土地利用類型的變化情況。

1.4 生態(tài)系統(tǒng)服務(wù)價值評估

以謝高地等［10］制定的生態(tài)系統(tǒng)服務(wù)價值當(dāng)量表為基礎(chǔ)，確定了中國1個生態(tài)系統(tǒng)服務(wù)價值當(dāng)量的經(jīng)濟(jì)價值為3 406.50元·hm-2，同期全國和南充市糧食單產(chǎn)分別為4 973.24 kg·hm-2和5 589.00 kg·hm-2，即南充市1個生態(tài)系統(tǒng)服務(wù)價值當(dāng)量為3 828.27元·hm-2。參考雷金睿等［7］在海南島和冉玉菊等［20］在滇中城市群的研究成果，結(jié)合研究區(qū)實(shí)際情況構(gòu)建南充市單位面積生態(tài)系統(tǒng)服務(wù)價值系數(shù)表（表1），生態(tài)系統(tǒng)服務(wù)價值計算引用Costanza等［5］的生態(tài)系統(tǒng)服務(wù)價值系數(shù)法。

1.5 生態(tài)系統(tǒng)服務(wù)價值損益

結(jié)合土地利用轉(zhuǎn)移矩陣測算LUCC對ESV的影響，用1 km×1 km漁網(wǎng)計算并分析ESV的空間損益情況［21］。

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

全局空間自相關(guān)可探索要素的整體空間分布狀況［22］，通常用Moran’s I指數(shù)與Z度量［17］，本文采用Moran’s I對南充市生態(tài)系統(tǒng)服務(wù)價值格網(wǎng)數(shù)據(jù)進(jìn)行全局空間自相關(guān)檢驗(yàn)。局部空間自相關(guān)旨在研究鄰近區(qū)域間對某屬性值的集聚程度［23］，可檢驗(yàn)局部區(qū)域的空間關(guān)聯(lián)性和異質(zhì)性，本文采用Moran’s Ii對南充市生態(tài)系統(tǒng)服務(wù)價值進(jìn)行局部空間自相關(guān)檢驗(yàn)。采用熱點(diǎn)分析 （Getis-Ord Gi*）探索ESV空間變化的聚集和分異特征［24］，即“冷點(diǎn)”和“熱點(diǎn)”分布格局。

2 結(jié)果與分析

2.1 南充市土地利用動態(tài)變化特征

2.1.1 土地利用面積變化

根據(jù)2000—2019年南充市土地利用類型變化情況來看（表2，圖1），南充市土地利用以耕地為主，主要分布于地勢較低平區(qū)域；其次是林地，主要分布于北部、西部和西南部山地；其余地類面積占比較少。從面積變化看，2000—2019年耕地面積變化量最大，累計減少97 070.57 hm2；其次為林地，增加77 632.09 hm2；建設(shè)用地和水域面積變化也較大，分別增加12 780.84 hm2和6 570.18 hm2；灌木、裸地和草地面積變化不大。從面積變化率看，2000—2019年建設(shè)用地增長109.44%，林地增長81.43%，水域增長40.62%，耕地減少8.63%，灌木、草地和裸地占比極低，變化率不予重點(diǎn)考慮。從面積時間變化來看，耕地面積持續(xù)減少，其中2000—2010年迅速下降；林地面積在2000—2010年急速增加，2010—2019年有所減少；建設(shè)用地和水域持續(xù)增加，其中建設(shè)用地在2010—2019年增加迅猛；其余土地利用類型面積變化不大。

2.1.2 土地利用變化動態(tài)度

深入分析2000—2019年南充市土地利用動態(tài)度（表3），從時間段來看，2000—2010年林地、灌木、草地和水域動態(tài)度在4%以上，增長較快，耕地和裸地動態(tài)度小于0，建設(shè)用地動態(tài)度為1.82%；2010—2019年耕地、林地、灌木動態(tài)度小于0，草地、裸地和建設(shè)用地動態(tài)度在3%以上，其中建設(shè)用地擴(kuò)張迅速，達(dá)7.71%；從2000—2019年整體來看，除耕地和灌木動態(tài)度小于0外，其余地類動態(tài)度均呈增長趨勢，其中林地、草地、

水域和建設(shè)用地動態(tài)度分別為4.29%、6.47%、2.14%和5.76%。從主要地類面積變化動態(tài)度來看，耕地動態(tài)度均小于0，其面積先迅速減少后減幅變緩；2000—2010年和2010—2019年林地動態(tài)度分別為8.81%和-0.35%，林地先迅猛增加后小幅減少；2000—2010年和2010—2019年水域動態(tài)度分別為4.03%和0.03%，水域先快速增加后緩慢擴(kuò)張；2000—2010年和2010—2019年建設(shè)用地動態(tài)度分別為1.82%和7.71%， 2010—2019年建設(shè)用地快速擴(kuò)張。

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

從南充市2000—2019年土地利用轉(zhuǎn)移矩陣來看（表4），耕地轉(zhuǎn)出面積最大，轉(zhuǎn)為林地、建設(shè)用地和水域的面積分別為99 351.14 hm2、16 681.78 hm2和3 881.07 hm2；林地有21 965.10 hm2轉(zhuǎn)為耕地，建設(shè)用地有4 068.35 hm2轉(zhuǎn)為水域，水域有913.22 hm2轉(zhuǎn)為耕地，其余地類面積轉(zhuǎn)出不大。從轉(zhuǎn)入情況來看，林地轉(zhuǎn)入面積主要來自耕地和水域；耕地主要由林地和水域轉(zhuǎn)入；建設(shè)用地主要由耕地和水域轉(zhuǎn)入；水域則主要由建設(shè)用地和耕地轉(zhuǎn)入；灌木、草地和水域轉(zhuǎn)入面積較小。

從土地利用轉(zhuǎn)移空間分布來看（圖2），2000—2019年南充市土地利用轉(zhuǎn)變集中于嘉陵江流域、區(qū)縣城區(qū)周邊和山地地區(qū)，其中區(qū)縣城區(qū)周邊因城鎮(zhèn)建設(shè)導(dǎo)致土地利用變化迅速，主要表現(xiàn)為耕地轉(zhuǎn)為建設(shè)用地，市轄三區(qū)尤為劇烈。嘉陵江流域因生態(tài)建設(shè)的推進(jìn)，導(dǎo)致大量建設(shè)用地轉(zhuǎn)為水域和林地，山地地區(qū)則因退耕還林而促使大面積耕地轉(zhuǎn)為林地。

2.2 南充市生態(tài)系統(tǒng)服務(wù)價值動態(tài)變化特征

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

2000—2019年，南充市ESV先急速增長后有所回落（表5），累計增長738 004.79萬元，增長率達(dá)22.89%。從各地類ESV來看，林地、水域、草地和裸地ESV呈增長趨勢，與其面積變化一致，其中水域和林地ESV分別增長315 939.57萬元和629 462.52萬元，分別增長40.62%和81.43%；耕地和建設(shè)用地分別減少149 016.45萬元和58 763.12萬元，變化率分別為-8.63%和109.44%。林地和水域面積增加是2000—2019年南充市ESV總體增加的主要原因，建設(shè)用地的擴(kuò)張則使2010—2019年ESV有所回落。

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

從南充市單項(xiàng)生態(tài)系統(tǒng)服務(wù)價值可看出（表6），3個時期中水文調(diào)節(jié)占比均最大，約30%左右，土壤保持約為15%。2000—2019年，單項(xiàng)ESV變化量最大的為水文調(diào)節(jié)，增加36.19億元，增長率達(dá)40.65%；氣候調(diào)節(jié)增加17.67億元，增長率為44.69%。僅食物生產(chǎn)和水資源供給ESV呈下降趨勢，分別減少2.15億元和0.68億元，增長率為-5.65%和-17.76%。水文調(diào)節(jié)和氣候調(diào)節(jié)價值量增加推動南充市ESV總體增長，得益于水域和林地面積的增長，而過快的城鎮(zhèn)建設(shè)導(dǎo)致水資源供給價值量下降，對ESV增長起阻礙作用。

2.2.3 生態(tài)系統(tǒng)服務(wù)價值空間分布和損益變化

為深入分析南充市ESV及其損益的空間分布，運(yùn)用自然斷點(diǎn)法制作2000—2019年南充市ESV及其損益的空間分布圖（圖3）。

從2000年和2019年ESV分布圖來看，南充市嘉陵江和升鐘水庫等水體區(qū)域ESV極高以及區(qū)縣城鎮(zhèn)用地周邊極低外，呈現(xiàn)出中部低邊緣高的特點(diǎn)，這與南充市耕地和林地分布情況極為契合，即中部地區(qū)地勢較和緩，耕地面積廣，ESV較低，邊遠(yuǎn)地區(qū)地勢起伏較大，林地多，ESV較高。從2000—2019年ESV空間變化來看，高值區(qū)不斷擴(kuò)張，低值區(qū)逐漸萎縮，但區(qū)縣城區(qū)周邊隨建設(shè)用地的迅速增長，ESV極低值范圍在快速擴(kuò)展。

從2000—2019年ESV的損益來看，各縣市區(qū)城鎮(zhèn)建設(shè)區(qū)ESV減損較大，市轄三區(qū)形成連片極高損失區(qū)，源于水域和耕地大面積轉(zhuǎn)為建設(shè)用地；山地和水域覆蓋地區(qū)因林地和水域面積增加，成為ESV增益區(qū)，東河、西河與嘉陵江交匯處和高坪區(qū)北部嘉陵江河段等水域形成極高增益區(qū)，源于水域面積的增長。

2.2.4 各縣（市、區(qū)）生態(tài)系統(tǒng)服務(wù)價值變化

從行政單位尺度對南充市ESV進(jìn)行分析，從空間分布來看（圖4），各縣（市、區(qū)）ESV呈北高南低的特點(diǎn)；閬中、南部、儀隴、營山4地所轄面積較大，林地面積較多，2000—2019年ESV貢獻(xiàn)率達(dá)64%以上；市轄三區(qū)面積較小且建設(shè)用地集中，ESV總量較低。從時間變化來看，2000—2010年，各地ESV均明顯增加，閬中、南部增長迅猛；2010—2019年，閬中、南部、西充、儀隴4地有所回落，其中閬中ESV減少最多。

采用自然斷點(diǎn)法進(jìn)一步分析南充市各縣（市、區(qū)）單位面積ESV及其時空變化情況（圖5），從空間分布來看，2000—2019年南充市單位面積ESV分布大體呈北高南低的特點(diǎn)，閬中和南部為高值區(qū)域，西充和順慶為低值區(qū)域。從時間變化來看，2000—2010年嘉陵、高坪、蓬安、營山、儀隴5地單位面積ESV整體偏低，ESV上升速率低于整體水平；2010—2019年高坪、蓬安2地單位面積ESV有所上漲，生態(tài)改善明顯；整體而言，2010—2019年各地單位面積ESV均隨ESV增加而不同程度上升，生態(tài)在逐步改善。

2.3 南充市生態(tài)系統(tǒng)服務(wù)價值空間自相關(guān)分析

2.3.1 全局空間自相關(guān)分析

從南充市ESV全局自相關(guān)檢驗(yàn)結(jié)果來看（表7），2000—2019年3個時期Moran’s I均大于0，P值均小于0.001，說明南充市ESV空間分布具有較強(qiáng)的正相關(guān)性，空間上具有極顯著的聚集性。 I值先增后降，表明ESV的空間自相關(guān)性先增強(qiáng)后減弱，說明2000—2010年退耕還林等宏觀政策促使ESV空間聚集性增強(qiáng)，2010—2019年在城鎮(zhèn)建設(shè)等經(jīng)濟(jì)活動影響下，ESV聚集性有所回落。從2000—2010年和2010—2019年ESV動態(tài)變化的空間自相關(guān)性來看，全局Moran’s I呈下降趨勢，但也具有顯著的空間聚集性，說明南充市ESV動態(tài)變化過程的空間聚集性在減弱，2000—2019年ESV總體動態(tài)變化的I值為0.199（P＜0.001），具有極顯著空間聚集效應(yīng)。

2.3.2 局部空間自相關(guān)分析

從南充市ESV局部自相關(guān)檢驗(yàn)來看（圖6），2000—2019年ESV高值聚集區(qū)主要分布于嘉陵江沿岸、升鐘水庫等水域和邊緣山區(qū)，其范圍逐漸擴(kuò)張；ESV低值區(qū)分布于城鎮(zhèn)建設(shè)密集區(qū)，隨著城市化推進(jìn)，低值區(qū)范圍在不斷變大，ESV空間分異加劇。

2.3.3 熱點(diǎn)分析

從南充市ESV變化熱點(diǎn)圖來看（圖7），2000—2010年ESV增值熱點(diǎn)集中于市轄三區(qū)嘉陵江河段、蓬安縣西北部—閬中市東北部嘉陵江河段、升鐘水庫，退耕區(qū)域也有零星分布，主要由于耕地和建設(shè)用地轉(zhuǎn)化為水域和林地，局部ESV升高；ESV損失冷點(diǎn)區(qū)域分散于城鎮(zhèn)建設(shè)和林地退化區(qū)域，冷點(diǎn)面積小于熱點(diǎn)區(qū)。2010—2019年ESV增值熱點(diǎn)集中分布于閬中東北部解元水庫、北部嘉陵江河段、高坪區(qū)北部和順慶區(qū)東部嘉陵江河段、儀隴縣北部前進(jìn)水庫、同盟水庫等，水庫修建或生態(tài)建設(shè)推動水域面積擴(kuò)張區(qū)域；ESV損失冷點(diǎn)則分布于各地城鎮(zhèn)建設(shè)區(qū)域。從2000—2019年整體變化來看，與南充市土地利用轉(zhuǎn)移情況較為一致，ESV增值熱點(diǎn)分布于新修水庫庫區(qū)、水域面積擴(kuò)張河段和退耕區(qū)域等水域、林地面積增長地區(qū)；ESV損失冷點(diǎn)分布于城鎮(zhèn)化加速推進(jìn)，建設(shè)用地增加和林地減少的地區(qū)；ESV損失冷點(diǎn)面積小于增值熱點(diǎn)，印證了南充市土地利用的總體ESV在增加。

3 討 論

進(jìn)入新世紀(jì)以來，在國家退耕還林還草政策全面實(shí)施、成渝雙城經(jīng)濟(jì)圈建設(shè)推進(jìn)、人口流失等因素影響下，南充市土地利用結(jié)構(gòu)顯著變化，進(jìn)而引起生態(tài)系統(tǒng)服務(wù)價值的變化。林地和水域能創(chuàng)造極高的生態(tài)服務(wù)價值［22］，得益于水域和林地面積的增加，2000—2019年南充市ESV增加738 004.79萬元，此結(jié)果與羅芳等［17］對四川省生態(tài)系統(tǒng)服務(wù)價值的研究結(jié)論相符，應(yīng)鞏固退耕成果，推動南充市林草資源穩(wěn)步增長，加強(qiáng)對嘉陵江流域和升鐘水庫等水域的管理，改善水域生態(tài)結(jié)構(gòu)。建設(shè)用地雖對ESV起一定促進(jìn)作用［25］，但對于水資源調(diào)節(jié)、氣體調(diào)節(jié)和凈化環(huán)境等起阻礙作用，今后需科學(xué)規(guī)劃城市建設(shè)并整治違建濫建等現(xiàn)象，切實(shí)保護(hù)耕地和水域資源?？傮w而言，南充生態(tài)環(huán)境在持續(xù)改善，仍需落實(shí)生態(tài)文明建設(shè)，優(yōu)化土地利用結(jié)構(gòu)，促進(jìn)生態(tài)系統(tǒng)可持續(xù)健康發(fā)展。

除土地利用變化外，人口分布、經(jīng)濟(jì)發(fā)展水平、氣候條件和產(chǎn)業(yè)結(jié)構(gòu)布局等自然和社會經(jīng)濟(jì)因素也影響區(qū)域ESV。本文僅從土地利用變化來探討南充市ESV的時空變化具有一定局限性，今后開展相關(guān)研究時需重點(diǎn)關(guān)注自然和社會經(jīng)濟(jì)因子對ESV時空變化的影響、ESV的模型估算及預(yù)測模擬研究，以提高ESV估算的準(zhǔn)確性。

參考文獻(xiàn)：

［1］ 孫慧蘭，李衛(wèi)紅，陳亞鵬，等.新疆伊犁河流域生態(tài)服務(wù)價值對土地利用變化的響應(yīng)［J］.生態(tài)學(xué)報，2010，30（4）：887-894.

［2］ 樊宏，張建平.岷江上游半干旱河谷土地利用/土地覆蓋研究［J］.中國沙漠，2002（3）：72-77.

［3］ 韓蕊，孫思琦，郭濼，等.川東地區(qū)生態(tài)系統(tǒng)服務(wù)價值時空演變及其驅(qū)動力分析［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報，2019，35（9）：1136-1143.

［4］ BOYD J，BANZHAF S.What are ecosystem services？ The need for standardized environmental accounting units［J］.Ecological Economics，2007，63（2-3）：616-626.

［5］ COSTANZA R，GROOT R，SUTTON P，et al.Changes in the global value of ecosystem services［J］.Global Environmental Change，2014，26：152-158.

［6］ GROOT R，BRANDER L，PLOEG S，et al.Global estimates of the value of ecosystems and their services in monetary units［J］.Ecosystem Services，2012，1（1）：50-61.

［7］ 雷金睿，陳宗鑄，陳小花，等.1980—2018年海南島土地利用與生態(tài)系統(tǒng)服務(wù)價值時空變化［J］.生態(tài)學(xué)報，2020，40（14）：4760-4773.

［8］ 錢彩云，鞏杰，張金茜，等.甘肅白龍江流域生態(tài)系統(tǒng)服務(wù)變化及權(quán)衡與協(xié)同關(guān)系［J］.地理學(xué)報，2018，73（5）：868-879.

［9］ COSTANZA R，DARGE R，GROOT R，et al.The value of the world’s ecosystem services and natural capital［J］.Nature，1997，387：253-260.

［10］謝高地，張彩霞，張雷明，等.基于單位面積價值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價值化方法改進(jìn)［J］.自然資源學(xué)報，2015，30（8）：1243-1254.

［11］任鈺覃，吳勇.閬中市土地利用變化對生態(tài)系統(tǒng)服務(wù)價值影響［J］.西華師范大學(xué)學(xué)報（自然科學(xué)版），2020，41（4）：404-412.

［12］羅芳，潘安，陳忠升，等.龍門山斷裂帶土地利用變化及生態(tài)系統(tǒng)服務(wù)價值損益評估［J］.西華師范大學(xué)學(xué)報（自然科學(xué)版），2021，42（4）：417-425.

［13］景海超，劉穎慧，賀佩，等.青藏高原典型區(qū)生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性及其影響因素：以那曲市為例［J］.生態(tài)學(xué)報，2022，42（7）：2657-2673.

［14］廖珍梅，楊薇，蔡宴朋，等.大清河-白洋淀流域生態(tài)功能評價及分區(qū)初探［J］.環(huán)境科學(xué)學(xué)報，2022，42（1）：131-140.

［15］孫夢華，牛文浩，張蚌蚌，等.黃河流域土地利用變化下生態(tài)系統(tǒng)服務(wù)價值時空演變及其響應(yīng)：以陜甘寧地區(qū)為例［J］.應(yīng)用生態(tài)學(xué)報，2021，32（11）：3913-3922.

［16］李魁明，王曉燕，姚羅蘭，等.京津冀地區(qū)生態(tài)系統(tǒng)服務(wù)價值時空變化及驅(qū)動因子分析［J］.環(huán)境工程技術(shù)學(xué)報，2022，12（4）：1114-1122.

［17］羅芳，潘安，陳忠升，等.四川省土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響研究［J］.云南農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)），2021，36（4）：734-744.

［18］YANG J ，HUANG X.The 30 m annual land cover dataset and its dynamics in China from 1990 to 2019［J］.Earth System Science Data，2021，13（8）：3907-3925.

［19］王超，常勇，侯西勇，等.基于土地利用格局變化的膠東半島生境質(zhì)量時空演變特征研究［J］.地球信息科學(xué)學(xué)報，2021，23（10）：1809-1822.

［20］冉玉菊，雷冬梅，劉林，等.滇中城市群2000—2020年土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響［J］.水土保持通報，2021，41（4）：310-322.

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

［22］喬斌，祝存兄，曹曉云，等.格網(wǎng)尺度下青海瑪多縣土地利用及生態(tài)系統(tǒng)服務(wù)價值空間自相關(guān)分析［J］.應(yīng)用生態(tài)學(xué)報，2020，31（5）：1660-1672.

［23］陳睿，楊燦，楊艷，等.洞庭湖生態(tài)經(jīng)濟(jì)區(qū)生態(tài)系統(tǒng)服務(wù)價值的時空演繹及其驅(qū)動因素［J］.應(yīng)用生態(tài)學(xué)報，2021，33（1）：169-179.

［24］郭椿陽，高尚，周伯燕，等.基于格網(wǎng)的伏牛山區(qū)土地利用變化對生態(tài)服務(wù)價值影響研究［J］.生態(tài)學(xué)報，2019，39（10）：3482-3493.

［25］董家華，舒廷飛，謝慧，等.城市建設(shè)用地生態(tài)服務(wù)功能價值計算與應(yīng)用［J］.同濟(jì)大學(xué)學(xué)報（自然科學(xué)版），2007，35（5）：636-640.

Spatio-temporal Changes of Land Use and Ecosystem Service Valueof Nanchong City in Recent 20 Years

ZHOU Jian-wei1ac，2，WANG Lei1bc，2，YANG Hai-qing1bc，2，DENG Qing-chun1bc，2，LIU Hui1bc，2，LUO Jun1bc，2

（1.a.College of Life Science，b.School of Geographical Sciences，c.Sichuan Provincial Engineering Laboratory of Monitoring and

Control for Soil Erosion on Dry Valleys，China West Normal University，Nanchong Sichuan 637009，China；2.Liangshan Soil Erosion and Ecological Restoration in Dry Valleys Observation and Research Station，Xide Sichuan 616753，China）

Abstract：Under the policy influence of urbanization promotion and returning farmland to forest and grassland，the land use structure of Nanchong has changed significantly.Based on the land use data of three periods from 2000 to 2019，the methods of value equivalent factor and spatial statistical analysis are employed to explore the spatio-temporal changes of land use and ecosystem service value （ESV） of Nanchong.The results show as follows：（1） the main land use type of Nanchong is cultivated land，accounting for more than 82%；the cultivated land and shrub had shown a decreasing trend from 2000 to 2019 whereas woodland，water area，grassland and construction land had shown an increasing trend，among which the growth rates of construction land （109.44%） and water area（40.62%） was the most significant；（2） from 2000 to 2019，the total ESV of Nanchong had increased rapidly at first and then decreased slightly，with a net increase of 7.38 billion yuan，a growth rate of 22.89%；the ESV of cultivated land had accounted for more than 38%； hydrological regulation and soil conservation were dominant in single ESV，accounting for more than 44% of the total ESV；the total ESV and ESV per unit area of each county （city，district） were high in the north and low in the south；（3） ESV and its changes of Moran’s I index in the three periods from 2000 to 2019 were significantly autocorrelated； the high value region and low value region of ESV had continued to expand，and the spatial aggregation was gradually enhanced；ESV loss area and variation cold spot area were concentrated in the urban construction area，while gain area and variation hot spot area were centralized in water area and woodland covered area.The research results have provided references for the land planning and ecological civilization construction in Nanchong.

Keywords：ecosystem service value；spatio-temporal changes；policy driving；value equivalent factor method； spatial statistical analysis