2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化評(píng)估

摘要：近年來，廣西壯族自治區(qū)（簡(jiǎn)稱廣西）社會(huì)經(jīng)濟(jì)穩(wěn)步發(fā)展，但部分因素依然制約區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值的持續(xù)提升。分析2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的時(shí)空變化特征及二者關(guān)聯(lián)性，可為廣西生態(tài)系統(tǒng)保護(hù)和修復(fù)提供參考。本文運(yùn)用單位面積當(dāng)量因子法、生態(tài)風(fēng)險(xiǎn)核算模型、雙變量空間自相關(guān)分析，核算2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)指數(shù)，并分析其時(shí)空分異特征及其相關(guān)性。結(jié)果表明，2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值增減起伏，從2000年的25 240.24億元增加到2020年的25 251.56億元，總體增幅為0.05%。生態(tài)風(fēng)險(xiǎn)等級(jí)總體呈現(xiàn)南高北低的空間格局，各區(qū)域面積表現(xiàn)出增減起伏的走向。從數(shù)量相關(guān)性來看，生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)具有較強(qiáng)的相關(guān)性；從空間相關(guān)性來看，生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)在空間上呈聚集分布，呈顯著的負(fù)相關(guān)。廣西生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)呈負(fù)相關(guān)，二者的結(jié)合研究可以更好地體現(xiàn)人類活動(dòng)造成的生態(tài)環(huán)境變化，為生態(tài)系統(tǒng)修復(fù)及土地可持續(xù)利用提供支持。

關(guān)鍵詞：生態(tài)系統(tǒng)服務(wù)價(jià)值；生態(tài)風(fēng)險(xiǎn)；時(shí)空變化；評(píng)估；廣西

中圖分類號(hào)：X171.1 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1008-9500（2024）03-0-08

DOI：10.3969/j.issn.1008-9500.2024.03.032

Abstract： In recent years， the social economy of Guangxi Zhuang Autonomous Region （referred to as Guangxi） has steadily developed， but some factors still constrain the continuous improvement of regional ecosystem service value. Analyzing the spatiotemporal changes of ecosystem service value and ecological risk in Guangxi from 2000 to 2020， as well as their correlation， can provide reference for the protection and restoration of ecosystems in Guangxi. This paper uses the unit area equivalent factor method， ecological risk accounting model， and bivariate spatial autocorrelation analysis to calculate the ecosystem service value and ecological risk index in Guangxi from 2000 to 2020， and analyzes their spatiotemporal differentiation characteristics and correlation. The results show that the ecosystem service value in Guangxi fluctuated from 2000 to 2020， increasing from 2 524.024 billion yuan in 2000 to 2 525.156 billion yuan in 2020， with an overall increase of 0.05%. The overall ecological risk level presents a spatial pattern of high in the south and low in the north， and the area of each region shows a fluctuating trend of increase and decrease. From the perspective of quantitative correlation， there is a strong correlation between ecosystem service value and ecological risk; from the perspective of spatial correlation， the value of ecosystem services and ecological risks exhibit a clustered distribution in space， showing a significant negative correlation. The ecosystem service value in Guangxi is negatively correlated with ecological risks， and the combination study of the two can better reflect the ecological environment changes caused by human activities， providing support for ecosystem restoration and sustainable land use.

Keywords： ecosystem service value; ecological risks; spatiotemporal changes; assessment; Guangxi

近幾十年來，由于生產(chǎn)和生活的需要，人類對(duì)土地等自然資源的開發(fā)利用強(qiáng)度不斷變大，不斷改變自然環(huán)境，各種生態(tài)系統(tǒng)受到破壞，面臨退化。隨著人類活動(dòng)范圍的不斷擴(kuò)大，氣候變化異常頻繁出現(xiàn)，全球環(huán)境受到威脅，生態(tài)和氣候成為世界關(guān)注的焦點(diǎn)[1]。因此，有必要研究土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響，分析生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的相關(guān)性，為實(shí)現(xiàn)可持續(xù)發(fā)展提供科學(xué)依據(jù)。目前，國內(nèi)外學(xué)者在生態(tài)系統(tǒng)服務(wù)價(jià)值研究方面取得一定進(jìn)展。研究內(nèi)容包括單個(gè)生態(tài)系統(tǒng)服務(wù)價(jià)值（森林等生態(tài)系統(tǒng)）[2-3]、草原與海洋等生態(tài)系統(tǒng)的功能價(jià)值（涵養(yǎng)水源、固碳釋氧等）[4-5]；研究區(qū)域大多集中在典型山區(qū)和高原地區(qū)[6-7]、邊境地區(qū)和城市群[8-9]；研究尺度涵蓋全球[10]、國家[11-12]和地區(qū)[13-14]等。當(dāng)前研究注重把生態(tài)系統(tǒng)服務(wù)價(jià)值加入生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)，或者將生態(tài)系統(tǒng)服務(wù)價(jià)值與景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)合。生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的評(píng)價(jià)逐漸趨于結(jié)合，當(dāng)前對(duì)兩者時(shí)空關(guān)系的研究還不多。

