中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)效應及其障礙因素研究

［摘" " 要］在生態(tài)文明建設背景下，科學衡量旅游經(jīng)濟與生態(tài)環(huán)境的協(xié)調(diào)關系，厘清影響其協(xié)調(diào)發(fā)展的障礙因素，具有重要意義。文章首次引入生態(tài)系統(tǒng)服務測算，并在方法層面對現(xiàn)有的靜態(tài)測算模型進行動態(tài)化創(chuàng)新，定量測度2000—2020年中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)度，科學度量其空間演化格局，并識別影響其協(xié)調(diào)發(fā)展的主要障礙因素。結(jié)果顯示：1）中國旅游經(jīng)濟綜合指數(shù)總體呈現(xiàn)波動上升態(tài)勢，表現(xiàn)出東南高、西北低的分布格局；2）中國生態(tài)系統(tǒng)服務有一定的改善，整體呈現(xiàn)增長態(tài)勢，表現(xiàn)出東高西低，由東南向西北遞減的分布格局；3）中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務長期處于高水平磨合階段，耦合協(xié)調(diào)效應穩(wěn)步增強，呈東高西低、南高北低的分布格局，極核特征明顯；4）旅游經(jīng)濟效益、支持服務、旅游經(jīng)濟規(guī)模和調(diào)節(jié)服務是制約旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的主要因素。

［關鍵詞］旅游經(jīng)濟；生態(tài)系統(tǒng)服務；耦合協(xié)調(diào)；障礙度

［中圖分類號］F59

［文獻標識碼］A

［文章編號］1002-5006（2024）04-0093-16

DOI： 10.19765/j.cnki.1002-5006.2024.04.012

0 引言

自然是人類賴以生存發(fā)展的基本條件，人與自然和諧共生是馬克思主義生態(tài)文明觀的核心[1]。旅游業(yè)資源消耗與環(huán)境依賴的雙重屬性決定了旅游經(jīng)濟與生態(tài)環(huán)境間有著相互依存、相互促進的密切關系。一方面，旅游活動對自然生態(tài)系統(tǒng)造成了劇烈干擾，而旅游業(yè)的發(fā)展又促進了旅游地經(jīng)濟的發(fā)展，為生態(tài)保護和環(huán)境恢復提供了資金支持[2-3]。另一方面，友好的環(huán)境不僅能夠提升目的地的旅游吸引力[4]，而且是旅游業(yè)可持續(xù)發(fā)展的重要基礎[5]。黨的二十大報告明確指出“走生態(tài)優(yōu)先、綠色發(fā)展之路”，強調(diào)尊重自然、順應自然和保護自然，必須站在人與自然和諧共生的高度謀劃發(fā)展。旅游業(yè)作為國民經(jīng)濟戰(zhàn)略性支柱產(chǎn)業(yè)，是實現(xiàn)“綠水青山”生態(tài)服務價值的重要渠道[6]。在生態(tài)文明建設背景下，科學衡量旅游經(jīng)濟與生態(tài)系統(tǒng)服務的協(xié)調(diào)關系，厘清影響其協(xié)調(diào)發(fā)展的障礙因素，對于旅游業(yè)創(chuàng)造性地把生態(tài)優(yōu)先和綠色發(fā)展結(jié)合起來，實現(xiàn)人與自然和諧共生具有重要使命。

旅游與生態(tài)環(huán)境的復雜關系一直是旅游地理學、旅游生態(tài)學研究的前沿和熱點問題。自20世紀80年代起，學者們開始關注旅游活動與生態(tài)環(huán)境的辯證關系[7-8]，研究多圍繞旅游環(huán)境容量[9-10]、旅游環(huán)境承載力[11-12]、旅游業(yè)可持續(xù)發(fā)展[13-14]和旅游地生態(tài)安全[15-17]等問題展開。隨著生態(tài)環(huán)境危機日益凸顯，相關學者對粗放式的旅游經(jīng)濟增長模式進行理性反思，旅游經(jīng)濟與生態(tài)環(huán)境的協(xié)調(diào)關系逐漸成為研究的焦點。當前，相關研究的特點表現(xiàn)為以下幾個方面：一是研究方法由定性描述向定量研究轉(zhuǎn)變，耦合協(xié)調(diào)度模型應用最為廣泛；二是研究尺度以微觀和中觀尺度為主，如對曲靖市[18]、神農(nóng)架林區(qū)[19]、浙江省[20]、長江與珠江流域省市[21-22]展開研究；三是研究內(nèi)容從雙系統(tǒng)協(xié)調(diào)關系向多系統(tǒng)耦合協(xié)調(diào)拓展。一是相關研究均使用綜合指數(shù)法對旅游經(jīng)濟和生態(tài)環(huán)境進行間接評價[18-23]，缺乏對更科學、精準的直接評估方法的探索；二是多數(shù)研究關注旅游經(jīng)濟與生態(tài)環(huán)境的協(xié)調(diào)發(fā)展情況，忽視了對協(xié)調(diào)發(fā)展影響因素的進一步探究；三是具有更高戰(zhàn)略價值的全國尺度的研究也較為不足。

生態(tài)系統(tǒng)服務是指人類從自然中直接或者間接獲得的所有惠益[24]，是理解自然與社會之間關系的橋梁，而旅游業(yè)正是這種關系的重要表現(xiàn)形式[25]。生態(tài)系統(tǒng)服務已經(jīng)形成了一個相對成熟的研究框架，提供了一種具體、準確地量化生態(tài)環(huán)境的直接評估方法[26]，運用生態(tài)系統(tǒng)服務測量生態(tài)環(huán)境質(zhì)量，能夠解決現(xiàn)有研究使用綜合指數(shù)法帶來的問題。此外，使用生態(tài)系統(tǒng)服務測算還具有以下優(yōu)勢：首先，現(xiàn)有研究使用的生態(tài)環(huán)境測量指標建立在研究者的主觀判斷選取之上[27]，既存在未能全面測算生態(tài)環(huán)境的片面性問題，已有不同研究的討論效果也不夠精準，而生態(tài)系統(tǒng)服務是當前關于生態(tài)環(huán)境測算較為全面的成熟框架，能夠解決上述問題；其次，生態(tài)環(huán)境為人類社會提供的服務屬于公共物品，具有外部性，即非排他性和非競爭性[28-30]，已有研究傾向于評估生態(tài)環(huán)境的直接影響，忽視了如環(huán)境惡化導致的負外部性和空氣凈化帶來的正外部性等間接影響[31-32]，生態(tài)系統(tǒng)服務不僅可以識別生態(tài)環(huán)境的直接影響，還能量化生態(tài)環(huán)境被低估的外部性這一間接影響；最后，已有研究選取的生態(tài)環(huán)境評估指標具有不同的量綱[13-21]，不能實現(xiàn)同一計量單位下環(huán)境成本和財務效益的比較[33]，而生態(tài)系統(tǒng)服務測算可以解決這一問題。具體而言，生態(tài)系統(tǒng)服務價值是生態(tài)系統(tǒng)服務的幣值表現(xiàn)形式，通過經(jīng)濟學市場價格的方式對生態(tài)環(huán)境價值、質(zhì)量及其變化進行量化，將生態(tài)環(huán)境成本和效益貨幣化，從而實現(xiàn)不同指標間的比較以及不同區(qū)域評估結(jié)果的對接[34-35]。

生態(tài)系統(tǒng)服務價值測算可分為價值量和物質(zhì)量兩種方式[36-37]。價值量評估方式中的當量因子法認為，同種類型土地具有相同的價值當量[27]，其評估結(jié)果不僅易于被決策者理解，還能增進公眾對自然資本的認可[34-35]。Costanza等最早運用當量因子法測算了全球生態(tài)系統(tǒng)服務價值[38]。在此基礎上，國內(nèi)學者編制的中國生態(tài)系統(tǒng)服務價值當量因子表被廣泛應用[39-40]。然而生態(tài)系統(tǒng)服務受生態(tài)機制的調(diào)控，涉及生物體及其物理、化學環(huán)境之間復雜的相互作用[41]，可能受到氣候、土地利用和其他人類活動變化的影響[42]。因此，這一靜態(tài)系數(shù)表方法忽視了生態(tài)系統(tǒng)服務的動態(tài)性和多樣性[43-44]，無法顯示區(qū)域自然資源和社會經(jīng)濟特征的差異[45-46]，不能適用于所有地區(qū)。

有鑒于此，本文將生態(tài)系統(tǒng)服務引入旅游經(jīng)濟與生態(tài)環(huán)境協(xié)調(diào)發(fā)展研究中，更加科學、客觀地度量生態(tài)環(huán)境為人類社會提供的功能服務，推進二者協(xié)調(diào)發(fā)展的實證研究。進一步地，鑒于社會經(jīng)濟和自然資源的區(qū)域差異性，本文將社會經(jīng)濟水平和資源稀缺性引入生態(tài)系統(tǒng)服務價值測算模型，對傳統(tǒng)靜態(tài)測算模型進行了延伸，構建了動態(tài)的生態(tài)系統(tǒng)服務價值測算框架，嘗試在戰(zhàn)略價值更高的全國尺度對2000—2020年旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)關系及其時空格局展開研究。同時，利用障礙度模型識別影響旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的主要障礙因素，以期為科學、合理地制定旅游業(yè)綠色可持續(xù)發(fā)展路徑提供借鑒與參考。

1 旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)機理分析

