宜賓市近20年土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空效應(yīng)

何治學(xué)　王磊　羅君　劉輝　鄧青春　楊海青

摘 要：土地利用變化和生態(tài)系統(tǒng)服務(wù)價(jià)值一直是生態(tài)學(xué)研究的重點(diǎn)之一。本文基于2000、2010和2020年宜賓市遙感影像圖，采用GIS制圖和面積轉(zhuǎn)移矩陣等方法分析了2000—2020年宜賓市土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空效應(yīng)。結(jié)果表明：（1）近20年來，宜賓市土地利用類型轉(zhuǎn)換面積以耕地、林地和建設(shè)用地為主，建設(shè)用地面積增速呈現(xiàn)上升的趨勢(shì)，2000—2010年與2010—2020年土地利用綜合動(dòng)態(tài)度分別為0.32%和0.31%；（2）生態(tài)系統(tǒng)服務(wù)價(jià)值在時(shí)間序列上呈現(xiàn)先降低后回升的趨勢(shì)，總價(jià)值從490.64×108元增長(zhǎng)到491.24×108元，其中林地、耕地和水域的貢獻(xiàn)率最大；在空間分布上，單元格生態(tài)系統(tǒng)服務(wù)價(jià)值主要以中和較高等級(jí)區(qū)域?yàn)橹?，低等?jí)區(qū)域主要集中于宜賓市主城區(qū)和南溪區(qū)（建設(shè)用地為主）。（3）生態(tài)系統(tǒng)服務(wù)價(jià)值聚類分布具有顯著正相關(guān)性，H-H聚類與H-L聚類的生態(tài)系統(tǒng)服務(wù)價(jià)值增值區(qū)域主要分布在耕地轉(zhuǎn)為林地的區(qū)域，L-L聚類與L-H聚類的減值區(qū)域則主要分布在林地轉(zhuǎn)為耕地與建設(shè)用地的區(qū)域。本研究通過直觀清晰的展現(xiàn)土地利用與生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空變化，有助于揭示兩種變化之間的關(guān)系，并為指導(dǎo)生態(tài)環(huán)境建設(shè)提供科學(xué)依據(jù)。

關(guān)鍵詞：宜賓市；土地利用變化；生態(tài)服務(wù)價(jià)值；時(shí)空效應(yīng)；空間統(tǒng)計(jì)

中圖分類號(hào)：K903 文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1673-5072（2023）04-0350-08

土地利用和生態(tài)系統(tǒng)服務(wù)價(jià)值（Ecological Service Value，ESV）是對(duì)一個(gè)地區(qū)生態(tài)環(huán)境質(zhì)量最為直觀的評(píng)估［1-3］。ESV的形成受到多種因素的影響，其時(shí)空動(dòng)態(tài)變化過程與生態(tài)結(jié)構(gòu)和功能密切相關(guān)［4-5］。

土地是人類活動(dòng)最重要的載體，在進(jìn)行ESV評(píng)估時(shí)，土地利用的變化可以直觀地表現(xiàn)人類活動(dòng)的規(guī)模和強(qiáng)度［6］，是影響生態(tài)系統(tǒng)結(jié)構(gòu)和功能最直接的驅(qū)動(dòng)因素之一［3］。土地利用/覆被變化（Land-Use and Land-Cover Change，LUCC）會(huì)顯著影響區(qū)域ESV以及生態(tài)敏感性的波動(dòng)［7-8］。同時(shí)，不同地區(qū)的地形、地貌和氣候等自然條件決定了地區(qū)的經(jīng)濟(jì)發(fā)展和生態(tài)治理方向，從而直接影響ESV的變化［7-8］。Costanza等［9］基于16個(gè)生物群落估算出了17個(gè)生態(tài)系統(tǒng)服務(wù)的當(dāng)前經(jīng)濟(jì)價(jià)值并制定了生態(tài)服務(wù)價(jià)值化評(píng)估方法，但該評(píng)估方法過于廣泛，不同國(guó)家在直接運(yùn)用時(shí)經(jīng)常會(huì)出現(xiàn)低估或者忽視了某些ESV的情況，從而導(dǎo)致研究結(jié)果的誤差較大。于是，Xie等［10-11］在Costanza等［9］的基礎(chǔ)上，通過大量走訪和問卷調(diào)查得出了新的生態(tài)系統(tǒng)服務(wù)評(píng)估體系，并在中國(guó)得到了廣泛的應(yīng)用［12-14］。部分研究人員充分考慮生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性，采用InVEST模型進(jìn)行各項(xiàng)生態(tài)系統(tǒng)服務(wù)估算［15-16］。土地利用和ESV研究作為生態(tài)學(xué)研究中的重要一環(huán)受到研究者們的廣泛關(guān)注。這類研究在國(guó)家［17］、流域［18］、城市群［19-20］等各尺度區(qū)域上成果豐富，但對(duì)于建設(shè)用地的生態(tài)作用研究往往受到了忽視，建設(shè)用地的ESV評(píng)估需要進(jìn)一步深入。