廣西壯族自治區(qū)（簡(jiǎn)稱廣西）是我國唯一一個(gè)沿邊、沿江、沿海的少數(shù)民族自治區(qū)[15]。新時(shí)代，廣西是“一帶一路”有機(jī)銜接的重要門戶，是面向東盟的重要國際通道、西部陸海新通道和新的自由貿(mào)易區(qū)[16]。近幾十年來，廣西社會(huì)經(jīng)濟(jì)蓬勃發(fā)展，但人口增長、土地資源粗放利用、環(huán)境污染和生態(tài)破壞等問題制約生態(tài)系統(tǒng)服務(wù)價(jià)值的持續(xù)提升。因此，本文基于土地利用數(shù)據(jù)和網(wǎng)格地理信息系統(tǒng)（Geographic Information System，GIS），采用5 km×5 km網(wǎng)格作為統(tǒng)計(jì)單元，運(yùn)用單位面積當(dāng)量因子法、生態(tài)風(fēng)險(xiǎn)核算模型、雙變量空間自相關(guān)分析，核算2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)指數(shù)，并分析其時(shí)空分異特征，探究兩者的關(guān)聯(lián)性，為今后廣西生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)的研究以及生態(tài)系統(tǒng)保護(hù)提供支持。

1 研究區(qū)概況

廣西地處華南地區(qū)，地理坐標(biāo)為104°28′E～112°04′E、20°54′N～26°24′N，它是我國西部唯一的沿海地區(qū)，背靠西南地區(qū)，毗鄰廣東省，并連接?xùn)|南亞。廣西在海上絲綢之路中處于極其重要的位置，同時(shí)在國家對(duì)外開放和西部大開發(fā)中具有特殊地位。廣西面積為23.76萬km2，占全國面積的2.6%。廣西礦產(chǎn)資源十分豐富，尤其是鋁、錫等有色金屬。廣西河流眾多，水資源充足，總量達(dá)到2 377.5億m3左右。境內(nèi)河流分為四大水系，即長江、珠江、百都河、桂南獨(dú)流入海。廣西有豐富的海洋資源，海岸線長約為1 788 km，大陸海岸線曲折，有很多溺谷且面積寬廣，可以開發(fā)的沿海港口就有20余個(gè)。廣西灘涂面積為10萬hm2左右，紅樹林面積居全國首位。廣西位于低緯度地區(qū)，中部穿過北回歸線，南部與北部灣相鄰，西部和云貴高原相接，北部與南嶺接壤，它是典型的亞熱帶季風(fēng)氣候區(qū)。廣西氣候暖和，多雨水、光照，全年平均氣溫介于17.6～23.6 ℃。

2 數(shù)據(jù)來源與研究方法

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

2000年、2005年、2010年、2015年和2020年廣西土地利用數(shù)據(jù)來自中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心，空間分辨率為30 m。廣西行政區(qū)劃矢量數(shù)據(jù)來源于國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心。數(shù)字高程模型（Digital Elevation Model，DEM）數(shù)據(jù)來自地理空間數(shù)據(jù)云，空間分辨率為30 m。其他基礎(chǔ)統(tǒng)計(jì)數(shù)據(jù)來自歷年（2000—2021年）發(fā)布的《中國統(tǒng)計(jì)年鑒》與《廣西統(tǒng)計(jì)年鑒》，相關(guān)數(shù)據(jù)根據(jù)基礎(chǔ)數(shù)據(jù)計(jì)算得到。

2.2 研究方法

2.2.1 生態(tài)系統(tǒng)服務(wù)價(jià)值核算模型