生態(tài)系統(tǒng)服務是連接社會經(jīng)濟系統(tǒng)和自然生態(tài)系統(tǒng)的重要橋梁，在判斷旅游活動與自然環(huán)境是否和諧共生中起到關鍵作用[38]。首先，人與自然互動過程中，人類從自然中獲得的惠益包含旅游體驗[26]，就這個意義而言，旅游業(yè)是生態(tài)系統(tǒng)服務對人類福祉做出貢獻的經(jīng)濟受益者。其次，旅游業(yè)作為幸福產(chǎn)業(yè)不僅僅受生態(tài)系統(tǒng)服務的影響，還創(chuàng)造性地重構了生態(tài)系統(tǒng)服務的產(chǎn)生過程，是生態(tài)系統(tǒng)服務實現(xiàn)人類福祉的重要渠道[47]。最后，旅游業(yè)發(fā)展依賴于生態(tài)系統(tǒng)服務，但旅游經(jīng)濟的快速發(fā)展不可避免地驅(qū)動生態(tài)系統(tǒng)服務產(chǎn)生積極或消極的變化[48-49]。因此，有必要正確認識旅游經(jīng)濟與生態(tài)系統(tǒng)服務之間的復雜關系。基于一般系統(tǒng)理論，耦合是指兩個或多個子系統(tǒng)相互作用與相互影響的現(xiàn)象[50]，同時，子系統(tǒng)間錯綜復雜的交互作用隨時間推移而動態(tài)變化[51]。耦合協(xié)調(diào)是指各子系統(tǒng)或要素之間相互影響的強弱程度，以及反映子系統(tǒng)間如何從無序演化為有序發(fā)展的過程[52]。旅游經(jīng)濟與生態(tài)系統(tǒng)服務子系統(tǒng)通過各要素耦合互饋，削弱系統(tǒng)間的權衡效應，加強協(xié)同效應，最終實現(xiàn)系統(tǒng)間由無序向有序發(fā)展的演變[53]。因此，本文揭示旅游經(jīng)濟與生態(tài)系統(tǒng)服務的互動機制具有重要意義。

旅游經(jīng)濟系統(tǒng)與生態(tài)系統(tǒng)服務系統(tǒng)協(xié)調(diào)發(fā)展的本質(zhì)是旅游經(jīng)濟發(fā)展合理配置和利用生態(tài)系統(tǒng)服務功能，釋放生態(tài)系統(tǒng)服務改善人類福祉的效用，實現(xiàn)旅游經(jīng)濟發(fā)展與生態(tài)系統(tǒng)服務效益的雙贏。從旅游經(jīng)濟對生態(tài)系統(tǒng)服務的作用來看，一方面，旅游經(jīng)濟發(fā)展能促進旅游目的地的經(jīng)濟發(fā)展，不僅有助于基礎設施升級以及產(chǎn)業(yè)結(jié)構優(yōu)化，還能間接改善生態(tài)系統(tǒng)服務水平，提升生態(tài)系統(tǒng)服務的整體能力[49， 54]。另一方面，旅游經(jīng)濟的快速發(fā)展不可避免地造成生態(tài)系統(tǒng)服務能力的衰退，如旅游基礎設施的建設改變了旅游目的地土地覆蓋和土地利用類型，導致水體污染、大氣污染、生物多樣性減少和景觀破碎化等問題[33， 55-56]。此外，旅游經(jīng)濟效益驅(qū)動旅游企業(yè)盲目發(fā)展，忽視旅游開發(fā)推動土地利用變化帶來的消極影響，對生態(tài)系統(tǒng)服務產(chǎn)生更嚴重的脅迫效應[57]。從生態(tài)系統(tǒng)服務對旅游經(jīng)濟的作用來看，一方面，旅游經(jīng)濟的發(fā)展以生態(tài)系統(tǒng)服務為基礎，如供給服務為游客提供食物、淡水和能源等[34]，水域、林地和草地等為旅游經(jīng)濟發(fā)展提供生態(tài)空間[58]，良好的大氣環(huán)境、生物多樣性等調(diào)節(jié)服務提高旅游目的地形象[59]，從而吸引更多的游客開展旅游活動，促進旅游經(jīng)濟增長。另一方面，良好的生態(tài)系統(tǒng)服務是旅游經(jīng)濟可持續(xù)發(fā)展的先決條件，而生態(tài)系統(tǒng)服務能力的喪失和退化又會對旅游經(jīng)濟發(fā)展產(chǎn)生制約效應，如大氣環(huán)境、水資源供給和美學價值等調(diào)節(jié)服務水平、供給服務水平和文化服務水平的下降導致游客滿意度的降低[60]，從而制約旅游經(jīng)濟發(fā)展。因此，旅游經(jīng)濟與生態(tài)系統(tǒng)服務各系統(tǒng)要素之間相互協(xié)調(diào)、共生和良性耦合，才能實現(xiàn)旅游經(jīng)濟綠色可持續(xù)發(fā)展。

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

2.1 旅游經(jīng)濟評估方法

2.1.1" " 指標體系構建

本文遵循系統(tǒng)性、科學性、代表性和數(shù)據(jù)可獲得性的原則，以動態(tài)反映旅游經(jīng)濟系統(tǒng)要素和結(jié)構為目標，構建旅游經(jīng)濟指標體系。首先，對相關文獻進行比較和借鑒[18-21]，選擇已有研究中高頻使用的相關指標。其次，在旅游經(jīng)濟研究課題組展開討論，構建旅游經(jīng)濟的初始指標庫。最后，運用Stata軟件剔除高度相關的指標，但保留對指標體系影響較大的指標，構建出由旅游經(jīng)濟效益、旅游經(jīng)濟規(guī)模和旅游經(jīng)濟強度3個準則層一級指標和18個指標層二級指標組成的旅游經(jīng)濟評價指標體系（表1）。

2.1.2" " 測算方法

采用綜合評價模型對旅游經(jīng)濟進行測算。主要步驟：1）運用極值法對指標數(shù)據(jù)進行標準化處理；2）通過面板熵值法計算旅游經(jīng)濟各指標權重；3）運用多重線性加權法測算我國各省旅游經(jīng)濟的綜合水平。具體公式參見相關文獻[18-19]。

2.2 生態(tài)系統(tǒng)服務評估方法

2.2.1" " 標準單位當量因子價值的確定

標準單位生態(tài)系統(tǒng)服務價值當量因子，是指中國農(nóng)田年均單位面積自然糧食產(chǎn)量的經(jīng)濟價值[39-40]。自然糧食產(chǎn)量一般被確定為每公頃平均實際糧食產(chǎn)量的1/7，公式如下：

[E=17×121s=20002020rsk×k=1nPsk×Qsk×Fsk] （1）

式（1）中，[E]為標準單位當量因子價值（元/hm2）；[rsk]指農(nóng)產(chǎn)品價格指數(shù)，將每年的價格轉(zhuǎn)換為2020年的固定價格；[Psk]為第[k]種作物在第[s]年種植的面積比重（%）；[n]為研究區(qū)主要作物數(shù)量；[Qsk]是第[k]種作物在第[s]年的每公頃產(chǎn)量（kg/hm2）；[Fsk]為第[k]種作物在第[s]年的平均市場價格（元/kg）。由于種植面積和農(nóng)產(chǎn)品價格易受市場因素和自然因素的影響，本文選取2000—2020年的平均值作為標準單位當量因子值。

2.2.2" " 生態(tài)系統(tǒng)服務價值系數(shù)的調(diào)整

1）社會經(jīng)濟調(diào)整系數(shù)

本文將支付意愿納入生態(tài)系統(tǒng)服務價值的評估體系，以更好地將生態(tài)系統(tǒng)服務價值的實際變化與社會經(jīng)濟發(fā)展水平的變化相適應。生長曲線模型通常用于表達人們的支付意愿與社會經(jīng)濟發(fā)展階段之間的關系[43]，公式如下：

[l=L1+e-1En-3] （2）

式（2）中，[l]為社會經(jīng)濟發(fā)展水平；[L]指社會發(fā)展水平較高時人們的支付意愿，設置為1；[En]為恩格爾系數(shù)；[e]為自然對數(shù)。

然而，恩格爾系數(shù)并不能很好地反映不同區(qū)域發(fā)展的差異[44]。因此，本文在已有研究的基礎上，采用區(qū)域收入水平和城市化水平對調(diào)整系數(shù)進行進一步修正。調(diào)整后的社會經(jīng)濟發(fā)展水平計算公式如下：

[l′=l×uU×diDI] （3）

式（3）中，[l′]為調(diào)整后的社會經(jīng)濟發(fā)展水平；[u]和[U]分別指研究區(qū)域城市化水平和全國城市化水平；[di]為研究地區(qū)的人均可支配收入（元/人）；[DI]是全國人均可支配收入（元/人）。

最后，社會經(jīng)濟調(diào)整系數(shù)計算公式如下：

[D=l′slC] （4）

式（4）中，D為社會經(jīng)濟調(diào)整系數(shù)；[l′s]是研究區(qū)域的社會經(jīng)濟發(fā)展水平；[lC]為中國社會經(jīng)濟發(fā)展水平。

2）資源稀缺性調(diào)整系數(shù)

資源稀缺程度，是指區(qū)域社會經(jīng)濟發(fā)展過程中生態(tài)資源的供需關系[61]。當生態(tài)資源的供給小于需求時，資源稀缺出現(xiàn)。生態(tài)資源稀缺性越高，人們?yōu)樯鷳B(tài)資源付費的意愿就越高。由于人是生態(tài)系統(tǒng)的服務對象，因此選擇人口密度來建立資源稀缺系數(shù)[62]。此外，為便于計算，資源稀缺系數(shù)由人口密度的對數(shù)構建，公式如下：

[S=logplogP] （5）

式（5）中，[S]是資源稀缺系數(shù)；[p]為研究區(qū)人口密度；[P]是全國平均人口密度。

2.2.3" "生態(tài)系統(tǒng)服務價值測算

1）非建設用地生態(tài)系統(tǒng)服務價值估算

基于當量因子法和上述調(diào)整系數(shù)[38， 63]，構建非建設用地生態(tài)系統(tǒng)服務價值動態(tài)評估模型：

