周建偉　羅君

摘要：基于2000～2020年高精度土地利用數(shù)據(jù)，利用生態(tài)系統(tǒng)服務(wù)價值（ESV）和景觀生態(tài)風(fēng)險指數(shù)（LERI）評估模型、雙變量空間自相關(guān)模型、地理探測器，探究安寧河流域和白龍江流域ESV和LERI的時空演變、關(guān)聯(lián)和分異特征。結(jié)果表明：安寧河流域和白龍江流域土地利用以林草為主；城鎮(zhèn)擴張影響下兩地ESV和LERI 均呈降低趨勢，ESV以低和中等等級為主，LERI以極低和低等級為主；ESV和LERI具有顯著空間負相關(guān)性且兩者關(guān)系以高價值-低風(fēng)險和低價值-高風(fēng)險為主；人類活動對兩地ESV和LERI空間分異起主導(dǎo)作用；林地和草地是兩流域ESV的重要來源，建設(shè)用地的激增使兩流域ESV不斷減少、LERI有所降低。需因地制宜完善用地結(jié)構(gòu)，加強林草及水域保護，合理控制城鎮(zhèn)擴張，推動流域土地利用結(jié)構(gòu)可持續(xù)發(fā)展。

關(guān) 鍵 詞：生態(tài)系統(tǒng)服務(wù)價值； 景觀生態(tài)風(fēng)險指數(shù)； 干旱河谷； 時空演變； 時空分異

中圖法分類號： X171

文獻標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2023.04.013

0 引 言

隨著經(jīng)濟社會的快速發(fā)展，人類活動對自然環(huán)境的影響加?。?］，致使生態(tài)系統(tǒng)結(jié)構(gòu)與功能惡化［2］、區(qū)域生態(tài)風(fēng)險不斷上升［3］等一系列生態(tài)問題日益嚴重。全面準(zhǔn)確把握生態(tài)系統(tǒng)功能的變化、科學(xué)防范生態(tài)風(fēng)險是區(qū)域生態(tài)環(huán)境保護和生態(tài)安全格局構(gòu)建的前提和關(guān)鍵［4］。

生態(tài)系統(tǒng)服務(wù)價值（Ecosystem Service Value，ESV）和景觀生態(tài)風(fēng)險指數(shù)（Landscape Ecological Risk Index，LERI）評估是生態(tài)環(huán)境評價的重要內(nèi)容，同生態(tài)安全評價息息相關(guān)［1］。生態(tài)系統(tǒng)服務(wù)指人類通過生態(tài)系統(tǒng)的結(jié)構(gòu)、過程和功能直接或間接獲得的生命支持產(chǎn)品和服務(wù)［5］，是人類賴以生存和發(fā)展的自然條件和效用［6］，對其進行價值評估是聯(lián)系生態(tài)系統(tǒng)研究與管理決策的關(guān)鍵環(huán)節(jié)［7］，有助于使自然資本在公共決策中得到認可。當(dāng)前，ESV評估方法主要分為單位面積服務(wù)價值功能價格法［8］和單位面積價值當(dāng)量因子法［9］，當(dāng)量因子法因適用于區(qū)域ESV評估［10］而被學(xué)者廣泛應(yīng)用于各級尺度研究［11-12］。景觀生態(tài)風(fēng)險指自然或人為因素干擾下景觀格局和生態(tài)過程相互作用可能產(chǎn)生的不利后果［13-14］，LERI評價可為區(qū)域綜合風(fēng)險防范提供決策依據(jù)，能有效引導(dǎo)區(qū)域景觀格局優(yōu)化與管控［15］。ESV和LERI的綜合運用，有利于使人類福祉與生態(tài)環(huán)境變化相聯(lián)系，可為區(qū)域生態(tài)環(huán)境保護提供決策支持［16］。