基于謝高地等[17]修改的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表，考慮廣西實(shí)際情況和相關(guān)修正規(guī)則，得到研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表。通過糧食產(chǎn)量、種植面積和平均收購價(jià)格，計(jì)算地均生態(tài)系統(tǒng)服務(wù)價(jià)值和修正后糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值（生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子）。將糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值和生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子相乘，能夠得出各土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，如表1所示。采用ArcGIS軟件計(jì)算每個(gè)網(wǎng)格的各土地利用類型面積，然后采用式（1）計(jì)算各土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值，再采用式（2）計(jì)算每個(gè)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值。

式中：ESVi為第i種土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值；ESV為每個(gè)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值；Ai為第i種土地利用類型的面積；VCi為第i種土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)。

2.2.2 土地利用生態(tài)風(fēng)險(xiǎn)指數(shù)核算模型

土地利用類型的空間組合及比例不同，土地利用生態(tài)風(fēng)險(xiǎn)也存在差異，因此按照不同土地利用類型的面積比例，采用式（3）計(jì)算廣西土地利用的生態(tài)風(fēng)險(xiǎn)指數(shù)。其間采用反正切函數(shù)對(duì)相關(guān)數(shù)據(jù)進(jìn)行歸一化處理，并結(jié)合相關(guān)研究[18-20]和研究區(qū)特點(diǎn)，得到研究區(qū)各土地利用類型的生態(tài)風(fēng)險(xiǎn)系數(shù)Ci。其中，耕地取0.326 5，林地取0.145 5，草地取0.248 1，水體取0.282 1，建設(shè)用地取0.398 5，未利用地取0.318 2。

式中：ERIk為第k個(gè)網(wǎng)格內(nèi)土地利用生態(tài)風(fēng)險(xiǎn)指數(shù)；Ski為第k個(gè)網(wǎng)格內(nèi)土地利用類型i的面積；Si為研究區(qū)第i種土地利用類型的面積；Ci為第i種土地利用類型的生態(tài)風(fēng)險(xiǎn)系數(shù)。

2.2.3 雙變量空間自相關(guān)分析

雙變量空間自相關(guān)是單變量莫蘭指數(shù)的延伸和拓展，全局莫蘭指數(shù)可反映生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)的空間相關(guān)性和差異[21]，采用式（4）進(jìn)行計(jì)算。

式中：Iab為單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值a和生態(tài)風(fēng)險(xiǎn)指數(shù)b的全局莫蘭指數(shù)；n為評(píng)價(jià)單元的個(gè)數(shù)；Wij為評(píng)價(jià)單元i與j的空間權(quán)重；yia和yib分別為第i個(gè)評(píng)價(jià)單元單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)指數(shù)；ya與yb分別為單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)指數(shù)平均值；σa和σb分別為生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)指數(shù)的方差。

3 結(jié)果與分析

3.1 生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化分析

運(yùn)用ArcGIS軟件，并結(jié)合廣西生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，可以得出2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值及其變化，如表2、表3所示。數(shù)據(jù)顯示，2000—2020年，廣西生態(tài)系統(tǒng)服務(wù)價(jià)值由2000年的25 240.24億元增長到2020年的25 251.56億元，增長率為0.05%，其表現(xiàn)出先增加、后減少、再增加的變化趨勢(shì)。單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值由106 229.97元/hm2增長到106 277.61元/hm2。6種土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值占總價(jià)值的比重排序?yàn)椋毫值兀静莸兀舅颍靖兀疚蠢玫兀窘ㄔO(shè)用地。總體來看，耕地生態(tài)系統(tǒng)服務(wù)價(jià)值不斷減少，林地生態(tài)系統(tǒng)服務(wù)價(jià)值先增后減，草地生態(tài)系統(tǒng)服務(wù)價(jià)值增減起伏，水域生態(tài)系統(tǒng)服務(wù)價(jià)值持續(xù)增長，建設(shè)用地和未利用地生態(tài)系統(tǒng)服務(wù)價(jià)值變化不明顯。