[Vkf=E×D×S×Ekf]" （6）

[ESVk=f=1mAk×Vkf] （7）

[ESVf=k=1nAk×Vkf] （8）

[ESV=k=1mf=1nAk×Vkf] （9）

式（6）中，[Vkf]是第[k]種土地利用類型的第[f]項生態(tài)系統(tǒng)服務功能價值系數(shù)；[E]為標準單位價值系數(shù)（元/hm2）；[D]為研究區(qū)域社會經(jīng)濟調(diào)整系數(shù)；[S]為資源稀缺調(diào)整系數(shù)；[Ekf]為第[k]種土地利用類型的第[f]項服務功能的當量因子。式（7）中，[ESVk]為第[k]種土地利用類型的生態(tài)系統(tǒng)服務價值；[Ak]是第[k]種土地利用類型的面積（hm2）。式（8）中，[ESVf]為第[f]項生態(tài)系統(tǒng)服務功能的生態(tài)系統(tǒng)服務價值。式（9）中，[ESV]為生態(tài)系統(tǒng)服務總價值。

2）建設用地生態(tài)系統(tǒng)服務價值估算

建設用地提供的生態(tài)系統(tǒng)服務在許多研究中被忽視，而建設用地提供的服務價值既有積極的一面，也有消極的一面。本文將建設用地生態(tài)系統(tǒng)服務價值分為正向價值和負向價值進行估算。旅游景區(qū)、公園、住宅綠地、綠化帶廊道等城市綠地多種植草木，與草地、林地類似[64]。因此，建設用地正向價值系數(shù)由這兩個土地類型價值系數(shù)的平均值決定。公式如下：

[Ecj=GC×12（Ewj+Egj）] （10）

[Vcj=Ecj×E] （11）

式（10）中，[Ecj]、[Ewj]和[Egj]分別為建設用地、林地、草地的第[j]項生態(tài)系統(tǒng)服務功能的當量因子；[GC]為城市綠地覆蓋率（%）。式（11）中，[Vcj]是單位建設用地正向生態(tài)系統(tǒng)服務價值。

建設用地對生態(tài)環(huán)境造成諸多負面影響，主要表現(xiàn)為：廢氣、廢水和固體廢物的排放會破壞供水、廢物處理、水文調(diào)節(jié)和氣體調(diào)節(jié)等生態(tài)系統(tǒng)服務。本文采用預防成本法和替代成本法，對建設用地負向生態(tài)系統(tǒng)服務進行評價[65-66]。其中，采用預防成本法對水文調(diào)節(jié)、廢物處理和氣體調(diào)節(jié)的生態(tài)系統(tǒng)服務價值進行估算，替代成本法對供水的生態(tài)系統(tǒng)服務價值進行估算。

單位建設用地供水負向生態(tài)系統(tǒng)服務價值計算公式如下：

[Vws=-W×PwA] （12）

式（12）中，[Vws]為單位建設用地供水的生態(tài)系統(tǒng)服務價值（元/hm2/年）；[W]為用水總量（m3/年）；[Pw]是耗水單價（元/m3）；[A]為建設用地總面積。

單位建設用地氣體調(diào)節(jié)負向生態(tài)系統(tǒng)服務價值計算公式如下：

[Vga=-g=1xQgCgA] （13）

式（13）中，[Vga]為單位建設用地氣體調(diào)節(jié)的生態(tài)系統(tǒng)服務價值（元/ hm2/年）；[Qg]是第[g]型廢氣量（m3/年）；[Cg]為第[g]型廢氣處理成本（元/m3）。

單位建設用地水文調(diào)節(jié)負向生態(tài)系統(tǒng)服務價值計算公式如下：

[Vha=-w=1tQwCwA] （14）

式（14）中，[Vha]為單位建設用地水文調(diào)整的生態(tài)系統(tǒng)服務價值（元/ hm2/年）；[Qw]為第[w]型廢水量（t/年）；[Cw]為第[w]型廢水的處理費用（元/t）。

單位建設用地固體廢物處理負向生態(tài)系統(tǒng)服務價值計算公式如下：

[Vpe=-r=1mQrCrA] （15）

式（15）中，[Vpe]為單位建設用地廢物處理的生態(tài)系統(tǒng)服務價值（元/hm2/年）；[Qr]為第[r]類固體廢物數(shù)量（噸/年）；[Cr]是第[r]類固體廢物處理費用（元/噸）。

2.3 耦合協(xié)調(diào)度模型和障礙度模型

本文運用耦合協(xié)調(diào)度模型計算旅游經(jīng)濟與生態(tài)系統(tǒng)服務的協(xié)調(diào)發(fā)展水平，具體公式參見相關文獻[18-20]，耦合協(xié)調(diào)度劃分標準參照廖重斌的研究[22]。本文采用障礙度模型，對旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的主導障礙因素進行測算與排序，具體公式參見相關文獻[67]。

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

本文數(shù)據(jù)主要包括土地利用數(shù)據(jù)、環(huán)境污染數(shù)據(jù)、社會經(jīng)濟數(shù)據(jù)和旅游經(jīng)濟數(shù)據(jù)（不含港澳臺）。其中，2000年、2005年、2010年、2015年和2020年5期土地利用數(shù)據(jù)來源于中國土地利用監(jiān)測數(shù)據(jù)庫，分辨率為[100m×100m]；環(huán)境污染數(shù)據(jù)來源于《中國生態(tài)環(huán)境統(tǒng)計年報（2000—2020）》和《中國城市統(tǒng)計年鑒（2000—2020）》；社會經(jīng)濟數(shù)據(jù)如國內(nèi)生產(chǎn)總值、城市化率、糧食產(chǎn)量和農(nóng)產(chǎn)品銷售價格等來源于中國經(jīng)濟社會大數(shù)據(jù)統(tǒng)計平臺；旅游經(jīng)濟數(shù)據(jù)主要來源于《中國旅游統(tǒng)計年鑒（2000—2020）》《中國區(qū)域經(jīng)濟統(tǒng)計年鑒（2000—2020）》《國民經(jīng)濟和社會發(fā)展統(tǒng)計公報（2000—2020）》。采用線性插值法補全缺失值。根據(jù)中國科學院土地利用覆蓋分類體系，建立林地、草地、耕地、濕地、水域、建設用地（交通用地、居民區(qū)等）和未利用土地（裸地和鹽堿地）7個土地類型。研究區(qū)域行政邊界數(shù)據(jù)來源于國家測繪地理信息局標準地圖服務網(wǎng)站（http：//bzdt.ch.mnr.gov.cn/）。

3 結(jié)果分析

3.1 旅游經(jīng)濟測度結(jié)果

2000—2020年中國省際旅游經(jīng)濟總體呈波動上升態(tài)勢，旅游經(jīng)濟指數(shù)由0.043上升至0.162，累計增長12%，表明中國旅游經(jīng)濟在過去20年穩(wěn)步提升（圖1），與王兆峰和陳青青[68]、王勝鵬等[69]研究結(jié)論相似。同時，中國省際旅游經(jīng)濟演化過程具有明顯的階段性特征，大致可分為5個階段。具體來看，2000—2002為平穩(wěn)上升期，年均增幅為7.18%；2002—2003年由于“非典”的暴發(fā)呈下降趨勢，特別是廣東、北京出現(xiàn)明顯下跌；2003—2012年為快速上升期，年均增幅為18.50%；2012—2019年呈先下降而后持續(xù)增長趨勢；2019—2020年受新冠疫情的沖擊，旅游經(jīng)濟指數(shù)呈現(xiàn)斷崖式下跌。借鑒已有研究的劃分標準，進一步將樣本省份劃分為東部、中部和西部3個區(qū)域子樣本[70]。東部、中部和西部3個區(qū)域的旅游經(jīng)濟指數(shù)演化軌跡與整體趨勢大體保持一致，均呈總體上升態(tài)勢，東部地區(qū)旅游經(jīng)濟指數(shù)高于全國平均水平。分區(qū)域看，東部地區(qū)旅游經(jīng)濟漲幅最大（累計增長13.34%），其次為中部地區(qū)（累計增長12.54%），而西部地區(qū)的漲幅最低（累計增長9.55%），與李維維等[71]、王兆峰和李靜怡[72]的研究結(jié)論相似。上述結(jié)果表明我國旅游經(jīng)濟發(fā)展面臨著突出的不均衡現(xiàn)象，但中西部地區(qū)后發(fā)增長優(yōu)勢明顯。

采用自然斷點法將研究時序內(nèi)旅游經(jīng)濟發(fā)展水平分為高水平區(qū)、較高水平區(qū)、中等水平區(qū)、較低水平區(qū)和低水平區(qū)5個等級，選擇2000年、2005年、2010年、2015年和2020年5個截面進行可視化處理，展現(xiàn)中國旅游經(jīng)濟演化特征（圖2）?？梢?，2000—2020年中國旅游經(jīng)濟變化比較明顯，總體呈東南高、西北低的分布格局?？赡軞w因于：其一，整體而言東南地區(qū)社會經(jīng)濟發(fā)達，旅游業(yè)起步較早，旅游經(jīng)濟發(fā)展的空間溢出效應使東南地區(qū)旅游經(jīng)濟空間集聚性增強，旅游經(jīng)濟綜合水平整體較高；其二，旅游經(jīng)濟快速發(fā)展導致人地矛盾凸顯，倒逼旅游產(chǎn)業(yè)轉(zhuǎn)型升級、綠色發(fā)展，進而促進東南地區(qū)旅游經(jīng)濟的空間集聚效應進一步增強。具體來看，2000—2010年低水平區(qū)和中等水平區(qū)趨于收縮，內(nèi)蒙古、山西、重慶、貴州和江西由低水平區(qū)轉(zhuǎn)移為較低水平區(qū)，黑龍江、河北和云南由中等水平區(qū)轉(zhuǎn)移為較低水平區(qū)。2010—2020年高水平區(qū)和較高水平區(qū)明顯擴張并且趨于集中，四川、山東、江蘇、上海、浙江和湖南由中等和較高水平區(qū)轉(zhuǎn)移為高水平區(qū)，較低和中等水平區(qū)略有收縮。

