拉薩市城關區近50年城市擴展過程對自然生境質量的綜合影響

聶 宇,楊彥敏,王一航,劉志鋒,*,何春陽,陳 歆

1 北京師范大學地表過程與資源生態國家重點實驗室人與環境系統可持續研究中心,北京 100875

2 北京師范大學地理科學學部自然資源學院,北京 100875

3 北京大學城市與環境學院,北京 100871

4 西藏自治區生態環境監測中心, 拉薩 850000

城市擴展過程是指由于土地利用方式的改變,導致空間上城市規模的擴大,使得非城市土地轉化為城市土地的過程[1]。自然生境質量是指森林、草地和濕地等自然生境能夠提供適合于物種與種群的資源與環境的能力,對于維持和提高區域生物多樣性至關重要[2]。伴隨著全球正在經歷的快速城市化進程,城市擴展過程已經成為影響自然生境質量的重要因素[3—4]。一方面,城市擴展過程會通過改變當地的土地利用/覆蓋情況直接導致自然生境的喪失和破碎化[5]。另一方面,城市擴展過程還會通過加劇環境污染與人類活動干擾、增加自然資源占用量和促進外來物種入侵間接造成自然生境質量下降[6—7]。近年來,已有研究在全球、國家和局地尺度上評價了城市擴展過程對自然生境質量的影響[8]。這些研究均表明,雖然在全球和國家尺度上城市擴展過程對自然生境質量的影響還較小,但在部分快速城市化熱點地區和生物多樣性保護熱點地區,城市擴展過程對自然生境質量的影響已不可忽視[3,5,9]。因此,深入理解這些熱點地區城市擴展過程對自然生境質量影響的格局、過程和機制對于當地的生物多樣性保護和可持續發展具有重要意義。

拉薩市城關區不僅是青藏高原主要的城市化地區,而且是至關重要的高寒濕地保護地[10—12]。作為西藏自治區行政中心,拉薩市城關區在西部大開發、新型城鎮化和“一帶一路”倡議的帶動下,正在經歷快速的城市化進程[13—14]。城鎮人口從1980年的38.63萬增長到了2018年的55.44萬,增長了近50%；城鎮建設用地面積也隨之快速增加[15]。該區域內的拉魯濕地國家級自然保護區是青藏高原重要的高寒濕地生態系統保護區,蘊藏著大量珍貴的野生動植物資源[16]。隨著城鎮建設用地的不斷擴張,拉魯濕地面積不斷萎縮,已經從歷史上超過1000 hm2下降到2005年設立國家級自然保護區時的核心面積660 hm2[17]。目前,區域內建設用地已經緊鄰保護區,致使該地區成為青藏高原城市擴展與生物多樣性保護之間矛盾最為突出的地區。量化該地區城市擴展過程對自然生境質量的綜合影響是優化城市擴展格局和生物多樣性保護策略的基礎,有助于實現當地乃至整個青藏高原生態文明建設。

目前,已有學者在拉薩市城關區開展了城市擴展過程對自然生境影響評估的相關研究工作。例如,張鐿鋰等[12]基于地形圖與多源遙感數據獲取了該地區1951—1999年土地利用信息,分析了城市擴展過程對自然生境的占用情況。次仁等[18]分析了距城市距離與拉魯濕地的鳥類生物多樣性之間的關系。稅燕萍等[19]評估了該地區1995—2015年城鎮建設用地增加對自然生境質量的影響。已有研究為認識該地區城市擴展過程對自然生境的影響提供了良好的基礎。但是,現有研究在獲取長時間序列的城市擴展過程信息以及綜合評估城市擴展過程對自然生境的綜合影響等方面還有待進一步完善。

