重慶地區(qū)極端高溫事件的模擬研究

魏麟驍，李永華，何卷雄，唐紅玉，辛?xí)愿?/p>

(1.重慶市氣候中心，重慶　401147；2.中國(guó)科學(xué)院大氣物理研究所國(guó)際氣候與環(huán)境科學(xué)中心，北京　100029；3.中國(guó)氣象局國(guó)家氣候中心，北京　100081)

1　引言

重慶位于西南地區(qū)東部，屬亞熱帶季風(fēng)性濕潤(rùn)氣候，地處長(zhǎng)江上游，氣候?yàn)?zāi)害頻發(fā)，尤其是與高溫有關(guān)的災(zāi)害較多。近幾十年來(lái)，由于人類活動(dòng)和自然因素的綜合影響，全球極端降水和高溫事件頻發(fā)，而重慶地區(qū)的高溫災(zāi)害尤其嚴(yán)重，對(duì)經(jīng)濟(jì)和社會(huì)的發(fā)展造成了很大影響。20世紀(jì)中期以來(lái)，中國(guó)地區(qū)的暖事件增加而冷事件減少，并且暖事件對(duì)平均態(tài)改變的敏感程度比冷事件高[1]。例如2003年夏季，高溫?zé)崂讼須W洲，打破了近500 a的記錄，造成上萬(wàn)人死亡[2]，同年，我國(guó)長(zhǎng)江流域也出現(xiàn)了大范圍的高溫事件。2006年夏季，重慶和四川地區(qū)遭受特大伏旱，造成的災(zāi)害十分嚴(yán)重[3]。2013年夏季，我國(guó)南方地區(qū)遭受大范圍的熱浪襲擊，影響的范圍、持續(xù)時(shí)間和強(qiáng)度均為歷史罕見。鑒于此，我國(guó)很多學(xué)者對(duì)極端高溫開展了大量的研究。況雪源等[4]對(duì)近50多年來(lái)全國(guó)群發(fā)性高溫事件進(jìn)行了識(shí)別，發(fā)現(xiàn)2000年以后的群發(fā)性高溫事件都是5個(gè)年代中最顯著的，并且長(zhǎng)江以南地區(qū)發(fā)生頻率最高。唐恬等[5]對(duì)2013年我國(guó)南方高溫天氣進(jìn)行了分析，發(fā)現(xiàn)這次高溫事件的發(fā)生次數(shù)突破了歷史記錄，其中，8月的極端高溫十分突出。因此，提前對(duì)高溫事件做出準(zhǔn)確預(yù)測(cè)對(duì)國(guó)計(jì)民生意義重大，也是氣象部門的服務(wù)重點(diǎn)之一。

隨著我國(guó)短期氣候預(yù)測(cè)系統(tǒng)的不斷發(fā)展，月動(dòng)力延伸期預(yù)測(cè)(DERF2.0)產(chǎn)品在氣候預(yù)測(cè)業(yè)務(wù)中得到廣泛應(yīng)用[6,7]，預(yù)測(cè)水平逐年提升，但是全球模式對(duì)于區(qū)域氣候的模擬能力有限，在業(yè)務(wù)應(yīng)用中需要使用更高分辨率的區(qū)域氣候模式對(duì)地表氣溫和極端高溫事件進(jìn)行模擬和預(yù)測(cè)。高學(xué)杰等[8]使用區(qū)域氣候模式單向嵌套全球模式的輸出結(jié)果對(duì)東亞地區(qū)的氣候進(jìn)行了模擬和預(yù)估，發(fā)現(xiàn)區(qū)域模式提供了氣溫和降水更詳細(xì)的空間分布信息，對(duì)夏季的季風(fēng)降水型模擬也有了很大改進(jìn)。李東歡等[9]利用中國(guó)科學(xué)院大氣物理研究所全球模式 FGOALS-g2 的模擬數(shù)據(jù)驅(qū)動(dòng)區(qū)域氣候模式 RegCM3，發(fā)現(xiàn)相對(duì)于全球模式，區(qū)域模式由于水平分辨率較高，能在刻畫地表氣溫分布的細(xì)節(jié)上體現(xiàn)出優(yōu)勢(shì)。余正[10]通過(guò)對(duì)1961—2005年的模擬研究表明，IAP AGCM4.1模式可以較好的模擬出高溫日數(shù)的空間分布特征，但是高估了華北、江淮流域和川渝地區(qū)等地的高溫日數(shù)。許建偉等[11]利用NCEP/DOE再分析資料驅(qū)動(dòng)WRF模式對(duì)1999—2008年夏季黑河流域及其周邊地區(qū)的氣溫和降水進(jìn)行了模擬，結(jié)果表明，在復(fù)雜地形條件下，WRF模式較再分析資料能夠更精細(xì)地模擬出復(fù)雜地形條件下氣溫和降水的分布特征。由此可見采用動(dòng)力降尺度的方法進(jìn)行短期氣候預(yù)測(cè)是具有一定優(yōu)勢(shì)的，本文的目的是分析重慶短期氣候預(yù)測(cè)系統(tǒng)(IAP AGCM4.1和WRF的耦合模式)對(duì)重慶夏季地表氣溫和極端高溫事件的模擬能力，并對(duì)使用該系統(tǒng)預(yù)測(cè)極端高溫事件提出建議。