前人針對ESV和LERI各自進行了大量研究［13，17］，也開展了ESV和LERI空間關(guān)聯(lián)的相關(guān)研究［1，4］。但當(dāng)前ESV和LERI空間關(guān)聯(lián)和空間分異研究局限于單一區(qū)域或單一時段，缺乏不同區(qū)域及多時段動態(tài)對比研究，不僅無法對比不同區(qū)域ESV和LERI的空間關(guān)聯(lián)性，更無法全面表征ESV和LERI 空間分異性的動態(tài)變化。沿安寧河和白龍江分布的狹長干旱地段作為中國干旱河谷的重要組成部分，是區(qū)域社會、經(jīng)濟和文化的中心。安寧河流域是攀西城市群經(jīng)濟發(fā)展的核心地帶，白龍江流域是成渝城市群與城市群間的重要能源樞紐，兩地同為中國西部重要的多民族聚居區(qū)，對區(qū)域經(jīng)濟可持續(xù)發(fā)展及國家長治久安意義重大。在城鎮(zhèn)化進程加快、西部大開發(fā)、退耕還林還草等因素影響下，安寧河流域和白龍江流域經(jīng)濟增長迅速，但隨之而來的人口激增、城鎮(zhèn)過度擴張等問題不僅加劇區(qū)域人地矛盾突出及生態(tài)退化，也威脅中國生態(tài)文明建設(shè)和生態(tài)安全格局構(gòu)建。研究安寧河流域和白龍江流域ESV和LERI時空演變、二者空間關(guān)聯(lián)性及分異性特征，可為中國西部干旱河谷區(qū)土地利用規(guī)劃、生態(tài)修復(fù)、生態(tài)安全格局構(gòu)建及人地關(guān)系協(xié)調(diào)發(fā)展提供科學(xué)參考。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

安寧河為雅礱江一級支流，流經(jīng)四川省涼山州和攀枝花市（見圖1），流域面積約11 053.59 km2，海拔997～5 081 m，地處橫斷山脈東緣，屬亞熱帶季風(fēng)氣候區(qū)，夏季高溫多雨、冬季干燥少雨［18］；流域內(nèi)安寧河平原為四川省第二大平原，是川西南地區(qū)重要的經(jīng)濟增長極，2020年末GDP約960.51億元、人口約192.61萬人、城鎮(zhèn)化率約52.65%。白龍江為嘉陵江一級支流，流經(jīng)甘肅省甘南州、隴南市及四川省阿壩州和廣元市（見圖1），流域面積約32 113.07 km2，海拔450～4 865 m，處青藏高原向黃土高原、秦巴山區(qū)及暖溫帶向北亞熱帶過渡地帶［19］，夏季高溫多雨、冬季溫暖少雨，2020年末GDP約510.07億元、人口約182.11萬人、城鎮(zhèn)化率約37.72%。白龍江流域為川、甘交界處的重要多民族聚居區(qū)，流域內(nèi)梯級水電開發(fā)對中國能源安全起重要作用。

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

2000～2020年3期30 m分辨率土地利用數(shù)據(jù)來源于地球大數(shù)據(jù)科學(xué)工程數(shù)據(jù)共享服務(wù)系統(tǒng)數(shù)據(jù)（http：∥data.casearth.cn/），該數(shù)據(jù)由中國科學(xué)院空天信息創(chuàng)新研究院利用Landsat衛(wèi)星（Lnadsat TM，ETM+和OLI）所生產(chǎn)，總體精度為82.5%［20］，Kappa系數(shù)達0.784，滿足研究需求。地理探測因子數(shù)據(jù)包括源自中國科學(xué)院資源環(huán)境數(shù)據(jù)中心（http：∥www.resdc.cn/）的歸一化植被指數(shù)、年降水、年均溫、地均GDP和人口密度數(shù)據(jù)；高程數(shù)據(jù)來自地理空間數(shù)據(jù)云（http：∥www.gscloud.cn/），并由高程計算得到坡度和坡向數(shù)據(jù)；人為影響指數(shù)利用公式計算而得。糧食產(chǎn)量和價格數(shù)據(jù)源于研究期內(nèi)兩流域主要縣級行政區(qū)統(tǒng)計年鑒、四川省糧食和物資儲備局（https：∥lwj.sc.gov.cn/）和甘肅省糧食和物資儲備局（http：∥lscb.gansu.gov.cn/）。

1.3 研究方法

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

基于ESV評估模型［5］，參照謝高地等［21］提出的中國陸地生態(tài)系統(tǒng)單位面積ESV當(dāng)量（見表1），依據(jù)研究期內(nèi)安寧河流域和白龍江流域主要縣級行政區(qū)平均糧食產(chǎn)量（分別為5 389.86，4 088.20 kg/hm2）及2020年四川省和甘肅省主要農(nóng)作物收購價格（分別為2.41，3.15元/kg），按“單位面積農(nóng)田ESV等于平均糧食單產(chǎn)市場經(jīng)濟價值的1/7”［21］，得出安寧河流域和白龍江流域ESV當(dāng)量因子經(jīng)濟價值（分別為1 855.65，1 839.69元/hm2），最終得到研究區(qū)ESV總量，公式為