3.2 生態(tài)系統(tǒng)服務測度結(jié)果

2000—2020年中國生態(tài)系統(tǒng)服務價值呈波動上升態(tài)勢（表2），由2000年的290 669.3億元上升到2020年的303 179.1億元，共計增長12509.8億元，年均增長625.5億元，增長率為4.3%。具體來看，2000—2020年生態(tài)系統(tǒng)服務總值歷經(jīng)了前緩增后驟減再增加的“N”形變化過程。由2000年的290 669.3億元增加到2010年的291 422.3億元，共計增加752.97億元，增加率為0.26%，而后階段性驟減到2015年的290 681.6億元，減少率為0.25%，幾乎回落到2000年的水平，張中浩等[73]、劉志濤等[74]也發(fā)現(xiàn)相似特征。究其原因，中國伴隨著改革開放環(huán)境質(zhì)量經(jīng)歷了從良好、惡化到總體好轉(zhuǎn)的演化進程。2010年環(huán)境惡化的癥狀開始急劇顯現(xiàn)，2010—2015年間，生態(tài)系統(tǒng)服務價值下降740.7億元。中國政府痛定思痛之后，走上了艱難的轉(zhuǎn)型之路。2013年，我國發(fā)布《大氣污染防治行動計劃》，要求重點區(qū)域的省市削減煤炭消費量，在2017年底實現(xiàn)全國空氣質(zhì)量明顯改善。2014年，中國公開做出關于碳排放碳達峰的國際承諾，有效遏制了生態(tài)環(huán)境惡化。2018年3月，“生態(tài)文明”理念被寫入憲法，有助于進一步塑造中國綠色發(fā)展進程。因此，2015—2020年間，中國生態(tài)系統(tǒng)服務價值扭虧為盈，累計增長12 497.5億元，占總增長量的99.9%，表明我國的生態(tài)保護措施卓有成效，人地關系正逐步改善。

從生態(tài)系統(tǒng)服務類型來看，近20年來4類生態(tài)系統(tǒng)服務相對比例沒有明顯變化，生態(tài)系統(tǒng)服務價值結(jié)構較為穩(wěn)定，支持服務價值占比最高，達到54%左右，調(diào)節(jié)服務價值次之，達到40%左右。以上結(jié)果表明，我國自然生態(tài)系統(tǒng)主要發(fā)揮維持人類賴以生存的自然環(huán)境與效用、消解人類污染與保障生態(tài)系統(tǒng)穩(wěn)定的作用。具體來看，2000—2020年我國4類生態(tài)系統(tǒng)服務價值均有不同程度的增加，其中休閑服務增加幅度最大，共計增加6.22%。階段性變化中，支持服務與休閑服務價值歷經(jīng)了先減少后略增，再驟降后提升的“W”形波動變化；供給服務和調(diào)節(jié)服務價值則均呈現(xiàn)前增后減、再增加的“N”形變動。

將研究時序內(nèi)各省份生態(tài)系統(tǒng)服務價值分為5個等級，繪制生態(tài)系統(tǒng)服務價值空間分布情況（圖3）。從空間分布來看，總體呈東高西低、由東南向西北遞減的分布格局，高值區(qū)域主要分布于我國的東北和西南地區(qū)，不同時期生態(tài)系統(tǒng)服務價值存在明顯的空間分異（圖3）。原因可能源于：其一，林地生態(tài)價值當量高于耕地、草地和水體等[75]，我國東北地區(qū)和西南地區(qū)多為山地丘陵地形，森林覆蓋率較高，如大小興安嶺是我國重要的生態(tài)屏障，因此生態(tài)系統(tǒng)服務價值較高；其二，隨著人口遷移、土地利用變化、城市擴展和經(jīng)濟發(fā)展等人類活動增加，導致生態(tài)系統(tǒng)服務發(fā)生變化且空間差異顯著?？傮w來看，2000—2020年高值區(qū)和低值區(qū)均趨于收縮，空間分布格局整體變化較小。其中，2000—2010年高值區(qū)主要分布于黑龍江、廣西和云南，次高值區(qū)分布于長江中下游省份，低值區(qū)主要分布于黃河沿線省份；2010—2020年高值區(qū)和低值區(qū)處于收縮狀態(tài)，廣西由高值區(qū)轉(zhuǎn)為為次高值區(qū)，河北和海南省分別由低值區(qū)轉(zhuǎn)為為次低值區(qū)和次高值區(qū)域。

3.3 旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)關系

基于耦合協(xié)調(diào)度模型，測算出中國省際2000—2020年旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)度，并分為東部、中部和西部3個地區(qū)子樣本進一步比較（表3）。

首先，從耦合度指數(shù)（C）來看，中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合度在研究期內(nèi)整體呈緩慢下降趨勢，耦合度指數(shù)在0.85至0.88之間，長期處于高水平磨合階段，表明中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務具有較強的關聯(lián)性，相互作用得到有效激發(fā)，耦合發(fā)展相對穩(wěn)定。分區(qū)域來看，東部、中部和西部3個區(qū)域的耦合度在研究時序內(nèi)均呈現(xiàn)下降態(tài)勢，中部地區(qū)耦合度水平最高，其次是西部地區(qū)，而東部地區(qū)相對較低。究其原因，隨著外部經(jīng)濟環(huán)境和自然資源狀況的不斷變化，我國旅游經(jīng)濟發(fā)展與生態(tài)系統(tǒng)服務之間的矛盾也不斷增加，旅游經(jīng)濟發(fā)展在一定程度上改變了土地利用和土地覆蓋類型，導致生態(tài)系統(tǒng)服務水平的下降。同時，生態(tài)系統(tǒng)服務水平衰退抑制了旅游經(jīng)濟的發(fā)展。中部和西部地區(qū)擁有良好的自然生態(tài)系統(tǒng)本底，生態(tài)系統(tǒng)服務水平整體處于較高水平，旅游經(jīng)濟后發(fā)優(yōu)勢明顯，從而旅游經(jīng)濟與生態(tài)系統(tǒng)服務二者的互動效應較高。東部地區(qū)旅游經(jīng)濟發(fā)展水平較高，旅游經(jīng)濟的快速發(fā)展導致生態(tài)系統(tǒng)服務與旅游經(jīng)濟二者耦合度不容樂觀，交互效應和作用力低于整體水平且有所下滑。因而，旅游經(jīng)濟在發(fā)展過程中需要與時俱進地調(diào)整發(fā)展模式和系統(tǒng)結(jié)構。

其次，從耦合協(xié)調(diào)度指數(shù)（D）來看，2000—2020年旅游經(jīng)濟與生態(tài)系統(tǒng)服務的耦合協(xié)調(diào)度呈穩(wěn)步上升態(tài)勢，從0.54上升至0.60，耦合協(xié)調(diào)等級由瀕臨失調(diào)躍升為勉強協(xié)調(diào)階段，但研究時序內(nèi)二者整體的耦合協(xié)調(diào)度水平不高，系統(tǒng)間協(xié)調(diào)發(fā)展仍待進一步增強?？赡艿脑蚴牵何覈糜谓?jīng)濟發(fā)展水平還有較大的提升空間，旅游經(jīng)濟的快速發(fā)展在一定程度上以犧牲目的地生態(tài)系統(tǒng)服務效益為代價，在發(fā)展過程中缺乏科學合理規(guī)劃，導致系統(tǒng)間協(xié)調(diào)發(fā)展水平不高。分區(qū)域來看，中部地區(qū)耦合協(xié)調(diào)度高于全國平均水平，實現(xiàn)了初級協(xié)調(diào)向中級協(xié)調(diào)階段的躍遷，東部地區(qū)趨近于全國平均水平，而西部地區(qū)距離全國平均水平仍有較大差距。由此可見，隨著旅游經(jīng)濟和生態(tài)系統(tǒng)服務綜合水平的不斷提高，二者實現(xiàn)由無序發(fā)展向接近有序發(fā)展的轉(zhuǎn)變，這與“綠色發(fā)展”“美麗中國”“生態(tài)文明建設”等政策指導息息相關，特別是綠色可持續(xù)發(fā)展理念推動旅游業(yè)的結(jié)構性改革，旅游經(jīng)濟高質(zhì)量發(fā)展初有成效，旅游經(jīng)濟對生態(tài)系統(tǒng)服務的脅迫效應有所下降，旅游經(jīng)濟與生態(tài)系統(tǒng)服務的協(xié)同效應趨于良好，系統(tǒng)間耦合協(xié)調(diào)發(fā)展有所改善。中部地區(qū)雖然旅游業(yè)起步較晚，旅游經(jīng)濟綜合水平相對較低，但基于良好自然生態(tài)狀況的生態(tài)旅游、康養(yǎng)旅游等新型旅游形式的發(fā)展，以綠色發(fā)展為契機促進了生態(tài)系統(tǒng)服務水平的不斷增強，加快了系統(tǒng)間協(xié)調(diào)發(fā)展，實現(xiàn)中部地區(qū)由初級協(xié)調(diào)階段向中級協(xié)調(diào)發(fā)展階段的轉(zhuǎn)變，二者耦合協(xié)調(diào)發(fā)展成效顯著。

