• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2018-07-16 11:54:46魏麟驍李永華何卷雄唐紅玉辛?xí)愿?/span>
    中低緯山地氣象 2018年3期
    關(guān)鍵詞:模擬出氣溫偏差

    魏麟驍,李永華,何卷雄,唐紅玉,辛?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)J降哪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í)需要格外注意。

    猜你喜歡
    模擬出氣溫偏差
    氣溫的神奇“魔法”
    基于FY-3D和FY-4A的氣溫時(shí)空融合
    如何走出文章立意偏差的誤區(qū)
    兩矩形上的全偏差
    春 夜
    放射夕陽(yáng)之光
    中華手工(2016年4期)2016-04-20 03:10:35
    關(guān)于均數(shù)與偏差
    與氣溫成反比的東西
    聲音從哪里來(lái)
    能模擬真實(shí)天氣的可視化預(yù)報(bào)器
    大眾健康(2014年12期)2014-01-19 03:15:54
    午夜老司机福利片| 大型av网站在线播放| 精品午夜福利视频在线观看一区| 亚洲一区中文字幕在线| av在线天堂中文字幕 | 这个男人来自地球电影免费观看| 色哟哟哟哟哟哟| 午夜精品久久久久久毛片777| 国产男靠女视频免费网站| av天堂久久9| 国产一区二区在线av高清观看| 99国产精品99久久久久| 久久精品成人免费网站| 国产真人三级小视频在线观看| 高潮久久久久久久久久久不卡| 午夜福利在线免费观看网站| 日本黄色视频三级网站网址| 亚洲五月色婷婷综合| 丰满迷人的少妇在线观看| aaaaa片日本免费| 日韩人妻精品一区2区三区| 亚洲男人天堂网一区| 午夜a级毛片| 亚洲成人精品中文字幕电影 | 激情视频va一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成av片中文字幕在线观看| 国产欧美日韩一区二区三区在线| 国产精品乱码一区二三区的特点 | 别揉我奶头~嗯~啊~动态视频| 欧美+亚洲+日韩+国产| 成人永久免费在线观看视频| 欧美中文综合在线视频| 99久久国产精品久久久| 一级片免费观看大全| 国产精品九九99| 亚洲国产精品999在线| 一级黄色大片毛片| 精品高清国产在线一区| 国产成人啪精品午夜网站| 欧美色视频一区免费| 亚洲久久久国产精品| 色哟哟哟哟哟哟| 亚洲三区欧美一区| 黄色成人免费大全| 岛国在线观看网站| a级毛片黄视频| 曰老女人黄片| 亚洲国产精品合色在线| 搡老熟女国产l中国老女人| 黑人操中国人逼视频| 久久国产亚洲av麻豆专区| 中出人妻视频一区二区| 国产精品爽爽va在线观看网站 | 午夜福利免费观看在线| 欧美老熟妇乱子伦牲交| 亚洲第一av免费看| 亚洲伊人色综图| 色哟哟哟哟哟哟| 国产精品1区2区在线观看.| 精品国产国语对白av| 50天的宝宝边吃奶边哭怎么回事| 成人三级做爰电影| 国产精品亚洲av一区麻豆| 久9热在线精品视频| 一个人观看的视频www高清免费观看 | 超碰97精品在线观看| 欧美日韩av久久| 日韩国内少妇激情av| 50天的宝宝边吃奶边哭怎么回事| 两人在一起打扑克的视频| 午夜福利,免费看| 人人妻,人人澡人人爽秒播| 黄网站色视频无遮挡免费观看| 咕卡用的链子| 久久国产精品人妻蜜桃| 少妇的丰满在线观看| 99热国产这里只有精品6| 日韩人妻精品一区2区三区| 视频在线观看一区二区三区| 国产午夜精品久久久久久| 少妇被粗大的猛进出69影院| 在线观看免费视频网站a站| 久久精品aⅴ一区二区三区四区| 久久狼人影院| 日韩欧美一区视频在线观看| 亚洲片人在线观看| 国产一卡二卡三卡精品| 十八禁网站免费在线| www.999成人在线观看| 国产精品香港三级国产av潘金莲| 午夜免费观看网址| 亚洲熟妇中文字幕五十中出 | 国产单亲对白刺激| 校园春色视频在线观看| 国产成人啪精品午夜网站| 高清在线国产一区| 国产人伦9x9x在线观看| 精品一区二区三区av网在线观看| 搡老乐熟女国产| 一级毛片女人18水好多| 日韩免费高清中文字幕av| 又紧又爽又黄一区二区| 亚洲五月色婷婷综合| 精品国产超薄肉色丝袜足j| ponron亚洲| 亚洲一区二区三区色噜噜 | 新久久久久国产一级毛片| 久久人人精品亚洲av| 一区福利在线观看| 亚洲,欧美精品.| 国产一区二区激情短视频| 中文字幕人妻丝袜制服| 女同久久另类99精品国产91| 亚洲欧美精品综合一区二区三区| 99国产精品一区二区三区| 高清av免费在线| 性少妇av在线| 国产精品秋霞免费鲁丝片| 久久久国产成人精品二区 | 人人妻,人人澡人人爽秒播| 波多野结衣av一区二区av| 亚洲成人免费av在线播放| 久久人人97超碰香蕉20202| av天堂在线播放| 成人三级黄色视频| 亚洲美女黄片视频| 婷婷精品国产亚洲av在线| 国产成人系列免费观看| 免费在线观看完整版高清| 久久香蕉国产精品| 黄色视频,在线免费观看| 久久人妻福利社区极品人妻图片| 狠狠狠狠99中文字幕| 久9热在线精品视频| 精品一区二区三区四区五区乱码| 日本一区二区免费在线视频| av福利片在线| 国产成人一区二区三区免费视频网站| 日韩一卡2卡3卡4卡2021年| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看黄色视频的| 亚洲成a人片在线一区二区| 老司机午夜福利在线观看视频| 十分钟在线观看高清视频www| 成人国语在线视频| 女人被狂操c到高潮| 人妻久久中文字幕网| 正在播放国产对白刺激| 日日爽夜夜爽网站| 久久精品亚洲精品国产色婷小说| 国产亚洲精品综合一区在线观看 | 黑人猛操日本美女一级片| 亚洲精品国产色婷婷电影| 国产成人av激情在线播放| 亚洲av成人一区二区三| 超碰成人久久| 国产熟女xx| 成人精品一区二区免费| 亚洲一区中文字幕在线| 精品少妇一区二区三区视频日本电影| 每晚都被弄得嗷嗷叫到高潮| 巨乳人妻的诱惑在线观看| 黄片小视频在线播放| 免费女性裸体啪啪无遮挡网站| 两性夫妻黄色片| 操美女的视频在线观看| 搡老乐熟女国产| 国产亚洲精品综合一区在线观看 | 999精品在线视频| 欧美+亚洲+日韩+国产| 日本黄色日本黄色录像| 免费搜索国产男女视频| 如日韩欧美国产精品一区二区三区| 午夜精品在线福利| 久久这里只有精品19| 丰满人妻熟妇乱又伦精品不卡| 热99国产精品久久久久久7| 国产亚洲av高清不卡| 悠悠久久av| 搡老岳熟女国产| 国产成年人精品一区二区 | 可以在线观看毛片的网站| 97人妻天天添夜夜摸| 久久国产乱子伦精品免费另类| 免费日韩欧美在线观看| 国产精品99久久99久久久不卡| 激情视频va一区二区三区| 久久亚洲精品不卡| 国产亚洲欧美在线一区二区| 18美女黄网站色大片免费观看| 热99re8久久精品国产| 色尼玛亚洲综合影院| 一级毛片女人18水好多| 搡老熟女国产l中国老女人| 精品一区二区三区视频在线观看免费 | 亚洲va日本ⅴa欧美va伊人久久| 国产精品1区2区在线观看.| 999精品在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美老熟妇乱子伦牲交| 99国产综合亚洲精品| a级毛片在线看网站| 两性夫妻黄色片| 国产成人欧美在线观看| av天堂久久9| 国产熟女xx| 天堂俺去俺来也www色官网| 亚洲精品在线观看二区| av视频免费观看在线观看| 亚洲中文av在线| 交换朋友夫妻互换小说| 又大又爽又粗| 国产片内射在线| 亚洲一卡2卡3卡4卡5卡精品中文| 狂野欧美激情性xxxx| 久久人妻av系列| 久久精品国产综合久久久| 丰满饥渴人妻一区二区三| 欧美日本亚洲视频在线播放| 久久中文看片网| 国产三级黄色录像| 亚洲精品中文字幕一二三四区| 中亚洲国语对白在线视频| 在线av久久热| 丝袜美足系列| 亚洲中文字幕日韩| 久久欧美精品欧美久久欧美| 婷婷六月久久综合丁香| 亚洲精品粉嫩美女一区| 老汉色av国产亚洲站长工具| 少妇裸体淫交视频免费看高清 | 国产在线观看jvid| 成年人免费黄色播放视频| 午夜福利在线观看吧| 狂野欧美激情性xxxx| 亚洲av片天天在线观看| 999精品在线视频| 精品国产一区二区三区四区第35| 久久香蕉国产精品| 亚洲人成电影观看| www.熟女人妻精品国产| 九色亚洲精品在线播放| 亚洲国产欧美日韩在线播放| 欧美黑人精品巨大| 日韩视频一区二区在线观看| 中国美女看黄片| 欧美一级毛片孕妇| 咕卡用的链子| 不卡一级毛片| 99re在线观看精品视频| 女性生殖器流出的白浆| 丝袜人妻中文字幕| 亚洲人成伊人成综合网2020| 少妇的丰满在线观看| 91麻豆av在线| 怎么达到女性高潮| 久久国产精品影院| 国产主播在线观看一区二区| 欧美在线黄色| 99精国产麻豆久久婷婷| 色精品久久人妻99蜜桃| 波多野结衣一区麻豆| 国产区一区二久久| 天堂影院成人在线观看| 一级毛片高清免费大全| 超碰97精品在线观看| 久久中文看片网| 高清毛片免费观看视频网站 | 一区二区三区国产精品乱码| 12—13女人毛片做爰片一| 免费少妇av软件| 极品教师在线免费播放| 免费不卡黄色视频| 又黄又爽又免费观看的视频| 亚洲伊人色综图| 巨乳人妻的诱惑在线观看| 露出奶头的视频| 久久午夜亚洲精品久久| 国产成人啪精品午夜网站| 欧美乱色亚洲激情| 自线自在国产av| 88av欧美| 国产黄色免费在线视频| 亚洲欧美精品综合久久99| 无限看片的www在线观看| 精品人妻1区二区| 一夜夜www| 欧美中文日本在线观看视频| 美国免费a级毛片| 12—13女人毛片做爰片一| 88av欧美| 两性夫妻黄色片| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 久久久久国产精品人妻aⅴ院| 悠悠久久av| 99精品欧美一区二区三区四区| 欧美日韩乱码在线| 日韩人妻精品一区2区三区| 国产色视频综合| 久久午夜亚洲精品久久| 校园春色视频在线观看| 亚洲中文日韩欧美视频| 美女国产高潮福利片在线看| 免费在线观看影片大全网站| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久人妻精品电影| 人妻久久中文字幕网| 午夜福利影视在线免费观看| 国产三级在线视频| 日日干狠狠操夜夜爽| 青草久久国产| 日韩大尺度精品在线看网址 | 国产亚洲欧美98| 国产1区2区3区精品| 国产精品99久久99久久久不卡| 久久国产乱子伦精品免费另类| 亚洲av电影在线进入| 国产欧美日韩一区二区精品| 精品久久久久久久久久免费视频 | 日韩 欧美 亚洲 中文字幕| 日本a在线网址| 欧美日韩亚洲国产一区二区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产欧美日韩av| a在线观看视频网站| 精品久久久久久电影网| 淫妇啪啪啪对白视频| 黄色成人免费大全| 国产精品久久久人人做人人爽| 欧美成人性av电影在线观看| 免费久久久久久久精品成人欧美视频| 色老头精品视频在线观看| 在线观看一区二区三区| 中文字幕人妻丝袜制服| av片东京热男人的天堂| 国产亚洲精品综合一区在线观看 | 乱人伦中国视频| 久久久久久大精品| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品第一综合不卡| 久久中文字幕一级| 操美女的视频在线观看| 欧美人与性动交α欧美软件| 成人三级黄色视频| 男女下面插进去视频免费观看| 欧美日韩亚洲国产一区二区在线观看| 一边摸一边抽搐一进一出视频| 国产一区在线观看成人免费| 成人免费观看视频高清| 日韩有码中文字幕| 久久人妻av系列| 麻豆成人av在线观看| 亚洲激情在线av| 麻豆国产av国片精品| 9191精品国产免费久久| 久久久国产欧美日韩av| 叶爱在线成人免费视频播放| e午夜精品久久久久久久| 精品国产国语对白av| 亚洲美女黄片视频| 亚洲av电影在线进入| 性欧美人与动物交配| 亚洲国产精品合色在线| 黄色女人牲交| 如日韩欧美国产精品一区二区三区| 黄色怎么调成土黄色| 18禁观看日本| videosex国产| www.熟女人妻精品国产| 成人av一区二区三区在线看| 99精品欧美一区二区三区四区| 在线观看午夜福利视频| 国产aⅴ精品一区二区三区波| 在线观看66精品国产| 免费高清在线观看日韩| 午夜久久久在线观看| 亚洲精品中文字幕在线视频| 日韩视频一区二区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 成人av一区二区三区在线看| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 少妇裸体淫交视频免费看高清 | 久久人妻熟女aⅴ| 女人被狂操c到高潮| 在线观看免费高清a一片| 国产野战对白在线观看| 国产有黄有色有爽视频| 国产精品美女特级片免费视频播放器 | 日本三级黄在线观看| 久久久久国产精品人妻aⅴ院| 在线永久观看黄色视频| 午夜精品在线福利| 国产熟女xx| 免费少妇av软件| 久久午夜亚洲精品久久| 亚洲一区二区三区不卡视频| 日韩欧美一区二区三区在线观看| 在线观看免费午夜福利视频| 在线十欧美十亚洲十日本专区| 免费在线观看完整版高清| 三级毛片av免费| 亚洲国产精品合色在线| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 国产免费男女视频| 亚洲视频免费观看视频| 免费在线观看亚洲国产| 在线观看66精品国产| 国产一区二区在线av高清观看| 99久久精品国产亚洲精品| 成年版毛片免费区| 99精品欧美一区二区三区四区| 国产区一区二久久| 午夜两性在线视频| 国产成年人精品一区二区 | 欧美日韩黄片免| 久久人妻福利社区极品人妻图片| 国产aⅴ精品一区二区三区波| 午夜影院日韩av| 欧美日韩视频精品一区| 91九色精品人成在线观看| 婷婷六月久久综合丁香| 国产97色在线日韩免费| 亚洲专区国产一区二区| 国产亚洲欧美在线一区二区| 高清av免费在线| 黄色 视频免费看| 国产亚洲精品久久久久5区| 搡老乐熟女国产| 亚洲成人免费av在线播放| 久久精品成人免费网站| 在线观看免费高清a一片| 一区二区日韩欧美中文字幕| 操出白浆在线播放| 一二三四在线观看免费中文在| 淫秽高清视频在线观看| 伦理电影免费视频| av天堂在线播放| 桃红色精品国产亚洲av| 热re99久久国产66热| 男人舔女人下体高潮全视频| 国产成年人精品一区二区 | 亚洲第一av免费看| 久99久视频精品免费| 国产精品免费视频内射| 国产精品爽爽va在线观看网站 | 美女国产高潮福利片在线看| 日本免费a在线| 韩国av一区二区三区四区| 欧美日本中文国产一区发布| www.熟女人妻精品国产| 村上凉子中文字幕在线| 中文欧美无线码| 人人澡人人妻人| 午夜精品国产一区二区电影| 亚洲性夜色夜夜综合| 9色porny在线观看| 久久人人精品亚洲av| 欧美成人性av电影在线观看| 人人澡人人妻人| 日韩成人在线观看一区二区三区| cao死你这个sao货| 欧美成人免费av一区二区三区| 女性被躁到高潮视频| 国产成年人精品一区二区 | 国产高清视频在线播放一区| 1024视频免费在线观看| 黑人操中国人逼视频| 亚洲,欧美精品.| 中出人妻视频一区二区| 日韩免费av在线播放| 高清毛片免费观看视频网站 | 亚洲成人久久性| 大陆偷拍与自拍| 精品国产一区二区久久| 中文字幕精品免费在线观看视频| 美女大奶头视频| 一个人观看的视频www高清免费观看 | 精品国产一区二区久久| 777久久人妻少妇嫩草av网站| 免费观看人在逋| 免费在线观看视频国产中文字幕亚洲| 亚洲精品粉嫩美女一区| 日本三级黄在线观看| 欧美av亚洲av综合av国产av| 99香蕉大伊视频| 欧美精品亚洲一区二区| 日本免费a在线| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 中亚洲国语对白在线视频| 日日夜夜操网爽| 精品国产一区二区三区四区第35| 久久精品亚洲av国产电影网| av电影中文网址| 两个人免费观看高清视频| 免费在线观看黄色视频的| 在线观看66精品国产| 亚洲av第一区精品v没综合| 啦啦啦在线免费观看视频4| 男女床上黄色一级片免费看| 视频区图区小说| 一边摸一边做爽爽视频免费| 亚洲性夜色夜夜综合| 日日爽夜夜爽网站| 80岁老熟妇乱子伦牲交| 久久精品国产99精品国产亚洲性色 | 麻豆av在线久日| 亚洲av五月六月丁香网| 亚洲av第一区精品v没综合| 中文字幕高清在线视频| 99在线人妻在线中文字幕| 精品久久蜜臀av无| 久久婷婷成人综合色麻豆| a在线观看视频网站| 一边摸一边抽搐一进一出视频| av有码第一页| 久久人人爽av亚洲精品天堂| 午夜福利欧美成人| 亚洲一区二区三区欧美精品| 岛国在线观看网站| 无遮挡黄片免费观看| 精品欧美一区二区三区在线| 日韩精品免费视频一区二区三区| 成人亚洲精品av一区二区 | 久久九九热精品免费| 少妇 在线观看| 男女之事视频高清在线观看| 久久国产亚洲av麻豆专区| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 欧美av亚洲av综合av国产av| e午夜精品久久久久久久| 91成人精品电影| 国产精品美女特级片免费视频播放器 | 麻豆av在线久日| 少妇 在线观看| 中文欧美无线码| 国产一区二区三区视频了| 99香蕉大伊视频| 嫩草影视91久久| 十八禁网站免费在线| 欧美日韩一级在线毛片| 国产精品国产高清国产av| 黑丝袜美女国产一区| 亚洲国产欧美网| 国产蜜桃级精品一区二区三区| 国产高清国产精品国产三级| 91麻豆精品激情在线观看国产 | 女人被躁到高潮嗷嗷叫费观| 中文欧美无线码| 91字幕亚洲| 99久久精品国产亚洲精品| 国产又爽黄色视频| 黄色女人牲交| 久久香蕉国产精品| 制服诱惑二区| 一级毛片精品| ponron亚洲| 午夜老司机福利片| 亚洲欧美激情在线| 交换朋友夫妻互换小说| 美女国产高潮福利片在线看| 满18在线观看网站| 国产伦一二天堂av在线观看| 国产av精品麻豆| 欧美 亚洲 国产 日韩一| 国产精品自产拍在线观看55亚洲| 亚洲av片天天在线观看| 淫秽高清视频在线观看| 一区二区三区国产精品乱码| 女生性感内裤真人,穿戴方法视频| 99国产精品99久久久久| 久久久国产精品麻豆| 香蕉久久夜色| 一个人免费在线观看的高清视频| 亚洲av成人一区二区三| 中文字幕人妻丝袜制服| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 午夜两性在线视频| 国产成人精品久久二区二区91| 午夜91福利影院| 一个人免费在线观看的高清视频| 91精品国产国语对白视频| 老汉色av国产亚洲站长工具| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区三区四区久久 | 亚洲国产精品合色在线| 欧美激情高清一区二区三区| 欧美乱色亚洲激情| 日本免费a在线| 免费在线观看视频国产中文字幕亚洲| 狠狠狠狠99中文字幕| 亚洲一区中文字幕在线| 亚洲一区二区三区欧美精品| 国产亚洲av高清不卡| 男人的好看免费观看在线视频 | 日本免费一区二区三区高清不卡 | 三级毛片av免费| 嫩草影视91久久| 久99久视频精品免费| 亚洲伊人色综图| aaaaa片日本免费| 嫩草影视91久久| 桃色一区二区三区在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久香蕉激情| 亚洲全国av大片| 国产97色在线日韩免费|