李洪慶，張俊紅，楊 瑀

（河海大學(xué) 公共管理學(xué)院，江蘇 南京 210000）

0 引言

中國(guó)是世界上生態(tài)脆弱區(qū)分布面積最大、脆弱生態(tài)類型最多、生態(tài)脆弱性表現(xiàn)最明顯的國(guó)家之一，中度以上生態(tài)脆弱區(qū)占國(guó)土總面積的55%左右，生態(tài)脆弱區(qū)的保護(hù)與發(fā)展壓力巨大。由于自身生態(tài)系統(tǒng)結(jié)構(gòu)穩(wěn)定性較差，抵抗外來(lái)干擾及維持自身穩(wěn)定性能力較弱，生態(tài)脆弱區(qū)極易產(chǎn)生生態(tài)退化問題且難以自我修復(fù)［1］，表現(xiàn)為對(duì)外部自然環(huán)境變化和人類活動(dòng)具有高度的敏感性［2］，在極端情況下容易發(fā)生系統(tǒng)性突變［3］。當(dāng)今全球氣候變暖、二氧化碳濃度升高、極端天氣頻發(fā)和人類活動(dòng)頻繁，這些因素對(duì)生態(tài)脆弱區(qū)的生態(tài)系統(tǒng)結(jié)構(gòu)與分布、演替與自我調(diào)節(jié)能力［4］、生產(chǎn)力［5］和碳存儲(chǔ)［6］、區(qū)域生態(tài)平衡［7］，甚至全球社會(huì)經(jīng)濟(jì)發(fā)展與生態(tài)安全產(chǎn)生了極大的影響［8］。在生態(tài)脆弱區(qū)開展生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估、生態(tài)風(fēng)險(xiǎn)評(píng)估、識(shí)別生態(tài)風(fēng)險(xiǎn)區(qū)對(duì)區(qū)域生態(tài)保護(hù)、資源管理、國(guó)土空間修復(fù)非常必要［9-10］。

生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)能夠分析評(píng)價(jià)生態(tài)環(huán)境狀況與時(shí)空分布特征，而風(fēng)險(xiǎn)區(qū)識(shí)別則可精確地展示生態(tài)脆弱區(qū)的脆弱類型、存在問題、高敏感性區(qū)域［11］。隨著對(duì)生態(tài)脆弱區(qū)開發(fā)和保護(hù)的重視，生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)與識(shí)別受到了學(xué)者們?cè)絹?lái)越多的關(guān)注［12］；在評(píng)價(jià)尺度上從全國(guó)、區(qū)域尺度逐步精確到小流域尺度［13］；識(shí)別對(duì)象從單一風(fēng)險(xiǎn)源轉(zhuǎn)向多重風(fēng)險(xiǎn)源［14］；評(píng)價(jià)對(duì)象也從單一類型的生態(tài)系統(tǒng)轉(zhuǎn)到多種生態(tài)系統(tǒng)的空間鑲嵌［15］；研究區(qū)域不僅涵蓋了東部海岸帶［16］、西藏［17］、貴州喀斯特［18］和新疆干旱區(qū)［19］等典型生態(tài)脆弱區(qū)，而且還向城市內(nèi)部空間［20］、城市群之間拓展［21］。生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)方法以風(fēng)險(xiǎn)源—匯評(píng)價(jià)［22］和景觀指數(shù)評(píng)價(jià)方法［23］為主，前者側(cè)重于風(fēng)險(xiǎn)源、風(fēng)險(xiǎn)受體的識(shí)別，評(píng)價(jià)過程受主觀意識(shí)的影響較大，后者更加注重區(qū)域時(shí)空異質(zhì)性和尺度性，具有空間可視化的突出優(yōu)勢(shì)［12，24］；景觀指數(shù)評(píng)價(jià)方法結(jié)合空間疊加方法可以完成生態(tài)風(fēng)險(xiǎn)的空間評(píng)價(jià)，因此用景觀干擾度和脆弱度指數(shù)構(gòu)建的景觀生態(tài)風(fēng)險(xiǎn)評(píng)估模型的應(yīng)用較為廣泛［25］。