運(yùn)用ArcGIS軟件的漁網(wǎng)工具，借助面積制表、字段計(jì)算器等工具，計(jì)算每個(gè)評(píng)價(jià)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值。經(jīng)分析，評(píng)價(jià)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值呈正態(tài)分布，因此使用普通克里金法對(duì)其進(jìn)行插值。為了清晰地展現(xiàn)生態(tài)系統(tǒng)服務(wù)價(jià)值在空間上的異質(zhì)性，使用自然斷裂法將其劃分為低、較低、中、較高、高五個(gè)等級(jí)。生態(tài)系統(tǒng)服務(wù)價(jià)值等級(jí)為高和較高的占比最大，其次為較低和中。總體上看，2000—2020年廣西生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布變化不顯著，呈相對(duì)穩(wěn)定的狀態(tài)。生態(tài)系統(tǒng)服務(wù)價(jià)值等級(jí)為高和較高的區(qū)域集中分布在廣西的北部、東部和西南部，部分分布在東南部，即分布在河池市、百色市、柳州市、桂林市、梧州市、賀州市、玉林市、防城港市等城市。廣西北部、西南部以及東部多山地和丘陵，如十萬大山、九萬大山、金鐘山等，這些區(qū)域森林資源較為豐富且林地、草地比較多，生態(tài)環(huán)境良好，因此生態(tài)系統(tǒng)服務(wù)價(jià)值較高。生態(tài)系統(tǒng)服務(wù)價(jià)值等級(jí)為中的區(qū)域主要分布在桂林市東北部、柳州市中部和百色市東南部，南寧市、崇左市、欽州市和賀州市等城市的部分地區(qū)也有分布。生態(tài)系統(tǒng)服務(wù)價(jià)值等級(jí)為較低的區(qū)域則分布在南寧市、桂林市、柳州市、來賓市、賀州市和欽州市等城市的部分區(qū)域。生態(tài)系統(tǒng)服務(wù)價(jià)值等級(jí)為低的區(qū)域多分布在南寧市、貴港市、來賓市和北海市，主要是南寧市及其周邊地區(qū)，廣西沿海地區(qū)也分布較多。相比其他生態(tài)系統(tǒng)服務(wù)價(jià)值高的區(qū)域，這些城市的山地和林地面積較少，人類對(duì)自然資源的開發(fā)強(qiáng)度較大，各種建設(shè)項(xiàng)目占用林地、草地，林地、草地轉(zhuǎn)為耕地，使得生態(tài)系統(tǒng)服務(wù)價(jià)值較低。

3.2 生態(tài)風(fēng)險(xiǎn)時(shí)空變化分析

借助ArcGIS軟件的漁網(wǎng)工具，使用5 km×5 km的網(wǎng)格劃分評(píng)價(jià)單元，得到9 940個(gè)評(píng)價(jià)單元，在格網(wǎng)尺度下計(jì)算土地利用生態(tài)風(fēng)險(xiǎn)指數(shù)。經(jīng)分析，生態(tài)風(fēng)險(xiǎn)指數(shù)服從正態(tài)分布，因此采用普通克里金法對(duì)其進(jìn)行空間插值。從插值結(jié)果看，2000—2020年，土地利用生態(tài)風(fēng)險(xiǎn)指數(shù)均介于0.024～0.357。為了更好地反映土地利用生態(tài)風(fēng)險(xiǎn)指數(shù)的空間分布，使用自然斷裂法將其劃分為低風(fēng)險(xiǎn)、較低風(fēng)險(xiǎn)、中風(fēng)險(xiǎn)、較高風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)五個(gè)風(fēng)險(xiǎn)等級(jí)。

2000—2020年，生態(tài)風(fēng)險(xiǎn)等級(jí)為高風(fēng)險(xiǎn)和較高風(fēng)險(xiǎn)的區(qū)域主要集中在南寧市、貴港市、玉林市、欽州市、北海市和來賓市，而河池市、百色市、柳州市、桂林市、賀州市和梧州市則分布相對(duì)較少。經(jīng)分析，南寧市、貴港市和玉林市等城市經(jīng)濟(jì)不斷發(fā)展，大量的房地產(chǎn)、商業(yè)、工業(yè)等項(xiàng)目需要建設(shè)用地，而建設(shè)用地存量不足，因此其占用林地、草地等生態(tài)風(fēng)險(xiǎn)較低的土地利用類型，使這些區(qū)域的生態(tài)風(fēng)險(xiǎn)處在較高水平。另外，因?yàn)榻?jīng)濟(jì)發(fā)展和糧食安全的需要，這些地區(qū)的林地、草地轉(zhuǎn)為耕地，也會(huì)使得生態(tài)風(fēng)險(xiǎn)增大。生態(tài)風(fēng)險(xiǎn)等級(jí)為中風(fēng)險(xiǎn)的區(qū)域主要集中在廣西的北部和東部，少量分布在東南部和西南部。生態(tài)風(fēng)險(xiǎn)等級(jí)為較低風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)的區(qū)域主要分布在河池市、百色市、柳州市和桂林市，少部分集中分布在崇左市、欽州市、防城港市和北海市的部分區(qū)域。經(jīng)分析，河池市、百色市、柳州市和桂林市都位于廣西北部，多山地和丘陵，有林地、草地等，生態(tài)環(huán)境較好，自然災(zāi)害較少，所以這些地區(qū)的生態(tài)風(fēng)險(xiǎn)比其他地區(qū)低。總的來說，南高北低是廣西生態(tài)風(fēng)險(xiǎn)空間分布的基本特征。經(jīng)分析，廣西北部多山地和丘陵，植被種類較為豐富，大規(guī)模的自然災(zāi)害較少，生態(tài)系統(tǒng)處在比較健康和穩(wěn)定的狀態(tài)，其生態(tài)風(fēng)險(xiǎn)較低；相對(duì)北部而言，廣西南部地勢(shì)比較平緩，經(jīng)濟(jì)發(fā)展開發(fā)區(qū)域集中，建設(shè)用地和耕地較多，生態(tài)風(fēng)險(xiǎn)處在比較高的水平。

