青藏高原旅游活動強度時空特征及影響因素分析

鄧婷 ，徐佩，逯亞峰，李政旸

（1.中國科學院 水利部成都山地災害與環(huán)境研究所，四川 成都 610041；2.中國科學院大學 資源與環(huán)境學院，北京 100049）

近半個世紀以來，不斷加劇的氣候變化和人類活動嚴重威脅著全球自然生態(tài)系統[1]。青藏高原是中國乃至亞洲重要的生態(tài)安全屏障，高寒條件使得青藏高原生態(tài)系統敏感且環(huán)境極其脆弱。在全球變化和人類活動綜合影響下，青藏高原生態(tài)系統退化明顯，生態(tài)安全屏障面臨嚴峻挑戰(zhàn)[2]。同時，壯麗的自然風光和獨特的人文景觀讓青藏高原成為國內外旅游者青睞的旅游目的地。旅游業(yè)已經成為青藏高原各省（區(qū)）的支柱產業(yè)和重點產業(yè)[3]。2019年，西藏自治區(qū)全年接待游客4 012.15萬人次，旅游收入559.28億元，旅游收入占GDP比重達32.94%。旅游活動是青藏高原重要的人類活動之一，適當強度的旅游活動能優(yōu)化當地生態(tài)環(huán)境，促進文化保護傳承，帶動欠發(fā)達地區(qū)發(fā)展，但過度的旅游活動將增大資源環(huán)境壓力，降低生態(tài)系統穩(wěn)定性，威脅青藏高原生態(tài)安全[4-6]。開展青藏高原旅游活動強度評估，明晰青藏高原旅游活動強度的時空特征和影響因素，將支撐青藏高原生態(tài)保護與旅游經濟協同促進，是青藏高原實現旅游可持續(xù)發(fā)展的前提。

伴隨著全球旅游業(yè)快速發(fā)展，越來越多的學者開始關注對旅游活動強度的研究。旅游活動強度的研究從早期的定性分析轉變?yōu)楸O(jiān)測、調查、統計等定量評估。目前，監(jiān)測方面主要通過監(jiān)測客流量、車流量等手段來動態(tài)反映旅游活動強度[7]；調查方面根據調查旅游地步道附近土壤理化性質[8]、植被群落結構[9]、生物多樣性[10-12]等側面反映旅游活動強度；而統計方面則依靠旅游人口規(guī)模和旅游經濟效益相關統計數據來直接反映旅游活動強度[13-15]。

然而，上述研究主要集中在經濟發(fā)展水平較高、旅游業(yè)發(fā)展較好的地區(qū)，對于經濟基礎相對薄弱、旅游業(yè)發(fā)展迅速的高原地區(qū)關注不足。同時，現有研究的研究單元主要包括省區(qū)[16]、城市[17-18]、景區(qū)[19-21]等類型，且多針對某一具體景區(qū)開展研究，而以縣級尺度開展的研究較少。有學者以青藏高原東南緣香格里拉市積雪動態(tài)變化反映旅游活動[22]，利用騰訊位置大數據分析青藏高原旅游活動[23]。但都只強調不同旅游活動強度下單一指標的變化，對不同區(qū)域的指標選取缺乏系統性。此外，有研究指出青藏高原旅游資源破壞[22]、水土流失[24]、濕地減少[25]和草地退化[26]等問題與人類活動有關，卻少有研究直接明確青藏高原旅游活動強度的時空特征及影響因素。

基于此，本文通過構建青藏高原旅游活動強度綜合評價指標體系，開展定量評估，借助空間自相關方法揭示青藏高原旅游活動強度的時空特征，并使用地理探測器進行影響因素解析，以期能為協調青藏高原生態(tài)保護和旅游經濟發(fā)展提供科學參考。

1 研究區(qū)概況