當(dāng)前大部分研究是在風(fēng)險(xiǎn)脆弱性評(píng)價(jià)的基礎(chǔ)上直接完成風(fēng)險(xiǎn)空間布局，重點(diǎn)突出脆弱性。山南市屬于青藏高原典型的復(fù)合侵蝕生態(tài)脆弱區(qū)，在全球氣候變化的背景下該市的生態(tài)系統(tǒng)表現(xiàn)出高脆弱性和敏感性，且隨著該市人口迅速增長(zhǎng)、經(jīng)濟(jì)快速發(fā)展、城鎮(zhèn)化加速推進(jìn)，當(dāng)?shù)厝祟惢顒?dòng)對(duì)區(qū)域生態(tài)安全的影響愈加強(qiáng)烈，對(duì)構(gòu)建西部生態(tài)安全屏障構(gòu)成了威脅［26］。生態(tài)系統(tǒng)服務(wù)在一定程度上可真實(shí)地反映生態(tài)功能的強(qiáng)弱，并正相關(guān)地表征生態(tài)安全的狀況［27-28］，因此本研究將生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)相結(jié)合，分析了兩者的空間關(guān)系，不僅突出了生態(tài)敏感性高的區(qū)域，而且揭示了風(fēng)險(xiǎn)區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值貢獻(xiàn)等級(jí)，可以為風(fēng)險(xiǎn)區(qū)的保護(hù)、修復(fù)和利用提供科學(xué)依據(jù)。

1 研究方法

1.1 研究區(qū)概況

山南市屬西藏自治區(qū)下轄地級(jí)市（圖1），位于岡底斯山至念青唐古拉山以南、雅魯藏布江干流中下游地區(qū)，北接拉薩，西與日喀則毗鄰，東與林芝相連，南與印度、不丹兩國(guó)接壤，地處東經(jīng)90°14′～94°22′、北緯27°08′～29°47′之間，面積7.93萬(wàn)km2。山南市轄1個(gè)市轄區(qū)（乃東區(qū)）和11個(gè)縣（扎囊縣、貢嘎縣、桑日縣、瓊結(jié)縣、曲松縣、措美縣、洛扎縣、加查縣、隆子縣、錯(cuò)那縣、浪卡子縣），屬于典型的藏南谷地，地勢(shì)自西向東逐漸降低，平均海拔約3700 m。氣候?yàn)闇貛Ц珊敌詺夂颍昃邓坎坏?50 mm，雨季多集中在6—9月。在2020年，山南市實(shí)現(xiàn)地區(qū)生產(chǎn)總值215.4億元，同比增長(zhǎng)了7.9%，其中第一產(chǎn)業(yè)增加值8.14億元，同比增長(zhǎng)了7.8%；第二產(chǎn)業(yè)增加值115.51億元，同比增長(zhǎng)了17.1%；第三產(chǎn)業(yè)增加值91.75億元，同比下降了0.5%。

圖1 研究區(qū)概況及評(píng)價(jià)單元?jiǎng)澐?/p>

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

本文所采用的1980、2000、2020年3期土地利用矢量數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心，空間分辨率為30 m，包括耕地、林地、草地、水域、建設(shè)用地、未利用地共6種一級(jí)地類和25種二級(jí)地類。計(jì)算各生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值所用的糧食產(chǎn)量、糧食價(jià)格數(shù)據(jù)來(lái)源于《山南市統(tǒng)計(jì)年鑒》《全國(guó)農(nóng)產(chǎn)品成本收益資料匯編》。

利用ArcGIS軟件，采用5 km×5 km的單元網(wǎng)格對(duì)山南市土地利用景觀格局進(jìn)行網(wǎng)格化重采樣，得到3264個(gè)評(píng)價(jià)單元（圖1）；基于評(píng)價(jià)單元計(jì)算生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)指數(shù)，并將數(shù)值賦值于評(píng)價(jià)單元中心點(diǎn)，進(jìn)行克里格插值，分析其空間格局分布特征。

1.3 研究方法