長時間序列遙感數據和生境模型的結合應用,為全面評估該地區城市擴展過程對自然生境的影響提供了新途徑。鎖眼衛星提供1960—1975年期間全色波段、空間分辨率為1.8 m至12.2 m的高分辨率遙感影像數據,能反映清晰的地物光譜信息、紋理特征和空間分布情況[20]。在1995年免費向公眾開放后,該套數據已經在20世紀60年代濕地的空間恢復、土地利用/覆蓋信息提取等方面得到了較好的應用[21—23]。20世紀70年代至今的Landsat系列遙感影像和近期具有米級空間分辨率的Google Earth遙感影像也可用于準確識別城市與自然生境范圍[24—25]。此外,InVEST模型的生境質量模塊也被廣泛用于城市擴展過程對自然生境質量綜合影響的評價研究中[9,26]。因此結合鎖眼衛星遙感數據、Landsat系列遙感數據、Google Earth遙感數據與InVEST模型能實現近50年城市擴展過程對自然生境質量的綜合影響評估。

本研究目的是揭示拉薩市城關區近50年城市擴展過程對自然生境質量的綜合影響。擬回答以下3個科學問題:(1)該地區近50年城市擴展過程對自然生境質量影響的時空動態是怎樣的？(2)不同自然生境受到的影響有何差異？(3)該地區城市擴展過程導致的自然生境質量下降會對哪些物種的生存構成威脅？為了回答上述問題,首先結合鎖眼衛星、Landsat和Google Earth遙感影像量化區域土地利用/覆蓋變化情況。其次,基于InVEST模型中的生境質量模塊模擬量化城市擴展過程對自然生境質量的直接和間接影響。最后,進一步分析了自然生境質量下降對區域生物多樣性帶來的威脅。

1 研究區和數據

1.1 研究區

本文研究區為拉薩市城關區(東經91°06′至91°56′,北緯29°24′至29°36′),位于青藏高原腹地拉薩河下游,總面積520.51 km2。當地平均海拔3650 m,是世界上海拔最高的城市之一(圖1)。區域地勢北南高、中間低,南北為山地；中部為雅魯藏布江支流拉薩河中游河谷平原,地勢平坦。該地區為高原溫帶半干旱季風氣候,年平均氣溫7.8 ℃,年均降水量200—510 mm,全年日照時長3000小時以上,被稱為“日光城”[12]。拉魯濕地國家級自然保護區位于研究區西北部,總面積1220 hm2,核心面積660 hm2,是當地重要的自然生境[17]。1980—2018年,拉薩市城關區總人口從38.63萬增長到55.44萬,增長了43.52%[15],是青藏高原地區主要的人口聚集區。

圖1 研究區

1.2 數據