本文選擇2000年、2005年、2010年、2015年和2020年作為橫斷面，將耦合協(xié)調(diào)度劃分為10個等級，從縱向維度比較分析旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)關系的空間分異和演化特征（圖4）。結(jié)果顯示，中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)度呈現(xiàn)“東高西低、南高北低”的分布特征，極核特征明顯?？傮w來看，2000—2020年高值區(qū)主要集中于珠江流域和長江流域省份，并向區(qū)域外緣延伸；低值區(qū)則多分布于西北地區(qū)；高值區(qū)趨于擴張，低值區(qū)趨于收縮，空間分布格局變化較大。

具體來看，2000年耦合協(xié)調(diào)度主要表現(xiàn)為5種類型：25.81%（8個）的省份發(fā)展初級失調(diào)，22.58%（7個）的省份處于勉強協(xié)調(diào)階段，16.13%（5個）的省份發(fā)展瀕臨失調(diào)，12.90%（4個）的省份已輕度失調(diào)，而良好協(xié)調(diào)的省份僅占9.68%（3個），耦合協(xié)調(diào)階段整體呈現(xiàn)出由東向西逐漸遞減的規(guī)律，與這一時期我國的旅游經(jīng)濟發(fā)展水平現(xiàn)狀基本相符。這主要由該時期我國旅游經(jīng)濟和生態(tài)系統(tǒng)服務均處于較低水平所致，尤其是寧夏、青海、西藏和甘肅等省份處于西北經(jīng)濟塌陷區(qū)，雖然旅游經(jīng)濟發(fā)展一定程度上改善了區(qū)域經(jīng)濟狀態(tài)，但粗放型發(fā)展模式使得生態(tài)系統(tǒng)服務嚴重衰退，系統(tǒng)間耦合協(xié)調(diào)效應較低。2000—2010年良好協(xié)調(diào)區(qū)域和中度失調(diào)區(qū)域均趨于收縮，而中級協(xié)調(diào)區(qū)域趨于擴張，總體空間演變較?。辉颇虾秃邶埥闪己脜f(xié)調(diào)階段轉(zhuǎn)為為中級協(xié)調(diào)階段，青海由中度失調(diào)階段躍升為輕度失調(diào)階段，湖南由初級協(xié)調(diào)階段躍升為中級協(xié)調(diào)階段。2010—2020年中級協(xié)調(diào)區(qū)域進一步擴張，內(nèi)蒙古、浙江、福建、江西和湖北均由初級協(xié)調(diào)階段躍升為中級協(xié)調(diào)階段；良好協(xié)調(diào)區(qū)域和中度失調(diào)區(qū)域均趨于擴張，廣西和湖南分別由初級協(xié)調(diào)和中級協(xié)調(diào)階段躍升為良好協(xié)調(diào)階段，而重慶和上海分別由瀕臨失調(diào)和初級失調(diào)階段轉(zhuǎn)移為中度失調(diào)階段；瀕臨失調(diào)、輕度失調(diào)、勉強協(xié)調(diào)和初級協(xié)調(diào)區(qū)域均為收縮態(tài)勢；特別是，四川和云南由中級協(xié)調(diào)階段躍升為優(yōu)質(zhì)協(xié)調(diào)階段。這得益于，近年來我國旅游經(jīng)濟水平不斷提升，旅游經(jīng)濟可持續(xù)發(fā)展得到重視，如內(nèi)蒙古借助草地資源優(yōu)勢發(fā)展環(huán)境友好型生態(tài)旅游；特別是，在西部開發(fā)政策扶持下，四川和云南借助優(yōu)渥的自然資源和豐富的旅游資源，實現(xiàn)旅游產(chǎn)業(yè)轉(zhuǎn)型升級，旅游經(jīng)濟發(fā)展綠色化和生態(tài)化。同時，各省紛紛開展大氣污染整治行動，削減煤炭消費量，實現(xiàn)“雙碳”目標，有效遏制了生態(tài)系統(tǒng)服務下降，推動多個省份耦合協(xié)調(diào)階段實現(xiàn)更高水平躍遷，耦合協(xié)調(diào)度優(yōu)化效果明顯。然而，上海和重慶作為熱門旅游城市和常住人口超過1000萬的超大城市，人口大量集聚，旅游活動對自然系統(tǒng)改造劇烈，導致旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)等級逐漸下降。

3.4 旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)發(fā)展的主導障礙因素

本文利用障礙度模型，進一步診斷影響2000—2020年中國省際旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的主要因素。2000—2020年旅游經(jīng)濟與生態(tài)系統(tǒng)服務兩個系統(tǒng)各準則層障礙度的結(jié)構相對穩(wěn)定。在旅游經(jīng)濟系統(tǒng)中，主導障礙因素為旅游經(jīng)濟效益，占整個系統(tǒng)的41.70%。在生態(tài)系統(tǒng)服務中，障礙因素以支持服務和調(diào)節(jié)服務為主，占比分別為37.89%和34.21%。障礙度平均值大小排序為：旅游經(jīng)濟效益＞支持服務＞旅游經(jīng)濟規(guī)模＞調(diào)節(jié)服務＞旅游經(jīng)濟強度＞供給服務＞文化服務。另外，研究期內(nèi)旅游經(jīng)濟系統(tǒng)中旅游經(jīng)濟效益和旅游經(jīng)濟強度的障礙度呈上升趨勢，對兩系統(tǒng)耦合協(xié)調(diào)發(fā)展的制約作用逐步加劇。生態(tài)系統(tǒng)服務中支持服務和供給服務的障礙度也不斷增加，說明在旅游經(jīng)濟騰飛和生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展過程中，旅游發(fā)展與生態(tài)環(huán)境的競爭日益加劇，對生態(tài)系統(tǒng)支持服務和供給服務能力產(chǎn)生了嚴重的制約作用。

隨著旅游經(jīng)濟的不斷發(fā)展，旅游活動對生態(tài)環(huán)境質(zhì)量和生態(tài)平衡的干擾，導致中國省際旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)關系較為復雜、區(qū)域差異顯著。為確定具體的障礙指標，對2000—2020年各子系統(tǒng)排名前5的障礙因子進行提取（表4），發(fā)現(xiàn)在旅游經(jīng)濟系統(tǒng)中，雖然各項障礙因子的占比隨時間發(fā)生了一定波動，但歷年排名前5的因素均為旅游產(chǎn)值密度、入境旅游收入、入境旅游者人數(shù)、旅游空間密度和民航客運量，且5項指標內(nèi)部排序相對穩(wěn)定。上述5項指標制約旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的原因可能是：其一，在旅游經(jīng)濟發(fā)展過程中，雖然各省份旅游經(jīng)濟規(guī)模不斷擴大，但旅游規(guī)劃與資源管理方面仍存在一定缺陷，導致土地和生態(tài)資源的利用程度不足，旅游產(chǎn)業(yè)的增長未能與生態(tài)資源投入相匹配，未能實現(xiàn)旅游經(jīng)濟的高效綠色增長；其二，各省份旅游經(jīng)濟效益中占主要地位的是國內(nèi)旅游收入，多數(shù)旅游目的地在國際市場上的知名度不足、擁有的高產(chǎn)值，但世界級旅游經(jīng)濟要素相對匱乏，導致入境旅游市場競爭力較弱，因此，入境旅游收入和入境旅游者人數(shù)限制了旅游經(jīng)濟與生態(tài)系統(tǒng)服務間的協(xié)調(diào)發(fā)展；其三，旅游資源分布不均衡，以及旅游經(jīng)濟利益驅(qū)動旅游企業(yè)過度開發(fā)旅游資源，進一步加劇了生態(tài)環(huán)境的壓力，影響了生態(tài)系統(tǒng)的健康；其四，盡管旅游目的地的交通可達性不斷提高、整體交通網(wǎng)絡密度進一步提升，但部分省份旅游地民航通達性較低、交通網(wǎng)絡和航空基礎設施建設相對滯后，制約了旅游地的可達性和跨區(qū)域協(xié)調(diào)發(fā)展。因此，加強旅游資源規(guī)劃與管理、提升旅游經(jīng)濟增長效率、改善旅游目的地吸引力和增加航空運輸基礎設施建設等，成為未來提升旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的重要著力點。

在生態(tài)系統(tǒng)服務指標層中，每年處于主導障礙因素排名前5的均為生物多樣性、氣候調(diào)節(jié)、美學景觀、土壤保持和凈化環(huán)境，各因子在不同年份的占比變化較小，排名維持了相對的穩(wěn)定性。上述5項指標制約旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的原因可能是：其一，旅游業(yè)依賴自然生態(tài)系統(tǒng)，保持生物多樣性對于提供生態(tài)系統(tǒng)服務至關重要，生物多樣性的喪失也會導致旅游目的地形象和吸引力下降，從而對旅游經(jīng)濟產(chǎn)生負面影響；其二，氣候調(diào)節(jié)和凈化環(huán)境是生態(tài)系統(tǒng)為人類社會系統(tǒng)提供的重要服務，氣候變化、大氣污染和極端天氣事件等可能影響旅游目的地形象和舒適性，導致旅游經(jīng)濟受損；其三，美學景觀是生態(tài)系統(tǒng)文化服務中最重要的社會價值，保持和維護美學景觀對于游客體驗和目的地可持續(xù)發(fā)展至關重要；其四，旅游活動可能導致土壤侵蝕和破壞，對生態(tài)系統(tǒng)造成負面影響，從而導致旅游經(jīng)濟對生態(tài)系統(tǒng)服務脅迫效應的增強。因此，需要加強對森林、水域和濕地等具有較高休閑文化服務和支持服務的土地生態(tài)系統(tǒng)的保護和發(fā)展，注重發(fā)揮森林資源、天然濕地等優(yōu)質(zhì)的氣候調(diào)節(jié)、土壤保持和凈化環(huán)境功能，筑牢旅游經(jīng)濟綠色發(fā)展的生態(tài)本底。