1.3.1 “三生空間”分類體系 同一用地類型可能兼顧多種功能，從而展示出不同生態(tài)系統(tǒng)服務(wù)價(jià)值，例如草地資源涵蓋了高覆蓋、低覆蓋類型，很難從利用功能視角開展土地分類工作。為了更科學(xué)精確地計(jì)算生態(tài)系統(tǒng)服務(wù)價(jià)值，同時(shí)考慮國(guó)土空間開發(fā)與利用，本研究從“三生空間”角度出發(fā)分析了山南市的關(guān)鍵生態(tài)風(fēng)險(xiǎn)區(qū)?！叭臻g”具有空間功能復(fù)合性、空間尺度差異性及空間用地異質(zhì)性等特征。借鑒前人的研究成果［29-30］，本研究建立的“三生空間”分類體系見表1。

表1 “三生空間”分類體系

1.3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估 參考謝高地等［31］修訂的中國(guó)生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量表，結(jié)合研究區(qū)的實(shí)際情況對(duì)部分當(dāng)量因子進(jìn)行修正，完成了適用于山南市“三生空間”生態(tài)系統(tǒng)服務(wù)價(jià)值的當(dāng)量因子表。1個(gè)標(biāo)準(zhǔn)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值量為平均糧食單產(chǎn)市場(chǎng)價(jià)格的1/7［32］，根據(jù)山南市多年來(lái)糧食平均產(chǎn)量5642.27 kg/hm2和平均糧食價(jià)格1.38元/kg，計(jì)算出山南市區(qū)單位面積農(nóng)田生態(tài)系統(tǒng)糧食作物的經(jīng)濟(jì)價(jià)值為1112.33元/hm2，完成了適合山南市的生態(tài)系統(tǒng)服務(wù)機(jī)制當(dāng)量表，涉及供給服務(wù)、調(diào)節(jié)服務(wù)、支持服務(wù)和文化服務(wù)四大類，共計(jì)11個(gè)二級(jí)類。

1.3.3 生態(tài)風(fēng)險(xiǎn)評(píng)價(jià) 本研究采用指數(shù)評(píng)價(jià)法進(jìn)行生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)。選取了景觀干擾度指數(shù)、脆弱度指數(shù)、損失度指數(shù)計(jì)算景觀生態(tài)風(fēng)險(xiǎn)指數(shù)，計(jì)算公式為：

式（1）中：ERI為景觀生態(tài)風(fēng)險(xiǎn)指數(shù)；Ri為景觀損失度指數(shù)；Fi為景觀脆弱度指數(shù)；Si為景觀干擾度指數(shù)；Ci為景觀破碎度指數(shù)；Ni為景觀分離度指數(shù)；Di為景觀優(yōu)勢(shì)度指數(shù)；Aki為第k個(gè)風(fēng)險(xiǎn)小區(qū)內(nèi)景觀i的面積；Ak為第k個(gè)風(fēng)險(xiǎn)小區(qū)的總面積；ni為景觀i的斑塊數(shù)；n為景觀斑塊總數(shù)量；Ai為景觀類型i的斑塊面積；A為景觀總面積。

借鑒相關(guān)學(xué)者的研究成果［20，33］，將式（1）中的a、b、c、d、e值分別設(shè)為0.5、0.3、0.2、0.6、0.4。不同空間的脆弱度賦值如下：其他生態(tài)空間6，水域生態(tài)空間5，農(nóng)業(yè)生產(chǎn)空間4，草地生態(tài)空間3，林地生態(tài)空間2，工業(yè)生產(chǎn)空間1，城鎮(zhèn)生活空間1，農(nóng)村生活空間1。

1.3.4 雙變量的空間自相關(guān)分析 采用雙變量空間分析模型，利用全局自相關(guān)系數(shù)度量生態(tài)系統(tǒng)服務(wù)和生態(tài)風(fēng)險(xiǎn)整體空間相關(guān)特征，計(jì)算公式為：

式（2）中：I為生態(tài)系統(tǒng)服務(wù)價(jià)值與景觀生態(tài)風(fēng)險(xiǎn)指數(shù)的雙變量全局自相關(guān)系數(shù)；yi，s和yi，r分別為第i評(píng)價(jià)單元單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值和景觀生態(tài)風(fēng)險(xiǎn)指數(shù)；σs和σr為方差；Wij為空間權(quán)重。

2 結(jié)果與分析

2.1 山南市“三生空間”的時(shí)空演變特征

結(jié)合3期土地利用數(shù)據(jù)，根據(jù)“三生空間”分類體系，得出了山南市1980、2000、2020年“三生空間”的時(shí)間演變圖（圖2）和空間格局變化特征（表2）。