式中：ESV表示生態(tài)系統(tǒng)服務(wù)價值總量；Ai表示評價格網(wǎng)第i種地類的面積；Sij表示第i種地類第j種生態(tài)系統(tǒng)服務(wù)類型單位面積所提供的ESV當(dāng)量；i為地類數(shù)量；j為生態(tài)系統(tǒng)服務(wù)類型。

1.3.2 景觀生態(tài)風(fēng)險指數(shù)評估模型

從區(qū)域生態(tài)系統(tǒng)的景觀結(jié)構(gòu)出發(fā)，選取景觀干擾度指數(shù)、景觀脆弱度指數(shù)、景觀損失度指數(shù)構(gòu)建綜合生態(tài)風(fēng)險指數(shù)，以此分析研究區(qū)生態(tài)風(fēng)險大小及其變化情況，計算公式［22］如下：

式中：LERIk為第k個風(fēng)險小區(qū)的景觀生態(tài)風(fēng)險指數(shù)；Aki為第k個風(fēng)險小區(qū)內(nèi)第i種地類的面積；Ak為第k個風(fēng)險小區(qū)的面積；n為評價格網(wǎng)數(shù)；Ri為第i類地類的景觀損失度指數(shù)。各參數(shù)計算公式及含義［23-25］詳見表2。

1.3.3 人為影響指數(shù)

人為影響指數(shù)（Human Active Index，HAI）能表征人為活動對土地利用和景觀結(jié)構(gòu)變化的影響［26］，故利用HAI分析研究區(qū)ESV和LERI同人類干擾強度間的相互關(guān)系。公式如下：

式中：HAI為人為影響指數(shù)；Ai為評價格網(wǎng)內(nèi)第i種地類的總面積；Si為第i種地類的人為影響強度參數(shù)，基于前人研究［26-27］對Si進行賦值；n為地類數(shù)量；AT為評價格網(wǎng)總面積。

1.3.4 雙變量空間自相關(guān)模型

為探索研究區(qū)ESV和LERI的空間關(guān)聯(lián)特征，運用GeoDa雙變量空間分析模型，利用全局空間自相關(guān)Morans I指數(shù)反映兩者整體空間關(guān)聯(lián)和差異狀況，計算公式［28］如下：

式中：Iel為單位面積生態(tài)系統(tǒng)服務(wù)價值e和生態(tài)風(fēng)險指數(shù)l的雙變量全局自相關(guān)系數(shù)；yie和yil分別為第i個評價小區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價值和生態(tài)風(fēng)險指數(shù)；σe和σl為方差；n為評價小區(qū)數(shù)；wij為基于空間鄰接關(guān)系建立的權(quán)重矩陣。

1.3.5 地理探測器

地理探測器（Geodetector）可探索要素的空間分異特征［29］，采用其因子探測器對研究區(qū)ESV和LERI空間分異的驅(qū)動因子進行探測，公式為

式中：q表示自變量對因變量空間分異的解釋力；σh2和σ2分別為評價小區(qū)和全區(qū)Y值的方差；Nh和N分別為評價小區(qū)h和全區(qū)的樣本數(shù)。

2 結(jié)果與分析

2.1 土地利用動態(tài)變化分析

2000～2020年兩個典型干旱河谷土地利用均以林地為主，其次為草地和耕地，其余地類占比較小（見圖2），安寧河流域與白龍江流域的林草地占比分別在89.74%和84.78%以上、耕地占9.26%和13.15%以上。20 a來除安寧河流域林地和草地、白龍江流域耕地和草地有所減少外，其余地類均呈增長趨勢（見表3），其中安寧河流域草地縮減3.58%（113.95 km2）、建設(shè)用地增長激增85.73%（69.02 km2），白龍江流域耕地減少4.63%（144.31 km2）、建設(shè)用地擴張120.59%（116.51 km2）。

2.2 ESV動態(tài)變化分析

2000～2020年安寧河流域和白龍江流域ESV呈緩慢下降趨勢，分別減少0.61%和0.07%，兩地林草地對ESV的貢獻率突出，分別達91.58%和95.19%以上（見表3）。近20 a，安寧河流域林地、白龍江流域耕地以及兩地草地和建設(shè)用地ESV呈減少趨勢，其余地類ESV不斷增加，其中安寧河流域水域和耕地ESV增長11.52%和3.16%、草地和建設(shè)用地ESV變化率為-3.58%和85.73%，白龍江流域水域和林地ESV增長12.46%和0.22%、建設(shè)用地和耕地變化率達120.59%和-4.63%?？梢姲矊幒恿饔蛄植莸孛娣e的銳減、白龍江流域耕地和草地的縮減以及兩地城鎮(zhèn)建設(shè)的快速推進，使區(qū)域生態(tài)環(huán)境逐步惡化。

