周 斌，趙 靜，青 泉，楊淑群，劉 佳

(1.四川省氣候中心，成都 610072；2.四川省成都市氣象局，成都 611130；3.四川省氣象臺(tái)，成都 610072)

引言

氣候是自然環(huán)境的重要組成部分，也是影響著人們居住適宜性的重要環(huán)境因子[1]。氣候條件對(duì)旅游影響甚廣，任何一個(gè)旅游景區(qū)均難免不受影響[2]。西嶺雪山是國(guó)家AAAA級(jí)風(fēng)景名勝區(qū)，位于四川省成都市大邑縣境內(nèi)，距成都僅95km，總面積483km2，海拔1260～5364m，景區(qū)屬于立體氣溫帶, 雨量充沛、日照少、云霧多、氣候景觀獨(dú)特。年平均氣溫為16℃，1月平均氣溫為4℃，7月平均氣溫為23℃，景區(qū)氣溫較平原地區(qū)低2℃；年降水量>1500mm，年均日照時(shí)數(shù)680h[3]。該區(qū)域夏季最高氣溫25℃，冬季最低氣溫為-6℃，每年11月底到次年3月底為積雪期，積雪厚度為60cm以上，適宜高山滑雪旅游[4]。賞雪滑雪旅游資源屬于氣候旅游資源，與氣象條件有著密不可分的聯(lián)系，例如氣溫以影響降水的形式來(lái)影響、決定冰雪的融化速率[5]。因此降雪積雪時(shí)段的氣象條件預(yù)報(bào)，天氣變化導(dǎo)致的氣象災(zāi)害直接影響外出旅游活動(dòng)[6]。長(zhǎng)期以來(lái)，國(guó)內(nèi)學(xué)者對(duì)冰雪產(chǎn)業(yè)的研究抱有較大興趣，但在滑雪場(chǎng)降雪預(yù)測(cè)這一研究領(lǐng)域文獻(xiàn)較少。目前，降雪預(yù)報(bào)主要通過(guò)統(tǒng)計(jì)或診斷分析[7-10]探尋相關(guān)因子建立降雪預(yù)報(bào)模型，利用環(huán)流相似建立預(yù)報(bào)方法[11]以及對(duì)數(shù)值模式產(chǎn)品[12-13]的訂正等手段來(lái)建立降雪預(yù)報(bào)方法。以往研究偏重于東北[14-15]、華北[16]和西北[17-18]區(qū)域等大范圍的降雪量級(jí)、起止時(shí)間和落區(qū)的預(yù)報(bào)，相關(guān)研究非常多，其預(yù)報(bào)準(zhǔn)確率也較高。而關(guān)于四川區(qū)域，特別是旅游景區(qū)和滑雪場(chǎng)區(qū)域的降雪預(yù)報(bào)，由于資料缺乏的原因，分析研究的非常少。本文在已有的研究基礎(chǔ)上，通過(guò)分析西嶺雪山旅游區(qū)周邊區(qū)域自動(dòng)站溫度和降水與海拔的關(guān)系，結(jié)合城鎮(zhèn)預(yù)報(bào)和EC數(shù)值預(yù)報(bào)產(chǎn)品，研究西嶺雪山旅游區(qū)最低氣溫、降水相態(tài)隨海拔高度的變化特征，進(jìn)而建立西嶺雪山降雪的訂正預(yù)報(bào)方法，以此滿足公眾對(duì)旅游氣象服務(wù)的精細(xì)化需求。

1 資料及方法

1.1 資料來(lái)源

文中使用的資料包括國(guó)家基本站大邑及區(qū)域自動(dòng)站霧山、花水灣、滑雪場(chǎng)、日月坪5站，2005～2017年每年10月～次年2月逐日最低溫度。大邑、霧山、花水灣、滑雪場(chǎng)、日月坪的海拔高度逐漸升高(表1)，特別是滑雪場(chǎng)和日月坪海拔高度在2000m以上，冬季的積雪主要在滑雪場(chǎng)及其以上的區(qū)域，以滑雪場(chǎng)作為重要的參考點(diǎn)，進(jìn)行溫度和降水的訂正預(yù)報(bào)。

表1 西嶺雪山各站海拔高度

1.2 方法

(1)候最低溫度差值訂正：以預(yù)報(bào)的最低溫度為基礎(chǔ)，加上逐候最低溫度差值后得到日最低溫度預(yù)報(bào)；