表2 1980—2020年山南市“三生空間”的結(jié)構(gòu)變化

圖2 1980—2020年山南市“三生空間”的時(shí)間演變

從空間上分析，2020年山南市“三生空間”以生態(tài)空間為主，約占區(qū)域總面積的98%；生產(chǎn)空間和生活空間僅占2%左右，主要分布于雅魯藏布江流域山南寬谷段。從時(shí)間動(dòng)態(tài)上分析，生產(chǎn)空間中農(nóng)業(yè)生產(chǎn)和工業(yè)生產(chǎn)空間均呈上升趨勢(shì)，其中農(nóng)業(yè)生產(chǎn)空間由1980年的608.42 km2增加到2020年的1389.59 km2，變化率為128.29%，主要由林地生態(tài)空間和草地生態(tài)空間分別轉(zhuǎn)入528.06和515.87 km2，分布在山南寬谷段和錯(cuò)那縣東南部。1980—2020年生活空間中城鎮(zhèn)生活空間和農(nóng)村生活空間分別增加了20.65和2.45 km2，變化率分別為882.48%、25.47%；生態(tài)空間內(nèi)部土地類型轉(zhuǎn)移較為明顯，主要由草地生態(tài)空間轉(zhuǎn)為林地生態(tài)空間，由林地生態(tài)空間轉(zhuǎn)為草地生態(tài)空間，以及由草地生態(tài)空間轉(zhuǎn)為其他生態(tài)空間，轉(zhuǎn)入面積分別為8023.86、5634.79和5588.74 km2，主要發(fā)生在雅魯藏布江山南寬谷北岸、隆子縣、錯(cuò)那縣、浪卡子縣和措美縣。生產(chǎn)和生活空間的擴(kuò)張擠壓了生態(tài)空間，導(dǎo)致40年間草地、水域和林地的生態(tài)空間分別減少了1321.09、382.77、2864.33 km2。上述結(jié)果說明隨著人口的迅速增加以及工業(yè)化、城鎮(zhèn)化的加速推進(jìn)，生產(chǎn)、生活用地需求不斷增長(zhǎng)，導(dǎo)致人地矛盾加劇，生態(tài)空間被不斷壓縮，對(duì)脆弱區(qū)生態(tài)系統(tǒng)造成了一定的壓力。

2.2 山南市生態(tài)系統(tǒng)服務(wù)價(jià)值的評(píng)價(jià)

依據(jù)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表，計(jì)算出1980—2020年山南市“三生空間”中各類型空間的生態(tài)系統(tǒng)服務(wù)價(jià)值（表3），山南市1980、2000、2020年生態(tài)系統(tǒng)服務(wù)價(jià)值總量分別為1294.32億、1281.34億、1129.26億元，總體呈下降態(tài)勢(shì)，總量減少了165.06億元，且2000—2020年生態(tài)系統(tǒng)服務(wù)價(jià)值的下降幅度遠(yuǎn)大于1980—2000年的下降幅度；生態(tài)系統(tǒng)服務(wù)價(jià)值主要由草地、林地和水域生態(tài)空間貢獻(xiàn)，占總額的98%以上。

表3 1980—2020年山南市生態(tài)系統(tǒng)服務(wù)價(jià)值的變化

在1980—2020年期間，農(nóng)業(yè)生產(chǎn)空間的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化率有所上升，為128.49%。根據(jù)時(shí)空格局演變分析，山南市城鎮(zhèn)化的快速發(fā)展和農(nóng)業(yè)、工業(yè)用地的增加侵占了林地和草地，工業(yè)、城鎮(zhèn)和農(nóng)村生產(chǎn)生活空間的增加對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生了更多的負(fù)向作用，而農(nóng)業(yè)生產(chǎn)空間的生態(tài)系統(tǒng)服務(wù)價(jià)值增加明顯。