通過1 km×1 km和2 km×2 km格網(wǎng)分別計算安寧河流域與白龍江流域單位面積ESV（萬元/hm2），并利用ArcGIS 10.7進行克里金插值，依據(jù)自然斷點法進行分級。其中，安寧河流域分為極低價值［-1.033 0，1.792 1）、低價值［1.792 1，2.929 1）、中等價值［2.929 1，3.740 2）、高價值［3.740 2，11.822 5）和極高價值［11.822 5，25.510 0）；白龍江流域分為極低價值［1.148 4，2.403 5）、低價值［2.403 5，3.194 0）、中等價值［3.194 0，3.845 0）、高價值［3.845 0，6.986 7）和極高價值［6.986 7，9.175 8）。得到2000～2020年兩地單位面積ESV的空間分布圖（見圖3），同時對不同等級ESV的面積占比進行統(tǒng)計（見圖4）。近20 a，安寧河流域與白龍江流域單位面積ESV以高和中等兩個等級為主，兩等級面積占比之和分別在69.6%和71.2%以上，此外兩地低等級單位面積ESV占比分別在22.2%和27.5% 以上，除安寧河流域低、中等和極高3個等級及白龍江流域中等等級單位面積ESV占比有所降低外，其余等級緩慢增長。高ESV集中于兩流域水域及林地密集帶，低ESV主要分布于耕地和建設(shè)用地集中的河谷地帶，中等ESV則集中于兩地草地及林草過渡帶。

近20 a，安寧河流域和白龍江流域單位面積ESV平均值分別下降0.6%和0.1%，表明經(jīng)濟快速發(fā)展、城鎮(zhèn)化加快推進等人為因素干擾下，兩地生態(tài)環(huán)境在不斷惡化，威脅區(qū)域生態(tài)可持續(xù)發(fā)展。

2.3 LERI動態(tài)變化分析

同樣借助ArcGIS 10.7自然斷點法將格網(wǎng)尺度下兩流域LERI進行分級（見圖5），并統(tǒng)計各等級面積占比變化（見圖6）。其中，安寧河流域分為極低風(fēng)險［0.016 1，0.026 4）、低風(fēng)險［0.026 4，0.041 2）、中等風(fēng)險［0.041 2，0.063 4）、高風(fēng)險［0.063 4，0.153 6）和極高風(fēng)險［0.153 6，0.332 1）；白龍江流域分為極低風(fēng)險［0.013 8，0.026 4）、低風(fēng)險［0.026 4，0.041 7）、中等風(fēng)險［0.041 7，0.061 9）、高風(fēng)險［0.061 9，0.089 1）和極高風(fēng)險［0.089 1，0.191 7）。2000～2020年安寧河流域和白龍江流域LERI以極低和低兩個等級為主，兩等級面積占比之和分別達66.01%和75.00%以上，兩地中等LERI占比分別在17.15%和12.38%以上，除安寧河流域極低、低兩等級及白龍江流域極低等級LERI占比有所增長外，其余等級緩慢增長。高LERI集中在兩地耕地、建設(shè)用地集中的河谷地帶和水域面積較大區(qū)，低LERI主要分布于兩地林地密集區(qū)域，中等LERI則處于河谷人類活動頻繁區(qū)與林草等植被過渡地帶。

近20 a，安寧河流域和白龍江流域LERI平均值分別下降8.24%和3.50%，表明隨著人類活動對流域生態(tài)環(huán)境的干擾，區(qū)域土地利用趨于集中，致使破碎度降低，兩地景觀生態(tài)風(fēng)險得到一定緩解。

2.4 ESV與LERI的空間關(guān)聯(lián)分析

運用GeoDa計算2000～2020年研究區(qū)ESV和LERI的雙變量空間自相關(guān)指數(shù)（見表4）。兩流域Morans I均小于0且逐步降低，同時|Z|均大于2.58，P＜0.001表明ESV和LERI存在顯著空間負相關(guān)性并呈增強趨勢。