2　資料和方法

2.1　模式簡(jiǎn)介

使用的耦合模式(IAP AGCM_WRF)是在中國(guó)科學(xué)院大氣物理研究所IAP AGCM4.1和美國(guó)的中尺度數(shù)值模式WRF3.2的基礎(chǔ)上建立的。IAP-AGCM4.1模式的水平分辨率為1.4°×1.4°，模式層頂為10 hPa。大氣環(huán)流模式的初始條件來(lái)自于NCEP/DOE再分析資料，海溫(SST)和海冰(SEA ICE)來(lái)自于美國(guó)國(guó)家海洋和大氣管理局(NOAA)的逐周OISST資料。中尺度模式的水平分辨率為20 km，中心點(diǎn)位于(29°N，106°E)，經(jīng)緯方向的格點(diǎn)數(shù)為111×111，范圍覆蓋了我國(guó)西南的大部分地區(qū)，模式的垂直方向?yàn)?1層。該系統(tǒng)的運(yùn)行步驟是首先將IAP-AGCM4.1每隔6 h輸出的氣候要素場(chǎng)進(jìn)行空間插值，生成驅(qū)動(dòng)WRF模式的初始和邊界條件，接下來(lái)再運(yùn)行WRF模式，得到預(yù)測(cè)結(jié)果。在正式開始試驗(yàn)前，我們進(jìn)行了大量的參數(shù)化方案敏感性試驗(yàn)(2000—2014年)，通過(guò)與再分析資料的對(duì)比，選出了模擬效果較好且運(yùn)行比較穩(wěn)定的參數(shù)化方案組合(見表1)。同時(shí)，也使用了NCEP/DOE再分析資料為驅(qū)動(dòng)場(chǎng)開展了模擬，參數(shù)化方案的配置、海溫和海冰覆蓋資料與上一個(gè)試驗(yàn)相同，2個(gè)試驗(yàn)的積分時(shí)間段都為1982年1月1日—2014年12月31日。

根據(jù)前人的研究，區(qū)域氣候模式對(duì)于將其驅(qū)動(dòng)的大尺度模式的模擬結(jié)果會(huì)有較大改進(jìn)，那么與觀測(cè)資料驅(qū)動(dòng)的區(qū)域氣候模式的模擬結(jié)果相比其差距如何？這是本文要解決的問(wèn)題。本文將分析IAP-AGCM4.1驅(qū)動(dòng)的WRF3.2模式對(duì)重慶夏季地表氣溫及極端高溫事件的模擬結(jié)果，并討論IAP-AGCM4.1驅(qū)動(dòng)的WRF模式(IWRF)與再分析資料(NCEP/DOE)驅(qū)動(dòng)的WRF模式(NWRF)在模擬結(jié)果上的差異。本文的結(jié)論對(duì)于探討重慶短期氣候預(yù)測(cè)系統(tǒng)在預(yù)測(cè)極端高溫的能力上具有一定的參考價(jià)值，并對(duì)如何使用該系統(tǒng)提出建議。