生態(tài)空間價(jià)值量由1980年的1291.62億元下降到2020年的1123.08億元，減少了13.05%。經(jīng)濟(jì)社會(huì)的發(fā)展侵占了生態(tài)用地空間，導(dǎo)致林地生態(tài)空間和草地生態(tài)空間的生態(tài)系統(tǒng)服務(wù)供給值下降，分別下降了32.12%、12.94%。其他生態(tài)空間面積的增加導(dǎo)致其生態(tài)系統(tǒng)服務(wù)價(jià)值上升較多。水域空間的生態(tài)服務(wù)價(jià)值不降反升，增加了21.06%，主要原因是雅魯藏布江干流河道治理工程的實(shí)施導(dǎo)致河渠湖泊以及灘地面積不斷增加，表明近20年干流河道治理效果十分顯著。

利用克里格插值法對(duì)3264個(gè)評(píng)價(jià)單元進(jìn)行空間插值，然后基于自然間斷點(diǎn)分級(jí)法結(jié)合研究區(qū)的實(shí)際情況，將生態(tài)系統(tǒng)服務(wù)價(jià)值分為低（＜1000萬(wàn)元）、較低（1000萬(wàn)～2000萬(wàn)元）、中（2000萬(wàn)～3000萬(wàn)元）、較高（3000萬(wàn)～4000萬(wàn)元）、高（＞4000萬(wàn)元）5個(gè)等級(jí)（圖3）。

圖3 1980—2020年山南市生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布

空間分析結(jié)果表明：山南市生態(tài)系統(tǒng)服務(wù)價(jià)值空間分異明顯，整體上呈現(xiàn)南高北低的空間分布格局。山南市單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值以高等級(jí)為主，分布在南部林區(qū)，面積占比達(dá)到40%以上；低、較低價(jià)值等級(jí)面積比例分別由1980年的2.34%、14.08%增加到2020年的11.98%、22.43%，主要分布在山南市中部，植被類型大部分為草地和未利用地，植被退化、森林面積較少、人類活動(dòng)頻繁是導(dǎo)致該區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值降低的主要原因。北部雅魯藏布江山南寬谷流域生態(tài)系統(tǒng)服務(wù)價(jià)值的范圍明顯擴(kuò)大，說明山南寬谷多年來(lái)植樹造林種草、土地沙化治理等一系列生態(tài)保護(hù)措施取得了一定的成效。

2.3 山南市生態(tài)風(fēng)險(xiǎn)的評(píng)價(jià)

利用Fragstats軟件計(jì)算出山南市各類景觀指標(biāo)；采用克里格插值法對(duì)3264個(gè)評(píng)價(jià)單元進(jìn)行空間插值，并基于自然間斷點(diǎn)分級(jí)法結(jié)合研究區(qū)的實(shí)際情況，將生態(tài)風(fēng)險(xiǎn)分為低（＜0.02）、較低（0.02～0.03）、中（0.03～0.04）、較 高（0.04～0.05）、高（＞0.05）5個(gè)等級(jí)（圖4）。

從圖4可以看出，山南市的生態(tài)風(fēng)險(xiǎn)總體上呈現(xiàn)出南低北高的特征。具體來(lái)說，1980和2000年山南市以中風(fēng)險(xiǎn)和低生態(tài)風(fēng)險(xiǎn)等級(jí)為主；2020年以較低風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)等級(jí)為主，但是高風(fēng)險(xiǎn)區(qū)明顯增多，與1980年相比，2020年高風(fēng)險(xiǎn)區(qū)面積急劇增加了2590.23 km2，增長(zhǎng)率為839.61%，高風(fēng)險(xiǎn)區(qū)主要分布于雅魯藏布江兩岸、羊卓雍措、乃東區(qū)及雅礱河、中部林草交錯(cuò)區(qū)，主要因?yàn)樵搮^(qū)域經(jīng)濟(jì)、社會(huì)發(fā)展迅速，人類活動(dòng)活躍，景觀類型主要為生活用地和生產(chǎn)用地。南部部分地區(qū)由低風(fēng)險(xiǎn)轉(zhuǎn)為較低風(fēng)險(xiǎn)，風(fēng)險(xiǎn)區(qū)面積增加，區(qū)域生態(tài)風(fēng)險(xiǎn)等級(jí)有所提升，可能存在林地、草地輕微退化的現(xiàn)象。2020年中部和北部地區(qū)的生態(tài)風(fēng)險(xiǎn)出現(xiàn)了明顯分化，除了高生態(tài)風(fēng)險(xiǎn)區(qū)面積增加外，大面積的中風(fēng)險(xiǎn)地區(qū)轉(zhuǎn)變?yōu)榈惋L(fēng)險(xiǎn)地區(qū)，表明山南市生態(tài)文明建設(shè)、圍欄育草、封山育林等一系列生態(tài)建設(shè)工程的大力開展和生態(tài)補(bǔ)償政策的實(shí)施促進(jìn)了生態(tài)恢復(fù)，但人類活動(dòng)的強(qiáng)烈干擾、全球氣候變暖、未利用地斑塊破碎化等使得部分區(qū)域的生態(tài)風(fēng)險(xiǎn)等級(jí)升高。

