李昌玉，王鏡，李姝彬，李京梅

(西寧市氣象局，青海 西寧 810016)

引言

從20世紀80年代開始，國內(nèi)外學者從不同區(qū)域、不同方法、不同評價體系對生態(tài)脆弱性進行了大量的研究，但針對西寧地區(qū)生態(tài)環(huán)境脆弱性的評價研究報道尚未發(fā)現(xiàn)，在查閱相關文獻資料的基礎上，借鑒前人的研究方法對西寧地區(qū)的生態(tài)環(huán)境脆弱性進行評價研究。

青藏高原生態(tài)系統(tǒng)是中國重要的碳素匯集生態(tài)系統(tǒng)，具有維系國家生態(tài)安全的重要地位和作用。是全國生態(tài)環(huán)境變化的“脆弱帶”[1]。西寧地處青藏高原的東北部，也是生態(tài)環(huán)境變化脆弱區(qū)。西寧地區(qū)經(jīng)過發(fā)展，已成為青藏高原人口超百萬的區(qū)域性現(xiàn)代化城市[2]，隨著城市規(guī)模的不斷擴大，城鎮(zhèn)人口的快速發(fā)展，人類活動的不斷增大，使得生態(tài)環(huán)境壓力越加緊張，生態(tài)環(huán)境趨向于越來越脆弱。因此西寧地區(qū)進行生態(tài)環(huán)境脆弱性評價項目的確定，可以認識西寧地區(qū)生態(tài)環(huán)境的脆弱狀況，并根據(jù)西寧地區(qū)生態(tài)環(huán)境評價結果，對不同脆弱區(qū)提出相應的應對措施和保護手段，為實現(xiàn)資源和經(jīng)濟社會的可持續(xù)發(fā)展提供技術支撐，促進西寧地區(qū)穩(wěn)定的生態(tài)系統(tǒng)，尋找出在生態(tài)脆弱欠發(fā)達地區(qū)走一條綠色發(fā)展的新路。

1 研究區(qū)概況

西寧地區(qū)位于青藏高原東北部湟水谷地，地理位置100°52′04″～101°54′54″E，36°13′42″～37°28′11″N[3]，面積7 665km2[4]，大體上由達坂山—湟水谷地—拉脊山組成，地勢西北高、東南低，最高海拔4877m，最低海拔2162m位于湟水河谷地，全市海拔高度差最大為2715m。水熱條件分布不均勻，經(jīng)西寧市氣象局統(tǒng)計1990—2020年30年氣象要素數(shù)據(jù)，年平均氣溫4.5℃，年降水量470.6mm，年平均風速1.3m/s，年日照時數(shù)2 582.6h，年平均氣壓751.0hPa。與同緯度地區(qū)相比，西寧地區(qū)具有氣壓和氣溫較低、晝夜溫差大、日照時數(shù)長、蒸發(fā)量大的氣候特點，土壤貧瘠、植被稀疏、地質結構復雜，兼有農(nóng)業(yè)和少部分牧業(yè)，生態(tài)環(huán)境先天脆弱。2019年西寧地區(qū)總人口238.71萬人，是青藏高原唯一過百萬人口的城市[5]，是黃河流域最上游和青藏高原上人類活動強度最大的地區(qū)[6]。

2 資料和方法

2.1 資料來源

使用的氣象數(shù)據(jù)2020年均氣溫和年降水量，來自青海省、市(州)、縣(區(qū))氣象臺業(yè)務平臺，選取西寧地區(qū)國家氣象站和區(qū)域氣象自動站共43個站點；人口密度和GDP等社會統(tǒng)計數(shù)據(jù)來源于《2020年西寧地區(qū)統(tǒng)計年鑒》；地形因素數(shù)據(jù)來自SRTM DEM，分辨率為90m；土地資源數(shù)據(jù)來自中科院資源環(huán)境數(shù)據(jù)中心2020年度植被指數(shù)數(shù)據(jù)集，分辨率為30m。

2.2 研究方法

生態(tài)脆弱性是生態(tài)系統(tǒng)在特定時空尺度相對于外界干擾所具有的敏感反應和自我恢復能力，是生態(tài)系統(tǒng)的固有屬性[7]。是指在外界干擾下生態(tài)系統(tǒng)產(chǎn)生的一系列反應能力。生態(tài)脆弱性指標選擇需要考慮研究區(qū)的抗干擾能力、氣候變化敏感性、時空波動性等多方面特征[8]。本項目在西寧地區(qū)生態(tài)環(huán)境最基本的特征基礎上，參照國內(nèi)研究人員已研究的脆弱性結果[9-11]構建指標，以及對現(xiàn)有研究數(shù)據(jù)的可獲得性和表向性的考慮，西寧地區(qū)生態(tài)環(huán)境脆弱性評價分析從地形因素、氣象因素、社會經(jīng)濟(人類活動)、土地資源四個方面選取植被覆蓋度、土地利用、土壤類型、高程、坡度、坡向、年平均氣溫、平均降水量、GDP、人口、工業(yè)總產(chǎn)值11個指標來進行研究。為得到西寧區(qū)域內(nèi)生態(tài)環(huán)境的特征細節(jié)，應用規(guī)則網(wǎng)格結合縣級行政區(qū)劃單元用海拔高程和坡度來表現(xiàn)西寧地形地貌脆弱情況；年溫度和年降水量來代表西寧氣候狀況；生存環(huán)境的脆弱度用土地利用狀況來表現(xiàn)；人口和經(jīng)濟壓力來代表西寧經(jīng)濟社會因素引起的脆弱度。按照以往研究人員的經(jīng)驗，各指標正負作用的性質，如社會經(jīng)濟越好，表示人類對生態(tài)環(huán)境干擾越大，生態(tài)脆弱度就越大，即呈正相關；植被覆蓋程度越好，生態(tài)脆弱程度越小[12]，即呈負相關。表1為西寧地區(qū)各因素評價指標和權重。