2.2　觀測(cè)資料和分析方法

使用CN05.1格點(diǎn)資料來(lái)檢驗(yàn)?zāi)Ｊ降哪M性能，該資料的分辨率為0.25°×0.25°，使用了中國(guó)境內(nèi)2 400個(gè)臺(tái)站的觀測(cè)資料，采用距平逼近法，由氣候場(chǎng)和距平場(chǎng)分別插值后疊加得到[12,13]。使用的平均氣溫、最高和最低氣溫資料，起始時(shí)間為1982年1月1日—2014年12月31日。為了方便與觀測(cè)資料進(jìn)行比較，先將2個(gè)模式的結(jié)果統(tǒng)一插值到與觀測(cè)資料相同的格點(diǎn)上。對(duì)模擬結(jié)果進(jìn)行檢驗(yàn)所用的統(tǒng)計(jì)量包括模擬結(jié)果與觀測(cè)場(chǎng)的偏差、均方根誤差、皮爾遜相關(guān)系數(shù)，以上3類檢驗(yàn)方法的計(jì)算公式分別為：

(1)

(2)

(3)

表1　模式選項(xiàng)設(shè)置Tab.1　Main physics options in the experiment

根據(jù)重慶地區(qū)地表氣溫的分布情形以及重慶市氣候中心預(yù)測(cè)業(yè)務(wù)的常用劃分方法，將重慶地區(qū)劃分為4個(gè)子區(qū)域，如圖1所示。

圖1　重慶地形以及文中的4個(gè)分區(qū)。黑色方框分別表示4個(gè)分區(qū)：西部(W:28.5～30.5°N，105～107°E)、中部(M:29～31°N，107～108°E)、東北部(NE:30～32°N，108～110°E)和東南部(SE:28～30°N，108～110°E)Fig.1　The topography of the east of the southwest area (shading,units:m) and distribution of 4 sub-regions in the east of the southwest area.The boxes illustrate the 4 sub-regions:western region(W:28.5～30.5°N，105～107°E),central region(M:29～31°N，107～108°E),northeastern region(NE:30～32°N，108～110°E) and southeastern region(SE:28～30°N，108～110°E)

名稱縮寫定義單位熱浪持續(xù)指數(shù)HWDI至少持續(xù)3 d,日最高溫度不低于基準(zhǔn)期第90個(gè)百分位數(shù)的天數(shù)d暖晝指數(shù)TX90p每年日最高氣溫大于或者等于基準(zhǔn)期內(nèi)90%分位數(shù)的天數(shù)占統(tǒng)計(jì)日數(shù)的百分比d日最高氣溫最大值TXx每年日最高氣溫最大值℃

3　結(jié)果分析

3.1　夏季地表氣溫氣候態(tài)模擬

從觀測(cè)來(lái)看，重慶地區(qū)的平均、最高和最低氣溫氣候態(tài)的空間分布特征較為相似：高溫區(qū)出現(xiàn)在盆地中部和東南部，低溫區(qū)主要出現(xiàn)在高海拔地區(qū)。

IWRF模式可以較好的模擬出3個(gè)地表氣溫指標(biāo)的空間分布，能夠再現(xiàn)盆地中部和東南部地區(qū)的高溫中心，但是模擬的高溫中心量級(jí)偏大(圖略)。IWRF模擬的平均氣溫在盆地中部和東北、東南大部地區(qū)主要表現(xiàn)為1～3 ℃的暖偏差，高海拔地區(qū)為冷偏差；模擬的最高氣溫的偏差分布與平均氣溫相似，盆地南部的最大偏差在3 ℃左右，高海拔地區(qū)的偏差也較平均氣溫大；相對(duì)來(lái)說(shuō)，IWRF對(duì)于最低氣溫模擬的較好，在盆地四周為1 ℃左右的暖偏差。