4 結(jié)論與討論

4.1 結(jié)論

在生態(tài)文明建設背景下，只有樹立尊重自然生態(tài)環(huán)境、堅持綠色發(fā)展的理念，才能實現(xiàn)旅游業(yè)可持續(xù)發(fā)展。本文將旅游經(jīng)濟和生態(tài)系統(tǒng)服務納入同一分析框架，定量測度2000—2020年中國省域旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)度，科學度量其空間演化格局，識別影響其發(fā)展的主要障礙因素，為探尋旅游業(yè)綠色發(fā)展的最優(yōu)路徑提供科學依據(jù)和實證支撐。主要研究結(jié)論如下。

第一，研究期內(nèi)，中國旅游經(jīng)濟綜合指數(shù)總體呈現(xiàn)波動上升態(tài)勢，旅游經(jīng)濟由2000年的0.043增長到2020年的0.162，累計增長12%，表明中國旅游經(jīng)濟不斷穩(wěn)步提升，旅游發(fā)展成效顯著。此外，研究期內(nèi)旅游經(jīng)濟總體呈東南高、西北低的空間分布格局。

第二，研究期內(nèi)，中國生態(tài)系統(tǒng)服務有一定的改善，整體呈現(xiàn)增長態(tài)勢，生態(tài)系統(tǒng)服務價值由2000年的290 669.3億元上升至2020年的303 179.1億元，年均增長625.5億元，生態(tài)系統(tǒng)服務價值結(jié)構較為穩(wěn)定，表明我國的生態(tài)保護措施卓有成效，人地關系正逐步改善。從空間分布來看，研究期內(nèi)生態(tài)系統(tǒng)服務價值呈東高西低、由東南向西北遞減的空間分布格局。

第三，研究期內(nèi)，中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)效應穩(wěn)步增強，且兩大系統(tǒng)長期處于高水平磨合階段。從時序維度來看，旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展經(jīng)歷了瀕臨失調(diào)階段（2000—2015）→勉強協(xié)調(diào)階段（2020）的演化軌跡。從空間分布來看，旅游經(jīng)濟與生態(tài)系統(tǒng)服務耦合協(xié)調(diào)度具有東高西低、南高北低的分布特征，極核特征明顯，表明中國旅游經(jīng)濟與生態(tài)系統(tǒng)服務存在正向協(xié)調(diào)關系。

第四，研究期內(nèi)，旅游經(jīng)濟效益、支持服務、旅游經(jīng)濟規(guī)模和調(diào)節(jié)服務是制約旅游經(jīng)濟和生態(tài)系統(tǒng)服務耦合協(xié)調(diào)發(fā)展的主要因素。為避免旅游經(jīng)濟快速發(fā)展對生物多樣性、氣候調(diào)節(jié)、美學景觀、土壤保持和凈化環(huán)境等生態(tài)系統(tǒng)服務功能產(chǎn)生脅迫作用，推進旅游經(jīng)濟綠色發(fā)展勢在必行。在旅游經(jīng)濟發(fā)展中，應進一步提升旅游目的地的交通可達性，合理規(guī)劃土地等生態(tài)資源，提高旅游經(jīng)濟發(fā)展質(zhì)量，加強對森林資源、天然濕地等具有較高休閑文化服務和支持服務的土地生態(tài)系統(tǒng)的保護和發(fā)展，推動耦合協(xié)調(diào)關系的躍升。

4.2 討論

本文關于中國旅游經(jīng)濟發(fā)展水平的時空演變結(jié)論與以往研究[21， 68， 72]一致，符合實際發(fā)展規(guī)律。相較于已有旅游與生態(tài)環(huán)境協(xié)調(diào)發(fā)展研究，本文創(chuàng)新性地引入生態(tài)系統(tǒng)服務框架對生態(tài)環(huán)境質(zhì)量進行測算，獲得了一個更全面、更準確的度量結(jié)果，為相關研究提供了參考。若未來研究能廣泛使用這一公認的、較為系統(tǒng)的測算方式，那么以微觀和中觀為主要尺度的相關研究，就可以直接對研究區(qū)域間的研究結(jié)果進行對比，有效地提高領域的研究效率。此外，生態(tài)系統(tǒng)服務測算是以經(jīng)濟學市場價格為表現(xiàn)形式的，這種方式能夠喚醒公眾的環(huán)保意識，也為政府部門的決策工作提供了一個有效的工具。本文對當前生態(tài)系統(tǒng)服務測算進行了動態(tài)化改進，綜合考慮了現(xiàn)實中區(qū)域社會經(jīng)濟發(fā)展水平和資源條件的異質(zhì)性，從方法論層面創(chuàng)新了生態(tài)系統(tǒng)服務測算方式。

在此基礎上，本文對旅游經(jīng)濟與生態(tài)系統(tǒng)服務關系展開耦合協(xié)調(diào)分析，識別影響旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展的主要障礙因素，具有重要的理論和實踐意義。在理論上，可以豐富新時代旅游業(yè)綠色可持續(xù)發(fā)展的研究文獻，將旅游經(jīng)濟與生態(tài)系統(tǒng)服務有機整合，詮釋“兩山理論”的發(fā)展觀，反哺生態(tài)文明建設。在實踐上，可以指導旅游發(fā)展綠色化和生態(tài)化，實現(xiàn)旅游經(jīng)濟增長方式的轉(zhuǎn)變和發(fā)展理念的超越，探索旅游經(jīng)濟綠色可持續(xù)的發(fā)展路徑。然而，本文作為旅游經(jīng)濟與生態(tài)系統(tǒng)服務協(xié)調(diào)發(fā)展研究的初步探索，存在一些局限性。例如囿于生態(tài)系統(tǒng)文化服務數(shù)據(jù)可獲取性較低，本文僅選取了美學景觀作為其測量指標，后續(xù)研究將結(jié)合陳述偏好法或顯示偏好法[35]，豐富文化服務的評估數(shù)據(jù)，進一步豐富這一動態(tài)生態(tài)系統(tǒng)服務價值分析框架的穩(wěn)健性與適用性。

參考文獻（References）

[1] 唐葉萍. 論人與自然和諧發(fā)展的生態(tài)倫理意蘊[J]. 道德與文明， 2009（3）： 100-102. [TANG Yeping. On the ecological ethics implication of harmonious development between man and nature[J]. Morality and Civilization， 2009（3）： 100-102.]

[2] 童昀， 馬勇， 劉軍. 入境旅游提升了城市綠色全要素生產(chǎn)率嗎？——基于非線性視角[J]. 旅游學刊， 2021， 36（9）： 120-133. [TONG Yun， MA Yong， LIU Jun. Does the inbound tourism increase the city’s green total factor productivity？ —Based on nonlinear perspective[J]. Tourism Tribune， 2021， 36（9）： 120-133.]

[3] WANG Q， MAO Z， XIAN L， et al. A study on the coupling coordination between tourism and the low-carbon city[J]. Asia Pacific Journal of Tourism Research， 2019， 24（6）： 550-562.

[4] TAN S， YANG J， YAN J， et al. A holistic low carbon city indicator framework for sustainable development[J]. Applied Energy， 2017， 185： 1919-1930.

[5] SHI C， PENG J. Construction of low-carbon tourist attractions based on low-carbon economy[J]. Energy Procedia， 2011， 5： 759-762.

[6] 馬勇， 郭田田. 踐行“兩山理論”生態(tài)旅游發(fā)展的核心價值與實施路徑[J]. 旅游學刊， 2018， 33（8）： 16-18. [MA Yong， GUO Tianian. Practicing the “two mountains theory”： The core value and implementation path of ecotourism development[J]. Tourism Tribune， 2018， 33（8）： 16-18.]

[7] EDWARDS J R. The UK heritage coasts： An assessment of the ecological impacts of tourism[J]. Annals of Tourism Research， 1987， 14（1）： 71-87.

[8] MORRISON P， SELMAN P. Tourism and the environment： A case study from Turkey[J]. Environmentalist， 1991， 11（2）： 113-129.

[9] COLE D N， MONZ C A. Impacts of camping on vegetation： Response and recovery following acute and chronic disturbance[J]. Environmental Management， 2003， 32（6）： 693-705.

[10] NAVARRO J E， TEJADA T M， ALMEIDA G F， et al. Carrying capacity assessment for tourist destinations. Methodology for the creation of synthetic indicators applied in a coastal area[J]. Tourism Management， 2012， 33（6）： 1337-1346.

[11] 劉佳， 于水仙， 王佳. 濱海旅游環(huán)境承載力評價與量化測度研究——以山東半島藍色經(jīng)濟區(qū)為例[J]. 中國人口·資源與環(huán)境， 2012， 22（9）： 163-170. [LIU Jia， YU Shuixian， WANG Jia. Evaluation and quantitative measurement on coastal tourism environmental carrying capacity： A case study on the Blue Shandong Peninsula Economic Zone[J]. China Population， Resources and Environment， 2012（9）： 163-170.]

[12] 馬駿. 基于生態(tài)環(huán)境閾限與旅游承載力背景下生物多樣性保護策略研究——以世界自然遺產(chǎn)武陵源核心景區(qū)為例[J]. 經(jīng)濟地理， 2016， 36（4）： 195-202. [MA Jun. Biodiversity protection strategy under the background of ecological environment threshold and the tourism carrying capacity—According to the core of the World Natural Heritage Wulingyuan Scenic Area as an example[J]. Economic Geography， 2016， 36（4）： 195-202.]

[13] HEDLUND T. The impact of values， environmental concern， and willingness to accept economic sacrifices to protect the environment on tourists’ intentions to buy ecologically sustainable tourism alternatives[J]. Tourism and Hospitality Research， 2011， 11（4）： 278-288.

