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

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

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

［中圖分類號］F59

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

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

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

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

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

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

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

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

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

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

2.1.1" " 指標(biāo)體系構(gòu)建

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

2.1.2" " 測算方法

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

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

2.2.1" " 標(biāo)準(zhǔn)單位當(dāng)量因子價(jià)值的確定

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

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

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

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

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

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

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

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

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

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

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

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

[D=l′slC] （4）

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

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

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

[S=logplogP] （5）

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

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

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

基于當(dāng)量因子法和上述調(diào)整系數(shù)[38， 63]，構(gòu)建非建設(shè)用地生態(tài)系統(tǒng)服務(wù)價(jià)值動態(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]項(xiàng)生態(tài)系統(tǒng)服務(wù)功能價(jià)值系數(shù)；[E]為標(biāo)準(zhǔn)單位價(jià)值系數(shù)（元/hm2）；[D]為研究區(qū)域社會經(jīng)濟(jì)調(diào)整系數(shù)；[S]為資源稀缺調(diào)整系數(shù)；[Ekf]為第[k]種土地利用類型的第[f]項(xiàng)服務(wù)功能的當(dāng)量因子。式（7）中，[ESVk]為第[k]種土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值；[Ak]是第[k]種土地利用類型的面積（hm2）。式（8）中，[ESVf]為第[f]項(xiàng)生態(tài)系統(tǒng)服務(wù)功能的生態(tài)系統(tǒng)服務(wù)價(jià)值。式（9）中，[ESV]為生態(tài)系統(tǒng)服務(wù)總價(jià)值。

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

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

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

[Vcj=Ecj×E] （11）

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

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

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

[Vws=-W×PwA] （12）

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3 結(jié)果分析

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

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

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

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

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

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

將研究時(shí)序內(nèi)各省份生態(tài)系統(tǒng)服務(wù)價(jià)值分為5個(gè)等級，繪制生態(tài)系統(tǒng)服務(wù)價(jià)值空間分布情況（圖3）。從空間分布來看，總體呈東高西低、由東南向西北遞減的分布格局，高值區(qū)域主要分布于我國的東北和西南地區(qū)，不同時(shí)期生態(tài)系統(tǒng)服務(wù)價(jià)值存在明顯的空間分異（圖3）。原因可能源于：其一，林地生態(tài)價(jià)值當(dāng)量高于耕地、草地和水體等[75]，我國東北地區(qū)和西南地區(qū)多為山地丘陵地形，森林覆蓋率較高，如大小興安嶺是我國重要的生態(tài)屏障，因此生態(tài)系統(tǒng)服務(wù)價(jià)值較高；其二，隨著人口遷移、土地利用變化、城市擴(kuò)展和經(jīng)濟(jì)發(fā)展等人類活動增加，導(dǎo)致生態(tài)系統(tǒng)服務(wù)發(fā)生變化且空間差異顯著。總體來看，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)濟(jì)與生態(tài)系統(tǒng)服務(wù)耦合協(xié)調(diào)關(guān)系

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

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