NWRF模擬的T2m和Tmax的高值中心在量級(jí)上要好于IWRF，模擬的Tmin的分布特征與IWRF相近(圖略)。在盆地內(nèi)部和東南部地區(qū)，NWRF模擬的T2m和Tmax的偏差均小于IWRF，尤其對(duì)于盆地南部改進(jìn)最大，這與使用了較為準(zhǔn)確的強(qiáng)迫場(chǎng)有很大關(guān)系，但是NWRF對(duì)于Tmin的模擬效果沒(méi)有IWRF的效果好。通過(guò)IWRF與NWRF之間的差值場(chǎng)可以看到，IWRF主要在盆地中部出現(xiàn)暖偏差，偏差幅度在1～2 ℃之間。

圖2　1982—2014年夏季重慶地表氣溫的氣候態(tài)和模式模擬偏差(單位：℃)，從上到下依次為日平均氣溫T2m、日最高氣溫Tmax、日最低氣溫Tmin；左列為觀測(cè)資料的氣候態(tài)分布，第2、第3列分別為IWRF和NWRF的模擬偏差,右列為IWRF減去NWRF的模擬結(jié)果Fig.2　Climatology and simulation bias of summer surface air temperature during 1982-2014 (units:℃).The top,middle and bottom row indicates the results of t2m,tmax,and tmin,respectively.The first column indicates the results from CN05.1,the second column indicates the bias between IWRF and CN05.1,the third column indicates the bias between NWRF and CN05.1,and the last column indicates the bias between IWRF and NWRF.

模式時(shí)間空間相關(guān)系數(shù)TaveTmaxTmin差值TaveTmaxTmin均方根誤差TaveTmaxTminIWRFNWRF逐月逐月0.960.950.940.930.970.970.980.200.67-0.37-0.14-0.791.520.841.951.220.910.99

IWRF和NWRF模擬的3個(gè)地表氣溫指標(biāo)與觀測(cè)的空間相關(guān)系數(shù)均在0.9以上，說(shuō)明這2個(gè)試驗(yàn)結(jié)果均對(duì)夏季氣溫的空間分布有較好的表征。從偏差和均方根誤差來(lái)看，NWRF除最低氣溫較IWRF大外，平均氣溫和最高氣溫均優(yōu)于IWRF。總體來(lái)看，這2個(gè)試驗(yàn)對(duì)于重慶夏季氣溫的模擬表現(xiàn)的都較好。

圖3給出了2個(gè)試驗(yàn)對(duì)夏季重慶氣溫概率密度分布的模擬情況，可以看到，2個(gè)試驗(yàn)對(duì)地表氣溫3個(gè)指標(biāo)的概率密度分布都有較好的模擬，IWRF與NWRF的概率分布相似，對(duì)最低氣溫模擬的最好，對(duì)于其余2個(gè)指標(biāo)，IWRF與觀測(cè)相比右移，其模擬的極端高溫略偏多，與圖2f相對(duì)應(yīng)。

圖3　1982—2014年中國(guó)夏季地表氣溫概率密度分布(單位：℃)，(a)T2m、(b)Tmax、(c)Tmin。圖中的概率密度分布使用所選時(shí)段、所選區(qū)域中各個(gè)點(diǎn)的逐日數(shù)據(jù)(不是區(qū)域平均) 計(jì)算而成。黑色表示觀測(cè)資料CN05.1中的結(jié)果，紅色表示IWRF模擬的結(jié)果，藍(lán)色表示NWRF模擬的結(jié)果Fig.3　Frequency distribution of daily (a)t2m,(b)tmax and(c)tmin from individual grid (not average).Black indicates the results from CN05.1,red indicates the results from IWRF and green indicates the results from NWRF

根據(jù)上文所劃分的4個(gè)子區(qū)域，對(duì)地表氣溫進(jìn)行區(qū)域平均可以更直觀地考察2個(gè)試驗(yàn)對(duì)地表氣溫的模擬能力。從圖4(a,b,c)可以看到，除了最低氣溫，IWRF模擬的其余2個(gè)氣溫指標(biāo)普遍高于觀測(cè)。IWRF對(duì)西部和中部地區(qū)的模擬偏差較大。使用泰勒?qǐng)D可以更直觀的表征模擬和觀測(cè)的地表氣溫的空間相似程度，2個(gè)試驗(yàn)的空間分布?xì)w一化標(biāo)準(zhǔn)差大多介于0.5～2.5之間。從空間相關(guān)系數(shù)來(lái)看，對(duì)于各個(gè)子區(qū)域的空間分布都有較好的模擬。總體來(lái)看，IWRF模擬的西部和中部地區(qū)的平均氣溫和最高氣溫偏差較大，NWRF的表現(xiàn)好于IWRF，它們對(duì)最低氣溫的模擬效果較好。