由于特殊的地理區(qū)位，宜賓市成為國(guó)家“五縱七橫”交通規(guī)劃建設(shè)中長(zhǎng)江東西軸線與南北干線的交匯點(diǎn)，是建設(shè)長(zhǎng)江上游沿江經(jīng)濟(jì)帶和川南經(jīng)濟(jì)新增長(zhǎng)極的重要支撐城市。本研究以宜賓市作為研究區(qū)，研究目的包括：（1）揭示2000—2020年宜賓市土地利用和ESV時(shí)空變化規(guī)律；（2）探究土地利用變化與ESV變化的空間可視化特征，尤其是建設(shè)用地變化與ESV的關(guān)系。本研究以期為生態(tài)文明城市的建設(shè)和發(fā)展提供科學(xué)的參考依據(jù)，為制定合理的區(qū)域生態(tài)環(huán)境管理政策奠定理論和實(shí)踐基礎(chǔ)。

1 研究區(qū)概況

宜賓市（103°36′—105°20′E，27°50′—29°16′N，圖1）位于四川盆地南緣，地處川、滇、黔三省交界處，因金沙江和岷江在此處匯合成長(zhǎng)江，而被稱為“萬里長(zhǎng)江第一城”。宜賓市總面積達(dá)13 283 km2，下轄3區(qū)7縣，總?cè)丝诩s551.5萬人［21］。該市海拔236～2 008 m，地形呈現(xiàn)西南高、東北低的態(tài)勢(shì)，地貌以中低山地和丘陵為主；水系以長(zhǎng)江為主脈，支流共600余條；屬中亞熱帶濕潤(rùn)季風(fēng)氣候，年均氣溫18 ℃，年均降水量1 050～1 618 mm，雨季集中在5—10月，具有氣候溫和、雨量充沛、四季分明等特點(diǎn)。

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

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

宜賓市2000和2010年的土地利用類型數(shù)據(jù)來源于全國(guó)地理信息資源目錄服務(wù)系統(tǒng)（http：//www.webmap.cn）提供的30 m分辨率全球地表覆蓋數(shù)據(jù)，2020年土地利用數(shù)據(jù)來源于地理國(guó)情監(jiān)測(cè)云平臺(tái)（http：//www.dsac.cn），根據(jù)該數(shù)據(jù)集的分類體系，通過重分類將數(shù)據(jù)集劃分為耕地、林地、灌木、草地、水域、未利用地和建設(shè)用地（圖1）。糧食播種面積及其凈利潤(rùn)數(shù)據(jù)集來源于四川省糧食和物資儲(chǔ)備局（https：//lwj.sc.gov.cn/sclwj/index.shtml）和《中國(guó)統(tǒng)計(jì)年鑒》。

2.2 研究方法

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

土地利用轉(zhuǎn)移矩陣全面具體地反映某一區(qū)域土地利用變化的結(jié)構(gòu)特征［22］。本文基于宜賓市2000—2020年的土地利用類型變化數(shù)據(jù)，運(yùn)用ENVI 5.3制作土地利用轉(zhuǎn)移矩陣［23］。

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

土地利用動(dòng)態(tài)度是反應(yīng)土地利用類型變化程度的指數(shù)，可分為單一動(dòng)態(tài)度和綜合動(dòng)態(tài)度。單一動(dòng)態(tài)度指特定區(qū)域一定時(shí)間范圍內(nèi)某種土地利用類型變化情況［24］，綜合動(dòng)態(tài)度指整個(gè)區(qū)域內(nèi)所用土地利用的年際變化［25］。

2.2.3 ESV核算