(2)模式溫度訂正：用距滑雪場(chǎng)最近的格點(diǎn)溫度預(yù)報(bào)數(shù)據(jù)代替滑雪場(chǎng)站點(diǎn)溫度預(yù)報(bào)，選出最低溫度，以氣溫預(yù)報(bào)變化量，建立最低氣溫預(yù)報(bào)的模式溫度訂正方法；

(3)集成預(yù)報(bào)：以候最低溫度差值訂正和模式溫度訂正進(jìn)行加權(quán)平均的最低溫度集成訂正方法。

2 滑雪場(chǎng)最低溫度的訂正

2.1 候最低溫度差值訂正

整理花水灣、霧山、滑雪場(chǎng)、日月坪和大邑各站點(diǎn)2005～2017年每年10月～次年2月每天的最低溫度，求出各站點(diǎn)每候的候平均最低溫度。計(jì)算花水灣、霧山、滑雪場(chǎng)、日月坪和大邑最低溫度差值如表2。

表2 各站和大邑逐候最低溫度差值(℃)

以大邑本站預(yù)報(bào)的24、48、72小時(shí)的最低溫度為基礎(chǔ)，加上逐候最低溫度差值后得到霧山、花水灣、滑雪場(chǎng)、日月坪的日最低溫度預(yù)報(bào)。

Yitmin24=Y0tmin24-Ai

Yitmin48=Y0tmin48-Ai

Yitmin72=Y0tmin72-Ai

其中Y0tmin24、Y0tmin48、Y0tmin72分別是大邑24～72小時(shí)最低溫度預(yù)報(bào)值，i=1、2、3、4，分別表示霧山、花水灣、滑雪場(chǎng)、日月坪，Yitmin24、Yitmin48、Yitmin72分別是某站24～72小時(shí)最低溫度預(yù)報(bào)值，Ai為對(duì)應(yīng)時(shí)段的候最低溫度差值(表2)。

以2017年10月～2018年2月滑雪場(chǎng)24小時(shí)最低溫度預(yù)報(bào)為例，對(duì)其進(jìn)行了試報(bào)，剔除缺資料時(shí)段，試報(bào)結(jié)果如圖1所示，按照預(yù)報(bào)誤差≤2℃為正確，準(zhǔn)確率為：74/134=55.2%。

2.2 模式溫度訂正

因滑雪場(chǎng)海拔高度2152m(約800hPa高度)，以EC025模式800hPa溫度預(yù)報(bào)(103.25°E、30.75°N)格點(diǎn)(距滑雪場(chǎng)最近的格點(diǎn))數(shù)據(jù)代替滑雪場(chǎng)站點(diǎn)溫度預(yù)報(bào)，選出最低溫度，以氣溫預(yù)報(bào)變化量建立最低氣溫預(yù)報(bào)方程：

Y2=T0+T1-T2

其中Y2是滑雪場(chǎng)24～72小時(shí)最低溫度預(yù)報(bào)值；T0為當(dāng)天14時(shí)前的最低溫度。

當(dāng)T1為EC當(dāng)天08時(shí)起報(bào)的24小時(shí)預(yù)報(bào)溫度；T2為EC前1天08時(shí)起報(bào)的48小時(shí)預(yù)報(bào)溫度時(shí)，Y2是滑雪場(chǎng)24小時(shí)最低溫度預(yù)報(bào)值。

當(dāng)T1為EC當(dāng)天08時(shí)起報(bào)的48小時(shí)預(yù)報(bào)溫度；T2為EC前1天08時(shí)起報(bào)的72小時(shí)預(yù)報(bào)溫度時(shí)，Y2是滑雪場(chǎng)48小時(shí)最低溫度預(yù)報(bào)值。

當(dāng)T1為EC當(dāng)天08時(shí)起報(bào)的72小時(shí)預(yù)報(bào)溫度；T2為EC前1天08時(shí)起報(bào)的96小時(shí)預(yù)報(bào)溫度時(shí)，Y2是滑雪場(chǎng)72小時(shí)最低溫度預(yù)報(bào)值。

以2017年10月～2018年2月滑雪場(chǎng)24小時(shí)最低溫度預(yù)報(bào)為例，對(duì)其進(jìn)行了試報(bào)，剔除缺資料時(shí)段，試報(bào)結(jié)果如下圖2所示，按照預(yù)報(bào)誤差≤2℃為正確，準(zhǔn)確率為：86/134=64.2%。