圖4　區(qū)域平均的地表氣溫(單位：℃),(a)、(b)、(c)分別為T2m、Tmax和Tmin，黑色表示CN05.1，紅色表示IWRF，藍(lán)色表示NWRF；(d)T2m、(e)Tmax和(f)Tmin為泰勒?qǐng)D，紅色表示IWRF模擬結(jié)果，綠色表示NWRF模擬結(jié)果,1～4分別表示文中劃分的W、M、NE、SE 4個(gè)子區(qū)域Fig.4　(a),(b) and (c)indicates the regional averaged summer surface air temperature (units:℃),t2m,tmax and tmin,respectively,black indicates CN05.1,red indicates IWRF and blue indicates NWRF.(d),(e) and (f) indicates the Taylor diagram for displaying pattern statistics of T2m,Tmax,Tmin,respectively.Red indicates the results from IWRF,and blue indicates the results from NWRF.1-4 represents W,M,NE and SE.

3.2　夏季極端高溫事件氣候態(tài)模擬

從圖5可知，日最高氣溫最大值(TXx)與地表氣溫指標(biāo)的分布相似，都是在西部地區(qū)出現(xiàn)大值中心，IWRF和NWRF都能較好的模擬出TXx的空間分布和強(qiáng)度，兩者對(duì)于東北部和東南部地區(qū)的模擬存在較強(qiáng)的冷偏差，這可能是由于地形的原因造成的。觀測(cè)場(chǎng)上，熱浪持續(xù)指數(shù)(HWDI)的大值中心出現(xiàn)在四川東北部地區(qū)，并向東擴(kuò)展到重慶的中部和東北部地區(qū)，基本在6 d左右，IWRF和NWRF與實(shí)況相比模擬的HWDI偏少，就偏差圖來(lái)看，IWRF的結(jié)果優(yōu)于NWRF，偏差區(qū)域主要出現(xiàn)在四川的東北部，偏少1～2 d左右，而NWRF模擬的模擬結(jié)果在整個(gè)區(qū)域都偏少。觀測(cè)場(chǎng)上，暖晝指數(shù)(TX90p)的大值中心出現(xiàn)在重慶大部以及湖北和湖南部分地區(qū)，IWRF的模擬結(jié)果仍然好于NWRF，偏差主要出現(xiàn)在東北部部分地區(qū)，偏少0.5～1 d；NWRF的模擬結(jié)果在重慶大部地區(qū)都偏少，尤其在東北部地區(qū)偏少最多，偏少1～2 d。從極端高溫指數(shù)的空間分布來(lái)看，除了TXx，IWRF對(duì)HWDI和TX90p的模擬要好于NWRF。

圖5　1982—2014年夏季重慶極端高溫氣候態(tài)，從上到下依次為日最高氣溫最大值TXx(單位：℃)、熱浪持續(xù)指數(shù)HWDI(單位：d)、暖晝指數(shù)TX90p(單位：d)；左列為觀測(cè)場(chǎng)的氣候態(tài)結(jié)果，第2、第3列分別為IWRF和NWRF的模擬偏差,右列為IWRF減去NWRF的模擬結(jié)果Fig.5　Climatology and simulation bias of summer high temperature indices during 1982-2014.The top row indicates the results of TXx (units:℃),the middle row indicates the results of HWDI (units:day),and the bottom row indicates the results of TX90p (units:day).The left column indicates the results from CN05.1,the second column indicates the bias between IWRF and CN05.1,the third column indicates the bias between NWRF and CN05.1,and the last column indicates the bias between IWRF and NWRF