本研究參考謝高地等［26］所修訂的ESV當(dāng)量因子表，將單位面積（1 hm2）的稻谷、小麥和玉米三大糧食的凈利潤(rùn)作為1個(gè)標(biāo)準(zhǔn)當(dāng)量因子的ESV量。根據(jù)宜賓市實(shí)際情況進(jìn)行ESV系數(shù)修正，標(biāo)準(zhǔn)當(dāng)量因子計(jì)算公式為［26］：D=Sr×Fr×Sw×Fw+Sc×Fc，式中，D為1個(gè)標(biāo)準(zhǔn)當(dāng)量因子的ESV量（元·hm-2）；Sr、Sw、Sc分別為研究區(qū)稻谷、小麥和玉米播種面積占全國(guó)三種作物播種總面積百分比（%）；Fr、Fw、Fc分別為稻谷、小麥和玉米單位面積平均凈利潤(rùn)（元·hm-2）?；诠接?jì)算得到宜賓市2000—2020年的標(biāo)準(zhǔn)當(dāng)量因子為2 524.1元·hm-2。VCk=D×Ek，式中，VCk為k類土地利用類型單位面積的ESV（元·hm-2）；Ek為k類土地利用類型單位面積ESV當(dāng)量，各土地利用類型的值詳見文獻(xiàn)［26］的表1。通過以上公式，得到了宜賓市單位面積ESV表（表1）。

再運(yùn)用表1進(jìn)行ESV計(jì)算，公式為［27］：∑Ak×VCk，式中，ESV為生態(tài)系統(tǒng)服務(wù)價(jià)值（元）；Ak：第k類土地利用類型分布面積（hm2）。然后通過運(yùn)用GIS軟件中的漁網(wǎng)工具，以2 km×2 km格網(wǎng)計(jì)算單元格ESV，并結(jié)合幾何間隔法和自然間斷點(diǎn)法將研究區(qū)ESV劃分成低［0，5×104］、較低（5×104，10×104］、中（10×104，15×104］、較高（15×104，30×104］、高（30×104，∞］5個(gè)等級(jí)，從而制成單元格ESV空間分布圖。

2.3 空間自相關(guān)性分析

運(yùn)用全局莫蘭指數(shù)（I）［14］和局部莫蘭指數(shù)（Ii）［14］描述宜賓市不同單元格ESV分布格局的空間自相關(guān)性。

3 結(jié)果與分析

3.1 土地利用類型變化特征

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

2000—2020年宜賓市土地利用變化的統(tǒng)計(jì)結(jié)果圖2所示。在這20年間，退耕還林還草和開墾工程同時(shí)進(jìn)行，導(dǎo)致耕地和林地面積發(fā)生了巨大變化。耕地主要轉(zhuǎn)出為林地（94 654.98 hm2）和建設(shè)用地（12 015.90 hm2），占耕地轉(zhuǎn)出總面積的85.71%。耕地的轉(zhuǎn)入面積共107 957.97 hm2，主要源自林地（106 327.62 hm2），因?yàn)榉鲐毠ぷ鞯拇罅﹂_展，山區(qū)林地不斷被開墾，變化區(qū)域以屏山縣與筠連縣為主；其次源自水域（1 383.66 hm2），只有少數(shù)來源于灌木、草地和建設(shè)用地，但草地轉(zhuǎn)出占比最大的是耕地，占47.31%。

3.1.2 土地利用動(dòng)態(tài)度變化

由表2可以看出，2000—2010年，耕地變化量達(dá)42 216.84 hm2，林地達(dá)-48 456.72 hm2，土地利用單一動(dòng)態(tài)度中最大的是灌木，達(dá)到61.53%，最小為水域，僅0.40%，綜合動(dòng)態(tài)度為0.32%；2010—2020年。耕地變化量為-45 032.31 hm2，林地為36 398.70 hm2，土地利用單一動(dòng)態(tài)度中最大的是草地，達(dá)23.40%，其次是灌木，達(dá)-7.59%，綜合動(dòng)態(tài)度為0.31%。

3.2 ESV變化特征

3.2.1 ESV時(shí)空變化

由表3可知，2000—2020年宜賓市總ESV從490.64×108元增長(zhǎng)到491.24×108元，增長(zhǎng)率達(dá)到0.12%。2000—2010年，宜賓市開墾工程大規(guī)模進(jìn)行，導(dǎo)致林地和草地面積大量減少，生態(tài)服務(wù)價(jià)值也隨之急劇降低了21.32×108元。2010—2020年，宜賓市主打退耕還林還草和長(zhǎng)江流域沿岸生態(tài)修復(fù)戰(zhàn)略，促使生態(tài)環(huán)境得到了大幅度的恢復(fù)，ESV較2000年有一定的提升。

