羅 潤(rùn),周年興,2
(1.南京師范大學(xué)地理科學(xué)學(xué)院,江蘇 南京 210023) (2.江蘇省地理信息資源開(kāi)發(fā)與利用協(xié)同創(chuàng)新中心,江蘇 南京 210023)
旅游業(yè)的發(fā)展高度依賴于自然資源與生態(tài)環(huán)境[1],環(huán)境污染損害旅游景觀美感,影響旅游體驗(yàn),不利于旅游業(yè)的可持續(xù)發(fā)展. 與噪音污染、水污染、固體廢物污染相較而言,空氣污染通常難以被察覺(jué)和被人體感知,伴隨空氣污染媒體披露的網(wǎng)絡(luò)化,空氣質(zhì)量成為影響游客目的地選擇的潛在外在因素. 因此,厘清空氣質(zhì)量與景區(qū)客流量的關(guān)聯(lián)及空氣質(zhì)量對(duì)景區(qū)客流量的影響,對(duì)于精準(zhǔn)科學(xué)地防治污染、保障旅游者的身心健康、促進(jìn)旅游地可持續(xù)發(fā)展具有重要意義.
國(guó)內(nèi)外學(xué)者針對(duì)空氣污染和旅游的研究主要涵蓋以下方面:(1)旅游對(duì)空氣污染的影響[2]. 研究結(jié)論基本驗(yàn)證了旅游-環(huán)境庫(kù)茲涅茨曲線,旅游對(duì)PM2.5濃度產(chǎn)生“倒U型”影響[3]. (2)空氣污染對(duì)旅游經(jīng)濟(jì)的影響. 研究結(jié)論主要分為兩類:①空氣污染給入境旅游和國(guó)內(nèi)旅游造成經(jīng)濟(jì)損失[4-12];②PM10對(duì)旅游業(yè)整體產(chǎn)生的負(fù)向影響不顯著,但 PM10對(duì)媒體關(guān)注度較高的城市旅游發(fā)展產(chǎn)生顯著抑制效應(yīng)[13]. (3)空氣污染對(duì)旅游者決策、風(fēng)險(xiǎn)感知、情感體驗(yàn)的影響[14-17]. 主要結(jié)論包括兩方面:①空氣污染促使旅游者產(chǎn)生“逃離污染”的旅游行為;②空氣污染的程度對(duì)旅游者情感體驗(yàn)的影響具有差異性,只有當(dāng)空氣質(zhì)量達(dá)到或超過(guò)中度污染時(shí)才會(huì)對(duì)旅游者情感體驗(yàn)產(chǎn)生顯著的負(fù)向影響[18]. (4)空氣污染物與旅游流的時(shí)空關(guān)聯(lián)研究等. 研究結(jié)論主要為霧霾污染與入境旅游呈空間負(fù)相關(guān)[19-20]. 受數(shù)據(jù)限制,現(xiàn)有研究側(cè)重于分析空氣污染對(duì)入境旅游的年際影響,缺乏旅游者空氣質(zhì)量敏感度、不同等級(jí)空氣質(zhì)量與景區(qū)客流量的分析. 本文基于2018—2019年南京市空氣質(zhì)量指數(shù)(AQI)、4A和5A級(jí)景區(qū)客流量數(shù)據(jù),運(yùn)用HP濾波、小波分析法研究空氣質(zhì)量與景區(qū)客流量的時(shí)間關(guān)聯(lián)、空氣質(zhì)量對(duì)景區(qū)客流量的影響. 本研究揭示了空氣質(zhì)量、旅游者空氣質(zhì)量敏感度對(duì)旅游流的影響機(jī)理,豐富了旅游流影響因素研究,為相關(guān)部門制定環(huán)境和旅游政策提供決策參考,促進(jìn)空氣污染協(xié)同治理和旅游業(yè)高質(zhì)量發(fā)展.
1.1.1 HP濾波法
Hodrick-Prescott Filter(HP濾波)能將AQI、景區(qū)客流量時(shí)間序列中的周期性與趨勢(shì)性相剝離[21]. 與直接繪制折線圖并觀察發(fā)展趨勢(shì)相較而言,運(yùn)用HP濾波法分析發(fā)展趨勢(shì)更具科學(xué)性和客觀性.
1.1.2 小波分析法
小波分析可用于分析AQI、景區(qū)客流量的冷熱分布、周期長(zhǎng)度、周期峰谷等. 運(yùn)用MATLAB軟件計(jì)算AQI、景區(qū)客流量的小波實(shí)部和小波方差等. 基于小波方差值分析AQI和景區(qū)客流量的第一主周期對(duì)應(yīng)的時(shí)間尺度,基于小波實(shí)部圖分析冷熱分布等[22-25]. 顏色越深(淺)表示AQI或景區(qū)客流量越小(大).
1.1.3 多元回歸分析法
為分析空氣質(zhì)量對(duì)旅游景區(qū)客流量的影響,建立如下線性和非線性回歸方程,以調(diào)整R2最大為方程最優(yōu)解的判定原則:
lnTour=α1ln AQI+β1lnTemp+β2lnWind+β3Wea+β4Hol+β5Seas+β6Year+c,
(1)
lnTour=α1(ln AQI)2+α2ln AQI+β1lnTemp+β2lnWind+β3Wea+β4Hol+β5Seas+β6Year+c,
(2)
(3)
在相關(guān)研究的基礎(chǔ)上[26-28],選取AQI表征空氣質(zhì)量狀況. 選取日均溫、風(fēng)力、天氣、節(jié)假日、季節(jié)等作為控制變量,景區(qū)游客量作為因變量(表1). 為消除異方差影響,對(duì)變量取對(duì)數(shù). 將風(fēng)力級(jí)別作為風(fēng)力的評(píng)價(jià)指標(biāo),對(duì)日最高氣溫和日最低氣溫取平均數(shù)得到日均溫,刪除日均溫為非正數(shù)的天數(shù). 此外,構(gòu)建 4個(gè)虛擬變量:①節(jié)假日. 節(jié)假日是旅游者出游的決定性因素,研究將法定節(jié)假日和非工作日設(shè)置為1,將工作日設(shè)置為0. ②天氣狀況.良好的天氣狀況通常促進(jìn)旅游者出游,研究將雨雪天氣設(shè)置為1,多云和晴天設(shè)置為0. ③季節(jié). 旅游流具有明顯的季節(jié)性特征[29-31],春、秋季節(jié)旅游景觀審美價(jià)值和旅游者體感舒適度較高,游客量較多,夏、冬季節(jié)旅游者體感舒適度較低,游客量較少. 季節(jié)虛擬變量將春(3、4、5月)、秋(9、10、11月)設(shè)置為1,夏(6、7、8月)、冬(1、2、12月)設(shè)置為0. ④年份. 為控制年際變化對(duì)景區(qū)客流量造成的影響,將2018年、2019年分別設(shè)置為0、1.
表1 變量的描述性統(tǒng)計(jì)Table 1 Descriptive statistics of variables
南京市2018—2019年AQI日度數(shù)據(jù)來(lái)源于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心,4A和5A級(jí)旅游景區(qū)日接待客流量數(shù)據(jù)來(lái)源于江蘇省文化和旅游廳景區(qū)系統(tǒng). 2018—2019年日均溫、天氣狀況、風(fēng)力風(fēng)向等數(shù)據(jù)來(lái)源于天氣網(wǎng)站,節(jié)假日數(shù)據(jù)根據(jù)中國(guó)政府網(wǎng)的節(jié)假日安排整理獲得.
圖1 HP濾波和小波分析Fig.1 HP filtering and wavelet analysis
研究基于2018年1月—2019年12月南京市4A和5A級(jí)景區(qū)客流量和AQI數(shù)據(jù),運(yùn)用HP濾波、小波分析法從趨勢(shì)性、周期性和季節(jié)性視角分析AQI與旅游景區(qū)客流量的關(guān)聯(lián)特征.由圖1知,2018年1月—2019年12月南京市AQI呈下降趨勢(shì),景區(qū)客流量呈上升趨勢(shì),可見(jiàn)AQI與景區(qū)客流量在時(shí)間趨勢(shì)方面具有負(fù)向關(guān)聯(lián)特征.由小波實(shí)部等值圖知,在10月的時(shí)間尺度附近,AQI和景區(qū)客流量呈明顯的周期性.在該時(shí)間尺度,AQI與景區(qū)客流量呈現(xiàn)“春秋冬負(fù)向關(guān)聯(lián)、夏季正向關(guān)聯(lián)”的特征.其中,春秋季南京市AQI較低,而景區(qū)客流量規(guī)模較大;冬季AQI較高,而景區(qū)客流量規(guī)模較小;夏季AQI較低,且景區(qū)客流量規(guī)模較小.需要注意的是,這里的正向關(guān)聯(lián)指的是AQI與旅游景區(qū)客流量?jī)蓚€(gè)序列在發(fā)展趨勢(shì)、季節(jié)性和周期性方面的發(fā)展規(guī)律具有一致性或相似性,負(fù)向關(guān)聯(lián)指的是兩者的發(fā)展趨勢(shì)、季節(jié)性和周期性規(guī)律相反.
上述分析表明AQI和旅游景區(qū)客流量具有一定的關(guān)聯(lián)性,采用回歸模型進(jìn)一步判別AQI對(duì)旅游流的影響作用.經(jīng)檢驗(yàn),變量間不存在嚴(yán)重的共線性問(wèn)題且時(shí)間序列是平穩(wěn)的.AQI和景區(qū)客流量的協(xié)整檢驗(yàn)結(jié)果顯示兩者存在長(zhǎng)期協(xié)整關(guān)系.格蘭杰因果檢驗(yàn)結(jié)果顯示,AQI是景區(qū)客流量的格蘭杰原因,景區(qū)客流量不是AQI的格蘭杰原因.回歸結(jié)果(表2)顯示,線性模型的擬合優(yōu)度最高.遵循調(diào)整R2最大的原則,選取線性模型為最優(yōu)解釋模型.AQI對(duì)景區(qū)客流量的影響系數(shù)為-0.011 5,但不顯著,說(shuō)明AQI對(duì)景區(qū)客流量產(chǎn)生微弱的負(fù)向影響.在控制變量當(dāng)中,節(jié)假日、季節(jié)、日均溫對(duì)旅游景區(qū)客流量產(chǎn)生顯著正向影響,天氣對(duì)景區(qū)客流量產(chǎn)生了顯著負(fù)向影響,這與預(yù)期相符.風(fēng)力因素對(duì)景區(qū)客流量具有正向影響,但不顯著.
表2 回歸結(jié)果Table 2 Regression results
依據(jù)國(guó)家對(duì)AQI的分級(jí)標(biāo)準(zhǔn),將南京市空氣質(zhì)量劃分為6個(gè)等級(jí)(表3).
表3 AQI等級(jí)劃分標(biāo)準(zhǔn)Table 3 Criterion of classification of AQI
為進(jìn)一步分析不同等級(jí)空氣質(zhì)量對(duì)景區(qū)客流量的影響,建立如下方程:
lnTour=α1Exc+β1lnTemp+β2lnWind+β3Wea+β4Hol+β5Seas+β6Year+c,
(4)
lnTour=α1Poll+β1lnTemp+β2lnWind+β3Wea+β4Hol+β5Seas+β6Year+c,
(5)
式(4)和式(5)中,Exc(空氣優(yōu)或良為1,其他為0)和Poll(輕度、中度和重度污染為1,其他為0)分別為空氣優(yōu)良、空氣污染的虛擬變量. 如表4所示:優(yōu)良的空氣質(zhì)量對(duì)景區(qū)客流量的影響系數(shù)為0.005 2,但不顯著,表明優(yōu)良的空氣質(zhì)量對(duì)景區(qū)客流量具有微弱的促進(jìn)作用. 輕度、中度和重度污染對(duì)景區(qū)客流量的影響系數(shù)為-0.005 2,但不顯著,表明空氣污染對(duì)景區(qū)客流量具有微弱的抑制作用.
表4 不同等級(jí)的空氣質(zhì)量對(duì)景區(qū)客流量的影響Table 4 The impact of different grades of air quality on tourist flow of scenic spots
為進(jìn)一步分析AQI對(duì)不同類型景區(qū)客流量的影響,按照室內(nèi)外標(biāo)準(zhǔn)將南京市20家景區(qū)劃分為9家室內(nèi)和11家戶外景區(qū). 如表5所示:AQI對(duì)室內(nèi)景區(qū)、戶外景區(qū)的影響系數(shù)分別為-0.028 3、0.019 0,表明空氣質(zhì)量指數(shù)對(duì)室內(nèi)景區(qū)客流量產(chǎn)生的負(fù)向影響大于對(duì)戶外景區(qū)客流量產(chǎn)生的影響. 室內(nèi)景區(qū)一般為博物館等人文類景區(qū),戶外景區(qū)多為森林、山岳、湖泊等自然類景區(qū). 節(jié)假日、季節(jié)、溫度、天氣狀況等因素對(duì)戶外景區(qū)客流量的影響顯著大于對(duì)室內(nèi)景區(qū)客流量的影響. 與室內(nèi)景區(qū)相較而言,戶外景區(qū)客流量更多受到閑暇時(shí)間和氣候條件的影響,旅游者對(duì)空氣質(zhì)量的敏感度較低,因此戶外景區(qū)客流量受到空氣質(zhì)量的影響較小. 由于南京市空氣質(zhì)量以優(yōu)良為主,AQI的小幅度上升不會(huì)對(duì)戶外景區(qū)客流量產(chǎn)生顯著的抑制作用.
表5 空氣質(zhì)量對(duì)不同類型景區(qū)客流量的影響Table 5 The impact of air quality on tourist flow of different types of scenic spots
為進(jìn)一步探究旅游者空氣質(zhì)量敏感度對(duì)AQI與景區(qū)客流量關(guān)系的調(diào)節(jié)效應(yīng),引入旅游者空氣質(zhì)量敏感度指標(biāo). “南京空氣質(zhì)量”百度搜索指數(shù)是潛在旅游者通過(guò)百度主動(dòng)關(guān)注與搜索南京空氣質(zhì)量的數(shù)值,反映了潛在旅游者對(duì)南京空氣質(zhì)量的敏感度和關(guān)注度,因此研究用其表征旅游者空氣質(zhì)量敏感度并建立如下方程:
lnTour=α1ln AQI+α2Sens+α3ln AQI_Sens+β1lnTemp+β2lnWind+
β3Wea+β4Hol+β5Seas+β6Year+c,
(6)
式中,Sens為“南京空氣質(zhì)量”百度指數(shù)的虛擬變量,分為高(百度指數(shù)高于均值,Sens=1)和低(百度指數(shù)等于或低于均值,Sens=0)兩類.ln AQI_Sens為AQI與旅游者空氣質(zhì)量敏感度的交互項(xiàng). 回歸結(jié)果如表6所示:AQI對(duì)景區(qū)客流量的影響系數(shù)為-0.004 8,但不顯著,Sens對(duì)景區(qū)客流量的影響系數(shù)為0.105 4,但不顯著,說(shuō)明旅游者的空氣質(zhì)量意識(shí)對(duì)景區(qū)客流量產(chǎn)生微弱的促進(jìn)作用,這是由于南京市空氣質(zhì)量位居全國(guó)前列,旅游者對(duì)南京市空氣質(zhì)量的關(guān)注度越高,優(yōu)良的空氣質(zhì)量越促進(jìn)旅游者到南京旅游.ln AQI_Sens對(duì)景區(qū)客流量的影響系數(shù)為-0.019 7,但不顯著,表明當(dāng)旅游者對(duì)南京市空氣質(zhì)量敏感度較高時(shí),AQI對(duì)景區(qū)客流量的影響系數(shù)為-0.024 5(-0.004 8+(-0.019 7)),AQI每提高1個(gè)單位,景區(qū)客流量減少2.45%. 因此,旅游者對(duì)空氣質(zhì)量的高敏感度會(huì)增強(qiáng)AQI對(duì)景區(qū)客流量的負(fù)向影響.
表6 旅游者空氣質(zhì)量敏感度對(duì)空氣質(zhì)量與景區(qū)客流量關(guān)系的調(diào)節(jié)作用Table 6 The moderating effect of tourists’ air quality sensitivityon the relationship between AQI and tourist flow
借鑒王兆華等[32]的內(nèi)生性處理方法,選用距離南京市最近的上風(fēng)向城市的AQI作為工具變量,設(shè)為lnUpwind_AQI. 由于風(fēng)向是自然因素,符合工具變量外生性要求,同時(shí)上風(fēng)向城市的AQI與南京市的AQI 的相關(guān)性為0.883 6,符合工具變量與核心解釋變量相關(guān)性的要求. 運(yùn)用2SLS的回歸結(jié)果(表7)顯示,AQI對(duì)景區(qū)客流量的影響系數(shù)為-0.019 3,但不顯著,回歸結(jié)果基本不變. 借鑒沈永建等[33]的穩(wěn)健性檢驗(yàn)方法,將核心解釋變量滯后一期,滯后一期的AQI對(duì)景區(qū)客流量的影響系數(shù)為-0.031 0,但不顯著,回歸結(jié)果基本不變,研究結(jié)論較為穩(wěn)健.
表7 穩(wěn)健性與內(nèi)生性檢驗(yàn)結(jié)果Table 7 Results of robustness and endogeneity test
(1)南京市AQI與旅游景區(qū)客流量的總體發(fā)展趨勢(shì)存在負(fù)向關(guān)聯(lián),AQI呈下降趨勢(shì),旅游景區(qū)客流量呈上升趨勢(shì).AQI和景區(qū)客流量在不同季節(jié)的關(guān)聯(lián)性存在差異. (2)南京市AQI對(duì)景區(qū)客流量產(chǎn)生微弱的負(fù)向影響,不同等級(jí)的空氣質(zhì)量對(duì)景區(qū)客流量產(chǎn)生的影響具有差異,優(yōu)良的空氣質(zhì)量對(duì)景區(qū)客流量產(chǎn)生微弱的促進(jìn)作用,輕度污染、中度污染和重度污染對(duì)旅游景區(qū)客流量產(chǎn)生微弱的抑制作用. (3)空氣質(zhì)量對(duì)不同類型旅游景區(qū)客流量的影響存在差異,AQI對(duì)室內(nèi)景區(qū)客流量的負(fù)向影響大于對(duì)戶外景區(qū)客流量的影響. (4)旅游者對(duì)南京市空氣質(zhì)量的高敏感度和關(guān)注度會(huì)增強(qiáng)AQI對(duì)旅游景區(qū)客流量的抑制作用.