2.3 集成預(yù)報(bào)

對(duì)候最低溫度差值訂正和模式溫度訂正方法進(jìn)行集成(兩種方法不分權(quán)重)：

Y3=(Yitmin24+Y2)/2

Y3=(Yitmin48+Y2)/2

Y3=(Yitmin72+Y2)/2

對(duì)2017年10月～2018年2月滑雪場(chǎng)24小時(shí)最低溫度預(yù)報(bào)進(jìn)行了試報(bào)，剔除缺資料時(shí)段，試報(bào)結(jié)果如下圖3所示，按照預(yù)報(bào)誤差≤2℃為正確，準(zhǔn)確率為：89/134=66.4%。

由于預(yù)報(bào)方法涉及大邑城鎮(zhèn)預(yù)報(bào)，預(yù)報(bào)方法準(zhǔn)確率和大邑城鎮(zhèn)預(yù)報(bào)準(zhǔn)確率有很大關(guān)系。如果大邑城鎮(zhèn)預(yù)報(bào)存在較大誤差，會(huì)引起滑雪場(chǎng)的預(yù)報(bào)誤差增大。剔除大邑城鎮(zhèn)最低溫度預(yù)報(bào)誤差大于2℃的天數(shù)，試報(bào)結(jié)果如下圖4所示，按照預(yù)報(bào)誤差≦2℃為正確，試報(bào)準(zhǔn)確率為：63/77=81.8%。

3 滑雪場(chǎng)降雪預(yù)報(bào)

3.1 閾值的確定

根據(jù)觀測(cè)數(shù)據(jù)，結(jié)合實(shí)地調(diào)查，當(dāng)滑雪場(chǎng)氣溫在5℃以下時(shí)，就可能降雪。因此選取2015～2017年11月～次年2月降水日，滑雪場(chǎng)小時(shí)降水≥0.1mm、氣溫在5℃以下的個(gè)例，分析降水時(shí)對(duì)應(yīng)的平均溫度。統(tǒng)計(jì)結(jié)果如表3所示：滑雪場(chǎng)從10月～次年1月降水時(shí)溫度為下降趨勢(shì)，10月有降水時(shí)的平均溫度為2.6℃，11月為0.93℃，12月為0.16℃，1月為-0.39℃，2月為-0.18℃。據(jù)此，可以把滑雪場(chǎng)出現(xiàn)降雪的溫度閾值確定為2.6℃左右。

表3 滑雪場(chǎng)小時(shí)降水≥0.1mm時(shí)對(duì)應(yīng)的平均氣溫

3.2 滑雪場(chǎng)降雪預(yù)報(bào)

以最低溫度集成訂正方法預(yù)報(bào)西嶺雪山滑雪場(chǎng)的最低溫度，結(jié)合智能網(wǎng)格降水預(yù)報(bào)結(jié)果和降雪最低溫度閾值判別，預(yù)報(bào)滑雪場(chǎng)是否有降雪。

如果智能網(wǎng)格預(yù)報(bào)有降水時(shí)，且0℃≤Tmin24、Tmin48、Tmin72≤2.6℃，預(yù)報(bào)有雨夾雪；Tmin24、Tmin48、Tmin72≤0℃，預(yù)報(bào)有雪(Tmin24、Tmin48、Tmin72分別為24、48、72小時(shí)最低溫度預(yù)報(bào)值)。

當(dāng)智能網(wǎng)格預(yù)報(bào)的降水量R為0～2.4mm，預(yù)報(bào)有小雪；R為2.5～4.9mm，預(yù)報(bào)有中雪；R為5～9.9mm，預(yù)報(bào)有大雪；R≥10mm，預(yù)報(bào)有暴雪。

3.3 滑雪場(chǎng)降雪預(yù)報(bào)檢驗(yàn)