從不同土地利用類型來看，耕地和灌木的總ESV和各項(xiàng)服務(wù)的ESV皆呈現(xiàn)先增后減的趨勢(shì)，年均變化率分別是-0.02%和2.76%；林地和草地呈現(xiàn)先減后增的態(tài)勢(shì)，年均變化率分別為-0.10%和1.70%；水域總ESV持續(xù)增長(zhǎng)，年均增長(zhǎng)率為0.72%。從不同服務(wù)功能來看，ESV中占比最高的是調(diào)節(jié)服務(wù)，達(dá)到61.88%，最低的是文化服務(wù)，僅占3.90%。各項(xiàng)生態(tài)系統(tǒng)服務(wù)類型總價(jià)值皆呈現(xiàn)先減后增的趨勢(shì)，其中，調(diào)節(jié)服務(wù)和支持服務(wù)價(jià)值變化最明顯，年均變化率分別為0.04%和-0.07%。調(diào)節(jié)服務(wù)主要來源于林地，其次是水域，這說明調(diào)節(jié)服務(wù)在以林地和水域?yàn)橹鞯囊速e市生態(tài)系統(tǒng)服務(wù)中占主導(dǎo)地位。

由圖 3可知，宜賓市單元格ESV主要以中和較高等級(jí)區(qū)域?yàn)橹?，面積占比達(dá)到60%以上。2000—2010年，中和較高等級(jí)面積比例減小，低和較低等級(jí)面積占比增加；2010—2020年，中、較高和高等級(jí)面積比例皆增加。西北部和南部的山區(qū)以及三江（長(zhǎng)江、金沙江和岷江）流域沿岸地區(qū)的ESV較高，以較高和高等級(jí)區(qū)域?yàn)橹?；中等?jí)區(qū)域主要分布在高縣、珙縣和興文縣；宜賓市主城區(qū)和南溪區(qū)是低等級(jí)的主要集中區(qū)域。

3.2.2 ESV空間自相關(guān)性

由表4可知，宜賓市2000—2020年各年份及年份差的ESV全局莫蘭指數(shù)均大于0（P<0.01），說明宜賓市ESV空間分布上具有顯著正相關(guān)性。2000—2020年莫蘭指數(shù)呈現(xiàn)持續(xù)下降的趨勢(shì)，說明2000—2020年宜賓市ESV的空間相關(guān)性在不斷減小。原因在于近年來宜賓市為提升自身的競(jìng)爭(zhēng)力對(duì)土地利用的結(jié)構(gòu)進(jìn)行了整改，加快宜賓市區(qū)南岸、三江新區(qū)和南溪區(qū)的城市化建設(shè)進(jìn)度以及羅龍工業(yè)園區(qū)的建設(shè)進(jìn)度，調(diào)整了各土地利用類型的面積與分布格局，從而導(dǎo)致ESV分散化。

如圖4所示：總體來說，宜賓市以高-高值（H-H）和低-低值（L-L）聚類為主。H-H聚類主要分布在北部和西北部，以屏山縣境內(nèi)金沙江流域和翠屏區(qū)為主；L-L聚類主要分布在南部和東部，以敘州區(qū)、高縣和珙縣行政區(qū)交界區(qū)域、筠連縣、長(zhǎng)寧縣、江安縣和興文縣交界區(qū)域?yàn)橹鳎桓?低值（H-L）聚類零散分布在L-L聚類范圍中或附近；低-高值（L-H）聚類主要分布在H-H聚類范圍內(nèi)，且以屏山縣為主。

3.3 土地利用變化與ESV變化聚類的疊加特征

將宜賓市2000—2020年土地利用變化（圖1）與ESV變化聚類（圖4）結(jié)合分析，結(jié)果如圖5所示：HH聚類與HL聚類的ESV增值區(qū)域具有相似的分布情況，主要分布在耕地轉(zhuǎn)為林地的區(qū)域。ESV的增長(zhǎng)主要來源于耕地轉(zhuǎn)為林地，增長(zhǎng)了45.59×108元；其次來源于耕地轉(zhuǎn)為水域，增長(zhǎng)了10.55×108萬元，水域面積的增加對(duì)生態(tài)系統(tǒng)服務(wù)功能起到極大的改善作用，僅耕地和林地轉(zhuǎn)為水域就促進(jìn)ESV增長(zhǎng)了11.02×108元。該區(qū)域具有優(yōu)質(zhì)的生態(tài)環(huán)境質(zhì)量，應(yīng)減少進(jìn)行規(guī)劃建設(shè)活動(dòng)，合理控制退耕還林還草工程的推行進(jìn)度，促使ESV保持穩(wěn)定的增長(zhǎng)速度。