圖4 1980—2020年山南市生態(tài)風(fēng)險(xiǎn)的空間分布

從“三生空間”視角分析：生產(chǎn)空間的風(fēng)險(xiǎn)性呈現(xiàn)升高趨勢(shì)，其中農(nóng)業(yè)生產(chǎn)空間中較高風(fēng)險(xiǎn)區(qū)面積明顯增多，且高風(fēng)險(xiǎn)區(qū)面積減少不明顯；工業(yè)生產(chǎn)空間中不僅風(fēng)險(xiǎn)區(qū)面積增加，而且較高和高風(fēng)險(xiǎn)區(qū)面積增加顯著。生活空間的總面積增加，風(fēng)險(xiǎn)等級(jí)總體上略有升高，其中城鎮(zhèn)和農(nóng)村生活空間中的低風(fēng)險(xiǎn)區(qū)面積占比大，但城鎮(zhèn)生活空間中較高和高風(fēng)險(xiǎn)區(qū)面積增加不明顯，而農(nóng)村生活空間中高風(fēng)險(xiǎn)區(qū)轉(zhuǎn)為較低風(fēng)險(xiǎn)區(qū)的面積增多。生態(tài)空間中不同生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)變較多，風(fēng)險(xiǎn)性總體偏好，但是高風(fēng)險(xiǎn)區(qū)面積明顯增加，其中林地和草地均以低風(fēng)險(xiǎn)區(qū)面積增加為主，林地的較高風(fēng)險(xiǎn)區(qū)面積有所增加但占比不大，草地的較高和高風(fēng)險(xiǎn)區(qū)面積變化較?。凰蚝推渌鷳B(tài)空間的較高和高風(fēng)險(xiǎn)區(qū)面積明顯增多，由0.26%上升到12.19%。綜上所述，生產(chǎn)空間的風(fēng)險(xiǎn)等級(jí)提高；生活空間的面積增加，但風(fēng)險(xiǎn)等級(jí)變化不明顯；生態(tài)空間中的水域和其他生態(tài)空間的風(fēng)險(xiǎn)等級(jí)較高。

2.4 山南市生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)相關(guān)性分析

基于GeoDa軟件，采用雙變量空間自相關(guān)模型，建立空間權(quán)重矩陣，對(duì)山南市生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)進(jìn)行了雙變量空間關(guān)聯(lián)特征分析。1980、2000、2020年全局自相關(guān)Moran I指數(shù)分別為-0.228、-0.228、-0.096，表明生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)存在一定程度的負(fù)向空間關(guān)聯(lián)性，即生態(tài)風(fēng)險(xiǎn)越高的區(qū)域，生態(tài)系統(tǒng)服務(wù)價(jià)值越低，在空間分布上也表現(xiàn)出高風(fēng)險(xiǎn)低價(jià)值的區(qū)域范圍大。

基于生態(tài)系統(tǒng)服務(wù)價(jià)值，將生態(tài)風(fēng)險(xiǎn)區(qū)分為不顯著、高價(jià)值—高風(fēng)險(xiǎn)、高價(jià)值—低風(fēng)險(xiǎn)、低價(jià)值—高風(fēng)險(xiǎn)、低價(jià)值—低風(fēng)險(xiǎn)5種模式（圖5）。從總體上分析，低價(jià)值—高風(fēng)險(xiǎn)區(qū)面積占比由1980年的37%減少到2020年的21%；高價(jià)值—高風(fēng)險(xiǎn)區(qū)面積明顯增多，占比由1980年的2%增加到2020年的4%；高價(jià)值—低風(fēng)險(xiǎn)區(qū)面積呈下降趨勢(shì)，占比由1980年的26%減少到2020年的22%。說明生態(tài)系統(tǒng)服務(wù)低價(jià)值區(qū)的生態(tài)風(fēng)險(xiǎn)明顯降低，但是高價(jià)值區(qū)的生態(tài)風(fēng)險(xiǎn)不斷上升，山南市的生態(tài)安全面臨較大的壓力。