圖6為2個(gè)試驗(yàn)對(duì)于極端高溫事件區(qū)域平均的模擬結(jié)果，IWRF和NWRF模擬的TXx(圖6a)在東北部地區(qū)偏低較為明顯，偏差幅度在2～3 ℃之間，但是，對(duì)西部地區(qū)模擬的較好。IWRF模擬的HWDI在各區(qū)域都要好于NWRF，在西部和東南部地區(qū)模擬的較好，NWRF在各區(qū)域模擬的HWDI都偏少2～4 d。在這4個(gè)區(qū)域，它們對(duì)TX90p都模擬的很好。本文也計(jì)算了2個(gè)試驗(yàn)?zāi)M的極端高溫事件的空間相關(guān)系數(shù)和均方根誤差，IWRF模擬的TXx 4個(gè)區(qū)域的空間相關(guān)系數(shù)均在0.8以上，但是均方根誤差較大，最大出現(xiàn)在東北地區(qū)，NWRF的模擬結(jié)果與IWRF相似。它們對(duì)于HWDI和TX90p的空間分布模擬的較差，但是在西部地區(qū)和東南部地區(qū)有一定的模擬能力，相關(guān)系數(shù)在0.3左右，均方根誤差也較小?？傊?，這2個(gè)試驗(yàn)對(duì)于TX90p和HWDI的空間分布模擬能力弱于對(duì)TXx的模擬。

圖6　區(qū)域平均的極端高溫指數(shù)，(a)、(b)、(c)分別為TXx(單位：℃)、HWDI(單位：d)和TX90p(單位：d)，黑色表示CN05.1，紅色表示IWRF，藍(lán)色表示NWRFFig.6　(a),(b) and (c)indicates the regional averaged extreme temperature indices,TXx,HWDI and TX90p,respectively,black indicates CN05.1,red indicates IWRF and blue indicates NWRF

3.3　極端高溫事件演變的模擬

按照上文劃分的4個(gè)區(qū)域，計(jì)算了它們區(qū)域平均的模擬和觀測(cè)的氣候指標(biāo)距平的年際變化情況(相較于1982—2011年)。圖7為觀測(cè)和模擬的重慶極端氣溫指標(biāo)的年際變化曲線，在西部地區(qū)，NWRF模擬出的TXx與觀測(cè)差距較大，相關(guān)系數(shù)只有0.23，但是對(duì)于東南部地區(qū)模擬的較好，相關(guān)系數(shù)達(dá)到0.51。對(duì)于熱浪持續(xù)指數(shù)和暖晝指數(shù)，NWRF與觀測(cè)的相關(guān)系數(shù)都在0.6～0.8之間，超過(guò)了0.01信度的顯著性檢驗(yàn)，2006年是重慶地區(qū)有氣象資料記錄以來(lái)最嚴(yán)重的高溫干旱年，無(wú)論是熱浪持續(xù)指數(shù)還是暖晝指數(shù)，NWRF在這4個(gè)區(qū)域都能模擬出該極值，這與使用了較為準(zhǔn)確的強(qiáng)迫場(chǎng)有關(guān)。IWRF不能模擬出極端高溫事件的年際變化特征。

圖7　1982—2014年觀測(cè)和模擬的極端氣溫指標(biāo)距平的年際變化(相對(duì)于1982—2011年)：從上到下依次為日最高氣溫最大值TXx、極端高溫日數(shù)HWDI、暖晝指數(shù)TX90pFig.7　Interannual variations of observed and simulated extreme temperature indicies for 1982—2014 with reference period 1982—2011:(a,b,c,d) TXx; (e,f,g,h) HWDI,(i,j,k,l) TX90p

圖8　1982—2014年夏季重慶極端高溫指數(shù)的趨勢(shì)變化(擴(kuò)大了100倍)，從上到下依次為日最高氣溫最大值TXx(單位：℃)、熱浪持續(xù)指數(shù)(單位：d)、暖晝指數(shù)TX90p(單位：d)；左列為觀測(cè)場(chǎng)的趨勢(shì)分布，第2、3列分別為IWRF和NWRF的趨勢(shì)分布Fig.8　The geographic distribution of the observed (a,d,g) and simulated by IWRF (b,e,h) and NWRF (c,f,i) long-term trends per year for 1982—2014 for extreme temperature indices in the eastern part of Southwest China,(a,b,c),(d,e,f),(g,h,i) represent TXx,HWDI,TX90p,respectively