進一步借助GeoDa對2000～2020年研究區(qū)ESV和LERI進行雙變量局部空間自相關(guān)分析（見圖7），并統(tǒng)計各自相關(guān)等級的面積占比及變化（見圖8）。安寧河流域和白龍江流域高價值-高風(fēng)險區(qū)分別占 1.77% 和1.91%以上，前者呈增長趨勢，主要分布于流域中上游邛海、大橋水庫（1999年蓄水）、冶勒水庫（2001年動工）周邊及支干流河面寬闊處；后者呈降低趨勢，集中于流域中下游碧口（1997年建成）、寶珠寺（1998年建成）和苗家壩（2014年建成）3個大型水電站庫區(qū)及支干流河面寬闊處，這些區(qū)域涵蓋面積較大且能提供高價值ESV的水域，但被人類活動侵占的風(fēng)險極高，故需加強對水域的保護。安寧河流域和白龍江流域高價值-低風(fēng)險區(qū)占27.47%和27.96%以上，前者呈擴張趨勢，集中于北部小相嶺山脈、西部牦牛山脈和東南部螺髻-魯南山脈；后者呈減少趨勢，主要分布于岷山-龍門山脈，上述區(qū)域以能創(chuàng)造較高ESV的林地為主，因面積極大，受經(jīng)濟活動影響有限，LERI較低。低價值-高風(fēng)險區(qū)主要分布于耕地和建設(shè)用地集中的寬闊河谷帶，ESV較低的同時極易受人類活動干擾，故LERI較高；其中安寧河冕寧縣-米易縣段、白龍江宕昌縣-文縣和劍閣縣段有大面積分布，兩地占比在12.93%和12.97%以上，均呈增長態(tài)勢。低價值-低風(fēng)險區(qū)安寧河流域和白龍江流域分別占1.24%和1.17%以上，均緩慢增長，零散分布在兩流域高山林草過渡帶，ESV較低且較少受人類活動影響，LERI較低。不顯著區(qū)主要分布于安寧河流域的東部、南部和白龍江流域的西北部、中部以及兩地河谷邊緣帶，多以草地、耕地向林地過渡帶為主，故ESV和LERI適中，兩地占55.73%和55.77%以上，前者有所降低，后者緩慢增長。

2.5 ESV與LERI空間分異的驅(qū)動分析

要素空間分異受自然條件和經(jīng)濟條件的復(fù)合作用影響，通過地理探測器對選取的高程（X1）、坡度（X2）、年均溫（X3）、年降水（X4）、歸一化植被指數(shù)（X5）、人口密度（X6）、地均GDP（X7）、HAI（X8）、ESV/LERI（X9）驅(qū)動因子進行探測。從ESV空間分異探測結(jié)果看（見圖9），安寧河流域ESV空間分異驅(qū)動因子按q統(tǒng)計量大小進行整體排序為X8＞X1＞X2＞X6＞X3＞X4＞X7＞X5，白龍江流域為X8＞X5＞X4＞X7＞X1＞X6＞X3＞X2?？梢奨8是兩地ESV分異的主導(dǎo)因子，解釋力分別在58.31%和42.29%以上。此外，X1、X2、X6是安寧河流域ESV分異的重要因素，解釋力在21.94%以上；X5和X4對白龍江流域ESV分異的影響較大，解釋力在33.31%和20.45%以上；其余因子對兩地ESV空間分異的影響力不顯著。從時間變化看，2000～2020年除安寧河流域X1、X5和X7以及白龍江流域X1、X4和X6的影響力有一定增強外，其余因子的影響力不同程度減弱。

從2000～2020年LERI空間分異探測結(jié)果看（見圖10），安寧河流域按q統(tǒng)計量由高到低排序為X8＞X1＞X2＞X6＞X3＞X4＞X7＞X5，白龍江流域為X8＞X5＞X4＞X6＞X7＞X3＞X1＞X2?？芍猉8是兩地LERI空間分異的主導(dǎo)因素，解釋力分別達67.80%和67.27% 以上。X1、X2、X3和X6對安寧河流域LERI分異產(chǎn)生較大影響，解釋力達25.09%以上；X4、X5和X6對白龍江流域LERI分異影響較大，解釋力分別在7.52%，31.75%和9.33%以上；其他因子對兩地LERI 空間分異的影響較小。在時間上，除安寧河流域X1、X5 、X7和X8以及白龍江流域X6的影響力有所增強外，其他因子影響力呈降低趨勢。