LL聚類與LH聚類的ESV減值區(qū)域主要分布在林地轉(zhuǎn)為耕地和建設(shè)用地的區(qū)域。ESV的減少主要來源于林地轉(zhuǎn)為耕地，減少了51.21×108元；其次來源于水域轉(zhuǎn)為耕地，減少了4.25×108元，建設(shè)用地面積的增加以及除建設(shè)用地外的其他土地利用類型轉(zhuǎn)為耕地時(shí)則會(huì)造成生態(tài)系統(tǒng)服務(wù)功能的降低。該區(qū)域生態(tài)環(huán)

境質(zhì)量處于持續(xù)劣化的階段，應(yīng)盡快控制因城市化進(jìn)度過快而導(dǎo)致建設(shè)用地面積迅速增長(zhǎng)的情況，大力開展植樹造林等活動(dòng)，推動(dòng)生態(tài)修復(fù)工程的進(jìn)行。

4 討 論

近年來，宜賓市人口迅速增長(zhǎng)、鄉(xiāng)村振興、退耕還林還草和城市化戰(zhàn)略大力推進(jìn)等社會(huì)經(jīng)濟(jì)因素共同造成了土地利用類型的極速變化［4］。2000—2020年宜賓市土地利用類型的變化主要集中在耕地、林地和建設(shè)用地，平均土地利用綜合動(dòng)態(tài)度為0.32%。2000—2020年宜賓市ESV整體呈現(xiàn)先降低后上升的趨勢(shì)，ESV增長(zhǎng)了6.12×107元，林地、水域和耕地的貢獻(xiàn)率占主要地位，總計(jì)在80%以上。ESV的變化與耕地、林地、水域和建設(shè)用地面積變化具有顯著的相關(guān)關(guān)系，林地面積的增加是推動(dòng)ESV上升的主要?jiǎng)恿Γ?9，27］。雖然林地和水域面積的增加促使ESV的上升，但同時(shí)也導(dǎo)致建筑和耕地等土地類型面積受到了擠壓，導(dǎo)致被擠壓區(qū)域生態(tài)環(huán)境受到更嚴(yán)重的破壞，使得宜賓市生態(tài)環(huán)境的平衡性受到了沖擊，不利于生態(tài)環(huán)境的可持續(xù)發(fā)展，從而形成ESV的空間異質(zhì)性［20］。因此，應(yīng)該注重建設(shè)用地與各土地類型的協(xié)調(diào)發(fā)展［28］，充分考慮建設(shè)用地對(duì)ESV的積極和消極作用，和諧各土地類型間的相互影響，減緩宜賓市生態(tài)環(huán)境受到的沖擊，從而有效的改善宜賓市生態(tài)環(huán)境，提高ESV，推動(dòng)生態(tài)文明城市建設(shè)進(jìn)程。

從ESV時(shí)空分布圖可以看出，在空間分布上，宜賓市ESV在2000—2010年呈現(xiàn)東北部和南部地區(qū)降低的趨勢(shì)，在2010—2020年呈現(xiàn)全市上升的態(tài)勢(shì)，且西部和南部地區(qū)增長(zhǎng)量最大。在時(shí)間序列上，ESV高值區(qū)域面積先減少后增加。因此，在制定生態(tài)保護(hù)政策時(shí)需根據(jù)各區(qū)域?qū)嶋H情況因地制宜，把握各區(qū)域政策力度，從而最大限度提升各區(qū)域生態(tài)環(huán)境質(zhì)量。且宜賓市ESV在空間聚類分布上表現(xiàn)出以僅增加或僅減少為主，導(dǎo)致宜賓市生態(tài)環(huán)境的穩(wěn)定性在城市發(fā)展過程中受到了沖擊，ESV主要源于調(diào)節(jié)服務(wù)和支持服務(wù)，貢獻(xiàn)率總計(jì)70%以上，各項(xiàng)服務(wù)功能間存在著重要的相互作用力［2，8］。因此，后續(xù)研究應(yīng)進(jìn)一步探究各項(xiàng)服務(wù)功能之間的關(guān)系，促進(jìn)各項(xiàng)服務(wù)間協(xié)調(diào)可持續(xù)發(fā)展，從而保護(hù)和維持生態(tài)環(huán)境的穩(wěn)定性。