[14] ARROBAS F， FERREIRA J， BRITO-HENRIQUES E， et al. Measuring tourism and environmental sciences students’ attitudes towards sustainable tourism[J]. Journal of Hospitality， Leisure， Sport amp; Tourism Education， 2020， 27： 100273.

[15] 翁鋼民， 潘越， 李凌雁. 基于改進DPSIR-DS模型的旅游生態(tài)安全等級測度及時空演變分析——以“絲綢之路”沿線五省區(qū)為例[J]. 旅游科學， 2018， 32（6）： 17-32. [WENG Gangmin， PAN Yue， LI Lingyan. The eco-security grading and spatial-temporal evolution of tourism based on improved DPSIR-DS model： A case study of the five northwestern provinces along the Silk Road[J]. Tourism Science， 2018， 32（6）： 17-32.]

[16] 穆學青， 郭向陽， 明慶忠， 等. 黃河流域旅游生態(tài)安全的動態(tài)演變特征及驅(qū)動因素[J]. 地理學報， 2022， 77（3）： 714-735. [MU Xueqing， GUO Xiangyang， MING Qingzhong， et al. Dynamic evolution characteristics and driving factors of tourism ecological security in the Yellow River Basin[J]. Acta Geographica Sinica， 2022， 77（3）： 714-735.]

[17] 袁宏瑞， 王群. 旅游與生態(tài)共生演進模式與生態(tài)安全判定——以國家重點生態(tài)功能區(qū)安徽省黃山區(qū)為例[J]. 旅游學刊， 2022， 37（12）： 25-38. [YUAN Hongrui， WANG Qun. The evolution model of the tourism-ecology symbiosis and the determination of ecological security—A case study of Huangshan District， Anhui province， a key ecological functional zone of China[J].Tourism Tribune， 2022， 37（12）： 25-38.]

[18] 郭向陽， 穆學青， 丁正山， 等. 城市生態(tài)環(huán)境與旅游經(jīng)濟協(xié)調(diào)效應及動態(tài)關系——以曲靖為例[J]. 經(jīng)濟地理， 2020， 40（7）： 231-240. [GUO Xiangyang， MU Xueqing， DING Zhengshan， et al. Coordination effect and dynamic relationship of urban ecological environment and tourism economy： A case study of Qujing[J]. Economic Geography， 2020， 40（7）： 231-240.]

[19] 馬勇， 李麗霞， 任潔. 神農(nóng)架林區(qū)旅游經(jīng)濟-交通狀況-生態(tài)環(huán)境協(xié)調(diào)發(fā)展研究[J]. 經(jīng)濟地理， 2017， 37（10）： 215-220; 227. [MA Yong， LI Lixia， REN Jie. Coordination development research among the tourism economy-traffic condition-ecological environment in Shengnongjia forest district[J]. Economic Geography， 2017， 37（10）： 215-220; 227.]

[20] FEI J， LIN Y， JIANG Q， et al. Spatiotemporal coupling coordination measurement on islands’ economy-environment-tourism system[J]. Ocean amp; Coastal Management， 2021， 212： 105793.

[21] 王兆峰， 劉慶芳. 長江經(jīng)濟帶旅游生態(tài)效率時空演變及其與旅游經(jīng)濟互動響應[J]. 自然資源學報， 2019， 34（9）： 1945-1961. [WANG Zhaofeng， LIU Qingfang. The spatio-temporal evolution of tourism eco-efficiency in the Yangtze River Economic Belt and its interactive response with tourism economy[J]. Journal of Natural Resources， 2019， 34（9）： 1945-1961.]

[22] 廖重斌. 環(huán)境與經(jīng)濟協(xié)調(diào)發(fā)展的定量評判及其分類體系——以珠江三角洲城市群為例[J]. 熱帶地理， 1999， 19（2）： 171-177. [LIAO Chongbin. Quantitaitve judgement and classification system for coordinated development of environment and economy—A case study of the city group in the Pearl River Delta[J]. Tropical Geography， 1999， 19（2）： 171-177.]

[23] LIU Y， SUK S， CAI Y. Spatial and temporal changes in the coupling of ecological environment and tourism development： The case of Kyushu， Japan[J]. Environmental Research Letters， 2022， 18： 014004.

[24] CHAPLIN-KRAMER R， SHARP R P， WEIL C， et al. Global modeling of nature’s contributions to people[J]. Science， 2019， 366（6462）： 255-258.

[25] DRIUS M， BONGIORNI L， DEPELLEGRIN D， et al. Tackling challenges for Mediterranean sustainable coastal tourism： An ecosystem service perspective[J]. Science of the Total Environment， 2019， 652： 1302-1317.

[26] CARPENTER S R， MOONEY H A， AGARD J， et al. Science for managing ecosystem services： Beyond the Millennium Ecosystem Assessment[J]. Proceedings of the National Academy of Sciences， 2009， 106（5）： 1305-1312.

[27] 黃木易， 岳文澤， 方斌， 等. 1970—2015年大別山區(qū)生態(tài)服務價值尺度響應特征及地理探測機制[J]. 地理學報， 2019， 74（9）： 1904-1920. [HUANG Muyi， YUE Wenze， FANG Bin， et al. Scale response characteristics and geographic exploration mechanism of spatial differentiation of ecosystem service values in Dabie mountain area， central China from 1970 to 2015[J]. Acta Geographica Sinica， 2019， 74（9）： 1904-1920.]

[28] BAKER J， SHEATE W R， PHILLIPS P， et al. Ecosystem services in environmental assessment—Help or hindrance？[J]. Environmental Impact Assessment Review， 2013， 40： 3-13.

[29] KARJALAINEN T P， MARTTUNEN M， SARKKI S， et al. Integrating ecosystem services into environmental impact assessment： An analytic-deliberative approach[J]. Environmental Impact Assessment Review， 2013， 40： 54-64.

[30] HONRADO J P， VIEIRA C， SOARES C， et al. Can we infer about ecosystem services from EIA and SEA practice？ A framework for analysis and examples from Portugal[J]. Environmental Impact Assessment Review， 2013， 40： 14-24.

[31] COSTANZA R， DEGROOT R， SUTTON P， et al. Changes in the global value of ecosystem services[J]. Global Environmental Change， 2014， 26： 152-158.

[32] KUMAR P， ESEN S E， YASHIRO M. Linking ecosystem services to strategic environmental assessment in development policies[J]. Environmental Impact Assessment Review， 2013， 40： 75-81.

[33] CHEN H. Land use trade-offs associated with protected areas in China： Current state， existing evaluation methods， and future application of ecosystem service valuation[J]. Science of the Total Environment， 2020， 711： 134688.

[34] CHEN H. Complementing conventional environmental impact assessments of tourism with ecosystem service valuation： A case study of the Wulingyuan scenic area， China[J]. Ecosystem Services， 2020， 43： 101100.

[35] 陳東軍， 鐘林生. 生態(tài)系統(tǒng)服務價值評估與實現(xiàn)機制研究綜述[J]. 中國農(nóng)業(yè)資源與區(qū)劃， 2023， 44（1）： 84-94. [CHEN Dongjun， ZHONG Linsheng. Review of the value evaluation and realization mechanism of ecosystem services[J]. Chinese Journal of Agricultural Resources and Regional Planning， 2023， 44（1）： 84-94.]

[36] 杜金鴻， 劉方正， 周越， 等. 自然保護地生態(tài)系統(tǒng)服務價值評估研究進展[J]. 環(huán)境科學研究， 2019， 32（9）： 1475-1482. [DU Jinhong， LIU Fangzheng， ZHOU Yue， et al. A review of ecosystem services assessment and valuation of protected areas[J]. Research of Environmental Sciences， 2019， 32（9）： 1475-1482.]

[37] 趙海蘭. 生態(tài)系統(tǒng)服務分類與價值評估研究進展[J]. 生態(tài)經(jīng)濟， 2015， 31（8）： 27-33. [ZHAO Hailan. Research progress of classification and value evaluation of ecosystem services[J]. Ecological Economy， 2015， 31（8）： 27-33.]

[38] COSTANZA R， D’ARGE R， DEGROOT R， et al. The value of the world’s ecosystem services and natural capital[J]. Nature， 1997， 387（6630）： 253-260.

[39] 謝高地， 魯春霞， 冷允法， 等. 青藏高原生態(tài)資產(chǎn)的價值評估[J]. 自然資源學報， 2003（2）： 189-196. [XIE Gaodi， LU Chunxia， LENG Yunfa， et al. Ecological assets valuation of the Tibetan Plateau[J]. Journal of Natural Resources， 2003， 18（2）： 189-196.]

[40] 謝高地， 甄霖， 魯春霞， 等. 一個基于專家知識的生態(tài)系統(tǒng)服務價值化方法[J]. 自然資源學報， 2008（5）： 911-919. [XIE Gaodi， ZHEN Lin， LU Chunxia， et al. Expert knowledge based valuation method of ecosystem services in China[J]. Journal of Natural Resources， 2008， 23（5）： 911-919.]

[41] BENNETT E M， PETERSON G D， GORDON L J. Understanding relationships among multiple ecosystem services[J]. Ecology Letters， 2009， 12（12）： 1394-1404.

[42] BATEMAN I J， HARWOOD A R， MACE G M， et al. Bringing ecosystem services into economic decision-making： Land use in the United Kingdom[J]. Science， 2013， 341（6141）： 45-50.

[43] XUE M， LUO Y. Dynamic variations in ecosystem service value and sustainability of urban system： A case study for Tianjin City， China[J]. Cities， 2015， 46： 85-93.

[44] ZANG Z， ZOU X， ZUO P， et al. Impact of landscape patterns on ecological vulnerability and ecosystem service values： An empirical analysis of Yancheng Nature Reserve in China[J]. Ecological Indicators， 2017， 72： 142-152.