研究表明，自然因素和社會經(jīng)濟因素共同對兩個典型干旱河谷ESV和LERI空間分異產(chǎn)生影響。X8均為兩地ESV和LERI分異的主導(dǎo)因子，說明人類活動對兩流域生態(tài)環(huán)境的干擾程度極大，需規(guī)范土地利用開發(fā)，加強對人類干擾所帶來的生態(tài)風(fēng)險管控。此外，安寧河流域X1、X2和白龍江流域X4和X5等自然因子較經(jīng)濟因子更易對流域ESV和LERI產(chǎn)生影響，故需加強流域生態(tài)環(huán)境保護，遵循因地制宜的資源開發(fā)模式，促進人地關(guān)系和諧健康發(fā)展。

3 討 論

本文借鑒近年來ESV和LERI相關(guān)研究［1，30］，基于土地利用數(shù)據(jù)探討了兩個典型干旱河谷ESV和LERI的時空演變、二者關(guān)聯(lián)性及空間分異性特征。結(jié)果顯示，2000～2020年兩個典型干旱河谷ESV和LERI呈持續(xù)縮減趨勢，與多數(shù)干旱河谷ESV預(yù)測結(jié)果相悖［31］。主要與以下因素有關(guān)：21世紀以來，因作為西部邊遠山區(qū)的區(qū)域經(jīng)濟中心，兩地在西部大開發(fā)、退耕還林還草及城鎮(zhèn)化建設(shè)加快和人口激增的影響下，土地利用結(jié)構(gòu)發(fā)生劇烈變化。建設(shè)用地的快速擴張，使安寧河流域草地、白龍江流域草地和耕地明顯減少，雖安寧河流域耕地、白龍江流域林地及兩地水域的增長對ESV起一定促進作用，但兩地快速城鎮(zhèn)化所產(chǎn)生的極大阻礙作用使ESV不斷下降，而當(dāng)前建設(shè)用地對ESV的巨大消極影響往往被忽視［12］。20 a 來，兩地LERI緩慢降低，景觀生態(tài)風(fēng)險狀況呈轉(zhuǎn)好趨勢，與其他經(jīng)濟較發(fā)達區(qū)結(jié)果一致［4，16］，在兩地主要地類斑塊數(shù)不斷減少影響下，分離度和破碎度降低，但人類活動頻繁區(qū)LERI極高，需格外注意生態(tài)風(fēng)險防范。因此，需根據(jù)ESV和LERI主要集聚關(guān)系，對干旱河谷不同分區(qū)因地制宜地采取針對性措施，以促進區(qū)域協(xié)調(diào)發(fā)展：河流支干流河谷地帶為低價值-高風(fēng)險區(qū)，人口密集、經(jīng)濟建設(shè)活動活躍，要格外注重生態(tài)環(huán)境保護，減少對耕地的侵占；水庫、湖泊及大型水電站庫區(qū)為高價值-高風(fēng)險區(qū)，對局部生態(tài)至關(guān)重要，需加強水域保護力度；兩地占主導(dǎo)的林地和草地為高價值-低風(fēng)險區(qū)，是流域生態(tài)安全及可持續(xù)發(fā)展的基石，需加大對林草生態(tài)系統(tǒng)的保護。

將ESV和LERI二者進行結(jié)合更有助于區(qū)域生態(tài)環(huán)境保護決策的制定與實施［16］。但本研究仍有所不足：ESV和LERI均通過土地利用進行評估，無法全面準(zhǔn)確反映區(qū)域?qū)嶋H；通過研究區(qū)糧食產(chǎn)量及價格來修正ESV的方法有待完善；ESV和LERI空間分異的驅(qū)動因子選取難以深刻反映區(qū)域內(nèi)部差異。因此，今后需探索更為科學(xué)的ESV和LERI評估方法，從自然和經(jīng)濟等多角度深度分析二者的空間分異性。

4 結(jié) 論

2000～2020年安寧河流域和白龍江流域ESV均呈減少趨勢，前者由370.597 3億元減少到368.337 2億元、后者由1 137.410 1億元減少到1 136.584 4億元，林地和草地ESV貢獻率分別達91.58%和95.19%以上，ESV以中等和高兩個等級為主，建設(shè)用地的快速擴張是兩地ESV縮減的主要原因；兩地LERI呈降低趨勢，分別下降8.24%和3.50%，以極低和低兩個等級為主，LERI降低主要緣于土地利用集約化加強引發(fā)的破碎度下降；兩地ESV和LERI具有顯著空間負相關(guān)性且呈增強態(tài)勢，高價值-低風(fēng)險和低價值-高風(fēng)險為兩者主要關(guān)系；人類活動對兩地ESV和LERI空間分異的干擾極大，其中HAI是兩地ESV和LERI空間分異的主導(dǎo)因子。

