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

    MRI-CGCM模式對(duì)東亞夏季風(fēng)的模擬評(píng)估及訂正

    2016-11-28 03:56:43羅連升段春鋒楊瑋徐敏程智丁小俊
    大氣科學(xué) 2016年6期
    關(guān)鍵詞:反氣旋雨帶年際

    羅連升 段春鋒 楊瑋 徐敏 程智 丁小俊

    ?

    MRI-CGCM模式對(duì)東亞夏季風(fēng)的模擬評(píng)估及訂正

    羅連升 段春鋒 楊瑋 徐敏 程智 丁小俊

    安徽省氣候中心,合肥230031

    本文利用CMAP月降水資料、NCEP再分析資料、NOAA的ERSST資料和日本氣象廳海氣耦合模式(MRI-CGCM)的輸出結(jié)果,從東亞夏季風(fēng)氣候態(tài)、主模態(tài)和年際變率等方面分析了MRI-CGCM模式對(duì)東亞夏季風(fēng)的預(yù)測(cè)性能,并且利用觀測(cè)的東亞夏季風(fēng)指數(shù)(EASMI)與模擬PC(principal component)的關(guān)系建立多元線性回歸方程來訂正EASMI(簡(jiǎn)稱PC訂正法)。結(jié)果表明:MRI-CGCM模式能夠較好再現(xiàn)東亞夏季風(fēng)降水和低層風(fēng)場(chǎng)的氣候態(tài),但模擬的西北太平洋反氣旋偏弱、偏東,使得模擬的副熱帶地區(qū)降水量偏小。模式較好地模擬出東亞夏季風(fēng)降水第一模態(tài)(EOF1)及相應(yīng)的低層風(fēng)場(chǎng),能夠較好再現(xiàn)出EOF1對(duì)應(yīng)El Ni?o衰減位相;模擬降水的EOF1與觀測(cè)之間的空間相關(guān)系數(shù)(ACC)為0.72,且能較好地再現(xiàn)其對(duì)應(yīng)的年際變率,其時(shí)間系數(shù)PC1與觀測(cè)之間的相關(guān)系數(shù)為0.41,能模擬出觀測(cè)EOF1的2 a和5 a主導(dǎo)周期;但模擬的我國(guó)以東梅雨鋒區(qū)雨帶位置偏南,這與模擬的西北太平洋反氣旋位置偏南有關(guān)。模式對(duì)降水第二模態(tài)EOF2的模擬能力比EOF1明顯下降,模擬EOF2與觀測(cè)之間的ACC降到0.36;雖然模式能較好地再現(xiàn)出EOF2對(duì)應(yīng)El Ni?o發(fā)展位相,但模擬的西太平洋反氣旋位置偏南,使得雨帶位置偏南,模擬的我國(guó)梅雨鋒區(qū)雨帶位于江南,與觀測(cè)場(chǎng)上江南少雨相反。模式較好地模擬出我國(guó)東部夏季降水和氣溫空間異常分布和年際變化,模擬與觀測(cè)夏季降水和氣溫的多年平均ACC分別為0.74和0.68。模式模擬我國(guó)東部、江淮流域和華南地區(qū)夏季降水多年平均PS評(píng)分分別為69、70和68分,略高于我國(guó)夏季降水業(yè)務(wù)預(yù)測(cè)多年平均評(píng)分(65分)。模擬的我國(guó)東部夏季氣溫與觀測(cè)多年平均PS評(píng)分為74分。PC訂正后EASMI與實(shí)況的相關(guān)系數(shù)由0.51提高到0.65、符號(hào)一致率由84%升到91%、標(biāo)準(zhǔn)差由0.75增大到1.4、大于1個(gè)標(biāo)準(zhǔn)差年數(shù)由6年變?yōu)?2年,訂正后在模擬變幅偏小和梅雨鋒區(qū)雨帶偏南等方面均有一定的改善,對(duì)應(yīng)西太平洋反氣旋位置和梅雨鋒區(qū)雨帶位置與實(shí)況較為吻合。

    東亞夏季風(fēng) 主模態(tài)分析 海氣耦合模式 模式訂正

    1 引言

    東亞夏季風(fēng)對(duì)我國(guó)東部、朝鮮半島以及日本等東亞地區(qū)的天氣和氣候有直接作用(Tao and Chen,1987),季風(fēng)異常引發(fā)的旱澇嚴(yán)重影響著該地區(qū)的經(jīng)濟(jì)發(fā)展(黃榮輝等,2006)。做好東亞夏季風(fēng)的預(yù)測(cè)是提高我國(guó)夏季降水預(yù)測(cè)水平的關(guān)鍵問題之一(陳麗娟等,2013)。目前,動(dòng)力氣候模式是世界各大氣候預(yù)測(cè)業(yè)務(wù)部門的主要工具(賈小龍等,2013),評(píng)估和檢驗(yàn)氣候模式對(duì)東亞夏季風(fēng)的預(yù)測(cè)能力,對(duì)于提高模式產(chǎn)品的業(yè)務(wù)應(yīng)用水平具有重要意義。

    目前已有很多動(dòng)力氣候模式對(duì)東亞夏季風(fēng)模擬能力進(jìn)行評(píng)估和檢驗(yàn)的研究。例如CMIP3和CMIP5的多個(gè)海氣耦合模式。研究表明CMIP3的多個(gè)海氣耦合模式基本上都能模擬出東亞夏季風(fēng)降水的氣候態(tài),但多數(shù)模式對(duì)東亞夏季風(fēng)降水年際變化的預(yù)測(cè)能力有待于提高,且普遍存在模擬值變幅偏小的缺點(diǎn)(姜大膀等,2004;張莉等,2008;孫穎和丁一匯,2008)。He and Zhou(2014)、何超等(2012)系統(tǒng)分析了中國(guó)科學(xué)院大氣物理研究所氣候系統(tǒng)模式FGOALS和CMIP5氣候模式對(duì)西太平洋副熱帶高壓(簡(jiǎn)稱西太副高)及其年際變率的模擬能力,研究表明模式能夠很好地模擬出西太副高年際變率第一模態(tài)及相應(yīng)的海溫強(qiáng)迫,但對(duì)第二模態(tài)的模擬能力明顯下降,且模擬的第二模態(tài)年際變幅偏小和位置偏南,對(duì)第二模態(tài)模擬不好的原因在于模式低估了西太副高與西太平洋地區(qū)海溫(SST)之間的局地海氣相互作用。Song and Zhou(2014a, 2014b)通過對(duì)CMIP5單獨(dú)大氣模式和海氣耦合模式的比較,表明海氣耦合模式對(duì)東亞夏季風(fēng)氣候態(tài)、年際變率和梅雨鋒區(qū)雨帶的模擬能力好于大氣模式,指出了海氣耦合過程在模擬東亞夏季風(fēng)年際變率中的重要作用。Zou and Zhou(2013)借助區(qū)域海氣耦合模式也證明了海氣耦合過程在模擬東亞夏季風(fēng)年際變率中的重要作用。Song et al.(2014)研究表明CMIP5模式能夠合理再現(xiàn)東亞夏季風(fēng)的年代際變化。上述研究也表明多模式集合平均預(yù)報(bào)效果普遍好于單個(gè)模式。

    雖然多模式集合預(yù)報(bào)是目前公認(rèn)的提高模式預(yù)測(cè)能力的有效方法,但前提是參加多模式集合的各氣候模式要有較好的預(yù)測(cè)性能(丁一匯,2011)。因此必須要充分了解氣候業(yè)務(wù)中應(yīng)用的各個(gè)氣候模式的預(yù)測(cè)性能。目前世界幾大氣候業(yè)務(wù)單位各自研發(fā)的海氣耦合模式,如國(guó)家氣候中心的BCC_ CM1.0、NCEP的CFSv2和日本氣象廳的MRI- CGCM廣泛應(yīng)用于我國(guó)短期氣候預(yù)測(cè)業(yè)務(wù)。研究表明這三個(gè)海氣耦合模式基本上能模擬出東亞夏季風(fēng)的氣候態(tài),對(duì)東亞夏季風(fēng)降水年際變化也要一定的預(yù)測(cè)能力(Rajendran et al.,2004; 司東等,2009;Gao et al.,2011;董敏等,2013;吳統(tǒng)文等,2013;Jiang et al.,2013),其中MRI-CGCM模式能較好地再現(xiàn)印度和西北太平洋地區(qū)季風(fēng)進(jìn)程。多數(shù)學(xué)者主要評(píng)估了BCC_CM1.0模式和CFSv2模式對(duì)東亞夏季風(fēng)的模擬能力,而對(duì)MRI-CGCM模式評(píng)估較少。因此,有必要開展MRI-CGCM模式對(duì)東亞夏季風(fēng)的預(yù)報(bào)性能的評(píng)估,給出該模式作為夏季預(yù)測(cè)的可靠性程度,有利于該模式在夏季預(yù)測(cè)中更好地使用。本文通過對(duì)比分析實(shí)況和MRI-CGCM模式模擬的東亞夏季風(fēng)降水的主要模態(tài)及其年際變化,評(píng)估模式對(duì)東亞夏季風(fēng)的模擬能力,在此基礎(chǔ)上,探討東亞夏季風(fēng)模式預(yù)測(cè)訂正的可能性,希望能夠?yàn)闁|亞夏季風(fēng)的預(yù)測(cè)提供有用的信息。

    2 模式、資料和方法簡(jiǎn)介

    本文所用的氣候模式是日本氣象廳海氣耦合模式MRI-CGCM(Yukimoto et al., 2001),由全球大氣環(huán)流模式(AGCM)和全球海洋環(huán)流模式(OGCM)耦合形成。AGCM模式水平方向采用三角形截?cái)啵═L95),水平分辨率近似于1.875°×1.875°,垂直方向分成40層(模式頂層為0.4 hPa)。OGCM模式是日本氣象研究所區(qū)域海洋模式(MRI.COM),在75°S~75°N地區(qū)水平分辨率在經(jīng)向?yàn)?.0°、緯向?yàn)?.3°~1.0°,垂直方向分成50層。模式輸出產(chǎn)品統(tǒng)一處理成網(wǎng)格點(diǎn)數(shù)據(jù),水平分辨率為2.5°×2.5°。本文利用模式以3月為初始場(chǎng)來回報(bào)的3~9月數(shù)據(jù),文中用到夏季(6~8月)的數(shù)據(jù),回報(bào)時(shí)段為1979~2010年,共有10個(gè)成員,文中提到的模擬場(chǎng)均為10個(gè)成員的平均。

    用來評(píng)估模式模擬能力的對(duì)比資料有:(1)CMAP月降水資料(Xie and Arkin,1997),水平分辨率為2.5°×2.5°;(2)GPCP月降水資料(Adler et al.,2003),水平分辨率為2.5°×2.5°;(3)NCEP/ NCAR再分析資料,水平分辨率為2.5°×2.5°(Kistler et al.,1996);(4)NOAA月平均的ERSST海溫資料,水平分辨率為2°×2°(Smith et al.,2008);(5)Ni?o3.4指數(shù),即(5°S~5°N,120°~170°W)范圍內(nèi)海表面溫度異常(SSTA)的區(qū)域平均值,由ERSST海溫資料計(jì)算得到。所有資料均使用1979~2010年時(shí)段進(jìn)行分析,為了便于討論,本文將再分析資料和觀測(cè)資料統(tǒng)稱為“觀測(cè)資料”。為了便于分析,將CMAP、GPCP和ERSST插值到全球144×73個(gè)(2.5°×2.5°)網(wǎng)格上。夏季數(shù)據(jù)均由月數(shù)據(jù)處理得到。文中冬季指的是前一年12月到當(dāng)年2月。

    本文所用的東亞季風(fēng)區(qū)為10°N~50°N、110°E~150°E,區(qū)域內(nèi)總格點(diǎn)數(shù)為289。文中評(píng)估用到的我國(guó)東部地區(qū)為(20°N~45°N,110°E~120°E),區(qū)域內(nèi)總格點(diǎn)數(shù)為55;華南地區(qū)為(20°N~27.5°N,110°E~120°E),區(qū)域內(nèi)總格點(diǎn)數(shù)為20;江淮流域?yàn)椋?7.5°N~35°N,110°E~120°E),區(qū)域內(nèi)格點(diǎn)數(shù)為20。

    本文用來評(píng)估模式的模擬性能指標(biāo)有距平符號(hào)一致率、時(shí)間相關(guān)系數(shù)TCC、空間相關(guān)系數(shù)ACC和國(guó)家氣候中心的預(yù)測(cè)技巧評(píng)分(簡(jiǎn)稱PS評(píng)分)。

    文中采用的PS評(píng)分為短期氣候預(yù)測(cè)業(yè)務(wù)常用方法(劉長(zhǎng)征等,2013),在距平符號(hào)預(yù)報(bào)準(zhǔn)確百分率的基礎(chǔ)上增加異常級(jí)加權(quán)得分,能夠較好地反映整體趨勢(shì)預(yù)測(cè)能力。

    本文用經(jīng)驗(yàn)正交函數(shù)(EOF)分解方法得到模擬和觀測(cè)的東亞地區(qū)夏季降水主要模態(tài)。

    3 結(jié)果分析

    能夠合理地模擬出觀測(cè)氣候態(tài)的特征是衡量模式性能的重要條件。由于觀測(cè)降水資料本身有很大的不確定性,尤其是在海洋上(Collins et al., 2013; He and Zhou, 2015),因此本文利用CMAP和GPCP兩種不同觀測(cè)降水資料來評(píng)估模式對(duì)東亞地區(qū)降水的預(yù)測(cè)性能。圖1是東亞季風(fēng)區(qū)1979~2010年32年平均的觀測(cè)和模擬的夏季平均日降水量的比較。CMAP觀測(cè)降水量在東亞季風(fēng)區(qū)從北到南呈緯向遞增,在40°N以南大部地區(qū)降水量大于4 mm d?1,大值中心位于熱帶地區(qū)(圖1a)。GPCP觀測(cè)降水量也表現(xiàn)出與CMAP相似的特征(圖1b),兩者在大陸上差異不大,但在海洋上差異較大,尤其在熱帶西太平洋地區(qū),CMAP降水量比GPCP偏大1~5 mm d?1,與Zhou et al.(2008)分析結(jié)果一致。模式較好地模擬出了東亞夏季風(fēng)降水的分布特征(圖1c)。從模擬降水量與CMAP(圖1d)和GPCP觀測(cè)降水量(圖1e)的差值來看,也是在熱帶西太平洋地區(qū)差異最大,模擬降水量與CMAP觀測(cè)降水量差值明顯小于與GPCP觀測(cè)降水量的差值,模式在海洋上模擬的降水量與CMAP觀測(cè)更為接近,因此本文以下分析用到的實(shí)況降水資料均為CMAP降水資料。從模擬降水量與實(shí)況降水量的差值來看(圖1d),東亞大部分地區(qū)降水量差異在1 mm d?1范圍內(nèi),而且東亞地區(qū)降水量差值大致呈“? + ?”緯向分布,我國(guó)華南地區(qū)、菲律賓附近及其以東熱帶西太平洋地區(qū)和我國(guó)黃淮地區(qū)到韓國(guó)、日本南部模擬降水量偏小1~5 mm d?1,而我國(guó)江南和20°N附近西太平洋地區(qū)模擬值偏大1~3 mm d?1。觀測(cè)850 hPa風(fēng)場(chǎng)在南海盛行偏西風(fēng),西北太平洋為反氣旋控制,我國(guó)東部主要受該反氣旋西北側(cè)西南氣流控制。模式較好地模擬出了低層風(fēng)場(chǎng)的上述主要特征(圖1c)。在西北太平洋地區(qū),模擬850 hPa風(fēng)場(chǎng)與實(shí)況差值場(chǎng)上表現(xiàn)為氣旋差異,即模擬的西北太平洋反氣旋偏東、偏弱,導(dǎo)致我國(guó)大陸到朝鮮半島、日本一帶水汽輸送偏弱,降水量偏少,同時(shí),南海和西太平洋這些降水量偏多的區(qū)域均位于氣旋性差異的南側(cè)。

    圖1 1979~2010年平均的夏季日降水量(陰影,單位:mm d?1)與850 hPa風(fēng)場(chǎng)(箭頭,單位:m s?1):(a)CMAP資料;(b)GPCP資料;(c)MRI-CGCM模式模擬;(d)MRI-CGCM模式模擬與CMAP資料差值;(e)MRI-CGCM模式模擬與GPCP資料差值

    由以上分析可知,MRI-CGCM模式能較好地再現(xiàn)東亞夏季風(fēng)降水和低層風(fēng)場(chǎng)的氣候態(tài)特征。

    3.2 東亞夏季降水主要模態(tài)的模擬

    對(duì)東亞季風(fēng)區(qū)1979~2010年觀測(cè)夏季降水進(jìn)行EOF分析,前2個(gè)模態(tài)的方差貢獻(xiàn)分別為37.7%和10.2%,占總方差貢獻(xiàn)47.9%。其中第一模態(tài)(EOF1)的空間分布特征是在西太平洋地區(qū)和副熱帶地區(qū)(25°N~40°N)夏季降水為反位相變化,當(dāng)西太平洋夏季降水量偏少(多)時(shí),從長(zhǎng)江中下游到韓國(guó)及日本一帶副熱帶地區(qū)降水量偏多(少)(圖2a)。這兩個(gè)地區(qū)降水是東亞夏季風(fēng)降水年際變率主要模態(tài)即偶極子模態(tài)降水的中心(Wang et al., 2008)。EOF1對(duì)應(yīng)的時(shí)間系數(shù)表現(xiàn)為明顯的年際變化特征(圖2e虛線),其中準(zhǔn)5年周期最為明顯,其次是準(zhǔn)2年周期,均通過了紅噪音譜的95%置信上限(圖2g)。第二模態(tài)(EOF2)的特點(diǎn)是東亞地區(qū)從南到北大致為“+?+”三極子分布,當(dāng)15°N以南熱帶地區(qū)和南海、我國(guó)華南降水量偏多(少)時(shí),同時(shí)我國(guó)江淮、黃淮到韓國(guó)、日本一帶副熱帶地區(qū)降水量也偏多(少),而這兩個(gè)多雨區(qū)中間的西北太平洋大部及我國(guó)江南降水量偏少(多)(圖2c)。EOF2對(duì)應(yīng)的時(shí)間系數(shù)也呈現(xiàn)出明顯的年際變化(圖2f虛線),也是準(zhǔn)5年和準(zhǔn)2年振蕩明顯(圖2h)。本文EOF分析結(jié)果與Wu et al.(2009)對(duì)東亞夏季降水進(jìn)行季節(jié)EOF得到的模態(tài)基本一致。

    同樣對(duì)同一時(shí)段模式模擬東亞夏季降水進(jìn)行EOF分析,前2個(gè)模態(tài)的方差貢獻(xiàn)分別為39%和19%,占總方差貢獻(xiàn)58%,模式模擬的EOF1方差貢獻(xiàn)與實(shí)況幾乎一樣,EOF2方差貢獻(xiàn)比實(shí)況大9%。模式能夠較好地抓住東亞夏季降水第一個(gè)模態(tài)的空間分布特征,模擬的EOF1與觀測(cè)之間的空間相關(guān)系數(shù)為0.72,通過了0.01的顯著水平檢驗(yàn)。在CMIP5的34個(gè)海氣耦合模式中,MRI-CGCM對(duì)東亞夏季降水EOF1的模擬能力排在前6名(Song and Zhou, 2014b)。模式很好地模擬出EOF1我國(guó)梅雨鋒區(qū)雨帶的位置和強(qiáng)度,只是我國(guó)以東梅雨鋒區(qū)雨帶偏窄和位置偏南(圖2c)。模式能較好地模擬出EOF1對(duì)應(yīng)的年際變率,對(duì)應(yīng)的時(shí)間系數(shù)PC1與觀測(cè)之間的相關(guān)系數(shù)為0.41(圖2e),通過了0.05的顯著水平檢驗(yàn),且能較好地再現(xiàn)出觀測(cè)EOF1的5 a和2 a主導(dǎo)周期(圖略)。相對(duì)于EOF1,模式對(duì)EOF2的模擬能力明顯下降,模擬的EOF2與觀測(cè)之間的ACC降到0.36,其對(duì)應(yīng)PC2與觀測(cè)之間的相關(guān)系數(shù)為0.35,通過了0.05的顯著水平檢驗(yàn),但模式僅能再現(xiàn)出觀測(cè)EOF2的2 a主導(dǎo)周期,未能模擬出5 a主導(dǎo)周期,表明模式對(duì)EOF2的年際變率模擬能力也差于EOF1。模式模擬的EOF2三極子模態(tài)明顯變窄且偏南,西北太平洋雨帶和副熱帶地區(qū)雨帶位置比觀測(cè)偏南,模擬的梅雨鋒區(qū)雨帶位置南移到江南,與觀測(cè)場(chǎng)上江南降水量偏少正好相反,即模擬的雨帶位置與觀測(cè)差異較大(圖2b、d)。模式模擬的東亞副熱帶地區(qū)雨帶偏南是目前大部分海氣耦合模式普遍存在的缺點(diǎn)(Song and Zhou, 2014a)。

    圖2 1979~2010年(a、b)觀測(cè)東亞季風(fēng)區(qū)夏季降水量與(c、d)模擬的夏季降水量EOF分解的前2個(gè)模態(tài)EOF1、EOF2及其(e、f)對(duì)應(yīng)的時(shí)間系數(shù)PC1、PC2和(g、h)相應(yīng)的功率譜。圖g、h中虛線為紅噪音的95%信度上限

    為了分析東亞夏季降水前兩個(gè)模態(tài)對(duì)應(yīng)的環(huán)流和SST分布,文中對(duì)模擬和觀測(cè)降水EOF分解對(duì)應(yīng)的時(shí)間系數(shù)PC1、PC2與同期低層850 hPa風(fēng)場(chǎng)及SST做了相關(guān)分析,空間分布見圖3。觀測(cè)場(chǎng)EOF1在熱帶印度洋和西太平洋低層夏季盛行東風(fēng),東亞地區(qū)低層風(fēng)場(chǎng)表現(xiàn)為西北太平洋、日本中部和鄂霍次克海附近的“反氣旋—?dú)庑礆庑苯?jīng)向波列(圖3a),這正是典型的東亞—太平洋(EAP)遙相關(guān)型波列(黃榮輝,1990)。其中西北太平洋反氣旋決定了東亞夏季降水的雨帶分布,該反氣旋偏強(qiáng)導(dǎo)致梅雨鋒區(qū)夏季降水量偏多、西北太平洋地區(qū)降水量偏少。西北太平洋反氣旋異常是聯(lián)系ENSO與東亞季風(fēng)環(huán)流的重要紐帶。在夏季觀測(cè)海溫場(chǎng)上,印度洋和西太平洋SST偏暖,而中東太平洋海溫趨于變冷。已有研究表明冬季El Ni?o事件對(duì)印度洋SST有增暖作用,印度洋是聯(lián)系冬季ENSO與夏季西北太平洋反氣旋異常的橋梁(Xie et al., 2009)。本文分析了觀測(cè)PC1與Ni?o3.4指數(shù)的超前滯后相關(guān)(圖4a虛線),圖4中橫坐標(biāo)(?1)、(0)和(1)分別表示前一年、當(dāng)年和后一年,由圖可見,PC1與前一年夏秋季以及前期冬春季Ni?o3.4指數(shù)的相關(guān)系數(shù)均超過0.4,通過了0.05的顯著水平檢驗(yàn),然而與同期夏季及后期秋冬季Ni?o3.4指數(shù)為弱負(fù)相關(guān)(未通過信度檢驗(yàn)),表明EOF1模態(tài)出現(xiàn)在El Ni?o事件峰值之后,也就是說EOF1模態(tài)反映的是El Ni?o衰減年東亞夏季風(fēng)環(huán)流的情況,這與Wu et al.(2009)等分析結(jié)果是一致的。夏季熱帶印度洋海溫增暖激發(fā)的大氣Kelvin波導(dǎo)致了西北太平洋反氣旋的產(chǎn)生,從而影響東亞夏季風(fēng)降水的雨帶分布(Wu et al., 2009; Song and Zhou, 2014a, 2014b)。

    模式能較好地模擬EOF1低層風(fēng)場(chǎng)EAP波列特征,但模擬的西北太平洋和日本中部的反氣旋—?dú)庑恢闷希▓D3c),導(dǎo)致模擬的我國(guó)以東梅雨鋒區(qū)雨帶位置偏南。模式能合理地再現(xiàn)EOF1在印度洋和太平洋SST分布(圖3c),模擬PC1與Ni?o3.4相關(guān)曲線變化與觀測(cè)PC1基本一致(圖4),只是模擬相關(guān)曲線變幅大于觀測(cè),表明模式對(duì)SST的模擬更為理想化,模擬的夏季中東太平洋海溫變冷幅度大于實(shí)況(圖3c)。

    觀測(cè)EOF2在赤道太平洋低層盛行西風(fēng),東亞地區(qū)從南到北在南海、西北太平洋和日本海附近為“氣旋—反氣旋—?dú)庑苯?jīng)向配置,南海氣旋導(dǎo)致華南地區(qū)夏季降水量偏多,同時(shí)看出西北太平洋反氣旋位置比EOF1偏北偏東,中心位置位于27°N附近,控制我國(guó)江南地區(qū),使得梅雨鋒區(qū)雨帶位于我國(guó)江淮、黃淮及其以東地區(qū),雨帶西邊界比EOF1偏東,位置也比EOF1偏北,而江南降水量偏少(圖2b)。在夏季觀測(cè)海溫場(chǎng)上最明顯的特征是熱帶中東太平洋SST偏暖(圖3b)。從觀測(cè)PC2與Ni?o3.4指數(shù)的相關(guān)來看(圖4),PC2與前期秋冬季為弱負(fù)相關(guān),到前期春季轉(zhuǎn)為正相關(guān),但未通過顯著性檢驗(yàn),到同期夏季相關(guān)系數(shù)陡升到0.58,與后期秋冬季及下一年春季均為明顯的正相關(guān),最大相關(guān)系數(shù)出現(xiàn)在同期夏季,表明EOF2對(duì)應(yīng)的是El Ni?o發(fā)展位相,這與黃榮輝和陳文(2002)指出在El Ni?o發(fā)展階段我國(guó)江淮流域到韓國(guó)日本一帶夏季降水量偏多而江南降水量偏少的結(jié)論相一致。研究表明對(duì)應(yīng)東亞夏季風(fēng)年際變率EOF分解第二模態(tài),赤道中東太平洋正的加熱異常在熱帶西太平洋激發(fā)出異常氣旋性環(huán)流(Wu el al., 2009),從而導(dǎo)致西北太平洋反氣旋位置比EOF1模態(tài)偏北,使得我國(guó)江淮流域到韓國(guó)、日本一帶降水量偏多,而江南降水量偏少。

    模式基本能模擬出了EOF2在熱帶太平洋和東亞地區(qū)的低層風(fēng)場(chǎng)特征(圖3d),但模擬的西北太平洋反氣旋位置比觀測(cè)明顯偏南,使得模擬的華南雨帶南移到南海地區(qū),相應(yīng)地模擬的梅雨鋒區(qū)雨帶南移到江南地區(qū)。模式能較好地再現(xiàn)夏季熱帶中東太平洋海溫偏暖的特征(圖3d),模擬PC2與Ni?o3.4相關(guān)變化曲線與觀測(cè)基本一致(圖4b),只是前期冬春季相關(guān)系數(shù)明顯大于觀測(cè),即模擬夏季熱帶中東太平洋海溫偏暖幅度大于觀測(cè)(圖3d)。

    圖3 圖2中的PC1、PC2與(a、b)同期觀測(cè)的、(c、d)模擬的850 hPa風(fēng)場(chǎng)和SST的相關(guān)。箭頭表示850 hPa風(fēng)場(chǎng)異常,陰影區(qū)表示與SST的相關(guān)系數(shù)(通過了0.1的顯著水平t檢驗(yàn)),字母“C”、“A”分別表示“氣旋”、“反氣旋”

    圖4 圖2中的(a)PC1和(b)PC2與Ni?o3.4指數(shù)超前滯后相關(guān)系數(shù)。橫坐標(biāo)中(?1)、(0)和(1)分別表示前一年、當(dāng)年和后一年,兩條水平虛線表示95%信度水平線,JJA、SON、DJF、MAM分別表示6、7、8月,9、10、11月,12、1、2月,3、4、5月

    3.3 不同區(qū)域夏季降水和氣溫的預(yù)測(cè)檢驗(yàn)

    從1979~2010年觀測(cè)與模擬的夏季降水量的時(shí)間相關(guān)系數(shù)(TCC)分布(圖5a)來看,在東亞熱帶地區(qū)TCC大于0.2,大部分地區(qū)通過了0.1的顯著水平檢驗(yàn),而副熱帶地區(qū)TCC明顯小于熱帶地區(qū),模式對(duì)熱帶地區(qū)的預(yù)報(bào)能力明顯高于副熱帶地區(qū)。在我國(guó)東部地區(qū),黃河以南大部地區(qū)TCC為正值,其中黃淮南部地區(qū)通過了0.1的顯著水平檢驗(yàn)。MRI-CGCM模式在我國(guó)東部的TCC分布與4個(gè)海氣耦合模式(CFSv2、GFDL CM2.1、ABOM POAMA2.4、FRCGC SINTEX-F)集合平均回報(bào)同一時(shí)段夏季降水在我國(guó)東部的TCC分布(Wang et al., 2015)基本一致,表明MRI-CGCM模式在我國(guó)東部夏季降水的預(yù)報(bào)能力與多模式集合平均預(yù)報(bào)能力相當(dāng)。從1979~2010年觀測(cè)與模擬的夏季氣溫的TCC分布(圖5b)來看,東亞地區(qū)TCC均為正值,其中熱帶地區(qū)大部大于0.3,通過了0.1的顯著水平檢驗(yàn);在副熱帶地區(qū),120°E以東大部TCC也通過了0.1的顯著水平檢驗(yàn),表明模式對(duì)東亞地區(qū)夏季氣溫具有較好的預(yù)測(cè)能力。由圖5也可看出,MRI-CGCM對(duì)東亞地區(qū)夏季氣溫的模擬能力明顯好于夏季降水。

    圖5 模擬的1979~2010年夏季(a)降水量、(b)氣溫與觀測(cè)的降水量、氣溫之間的時(shí)間相關(guān)系數(shù)分布。陰影區(qū)表示通過了0.1的顯著水平t檢驗(yàn)

    我國(guó)東部1979~2010年觀測(cè)與模擬的夏季降水量的ACC多年平均為0.74(圖6),通過了0.01的顯著水平檢驗(yàn),高于目前國(guó)家氣候中心預(yù)測(cè)業(yè)務(wù)上應(yīng)用的多模式解釋應(yīng)用集成方法(包含BCC_CM1.0、CFSv2、MRI-CGCM和歐洲中期天氣預(yù)報(bào)中心的System4等4個(gè)氣候業(yè)務(wù)模式)回報(bào)2003~2012年我國(guó)夏季降水最佳ACC(0.05)(劉長(zhǎng)征等,2013)。從每年降水ACC變化來看,每年ACC均通過了0.01的顯著水平檢驗(yàn),其中最小值為0.58(2003年),最大值為0.89(1999年)。模式模擬我國(guó)東部夏季氣溫與觀測(cè)之間的32年平均ACC為0.68(圖6),通過了0.01的顯著水平檢驗(yàn),其中ACC最小是0.48(2007年),最大是0.79(1979年)。表明模式對(duì)我國(guó)東部夏季降水和氣溫的空間異常分布型具有較好的模擬能力。

    圖6 模擬的1979~2010年我國(guó)東部夏季降水量、氣溫與觀測(cè)之間的空間相關(guān)系數(shù),圖中虛線表示99%信度水平線

    我國(guó)東部、華南地區(qū)和江淮流域模擬與觀測(cè)降水多年平均距平符號(hào)一致率分別為53%、53%和58%(圖7a),其中模式對(duì)江淮流域的預(yù)報(bào)技巧最高。1979~2010年中江淮流域距平符號(hào)一致率大于75%有9年(圖7a),分別為:1981、1985、1986、1993、1996、1997、1998、2002和2007年,這9年中除了1981年和1985年為正常年外,其他年份均為ENSO年,也就是說ENSO年模式對(duì)江淮流域夏季降水的預(yù)測(cè)能力較強(qiáng)。我國(guó)東部、華南地區(qū)和江淮流域模擬夏季降水32年平均PS評(píng)分分別為69、68和70分(圖7b),略高于我國(guó)夏季降水業(yè)務(wù)預(yù)測(cè)多年平均(65分)(艾秀等,2008;2015年國(guó)家氣候中心氣象現(xiàn)代化實(shí)施方案)。從各年P(guān)S評(píng)分來看,我國(guó)東部、華南地區(qū)和江淮流域PS評(píng)分超過70分的分別有16年、20年和19年,即各個(gè)區(qū)域32年中有一半以上年份預(yù)測(cè)評(píng)分超過70分,MRI-CGCM模式對(duì)我國(guó)各區(qū)域夏季降水年際變化具有較好的預(yù)測(cè)能力。模式模擬的1979~2010年我國(guó)東部夏季氣溫平均符號(hào)一致率和PS評(píng)分為62%和74分(圖略)。從各年P(guān)S評(píng)分來看,PS評(píng)分超過70分的有22年,即32年中有三分之二的年份預(yù)測(cè)評(píng)分超過70分(圖略),模式對(duì)我國(guó)東部夏季氣溫年際變化也具有較好的預(yù)測(cè)能力。

    圖7 模式模擬中國(guó)各區(qū)域夏季降水量和觀測(cè)降水量的(a)距平符號(hào)一致率和(b)PS評(píng)分

    3.4 模擬的東亞夏季風(fēng)指數(shù)的訂正

    Wang et al.(2008)對(duì)目前25種東亞夏季風(fēng)指數(shù)進(jìn)行了評(píng)估,指出Wang定義的東亞夏季風(fēng)指數(shù)(Wang and Fan,1999)和張慶云定義的東亞夏季風(fēng)指數(shù)(張慶云等,2003)與東亞夏季風(fēng)的多變量EOF分解PC1的相關(guān)系數(shù)分別為?0.97和?0.93,是25種指數(shù)中相關(guān)系數(shù)最高的前兩種。這兩個(gè)指數(shù)定義相似,均用東亞熱帶地區(qū)和副熱帶地區(qū)的緯向風(fēng)南北梯度來定義,反映東亞熱帶輻合帶(季風(fēng)槽)與副熱帶輻合帶(梅雨鋒)的強(qiáng)度存在反位相變化特征。兩者的差異只是選取的區(qū)域不同而已。這兩個(gè)指數(shù)物理意義明確且計(jì)算簡(jiǎn)單,兩者1979~2010年相關(guān)系數(shù)達(dá)0.93。我國(guó)短期氣候預(yù)測(cè)業(yè)務(wù)中普遍采用張慶云定義的指數(shù)(簡(jiǎn)稱EASMI)來表征東亞夏季風(fēng)強(qiáng)弱,因此本文也采用EASMI來評(píng)估模式對(duì)東亞夏季風(fēng)年際變化的模擬能力。EASMI是用東亞熱帶地區(qū)(10°~20°N,100°~150°E)與副熱帶地區(qū)(25°~35°N,100°~150°E)夏季平均的850 hPa風(fēng)場(chǎng)的緯向風(fēng)距平之差來表示。當(dāng)指數(shù)小(大)表明東亞夏季風(fēng)弱(強(qiáng)),東亞熱帶季風(fēng)槽區(qū)夏季降水量偏少(多),而梅雨鋒區(qū)降水量偏多(少)。

    從觀測(cè)EASMI與同期降水和850 hPa風(fēng)場(chǎng)的相關(guān)分布來看,當(dāng)東亞夏季風(fēng)偏弱(強(qiáng))時(shí),西太平洋地區(qū)為異常反氣旋(氣旋),受其影響,東亞夏季降水從南到北為偶極子經(jīng)向分布,西太平洋地區(qū)夏季降水量偏少(多),而我國(guó)江淮流域到韓國(guó)、日本一帶梅雨鋒區(qū)降水量偏多(少)(圖8a)。比較圖8a和圖2a可見,觀測(cè)EASMI強(qiáng)弱對(duì)應(yīng)的東亞夏季降水分布與觀測(cè)東亞夏季降水EOF1模態(tài)基本一致,1979~2010年觀測(cè)EASMI與PC1相關(guān)系數(shù)為0.84,即EASMI反映的是東亞夏季降水的EOF1模態(tài)。

    圖8 1979~2010年(a)觀測(cè)的EASMI與觀測(cè)的夏季降水量、(b)模擬的EASMI與模擬的降水量、(c)訂正的EASMI與觀測(cè)的夏季降水量和850 hPa風(fēng)場(chǎng)異常(箭頭,單位:m s?1)的相關(guān)關(guān)系分布。陰影區(qū)表示通過0.1的顯著水平t檢驗(yàn)

    1979~2010年模擬和觀測(cè)的EASMI變化曲線基本一致(圖9),相關(guān)系數(shù)為0.51,通過了0.01的顯著水平檢驗(yàn),表明MRI-CGCM模式對(duì)東亞夏季風(fēng)預(yù)測(cè)技巧較高,較好地模擬出東亞夏季風(fēng)的年際變化,但模擬的東亞夏季風(fēng)指數(shù)變幅偏小。模式EASMI能模擬出西太平洋反氣旋和東亞夏季風(fēng)降水偶極子經(jīng)向分布型,當(dāng)模擬東亞夏季風(fēng)偏弱(強(qiáng))時(shí),西太平洋地區(qū)降水量偏少(多),而梅雨鋒區(qū)降水量偏多(少),但是模擬的西太平洋反氣旋比實(shí)況偏南,導(dǎo)致模擬的梅雨鋒區(qū)雨帶位置偏南(圖8b)。

    圖9 1979~2010年觀測(cè)、模擬和訂正的EASMI

    模擬EASMI對(duì)應(yīng)的梅雨鋒區(qū)雨帶位置偏南,那么是否可以對(duì)模擬的EASMI進(jìn)行訂正,使得模擬的梅雨鋒區(qū)雨帶位置與實(shí)況更為接近呢。本文嘗試?yán)糜^測(cè)的EASMI與模擬模態(tài)對(duì)應(yīng)的PC建立多元回歸方程來訂正EASMI。把觀測(cè)的EASMI作為預(yù)報(bào)變量,東亞夏季降水前4個(gè)模態(tài)對(duì)應(yīng)的PC作為預(yù)報(bào)因子,建立多元線性回歸方程,文中把這種訂正方法稱為PC訂正法。文中選前4個(gè)PC的原因是前4個(gè)PC與觀測(cè)的EASMI相關(guān)性較好,相關(guān)系數(shù)分別為0.57、?0.16、0.22和?0.28,而第5模態(tài)之后的PC與觀測(cè)的EASMI的相關(guān)性明顯下降。根據(jù)PC訂正方法利用交叉檢驗(yàn)計(jì)算出1979~2010年訂正的EASMI。下面簡(jiǎn)要介紹訂正方案:例如,要訂正第年的EASMI,通過多元線性回歸方程:

    表1是1979~2010年訂正前后EASMI與觀測(cè)EASMI的比較,訂正指數(shù)與實(shí)況相關(guān)系數(shù)明顯增大,由訂正前0.51增到0.65。訂正EASMI變幅也明顯增大,標(biāo)準(zhǔn)差從訂正前0.75增大到訂正后1.4,增幅為87%;大于1個(gè)標(biāo)準(zhǔn)差年數(shù)從6年升到12年,偏多了一倍(表1和圖9)。訂正EASMI與觀測(cè)距平符號(hào)一致率也從訂正前84%升到91%,上升了7%,訂正EASMI年際波動(dòng)與實(shí)況更為吻合(表1和圖9)。訂正EASMI有1980年、1982年、1988年和2000年等4年,由訂正前與觀測(cè)符號(hào)相反變?yōu)榉?hào)一致,但是訂正后EASMI在1996年和2008年,由訂正前與觀測(cè)符號(hào)一致轉(zhuǎn)為符號(hào)相反。訂正后EASMI變幅明顯增大的年份有14年,分別為1981、1982、1986、1988、1989、1992、1995、1996、1998、1999、2003、2008、2009和2010年,這14年中除了1981年外,其他13年都是ENSO年,表明利用PC訂正法對(duì)ENSO年訂正效果更好,也就是說海溫信號(hào)較強(qiáng)年訂正效果更顯著。

    表1 1979~2010年觀測(cè)的和訂正前、后的EASMI比較

    從訂正EASMI與觀測(cè)夏季降水量和850 hPa風(fēng)場(chǎng)的相關(guān)分布(圖8b)來,訂正后對(duì)應(yīng)西太平洋反氣旋位置比訂正前偏北,與實(shí)況較為接近,使得訂正后我國(guó)梅雨鋒區(qū)雨帶位置也偏北,與實(shí)況更為吻合,表明訂正EASMI對(duì)模擬梅雨鋒區(qū)雨帶偏南的缺點(diǎn)也有一定的改善。

    4 結(jié)論

    本文利用日本氣象廳海氣耦合模式MRI- CGCM的輸出結(jié)果和觀測(cè)資料進(jìn)行比較,從東亞夏季風(fēng)氣候態(tài)、主模態(tài)和年際變化等三個(gè)方面,來評(píng)估模式對(duì)東亞夏季風(fēng)的預(yù)測(cè)性能,在此基礎(chǔ)上,利用PC訂正法對(duì)東亞夏季風(fēng)指數(shù)進(jìn)行了訂正,得到以下結(jié)論:

    (1)MRI-CGCM模式對(duì)東亞夏季風(fēng)氣候態(tài)有較好的模擬能力,能模擬出東亞地區(qū)夏季降水從北到南遞增、西北太平洋反氣旋以及我國(guó)東部受該氣旋的西北側(cè)西南氣流控制的特征。但模擬的西北太平洋反氣旋偏弱、偏東,導(dǎo)致模擬的副熱帶地區(qū)降水量偏小。

    (2)模式能夠較好地模擬出東亞夏季風(fēng)降水的第一模態(tài)空間分布特征,模擬的EOF1與觀測(cè)之間的ACC為0.72,通過了0.01的顯著水平檢驗(yàn);模式也能較好地再現(xiàn)其對(duì)應(yīng)的年際變率,其時(shí)間系數(shù)PC1與觀測(cè)之間的相關(guān)系數(shù)為0.41,能模擬出觀測(cè)EOF1的2 a和5 a主導(dǎo)周期;但模擬的降水量年際變率幅度偏小和我國(guó)以東梅雨鋒區(qū)雨帶位置偏南,這與模擬的西北太平洋反氣旋位置偏南偏弱有關(guān)。MRI-CGCM較好地再現(xiàn)出EOF1對(duì)應(yīng)El Ni?o衰減位相,熱帶印度洋和西太平洋夏季SST偏暖,而熱帶中東太平洋夏季SST偏冷,與觀測(cè)相符。

    (3)模式對(duì)降水EOF2的模擬能力明顯差于EOF1,模擬的EOF2與觀測(cè)之間的ACC為0.36,相應(yīng)的模擬PC2與觀測(cè)之間的相關(guān)系數(shù)為0.35,模式能再現(xiàn)觀測(cè)EOF2的2 a主導(dǎo)周期,但未能模擬出EOF2的5 a主導(dǎo)周期。模式能再現(xiàn)出EOF2對(duì)應(yīng)El Ni?o發(fā)展位相,但模擬的西太平洋反氣旋明顯偏南,導(dǎo)致模擬的雨帶位置偏南,模擬的梅雨鋒區(qū)雨帶位置位于江南,與觀測(cè)場(chǎng)上江南降水量偏少相反,模擬的EOF2雨帶位置與觀測(cè)差異較大。

    (4)MRI-CGCM較好地模擬出我國(guó)夏季降水和氣溫的空間異常分布和年際變化。模擬與觀測(cè)夏季降水和氣溫多年平均ACC分別為0.74和0.68,均通過了0.01的顯著水平檢驗(yàn),其中降水平均ACC高于目前國(guó)家氣候中心預(yù)測(cè)業(yè)務(wù)上應(yīng)用的多模式解釋應(yīng)用集成預(yù)測(cè)。模式模擬我國(guó)東部、江淮流域和華南地區(qū)夏季降水多年平均PS評(píng)分分別為69、70和68分,略高于我國(guó)夏季降水業(yè)務(wù)多年平均(65分)。模擬的我國(guó)東部夏季氣溫與觀測(cè)多年平均符號(hào)一致率和PS評(píng)分分別為62%和74分。

    (5)利用觀測(cè)EASMI與模擬PC的關(guān)系建立多元線性回歸方程來訂正EASMI,訂正后EASMI在相關(guān)系數(shù)、符號(hào)一致率、模擬EASMI變幅偏小和我國(guó)梅雨鋒區(qū)雨帶偏南等方面均有一定的改善。訂正EASMI與實(shí)況的相關(guān)系數(shù)由0.51提高到0.65,符號(hào)一致率由84%升到91%;訂正后EASMI標(biāo)準(zhǔn)差從0.75增大到1.4,大于1個(gè)標(biāo)準(zhǔn)差年數(shù)由6年變?yōu)?2年;訂正后對(duì)應(yīng)西太平洋反氣旋位置和我國(guó)梅雨鋒區(qū)雨帶位置與實(shí)況更為吻合;PC訂正法對(duì)ENSO年的訂正效果更為明顯。訂正EASMI在夏季預(yù)測(cè)業(yè)務(wù)中應(yīng)該具有實(shí)際的應(yīng)用價(jià)值。

    致謝 感謝兩位匿名審稿專家對(duì)論文修改提出的寶貴意見。

    Adler R F, Huffman G J, Chang A, et al. 2003. The version-2 Global Precipitation Climatology Project (GPCP) monthly precipitation analysis (1979–present) [J]. J. Hydrometeor., 4 (6): 1147–1167, doi:10.1175/1525- 7541(2003)004<1147:TVGPCP>2.0.CO;2.

    陳麗娟, 袁媛, 楊明珠, 等. 2013. 海溫異常對(duì)東亞夏季風(fēng)影響機(jī)理的研究進(jìn)展 [J]. 應(yīng)用氣象學(xué)報(bào), 24 (5): 521–532. Chen Lijuan, Yuan Yuan, Yang Mingzhu, et al. 2013. A review of physical mechanisms of the global SSTA impact on EASM [J]. J. Appl. Meteor. Sci. (in Chinese), 24 (5): 521–532, doi:10.3969/j.issn.1001-7313.2013.05.002.

    Collins M, AchutaRao K, Ashok K, et al. 2013. Observational challenges in evaluating climate models [J]. Nature Climate Change, 3 (11): 940–941, doi:10.1038/nclimate2012.

    丁一匯. 2011. 季節(jié)氣候預(yù)測(cè)的進(jìn)展和前景 [J]. 氣象科技進(jìn)展, 1 (3): 14–27. Ding Yihui. 2011. Progress and prospects of seasonal climate prediction [J]. Adv. Meteor. Sci. Technol. (in Chinese), 1 (3): 14–27.

    董敏, 吳統(tǒng)文, 王在志, 等. 2013. BCC_CSM1.0模式對(duì)20世紀(jì)降水及其變率的模擬 [J]. 應(yīng)用氣象學(xué), 24 (1): 1–11. Dong Min, Wu Tongwen, Wang Zaizhi, et al. 2013. Simulation of the precipitation and its variation during the 20th century using the BCC climate model (BCC_CSM1.0) [J]. J. Appl. Meteor. Sci. (in Chinese), 24 (1): 1–11, doi:10.3969/j.issn. 1001-7313.2013.01.001.

    汽輪機(jī)整個(gè)啟動(dòng)過程就是一個(gè)緩慢均勻的加熱過程,各部溫差及膨脹不正常的情況下就容易發(fā)生碰磨產(chǎn)生振動(dòng),升速率及暖機(jī)點(diǎn)的選擇和控制就至關(guān)重要。本次汽封更換間隙調(diào)整后,除臨界轉(zhuǎn)速區(qū)以外,升速率及機(jī)組帶負(fù)荷速度均應(yīng)比正常啟動(dòng)緩慢,以便汽輪機(jī)加熱和汽封磨合,通過脹差、汽缸膨脹、各部溫度綜合判斷暖機(jī)效果及汽封磨合狀況,再?zèng)Q定升速或加負(fù)荷,切忌單純憑時(shí)間或某一參數(shù)決定。

    Gao Hui, Yang Song, Kumar A, et al. 2011. Variations of the East Asian Meiyu and simulation and prediction by the NCEP climate forecast system [J]. J. Climate, 24 (1): 94–108, doi:10.1175/2010JCLI3540.1.

    He Chao, Zhou Tianjun. 2014. The two interannual variability modes of the western North Pacific subtropical high simulated by 28 CMIP5-AMIP models [J]. Climate Dyn., 43 (9–10): 2455–2469, doi:10.1007/s00382- 014-2068-x.

    He Chao, Zhou Tianjun. 2015. Responses of the western North Pacific subtropical high to global warming under RCP4.5 and RCP8.5 scenarios projected by 33 CMIP5 models: The dominance of tropical Indian Ocean–tropical western Pacific SST gradient [J]. J. Climate, 28 (1): 365–380, doi:10.1175/JCLI-D-13-00494.1.

    何超, 周天軍, 鄒立維, 等. 2012. 夏季西北太平洋副熱帶高壓的兩種年際變率模態(tài) [J]. 中國(guó)科學(xué): 地球科學(xué), 42 (12): 1923–1936. He Chao, Zhou Tianjun, Zou Liwei, et al. 2013. Two interannual variability modes of the northwestern Pacific subtropical anticyclone in boreal summer [J]. Sci. China Earth Sci., 56 (7): 1254–1265.

    黃榮輝. 1990. 引起我國(guó)夏季旱澇的東亞大氣環(huán)流異常遙相關(guān)及其物理機(jī)制的研究 [J]. 大氣科學(xué), 14 (1): 108–117. Huang Ronghui. 1990. Studies on the teleconnections of the general circulation anomalies of East Asia causing the summer drought and flood in China and their physical mechanism [J]. Chinese Journal of Atmospheric Sciences (Scientia Atmos. Sinica) (in Chinese), 14 (1): 108–117, doi:10.3878/j. issn.1006-9895.1990.01.14.

    黃榮輝, 陳文. 2002. 關(guān)于亞洲季風(fēng)與ENSO循環(huán)相互作用研究最近的進(jìn)展 [J]. 氣候與環(huán)境研究, 7 (2): 146–159. Huang Ronghui, Chen Wen. 2002. Recent progresses in the research on the interaction between Asian monsoon and ENSO cycle [J]. Climatic and Environmental Research (in Chinese), 7 (2): 146–159, doi:10.3878/j.issn.1006-9585. 2002.02.03.

    黃榮輝, 蔡榕碩, 陳際龍, 等. 2006. 我國(guó)旱澇氣候?yàn)?zāi)害的年代際變化及其與東亞氣候系統(tǒng)變化的關(guān)系 [J]. 大氣科學(xué), 30 (5): 730–743. Huang Ronghui, Cai Rongshuo, Chen Jilong, et al. 2006. Interdecaldal variations of drought and flooding disasters in China and their association with the East Asian climate system [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 30 (5): 730–743, doi:10.3878/j.issn.1006-9895.2006.05.02.

    賈小龍, 陳麗娟, 高輝, 等. 2013. 我國(guó)短期氣候預(yù)測(cè)技術(shù)進(jìn)展 [J]. 應(yīng)用氣象學(xué)報(bào), 24 (6): 641–655. Jia Xiaolong, Chen Lijuan, Gao Hui, et al. 2013. Advances of the short-range climate prediction in China [J]. J. Appl. Meteor. Sci. (in Chinese), 24 (6): 641–655, doi:10.11898/1001- 7313.20130601.

    姜大膀, 王會(huì)軍, 郎咸梅. 2004. 全球變暖背景下東亞氣候變化的最新情景預(yù)測(cè) [J]. 地球物理學(xué)報(bào), 47 (4): 590–596. Jiang Dabang, Wang Huijun, Lang Xianmei. 2004. East Asian climate change trend under global warming background [J]. Chinese J. Geophys. (in Chinese), 47 (4): 590–596, doi:10.3321/j.issn:0001-5733.2004.04.007.

    Kistler R, Collins W, Saha A, et al. 1996. The NCEP–NCAR 50-year reanalysis: Monthly means CD-ROM and documentation [J]. Bull. Amer. Meteor. Soc., 82 (2): 247–268, doi:10.1175/1520-0477(2001)082<0247: TNNYRM>2.3.CO;2.

    劉長(zhǎng)征, 杜良敏, 柯宗建, 等. 2013. 國(guó)家氣候中心多模式解釋應(yīng)用集成預(yù)測(cè) [J]. 應(yīng)用氣象學(xué)報(bào), 24 (6): 677–685. Liu Changzheng, Du Liangmin, Ke Zongjian, et al. 2013. Multi-model downscaling ensemble prediction in National Climate Center [J]. J. Appl. Meteor. Sci. (in Chinese), 24 (6): 677–685, doi:10.3969/j.issn.1001-7313.2013.06.004.

    Rajendran K, Kitoh A, Yukimoto S. 2004. South and East Asian summer monsoon climate and variation in the MRI Coupled Model (MRI-CGCM2) [J]. J. Climate, 17 (4): 763–782, doi:10.1175/1520- 0442(2004)017<0763:SAEASM>2.0.CO;2.

    司東, 丁一匯, 柳艷菊. 2009. 全球海氣耦合模式(BCC_CM1.0)對(duì)江淮梅雨降水預(yù)報(bào)的檢驗(yàn) [J]. 氣象學(xué)報(bào), 67 (6): 947–960. Si Dong, Ding Yihui, Liu Yanju. 2009. Evaluation of Meiyu prediction in the Yangtze–Huaihe region by coupled ocean–atmosphere general circulation model (BCC_CM1.0) [J]. Acta Meteor. Sinica (in Chinese), 67 (6): 947–960, doi:10.11676/qxxb2009.092.

    Smith T M, Reynolds R W, Peterson T C, et al. 2008. Improvements to NOAA’s historical merged land–ocean surface temperature analysis (1880–2006) [J]. J. Climate, 21 (10): 2283–2296, doi:10.1175/ 2007JCLI2100.1.

    Song Fengfei, Zhou Tianjun. 2014a. Interannual variability of East Asian summer monsoon simulated by CMIP3 and CMIP5 AGCMs: Skill dependence on Indian Ocean–western Pacific anticyclone teleconnection [J]. J. Climate, 27 (4): 1679–1697, doi:10.1175/JCLI-D-13-00248.1.

    Song Fengfei, Zhou Tianjun. 2014b. The climatology and interannual variability of East Asian summer monsoon in CMIP5 coupled models: Does air–sea coupling improve the simulations? [J] J. Climate, 27 (23): 8761–8777, doi:10.1175/JCLI-D-14-00396.1.

    Song Fengfei, Zhou Tianjun, Qian Yun. 2014. Responses of East Asian summer monsoon to natural and anthropogenic forcings in the 17 latest CMIP5 models [J]. Geophys. Res. Lett., 41 (2), doi:10.1002/ 2013GL058705.

    孫穎, 丁一匯. 2008. IPCC AR4氣候模式對(duì)東亞夏季風(fēng)年代際變化的模擬性能評(píng)估 [J]. 氣象學(xué)報(bào), 66 (5): 765–780. Sun Ying, Ding Yihui. 2008. Validation of IPCC AR4 climate models in simulating interdecadal change of East Asian summer monsoon [J]. Acta Meteor. Sinica (in Chinese), 66 (5): 765–780, doi:10.11676/qxxb2008.070.

    Tao Shiyan, Chen Longxun. 1987. A review of recent research on the East Asian summer monsoon in China [M]//Monsoon Meteorology. Chang C P, Krishnamurti T N, Eds. Oxford: Oxford University Press, 62–92.

    Wang Bin, Fan Zhen. 1999. Choice of South Asian summer monsoon indices [J]. Bull. Amer. Meteor. Soc., 80 (4): 629–638, doi:10.1175/1520-0477(1999)080<0629:COSASM>2.0.CO;2.

    Wang Bin, Wu Zhiwei, Li Jianping, et al. 2008. How to measure the strength of the East Asian summer monsoon [J]. J. Climate, 21 (17): 4449–4463, doi:10.1175/2008JCLI2183.1.

    Wang Bin, Lee J Y, Xiang Baoqiang. 2015. Asian summer monsoon rainfall predictability: A predictable mode analysis [J]. Climate Dyn., 44 (1–2): 61–74, doi:10.1007/s00382-014-2218-1.

    Wu Bo, Zhou Tianjun, Li T. 2009. Seasonally evolving dominant interannual variability modes of East Asian climate [J]. J. Climate, 22 (11): 2992–3005.

    吳統(tǒng)文, 宋連春, 李偉平, 等. 2013. 北京氣候中心氣候系統(tǒng)模式研發(fā)進(jìn)展——在氣候變化研究中的應(yīng)用 [J]. 氣象學(xué)報(bào), 72 (1): 12–29. Wu Tongwen, Song Lianchun, Li Weiping, et al. 2013. An overview on progress in Beijing Climate Center climate system model—Its development and application to climate change studies [J]. Atca Meteor. Sinaca (in Chinese), 72 (1): 12–29, doi:10.11676/qxxb2013.084.

    Xie Pingping, Arkin P A. 1997. Global precipitation: A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs [J]. Bull. Amer. Meteor. Soc., 78 (11): 2539–2558, doi:10.1175/1520-0477(1997)078<2539:GPAYMA>2.0.CO;2.

    Xie S P, Hu K M, Hafner J, et al. 2009. Indian Ocean capacitor effect on Indo-Western Pacific climate during the summer following El Ni?o [J]. J. Climate, 22 (3): 730–747, doi:10.1175/2008JCLI2544.1.

    Yukimoto S, Noda A, Kitoh A, et al. 2001. The new meteorological research institute coupled GCM (MRI-CGCM2): Model climate and variability [J]. Papers Meteor. Geophys., 51 (2): 47–88, doi:10.2467/mripapers.51.47.

    張莉, 丁一匯, 孫穎. 2008. 全球海氣耦合模式對(duì)東亞季風(fēng)降水模擬的檢驗(yàn) [J]. 大氣科學(xué), 32 (2): 261–276. Zhang Li, Ding Yihui, Sun Ying. 2008. Evaluation of precipitation simulation in East Asian monsoon areas by coupled ocean–atmosphere general circulation models [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 32 (2): 261–276, doi:10.3878/j.issn.1006-9895.2008.02.06.

    張慶云, 陶詩言, 陳烈庭. 2003. 東亞夏季風(fēng)指數(shù)的年際變化與東亞大氣環(huán)流 [J]. 氣象學(xué)報(bào), 61 (4): 559–568. Zhang Qingyun, Tao Shiyan, Chen Lieting. 2003. The interannual variability of East Asian summer monsoon indices and its association with the pattern of general circulation over East Asia [J]. Acta Meteor. Sinica (in Chinese), 61 (4): 559–568, doi:10.3321/j.issn:0577-6619.2003.05.005.

    Zhou Tianjun, Yu Rucong, Li Hongmei, et al. 2008. Ocean forcing to changes in global monsoon precipitation over the recent half-century [J]. J. Climate, 21 (15): 3833–3852, doi:10.1175/2008JCLI2067.1.

    Zou Liwei, Zhou Tianjun. 2013. Can a regional ocean–atmosphere coupled model improve the simulation of the interannual variability of the western North Pacific summer monsoon? [J] J. Climate, 26 (7): 2353–2367, doi:10.1175/JCLI-D-11-00722.1.

    羅連升, 段春峰, 楊瑋, 等. 2016.MRI-CGCM模式對(duì)東亞夏季風(fēng)的模擬評(píng)估及訂正[J].大氣科學(xué), 40 (6): 1320–1332. Luo Liansheng, Duan Chunfeng, Yang Wei, et al.2016.Evaluation and correction of the East Asian summer monsoon simulated by MRI coupled ocean–atmosphere general circulation model[J]. Chinese Journal of Atmospheric Sciences (in Chinese), 40 (6): 1320–1332, doi:10.3878/j.issn.1006-9895.1602.15188.

    Evaluation and Correction of the East Asian Summer Monsoon Simulated by MRI Coupled Ocean–Atmosphere General Circulation Model

    LUO Liansheng, DUAN Chunfeng, YANG Wei, XU Min, CHENG Zhi, and DING Xiaojun

    ,230031

    Based on the monthly precipitation data of Climate Prediction Center Merged Analysis of Precipitation (CMAP), the NCEP reanalysis data, the Extended Reconstruction Sea Surface Temperature (ERSST) from National Oceanic and Atmospheric Administration (NOAA), and the output of the MRI-CGCM (Meteorological Research Institute Coupled Ocean–Atmosphere General Circulation Model) of Japan Meteorological Agency, the ability of MRI-CGCM to simulate precipitation in the East Asian Summer Monsoon (EASM) region is evaluated from the perspective of the climatology, the primary modes of Empirical Orthogonal Function (EOF) and the interannual variation of the EASM. The multiple linear regression equation is established by the relationship between the observed East Asian Summer Monsoon index (EASMI) and the simulated principal component (PC), and applied to correct the EASMI (PC correction method). Results show that the MRI-CGCM can reasonably reproduce the basic EASM rainfall and low-level wind fields. However, the simulated western North Pacific anticyclone is weaker and eastward-shifted than normal, which leads to the underestimation of rainfall over the subtropical region. MRI-CGCM can capture the first leading EOF mode (EOF1) of the EASM rainfall and the corresponding wind fields in lower levels over the EASM region and the decaying phase of El Ni?o. The EOF1 space Correlation Coefficient (ACC) between the simulation and observation is 0.72. The interannual variability of EOF1 is reasonably simulated, and the correlation coefficient between the first component (PC1) of MRI-CGCM simulation and observation is 0.41. The simulated EOF1 well reflects the observed characteristics of EOF1. However, the simulated Mei-yu rainbelt over eastern China shifts southward, which is closely related to the southward shift of the western North Pacific anticyclone. The model ability for the simulation of the second leading EOF mode (EOF2) decreases significantly compared to that for the EOF1. The EOF2 ACC between simulation and observation is 0.36. MRI-CGCM can well reproduce the EOF2 that corresponds to the developing phase of El Ni?o. However, the simulated western North Pacific anticyclone shifts southward abnormally, which leads to the southward shift of the rain belt. The simulated Mei-yu rain belt over China is located in the middle and lower reaches of the Yangtze River , which is contrary to the observation that little rain occurs over this region. The spatial distribution and interannual variability of summer precipitation and temperature anomalies in eastern China are well simulated. The mean ACC between the simulated and observed summer precipitation (temperature) is 0.74 (0.68). The mean predict score (PS) of simulated summer precipitation over eastern China, the Yangtze-Huaihe River valley and southern China are 69, 70, and 68, respectively, which are higher than the average PS (65) of the operational summer precipitation prediction. The mean PS of summer mean temperature in eastern China is 74. Improvements in the PC corrected EASMI are reflected in the correlation coefficient, the anomaly sign consistency rate, the weaker magnitude of the simulated EASMI, and the southward shift of the Mei-yu rain belt. The correlation coefficient between the corrected EASMI and observed EASMI increases from 0.51 to 0.65, the anomaly sign consistency rate changes from 84% to 91%, the standard deviation increases from 0.75 to 1.3, the number of years with greater than one standard deviation changes from 6 to 12, and the locations of the western North Pacific anticyclone and Mei-yu rain belt corresponding to the corrected EASMI are more consistent with observations.

    East Asian Summer Monsoon (EASM), EOF, Coupled ocean–atmosphere general circulation model (CGCM), Model correction

    1006-9895(2016)06-1320-13

    P467

    A

    10.3878/j.issn.1006-9895.1602.15188

    2015-04-27;網(wǎng)絡(luò)預(yù)出版日期 2016-02-02

    羅連升,女,1977年出生,高工,主要研究短期氣候預(yù)測(cè)。E-mail:luolsh_gx@163.com

    公益性行業(yè)(氣象)科研專項(xiàng)201406021,安徽省自然科學(xué)基金項(xiàng)目1308085QD69、1408085MD73

    Public Science and Technology Research Founds Projects of Meteorology (Grant 201406021), Nature Science Foundation of Anhui Province (Grants 1308085QD69,1408085MD73)

    猜你喜歡
    反氣旋雨帶年際
    王金成:秋風(fēng)
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    案例教學(xué)法在“氣旋、反氣旋與天氣”教學(xué)設(shè)計(jì)中的應(yīng)用
    考試周刊(2019年2期)2019-01-28 10:08:56
    《低壓(氣旋)、高壓(反氣旋)與天氣》教學(xué)設(shè)計(jì)
    考試周刊(2017年66期)2018-01-30 22:20:31
    關(guān)于反氣旋(高壓)的若干問題解析
    關(guān)于反氣旋(高壓)的若干問題解析
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測(cè)
    與北大西洋接壤的北極海冰和年際氣候變化
    各種各樣的雨
    最大熵算法在氣象雨量預(yù)測(cè)中的應(yīng)用分析
    免费观看的影片在线观看| 少妇人妻一区二区三区视频| 欧美日韩综合久久久久久 | 国产高清激情床上av| 日韩 欧美 亚洲 中文字幕| 小说图片视频综合网站| 三级男女做爰猛烈吃奶摸视频| 白带黄色成豆腐渣| 黄色成人免费大全| 国产乱人伦免费视频| 国产一区在线观看成人免费| 成人亚洲精品av一区二区| 在线视频色国产色| 亚洲成人中文字幕在线播放| 久久精品国产自在天天线| 国产欧美日韩一区二区三| 91久久精品电影网| 亚洲美女黄片视频| 美女高潮喷水抽搐中文字幕| 国产av不卡久久| 亚洲av熟女| 国产成人aa在线观看| 最新在线观看一区二区三区| 午夜亚洲福利在线播放| 精品久久久久久久末码| 最近在线观看免费完整版| 在线免费观看不下载黄p国产 | 亚洲精品成人久久久久久| 久久国产乱子伦精品免费另类| avwww免费| 91麻豆精品激情在线观看国产| 亚洲一区高清亚洲精品| 亚洲最大成人手机在线| 操出白浆在线播放| 操出白浆在线播放| 法律面前人人平等表现在哪些方面| x7x7x7水蜜桃| 国产精品 欧美亚洲| 国产v大片淫在线免费观看| 久久久久国内视频| 18禁裸乳无遮挡免费网站照片| 欧美成人a在线观看| 我要搜黄色片| 91在线精品国自产拍蜜月 | 久久久国产成人精品二区| 日韩有码中文字幕| 欧美色欧美亚洲另类二区| 一级毛片女人18水好多| 91字幕亚洲| 久久久久久久久久黄片| av女优亚洲男人天堂| 国产精品一区二区免费欧美| 一区二区三区高清视频在线| 九九热线精品视视频播放| 亚洲人成电影免费在线| 日韩欧美在线乱码| 欧美av亚洲av综合av国产av| 91久久精品国产一区二区成人 | 精品电影一区二区在线| 成人特级av手机在线观看| 99精品久久久久人妻精品| 精品一区二区三区av网在线观看| 亚洲一区二区三区色噜噜| 99热只有精品国产| 一级作爱视频免费观看| 欧美高清成人免费视频www| 欧美又色又爽又黄视频| 久久久国产精品麻豆| 久久精品国产清高在天天线| 久久久久精品国产欧美久久久| 欧美日韩福利视频一区二区| 国产麻豆成人av免费视频| 少妇人妻精品综合一区二区 | 在线视频色国产色| 美女高潮喷水抽搐中文字幕| 淫秽高清视频在线观看| 嫩草影院精品99| 日韩欧美免费精品| 久久6这里有精品| 欧美成狂野欧美在线观看| 特大巨黑吊av在线直播| 久久国产精品影院| 国产三级黄色录像| 女人被狂操c到高潮| 老司机深夜福利视频在线观看| 美女免费视频网站| 欧美zozozo另类| 色av中文字幕| 手机成人av网站| 欧美日韩中文字幕国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 91九色精品人成在线观看| 国产精品影院久久| 国产伦人伦偷精品视频| 高清在线国产一区| 成人av一区二区三区在线看| 欧美性感艳星| 久99久视频精品免费| 桃红色精品国产亚洲av| 色哟哟哟哟哟哟| 免费av观看视频| 白带黄色成豆腐渣| 欧美又色又爽又黄视频| www.色视频.com| 久久久久久久午夜电影| 在线视频色国产色| 国产一区二区在线观看日韩 | 成年女人永久免费观看视频| 18禁美女被吸乳视频| 国产精华一区二区三区| 国产成人aa在线观看| 亚洲人成伊人成综合网2020| 人妻久久中文字幕网| svipshipincom国产片| 看黄色毛片网站| 欧美又色又爽又黄视频| 精品久久久久久久末码| 夜夜躁狠狠躁天天躁| 91久久精品电影网| 国产一区二区三区在线臀色熟女| 亚洲成a人片在线一区二区| 国产伦在线观看视频一区| 国产一区二区在线观看日韩 | 伊人久久精品亚洲午夜| 亚洲av美国av| 日韩有码中文字幕| 老熟妇乱子伦视频在线观看| 亚洲熟妇熟女久久| 国产精品久久久久久久久免 | 国产高潮美女av| 国产一区二区亚洲精品在线观看| 99久久无色码亚洲精品果冻| 亚洲av成人av| 一区二区三区激情视频| 在线观看免费视频日本深夜| 日韩欧美国产一区二区入口| 日韩欧美一区二区三区在线观看| 亚洲精品粉嫩美女一区| 制服丝袜大香蕉在线| 中文资源天堂在线| 搡老熟女国产l中国老女人| 99久久九九国产精品国产免费| 日韩精品青青久久久久久| 在线观看一区二区三区| 国产高清视频在线观看网站| bbb黄色大片| 国产精品99久久久久久久久| 宅男免费午夜| x7x7x7水蜜桃| 欧美日韩精品网址| 琪琪午夜伦伦电影理论片6080| 国产97色在线日韩免费| 午夜两性在线视频| 床上黄色一级片| 免费在线观看成人毛片| 亚洲天堂国产精品一区在线| 亚洲成人精品中文字幕电影| 精品欧美国产一区二区三| 亚洲中文日韩欧美视频| 欧美午夜高清在线| 老司机午夜十八禁免费视频| 国产老妇女一区| 久久香蕉国产精品| 蜜桃久久精品国产亚洲av| 日韩中文字幕欧美一区二区| 久久久久久国产a免费观看| 女人被狂操c到高潮| 国产乱人伦免费视频| 十八禁网站免费在线| 亚洲人与动物交配视频| 波野结衣二区三区在线 | 国产成人福利小说| 亚洲人成网站在线播| 全区人妻精品视频| 免费观看人在逋| 欧洲精品卡2卡3卡4卡5卡区| 国产三级中文精品| 最近最新免费中文字幕在线| 免费av观看视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品嫩草影院av在线观看 | 两人在一起打扑克的视频| 少妇高潮的动态图| 琪琪午夜伦伦电影理论片6080| 国产精品久久久久久久久免 | 日本五十路高清| 69av精品久久久久久| 怎么达到女性高潮| 男人舔女人下体高潮全视频| 在线观看免费午夜福利视频| 99视频精品全部免费 在线| 亚洲无线在线观看| 国产精品久久电影中文字幕| 久久精品人妻少妇| 亚洲精品色激情综合| 国产精品久久久久久人妻精品电影| 欧美色视频一区免费| 国产在视频线在精品| 国产精品1区2区在线观看.| www.999成人在线观看| e午夜精品久久久久久久| 亚洲专区国产一区二区| 观看美女的网站| 欧美大码av| 国产av在哪里看| 久久久久国内视频| 日本熟妇午夜| 中文字幕人妻丝袜一区二区| 午夜福利成人在线免费观看| 色哟哟哟哟哟哟| 久久久成人免费电影| 俄罗斯特黄特色一大片| 好男人在线观看高清免费视频| 在线播放国产精品三级| 老司机深夜福利视频在线观看| 91av网一区二区| 国产av不卡久久| 性色avwww在线观看| 久久人妻av系列| 小说图片视频综合网站| 每晚都被弄得嗷嗷叫到高潮| 国产精品,欧美在线| 99在线视频只有这里精品首页| 欧美中文日本在线观看视频| 国产伦在线观看视频一区| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 美女高潮的动态| 欧美最黄视频在线播放免费| 女人高潮潮喷娇喘18禁视频| av视频在线观看入口| 国产久久久一区二区三区| 亚洲精品在线美女| 国产高清视频在线观看网站| 1024手机看黄色片| 午夜日韩欧美国产| 亚洲精品一区av在线观看| 日日干狠狠操夜夜爽| 欧美黄色淫秽网站| 啦啦啦韩国在线观看视频| 91字幕亚洲| a级毛片a级免费在线| 国产成人欧美在线观看| 欧美不卡视频在线免费观看| 最新中文字幕久久久久| 欧美性感艳星| 亚洲男人的天堂狠狠| 99热精品在线国产| 无遮挡黄片免费观看| 欧美一区二区亚洲| 麻豆国产av国片精品| 18+在线观看网站| 久久久久亚洲av毛片大全| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| 欧美不卡视频在线免费观看| 国产精品一区二区三区四区久久| 啪啪无遮挡十八禁网站| 国产精品 欧美亚洲| 亚洲中文字幕一区二区三区有码在线看| 成年女人永久免费观看视频| 69av精品久久久久久| 一本精品99久久精品77| 国产高清视频在线播放一区| 欧美黑人欧美精品刺激| 1024手机看黄色片| 久久久久久久久中文| 国产午夜精品久久久久久一区二区三区 | 精品国产亚洲在线| 色综合婷婷激情| 校园春色视频在线观看| 九九在线视频观看精品| 亚洲性夜色夜夜综合| 欧美一级毛片孕妇| 亚洲人成网站在线播放欧美日韩| 国产国拍精品亚洲av在线观看 | 好看av亚洲va欧美ⅴa在| 色噜噜av男人的天堂激情| 亚洲av免费高清在线观看| 午夜日韩欧美国产| 黄色女人牲交| 中亚洲国语对白在线视频| 国产色婷婷99| 12—13女人毛片做爰片一| 有码 亚洲区| 无人区码免费观看不卡| 波多野结衣巨乳人妻| 亚洲精品在线观看二区| 欧美zozozo另类| 国产精品电影一区二区三区| 欧美一区二区精品小视频在线| 天美传媒精品一区二区| 亚洲国产精品久久男人天堂| 一进一出抽搐gif免费好疼| 精品免费久久久久久久清纯| 日本黄色片子视频| 国产免费男女视频| 黄色女人牲交| 亚洲国产色片| 亚洲va日本ⅴa欧美va伊人久久| 国产免费av片在线观看野外av| 性色avwww在线观看| 免费在线观看日本一区| 成人av一区二区三区在线看| 亚洲国产欧美人成| 亚洲av不卡在线观看| 两个人的视频大全免费| 乱人视频在线观看| 国产视频内射| 一个人看视频在线观看www免费 | 欧美中文日本在线观看视频| 91在线观看av| 国产高清视频在线观看网站| 法律面前人人平等表现在哪些方面| 亚洲国产高清在线一区二区三| 变态另类丝袜制服| 一进一出抽搐动态| 91九色精品人成在线观看| 国产精品久久久久久久久免 | 午夜免费激情av| 日本一本二区三区精品| 99国产极品粉嫩在线观看| 亚洲国产精品成人综合色| 国产成人av激情在线播放| 亚洲国产精品999在线| 一边摸一边抽搐一进一小说| 国产激情偷乱视频一区二区| 国产精品香港三级国产av潘金莲| 中文字幕av成人在线电影| 成人av一区二区三区在线看| 亚洲精品美女久久久久99蜜臀| aaaaa片日本免费| 手机成人av网站| 亚洲精品美女久久久久99蜜臀| 色视频www国产| av福利片在线观看| 在线免费观看的www视频| 性欧美人与动物交配| 午夜两性在线视频| 乱人视频在线观看| 亚洲精品乱码久久久v下载方式 | 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| 男女之事视频高清在线观看| 啦啦啦观看免费观看视频高清| 91字幕亚洲| 一个人观看的视频www高清免费观看| 99热这里只有是精品50| 悠悠久久av| 国产精品 国内视频| av专区在线播放| www日本黄色视频网| 色精品久久人妻99蜜桃| xxx96com| 成人国产一区最新在线观看| 国产一区二区三区在线臀色熟女| 国产中年淑女户外野战色| 色老头精品视频在线观看| 丰满的人妻完整版| 俺也久久电影网| 久久精品国产综合久久久| 国产免费av片在线观看野外av| 国内精品久久久久久久电影| 99视频精品全部免费 在线| 19禁男女啪啪无遮挡网站| 欧美极品一区二区三区四区| 亚洲专区国产一区二区| 国产成年人精品一区二区| 黄色片一级片一级黄色片| 国产精品国产高清国产av| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| 国产亚洲精品久久久久久毛片| 香蕉丝袜av| av专区在线播放| 日本免费一区二区三区高清不卡| 性色avwww在线观看| 天堂影院成人在线观看| 亚洲专区中文字幕在线| 91久久精品国产一区二区成人 | 亚洲,欧美精品.| 午夜久久久久精精品| 高清日韩中文字幕在线| 精品人妻一区二区三区麻豆 | 色av中文字幕| 九九久久精品国产亚洲av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 五月伊人婷婷丁香| 亚洲av电影不卡..在线观看| 少妇丰满av| 午夜福利成人在线免费观看| 午夜a级毛片| 久久伊人香网站| 精品国产美女av久久久久小说| 天堂影院成人在线观看| 免费高清视频大片| 在线国产一区二区在线| 性色av乱码一区二区三区2| 香蕉久久夜色| 天美传媒精品一区二区| 亚洲av二区三区四区| 69人妻影院| 午夜免费男女啪啪视频观看 | 在线观看一区二区三区| 国产真实乱freesex| 99久久九九国产精品国产免费| 中文字幕人妻丝袜一区二区| 久久久久性生活片| 亚洲专区国产一区二区| 一个人免费在线观看的高清视频| 国产成人av激情在线播放| 搡女人真爽免费视频火全软件 | 亚洲五月天丁香| 国产精品野战在线观看| 草草在线视频免费看| 日本撒尿小便嘘嘘汇集6| 夜夜看夜夜爽夜夜摸| a级毛片a级免费在线| 亚洲在线观看片| 亚洲精品色激情综合| 国产精品 欧美亚洲| 12—13女人毛片做爰片一| 桃红色精品国产亚洲av| 国语自产精品视频在线第100页| 天天添夜夜摸| 国产成人av激情在线播放| 18禁国产床啪视频网站| 两人在一起打扑克的视频| 久久久久久九九精品二区国产| 免费看十八禁软件| 国产熟女xx| 久久久成人免费电影| 欧美中文日本在线观看视频| 国产乱人视频| 亚洲精华国产精华精| 婷婷亚洲欧美| 精品免费久久久久久久清纯| 亚洲中文日韩欧美视频| 啦啦啦免费观看视频1| 久久久久亚洲av毛片大全| 99精品久久久久人妻精品| 少妇熟女aⅴ在线视频| 国产成人福利小说| 99热精品在线国产| 欧美黑人欧美精品刺激| 成人亚洲精品av一区二区| 中文字幕av成人在线电影| 日韩av在线大香蕉| 欧美绝顶高潮抽搐喷水| 噜噜噜噜噜久久久久久91| a级一级毛片免费在线观看| 欧美黑人欧美精品刺激| 亚洲国产中文字幕在线视频| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 亚洲 欧美 日韩 在线 免费| 亚洲av第一区精品v没综合| 看片在线看免费视频| 亚洲av免费在线观看| 制服人妻中文乱码| 亚洲片人在线观看| 在线看三级毛片| 婷婷精品国产亚洲av在线| 51国产日韩欧美| 国产午夜福利久久久久久| 国产精品久久久久久人妻精品电影| 757午夜福利合集在线观看| 丰满乱子伦码专区| 观看免费一级毛片| 男女做爰动态图高潮gif福利片| 欧美黄色淫秽网站| 成人国产一区最新在线观看| 1024手机看黄色片| 亚洲电影在线观看av| 夜夜爽天天搞| 亚洲av五月六月丁香网| 色播亚洲综合网| 欧美日韩精品网址| 国产黄色小视频在线观看| 国内精品久久久久精免费| 国产精品香港三级国产av潘金莲| 亚洲精品456在线播放app | www国产在线视频色| 亚洲美女黄片视频| 国产亚洲精品久久久com| 欧美乱色亚洲激情| 久久久久久久午夜电影| av黄色大香蕉| 久久精品国产自在天天线| 日本撒尿小便嘘嘘汇集6| ponron亚洲| 男人和女人高潮做爰伦理| av天堂中文字幕网| 亚洲国产精品成人综合色| 免费av不卡在线播放| 两个人视频免费观看高清| 国产亚洲精品综合一区在线观看| 九色成人免费人妻av| 国产色爽女视频免费观看| 久久这里只有精品中国| 亚洲国产中文字幕在线视频| 少妇丰满av| 欧美乱妇无乱码| 亚洲一区二区三区色噜噜| 欧美色视频一区免费| 国产高潮美女av| 久久久国产成人免费| 国产一区二区三区视频了| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 久久性视频一级片| 一级黄色大片毛片| 一边摸一边抽搐一进一小说| 99视频精品全部免费 在线| 给我免费播放毛片高清在线观看| 又紧又爽又黄一区二区| 久久精品91无色码中文字幕| 欧美日韩瑟瑟在线播放| 欧美一区二区亚洲| 久久久久免费精品人妻一区二区| 欧美一区二区亚洲| 人人妻人人看人人澡| 琪琪午夜伦伦电影理论片6080| 一区二区三区国产精品乱码| 香蕉丝袜av| 国产精品98久久久久久宅男小说| 毛片女人毛片| 亚洲精品亚洲一区二区| 国产日本99.免费观看| 亚洲国产精品sss在线观看| 国产精品永久免费网站| 亚洲熟妇中文字幕五十中出| 人人妻人人看人人澡| 国产高清有码在线观看视频| 小说图片视频综合网站| av女优亚洲男人天堂| 亚洲美女黄片视频| 久久伊人香网站| 老汉色av国产亚洲站长工具| 韩国av一区二区三区四区| 国产高清有码在线观看视频| 九九在线视频观看精品| 在线观看免费午夜福利视频| 日韩欧美精品v在线| 亚洲国产中文字幕在线视频| 欧美黑人欧美精品刺激| 午夜日韩欧美国产| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久久电影 | 在线观看av片永久免费下载| 少妇的丰满在线观看| 高清日韩中文字幕在线| 99久久精品一区二区三区| 欧美一区二区国产精品久久精品| 99热这里只有是精品50| 床上黄色一级片| 中文字幕高清在线视频| 18禁在线播放成人免费| 国产午夜精品论理片| 亚洲精华国产精华精| 又紧又爽又黄一区二区| 日本一二三区视频观看| 国产精品综合久久久久久久免费| 亚洲成a人片在线一区二区| 国内揄拍国产精品人妻在线| 亚洲av电影不卡..在线观看| 久久婷婷人人爽人人干人人爱| 国产综合懂色| 国产成人福利小说| 在线播放无遮挡| 日韩av在线大香蕉| 精品一区二区三区av网在线观看| 成人高潮视频无遮挡免费网站| 男女视频在线观看网站免费| 麻豆久久精品国产亚洲av| 亚洲最大成人手机在线| 久久久久久人人人人人| 国产97色在线日韩免费| 国产伦精品一区二区三区视频9 | 婷婷精品国产亚洲av| 色精品久久人妻99蜜桃| 国产精品亚洲一级av第二区| 日韩有码中文字幕| 亚洲国产精品sss在线观看| 欧美一级毛片孕妇| 国产蜜桃级精品一区二区三区| 亚洲在线观看片| 中文字幕熟女人妻在线| 18+在线观看网站| 少妇裸体淫交视频免费看高清| 国产毛片a区久久久久| 一边摸一边抽搐一进一小说| 我的老师免费观看完整版| 国产日本99.免费观看| 午夜a级毛片| 两个人看的免费小视频| 欧美性猛交黑人性爽| 一个人看视频在线观看www免费 | 国产毛片a区久久久久| 日韩亚洲欧美综合| 男女午夜视频在线观看| 精品99又大又爽又粗少妇毛片 | 国产成+人综合+亚洲专区| 久久久久久久精品吃奶| 亚洲一区二区三区色噜噜| 一级黄片播放器| 亚洲自拍偷在线| 法律面前人人平等表现在哪些方面| 久久99热这里只有精品18| 欧美大码av| 中文字幕熟女人妻在线| 日韩精品中文字幕看吧| a在线观看视频网站| 91久久精品电影网|