[45] LI F， ZHANG S， YANG J， et al. Effects of land use change on ecosystem services value in west Jilin since the reform and opening of China[J]. Ecosystem Services， 2018， 31： 12-20.

[46] XING L， XUE M， WANG X. Spatial correction of ecosystem service value and the evaluation of eco-efficiency： A case for China’s provincial level[J]. Ecological Indicators， 2018， 95： 841-850.

[47] CHURCH A， COLES T， FISH R. Tourism in sub-global assessments of ecosystem services[J]. Journal of Sustainable Tourism， 2017， 25（11）： 1529-1546.

[48] DONáZAR J A， CEBALLOS O， CORTéS-AVIZANDA A. Tourism in protected areas： Disentangling road and traffic effects on intra-guild scavenging processes[J]. Science of the Total Environment， 2018， 630： 600-608.

[49] 朱順順， 章錦河， 胡歡， 等. 旅游發(fā)展的生態(tài)系統(tǒng)服務價值增值效應研究[J]. 長江流域資源與環(huán)境， 2019， 28（3）： 603-613. [ZHU Shunshun， ZHANG Jinhe， HU Huan， et al. Study on the value-added effect of ecosystem service values in tourism development[J]. Resources and Environment in the Yangtze Basin， 2019， 28（3）： 603-613.]

[50] MA M， TANG J. Interactive coercive relationship and spatio-temporal coupling coordination degree between tourism urbanization and eco-environment： A case study in Western China[J]. Ecological Indicators， 2022， 142： 109149.

[51] GUAN D， GAO W， SU W， et al. Modeling and dynamic assessment of urban economy-resource-environment system with a coupled system dynamics-geographic information system model[J]. Ecological Indicators， 2011， 11（5）： 1333-1344.

[52] TANG Z. An integrated approach to evaluating the coupling coordination between tourism and the environment[J]. Tourism Management， 2015， 46： 11-19.

[53] SHAHEEN K， ZAMAN K， BATOOL R， et al. Dynamic linkages between tourism， energy， environment， and economic growth： Evidence from top 10 tourism-induced countries[J]. Environmental Science and Pollution Research， 2019， 26： 31273-31283.

[54] RIENSCHE M， CASTILLO A， FLORES-DíAZ A， et al. Tourism at Costalegre， Mexico： An ecosystem services-based exploration of current challenges and alternative futures[J]. Futures， 2015， 66： 70-84.

[55] WANG P， WANG J， ZHANG J， et al. Spatial-temporal changes in ecosystem services and social-ecological drivers in a typical coastal tourism city： A case study of Sanya， China[J]. Ecological Indicators， 2022， 145： 109607.

[56] LI J， BAI Y， ALATALO J M. Impacts of rural tourism-driven land use change on ecosystems services provision in Erhai Lake Basin， China[J]. Ecosystem Services， 2020， 42： 101081.

[57] MENDOZA-GONZáLEZ G， MARTíNEZ M L， LITHGOW D， et al. Land use change and its effects on the value of ecosystem services along the coast of the Gulf of Mexico[J]. Ecological Economics， 2012， 82： 23-32.

[58] LIU K， QIAO Y， SHI T， et al. Study on coupling coordination and spatiotemporal heterogeneity between economic development and ecological environment of cities along the Yellow River Basin[J]. Environmental Science and Pollution Research， 2021， 28： 6898-6912.

[59] GENG Y， WANG R， WEI Z， et al. Temporal-spatial measurement and prediction between air environment and inbound tourism： Case of China[J]. Journal of Cleaner Production， 2021， 287： 125486.

[60] GHERMANDI A， NUNES P A. A global map of coastal recreation values： Results from a spatially explicit meta-analysis[J]. Ecological Economics， 2013， 86： 1-15.

[61] SCHILLING M， CHIANG L. The effect of natural resources on a sustainable development policy： The approach of non-sustainable externalities[J]. Energy Policy， 2011， 39（2）： 990-998.

[62] LI T， CUI Y， LIU A. Spatiotemporal dynamic analysis of forest ecosystem services using “big data”： A case study of Anhui province， central-eastern China[J]. Journal of Cleaner Production， 2017， 142： 589-599.

[63] XIE G， ZHANG C， ZHEN L， et al. Dynamic changes in the value of China’s ecosystem services[J]. Ecosystem Services， 2017， 26： 146-154.

[64] CHUAI X， HUANG X， WU C， et al. Land use and ecosystems services value changes and ecological land management in coastal Jiangsu， China[J]. Habitat International， 2016， 57： 164-174.

[65] 董家華， 舒廷飛， 謝慧， 等. 城市建設用地生態(tài)服務功能價值計算與應用[J]. 同濟大學學報（自然科學版）， 2007， 35（5）： 636-640. [DONG Jiahua， SHU Tingfei， Xie Hui， et al. Calculative method for ecosystem services values of urban constructive lands and its application[J]. Journal of Tongji University （Natural Science Edition）， 2007， 35 （5）： 636-640.]

[66] HU X， WU C， HONG W， et al. Impact of land-use change on ecosystem service values and their effects under different intervention scenarios in Fuzhou city， China[J]. Geosciences Journal， 2013， 17： 497-504.

[67] 魯春陽， 文楓， 楊慶媛， 等. 基于改進TOPSIS法的城市土地利用績效評價及障礙因子診斷——以重慶市為例[J]. 資源科學， 2011， 33（3）： 535-541. [LU Chunyang， WEN Feng， YANG Qingyuan， et al. An evaluation of urban land use performance based on the improved TOPSIS method and diagnosis of its obstacle indicators： A case study of Chongqing[J]. Resources Science， 2011， 33（3）： 535-541.]

[68] 王兆峰， 陳青青. 長江經(jīng)濟帶旅游產(chǎn)業(yè)與生態(tài)環(huán)境交互脅迫關系驗證及協(xié)調(diào)效應研究[J]. 長江流域資源與環(huán)境， 2021， 30（11）： 2581-2593. [WANG Zhaofeng， CHEN Qingqing. Study on the verification and coordination effect of the interactive stress of tourism industry and ecological environment in the Yangtze River Economic Belt[J]. Resources and Environment in the Yangtze Basin， 2021， 30（11）： 2581-2593.]

[69] 王勝鵬， 喬花芳， 馮娟， 等. 黃河流域旅游生態(tài)效率時空演化及其與旅游經(jīng)濟互動響應[J]. 經(jīng)濟地理， 2020， 40（5）： 81-89. [WANG Shengpeng， QIAO Huafang， FENG Juan， et al. The spatio-temporal evolution of tourism eco-efficiency in the Yellow River Basin and its interactive response with tourism economy development level[J]. Economic Geography， 2020， 40（5）： 81-89.]

[70] 劉長生， 陳昀， 簡玉峰， 等. 中國旅游產(chǎn)業(yè)發(fā)展間接就業(yè)帶動能力測算及其時空差異[J]. 地理學報， 2022， 77（4）： 918-935. [LIU Changsheng， CHEN Yun， JIAN Yufeng， et al. Measurement of the driving capacity of tourism industry on indirect employment and its spatio-temporal differences in China[J]. Acta Geographica Sinica， 2022， 77（4）： 918-935.]

[71] 李維維， 陳田， 馬曉龍. 中國旅游化與城市化耦合協(xié)調(diào)關系的省際格局及形成機制分析[J]. 旅游科學， 2018， 32（5）： 14-29; 80. [LI Weiwei， CHEN Tian， MA Xiaolong. Analysis of the interprovincial pattern of the coupled and coordinated relationship and mechanism between touristization and urbanization in China[J]. Tourism Science， 2018， 32（5）： 14-29; 80.]

[72] 王兆峰， 李靜怡. 黃河流域旅游發(fā)展與生態(tài)環(huán)境耦合協(xié)調(diào)時空演變及交互脅迫關系驗證[J]. 長江流域資源與環(huán)境， 2022， 31（2）： 447-460. [WANG Zhaofeng， LI Jingyi. Verify and study the coupling coordination development and the interactive stress between tourism and eco-environment in the Yellow River Basin[J]. Resources and Environment in the Yangtze Basin， 2022， 31（2）： 447-460.]

[73] 張中浩， 聶甜甜， 高陽， 等. 長江經(jīng)濟帶生態(tài)系統(tǒng)服務與經(jīng)濟社會發(fā)展耦合協(xié)調(diào)關聯(lián)時空特征研究[J].長江流域資源與環(huán)境， 2022， 31（5）： 1086-1100. [ZHANG Zhonghao， NIE Tiantian， GAO Yang， et al. Study on temporal and spatial characteristics of coupling coordination correlation between ecosystem services and economic-social development in the Yangtze River Economic Belt[J]. Resources and Environment in the Yangtze Basin， 2022， 31（5）： 1086-1100.]

[74] 劉志濤， 王少劍， 方創(chuàng)琳. 粵港澳大灣區(qū)生態(tài)系統(tǒng)服務價值的時空演化及其影響機制[J]. 地理學報， 2021， 76（11）： 2797-2813. [LIU Zhitao， WANG Shaojian， FANG Chuanglin. Spatiotemporal evolution and influencing mechanism of ecosystem service value in the Guangdong-Hong Kong-Macao Greater Bay Area[J]. Acta Geographica Sinica， 2021， 76（11）： 2797-2813.]

[75] 朱治州， 鐘業(yè)喜. 長江三角洲城市群土地利用及其生態(tài)系統(tǒng)服務價值時空演變研究[J]. 長江流域資源與環(huán)境， 2019， 28（7）： 1520-1530. [ZHU Zhizhou， ZHONG Yexi. Spatio-temporal evolution of land use and ecosystem service value in Yangtze River Delta urban agglomeration[J]. Resources and Environment in the Yangtze Basin， 2019， 28（7）： 1520-1530.]