青藏高原位于亞洲中部，是世界上面積最大、海拔最高的高原，約占我國國土面積的四分之一，是我國重要的生態(tài)安全屏障。其周圍環(huán)繞著喜馬拉雅山、帕米爾高原、昆侖山等世界上最高的山脈，是地球上獨一無二的地貌景觀，有“世界屋脊”“地球第三極”“亞洲水塔”之稱。青藏高原行政范圍主要包括青海和西藏的全部，新疆南部、甘肅及云南的西北部和四川省西部（圖1），共計251個縣級行政單位。青藏高原自然和人文景觀豐富多樣，旅游資源本底良好，是旅游者青睞的旅游目的地。隨著交通等基礎設施條件明顯改善，青藏高原旅游業(yè)迅速發(fā)展[27]。快速發(fā)展的旅游業(yè)為青藏高原帶來了可觀的經濟收益，但較為粗放的旅游開發(fā)方式也給當地帶來一定的生態(tài)環(huán)境壓力[28]。

圖1 青藏高原范圍及位置

2 數據來源與研究方法

2.1 數據來源

統計數據來源于2015—2019年青藏高原各省（自治區(qū)）、市（州）、縣（區(qū)）行政單位統計年鑒、《國民經濟和社會發(fā)展統計公報》以及政府工作報告。2019年A級旅游景區(qū)名單來源于各省、自治區(qū)文化和旅游廳。2015—2019年道路數據下載自Open Street Map（https://download.geofabrik.de）。研究地域單元邊界數據下載自“國家青藏高原科學數據中心”（http://data.tpdc.ac.cn），以2019年底的行政區(qū)劃邊界為準。對于部分缺失數據，則使用其上級行政單元數據進行修正補充。

2.2 評價指標體系構建

為綜合評價青藏高原旅游活動強度，本文在參考馬勇等[13]、王振波等[14]、楊璇業(yè)等[29]學者構建的指標體系基礎上，遵循數據獲取的科學性、準確性和可獲取性原則，從旅游經濟強度、旅游人口強度兩個層面共選取7個指標構成青藏高原旅游活動強度綜合評價指標體系（表1）。其中，反映旅游人口強度的指標有3個，反映旅游經濟強度的指標有4個。采用縣級行政單元進行評價，按2019年的行政區(qū)劃，共計251個評價單元。

表1 青藏高原旅游活動強度綜合評價指標體系

旅游人口強度層面包括常住人口密度、旅游人口密度和旅游人口壓力3個指標。旅游者數量增加將直接導致旅游活動強度增強，因此旅游人口強度能很好地衡量各地區(qū)的旅游活動強度。常住人口密度為每平方千米居住的人口數量，以反映該區(qū)縣自身的人口壓力。旅游人口密度為每平方千米接待的旅游者人次，可衡量各區(qū)縣所承載的旅游人口。旅游人口壓力則為旅游人次與當地常住人口之比，反映旅游者給當地帶來的人口壓力。

旅游經濟強度層面包括旅游依存度、旅游貢獻率、旅游經濟密度和第三產業(yè)拉動率4個指標。為使旅游者在旅游地開展旅游活動，促進旅游者消費，帶來更好的旅游體驗和旅游經濟效益，旅游地往往會開發(fā)更多的旅游項目，增加旅游設施，擴寬旅游景區(qū)建設用地等。因此，旅游經濟強度能更全面地反映旅游活動強度，是衡量旅游活動強度的重要層面。旅游依存度為各區(qū)縣旅游收入與地區(qū)生產總值之比，以衡量某區(qū)縣對旅游的依賴程度。旅游貢獻率為當地旅游收入與第三產業(yè)增加值的比值，反映該地旅游業(yè)帶來的收入對第三產業(yè)的貢獻。旅游經濟密度則為每平方千米所產生的旅游經濟收益，以衡量旅游給各區(qū)縣帶來的旅游經濟收益狀況。第三產業(yè)拉動率為第三產業(yè)增加值與地區(qū)生產總值之比，從第三產業(yè)的角度反映旅游業(yè)對地區(qū)經濟的帶動強度。