表1 西寧地區(qū)生態(tài)環(huán)境脆弱性評價指標及權重

數(shù)據(jù)標準化：在參考各位學者研究生態(tài)脆弱性評價研究的基礎上，先對各指標數(shù)據(jù)分別用公式(1)、(2)進行標準化處理[13]。

(1)

(2)

計算各指標權重，此研究用層次分析法來確定。用各因素標準化值和該指標權重的乘積來計算各要素的脆弱性指數(shù)。最后計算綜合生態(tài)環(huán)境脆弱性指數(shù)。

3 特征和分析

3.1 地形因素脆弱性分析

地形因素是生態(tài)脆弱性評價中不可缺少的一部分，是生態(tài)環(huán)境的基本要素，地形因素生態(tài)脆弱性與海拔高度和坡度成正相關，海拔越高，坡度越大，生態(tài)脆弱性越大。從地形因素生態(tài)脆弱性空間分布來看(見圖1)，中度和重度脆弱占的面積較小兩者之和為19.66%，主要分布高海拔山區(qū)，即在大通北面達坂山、大通西面娘娘山、湟源西南面的日月山、湟中西南面的拉脊山等，為地勢起伏大、海拔高度高、山脈縱橫處。輕度脆弱區(qū)隨中度脆弱區(qū)依附在一起隨山體分布，面積占比為16.31%。潛在脆弱區(qū)面積占比最多為42.73%，主要分布在海拔相對較低的地帶，如大通東南部、西寧地區(qū)、湟中大部和湟源局地。

圖1 西寧地區(qū)地形系統(tǒng)生態(tài)脆弱性評價圖

3.2 氣象因素脆弱性分析

對處于高寒冷涼區(qū)的西寧地區(qū)而言，降水較少，蒸發(fā)量大，氣溫日較差大，冬季寒冷干燥，氣候對植被生長、土壤環(huán)境有嚴重的影響。因此，水熱條件越好生態(tài)脆弱性越小，水熱條件越差生態(tài)脆弱性越大。從表2可見，氣象因素脆弱性面積占比由大到小依次為中度占36.35%、重度占28.05%、輕度占26.75%、微度占8.34%、潛在脆弱占0.51%。從氣象因子的脆弱性空間分級圖(見圖2)可看出，西寧地區(qū)氣象因素的生態(tài)脆弱性，重度脆弱主要集中在湟源西部和西南部、大通西北部和湟中西北小部，這些地區(qū)相對而言降水較多但溫度較低。處于河谷的湟源縣城、大通縣城及東南部、湟中大部、西寧城區(qū)因降水相對較多，溫度較高，生態(tài)脆弱性逐步降小。

圖2 西寧地區(qū)氣候生態(tài)脆弱性評價圖

3.3 社會經(jīng)濟脆弱性分析

人口密度、人均GDP和工業(yè)總產(chǎn)值作為影響西寧地區(qū)社會經(jīng)濟因素的生態(tài)脆弱性，與生態(tài)脆弱性成正相關，人口多、GDP高和工業(yè)總產(chǎn)值高的地方，生態(tài)脆弱性越大。從社會經(jīng)濟因素評價圖(圖3)可見，西寧地區(qū)社會經(jīng)濟因素的生態(tài)脆弱性分布大小由中心城區(qū)向周邊擴展，重度脆弱集中在西寧城北區(qū)、城西區(qū)和湟中區(qū)，因這些區(qū)域人口密集、經(jīng)濟較發(fā)達給環(huán)境帶來的壓力相對較高故生態(tài)脆弱性較大；相反，研究區(qū)西北部和西南部社會經(jīng)濟因素影響較小的地區(qū)生態(tài)脆弱性較小。

圖3 西寧地區(qū)社會經(jīng)濟生態(tài)脆弱性評價圖

3.4 土地資源脆弱性分析

影響西寧地區(qū)土地資源生態(tài)脆弱性主要有植被覆蓋度和土地利用類型，其中植被覆蓋度與脆弱性成負相關，植被覆蓋度越大脆弱性越小。從西寧地區(qū)土地資源生態(tài)脆弱性空間分布圖(見圖4)可看出，中度脆弱、重度脆弱和輕度脆弱主要分布在湟水河谷地城區(qū)人造地表、裸地以及高海拔山脈植被少的地方，三者面積之和占13.2%，剩余大部分地區(qū)為潛在脆弱和微度脆弱區(qū)，兩者面積之和達87.14%，占比較大。