為了考察這2個(gè)試驗(yàn)對(duì)于重慶極端高溫事件長(zhǎng)期變化趨勢(shì)的模擬能力，分別計(jì)算了觀測(cè)與模擬的這3個(gè)極端氣溫指數(shù)的線性趨勢(shì)。圖8a、8b和8c為TXx的長(zhǎng)期變化趨勢(shì)，從觀測(cè)中可以看到，重慶大部地區(qū)每年的最高氣溫最大值都呈現(xiàn)出了增加趨勢(shì)，這2個(gè)試驗(yàn)均未能模擬出這種一致性的增加趨勢(shì)。圖8d、8e和8f是熱浪持續(xù)指數(shù)的變化趨勢(shì)，觀測(cè)中，四川和重慶東北部地區(qū)的增加趨勢(shì)十分顯著，其余地區(qū)也都呈現(xiàn)出增加的趨勢(shì)，IWRF未能模擬出這種增加趨勢(shì)，在四川東北部和重慶中西部地區(qū)呈減小趨勢(shì)，重慶東部地區(qū)為增加趨勢(shì)，NWRF對(duì)于這種全區(qū)一致的增加趨勢(shì)模擬的較好。圖8g、8h、8i為暖晝指數(shù)的趨勢(shì)分布圖，觀測(cè)中的趨勢(shì)分布和熱浪持續(xù)指數(shù)的分布較為一致，只是增加的趨勢(shì)更為顯著，IWRF同樣未能模擬出這種增加趨勢(shì)，但是可以模擬西部部分地區(qū)和東南、東北部分地區(qū)暖晝指數(shù)的增加趨勢(shì)，NWRF模擬得最好，只是趨勢(shì)值偏小。

4　結(jié)論

本文比較了重慶短期氣候預(yù)測(cè)系統(tǒng)和NCEP/DOE再分析資料驅(qū)動(dòng)的WRF模式對(duì)重慶夏季地表氣溫和極端高溫事件的模擬能力，主要結(jié)論如下：

①2個(gè)試驗(yàn)均能模擬出平均、最高和最低氣溫的空間分布，NWRF模擬的平均氣溫和最高氣溫高值中心的量級(jí)要好于IWRF，但是，IWRF模擬的最低氣溫要優(yōu)于NWRF。它們對(duì)3個(gè)地表氣溫指標(biāo)的概率密度分布都有較好的模擬。總體上，2個(gè)試驗(yàn)對(duì)于重慶夏季氣溫的空間分布模擬的很好。

②它們都能較好的模擬出西部地區(qū)TXx的空間分布，但對(duì)于東部地區(qū)存在較強(qiáng)的冷偏差。IWRF和NWRF模擬的熱浪持續(xù)指數(shù)和暖晝指數(shù)偏少，但是，IWRF的模擬結(jié)果優(yōu)于NWRF。

③NWRF對(duì)于極端氣溫指數(shù)的年際變化特征具有較強(qiáng)的模擬能力，IWRF不能模擬出極端高溫事件的年際變化。

④NWRF能夠模擬出HWDI和TX90p的增加趨勢(shì)，但是不能模擬出TXx的趨勢(shì)變化。IWRF對(duì)這3個(gè)極端高溫指數(shù)的長(zhǎng)期變化趨勢(shì)沒(méi)有模擬能力。

本文的結(jié)論表明IWRF對(duì)于重慶地表氣溫指標(biāo)具有較好的模擬能力，對(duì)于極端高溫事件模擬能力不足，但是，其模擬的熱浪持續(xù)指數(shù)和暖晝指數(shù)的空間分布要好于NWRF，并且模擬的暖晝指數(shù)在西部和東南部地區(qū)誤差較小，這說(shuō)明應(yīng)用該系統(tǒng)進(jìn)行極端高溫事件的預(yù)測(cè)還存在較大的局限性，應(yīng)用該系統(tǒng)預(yù)測(cè)極端高溫事件時(shí)需要格外注意。