參考文獻(xiàn)：

［1］ 魏慧，趙文武，張驍，等.基于土地利用變化的區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)價(jià)：以山東省德州市為例［J］.生態(tài)學(xué)報(bào)，2017，37（11）：3830-3839.

［2］ 嚴(yán)恩萍，林輝，王廣興，等.1990—2011年三峽庫(kù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值演變及驅(qū)動(dòng)力［J］.生態(tài)學(xué)報(bào)，2014，34（20）：5962-5973.

［3］ 陳佑啟，楊鵬.國(guó)際上土地利用/土地覆蓋變化研究的新進(jìn)展［J］.經(jīng)濟(jì)地理，2001，21（1）：95-100.

［4］ 劉桂林，張落成，張倩.長(zhǎng)三角地區(qū)土地利用時(shí)空變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.生態(tài)學(xué)報(bào)，2014，34（12）：3311-3319．

［5］ LI Z G，SUN Z S，TIAN Y J.Impact of land use/cover change on Yangtze River delta urban agglomeration ecosystem services value：temporal-spatial patterns and cold/hot spots ecosystem services value change brought by urbanization［J］.International Journal of Environmental Research and Public Health，2019，16（1）：123-123.

［6］ TASSER E，SCHURPKE U，ZODERER B M，et al.Towards an integrative assessment of land-use type values from the perspective of ecosystem services［J］.Ecosystem Services，2020，42：101082.

［7］ 吳后建，王學(xué)雷，寧龍梅，等.土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響：以武漢市為例［J］.長(zhǎng)江流域資源與環(huán)境，2006，15（2）：185-190.

［8］ 唐啟琳，劉方，劉秀明等.基于LUCC的喀斯特山區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)價(jià)［J］.環(huán)境科學(xué)與技術(shù)，2019，42（1）：170-177.

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

［10］XIE G D，LU C X，LENG Y F，et al.Ecological assets valuation of the Tibetan Plateau［J］.Journal of Natural Resources，2003，18（2）：189-196.

［11］XIE G D，ZHEN L，LU C X，et al.Expert knowledge based valuation method of ecosystem services in China［J］.Journal of Natural Resources，2008，23（5）：911-919.

［12］WANG W J，GUO H C，CHUAI X W，et al.The impact of land use change on the temporospatial variations of ecosystemsservices value in China and an optimized land use solution ［J］.Environmental Science & Policy，2014，44：62-72.

［13］趙永華，張玲玲，王曉峰.陜西省生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估及時(shí)空差異［J］.應(yīng)用生態(tài)學(xué)報(bào)，2011，22（10）：2662-2672.

［14］幸贊品，顏長(zhǎng)珍，馮坤，等.1975—2015年甘肅省白龍江流域自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值及其時(shí)空差異［J］.中國(guó)沙漠，2019，39（3）：172-182.

［15］李婷，劉康，胡勝，等.基于InVEST模型的秦嶺山地土壤流失及土壤保持生態(tài)效益評(píng)價(jià)［J］.長(zhǎng)江流域資源與環(huán)境，2014，23（9）：1242-1250.

［16］黃博強(qiáng)，黃金良，李迅，等.基于GIS和InVEST模型的海岸帶生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空動(dòng)態(tài)變化分析：以龍海市為例［J］.海洋環(huán)境科學(xué)，2015，34（6）：916-924.

［17］王兵，魯紹偉.中國(guó)經(jīng)濟(jì)林生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估［J］.應(yīng)用生態(tài)學(xué)報(bào)，2009，20（2）：417-425.

［18］粟曉玲，康紹忠，佟玲.內(nèi)陸河流域生態(tài)系統(tǒng)服務(wù)價(jià)值的動(dòng)態(tài)估算方法與應(yīng)用：以甘肅河西走廊石羊河流域?yàn)槔跩］.生態(tài)學(xué)報(bào)，2006，26（6）：2011-2019.

［19］劉桂林，張落成，張倩.長(zhǎng)三角地區(qū)土地利用時(shí)空變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.生態(tài)學(xué)報(bào)，2014，34（12）：3311-3319.