圖4 西寧地區(qū)土地資源生態(tài)脆弱性評價圖

3.5 生態(tài)脆弱性綜合評價分析

在計算分析了地形、社會經(jīng)濟、水文氣象和土地資源四個因素的生態(tài)脆弱性之后，用(3)式來計算綜合生態(tài)環(huán)境脆弱性[14]。

EVI=∑[A1Y1+A2Y2+A3Y3+A4YA]

(3)

式中，A1,A2,A3,A4,分別表示地形、社會經(jīng)濟、水文氣象和土地資源因素生態(tài)脆弱度，Y1,Y2,Y3,Y4分別為各因素對應的權重，此處權重用專家打分法確定。

西寧地區(qū)生態(tài)脆弱度使用Arcgis中的“自然間斷點”分為五個等級見表2。繪制出西寧地區(qū)生態(tài)環(huán)境脆弱性等級面積占比圖(見圖5)和西寧地區(qū)綜合生態(tài)脆弱性評價圖(見圖6)。由圖5可知，西寧地區(qū)生態(tài)脆弱性一般潛在區(qū)域面積占17.59%，微度脆弱占43.61%，輕度脆弱占28.86%，中度脆弱占8.0%，重度脆弱占1.94%，由此可見西寧地區(qū)微度脆弱面積占比最大。從圖6上分析，西寧地區(qū)生態(tài)脆弱度空間分布總體呈現(xiàn)高海拔向低海拔減小趨勢。重度和中度主要分布西寧地區(qū)北部和南部高海拔山脈區(qū)，這些地區(qū)因溫度低，以裸巖石質地景觀為主的高山，植被覆蓋度極低。輕度脆弱區(qū)與中度脆弱區(qū)依附在一起分布在高海拔山區(qū)，另外一塊輕度脆弱區(qū)聚集于主城區(qū)的城北區(qū)、城西區(qū)和城中區(qū)，其中城西區(qū)和城中區(qū)以微度脆弱之中鑲嵌輕度脆弱，而城北區(qū)是以微度和輕度脆弱混合并有一小塊中度脆弱。主要因城區(qū)人口密度大，建設開發(fā)用地增加，工業(yè)總產(chǎn)值高，尤其城北區(qū)工業(yè)總產(chǎn)值最高，人均GDP較高，故脆弱性較高。城東區(qū)因海拔最低，氣溫最高，加之工業(yè)總產(chǎn)值在城區(qū)中是最低的，故為微度脆弱區(qū)。大通、湟源和湟中海拔較低的谷地和緩坡地區(qū)為微度脆弱區(qū)，相對研究區(qū)而言微度脆弱區(qū)域分布范圍最廣，這些區(qū)域基本上坡度較緩，是以耕地、林地和草地為主的地區(qū)，植被好，人類活動少。

表2 西寧地區(qū)生態(tài)脆弱性評價占比(%)

圖5 西寧地區(qū)生態(tài)環(huán)境脆弱性等級面積占比圖

圖6 西寧地區(qū)綜合生態(tài)脆弱性評價圖

4 小結

(1)由西寧地區(qū)2020年生態(tài)脆弱性等級面積占比統(tǒng)計可知，微度脆弱面積占比最高達到43.61%，其次是輕度脆弱占28.86%，中度和重度脆弱占比僅為8.0%和1.94%，西寧地區(qū)整體處于微度脆弱性。

(2)生態(tài)脆弱性空間分布上整體呈現(xiàn)高海拔山區(qū)高，低海拔緩坡和河谷低的分布特征，重度和中度脆弱區(qū)主要分布在海拔較高的山區(qū)，潛在或微度脆弱區(qū)分布在海拔相對較低的谷地和緩坡，植被覆蓋好，工業(yè)總產(chǎn)值相對低的地區(qū)。大通、湟源和湟中海拔較低的谷地和緩坡地區(qū)為微度脆弱區(qū)，相對西寧地區(qū)而言微度脆弱區(qū)域分布范圍最廣。

(3)西寧城區(qū)內(nèi)人口集中，人類活動擾動大，工業(yè)總產(chǎn)值高，故生態(tài)脆弱性較高。其中城北區(qū)是以微度和輕度脆弱混合并有一小塊中度脆弱，是城區(qū)中最脆弱的區(qū)域，其次是城西區(qū)和城中區(qū)，城東區(qū)為微度脆弱區(qū)。

(4)生態(tài)環(huán)境脆弱性受自然、氣候和人類等各個方面的共同影響，相互關系錯綜復雜。限于收集數(shù)據(jù)的難度，在此次研究西寧地區(qū)脆弱性時只用了一年的資料，并且沒有滲透到更廣更細的層面選取更多的生態(tài)環(huán)境脆弱性指標，而是主要根據(jù)結構特點選取了四個有代表性因素的指標，對于本身復雜的生態(tài)環(huán)境系統(tǒng)，根據(jù)特定指標進行環(huán)境評價有一定的局限性，還有待于今后結合實證研究進一步完善。