由圖1和表4可知，2000—2020年廣西生態(tài)風(fēng)險(xiǎn)空間格局變化不大，但是各區(qū)域的生態(tài)風(fēng)險(xiǎn)空間差異比較明顯，總體呈南高北低的分布格局，以廣西中部為中心向四周分散，生態(tài)風(fēng)險(xiǎn)呈增減起伏的變化趨勢(shì)。其中，生態(tài)風(fēng)險(xiǎn)等級(jí)為低風(fēng)險(xiǎn)的區(qū)域面積先減少后增加；生態(tài)風(fēng)險(xiǎn)等級(jí)為較低風(fēng)險(xiǎn)的區(qū)域面積則持續(xù)增加；生態(tài)風(fēng)險(xiǎn)等級(jí)為中風(fēng)險(xiǎn)的區(qū)域面積則增減起伏，但總體趨于減少；生態(tài)風(fēng)險(xiǎn)等級(jí)為較高風(fēng)險(xiǎn)的區(qū)域面積先減少后增加，再減少，最后增加，總體呈增加趨勢(shì)；生態(tài)風(fēng)險(xiǎn)等級(jí)為高風(fēng)險(xiǎn)的區(qū)域面積先減少后增加，總體呈增加趨勢(shì)。2000—2020年，高風(fēng)險(xiǎn)區(qū)面積增加幅度最大，為1.93%，而中風(fēng)險(xiǎn)區(qū)面積減少幅度最大，達(dá)到2.61%；生態(tài)風(fēng)險(xiǎn)等級(jí)為較低風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)的區(qū)域面積變化相同，增加0.28%，較高風(fēng)險(xiǎn)區(qū)面積增加0.12%，總體而言，它們的變化幅度比高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)區(qū)小。由此可知，廣西生態(tài)風(fēng)險(xiǎn)變化較小，其生態(tài)系統(tǒng)位于比較穩(wěn)定的狀態(tài)。

3.3 相關(guān)性分析

3.3.1 數(shù)量相關(guān)性

為了揭示生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的聯(lián)系，借助SPSS軟件對(duì)各個(gè)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)進(jìn)行正態(tài)性檢驗(yàn)。結(jié)果表明，兩者不服從正態(tài)分布，不適合使用皮爾遜相關(guān)系數(shù)進(jìn)行相關(guān)性分析，而宜采用斯皮爾曼相關(guān)系數(shù)對(duì)相關(guān)性進(jìn)行檢驗(yàn)。研究區(qū)2000年、2005年、2010年、2015年和2020年的斯皮爾曼相關(guān)系數(shù)依次為-0.725 4、-0.734 6、-0.757 1、-0.764 5和-0.785 1，說明生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的相關(guān)性較強(qiáng)。為了探明生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的相關(guān)性，借助SPSS軟件繪制它們的散點(diǎn)圖，如圖2所示。散點(diǎn)圖呈斜三角分布，變化趨勢(shì)不明顯，二者在一定程度上呈負(fù)相關(guān)，即生態(tài)系統(tǒng)服務(wù)價(jià)值越大，生態(tài)風(fēng)險(xiǎn)越小。

3.3.2 空間相關(guān)性