圖5 山南市生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)的空間關(guān)系

高價(jià)值—高風(fēng)險(xiǎn)區(qū)面積占比由1980年的2%增加到2020年的4%，分布于雅魯藏布江山南寬谷流域及浪卡子縣的羊卓雍措，主要的景觀類型是水域以及水域與未利用地的交界處，大部分屬于生態(tài)空間中的水域生態(tài)空間；水域極易受到氣候變化的影響，在全球氣候變暖的大背景下，水域生態(tài)空間面臨著面積萎縮、景觀破碎化和向未利用地轉(zhuǎn)化的高風(fēng)險(xiǎn)。高價(jià)值—高風(fēng)險(xiǎn)區(qū)對(duì)山南市生態(tài)系統(tǒng)保護(hù)具有極為重要的作用，既有高生態(tài)系統(tǒng)服務(wù)價(jià)值表現(xiàn)，但同時(shí)又面臨著高度生態(tài)風(fēng)險(xiǎn)，該區(qū)域是未來(lái)生態(tài)環(huán)境保護(hù)和土地利用的重點(diǎn)關(guān)注區(qū)域。對(duì)于高價(jià)值—高風(fēng)險(xiǎn)區(qū)域，建議未來(lái)加強(qiáng)對(duì)水域生態(tài)景觀的保護(hù)，適當(dāng)減少人類活動(dòng)的干擾，減少景觀破碎化程度，同時(shí)加強(qiáng)水域周邊尤其是雅魯藏布江和羊卓雍措?yún)^(qū)域林地、草地的保護(hù)與修復(fù)，避免生態(tài)景觀向未利用地和建設(shè)用地轉(zhuǎn)化。

高價(jià)值—低風(fēng)險(xiǎn)區(qū)主要位于錯(cuò)那縣和加查縣南部，面積占比由1980年的26%減少到2020年的22%，減少范圍的空間分布十分明顯和集中，主要向不顯著區(qū)域轉(zhuǎn)化，具有向高價(jià)值—高風(fēng)險(xiǎn)轉(zhuǎn)化的趨勢(shì)。該區(qū)域的主要景觀類型為林地，其植被覆蓋度高，提供的生態(tài)系統(tǒng)服務(wù)價(jià)值較高，林地生態(tài)系統(tǒng)較為穩(wěn)定，生態(tài)風(fēng)險(xiǎn)較低。但是隨著人類活動(dòng)干擾程度的加劇，山南市南部部分區(qū)域土地利用類型由林地轉(zhuǎn)為草地或者耕地，導(dǎo)致該部分區(qū)域面積明顯減少。對(duì)于高價(jià)值—低風(fēng)險(xiǎn)區(qū)域，建議加強(qiáng)對(duì)林地生態(tài)空間的保護(hù)與管理，強(qiáng)化生態(tài)安全屏障，減少向草地或者耕地的轉(zhuǎn)化，盡量保持原有的生態(tài)安全空間。

低價(jià)值—高風(fēng)險(xiǎn)區(qū)主要位于山南市中部和北部，面積占比由1980年的37%減少到2020年的21%，主要為生態(tài)空間的未利用地和草地，以及生產(chǎn)空間的耕地，其提供的生態(tài)系統(tǒng)服務(wù)價(jià)值較低；結(jié)合山南市風(fēng)險(xiǎn)特征分布圖，減少的區(qū)域向無(wú)風(fēng)險(xiǎn)區(qū)轉(zhuǎn)移，生態(tài)安全等級(jí)有所提升。結(jié)合“三生空間”演變和生態(tài)系統(tǒng)服務(wù)價(jià)值分析，在氣候變暖和人類活動(dòng)加劇的雙重影響下，該區(qū)域的植被退化導(dǎo)致未利用地面積明顯增加，而人類生產(chǎn)生活占用了大量的耕地和草地，導(dǎo)致生態(tài)風(fēng)險(xiǎn)增大；但是隨著山南市圍欄、人工草地等生態(tài)保護(hù)措施和草地生態(tài)補(bǔ)償獎(jiǎng)勵(lì)政策的實(shí)施，草地景觀破碎度下降，草地風(fēng)險(xiǎn)指數(shù)逐漸降低，原低價(jià)值—高風(fēng)險(xiǎn)區(qū)面積明顯縮小。對(duì)于低價(jià)值—高風(fēng)險(xiǎn)區(qū)，建議注重草地質(zhì)量的提升，減緩草地退化速度，適當(dāng)控制生產(chǎn)和生活用地的擴(kuò)張，減少人類活動(dòng)的強(qiáng)烈干擾，提升其生態(tài)系統(tǒng)服務(wù)價(jià)值；同時(shí)也要注重對(duì)生態(tài)系統(tǒng)的保護(hù)，逐漸減少該區(qū)域的面積，提升生態(tài)安全等級(jí)。