熵值法能客觀確定各個指標權重，主要用來判斷某個指標的離散程度[30]。因此，本文采用熵值法計算青藏高原旅游活動強度綜合評價指標體系中各指標的權重（表1）。利用線性加權法和熵值法計算得出青藏高原旅游活動強度綜合評價函數f(x)：

式中：wi為指標Xi的權重，i為指標序號。

2.3 空間自相關分析

空間自相關分析結果可用于描述觀測值在整個研究區(qū)域的空間特征。本文采用全局Moran’s I指數檢驗青藏高原旅游活動強度在全局尺度的空間相關程度，其計算公式為[31]：

式中：Ii表示局部Moran’s I指數，在顯著性檢驗下，將顯著性水平達到閾值的空間單元分為高—高、低—低、低—高和高—低四種聚集類型。

2.4 地理探測器

地理探測器是用于探測空間分異性以及揭示其背后驅動力的一組統計學方法[33]。本文主要使用模型中的因子檢測器，用于檢測各因素對青藏高原旅游活動強度的貢獻度，解釋力用q值來度量，模型定義如公式（4）所示：

式中：q為影響因素的探測解釋力；L為因素的分類數；h=1, 2, …,N為研究區(qū)域內的樣本數量；σ2為研究區(qū)域整體樣本的方差值。q取值范圍為[0, 1]，數值越大表明該因素對旅游活動強度差異的解釋能力越強。

3 結果與分析

3.1 青藏高原旅游活動強度時空特征

根據評價指標體系和綜合評價函數f(x)計算得出青藏高原各區(qū)縣旅游活動強度綜合得分，以2015年綜合得分的四分之一分位數、中值、平均值、平均值與標準差之和作為分級標準，將青藏高原2015—2019年旅游活動強度由弱到強分為“弱”“較弱”“中等”“較強”和“強”五個等級，以反映青藏高原旅游活動強度的時間變化和空間格局，得到2015—2019年青藏高原縣域旅游活動強度空間分布（圖2）。

圖2 2015—2019年青藏高原旅游活動強度空間分布

從時間上看，2015—2019年青藏高原旅游活動強度呈增強趨勢。2015—2019年，青藏高原旅游活動強度為“強”等級的區(qū)縣由14個增至20個，面積占比增加；旅游活動強度為“中等”和“較強”等級的區(qū)縣由110個增至150個，面積占比明顯增大；而旅游活動強度為“弱”和“較弱”等級的區(qū)縣數量則從123個減至77個，面積占比逐漸減小。青藏高原旅游活動強度東部和中部變化明顯，大部分區(qū)縣旅游活動強度增強。西部變化微弱，僅個別區(qū)縣旅游活動強度增強。2015年，旅游活動強度為“中等”和“較強”等級的區(qū)縣僅在青藏高原東緣和南部有所分布。而到了2019年，除青藏高原東緣和南部外，中部和西北部的部分區(qū)縣旅游活動強度變?yōu)椤爸械取焙汀拜^強”等級。

從空間上看，青藏高原旅游活動強度呈現出“東南強、中西弱”的空間格局。青藏高原旅游活動強度為“強”等級的區(qū)縣零星分布在東南緣，“中等”和“較強”等級的區(qū)縣集中分布于青藏高原東南部。而旅游活動強度為“弱”和“較弱”等級的區(qū)縣主要分布在中部和西部。旅游活動強度“較強”等級區(qū)縣鄰近“強”等級區(qū)縣分布，旅游活動強度“較弱”等級區(qū)縣鄰近“弱”等級區(qū)縣分布，旅游活動強度“中等”等級區(qū)縣和“較弱”等級區(qū)縣分布最為廣泛。整體來看，青藏高原東南部的區(qū)縣旅游活動強度較強，而中西部的區(qū)縣旅游活動強度較弱。

3.2 青藏高原旅游活動強度空間關聯分析