為了進(jìn)一步揭示生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)在空間上是否存在相關(guān)性、是否呈現(xiàn)聚集狀態(tài)，借助GeoDa軟件對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)進(jìn)行雙變量空間自相關(guān)分析。生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的全局空間自相關(guān)分析結(jié)果如表5所示。從全局空間自相關(guān)來看，研究區(qū)2000年、2005年、2010年、2015年和2020年的莫蘭指數(shù)依次為-0.145 2、-0.173 6、-0.194 1、-0.223 4和-0.245 9，莫蘭指數(shù)均小于0，顯著性系數(shù)P均小于0.01，說明生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)在空間分布上呈現(xiàn)聚集狀態(tài)，且呈顯著的負(fù)相關(guān)。

從局部空間自相關(guān)來看，生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的相關(guān)性以不顯著為主，高價(jià)值-低風(fēng)險(xiǎn)和低價(jià)值-高風(fēng)險(xiǎn)次之，緊接著是低價(jià)值-低風(fēng)險(xiǎn)，相比較而言，高價(jià)值-高風(fēng)險(xiǎn)的占比最小。不顯著是指該區(qū)域的生態(tài)風(fēng)險(xiǎn)等級(jí)與生態(tài)系統(tǒng)服務(wù)價(jià)值不存在顯著的集聚關(guān)系。

高價(jià)值-低風(fēng)險(xiǎn)代表此區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值較高，而生態(tài)風(fēng)險(xiǎn)等級(jí)較低，主要分布在廣西的北部、西北部、西南部以及東北部的山地丘陵，如九萬大山、十萬大山等。這些區(qū)域地形復(fù)雜、降水較多、氣候濕潤、植被茂密、物種豐富、生境質(zhì)量較高，優(yōu)越的自然條件造就很高的生態(tài)系統(tǒng)服務(wù)價(jià)值，而人為活動(dòng)對(duì)此區(qū)域的影響較小，因此其生態(tài)風(fēng)險(xiǎn)較低。

低價(jià)值-高風(fēng)險(xiǎn)代表此區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值較低，而生態(tài)風(fēng)險(xiǎn)較高，主要分布在廣西的中部、西南部和沿海地帶。這些區(qū)域建設(shè)用地較多，生態(tài)系統(tǒng)服務(wù)價(jià)值較高的林地和草地被開發(fā)成建設(shè)用地，因此生態(tài)系統(tǒng)服務(wù)價(jià)值低。同時(shí)，這些地區(qū)大都分布耕地、未利用地，夾雜一些林地和草地，耕作活動(dòng)對(duì)生態(tài)系統(tǒng)的影響較大，因此生態(tài)風(fēng)險(xiǎn)較高。

低價(jià)值-低風(fēng)險(xiǎn)代表此區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值較低，同時(shí)生態(tài)風(fēng)險(xiǎn)較低，其分布與建設(shè)用地相似，主要分布在城市及其周邊地區(qū)，尤其是廣西的北部邊緣和沿海地帶。這些區(qū)域建設(shè)用地較多，人類活動(dòng)對(duì)自然環(huán)境的影響較大，廣西北部邊緣分布有疏林地，因此生態(tài)系統(tǒng)服務(wù)價(jià)值較低且生態(tài)風(fēng)險(xiǎn)較小。同時(shí)，在城鎮(zhèn)化與工業(yè)化的驅(qū)動(dòng)下，建設(shè)用地呈擴(kuò)張趨勢(shì)，低價(jià)值-低風(fēng)險(xiǎn)的區(qū)域也逐漸擴(kuò)張。未來規(guī)劃應(yīng)該重視土地利用結(jié)構(gòu)的調(diào)整，保護(hù)生態(tài)用地，遏制建設(shè)用地的無序擴(kuò)張，多利用存量用地，維護(hù)生態(tài)系統(tǒng)的可持續(xù)發(fā)展。

高價(jià)值-高風(fēng)險(xiǎn)表示此區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值較高，同時(shí)生態(tài)風(fēng)險(xiǎn)較高，主要分布在南寧市、賀州市、百色市、來賓市以及沿海地區(qū)，這些區(qū)域大都分布濕地及水域（左江、右江、郁江、賀江、紅水河、邕江等），且可能位于水域和未利用地的交界處。水域自然資源豐富，生態(tài)系統(tǒng)服務(wù)價(jià)值較高，但容易受到人類活動(dòng)或氣候的影響，生態(tài)敏感性較高，因此生態(tài)風(fēng)險(xiǎn)等級(jí)也較高。水域逐漸向建設(shè)用地轉(zhuǎn)變，水資源供給為這些區(qū)域創(chuàng)造高生態(tài)服務(wù)價(jià)值，同時(shí)人為干擾或自然因素產(chǎn)生較高的生態(tài)風(fēng)險(xiǎn)。