其次，從耦合協(xié)調(diào)度指數(shù)（D）來看，2000—2020年旅游經(jīng)濟(jì)與生態(tài)系統(tǒng)服務(wù)的耦合協(xié)調(diào)度呈穩(wěn)步上升態(tài)勢，從0.54上升至0.60，耦合協(xié)調(diào)等級由瀕臨失調(diào)躍升為勉強(qiáng)協(xié)調(diào)階段，但研究時(shí)序內(nèi)二者整體的耦合協(xié)調(diào)度水平不高，系統(tǒng)間協(xié)調(diào)發(fā)展仍待進(jìn)一步增強(qiáng)。可能的原因是：我國旅游經(jīng)濟(jì)發(fā)展水平還有較大的提升空間，旅游經(jīng)濟(jì)的快速發(fā)展在一定程度上以犧牲目的地生態(tài)系統(tǒng)服務(wù)效益為代價(jià)，在發(fā)展過程中缺乏科學(xué)合理規(guī)劃，導(dǎo)致系統(tǒng)間協(xié)調(diào)發(fā)展水平不高。分區(qū)域來看，中部地區(qū)耦合協(xié)調(diào)度高于全國平均水平，實(shí)現(xiàn)了初級協(xié)調(diào)向中級協(xié)調(diào)階段的躍遷，東部地區(qū)趨近于全國平均水平，而西部地區(qū)距離全國平均水平仍有較大差距。由此可見，隨著旅游經(jīng)濟(jì)和生態(tài)系統(tǒng)服務(wù)綜合水平的不斷提高，二者實(shí)現(xiàn)由無序發(fā)展向接近有序發(fā)展的轉(zhuǎn)變，這與“綠色發(fā)展”“美麗中國”“生態(tài)文明建設(shè)”等政策指導(dǎo)息息相關(guān)，特別是綠色可持續(xù)發(fā)展理念推動旅游業(yè)的結(jié)構(gòu)性改革，旅游經(jīng)濟(jì)高質(zhì)量發(fā)展初有成效，旅游經(jīng)濟(jì)對生態(tài)系統(tǒng)服務(wù)的脅迫效應(yīng)有所下降，旅游經(jīng)濟(jì)與生態(tài)系統(tǒng)服務(wù)的協(xié)同效應(yīng)趨于良好，系統(tǒng)間耦合協(xié)調(diào)發(fā)展有所改善。中部地區(qū)雖然旅游業(yè)起步較晚，旅游經(jīng)濟(jì)綜合水平相對較低，但基于良好自然生態(tài)狀況的生態(tài)旅游、康養(yǎng)旅游等新型旅游形式的發(fā)展，以綠色發(fā)展為契機(jī)促進(jìn)了生態(tài)系統(tǒng)服務(wù)水平的不斷增強(qiáng)，加快了系統(tǒng)間協(xié)調(diào)發(fā)展，實(shí)現(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個(gè)等級，從縱向維度比較分析旅游經(jīng)濟(jì)與生態(tài)系統(tǒng)服務(wù)耦合協(xié)調(diào)關(guān)系的空間分異和演化特征（圖4）。結(jié)果顯示，中國旅游經(jīng)濟(jì)與生態(tài)系統(tǒng)服務(wù)耦合協(xié)調(diào)度呈現(xiàn)“東高西低、南高北低”的分布特征，極核特征明顯。總體來看，2000—2020年高值區(qū)主要集中于珠江流域和長江流域省份，并向區(qū)域外緣延伸；低值區(qū)則多分布于西北地區(qū)；高值區(qū)趨于擴(kuò)張，低值區(qū)趨于收縮，空間分布格局變化較大。

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

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

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

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

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

4 結(jié)論與討論

4.1 結(jié)論

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

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

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

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

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

4.2 討論

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

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

參考文獻(xiàn)（References）

[1] 唐葉萍. 論人與自然和諧發(fā)展的生態(tài)倫理意蘊(yùn)[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]. 旅游學(xué)刊， 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ā)展的核心價(jià)值與實(shí)施路徑[J]. 旅游學(xué)刊， 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)境承載力評價(jià)與量化測度研究——以山東半島藍(lán)色經(jīng)濟(jì)區(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)境閾限與旅游承載力背景下生物多樣性保護(hù)策略研究——以世界自然遺產(chǎn)武陵源核心景區(qū)為例[J]. 經(jīng)濟(jì)地理， 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] 翁鋼民， 潘越， 李凌雁. 基于改進(jìn)DPSIR-DS模型的旅游生態(tài)安全等級測度及時(shí)空演變分析——以“絲綢之路”沿線五省區(qū)為例[J]. 旅游科學(xué)， 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] 穆學(xué)青， 郭向陽， 明慶忠， 等. 黃河流域旅游生態(tài)安全的動態(tài)演變特征及驅(qū)動因素[J]. 地理學(xué)報(bào)， 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)共生演進(jìn)模式與生態(tài)安全判定——以國家重點(diǎn)生態(tài)功能區(qū)安徽省黃山區(qū)為例[J]. 旅游學(xué)刊， 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] 郭向陽， 穆學(xué)青， 丁正山， 等. 城市生態(tài)環(huán)境與旅游經(jīng)濟(jì)協(xié)調(diào)效應(yīng)及動態(tài)關(guān)系——以曲靖為例[J]. 經(jīng)濟(jì)地理， 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)濟(jì)-交通狀況-生態(tài)環(huán)境協(xié)調(diào)發(fā)展研究[J]. 經(jīng)濟(jì)地理， 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)濟(jì)帶旅游生態(tài)效率時(shí)空演變及其與旅游經(jīng)濟(jì)互動響應(yīng)[J]. 自然資源學(xué)報(bào)， 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)濟(jì)協(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)服務(wù)價(jià)值尺度響應(yīng)特征及地理探測機(jī)制[J]. 地理學(xué)報(bào)， 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)服務(wù)價(jià)值評估與實(shí)現(xiàn)機(jī)制研究綜述[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] 杜金鴻， 劉方正， 周越， 等. 自然保護(hù)地生態(tài)系統(tǒng)服務(wù)價(jià)值評估研究進(jìn)展[J]. 環(huán)境科學(xué)研究， 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)服務(wù)分類與價(jià)值評估研究進(jìn)展[J]. 生態(tài)經(jīng)濟(jì)， 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)的價(jià)值評估[J]. 自然資源學(xué)報(bào)， 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] 謝高地， 甄霖， 魯春霞， 等. 一個(gè)基于專家知識的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法[J]. 自然資源學(xué)報(bào)， 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)服務(wù)價(jià)值增值效應(yī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] 董家華， 舒廷飛， 謝慧， 等. 城市建設(shè)用地生態(tài)服務(wù)功能價(jià)值計(jì)算與應(yīng)用[J]. 同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 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] 魯春陽， 文楓， 楊慶媛， 等. 基于改進(jìn)TOPSIS法的城市土地利用績效評價(jià)及障礙因子診斷——以重慶市為例[J]. 資源科學(xué)， 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)濟(jì)帶旅游產(chǎn)業(yè)與生態(tài)環(huán)境交互脅迫關(guān)系驗(yàn)證及協(xié)調(diào)效應(yīng)研究[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)效率時(shí)空演化及其與旅游經(jīng)濟(jì)互動響應(yīng)[J]. 經(jīng)濟(jì)地理， 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è)帶動能力測算及其時(shí)空差異[J]. 地理學(xué)報(bào)， 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)關(guān)系的省際格局及形成機(jī)制分析[J]. 旅游科學(xué)， 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)時(shí)空演變及交互脅迫關(guān)系驗(yàn)證[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)濟(jì)帶生態(tài)系統(tǒng)服務(wù)與經(jīng)濟(jì)社會發(fā)展耦合協(xié)調(diào)關(guān)聯(lián)時(shí)空特征研究[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)服務(wù)價(jià)值的時(shí)空演化及其影響機(jī)制[J]. 地理學(xué)報(bào)， 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)服務(wù)價(jià)值時(shí)空演變研究[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.]