［20］楊青，劉耕源.濕地生態(tài)系統(tǒng)服務(wù)價(jià)值能值評(píng)估：以珠江三角洲城市群為例［J］.環(huán)境科學(xué)學(xué)報(bào)，2018，38（11）：4527-4538.

［21］宜賓市年鑒編纂委員會(huì).宜賓市年鑒［M］.成都：四川科學(xué)技術(shù)出版社，2017.

［22］周書貴，邵全琴，曹?。?0年黃土高原土地利用/覆被變化特征分析［J］.地球信息科學(xué)學(xué)報(bào)，2016，18（2）：190-199．

［23］吳琳娜，楊勝天，劉曉燕，等.1976年以來北洛河流域土地利用變化對(duì)人類活動(dòng)程度的響應(yīng)［J］.地理學(xué)報(bào)，2014，69（1）：54-63.

［24］王秀蘭，包玉海.土地利用動(dòng)態(tài)變化研究方法探討［J］.地理科學(xué)進(jìn)展，1999，18（1）：83-89.

［25］宋開山，劉殿偉，王宗明，等.1954年以來三江平原土地利用變化及驅(qū)動(dòng)力［J］.地理學(xué)報(bào)，2008，63（1）：93-104.

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

［27］彭文甫，周介銘，羅懷良，等.城市土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值損益估算：以成都市為例［J］.水土保持研究，2011，18（4）：43-51+27.

［28］羅繼文，周禧，黃亞南，等.南沙區(qū)土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.華南師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，54（3）：100-110.

Spatial-temporal Effects of Land Use Change and Ecosystem Service Valueof Yibin City in Recent 20 Years

HE Zhi-xue1a，2 WANG Lei1bc，2，LUO Jun1bc，2，LIU Hui1bc，2，DENG Qing-chun1bc，2，YANG Hai-qing1c，2

（1.a.College of Life Science，b.Sichuan Provincial Engineering Laboratory of Monitoring and Control for Soil Erosion onDry Valleys，c.School of Geographical Sciences，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：Land use change and ecosystem service value have always been one of the ecological research focuses.In this paper，the spatial-temporal effects of land use change and ecosystem service value in Yibin city from 2000 to 2020 are analyzed by GIS mapping and area transfer matrix，which is based on remote sensing images of Yibin city in 2000，2010 and 2020.The results are as follows：（1）In the past 20 years，the transfer area of land use type in Yibin city has been dominated by cropland，forestland and construction land，with the growth rate of construction land showing an upward trend；the comprehensive dynamic values of land use from 2000 to 2010 and from 2010 to 2020 were 0.32% and 0.31%，respectively；（2） The value of ecosystem service has presented a trend of decteasing first and then increasing in the time series；the total value has increased from 490.64 × 108 yuan to 491.24 × 108 yuan，to which the forest land，water area and cultivated land have contributed the most；in the spatial distribution，the ecosystem service value of the unit is mostly centered on the middle and the higher grade areas and the low-grade areas are mainly concentrated in the main urban areas of Yibin and Nanxi District （mainly construction land）；（3） The cluster distribution of ecosystem service value has a significant positive correlation；the value-added areas of H-H Cluster and H-L Cluster are mainly distributed in the areas transfered from cropland to forestland，while the value-decreased areas of L-L cluster and L-H cluster are mainly distributed in the areas transfered from forestland to cropland and construction land.The visual and clear presentation of the temporal and spatial changes in this study is helpful to reveal the relationship between land use and ecosystem service value，and provide scientific basis for guiding the construction of ecological environment.

Keywords：Yinbin city；land use change；ecological service value；spatial-temporal effects；spatial statistics

收稿日期：2022-06-09基金項(xiàng)目：西華師范大學(xué)科研啟動(dòng)項(xiàng)目（18Q020，18Q033）

作者簡(jiǎn)介：何治學(xué)（1997—），男，碩士研究生，主要從事水土保持與荒漠化防治研究。

通信作者：羅君（1987—），男，博士，講師，主要從事水土保持與荒漠化防治研究。Email：luojunxmx@126.com

引文格式：何治學(xué)，王磊，羅君，等.宜賓市近20年土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空效應(yīng)［J］.西華師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2023，44（4）：350-357.［HE Z X，WANG L，LUO J，et al.Spatial-temporal effects of land use change and ecosystem service value of Yibin city in recent 20 years［J］.Journal of China West Normal University （Natural Sciences），2023，44（4）：350-357.］