參考文獻：

［1］李輝，周啟剛，李斌，等.近30年三峽庫區(qū)生態(tài)系統(tǒng)服務(wù)價值與生態(tài)風(fēng)險時空變化及相關(guān)性研究［J］.長江流域資源與環(huán)境，2021，30（3）：654-666.

［2］周忠學(xué)，仇立慧.城市化對生態(tài)系統(tǒng)服務(wù)功能影響的實證研究：以西安市南郊為例［J］.干旱區(qū)研究，2011，28（6）：974-979.

［3］李景剛，何春陽，李曉兵.快速城市化地區(qū)自然/半自然景觀空間生態(tài)風(fēng)險評價研究：以北京為例［J］.自然資源學(xué)報，2008，23（1）：33-47.

［4］朱潤苗，陳松林.1980～2020年福建省景觀生態(tài)風(fēng)險與生態(tài)系統(tǒng)服務(wù)價值的空間關(guān)系［J］.應(yīng)用生態(tài)學(xué)報，2022，33（6）：1599-1607.

［5］COSTANZA R，DARGE R，GROOT R，et al.The value of the world′s ecosystem services and natural capital［J］.Ecological Economics，1998，25（1）：3-15.

［6］COSTANZA R，GROOT R，BRAAT L，et al.Twenty years of ecosystem services：how far have we come and how far do we still need to go？［J］.Ecosystem Services，2017，28：1-16.

［7］虞依娜，彭少麟.生態(tài)系統(tǒng)服務(wù)價值評估的研究進展［J］.生態(tài)環(huán)境學(xué)報，2010，19（9）：2246-2252.

［8］歐陽志云，王如松，趙景柱.生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟價值評價［J］.應(yīng)用生態(tài)學(xué)報，1999，10（5）：635-640.

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

［10］ 于德永，郝蕊芳.生態(tài)系統(tǒng)服務(wù)研究進展與展望［J］.地球科學(xué)進展，2020，35（8）：804-815.

［11］張曉瑤，陸林，虞虎，等.青藏高原土地利用變化對生態(tài)系統(tǒng)服務(wù)價值影響的多情景模擬［J］.生態(tài)學(xué)雜志，2021，40（3）：887-898.

［12］朱相君，薛亮.關(guān)中地區(qū)生態(tài)系統(tǒng)服務(wù)價值變化及其與經(jīng)濟發(fā)展耦合關(guān)系［J］.生態(tài)學(xué)雜志，2022，41（4）：769-776.

［13］彭建，黨威雄，劉焱序，等.景觀生態(tài)風(fēng)險評價研究進展與展望［J］.地理學(xué)報，2015，70（4）：664-677.

［14］曹祺文，張曦文，馬洪坤，等.景觀生態(tài)風(fēng)險研究進展及基于生態(tài)系統(tǒng)服務(wù)的評價框架：ESRISK［J］.地理學(xué)報，2018，73（5）：843-855.

［15］賈艷艷，唐曉嵐.長江流域安徽段生態(tài)系統(tǒng)服務(wù)價值與景觀生態(tài)風(fēng)險時空演變及其關(guān)聯(lián)分析［J］.南京林業(yè)大學(xué)學(xué)報（自然科學(xué)版），2022，46（3）：31-40.

［16］李俊翰，高明秀.濱州市生態(tài)系統(tǒng)服務(wù)價值與生態(tài)風(fēng)險時空演變及其關(guān)聯(lián)性［J］.生態(tài)學(xué)報，2019，39（21）：7815-7828.

［17］殷楠，王帥，劉焱序.生態(tài)系統(tǒng)服務(wù)價值評估：研究進展與展望［J］.生態(tài)學(xué)雜志，2021，40（1）：233-244.

［18］廖成浩，姚昆，曾艷.安寧河流域土地生態(tài)安全動態(tài)評價與預(yù)測［J］.環(huán)境監(jiān)測管理與技術(shù)，2022，34（2）：27-31，37.

［19］范建容，楊超，包維楷，等.西南地區(qū)干旱河谷分布范圍及分區(qū)統(tǒng)計分析［J］.山地學(xué)報，2020，38（2）：303-313.

［20］ZHANG X，LIU L Y，CHEN X D，et al.GLC_FCS30：global land-cover product with fine classification system at 30 m using time-sLERIes Landsat imagery［J］.Earth System Science Data，2021，13（6）：2753-2776.