2000—2020年，高價(jià)值-高風(fēng)險(xiǎn)區(qū)的面積逐漸收縮，相關(guān)性轉(zhuǎn)變?yōu)椴伙@著。廣西積極實(shí)施封山育林、石漠化治理、退耕還林還草等政策，使得這些區(qū)域的土壤保持、生境質(zhì)量、碳固定、水源涵養(yǎng)等生態(tài)服務(wù)功能得以提高，因此生態(tài)風(fēng)險(xiǎn)緩慢降低，生態(tài)系統(tǒng)服務(wù)價(jià)值緩慢增高，表現(xiàn)為相關(guān)性不顯著。生態(tài)系統(tǒng)服務(wù)價(jià)值提升，表明生態(tài)環(huán)境變好。需要特別注意的是，生態(tài)風(fēng)險(xiǎn)等級(jí)為高風(fēng)險(xiǎn)，表明區(qū)域生態(tài)系統(tǒng)脆弱，要采取相應(yīng)的保護(hù)和修復(fù)措施。生態(tài)系統(tǒng)服務(wù)價(jià)值高且生態(tài)風(fēng)險(xiǎn)也較高的區(qū)域是廣西生態(tài)系統(tǒng)保護(hù)和治理的重點(diǎn)對(duì)象。這些區(qū)域呈現(xiàn)出較高的生態(tài)價(jià)值，同時(shí)面臨高生態(tài)風(fēng)險(xiǎn)的威脅，因此要重點(diǎn)關(guān)注生態(tài)系統(tǒng)和土地利用，降低其生態(tài)風(fēng)險(xiǎn)，保持較高的生態(tài)系統(tǒng)服務(wù)價(jià)值。

4 結(jié)論

2000—2020年，廣西生態(tài)系統(tǒng)服務(wù)價(jià)值表現(xiàn)出增減起伏的變化趨勢(shì)，從2000年的25 240.24億元增長到2020年的25 251.56億元，整體增幅只有0.05%，總體沒有出現(xiàn)大幅度的增長，而是偏向平緩。其間，高風(fēng)險(xiǎn)區(qū)面積增加幅度最大，為1.92%，而中風(fēng)險(xiǎn)區(qū)面積減少幅度最大，達(dá)到2.61%。2000—2020年，廣西生態(tài)風(fēng)險(xiǎn)空間格局變化不大，但是各區(qū)域生態(tài)風(fēng)險(xiǎn)空間差異比較明顯，總體呈南高北低的分布格局，生態(tài)風(fēng)險(xiǎn)區(qū)面積表現(xiàn)出增減起伏的變化。

從數(shù)量相關(guān)性來看，2000年、2005年、2010年、2015年和2020年的斯皮爾曼相關(guān)系數(shù)分別為-0.725 4、-0.734 6、-0.757 1、-0.764 5和-0.785 1，說明生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)具有較強(qiáng)的相關(guān)性。從空間相關(guān)性來看，2000年、2005年、2010年、2015年和2020年的莫蘭指數(shù)分別為-0.145 2、-0.173 6、-0.194 1、-0.223 4和-0.249 5，相關(guān)性系數(shù)P均小于0.01，說明生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)在空間上呈聚集分布，且呈顯著的負(fù)相關(guān)。從局部空間自相關(guān)分析結(jié)果來看，生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)的相關(guān)性以不顯著為主，高價(jià)值-低風(fēng)險(xiǎn)、低價(jià)值-高風(fēng)險(xiǎn)次之，而低價(jià)值-低風(fēng)險(xiǎn)、高價(jià)值-高風(fēng)險(xiǎn)占比相對(duì)較小。高價(jià)值-低風(fēng)險(xiǎn)主要分布在廣西的北部、西北部、西南部以及東北部的山地丘陵；低價(jià)值-高風(fēng)險(xiǎn)主要分布在廣西的中部、西南部和沿海地帶；低價(jià)值-低風(fēng)險(xiǎn)主要分布在城市及其周邊地區(qū)，尤其是廣西北部邊緣和沿海地區(qū)；高價(jià)值-高風(fēng)險(xiǎn)主要分布在南寧市、賀州市、百色市、來賓市以及沿海地區(qū)。

生態(tài)系統(tǒng)組成復(fù)雜，當(dāng)前，生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估還沒有形成規(guī)范的標(biāo)準(zhǔn)體系。所以，本文結(jié)合已有研究成果和修正規(guī)則，并基于廣西實(shí)際情況，修正生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估模型的相關(guān)系數(shù)。但是，研究區(qū)范圍較大，修正生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子時(shí)，可能考慮不夠周全，有待進(jìn)一步完善研究方法。