低價(jià)值—低風(fēng)險(xiǎn)區(qū)在1980、2000和2020年的面積占比均為5%，主要位于錯(cuò)那縣，主要的景觀類型為疏林地，屬于生態(tài)空間的林地生態(tài)空間，其提供的生態(tài)系統(tǒng)服務(wù)價(jià)值較低，同時(shí)生態(tài)風(fēng)險(xiǎn)較小，空間分布保持較為穩(wěn)定。建議低價(jià)值—低風(fēng)險(xiǎn)區(qū)優(yōu)化空間格局，通過多種生態(tài)建設(shè)措施增加植被覆蓋度，提升生態(tài)系統(tǒng)服務(wù)價(jià)值，促進(jìn)低價(jià)值—低風(fēng)險(xiǎn)區(qū)向高價(jià)值—低風(fēng)險(xiǎn)區(qū)轉(zhuǎn)變。

3 小結(jié)

（1）1980—2020年山南市的農(nóng)業(yè)生產(chǎn)空間面積增加了128.39%，生態(tài)空間依然是主要土地利用類型，約占98%，其中草地生態(tài)空間、林地生態(tài)空間、水域生態(tài)空間面積均有減少，草地生態(tài)空間部分區(qū)域退化且向其他生態(tài)空間轉(zhuǎn)移，“三生空間”的變化深刻影響了山南市生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)風(fēng)險(xiǎn)值。

（2）山南市生態(tài)系統(tǒng)服務(wù)價(jià)值空間分異明顯，整體上呈現(xiàn)南高北低的空間分布格局。40年間山南市生態(tài)系統(tǒng)服務(wù)價(jià)值由1294.32億元降低到1129.26億元，下降的主要原因是中部地區(qū)生活空間增加、草地生態(tài)空間退化、南部地區(qū)林地生態(tài)空間的生態(tài)系統(tǒng)服務(wù)價(jià)值下降。

（3）山南市生態(tài)風(fēng)險(xiǎn)等級(jí)呈現(xiàn)出南低北高的分布特征。40年間山南市生態(tài)風(fēng)險(xiǎn)由中風(fēng)險(xiǎn)和低生態(tài)風(fēng)險(xiǎn)為主向較低和低風(fēng)險(xiǎn)轉(zhuǎn)化，但是高風(fēng)險(xiǎn)區(qū)面積明顯增多，增加了2590.23 km2。山南市的生態(tài)環(huán)境總體上得到了改善，但局部區(qū)域生態(tài)環(huán)境有所惡化，生態(tài)安全面臨一定壓力。

（4）山南市生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)風(fēng)險(xiǎn)存在極顯著負(fù)相關(guān)關(guān)系。至2020年，高價(jià)值—高風(fēng)險(xiǎn)區(qū)面積增加了113.43%，主要分布在雅魯藏布江寬谷流域和其他水域；高價(jià)值—低風(fēng)險(xiǎn)區(qū)面積減少了12.73%，南部林地生態(tài)系統(tǒng)呈現(xiàn)減弱趨勢(shì)；低價(jià)值—高風(fēng)險(xiǎn)區(qū)面積減少了44.19%，中部草原生態(tài)空間雖然呈現(xiàn)退化狀態(tài)，但是生態(tài)安全風(fēng)險(xiǎn)逐步好轉(zhuǎn)。