［21］謝高地，張彩霞，張昌順，等.中國生態(tài)系統(tǒng)服務(wù)的價值［J］.資源科學(xué)，2015，37（9）：1740-1746.

［22］LIU Y C，LIU Y X，LI J L，et al.Evolution of landscape ecological risk at the optimal scale：a case study of the open coastal wetlands in Jiangsu，China［J］.International Journal of Environmental Research and Public Health，2018，15（8）：1691-1701.

［23］韓振華，李建東，殷紅，等.基于景觀格局的遼河三角洲濕地生態(tài)安全分析［J］.生態(tài)環(huán)境學(xué)報，2010，19（3）：701-705.

［24］胡和兵，劉紅玉，郝敬鋒，等.流域景觀結(jié)構(gòu)的城市化影響與生態(tài)風(fēng)險評價［J］.生態(tài)學(xué)報，2011，31（12）：3432-3440.

［25］高賓，李小玉，李志剛，等.基于景觀格局的錦州灣沿海經(jīng)濟開發(fā)區(qū)生態(tài)風(fēng)險分析［J］.生態(tài)學(xué)報，2011，31（12）：3441-3450.

［26］陳浮，葛小平，陳剛，等.城市邊緣區(qū)景觀變化與人為影響的空間分異研究［J］.地理科學(xué)，2001，21（3）：210-216.

［27］黃木易，方斌，岳文澤，等.近20 a來巢湖流域生態(tài)服務(wù)價值空間分異機制的地理探測［J］.地理研究，2019，38（11）：2790-2803.

［28］楊遠琴，任平，洪步庭.基于生態(tài)安全的三峽庫區(qū)重慶段土地利用沖突識別［J］.長江流域資源與環(huán)境，2019，28（2）：322-332.

［29］王勁峰，徐成東.地理探測器：原理與展望［J］.地理學(xué)報，2017，72（1）：116-134.

［30］JIA Y Y，TANG X L，LIU W.Spatial-temporal evolution and correlation analysis of ecosystem service value and landscape ecological risk in Wuhu City［J］.Sustainability，2020，12（7）：2803-2819.

［31］歐朝蓉，朱清科，孫永玉.干熱河谷區(qū)土地利用與生態(tài)系統(tǒng)服務(wù)價值：以云南省元謀縣為例［J］.農(nóng)業(yè)現(xiàn)代化研究，2018，39（3）：494-502.

（編輯：黃文晉）

Spatial-temporal evolution of ecosystem service value and landscape ecological risk in dry valleys

ZHOU Jianwei1，2，3，LUO Jun2，3

（1.College of Life Sciences，China West Normal University，Nanchong 637009，China； 2.Sichuan Provincial Engineering Laboratory of Monitoring and Control for Soil Erosion on Dry Valleys，Nanchong 637009，China； 3.Liangshan Soil Erosion and Ecological Restoration in Dry Valleys Observation and Research Station，Liangshan 616753，China）

Abstract：

Based on the high-precision land use data of Anning River Basin and Bailong River Basin from 2000 to 2020，the ecosystem service value （ESV） and landscape ecological riskindex （LERI） assessment model，bivariate spatial autocorrelation model and Geodetector were used to explore the spatio-temporal evolution，correlation and differentiation characteristics of ESV and LERI in different dry valleys.The results indicated that forest and grass dominated the Anning River Basin and Bailong River Basin.With the influence of urban expansion，ESV and LERI showed a decreasing trend.ESV was mainly in low and medium grade，while LERI was mainly in extremely low and low grade.ESV and LERI had significant negative spatial correlation，and the relationship between ESV and LERI were dominated by high value-low risk and low value-high risk.Human activities played a leading role in spatial differentiation of ESV and LERI.Forestland and grassland were important sources of ESV in the two basins，and the surge of construction land leaded to the continuous decreasing of ESV and LERI in the two basins.It is necessary to improve land use structure according to local conditions，strengthen protection of forest，grass and water area，control urban expansion reasonably，and promote sustainable development of land use structure in watershed.

Key words： ecosystem service value；landscape ecological risk index；dry valley；spatial-temporal evolution；spatio-temporal differentiation

收稿日期：2022-06-24

基金項目：國家自然科學(xué)基金項目（41971015）；生態(tài)安全與保護四川省重點實驗室開放基金項目（ESP201301）

作者簡介：周建偉，男，碩士研究生，主要從事水土保持與荒漠化防治研究。E-mail：gesang3650@163.com

通信作者：羅 君，男，講師，博士，主要從事地理綜合研究。E-mail：luojunxmx@126.com