2015—2019年，青藏高原旅游活動強度全局Moran’s I指數分別為0.58、0.61、0.66、0.62、0.38，數值均大于零且較高，說明青藏高原旅游活動強度在空間上有較強的正相關性，旅游活動強度在空間上呈現出較強的聚集效應，區(qū)縣之間存在相互影響。同時，通過局部空間自相關分析將具有顯著空間效應的區(qū)縣劃分為“高—高”“高—低”“低—高”“低—低”四種聚集類型（圖3）?？梢钥闯觯嗖馗咴糜位顒訌姸鹊木钟蚩臻g聚集特征明顯，總體以高—高聚集類型和低—低聚集類型為主，低—低聚集類型區(qū)縣數量遠多于高—高聚集類型區(qū)縣。其中，高—高聚集類型區(qū)縣集中分布在青藏高原東南緣，數量較為穩(wěn)定，而低—低聚集類型集中分布在青藏高原中部和西部地區(qū)，數量逐年減少。

圖3 2015—2019年青藏高原旅游活動強度局部空間關聯格局

值得注意的是，隨著時間推移，以青藏高原南部的拉薩市城關區(qū)為中心，其外圍低—低聚集類型的區(qū)縣逐漸減少。2015—2017年，拉薩市城關區(qū)周圍鄰近的區(qū)縣不再是低—低聚集類型，2018—2019年，拉薩市城關區(qū)以北、以西和以南的鄰近區(qū)縣也不再是低—低聚集類型，范圍進一步擴大。拉薩市城關區(qū)是青藏高原南部的旅游活動強度高值區(qū)，以其為中心，可以明顯看出旅游活動強度高值區(qū)對周圍區(qū)縣有著較強的輻射帶動作用。

3.3 青藏高原旅游活動強度影響因素解析

通過參考相關研究成果[34-36]，本文從資源本底、交通設施、經濟環(huán)境和服務條件4個因素來探測其對青藏高原旅游活動強度的影響。其中，資源本底因素屬于資源要素，變化較小，由2019年各區(qū)縣景區(qū)得分來反映，景區(qū)得分則由區(qū)縣各A級旅游景區(qū)與其等級之積加總得出（A～AAAAA級旅游景區(qū)賦值為1～5）；交通設施由各區(qū)縣各級道路里程來反映；經濟環(huán)境選取居民儲蓄存款余額表示；服務條件則由第三產業(yè)從業(yè)人員數量反映。以上指標均進行標準化處理后，采用K-means聚類法劃分為1～5級。再利用公式（4）探測2015—2019年各影響因素對青藏高原旅游活動強度的解釋力（表2）。

表2 青藏高原旅游活動強度影響因素地理探測結果

從結果（表2）來看，2015—2019年青藏高原旅游活動強度受到多因素綜合作用。資源本底和交通設施為主導因素，經濟環(huán)境持續(xù)發(fā)揮作用，服務條件的影響愈加明顯。2015—2019年，資源本底和交通設施因素的探測q值均相對較高。2015—2017年，青藏高原旅游活動強度受資源本底影響更大，2018—2019年則是交通設施的影響作用更強。青藏高原旅游資源本底豐富，特點鮮明，為青藏高原旅游業(yè)提供了發(fā)展基礎，是影響青藏高原旅游活動強度的主導因素。交通設施作為產生旅游活動的必要途徑，交通設施狀況和對旅游資源的可達性直接影響著旅游者對旅游目的地的選擇、開展旅游活動的出游范圍和旅游體驗等。穩(wěn)定良好的經濟環(huán)境是公眾產生旅游需求的前提，也是旅游地開發(fā)建設旅游景區(qū)的基礎。旅游業(yè)用較低的生態(tài)環(huán)境成本將旅游資源轉化為經濟效益，讓青藏高原在保護中實現發(fā)展，促進當地經濟良性發(fā)展。服務條件因素的探測q值由2015年的0.078增長至2019年的0.143，可見服務條件的影響因素愈加明顯。更好、更全面的服務條件能夠提升旅游體驗，日漸成為影響旅游活動強度的重要指標。