進(jìn)一步統(tǒng)計(jì)2017年2月中滑雪場(chǎng)的11個(gè)降雪個(gè)例(表4)，發(fā)現(xiàn)EC降水量預(yù)報(bào)值均較大，且滑雪場(chǎng)最低溫度均在0℃以下。利用EC格點(diǎn)降水預(yù)報(bào)，結(jié)合滑雪場(chǎng)最低溫度預(yù)報(bào)能預(yù)報(bào)出滑雪場(chǎng)的降雪(表5)。四川省氣象臺(tái)利用四川160個(gè)測(cè)站2019年1～3月24小時(shí)降水及對(duì)應(yīng)的EC24h有效時(shí)效降水預(yù)報(bào)資料，對(duì)EC晴雨預(yù)報(bào)閾值研究得到，在準(zhǔn)確率最高(85.22%)時(shí)對(duì)應(yīng)的閾值為1.8mm。如果以1.8mm 為晴雨閾值，對(duì)2017年2月中滑雪場(chǎng)的11個(gè)降雪個(gè)例降雪預(yù)報(bào)，發(fā)現(xiàn)除了2017年2月2日個(gè)例未預(yù)報(bào)出以外，其它10次降雪均預(yù)報(bào)出來(lái)，預(yù)報(bào)效果較好。

表4 2017年2月滑雪場(chǎng)降雪預(yù)報(bào)及檢驗(yàn)

表5 2017年2月EC降水量預(yù)報(bào)

4 結(jié)論

本文通過(guò)分析西嶺雪山旅游區(qū)周邊區(qū)域自動(dòng)站溫度和降水與海拔的關(guān)系，結(jié)合城鎮(zhèn)預(yù)報(bào)和EC數(shù)值預(yù)報(bào)產(chǎn)品，探討了西嶺雪山降雪的訂正預(yù)報(bào)方法：

(1) 用統(tǒng)計(jì)方法建立以大邑本站最低溫度預(yù)報(bào)值為基礎(chǔ)，加上逐候最低溫度差值后得到滑雪場(chǎng)的日最低溫度預(yù)報(bào)方程：Yitmin24=Y0tmin24-Ai，試報(bào)準(zhǔn)確率為55.2%。

(2) 結(jié)合滑雪場(chǎng)海拔特點(diǎn)，以氣溫預(yù)報(bào)變化量，建立滑雪場(chǎng)最低氣溫預(yù)報(bào)的模式溫度訂正方程：Y2=T0+T1-T2，試報(bào)準(zhǔn)確率為64.2%。

(3) 進(jìn)一步利用候最低溫度差值訂正和模式溫度訂正進(jìn)行加權(quán)平均，建立最低溫度集成訂正方程：Y3=(Yitmin24+Y2)/2，試報(bào)準(zhǔn)確率為66.4%。在此基礎(chǔ)上，剔除大邑城鎮(zhèn)最低溫度預(yù)報(bào)誤差>2℃的天數(shù)，按照預(yù)報(bào)誤差≤2℃為正確，試報(bào)準(zhǔn)確率為：63/77=81.8%。

(4) 以最低溫度集成訂正方法預(yù)報(bào)西嶺雪山滑雪場(chǎng)的最低溫度，結(jié)合智能網(wǎng)格降水預(yù)報(bào)結(jié)果和降雪最低溫度閾值判別，若智能網(wǎng)格預(yù)報(bào)有降水時(shí)，且0℃≤Tmin24、Tmin48、Tmin72≤2.6℃，預(yù)報(bào)有雨夾雪；Tmin24、Tmin48、Tmin72≤0℃，預(yù)報(bào)有雪(Tmin24、Tmin48、Tmin72分別為24、48、72小時(shí)最低溫度預(yù)報(bào)值)。 其中，當(dāng)智能網(wǎng)格預(yù)報(bào)的降水量R為0～2.4mm，預(yù)報(bào)有小雪；R為2.5～4.9mm，預(yù)報(bào)有中雪；R為5～9.9mm，預(yù)報(bào)有大雪；R≥10mm，預(yù)報(bào)有暴雪。

(5) 通過(guò)對(duì)以上三種最低溫度預(yù)報(bào)方法進(jìn)行檢驗(yàn)，發(fā)現(xiàn)最低溫度集成訂正方法預(yù)報(bào)準(zhǔn)確率最高。因此以該方法，結(jié)合智能網(wǎng)格降水預(yù)報(bào)結(jié)果和降雪最低溫度閾值判別，建立了西嶺雪山滑雪場(chǎng)降雪的訂正預(yù)報(bào)方法，該研究為西嶺雪山滑雪場(chǎng)的旅游景觀預(yù)報(bào)方法提供技術(shù)支撐。