參考文獻(xiàn)

1 壽飛云，李卓飛，黃 璐，等.基于生態(tài)系統(tǒng)服務(wù)供求評(píng)價(jià)的空間分異特征與生態(tài)格局劃分：以長三角城市群為例[J].生態(tài)學(xué)報(bào)，2020（9）：2813-2826.

2 鄔紫荊，曾 輝.基于Meta分析的中國森林生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估[J].生態(tài)學(xué)報(bào)，2021（14）：5533-5545.

3 周文昌，張 維，胡興宜，等.湖北省濕地生態(tài)系統(tǒng)的服務(wù)價(jià)值評(píng)估[J].水土保持通報(bào)，2021（3）：305-311.

4 趙 正，韓 鋒，侯一蕾.基于Meta回歸方法的中國城市森林生態(tài)系統(tǒng)服務(wù)功能價(jià)值再評(píng)估[J].長江流域資源與環(huán)境，2021（1）：64-75.

5 尤海舟，王 超，畢 君.河北省森林生態(tài)系統(tǒng)固碳釋氧服務(wù)功能價(jià)值評(píng)估[J].西部林業(yè)科學(xué)，2017（4）：121-127.

6 姜棟棟，楊 帆，馬偉波，等.大婁山區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值變化與情景預(yù)測(cè)[J].環(huán)境科學(xué)研究，2022（7）：1670-1680.

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

8 于 媛，李明玉，韓旭龍，等.邊境地區(qū)景觀格局演變及生態(tài)系統(tǒng)服務(wù)價(jià)值響應(yīng)：以延邊朝鮮族自治州為例[J].水土保持研究，2021（1）：315-322.

9 樂榮武，李 巍，周思楊，等.呼包鄂榆城市群生態(tài)系統(tǒng)服務(wù)價(jià)值驅(qū)動(dòng)因素及其交互效應(yīng)研究[J].生態(tài)學(xué)報(bào)，2023（23）：1-14.

10 Sannigrahi S，Bhatt S，Rahmat S，et al.Estimating global ecosystem service values and its response to land surface dynamics during 1995—2015[J].Journal of Environmental Management，2018，223：115-131.

11 包 妲，李 杰，劉 鋒，等.1995—2015年緬甸土地利用/覆被變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[J].生態(tài)學(xué)報(bào)，2021（17）：6960-6969.

12 Xu Y，Xie Y，Wu X，et al.Evaluating temporal-spatial variations of wetland ecosystem service value in China during 1990—2020 from the donor side based on cosmic exergy[J].Journal of Cleaner Production，2023，414：137485.

13 戚麗萍，欒兆擎，魏 勉，等.基于土地利用的江蘇省各市生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化研究[J].南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2023（4）：200-208.

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

15 丘海紅，胡寶清，張 澤.基于土地利用變化的廣西近20年生態(tài)系統(tǒng)服務(wù)價(jià)值研究[J].環(huán)境工程技術(shù)學(xué)報(bào)，2022（5）：1455-1465.

16 Qiu H，Hu B，Zhang Z.Impacts of land use change on ecosystem service value based on SDGs report：taking Guangxi as an example[J].Ecological Indicators，2021，133：108366.

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

18 彭文君，舒英格.基于GIS的石漠化山區(qū)縣域土地利用空間變化的生態(tài)風(fēng)險(xiǎn)測(cè)度[J].水土保持研究，2018（1）：342-348.

19 謝 玲，鄧曉軍，盧月燕，等.廣西石漠化地區(qū)土地利用空間變化的生態(tài)風(fēng)險(xiǎn)研究[J].中國農(nóng)業(yè)資源與區(qū)劃，2019（8）：113-121.

20 臧淑英，梁 欣，張思沖.基于GIS的大慶市土地利用生態(tài)風(fēng)險(xiǎn)分析[J].自然災(zāi)害學(xué)報(bào)，2005（4）：141-145.

21 賈占華，谷國鋒.東北地區(qū)人口分布的時(shí)空演變特征及影響因素[J].經(jīng)濟(jì)地理，2016（12）：60-68.

收稿日期：2024-01-02

作者簡(jiǎn)介：黃靈麗（2000—），女，廣西都安人，碩士研究生。研究方向：土地利用與區(qū)域發(fā)展。