研究所用數據主要包括遙感影像數據、物種數據以及地理信息輔助數據。遙感影像數據包括鎖眼衛星數據、Landsat系列遙感數據、Google Earth遙感數據。其中,鎖眼衛星數據來自北京攬宇方圓信息技術有限公司,空間分辨率為2.7 m,已進行幾何校正和正射校正,用于提取拉薩市城關區1968年土地利用信息。Landsat系列遙感數據來自于中國科學院計算機網絡信息中心地理空間數據云平臺(http://www.gscloud.cn),空間分辨率為30 m,用于提取區域1984、1990、2000、2010年土地利用信息。Google Earth遙感數據來自于谷歌地球(https://www.google.com/earth/),空間分辨率為1.91 m,用于提取2019年土地利用信息。此外,物種分布數據來自國際自然保護聯盟瀕危物種紅色名錄(IUCN Redlist)(https://www.iucnredlist.org/)。研究區的行政邊界、行政中心等基礎地理信息輔助數據來源于國家基礎地理信息中心(http://ngcc.sbsm.gov.cn/)。為了保證數據一致性,所有柵格數據采用Albers等積圓錐投影并重采樣至30 m分辨率。

2 方法

2.1 量化區域自然生境質量空間格局

首先,通過目視解譯,基于1968—2019年間的6期(1968、1984、1990、2000、2010、2019年)遙感數據獲取城鎮建設用地、農村建設用地、耕地、沼澤、林地、河漫灘、草地和水體8類土地利用/覆蓋信息(圖3)。解譯標準見圖2。

圖2 三類遙感數據的判讀標準

通過實地考察和混淆矩陣分析,對目視解譯提取的土地利用/覆蓋數據進行精度驗證。首先于2020年夏季進行了實地考察,分別在城鎮建設用地、林地、濕地(包括沼澤、河漫灘和水體)和草地等土地利用類型中選取了118個、40個、60個、43個樣本點,獲取了真實的土地利用/覆蓋信息。其次,在實地考察的基礎上對分類結果進行精度驗證。結果表明,全區2019年的數據分類總體精度為92.80%,Kappa系數為0.90(圖3)。其中,城鎮建設用地生產者精度為93.22%,用戶精度為95.65%；林地的生產者精度為90.00%,用戶精度為76.60%；濕地生產者精度為97.67%,用戶精度為97.67%；草地生產者精度為90.48%,用戶精度為98.28%。由此可見,獲取的拉薩市城關區土地利用/覆蓋數據準確性較高。

圖3 流程圖

圖4 土地利用/覆蓋分類結果的精度驗證

然后,參考國際自然保護聯盟發布的生境分類標準[27],選取林地、草地、沼澤、河漫灘、水體作為自然生境,城鎮建設用地作為威脅源[9]。采用InVEST模型的生境質量模塊量化1968—2019年區域自然生境質量空間格局動態。自然生境質量的計算公式為:

(1)

式中,Qxj為土地利用/覆蓋類型j中像元x上的生境質量；Hj為土地利用/覆蓋類型j的生境適宜性；Dxj為土地利用/覆蓋類型j中像元x受到總的脅迫水平；k為半飽和常數,通常取Dxj最大值的一半；z為歸一化常量。每個像元上生境受到總的脅迫水平計算公式為:

(2)

式中,R為威脅因子總數,wr為受威脅影響的權重,Yr為柵格總數,βx為可達性,Sjr為敏感性,irxy為像元y中威脅因子ry對像元x的威脅類型。城鎮的威脅類型為指數型,計算公式為:

(3)

式中,dxy為像元x與y的距離,drmax為威脅因子的最大影響距離。生境適宜性與威脅因子參數均參考相關研究[9,19,28]設置(表1)。

表1 InVEST模型生境質量模塊的參數設置

最后,參考Song等[9]和Moreria等[29]的研究,將自然生境質量進行分級評價。根據研究區實際情況,將自然生境質量由低到高劃分為0—0.3、0.3—0.6、0.6—0.8和0.8—1等4個等級,分別表示自然生境質量的低、中、較高和高水平。

2.2 評估區域城市擴展過程對自然生境質量的影響

通過量化城市擴展過程對自然生境質量的總影響、城市擴展過程對自然生境質量的直接與間接影響等三個方面,綜合評估拉薩市城關區城市擴展過程對區域自然生境質量的影響。城市擴展過程對自然生境質量的總影響是指城市擴展過程導致的自然生境質量變化率[30]。自然生境質量變化率計算公式為:

(4)

式中,ΔQi-1968表示區域1968年到第i年城市擴展過程導致的自然生境質量變化率,Q1968和Qi分別表示區域1968年和第i年的自然生境質量。根據Seto等[4]和McDonald等[31]開展的相關研究,城市擴展過程對自然生境質量的直接影響是指新增城市土地直接占用導致的生境質量變化情況；間接影響是指在城市最大影響范圍內未被直接占用的自然生境質量變化情況。

2.3 評估區域生境質量變化對生物多樣性的影響

通過統計區域內城鎮建設用地最大影響范圍內物種的分布情況評估拉薩市城關區生境質量變化對生物多樣性的影響。物種分布數據來自國際自然保護聯盟瀕危物種紅色名錄,其給出了全球12萬余種物種的潛在空間分布范圍、適宜生境和瀕危程度等信息。參考He等[3]的研究,通過空間疊加的方式篩選出拉薩市城關區城鎮建設用地最大影響范圍內可能存在的物種,統計了物種名稱、物種類型、受威脅情況、適宜的生境類型和是否為本地種等關鍵信息,進而量化了各物種潛在分布范圍內自然生境質量受城市擴展過程的影響情況。以上方法通過IDL語言編程、Python語言編程和ArcGIS 10.2軟件操作實現。

3 結果

3.1 區域1968年自然生境質量空間格局

拉薩市城關區自然生境面積大,生境質量整體水平中等。1968年,拉薩市城關區自然生境面積為476.58 km2,占區域總面積的91.57%；自然生境質量為0.60,屬于中等水平(圖5)。其中,高質量自然生境面積為20.12 km2,占全區自然生境總面積的4.22%,主要分布在區域城鎮建設用地的南北兩側,多為河流、沼澤和林地。較高質量的自然生境面積為305.31 km2,占全區自然生境總面積的64.10%,主要分布在遠離城市區域的南北兩側山地和拉薩河兩岸,多為草地與河漫灘。中等質量的自然生境面積為151.18 km2,占全區自然生境總面積的31.71%,分布在靠近城市的城關區南北兩側,以草地為主。

圖5 1968年自然生境質量空間格局

3.2 區域1968—2019年城市擴展過程

拉薩市城關區城鎮建設用地呈指數型增長(圖6)。1968—2019年,全區城鎮建設用地總面積從9.78 km2增長到的67.67 km2,增長了6.92倍,年均增長1.14 km2(圖7)。在1968—1984年期間,城市擴展過程相對緩慢。全區城鎮建設用地面積增加了10.58 km2,增長率為108.18%,年均增長0.66 km2。而在2010—2019年,該區城鎮建設用地迅速增加,增加了23.18 km2,增長率為52.10%,年均增長2.58 km2。相比1968—1984年,2010—2019年的城鎮建設用地增長速度增加了3.91倍。從空間上看,拉薩市城關區的城鎮建設用地主要在拉薩河北岸呈“倒L形”分布。1968—1984年城鎮建設用地主要向西、北方向擴展,1984—2000年城鎮建設用地轉為向西擴展,2000—2019年城鎮建設用地沿著拉薩河向東、北方向擴展。

圖6 1968—2019年土地利用/覆蓋變化過程

圖7 1968—2019年各地類面積的變化過程

3.3 區域1968—2019年城市擴展過程對自然生境的占用

拉薩市城關區城市擴展過程加速占用了區域自然生境。1968—2019年,城市擴展過程占用了39.57 km2的自然生境,占總自然生境面積的8.30%(表2)。在1968—1984年間由于城市擴展,城關區自然生境總面積減少了8.71 km2,年均下降0.55 km2。2010—2019年間的總面積減少了13.85 km2,年均下降1.54 km2,面積減少速度增加了1.8倍。從不同生境來看,草地、濕地與林地的總面積分別減少了29.69 km2、7.95 km2和1.93 km2(表2)。其中林地與草地的面積加速減少,分別從年均減少0.04 km2和0.39 km2增加到年均減少0.23 km2和1.37 km2。濕地的面積減速下降,從年均減少0.12 km2下降到年均減少0.03 km2(圖8、表2—7)。

圖8 1968—2019年自然生境年均面積損失情況

表2 1968—2019年土地利用/覆蓋轉移矩陣/km2

表3 1968—1984年土地利用/覆蓋轉移矩陣/km2

表4 1984—1990年土地利用/覆蓋轉移矩陣/km2

表5 1990—2000年土地利用/覆蓋轉移矩陣/km2

表6 2000—2010年土地利用/覆蓋轉移矩陣/km2

表7 2010—2019年土地利用/覆蓋轉移矩陣/km2

3.4 區域1968—2019年城市擴展過程對自然生境質量的影響

拉薩市城關區城市擴展過程導致區域自然生境質量加速下降。1968—2019年,城市擴展過程導致區域自然生境質量從0.57下降到0.50,下降了12.69%。在1968—1984年,自然生境質量下降了2.87%,而2010—2019年自然生境質量下降了4.22%,自然生境質量下降速度增加了1.5倍(圖9)。從空間上看,城鎮建設用地周邊地區自然生境質量下降明顯,自然生境質量變化率高達80%以上(圖10)。在不同自然生境類型中,林地和濕地受到城市擴展過程的影響較大。其中,林地生境質量從1968年的0.80下降到2019年的0.03,下降了96.59%。濕地的自然生境質量從0.88下降到0.45,下降了49.21%。(圖11、表8)。

圖9 1968—2019年自然生境質量變化趨勢

圖10 1968—2019年自然生境質量變化的空間格局

圖11 1968—2019年主要自然生境類型質量的變化趨勢

城市擴展過程對不同自然生境類型的直接和間接影響存在差異。1968—2019年,拉薩市城關區城市擴展過程直接導致全區自然生境質量下降7.06%,間接導致自然生境質量下降5.63%(圖12)。在不同自然生境類型中,林地、草地主要受到城市擴展過程的直接影響,濕地主要受到城市擴展過程的間接影響。其中,林地受到城市擴展過程的直接影響導致其生境質量下降80.73%,間接影響為15.86%,直接影響是間接影響的4.8倍；草地受到的直接影響使其生境質量下降6.64%,間接影響為1.27%,直接影響是間接影響的5.2倍(圖13,14,表8)。而城市擴展過程對濕地的直接影響導致其生境質量下降22.70%,間接影響使其生境質量下降26.51%,間接影響是直接影響的1.2倍(圖15,表8)。

表8 城市擴展過程對自然生境質量的影響/%

圖12 城市擴展過程對自然生境質量的直接影響與間接影響

圖13 城市擴展過程對林地生境質量的直接影響與間接影響

圖14 城市擴展過程對草地生境質量的直接影響與間接影響

圖15 城市擴展過程對濕地生境質量的直接影響與間接影響

3.5 區域1968—2019年自然生境質量下降對生物多樣性的影響

拉薩市城關區物種豐富。根據IUCN Redlist數據對區域物種分布進行統計后發現,區域內共包括鳥類、哺乳類和兩棲類等12類、161種物種(圖16),其中有鳥類生物48種,哺乳類生物28種,分別占總物種量的29.8%和16.8%。161種物種中有藏羚羊(Procaprapicticaudata)、馬麝(Moschuschrysogaster)、小黃黝魚(Micropercopsswinhonis)等7種本地種,以及青頭潛鴨(Aythyabaeri)、雪豹(Pantherauncia)、獵隼(Falcocherrug)等16種瀕危物種(表9)。拉薩市城關區自然生境質量下降對161種物種均有影響。值得注意的是,全區16種瀕危物種的適宜生境均為生境質量快速下降的林地、濕地或草地,且其中瀕危程度為易危及以上的物種有6種,占瀕危物種總量的37.50%(表9)。

圖16 拉薩市城關區生物多樣性狀況

表9 拉薩市城關區受威脅物種

4 討論

4.1 長時間序列遙感數據與InVEST模型為評估城市擴展過程對自然生境質量的綜合影響提供了可靠途徑

現有用于評估城市擴展過程對自然生境質量影響的數據主要為實地調查數據與短時間序列遙感數據[32—34]。其中,實地調查數據需要依賴可靠的數據獲取方式和獲取地點。因此,實地調查數據具有高花費、短時間、小范圍等缺陷。更重要的是,青藏高原地區缺乏大量的實地調查數據。短時間序列遙感數據只能反映近期的城市擴展過程,無法獲得長時間序列城市擴展信息與自然生境本底情況[8]。

結合鎖眼衛星數據、Landsat遙感數據和Google Earth遙感數據能夠彌補上述數據的不足。具有2.7 m空間分辨率的鎖眼衛星高分遙感影像可以清晰展現地物的光譜、紋理特征和空間分布情況,因此可以使用該數據準確識別區域的自然生境和城市建成區范圍[20,35]。例如,Brinkmann等[21]使用鎖眼衛星數據恢復了1960年西非四個城市地區的土地利用情況。Hepcan等[36]通過鎖眼衛星數據與航片恢復了1963年土耳其伊茲密爾地區的土地利用情況。劉志鋒等[23]基于鎖眼衛星數據提取了1964年圖們江下游地區濕地分布情況。Landsat系列遙感影像和具有米級空間分辨率的Google Earth遙感影像也可用于準確識別城市與自然生境范圍。例如,牛振國等[25]基于Landsat系列遙感影像和Google Earth遙感影像提取了2000年中國濕地的分布范圍。Gong等[24]基于Landsat系列遙感影像和Google Earth遙感影像提取出了中國1978—2017年城市不透水面分布范圍。結合以上三類遙感數據能夠準確量化長時間序列城市擴展過程和揭示歷史自然生境質量本底信息。

此外, InVEST模型生境質量模塊可以同時量化城市擴展過程對生境質量的直接影響和間接影響。模型通過結合土地覆被和生物多樣性威脅因素的信息生成生境質量地圖,土地利用/覆被變化會導致自然生境的減少和鄰近生境的退化[37],其中,自然生境減少為直接影響,鄰近生境退化的邊緣效應即為間接影響。因此,InVEST模型生境質量模塊可以通過評估生境適宜性變化和生境受威脅程度變化同時量化直接和間接影響[37]。例如,韓艷莉等[28]基于InVEST模型模擬結果發現2005—2015年青海湖流域內的建成區和交通用地面積增加直接占用了自然生境,最終導致自然生境質量不斷下降。稅燕萍等[19]以拉薩河流域為例,使用InVEST模型量化了城市擴展過程對自然生境的間接影響,研究發現,1995—2015年間拉薩河流域內社會經濟發展間接影響導致當地自然生境質量下降。

綜上所述,本文基于鎖眼衛星影像、Landsat系列遙感影像和Google Earth遙感影像長時間序列遙感數據,結合InVEST模型的生境質量模塊同時量化直接影響和間接影響,能夠全面地揭示拉薩市城關區長時間尺度的城市擴展過程對區域自然生境質量的影響。

4.2 拉薩市城關區城市擴展過程導致的自然生境質量下降已對生物多樣性造成威脅

拉薩市城關區的自然生境對于生物多樣性保護、當地生態系統的維持等具有重要作用。緊鄰拉薩城鎮建設用地的拉魯濕地國家級自然保護區是青藏高原十分重要的高寒濕地保護區。作為雅魯藏布江最大的支流,拉薩河有25 km長的河道在研究區內橫穿而過。研究表明,拉魯濕地內有針尾鴨(Anasacuta)、斑嘴鴨(Anaszonorhyncha)、斑頭雁(Anserindicus)等鳥類62種,拉薩河內有包括拉薩裂腹魚(Racomawaltoni)、尖裸鯉(Oxygymnocyprisstewarti)、雙須葉須魚(Ptychobarbusdipogon)等魚類24種[16,38]。

城市擴展過程導致的自然生境質量下降已經對當地物種生存造成了威脅。區域內鳥類和哺乳類生物主要的棲息地為濕地,而城市擴展過程對區域濕地造成了嚴重的影響。本文發現,1968—2019年,拉魯濕地生境質量從0.90下降到0.40,下降了55.42%。已有學者在拉魯濕地開展了相關研究,也得到了類似的發現。例如,劉曉曼等[17]研究表明發現,由于人類活動加劇等的影響,拉魯濕地面積在1989—2009年期間從942.45 hm2下降到685.33 hm2,減少了27.30%。范麗卿等[39]相關研究表明濕地內分布有5種土著魚類和7種外來魚類,隨著拉薩城市化過程與人類活動加劇,拉魯濕地內的土著魚類數量幾近滅絕,外來魚類成為優勢種。田瀚鑫等[40]對拉魯濕地水質開展的研究發現,區內枯水期水質主要為V類,豐水期水質主要為I和III類,水質污染區位于區內東北部和中南部,靠近城市道路,受到人類活動影響極大。周生靈等[41]研究發現拉魯濕地猛禽的觀察量減少,主要原因在于城市擴展過程導致的濕地周邊適宜猛禽營巢的生境持續減少。本研究與相關研究結果均表明拉魯濕地生境質量的下降對拉魯濕地內生物多樣性造成了嚴重的威脅。

4.3 政策啟示

拉薩市城關區的城市擴展過程對自然生境質量的影響具有三個特征:快速的城市擴展過程導致區域自然生境質量加速下降；在不同自然生境類型中,林地和濕地的生境質量下降幅度大；城市擴展過程對不同自然生境類型的直接和間接影響存在差異,林地主要受到直接影響,濕地主要受到間接影響。

本研究在拉薩市城關區發現的現象普遍存在于青藏高原地區。青藏高原自然生境本底相對脆弱,對城市擴展過程等人類活動具有高度敏感性[42—46]。Pan等[44]的研究發現,青藏高原地區的城市正在經歷加速的城市擴展過程。Wang等[10]研究結果表明,西寧、拉薩、香格里拉和格爾木等地是青藏高原地區快速城市擴展過程的熱點地區,這些地區大多分布于青藏高原上的河谷或綠洲地區。同時,這些地區也是濕地和森林廣泛分布、生物多樣性較為豐富的區域。因此,青藏高原上的城鎮發展勢必會對濕地和森林生態系統構成嚴重威脅(圖16)。

圖17 青藏高原主要城市的城市擴展過程與自然生境的空間關系

因此,亟待采取有效措施減少城市土地對自然生境的直接占用和間接影響。首先,加強土地管控,控制城市蔓延。其次,對拉魯濕地、拉薩河等高質量的自然生境地區劃定嚴格的保護紅線,嚴守生態保護紅線,減少城市土地對自然生境的直接占用[47—48]。最后,在保護區等重要自然生境周圍建立緩沖區,緩沖區內執行有效的管控措施,如禁止鳴笛,控制道路上的車流量與人流量,夜間減暗或關閉道路路燈等,以盡量減少城市對于自然生境的間接影響,從而緩解由于城市擴展過程導致的自然生境質量下降。

4.4 未來展望

拉薩市城關區位于青藏高原腹地,生態環境本底脆弱,快速的城市擴展過程對當地自然生境帶來嚴重威脅。本文結合長時間序列遙感數據與生境質量模型,全面評估了拉薩市城關區近50年的城市擴展過程對自然生境質量的綜合影響,對保護區域生物多樣性、促進城市景觀可持續、緩解人地矛盾具有參考價值。

研究仍存在不足。一方面,在生境質量模型參數設置方面,生境質量賦值與威脅因子權重等均參考現有研究確定,模型參數本地化還需要結合實地考察數據進一步完善。另一方面,沒有考慮自然生境斑塊破碎化等進一步導致的連通性下降等問題,導致評價結果可能存在低估。另外,在拉薩市城市擴展過程中,當地有關部門采取了植樹造林和濕地恢復等一系列措施來提高自然生境質量,在文中還未體現。

未來研究將通過多尺度融合與多技術結合,進一步聚焦青藏高原城市擴展過程對自然生境的綜合影響評估。通過實地調查、遙感技術、空間統計和模型模擬等多方法結合,結合青藏高原區域特征進行模型本地化處理,預估未來城市擴展過程以及城市周邊生態保護與恢復措施對區域自然生境的可能影響,為青藏高原城市景觀可持續與生態安全屏障建設提供科學建議。

5 結論

本文結合鎖眼衛星、Landsat和Google Earth在內的長時間序列遙感數據,應用InVEST模型生境質量模塊,有效地評估了拉薩市城關區近50年城市擴展過程對自然生境質量的綜合影響。研究區1968—2019年城市擴展過程導致區域自然生境質量加速下降,共下降了12.69%。自然生境質量同時受到城市擴展過程的直接影響和間接影響。在不同自然生境類型中,林地和濕地受到城市擴展過程影響明顯,生境質量分別下降了96.59%和49.21%。其中,林地主要被城鎮建設用地直接占用,濕地則主要受城市擴展過程的間接影響。該地區城市擴展過程導致的自然生境質量下降已經對16種瀕危物種構成了威脅。

因此,在青藏高原的城市擴展過程中,需同時關注新增城鎮建設用地對林地和濕地的直接占用和間接影響。建議在自然保護區等重要自然生境周邊建立緩沖區,在控制緩沖區內城市擴張和人類活動強度的同時積極推行生態恢復措施,協調城市建設與生物多樣性保護,以促進區域可持續發展。