4 .討論與結論

4.1 討論

青藏高原旅游活動強度呈現“東南強、中西弱”的空間分布受多種因素綜合影響，資源本底和交通設施仍是主導因素。與我國中部和東部地區(qū)相似，旅游活動強度受資源本底、交通設施、經濟環(huán)境和服務條件等因素綜合作用[37-38]。但中東部地區(qū)經濟發(fā)展水平較高，除旅游資源以外，城市發(fā)展、科技水平、網絡媒體和對外開放程度等均會影響地區(qū)旅游活動強度[39-40]。而青藏高原由于自身綜合發(fā)展水平較低，獨特的旅游資源仍是其旅游發(fā)展的核心吸引力。同時，交通設施作為產生旅游活動的必要途徑，對于可進入性相對較難的青藏高原來說也是關鍵的影響因素[41]。

青藏高原旅游活動強度時空特征與旅游業(yè)時空格局一致性較高。旅游業(yè)是產生旅游活動的中介體，為旅游活動提供服務，二者的空間分布可以相互印證。根據評價結果，青藏高原旅游活動強度在空間上有較強的聚集效應，旅游活動強度高等級區(qū)主要分布在青藏高原東南部。與本文旅游活動強度評價結果相一致，青藏高原旅游業(yè)空間分布同樣具有明顯的集聚特征[42]，以西寧和拉薩為核心，對周圍地區(qū)有著明顯的輻射帶動作用，形成西寧旅游帶和拉薩旅游帶[43]。隨著時間推移，高原旅游業(yè)企業(yè)數增長[42]，旅游經濟波動上升[35]，青藏高原旅游活動強度持續(xù)增強。

本文針對青藏高原全域開展定量評估，得到了更高精度的旅游活動強度縣級尺度評價結果，可將平衡生態(tài)保護與旅游經濟發(fā)展的科學參考落實到各個區(qū)縣。相關研究表明，青藏高原旅游生態(tài)環(huán)境壓力等級普遍處于較不安全狀態(tài)，旅游活動對青藏高原現有生態(tài)環(huán)境造成較大壓力[28]。因此，對于旅游活動強度過度的區(qū)縣，旅游開發(fā)時需堅持“生態(tài)優(yōu)先”原則，關注旅游地旅游活動強度，注重旅游地生態(tài)環(huán)境保護，不可超過當地資源環(huán)境承載力，以免造成不可逆的生態(tài)環(huán)境損傷。對于旅游活動強度適宜的區(qū)縣，可通過培育旅游景區(qū)集群，打造區(qū)域旅游景區(qū)增長極，帶動周圍景區(qū)發(fā)展，增加當地旅游經濟收益。而對于旅游活動強度較低的區(qū)縣來說，在穩(wěn)定良好的經濟環(huán)境背景下，通過挖掘特色旅游資源，完善基礎設施，優(yōu)化提升旅游相關服務，將有效促進產業(yè)結構優(yōu)化，增加新的經濟增長點。

4.2 結論

本文通過構建綜合評價指標體系開展了青藏高原旅游活動強度定量評估，在此基礎上借助空間自相關分析方法刻畫了青藏高原旅游活動強度時空特征，并運用地理探測器探討了旅游活動強度的影響因素，得出以下結論：

青藏高原旅游活動強度隨時間持續(xù)增強，在空間上表現出顯著的聚集效應，呈現出“東南強、中西弱”的空間格局，旅游活動強度“強”等級區(qū)對周圍地區(qū)有明顯的輻射帶動作用。青藏高原旅游活動強度受到多種因素綜合作用，資源本底和交通設施為主導因素，經濟環(huán)境因素的影響較為穩(wěn)定且持續(xù)發(fā)揮作用，服務條件因素的影響愈加明顯。