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

    基于前春對流層溫度和北大西洋濤動的中國夏季降水統(tǒng)計預(yù)測模型

    2011-12-15 07:06:00馮蕾魏鳳英朱艷峰
    大氣科學(xué) 2011年5期
    關(guān)鍵詞:對流層特征向量降水

    馮蕾 魏鳳英 朱艷峰

    1 中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實驗室,北京 100029

    2 中國科學(xué)院研究生院,北京 100049

    3 中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實驗室,北京 100081

    4 國家氣候中心,北京 100081

    基于前春對流層溫度和北大西洋濤動的中國夏季降水統(tǒng)計預(yù)測模型

    馮蕾1,2魏鳳英3朱艷峰4

    1 中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實驗室,北京 100029

    2 中國科學(xué)院研究生院,北京 100049

    3 中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實驗室,北京 100081

    4 國家氣候中心,北京 100081

    本文在分析中國夏季降水典型模態(tài)與前期春季對流層中上層溫度主要分布類型、北大西洋濤動 (簡稱NAO)之間關(guān)系的基礎(chǔ)上,提出基于前春對流層溫度和NAO的中國夏季降水統(tǒng)計預(yù)測模型,并對利用該模型預(yù)報的2004~2009年中國160站夏季降水進(jìn)行檢驗。結(jié)果表明:中國大范圍夏季降水多寡主要與5月NAO變化有關(guān),“南多 (少)北少 (多)型”降水異常分布主要與5月NAO異常以及前春中國東北、華北地區(qū)上空對流層溫度的冷異常有關(guān)。NAO對中國夏季降水的影響既包括年際尺度的,也包括年代際尺度的,而對流層溫度異常的影響主要體現(xiàn)在年代際尺度上。這兩個因子通過影響夏季大氣環(huán)流和季風(fēng)水汽輸送,從而影響中國夏季降水的趨勢分布?;谶@樣的聯(lián)系,本工作以中國夏季降水量原始場EOF分解前3個模態(tài)對應(yīng)的時間系數(shù)為預(yù)測對象,以前期春季平均的對流層中上層溫度距平REOF分解前4個模態(tài)對應(yīng)的時間系數(shù)以及逐月NAO指數(shù)序列作為預(yù)報因子,并結(jié)合降水序列本身、對流層溫度、NAO指數(shù)的多時間尺度分量建立了我國夏季降水的統(tǒng)計預(yù)測模型,取得了較好的預(yù)報效果。2004~2009年預(yù)報降水與觀測降水之間的6年平均距平相關(guān)系數(shù)達(dá)到0.335,表明該預(yù)報方案對我國夏季降水具有較高的預(yù)報技巧。

    對流層溫度 北大西洋濤動 中國夏季降水 統(tǒng)計預(yù)測模型

    1 引言

    中國位于東亞季風(fēng)區(qū),影響夏季降水的環(huán)流系統(tǒng)復(fù)雜,降水時空變率較大,夏季頻繁發(fā)生的旱澇災(zāi)害,給國民經(jīng)濟(jì)和人民生命安全帶來巨大損失。準(zhǔn)確預(yù)測夏季降水異常是業(yè)務(wù)部門和大氣科學(xué)界共同面臨的重要課題和難題 (丁一匯,2004;王斌,2008;Zhou et al.,2009a)。

    ENSO被認(rèn)為是季風(fēng)區(qū)乃至全球氣候可預(yù)報性的主要來源 (Wang et al.,2008;Wu et al.,2009,2010;Zhou et al.,2009b)。但對于東亞季風(fēng)區(qū),ENSO對降水變率的影響相對于其它地區(qū)要小得多(Wang et al.,2008)。因此,準(zhǔn)確預(yù)測中國夏季降水,還需要充分考慮ENSO以外的物理因子。目前,研究者們在做我國夏季降水預(yù)測時,考慮的預(yù)測因子除了ENSO外,還包括青藏高原積雪以及一些大氣環(huán)流變化特征量等 (趙振國,1996,1999;金祖輝和陶詩言,1999;Zhang et al.,1999;吳統(tǒng)文和錢正安,2000;黃嘉佑等,2004;彭京備等,2005;魏鳳英,2006;張慶云等,2007;武炳義等,2008;楊蓮梅和張慶云,2008;魏鳳英和黃嘉佑,2010;韓雪和魏鳳英,2010)。

    對流層溫度是一個綜合性的物理因子,其變化是各種動力因素和熱力因素共同作用的結(jié)果。各種環(huán)流因子、大氣活動中心對東亞氣候的影響,都繞不過對流層溫度變化這一環(huán)節(jié) (Zhou and Zou,2010)。近年來,圍繞對流層溫度變化對東亞氣候的影響,產(chǎn)生了一系列重要的研究成果:20世紀(jì)70年代末以來中國東部夏季降水的“南澇北旱”特征與東亞對流層中上層變冷存在直接的聯(lián)系 (Yu et al.,2004;Yu and Zhou,2007;宇如聰?shù)?2008),其物理概念模型可以表述為:對流層中上層出現(xiàn)的年代際尺度變冷,一方面通過其上層的氣旋式環(huán)流異常,使得東亞急流軸以南的西風(fēng)增強(qiáng),另一方面通過其下層的反氣旋式環(huán)流異常,導(dǎo)致東亞夏季風(fēng)減弱。西風(fēng)急流偏南和夏季風(fēng)減弱,最終導(dǎo)致中國東部“南澇北旱”型降水異常 (Xin et al.,2006)。中國夏季降水除了在東部有“南澇北旱”變化外,在中西部地區(qū) (90°E~110°E)也有“南澇(旱)北旱 (澇)”的變化,并且這種降水異常型與500 hPa等壓面上高原地區(qū)平均溫度和其東部平原平均溫度之差的異常變化有關(guān) (朱艷峰等,2010)。研究還表明,東亞對流層溫度的年代際變化是北半球?qū)α鲗訙囟饶甏H變化的組成部分 (Zhou and Zhang,2009)。早春季節(jié)歐亞大陸副熱帶地區(qū)的表層氣溫和對流層溫度的變冷趨勢,與冬季北大西洋濤動 (簡稱NAO)正位相增強(qiáng)所對應(yīng)的北半球副熱帶冷信號的準(zhǔn)正壓東傳有關(guān) (Yu and Zhou,2004)。東亞對流層上層變冷和青藏高原下游地面氣溫變冷開始于前春的3月,從3月到5月,冷中心逐漸南移并加強(qiáng),并且這種冷位相可以維持到9月份 (Yu and Zhou,2007)。前冬NAO的年代際增強(qiáng),使得東亞西風(fēng)急流增強(qiáng),通過影響高原積雪深度以及觸發(fā)高原東側(cè)獨(dú)特的云—輻射反饋過程,使得東亞對流層中上層溫度變冷 (Li et al.,2005,2008;辛?xí)愿璧?2009;Xin et al.,2010)。上述關(guān)系在年際尺度上同樣存在。盡管東亞對流層中上層溫度與東亞夏季降水的聯(lián)系受到前冬NAO的影響,但中國夏季降水對前冬NAO的響應(yīng)很弱,而與前春的NAO變化具有顯著的相關(guān) (王永波和施能,2001;龔道溢等,2002)。

    以上研究表明,前春季節(jié)的對流層溫度、NAO變化與中國夏季降水異常具有重要的影響,但目前將前春對流層溫度和NAO應(yīng)用于中國夏季降水預(yù)測的工作開展得很少。本工作的目的是定量分析前春對流層溫度和NAO變化對中國夏季降水的影響,并建立相應(yīng)的預(yù)測模型,對我國夏季 (6~8月)降水進(jìn)行預(yù)測試驗。

    本文第2節(jié)扼要介紹使用的資料和方法;第3節(jié)分析我國夏季降水以及前春對流層中上層溫度變化的主要分布類型,討論降水典型模態(tài)與對流層中上層溫度典型模態(tài)及NAO之間的聯(lián)系;第4節(jié)提出以前春對流層中上層溫度和NAO為主要預(yù)報因子的預(yù)測模型,并檢驗該模型對我國夏季降水趨勢分布的預(yù)測效果;第5節(jié)對全文進(jìn)行總結(jié)。

    2 資料和方法

    本文使用的降水資料為國家氣候中心氣候預(yù)測室提供的1951~2009年中國160站6~8月的降水資料,該資料被廣泛地應(yīng)用于東亞降水變率研究中(Zhou and Yu,2005)。所用的環(huán)流場資料為NCEP/NCAR逐月再分析資料 (Kalnay et al.,1996)。涉及到的物理量有氣溫、風(fēng)場、位勢高度場,分辨率為2.5°×2.5°。使用的 NAO指數(shù)來自http:∥www.cpc.noaa.gov/products/precip/CWlink/pna/nao.shtm l[2010-04-10]。

    本文使用的分析方法包括經(jīng)驗正交函數(shù)(EOF)、旋轉(zhuǎn)經(jīng)驗正交函數(shù) (REOF)、均值生成函數(shù)及回歸分析等方法 (魏鳳英,2007)。

    3 前春對流層中上層溫度及 NAO對中國夏季降水的影響

    3.1 中國夏季降水的主導(dǎo)模態(tài)

    在預(yù)測中,為了回答某一年中國地區(qū)總體降水趨勢和雨區(qū)分布問題,本文對中國160站1951~2009年夏季降水量原始場進(jìn)行 EOF分解。圖1給出前3個特征向量及其對應(yīng)的時間系數(shù)。由于是使用原始場進(jìn)行 EOF展開,那么,第1特征向量(EOF1)就是中國大范圍夏季降水的多年平均分布狀況,即自東南地區(qū)向西北地區(qū)遞減的分布型式(圖1a)(施能,2009)。第一特征向量對應(yīng)的時間系數(shù) (PC1)就是中國夏季大范圍降水總體多寡的逐年變化,PC1越大,意味著這一年全國大范圍降水是偏多的 (圖1b)。該模態(tài)能夠解釋總方差的90.90%。降水的第2特征向量 (EOF2)表現(xiàn)為華南—江淮—華北地區(qū)“-+-”的降水異常分布(圖1c),即以江淮地區(qū)為中心的三極型降水分布,這是我國最常出現(xiàn)的降水異常型,以下簡稱為“南北少 (多)中間多 (少)型”。從時間系數(shù) (PC2)來看 (圖1d),該類型降水年際振蕩比較顯著。第3特征向量 (EOF3)表現(xiàn)為長江以南地區(qū)降水與其以北地區(qū)相反的趨勢分布 (圖1e),這也是我國比較常見的降水分布型,以下簡稱為“南多 (少)北少 (多)型”。該模態(tài)對應(yīng)的時間系數(shù) (PC3)表明(圖1f),自1970年代末至2000年左右,“南澇北旱”型降水有明顯的增強(qiáng)趨勢,這與前人的結(jié)果(趙平和周秀驥,2006;Yu and Zhou,2007;Zhou et al.,2009a)一致。上述三個特征向量均通過了No rth準(zhǔn)則顯著性檢驗 (No rth et al.,1982),總計解釋了總方差的93.32%,決定著中國夏季降水的趨勢分布。

    3.2 前春對流層中上層溫度變化的主要模態(tài)

    為揭示春季 (3~5月)對流層溫度的變化模態(tài) ,對 (10°N~50°N,60°E~130°E) 區(qū)域內(nèi)的對流層平均溫度 (定義為200~500 hPa的厚度平均)距平做旋轉(zhuǎn)經(jīng)驗正交函數(shù) (REOF)分解。第1特征向量 (REOF1)解釋了總方差的22.7%(圖2a),其中心位于20°N以南地區(qū),以下稱為“亞洲低緯型”。第一模態(tài)對應(yīng)的時間系數(shù) (PC1′)表明 (圖2b),該模態(tài)具有顯著的年代際振蕩特征。第2特征向量 (REOF2)解釋了總方差的16.6%(圖2c),其中心主要位于我國東北、華北地區(qū)上空,以下稱為“亞洲中緯型”。對應(yīng)的時間系數(shù) (圖2d,PC2′)表明,1960年代以前,該模態(tài)無明顯的變化,1960年代中期到1990年代末期,該模態(tài)具有顯著的增強(qiáng)趨勢,2000年后又呈現(xiàn)出減弱趨勢。從時間演變上看,“亞洲中緯型”對流層中上層溫度的變化與我國“南多 (少)北少 (多)型”降水變化比較一致,對應(yīng)的溫度、降水模態(tài)的PC序列間的相關(guān)系數(shù)為0.3。第3特征向量 (REOF3)解釋了總方差的13.9%(圖2e),中心主要位于伊朗地區(qū)上空,以下簡稱“南亞型”。對應(yīng)的時間系數(shù) (圖 2f,PC3′)表明,該模態(tài)具有明顯的年代際和年際變化特征。第4特征向量 (REOF4)解釋了總方差的11.5%(圖2g),中心位于我國青藏高原地區(qū)上空,以下簡稱“高原型”。North準(zhǔn)則檢驗表明,上述4個特征向量是獨(dú)立的,總計解釋了總方差的64.7%。

    圖1 1951~2009年我國夏季降水量場EOF分解前3個特征向量 (a、c、e)及時間序列 (b、d、f)Fig.1 (a,c,e)The first three EOFmodesof summer precipitation over China for the period of 1951-2009 and(b,d,f)the corresponding time series

    3.3 中國夏季降水異常和前春對流層中上層溫度、

    NAO異常的聯(lián)系

    3.3.1 相關(guān)分析

    為揭示中國夏季降水異常與前春對流層中上層溫度、NAO活動異常之間的聯(lián)系,首先計算我國夏季降水量前3個特征向量時間系數(shù)序列 (PC1、PC2、PC3)與前春對流層溫度前4個特征向量時間系數(shù)序列 (PC1′、PC2′、PC3′、PC4′) 以及3、4、5月逐月的NAO指數(shù)之間的相關(guān) (如表1所示)。同時,為了考察降水與對流層中上層溫度、NAO在年代際和年際時間尺度上的聯(lián)系,本文使用權(quán)重系數(shù)長度為9,截斷頻率為0.1(即10年周期)的Lanczos濾波器對各序列進(jìn)行低通和高通濾波,分別得到其年代際變化 (具有10年以上時間尺度的周期)和年際變化 (具有10年以下時間尺度的周期)。濾波之后各序列的相關(guān)如表2和表3所示。

    圖2 1951~2009年前春對流層中上層溫度REOF分解前4個特征向量 (a、c、e、g)及時間序列 (b、d、f、h)Fig.2 (a,c,e,g)The first four rotated empiricalorthogonal function(REOF)modesof spring atmosphere temperature averaged over 500 hPa-200 hPa for the period of 1951-2009 and(b,d,f,h)the corresponding time series

    表1 1951~2009年我國夏季降水前3個特征向量時間系數(shù)序列與前春對流層中上層溫度前4個特征向量時間系數(shù)序列、逐月NAO指數(shù)的相關(guān)Table 1 The correlation coefficients between the first three summer precipitation time series(PC1,PC2,PC3)and the first four spring atmosphere temperature time series(PC1′,PC2′,PC3′,PC4′),NAO index during the period of 1951-2009

    考察表1到表3,首先,從總體上看,中國夏季降水的第1模態(tài),即大范圍降水趨勢的變化主要與前期5月的NAO指數(shù)變化有關(guān)。降水第3模態(tài),即“南多 (少)北少 (多)型”降水異常分布與前春對流層中上層溫度變化的第2模態(tài)以及5月NAO指數(shù)變化密切相關(guān),而降水的第2模態(tài),即“南北少 (多)中間多 (少)型”降水與前春對流層中上層溫度及NAO指數(shù)的變化均無明顯的關(guān)聯(lián);第二,前春5月NAO指數(shù)與中國夏季大范圍降水及“南多 (少)北少 (多)型”降水異常分布的聯(lián)系既有年際尺度上,也包含年代際尺度上的。前春 3月NAO變化與我國“南多 (少)北少 (多)型”夏季降水的相關(guān)在年代際尺度上很顯著,但在年際尺度上很弱。而前春3月NAO變化與中國大范圍夏季降水、以及前春4月NAO變化與我國“南多 (少)北少 (多)型”夏季降水的相關(guān)主要是年際尺度上的。第三,對流層中上層溫度第1模態(tài)與我國大范圍夏季降水變化具有較強(qiáng)的年際尺度上的相關(guān)。除此之外,對流層溫度變化各主要模態(tài)與中國夏季降水的聯(lián)系主要是年代際尺度上的。表1中原始序列包含了不同時間尺度的變化,使得降水與對流層溫度之間的相關(guān)未通過5%的顯著水平的顯著性檢驗。因此,有必要將各時間序列不同時間尺度的周期,即將可預(yù)報的分量提取出來,也作為預(yù)報因子,可以排除噪音,達(dá)到更好的預(yù)報效果。

    需要指出的是,不同季節(jié)的NAO活動對東亞氣候的影響不同。前冬NAO活動異常與前春東亞對流層溫度之間有著較好的對應(yīng)關(guān)系,即冬季NAO的正位相年,較強(qiáng)的西風(fēng)急流使得高原下游的云量增加,通過云—溫度反饋機(jī)制維持著東亞對流層溫度的冷異常 (Yu and Zhou,2004;Li et al.,2005)。計算表明,前冬NAO指數(shù)序列與前春對流層中上層溫度變化PC2′之間的相關(guān)系數(shù)為0.4,通過了5%的顯著水平的顯著性檢驗。而前春NAO的正位相與前春東亞對流層溫度的偏暖狀態(tài)有所對應(yīng),不過它們之間的相關(guān)系數(shù)很弱,均未通過顯著性檢驗,這表明前春NAO活動和對流層中上層溫度對我國夏季降水的影響是相互獨(dú)立的。春末和夏季NAO活動可能是通過羅斯貝波的傳播影響東亞夏季氣候 (Ding and Wang,2005;Bo the et al.,2010)。

    表2 同表1,但為去除年際變化(10年以下)后的相關(guān)Table 2 Same as table 1,but for the correlation coefficients between the time series with variation less than 10 years removed by Lanczos filter

    表3 同表1,但為去除年代際變化(10年以上)后的相關(guān)Table 3 Sameas table 1,but for the correlation coeff icients between the time series with variation more than 10 years removed by Lanczos filter

    3.3.2 前春NAO、對流層中上層溫度對我國夏季降水的影響分析

    相關(guān)分析表明,前春5月NAO、“亞洲中緯型”對流層中上層溫度變化與中國夏季降水 EOF分解第1、3模態(tài)的聯(lián)系最密切,這里僅分析這兩個因子對中國夏季降水的影響。

    中國位于東亞季風(fēng)區(qū),降水主要是由水汽的輻合造成的 (周天軍等,1999;Zhou and Yu,2005)。而水汽輸送與大尺度大氣環(huán)流形勢有著密切的聯(lián)系。研究表明,西太平洋副熱帶高壓和對流層高層西風(fēng)急流是影響我國夏季氣候的重要環(huán)流系統(tǒng)(Tao and Chen,1987;Zhang et al.,2006)。本文將分別計算前春5月NAO、“亞洲中緯型”對流層中上層溫度變化這兩個序列與夏季850 hPa水汽輸送場、500 hPa位勢高度場以及200 hPa緯向風(fēng)場的回歸,來探討其對大氣環(huán)流和我國夏季降水的影響。

    圖3 前春5月負(fù)NAO指數(shù)與夏季環(huán)流場的回歸:(a)850 hPa水汽輸送 (矢量,單位:kg·m-1·s-1);(b)500 hPa位勢高度(等值線,單位:gpm);(c)200 hPa緯向風(fēng) (等值線,單位:m/s)。陰影:氣候態(tài)下500 hPa位勢高度 (b)、200 hPa緯向風(fēng)(c)的形勢Fig.3 Summer circulation fields regressed upon the negative NAO index in May:(a)Water vapor transport(vector;units:kg·m-1·s-1)at 850 hPa;(b)500-hPa geopotential height(contour,units:gpm);(c)200-hPa zonal wind(contour,units:m/s).The shaded regions indicate the climatemean 500-hPa geopotential height(b)and 200-hPa zonal wind(c)

    前春5月NAO指數(shù)與我國夏季降水第1模態(tài)和第3模態(tài)的相關(guān)最顯著,而降水的第1、3模態(tài)中心均位于長江以南地區(qū),因此前春5月NAO對中國夏季降水的影響主要是對江南和華南地區(qū)的影響。同時,負(fù)相關(guān)意味著當(dāng)NAO處于負(fù)位相時,我國南方地區(qū)降水偏多;當(dāng)NAO處于正位相時,我國南方地區(qū)降水偏少。這里只給出我國南方地區(qū)降水偏多的情況。圖3為前春5月負(fù)NAO指數(shù)與水汽輸送及環(huán)流場的回歸。從圖3a可以看出,NAO的負(fù)位相年,850 hPa水汽輸送場上,西太平洋地區(qū)存在較強(qiáng)的異常反氣旋水汽環(huán)流,反氣旋西北側(cè)的異常西南水汽輸送有利于中國南方地區(qū)的夏季降水。500 hPa位勢高度場上 (圖3b),氣候平均的西太平洋副熱帶高壓位置與一異常高壓相重合,同時,異常高壓的中心位置偏向氣候平均西太平洋副高的西南側(cè),表明西太平洋副熱帶高壓強(qiáng)度偏強(qiáng),位置偏南、偏西,使得副高西側(cè)的西南水汽輸送偏強(qiáng)。200 hPa緯向風(fēng)場上 (圖3c),長江以南地區(qū)表現(xiàn)為西風(fēng)異常,長江以北地區(qū)為東風(fēng)異常,表明西風(fēng)急流位置較常年偏南,加強(qiáng)了南方地區(qū)高空的輻散,這種環(huán)流形勢有利于中國南方地區(qū)異常偏多的降水。

    對流層溫度變化的第2模態(tài)與我國夏季降水的第3模態(tài)之間具有顯著的正相關(guān),意味著“亞洲中緯型”溫度冷異常對應(yīng)著中國夏季降水“南澇北旱”型降水分布,而暖異常對應(yīng)著“南旱北澇”型異常降水分布。同樣,這里只討論“南澇北旱”的情況。圖4為對流層溫度變化第2模態(tài)對應(yīng)的時間系數(shù)與夏季水汽輸送及環(huán)流場的回歸分布。圖4a表明,850 hPa異常水汽輸送場上,我國長江以南地區(qū)表現(xiàn)為異常偏南水汽輸送,而其北方地區(qū)為異常反氣旋水汽輸送,使得我國華北、東北地區(qū)處于水汽輻散區(qū),不利于北方地區(qū)的夏季降水;500 hPa位勢高度場上 (圖4b)最明顯的特征是,東亞地區(qū)上空存在一個強(qiáng)大的異常高壓,中心位于青藏高原以北地區(qū)。該異常高壓由“北方型”溫度冷異常引起,但異常高壓中心比對流層中上層溫度冷中心位置 (40°N,120°E)略偏西。異常冷高壓阻止了西太平洋副高的北進(jìn),同時,西太平洋副熱帶高壓強(qiáng)度也偏強(qiáng),使得副高西北側(cè)較強(qiáng)的水汽輸送維持在我國南方地區(qū)。200 hPa緯向風(fēng)場上 (圖4c),西風(fēng)急流中心位置偏南,加強(qiáng)了我國長江及其以南地區(qū)的高空輻散,這種環(huán)流形勢使得我國南方地區(qū)降水偏多,北方地區(qū)降水偏少。

    以上分析表明前春對流層中上層溫度變化、NAO的異?;顒优c我國夏季降水之間存在著明確的物理上的聯(lián)系。下面主要利用前春NAO指數(shù)及對流層中上層溫度作為預(yù)報因子,建立我國夏季降水的統(tǒng)計預(yù)測模型。

    圖4 同圖3,但為前春對流層中上層溫度REOF分解第2模態(tài)對應(yīng)的時間序列與夏季環(huán)流場的的回歸Fig.4 Same as Fig.3,but for the summer circulation fields regressed on PC2′(the second REOF mode of sp ring atmosphere temperature averaged over 500 hPa-200 hPa)

    4 預(yù)報流程及預(yù)報效果檢驗

    依據(jù)氣候系統(tǒng)具有不同時間尺度的顯著周期振蕩,魏鳳英和張先恭 (1998)提出把氣象要素序列本身生成的顯著周期函數(shù)作為預(yù)報因子加入到預(yù)報方程中,使由于預(yù)報量和預(yù)報因子的時間趨勢引起的相關(guān)不穩(wěn)定性由選入的顯著性周期因子加以調(diào)整,目的是達(dá)到更好的降水預(yù)報效果。本文在魏鳳英和張先恭 (1998)提出的預(yù)報思路的基礎(chǔ)上進(jìn)行了改進(jìn),利用前春 (3~5月)NAO指數(shù)和對流層中上層 (500~200 hPa)溫度REOF分解前4個特征向量對應(yīng)的時間系數(shù)作為預(yù)報因子,同時考慮逐月NAO指數(shù)、對流層溫度各時間序列多時間尺度的周期振蕩,建立了圖5所示的降水統(tǒng)計預(yù)報流程。

    按照圖5所示的預(yù)報流程,首先提取我國夏季降水量原始場各PC序列不同時間尺度的均值生成函數(shù)作為預(yù)報因子的一部分。均生函數(shù)是由時間序列按一定的時間間隔計算均值而派生出來的。將均生函數(shù)定義域延拓到整個數(shù)軸上,即為周期性延拓(魏鳳英,2007)。這里以 PC1為例進(jìn)行說明。圖6a為PC1的功率譜分析,從圖中可以看出,我國夏季降水第1模態(tài)的變化存在顯著的14年左右的周期振蕩以及3年左右的年際振蕩周期。因此,將周期長度為3年和14年的均生函數(shù)延拓序列從PC1序列中提取出來 (如圖6b、c所示)。對于降水的PC2和PC3序列也作同樣的處理 (圖略)。功率譜分析表明,我國夏季降水的第2模態(tài)存在2年和5年左右的顯著周期振蕩,而降水的第3模態(tài)的周期性變化與第1模態(tài)類似,第1顯著周期為10年左右,次顯著周期為2年左右。

    預(yù)報依據(jù)的另外一部分來自前春對流層中上層溫度以及NAO信號。這里給出對流層中上層溫度距平進(jìn)行REOF分解得到的 PC1′序列的功率譜分析及均生函數(shù)延拓序列。該序列的功率譜分析 (圖7a)表明,周期長度為14年處,功率譜估計值為一峰值且大大超過標(biāo)準(zhǔn)譜,而年際周期振蕩不顯著。因此,提取周期長度為14年的均生函數(shù)延拓序列,如圖7b所示。對于對流層溫度的其它 PC′序列和逐月NAO指數(shù)序列也作同樣的處理 (圖略)。結(jié)果表明,對流層中上層溫度 REOF分解的 PC2′序列以及3月NAO指數(shù)變化周期與我國夏季降水EOF分解的第1、3模態(tài)的變化周期較為一致,均存在14年和3年左右的顯著周期。其它各預(yù)報因子的變化主要存在顯著的年際振蕩周期。

    圖5 預(yù)測模型的預(yù)報流程圖Fig.5 The process diagram of the forecasting model

    圖6 我國夏季降水EOF分解第1模態(tài)對應(yīng)時間序列的功率譜分析 (a)以及周期長度為3年 (b)和14年 (c)的均生函數(shù)延拓序列。(a)虛線表示顯著性水平為0.05的紅噪音標(biāo)準(zhǔn)譜Fig.6(a) The power spectral analysis on the first PC of summer precipitation over China by the EOF decomposition;(b,c)the extension time series by mean generating function with(b)3-year and(c)14-year period.The dashed line in(a)is the red noise power spectrum with the statistical significance at 0.05 level

    根據(jù)上述方法,共得到以下預(yù)報因子:夏季降水量EOF分解前3個PC序列本身的均生函數(shù)延拓序列;前春3~5月平均對流層中上層溫度REOF分解前4個PC序列、前春逐月NAO指數(shù)序列以及它們各自對應(yīng)的顯著周期的均生函數(shù)延拓序列。使用逐步回歸方法,分別建立中國夏季降水量前3個PC序列與這些預(yù)報因子之間的回歸預(yù)測模型。

    按照上述預(yù)報流程,分別對1951~2003、1951~2004、1951~2005 和 1951~2006、1951~2007、1951~2008年的夏季降水量原始場進(jìn)行 EOF分解,通過逐步回歸,建立降水的預(yù)報方程,預(yù)報出2004~2009年間逐年的夏季降水總量,然后轉(zhuǎn)化成降水距平百分率。表4為用獨(dú)立樣本資料做出的2004~2009年預(yù)報的效果檢驗,其中ACC表示預(yù)報場與實況場之間的距平相關(guān)系數(shù),這是國際上通用的、比較客觀的評定辦法,考慮了實況值和預(yù)測值的相似程度,包括數(shù)值的大小和符號的相似程度。從表4可以看出,2004~2009年預(yù)報場和觀測場之間平均的距平相關(guān)系數(shù)達(dá)到0.335。除了2009年外,其它5年都取得了比較好的預(yù)報效果,說明基于前春對流層中上層溫度變化和NAO活動異常所建立的預(yù)報方案,對我國夏季降水具有較好的預(yù)報技巧。

    圖8給出2004和2006年中國夏季降水的預(yù)報結(jié)果。2004年夏季我國降水異常分布呈南北走向,主要降水中心位于黃淮地區(qū)、長江中游以及西南、華北的部分地區(qū) (圖8a)。圖8c表明,該預(yù)測模型能夠較好地預(yù)報出黃淮地區(qū)和長江中游的多雨中心以及華南地區(qū)的少雨中心,但是江淮地區(qū)降水的預(yù)報與實況存在較大的誤差。2006年夏季,我國東部自南向北“+-+”的降水異常分布基本上預(yù)報出來了,如華南地區(qū)的多雨中心和長江中下游地區(qū)的少雨中心,但是預(yù)報的北方雨區(qū)的位置較實況更偏北一些。

    圖7 前春對流層中上層溫度REOF分解第1模態(tài)對應(yīng)時間序列的功率譜分析 (a)以及周期長度為14年的均生函數(shù)延拓序列 (b)。(a)虛線表示顯著性水平為0.05的紅噪音標(biāo)準(zhǔn)譜Fig.7 (a)The power spectral analysis on the first PC of sp ring atmosphere temperature averaged over 500 hPa-200 hPa by the REOF decomposition;(b)the extension time series by mean generating function with 14-year period.The dashed line in(a)is the red noise power spectrum with the statistical significance at 0.05 level

    表4 2004~2009年中國夏季降水預(yù)測效果檢驗Table 4 The forecasting verification of summer precipitation over China during 2004-2009

    對比2004~2009年我國夏季降水的觀測場和預(yù)報場還可以發(fā)現(xiàn),該預(yù)測模型基本上能夠把握我國夏季旱澇趨勢的分布特征,這與該統(tǒng)計預(yù)測模型中預(yù)報因子的選取有著密切的聯(lián)系。由于對流層溫度是一個比較綜合性的物理因子,前春對流層溫度異常已經(jīng)包含了來自海洋、高原積雪等各種強(qiáng)迫因子的變化,因此,只基于前春對流層中上層溫度和NAO變化的降水統(tǒng)計預(yù)測模型即可以預(yù)報出我國夏季降水分布的大部分信息。但是該預(yù)測模型還存在不足之處,比如預(yù)測的降水強(qiáng)度比觀測偏弱,并且未能預(yù)測出局地降水的更多細(xì)節(jié)特征。這是因為對流層溫度變化所描述的是大尺度海陸溫差的變化,故僅能反映與之聯(lián)系的大范圍降水異常。另外,前文的相關(guān)分析表明,我國夏季降水 EOF分解的第2模態(tài)與前春對流層溫度及NAO變化之間的相關(guān)很弱,導(dǎo)致該預(yù)測模型對降水第2模態(tài)的預(yù)報能力很低。這意味著要更好地預(yù)報我國夏季降水的分布,還需要在該預(yù)測模型中考慮影響降水第2模態(tài)的主要預(yù)報因子。

    圖8 2004年 (a、c)和2006年 (b、d)夏季降水距平百分率 (%):(a、b)實況;(c、d)預(yù)報Fig.8 The summer precipitation anomaly percentage(%)in(a,c)2004 and(b,d)2006:(a,b)Observation;(c,d)forecast

    5 結(jié)論和討論

    (1)我國夏季降水大范圍分布趨勢與前春5月NAO的變化密切相關(guān),“南多 (少)北少 (多)型”夏季降水主要與“亞洲中緯型”對流層中上層溫度變化有關(guān),而“南北少 (多)中間多 (少)型”夏季降水與前春NAO及對流層中上層溫度之間均無顯著相關(guān)。

    (2)前春NAO活動異常對我國夏季降水的影響既有年際尺度上的,也有年代際尺度上的,而前春對流層中上層溫度對我國夏季降水的影響主要是年代際尺度上的。

    (3)前春5月NAO活動對我國南方地區(qū)夏季降水的影響機(jī)理為:NAO的負(fù)位相年,夏季西風(fēng)急流位置偏南,西太平洋副熱帶高壓偏西、偏南,來自西太平洋的西南水汽輸送偏強(qiáng),造成我國夏季江南、華南地區(qū)降水偏多;相反,NAO正位相年,我國南方地區(qū)降水偏少。前春“亞洲中緯型”對流層溫度冷異常對我國夏季“南澇北旱”型降水格局的影響機(jī)理為:當(dāng)東北、華北地區(qū)上空異常偏冷時,我國北方大部分地區(qū)為異常高壓控制,使得西太平洋副高西側(cè)位置偏南,不利于副高西側(cè)水汽的向北輸送,同時,夏季西風(fēng)急流位置也偏南,使得南方地區(qū)降水偏多,而北方地區(qū)夏季降水偏少。相反,前春“亞洲中緯型”對流層溫度暖異常對應(yīng)著我國夏季降水“南旱北澇”型的異常分布。

    (4)本文使用的預(yù)測模型,除了考慮降水本身不同時間尺度的周期振蕩,還將前春NAO指數(shù)及對流層中上層溫度變化各時間序列不同顯著周期的均生函數(shù)延拓序列提取出來,加入預(yù)報方程,很大程度上排除了噪音,取得了較好的降水預(yù)報效果。

    本文基于前春NAO、對流層溫度變化與夏季降水異常的關(guān)系,建立了一個統(tǒng)計預(yù)測模型,對于我國夏季降水具有較好的預(yù)報能力。基于春季NAO和對流層溫度異常對夏季降水的統(tǒng)計預(yù)報,可以在汛期旱澇預(yù)測的滾動訂正中發(fā)揮作用。

    本文建立的降水統(tǒng)計預(yù)測模型主要基于再分析資料中的前春對流層中上層溫度變化,該模型的基礎(chǔ)是自春季到夏季對流層中上層溫度異常的持續(xù)性。而夏季降水與同期對流層中上層溫度之間的對應(yīng)關(guān)系無疑好于前期對流層溫度。如果能夠提前獲得理想的夏季對流層溫度資料,將會大大提高降水預(yù)報效果。數(shù)值模式,尤其是海氣耦合模式被認(rèn)為是目前最常用的季節(jié)預(yù)報工具。當(dāng)前的耦合氣候模式,對東亞表層溫度和降水的直接預(yù)報能力很低,但對包括溫度場在內(nèi)的東亞高層環(huán)流的預(yù)報能力則很強(qiáng) (Wang et al.,2008)。即使當(dāng)前的氣候模式難以直接準(zhǔn)確地預(yù)測降水的變化,但只要它對東亞對流層溫度的變化具備一定預(yù)測能力,就可以在此基礎(chǔ)上,進(jìn)一步基于“對流層溫度和降水變化”的統(tǒng)計關(guān)系,采用統(tǒng)計降尺度的辦法來預(yù)測降水變化。因此,根據(jù)“對流層溫度—表面降水”變化的統(tǒng)計關(guān)系,基于耦合模式預(yù)報的對流層溫度變化來預(yù)報降水的變化是一條切實可行的途徑。

    References)

    Bothe O,Fraedrich K,Zhu X H.2010.The large-scale circulations and summer drought and wetnesson the Tibetan plateau[J].International Journal of Climatology,30:844-855.

    Ding Q,Wang B.2005.Circumglobal teleconnection in the Northern Hemisphere summer[J].J.Climate,18:3483-3505.

    丁一匯.2004.我國短期氣候預(yù)測業(yè)務(wù)系統(tǒng) [J].氣象,30(12):11-17.Ding Y H.2004.China Operational Short-Range Climate Prediction System[J].Meteorological Monthly(in Chinese),30(12):11-17.

    龔道溢,朱錦紅,王紹武.2002.長江流域夏季降水與前期北極濤動的顯著相關(guān) [J].科學(xué)通報,47(7):546-549. Gong D Y,Zhu J H,Wang S W.2002.Significant relationship between sp ring AO and the summer rainfall along the Yangtze River[J].Chinese Science Bulletin(in Chinese),47(7):546-549.

    黃嘉佑,劉舸,趙昕奕.2004.副高、極渦因子對我國夏季降水的影響[J].大氣科學(xué),28(4):517-526. Huang Jiayou,Liu Ge,Zhao Xiyi.2004.The influence of subtropical high indexes and polar vortex indexes on the summertime precipitation in China[J].Chinese Journal of A tmospheric Sciences(in Chinese),28(4):517-526.

    韓雪,魏鳳英.2010.中國東部夏季降水與東亞垂直環(huán)流結(jié)構(gòu)及其預(yù)測試驗 [J].大氣科學(xué),34(3):533-547. Han Xue,Wei Fengying.2010.The influence of vertical atmospheric circulation pattern of over East Asia on summer precipitation in the east of China and its forecasting test[J].Chinese Journalof A tmospheric Sciences(in Chinese),34(3):533-547.

    金祖輝,陶詩言.1999.ENSO循環(huán)與中國東部地區(qū)夏季和冬季降水關(guān)系的研究[J].大氣科學(xué),23(6):663-671.Jin Zuhui,Tao Shiyan.1999.A study on the relationships between ENSO cycle and rainfalls during summer and winter in eastern China[J].Chinese Journal of A tmospheric Sciences(in Chinese),23(6):663-671.

    Kalnay E,Kanamitsu M,Kistler R,et al.1996.The NCEP/NCAR 40-year reanalysis project[J].Bull.Amer.Meteo r.Soc.,77:437-472.

    Li J,Yu R C,Zhou T J,et al.2005.Why is there an early sp ring cooling shift dow nstream of the Tibetan Plateau?[J].J.Climate,18(22):4660-4668.

    Li J,Yu R C,Zhou T J.2008.Teleconnection between NAO and climate dow nstream of the Tibetan Plateau[J].J.Climate,21:4680-4690.

    No rth G R,Bell T L,Cahalan R F,et al.1982.Samp ling errors in the estimation of empirical o rthogonal functions[J].Mon.Wea.Rea.,110:699-706.

    彭京備,陳烈庭,張慶云.2005.青藏高原異常雪蓋和 ENSO的多尺度變化及其與中國夏季降水的關(guān)系 [J].高原氣象,24(3):366-377. Peng Jingbei,Chen Lieting,Zhang Qingyun.2005.Multi-scale variations of snow cover over QXP and tropical Pacific SST and their influences on summer rainfall in China[J].Plateau Meteo rology(in Chinese),24(3):366-377.

    施能.2009.氣象統(tǒng)計預(yù)報 [M].北京:氣象出版社,128-143.Shi N.2009.Statistical Forecast of Meteorology(in Chinese)[M].Beijing:China Meteorological Press,128-143.

    Tao S,Chen L.1987.A review of recent research on the East A-sia summer monsoon inChina[M]//Chang C P,Krishnam urti T N.Monsoon Meteo rology.Oxfo rd:Oxfo rd university Press,60-92.

    王斌.2008.認(rèn)識和預(yù)報亞洲季風(fēng)氣候:前沿突破點(diǎn)和展望[J].氣象學(xué)報,66:653-669. Wang B.2008.Thrusts and prospects on understanding and predicting Asian monsoon climate[J].Acta Meteo rologica Sinica(in Chinese),66(5):653-669.

    Wang B,Lee J Y,Kang IS,et al.2008. Advance and prospectus of seasonal prediction:A ssessment of the APCC/CliPAS 14-model ensemble retrospective seasonal prediction(1980-2004)[J].Climate Dynamics,doi:10.1007/s00382-008-0460-0.

    王永波,施能.2001.近45a冬季北大西洋濤動異常與我國氣候的關(guān)系[J].南京氣象學(xué)院學(xué)報,24(3):315-322.Wang YB,Shi N.2001.Relation of North A tlantic Oscillation anomaly to China climate during 1951-1995[J].Journalof Nanjing Institute of Meteorology(in Chinese),24(3):315-322.

    魏鳳英,張先恭.1998.中國夏季降水趨勢分布的一個客觀預(yù)報方法[J].氣候與環(huán)境研究,3(3):218-226. Wei F Y,Zhang X G.1998.An objective predicting method for seasonal precipitation anomaly distribution in summer over China[J].Climatic and Environmental Research(in Chinese),3(3):218-226.

    魏鳳英.2006.長江中下游夏季降水異常變化與若干強(qiáng)迫因子的關(guān)系[J].大氣科學(xué),30(2):202-211. Wei Fengying.2006.Relationships between precipitation anomaly over the middle and lower reachesof the Changjiang River in summer and several forcing factors[J].Chinese Journal of A tmospheric Sciences(in Chinese),30(2):202-211.

    魏鳳英.2007.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù) (第2版)[M].北京:氣象出版社.Wei F Y.2007.Climate Statistical Diagnosing and Prediction(in Chinese)[M].2nd ed.Beijing:China Meteorological Press.

    魏鳳英,黃嘉佑.2010.大氣環(huán)流降尺度因子在中國東部夏季降水預(yù)測中的作用 [J].大氣科學(xué),34(1):202-212. Wei F Y,Huang J Y.2010.A study of dow nscaling factors of atmospheric circulations in the prediction model of summer precipitation in Easter China[J].Chinese Journal of A tmospheric Sciences(in Chinese),34(1):202-212.

    武炳義,張人禾,D'A rrigo Rosanne.2008.北極偶極子異常與中國東北夏季降水[J].科學(xué)通報,53(12):1422-1428. Wu Bingyi,Zhang Renhe,D'A rrigo Rosanne.2008.A rctic dipole anomaly and summer rainfall in no rtheast China[J].Chinese Science Bulletin(in Chinese),53(12):1422-1428.

    Wu B,Zhou T J,Li T.2009.Seasonally evolving dominant interannual variability modes of East A sian climate[J].J.Climate,22:2992-3005

    Wu B,Li T,Zhou T J.2010.Relative contributionsof the Indian O-cean and local SST anomalies to the maintenance of the western North Pacific anomalous anticyclone duringdecaying summer[J].J.Climate,23:2974-2986.

    吳統(tǒng)文,錢正安.2000.青藏高原冬春積雪異常與中國東部地區(qū)夏季降水關(guān)系的進(jìn)一步分析[J].氣象學(xué)報,58(5):570-581.Wu Tongw en,Qian Zheng'an.2000.Further analysis of the linkage between winter and sp ring snow dep th anomaly over Qinghai-Xizang Plateau and summer rainfall of eastern China[J].Acta Meteo rologica Sinica(in Chinese),58(5):570-581.

    Xin X G,Yu R C,Zhou T J,et al.2006.Drought in late sp ring of South China in recent decades[J].J.Climate,19(13):3197-3206.

    辛?xí)愿?宇如聰,周天軍.2009.中國東南部4~5月年代際干旱的南移特征及潛熱影響的數(shù)值模擬 [J].大氣科學(xué),33(6):1165-1173. Xin Xiaoge,Yu Rucong,Zhou Tianjun.2009.Southward movement of the decadal drought in southeastern China during Ap ril-May and numerical simulation of the effect of the condensation heating[J].Chinese Journal of A tmospheric Sciences(in Chinese),33(6):1165-1173.

    Xin X G,Zhou T J,Yu R C.2010.Increased Tibetan Plateau snow depth:An indicator of the connection between enhanced winter NAO and late-sp ring tropospheric cooling over East Asia[J].Advances in A tmospheric Sciences,27:788-794.

    楊蓮梅,張慶云.2008.北大西洋濤動對新疆夏季降水異常的影響[J].大氣科學(xué),32(5):1187-1196. Yang Lianmei,Zhang Qingyun.2008.Effects of the North A tlantic Oscillation on the summer rainfall anomalies in Xinjiang[J].Chinese Journal of A tmospheric Sciences(in Chinese),32(5):1187-1196.

    Yu R C,Wang B,Zhou T J.2004.Tropospheric cooling and summer monsoon weakening trend over East Asia[J].Geophys.Res.Lett.,31,L22212,doi:10.1029/2004GL021270.

    Yu R C,Zhou T J.2004.Impacts of winter-NAO on March cooling trends over subtropical Eurasia continent in the recent half century[J]. Geophys.Res.Lett.,31,L 12204,doi:10.1029/2004GL 019814.

    Yu R C,Zhou T J.2007.Seasonality and three-dimensional structure of interdecadal change in East Asian monsoon[J].J.Climate,20:5344-5355.

    宇如聰,周天軍,李建,等.2008.中國東部氣候年代際變化三維特征的研究進(jìn)展 [J].大氣科學(xué),32(4):893-905. Yu Rucong,Zhou Tianjun,Li Jian,et al.2008.Progress in the studies of three-dimensional structure of interdecadal climate change over eastern China[J].Chinese Journal of A tmospheric Sciences(in Chinese),32(4):893~905.

    張慶云,呂俊梅,楊蓮梅,等.2007.夏季中國降水型的年代際變化與大氣內(nèi)部動力過程及外強(qiáng)迫因子關(guān)系[J].大氣科學(xué),31(6):1290-1300.Zhang Q Y,LüJ M,Yang L M,et al.2007.The interdecadal variation of precipitation pattern over China during summer and its relationship with the atmospheric internal dynamic processes and extra-forcing factors[J].Chinese Journal of A tmospheric Sciences(in Chinese),31(6):1290-1300.

    Zhang R H,Sumi A,Kimoto M.1999.A diagnostic study of the impact ofon the precipitation in China[J].Adv.A tmos.Sci.,16:229-241.

    Zhang Y C,Kuang X Y,Guo WD,et al.2006.Seasonal evolution of the upper-tropospheric westerly jet co re over East Asia[J].Geophys. Res. Lett., 33, L 11708, doi: 10.1029/2006GL 026377.

    趙平,周秀驥.2006.近40年我國東部降水持續(xù)時間和雨帶移動的年代際變化[J].應(yīng)用氣象學(xué)報,17(5):548-556. Zhao Ping,Zhou Xiuji.2006.Decadal variability of rainfall persistence time and rain-belt shift over eastern China in recent 40 years[J].Journal of App lied Meteorological Sciences,17(5):548-556.

    趙振國.1996.厄爾尼諾現(xiàn)象對北半球大氣環(huán)流和中國降水的影響[J].大氣科學(xué),20(4):422-428. Zhan Zhenguo.1996.Impact ofevents on atmospheric circulations in the Northern Hemisphere and precipitation in China[J].Chinese Journal of A tmospheric Sciences(Scientia A tmospherica Sinica)(in Chinese),20(4):422-428.

    趙振國.1999.中國夏季旱澇及環(huán)境場 [M].北京:氣象出版社,297pp.Zhao Z G.1999.Summertime Floods and Droughts in China and the Associated Circulations(in Chinese)[M].Beijing:China Meteorological Press,297pp.

    周天軍,張學(xué)洪,王紹武.1999.全球水循環(huán)的海洋分量研究 [J].氣象學(xué)報,57:264-282. Zhou T J,Zhang X H,Wang SW.1999.The air-sea freshwater exchange derived from NCEP/NCAR reanalysis data[J].Acta Meteorologica Sinica(in Chinese),57:264-282.

    Zhou T J,Yu R C.2005.A tmospheric water vapor transport associated with typical anomalous summer rainfall patterns in China[J]. J. Geophys. Res.,110,D08104, doi: 10.1029/2004JD005413.

    Zhou T J,Gong D Y,Li J,et al.2009a.Detecting and understanding the multi-decadal variability of the East A sian summer monsoon:Recent progress and state of affairs[J].Meteo rologische Zeitschrift,18(4):455-467.

    Zhou T J,Wu B,Wang B.2009b.How well do A tmospheric General Circulation Models cap ture the leading modesof the interannual variability of Asian-Australian monsoon?[J].J.Climate,22:1159-1173.

    Zhou T J,Zhang J.2009.Harmonious inter-decadal changesof July-August upper tropospheric temperature across the North A tlantic,Eurasian continent,and North Pacific[J].Advances in A tmospheric Sciences,26,656-665.

    Zhou T J,Zou L W.2010.Understanding the predictability of East Asian summermonsoon from the rep roduction of land-sea thermal contrast change in AMIP-type simulation[J].J.Climate,23,6009-6026,doi:10.1175/2010JCL I3546.1.

    朱艷峰,張博,陳隆勛.2010.夏季青藏高原與其東部平原的熱力差異對中國降水的影響 [J].科學(xué)通報,55(6):483-489. Zhu Y F,Zhang B,Chen L X.2010.Thermal difference betw een the Tibetan Plateau and the plain east of Plateau and its influence on rainfall over China in the summer[J].Chinese Science Bulletin,55,doi:10.1007/s11434-009-0613-5.

    A Predictive Model for Summer Precipitation over China Based on Upper Tropospheric Temperature and North Atlantic Oscillation in the Preceding Spring

    FENG Lei1,2,WEI Fengying3,and ZHU Yanfeng4
    1StateKeyLaboratoryofNumericalModelingforAtmosphericSciencesandGeophysicalFluidDynamics,InstituteofAtmosphericPhysics,ChineseAcademyofSciences,Beijing100029
    2GraduateUniversityofChineseAcademyofSciences,Beijing100049
    3ChineseAcademyofMeteorologicalSciences,ChinaMeteorologicalAdministration,Beijing100081
    4NationalClimateCenter,ChinaMeteorologicalAdministration,Beijing100081

    Based on the analyses of relationship between summer precipitation and upper tropospheric temperature,No rth A tlantic Oscillation(NAO)in the preceding sp ring,a statistic predictivemodel for summer precipitation over China is proposed and the forecast results for years 2004-2009 are tested.The results show that the large scale precipitation over China mainly changeswith NAO in May,and themeridional dipole pattern of precipitation is associated with NAO in May and the anomalous cooling during 500 hPa-200 hPa in the preceding sp ring over No rth China and Northeast China.The NAO influences the summer precipitation over China on both interannual and interdecadal time scales,w hile the upper tropospheric temperature influences the precipitation mainly on interdecadal time scale.The anomalous condition of sp ring NAO and upper tropospheric temperature determines the pattern and magnitude of summer precipitation over China via large scale circulation and water vapo rtransport.A predictivemodel for summer precipitation over China is proposed based on both this statistic and physical relationship,with the first three princip le components(PCs)of summer precipitation as the predictands,and with the first four PCs of upper tropospheric temperature over Asia(10°N-50°N,60°E-130°E)and themonthly NAO index in the preceding sp ring as the predictors.The periodsof these series with different time scales are also considered to remove the noise and extract useful info rmation.The averaged anomaly correlation coefficient for years 2004-2009 is 0.335,indicting a good skill in fo recasting the summer precipitation over China.

    upper tropospheric temperature,North A tlantic Oscillation(NAO),summer precipitation over China,statistic predictive model

    1006-9895(2011)05-0963-14

    P461

    A

    馮蕾,魏鳳英,朱艷峰.2011.基于前春對流層溫度和北大西洋濤動的中國夏季降水統(tǒng)計預(yù)測模型 [J].大氣科學(xué),35(5):963-976. Feng Lei,Wei Fengying,Zhu Yanfeng.2011.A predictivemodel for summer precipitation over China based on upper tropospheric temperature and North A tlantic Oscillation in the preceding sp ring[J].Chinese Journal of A tmospheric Sciences(in Chinese),35(5):963-976.

    2010-09-20,2010-12-09收修定稿

    財政部/科技部公益類行業(yè)專項GYH Y200706010、GYHY200906016

    馮蕾,女,1982年出生,博士研究生,主要從事氣候診斷與預(yù)測研究。E-mail:leifeng@mail.iap.ac.cn

    猜你喜歡
    對流層特征向量降水
    二年制職教本科線性代數(shù)課程的幾何化教學(xué)設(shè)計——以特征值和特征向量為例
    克羅內(nèi)克積的特征向量
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    郴州地區(qū)對流層頂氣候概況
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    一類特殊矩陣特征向量的求法
    EXCEL表格計算判斷矩陣近似特征向量在AHP法檢驗上的應(yīng)用
    實時干涉測量中對流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    降水現(xiàn)象儀模擬軟件設(shè)計與實現(xiàn)
    成都地區(qū)2005~2015年對流層NO2柱濃度趨勢與時空分布
    成年女人看的毛片在线观看| 国产精品99久久久久久久久| 国产成人精品一,二区 | 欧美高清性xxxxhd video| 国产亚洲欧美98| 毛片女人毛片| 午夜福利在线观看吧| 国产视频首页在线观看| 亚洲综合色惰| 久久热精品热| 欧美日韩一区二区视频在线观看视频在线 | 欧洲精品卡2卡3卡4卡5卡区| 国产成人aa在线观看| 观看美女的网站| 热99在线观看视频| 中文欧美无线码| 亚洲精品久久久久久婷婷小说 | 可以在线观看的亚洲视频| 性插视频无遮挡在线免费观看| 深爱激情五月婷婷| 观看美女的网站| 欧美一区二区亚洲| a级毛片免费高清观看在线播放| 久久中文看片网| 日产精品乱码卡一卡2卡三| 国产爱豆传媒在线观看| 女同久久另类99精品国产91| 婷婷精品国产亚洲av| 91狼人影院| 又黄又爽又刺激的免费视频.| 亚洲av电影不卡..在线观看| 美女 人体艺术 gogo| 99久国产av精品国产电影| 亚洲aⅴ乱码一区二区在线播放| а√天堂www在线а√下载| 久久久久久国产a免费观看| 国产精品,欧美在线| 欧美bdsm另类| 久久精品夜夜夜夜夜久久蜜豆| 国产 一区精品| 深夜a级毛片| 少妇被粗大猛烈的视频| 91av网一区二区| 国产一区二区在线观看日韩| 亚洲激情五月婷婷啪啪| 亚洲精品日韩在线中文字幕 | 男人的好看免费观看在线视频| 国模一区二区三区四区视频| 亚洲av二区三区四区| 免费看光身美女| 婷婷亚洲欧美| 亚洲色图av天堂| 夜夜看夜夜爽夜夜摸| 男的添女的下面高潮视频| 国产一区二区三区av在线 | 久久精品夜夜夜夜夜久久蜜豆| 少妇人妻一区二区三区视频| 3wmmmm亚洲av在线观看| 国产精品一及| 我的女老师完整版在线观看| 久久国产乱子免费精品| 男人舔女人下体高潮全视频| 男女做爰动态图高潮gif福利片| 婷婷六月久久综合丁香| 亚洲自偷自拍三级| 亚洲av免费在线观看| 综合色av麻豆| 午夜亚洲福利在线播放| 卡戴珊不雅视频在线播放| 在线观看美女被高潮喷水网站| 亚洲最大成人av| av在线天堂中文字幕| 最近中文字幕高清免费大全6| 99riav亚洲国产免费| 麻豆成人午夜福利视频| 亚洲最大成人手机在线| 亚洲国产精品成人久久小说 | 不卡一级毛片| 久久99蜜桃精品久久| eeuss影院久久| 男女做爰动态图高潮gif福利片| 爱豆传媒免费全集在线观看| 国产一区二区在线av高清观看| 天美传媒精品一区二区| 精品久久久久久久末码| 国产真实乱freesex| 亚洲精品456在线播放app| 久久草成人影院| 亚洲第一电影网av| 乱码一卡2卡4卡精品| 国产av麻豆久久久久久久| 欧美人与善性xxx| 内射极品少妇av片p| 九九在线视频观看精品| 1000部很黄的大片| 2022亚洲国产成人精品| 五月伊人婷婷丁香| 美女国产视频在线观看| 久久久久久大精品| 最近中文字幕高清免费大全6| 午夜福利在线观看吧| 麻豆国产av国片精品| 免费在线观看成人毛片| 亚洲av电影不卡..在线观看| 亚洲在线自拍视频| 国产激情偷乱视频一区二区| 91精品一卡2卡3卡4卡| 日韩,欧美,国产一区二区三区 | 欧美激情久久久久久爽电影| 免费无遮挡裸体视频| 欧美日韩综合久久久久久| 边亲边吃奶的免费视频| 国产亚洲5aaaaa淫片| 少妇丰满av| 国产精品国产三级国产av玫瑰| 国产一区二区激情短视频| 成人漫画全彩无遮挡| 亚洲精品影视一区二区三区av| 禁无遮挡网站| 色视频www国产| 亚洲在线观看片| 麻豆成人av视频| 精品不卡国产一区二区三区| 日韩欧美 国产精品| av在线观看视频网站免费| 亚洲av中文av极速乱| 男人舔女人下体高潮全视频| 国产精品三级大全| 免费大片18禁| 亚洲aⅴ乱码一区二区在线播放| 亚洲av免费高清在线观看| 日韩一区二区三区影片| 在线播放国产精品三级| 赤兔流量卡办理| 黄色一级大片看看| 成人毛片60女人毛片免费| 97人妻精品一区二区三区麻豆| 免费电影在线观看免费观看| 97在线视频观看| av在线播放精品| 精品人妻熟女av久视频| 蜜桃久久精品国产亚洲av| 精品久久久噜噜| 久久久久国产网址| 天堂中文最新版在线下载 | 午夜福利在线观看免费完整高清在 | av免费在线看不卡| 97在线视频观看| 亚洲五月天丁香| 国产精品美女特级片免费视频播放器| 极品教师在线视频| 国产 一区 欧美 日韩| 欧美变态另类bdsm刘玥| 国产黄片美女视频| 人人妻人人澡欧美一区二区| av免费在线看不卡| 免费观看人在逋| 伦理电影大哥的女人| 国产高清不卡午夜福利| 欧美又色又爽又黄视频| 丰满的人妻完整版| 国产精品三级大全| 日本欧美国产在线视频| 99久国产av精品| 大又大粗又爽又黄少妇毛片口| 欧美+日韩+精品| 国产精品久久久久久久电影| 久久欧美精品欧美久久欧美| 久久精品国产清高在天天线| 我要看日韩黄色一级片| 最后的刺客免费高清国语| 日本撒尿小便嘘嘘汇集6| 中文在线观看免费www的网站| 国产伦一二天堂av在线观看| 亚洲色图av天堂| 天天躁夜夜躁狠狠久久av| 亚洲成人精品中文字幕电影| 有码 亚洲区| 男女视频在线观看网站免费| 乱系列少妇在线播放| 国产91av在线免费观看| kizo精华| 欧美一区二区精品小视频在线| 九九久久精品国产亚洲av麻豆| 成人漫画全彩无遮挡| 精品一区二区三区视频在线| 亚洲人成网站在线播| 亚洲三级黄色毛片| 99久久九九国产精品国产免费| 国产乱人视频| 色综合亚洲欧美另类图片| 中国国产av一级| 1024手机看黄色片| 亚洲久久久久久中文字幕| 在线国产一区二区在线| 国产高清不卡午夜福利| 免费观看a级毛片全部| 精品久久国产蜜桃| 久久久久九九精品影院| 嘟嘟电影网在线观看| 久久这里有精品视频免费| 美女脱内裤让男人舔精品视频 | 久久久久久久午夜电影| 18禁裸乳无遮挡免费网站照片| 国产一级毛片在线| 少妇的逼水好多| 人妻夜夜爽99麻豆av| 免费av毛片视频| 久久韩国三级中文字幕| 国产一级毛片在线| 中国美白少妇内射xxxbb| 亚洲欧美精品自产自拍| 日韩精品有码人妻一区| 综合色丁香网| 麻豆精品久久久久久蜜桃| 可以在线观看毛片的网站| 精品不卡国产一区二区三区| 26uuu在线亚洲综合色| 最近中文字幕高清免费大全6| 欧美3d第一页| 人人妻人人看人人澡| 一区二区三区高清视频在线| 国产精品一及| 毛片一级片免费看久久久久| 如何舔出高潮| 国产免费男女视频| 久久人人精品亚洲av| 嫩草影院精品99| 能在线免费观看的黄片| 啦啦啦韩国在线观看视频| 国产一级毛片七仙女欲春2| av免费在线看不卡| 国产精品久久视频播放| 成人二区视频| 麻豆成人av视频| 人人妻人人澡欧美一区二区| 九色成人免费人妻av| 亚洲国产精品成人久久小说 | 久久精品夜夜夜夜夜久久蜜豆| а√天堂www在线а√下载| 成年免费大片在线观看| 亚洲欧美成人精品一区二区| 91午夜精品亚洲一区二区三区| 一边亲一边摸免费视频| 久久国内精品自在自线图片| 真实男女啪啪啪动态图| 亚洲av二区三区四区| 在线观看一区二区三区| 一夜夜www| 欧美成人精品欧美一级黄| 麻豆成人av视频| 精品人妻视频免费看| 两个人视频免费观看高清| 晚上一个人看的免费电影| 精品少妇黑人巨大在线播放 | 亚洲真实伦在线观看| 97在线视频观看| 中文字幕精品亚洲无线码一区| 日韩在线高清观看一区二区三区| 日日啪夜夜撸| 黄色视频,在线免费观看| 少妇人妻精品综合一区二区 | 亚洲激情五月婷婷啪啪| 变态另类丝袜制服| 午夜福利视频1000在线观看| 中文在线观看免费www的网站| av天堂中文字幕网| 国产成人91sexporn| 毛片女人毛片| 天堂av国产一区二区熟女人妻| 亚洲精品日韩在线中文字幕 | 成人午夜精彩视频在线观看| 欧美一区二区精品小视频在线| 亚洲内射少妇av| 最近2019中文字幕mv第一页| 成人美女网站在线观看视频| av.在线天堂| 国产成人a∨麻豆精品| 免费观看的影片在线观看| 老女人水多毛片| 天堂中文最新版在线下载 | 亚洲成人久久性| 男女边吃奶边做爰视频| 男人和女人高潮做爰伦理| 午夜久久久久精精品| 精品不卡国产一区二区三区| 久久这里只有精品中国| 欧美成人a在线观看| 久久人人爽人人片av| 看片在线看免费视频| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 国产成人影院久久av| 国产av一区在线观看免费| 国产精品.久久久| 色噜噜av男人的天堂激情| 91久久精品国产一区二区成人| 99热这里只有精品一区| 国产精品国产高清国产av| 男女下面进入的视频免费午夜| 亚洲国产精品成人综合色| 亚洲国产精品国产精品| 国语自产精品视频在线第100页| 日韩精品有码人妻一区| 18禁在线无遮挡免费观看视频| 99久久成人亚洲精品观看| 在线观看免费视频日本深夜| 色吧在线观看| 亚洲最大成人手机在线| 精品国内亚洲2022精品成人| 亚洲一区二区三区色噜噜| 色哟哟·www| 美女 人体艺术 gogo| 欧美日韩国产亚洲二区| 亚洲精品日韩av片在线观看| 国产视频首页在线观看| 国产人妻一区二区三区在| 午夜免费男女啪啪视频观看| 91aial.com中文字幕在线观看| 久久6这里有精品| 最近最新中文字幕大全电影3| 麻豆久久精品国产亚洲av| a级毛片免费高清观看在线播放| 欧美3d第一页| 国产高清不卡午夜福利| 亚洲内射少妇av| 国产人妻一区二区三区在| 国产成人一区二区在线| 3wmmmm亚洲av在线观看| 国模一区二区三区四区视频| 日本一本二区三区精品| 人妻夜夜爽99麻豆av| 精品久久久久久久久av| 99热全是精品| a级一级毛片免费在线观看| 三级男女做爰猛烈吃奶摸视频| 国产av麻豆久久久久久久| av女优亚洲男人天堂| 国产国拍精品亚洲av在线观看| 日本免费一区二区三区高清不卡| 国内精品久久久久精免费| 99精品在免费线老司机午夜| 国产高潮美女av| 九九热线精品视视频播放| 亚洲精品456在线播放app| 亚洲人成网站在线播放欧美日韩| 国产亚洲91精品色在线| 亚洲一区高清亚洲精品| 人体艺术视频欧美日本| 我要看日韩黄色一级片| 亚洲第一电影网av| 神马国产精品三级电影在线观看| 久久6这里有精品| 联通29元200g的流量卡| 黄片wwwwww| 精品一区二区三区视频在线| 欧美激情国产日韩精品一区| 热99在线观看视频| 简卡轻食公司| 少妇高潮的动态图| 亚洲国产欧洲综合997久久,| 青春草国产在线视频 | 国产不卡一卡二| 99在线视频只有这里精品首页| 伊人久久精品亚洲午夜| 国产精品人妻久久久影院| 麻豆精品久久久久久蜜桃| 免费无遮挡裸体视频| 亚洲内射少妇av| 欧美极品一区二区三区四区| 亚洲中文字幕一区二区三区有码在线看| 久久久久久大精品| 欧美高清性xxxxhd video| 国产成人精品久久久久久| 嫩草影院精品99| 国内久久婷婷六月综合欲色啪| 一进一出抽搐gif免费好疼| 我要搜黄色片| 国产精品麻豆人妻色哟哟久久 | 高清在线视频一区二区三区 | 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 性欧美人与动物交配| av.在线天堂| 青青草视频在线视频观看| 真实男女啪啪啪动态图| 国产极品天堂在线| 亚洲精品久久国产高清桃花| 中国美白少妇内射xxxbb| 色吧在线观看| 国产精品伦人一区二区| 午夜福利在线观看免费完整高清在 | 99久久精品国产国产毛片| 国产伦在线观看视频一区| 婷婷六月久久综合丁香| 中文字幕精品亚洲无线码一区| 国产精品.久久久| 色综合站精品国产| 午夜精品国产一区二区电影 | 国产伦理片在线播放av一区 | 毛片一级片免费看久久久久| 大又大粗又爽又黄少妇毛片口| 久久精品久久久久久噜噜老黄 | 国产v大片淫在线免费观看| 1000部很黄的大片| 国产精品久久久久久av不卡| 欧美最黄视频在线播放免费| av国产免费在线观看| 国产人妻一区二区三区在| 看十八女毛片水多多多| 丰满乱子伦码专区| 亚洲人成网站在线播| 免费看a级黄色片| 好男人在线观看高清免费视频| 国产精品.久久久| 舔av片在线| 国产色爽女视频免费观看| 精华霜和精华液先用哪个| 老司机影院成人| 成年版毛片免费区| 一级二级三级毛片免费看| 国产精品久久久久久久久免| 久久精品91蜜桃| 日日摸夜夜添夜夜爱| 小蜜桃在线观看免费完整版高清| 亚洲欧洲国产日韩| 非洲黑人性xxxx精品又粗又长| 女同久久另类99精品国产91| 一级二级三级毛片免费看| 日韩大尺度精品在线看网址| 18禁黄网站禁片免费观看直播| 最近2019中文字幕mv第一页| 国产老妇女一区| 国产精品女同一区二区软件| 中文字幕久久专区| avwww免费| 国产人妻一区二区三区在| 欧美人与善性xxx| 床上黄色一级片| 精品久久久久久久久久久久久| 国产三级在线视频| 日韩亚洲欧美综合| av福利片在线观看| 最近视频中文字幕2019在线8| 亚洲第一电影网av| 亚洲乱码一区二区免费版| 最近手机中文字幕大全| 久久久久久久久久成人| 亚洲精品粉嫩美女一区| av黄色大香蕉| 亚洲国产精品sss在线观看| 熟女人妻精品中文字幕| av又黄又爽大尺度在线免费看 | 日韩制服骚丝袜av| 国产精品一二三区在线看| 我的老师免费观看完整版| 欧美激情在线99| 久久久久免费精品人妻一区二区| 一边亲一边摸免费视频| 在线观看免费视频日本深夜| 亚洲内射少妇av| 天堂影院成人在线观看| 少妇人妻一区二区三区视频| 亚洲欧美中文字幕日韩二区| 网址你懂的国产日韩在线| 免费看美女性在线毛片视频| 成人欧美大片| 日韩欧美三级三区| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久久电影| 一个人免费在线观看电影| 国产av麻豆久久久久久久| 国产精品久久久久久精品电影| 美女高潮的动态| 亚洲精品久久久久久婷婷小说 | 国产一区亚洲一区在线观看| 别揉我奶头 嗯啊视频| 婷婷色综合大香蕉| 欧美色视频一区免费| 国产精品一区二区在线观看99 | 精品久久久久久久久久久久久| 国产一区二区在线观看日韩| 国产亚洲欧美98| 久久欧美精品欧美久久欧美| 能在线免费观看的黄片| 国产极品精品免费视频能看的| 国产高清有码在线观看视频| 久久久久久久久中文| 亚洲自偷自拍三级| 久久精品国产鲁丝片午夜精品| 五月伊人婷婷丁香| 亚洲人与动物交配视频| 国产精品一区二区性色av| 免费观看的影片在线观看| 岛国毛片在线播放| 超碰av人人做人人爽久久| 小蜜桃在线观看免费完整版高清| АⅤ资源中文在线天堂| 日本一本二区三区精品| 26uuu在线亚洲综合色| 久久久国产成人精品二区| 男女啪啪激烈高潮av片| 日本欧美国产在线视频| 一级黄色大片毛片| 卡戴珊不雅视频在线播放| 亚洲欧美日韩东京热| 一级毛片电影观看 | 蜜桃亚洲精品一区二区三区| 亚洲欧美精品专区久久| 国产亚洲欧美98| av天堂在线播放| 人妻系列 视频| 日产精品乱码卡一卡2卡三| 日本撒尿小便嘘嘘汇集6| 亚洲精品久久久久久婷婷小说 | 青青草视频在线视频观看| 日日摸夜夜添夜夜添av毛片| 亚洲最大成人av| 97超碰精品成人国产| 淫秽高清视频在线观看| 国产在线男女| 一级毛片电影观看 | 亚洲成a人片在线一区二区| 亚洲激情五月婷婷啪啪| 午夜爱爱视频在线播放| 12—13女人毛片做爰片一| 一本一本综合久久| 国产一级毛片七仙女欲春2| 禁无遮挡网站| 精品熟女少妇av免费看| 亚洲欧美清纯卡通| 九九在线视频观看精品| 久久久久九九精品影院| 小说图片视频综合网站| 能在线免费看毛片的网站| 欧美区成人在线视频| 亚洲五月天丁香| 国产精品野战在线观看| 精品午夜福利在线看| 精品久久久久久久久久免费视频| 日韩欧美三级三区| 亚洲在线自拍视频| 国产男人的电影天堂91| 免费观看精品视频网站| 国产国拍精品亚洲av在线观看| 成人综合一区亚洲| 国国产精品蜜臀av免费| 国产激情偷乱视频一区二区| 亚洲国产欧美人成| 免费一级毛片在线播放高清视频| 色综合站精品国产| 91精品一卡2卡3卡4卡| 18禁在线播放成人免费| 国产伦理片在线播放av一区 | 蜜桃亚洲精品一区二区三区| 欧美bdsm另类| 国产乱人偷精品视频| 国产精品免费一区二区三区在线| av在线老鸭窝| 亚洲熟妇中文字幕五十中出| 久久这里有精品视频免费| 狂野欧美激情性xxxx在线观看| 久久久久久久久中文| 乱人视频在线观看| 国产精品99久久久久久久久| 国产黄片视频在线免费观看| 大香蕉久久网| 一个人免费在线观看电影| 天天一区二区日本电影三级| 日日摸夜夜添夜夜爱| 白带黄色成豆腐渣| 性色avwww在线观看| 永久网站在线| 久久久国产成人免费| 又粗又爽又猛毛片免费看| 人妻系列 视频| 黄片无遮挡物在线观看| 国产精品蜜桃在线观看 | 亚洲,欧美,日韩| 久久6这里有精品| 国产精品一区二区在线观看99 | 欧美高清成人免费视频www| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 久久韩国三级中文字幕| 女同久久另类99精品国产91| 麻豆成人午夜福利视频| 波多野结衣高清无吗| 91精品国产九色| 久久亚洲精品不卡| 久久久久久久亚洲中文字幕| 99在线人妻在线中文字幕| 国产在线男女| 天美传媒精品一区二区| 婷婷精品国产亚洲av| 日本欧美国产在线视频| 亚洲精品乱码久久久久久按摩| 婷婷色av中文字幕| 精品一区二区免费观看| 又黄又爽又刺激的免费视频.| 男女边吃奶边做爰视频| 国产成人一区二区在线| 91久久精品国产一区二区成人| 中文亚洲av片在线观看爽| 亚洲成人久久爱视频| 精品人妻视频免费看| 久久精品国产鲁丝片午夜精品| 国产高清视频在线观看网站| 国产av麻豆久久久久久久| 欧美性感艳星| 国产精品.久久久| 久久国内精品自在自线图片| 日韩精品青青久久久久久| 91精品国产九色| 我的女老师完整版在线观看|