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

    中國(guó)夏季降水異常EOF模態(tài)的時(shí)間穩(wěn)定性分析

    2014-12-14 09:13:48龐軼舒祝從文劉凱
    大氣科學(xué) 2014年6期
    關(guān)鍵詞:距平百分率降水

    龐軼舒 祝從文 劉凱

    中國(guó)氣象科學(xué)研究院,北京100081

    1 引言

    經(jīng)驗(yàn)正交分解(Empirical orthogonal function,EOF)又被稱為主成分分析(Principal component analysis,PCA),最早由Pearson提出(Pearson,1902),20世紀(jì)50年代Lorenz將該方法引入大氣科學(xué)領(lǐng)域(Lorenz,1956),隨后被廣泛應(yīng)用至今。EOF方法的優(yōu)點(diǎn)在于它能對(duì)分布不規(guī)則的氣象要素場(chǎng)進(jìn)行時(shí)空分解,得到的各特征向量(空間模態(tài)和時(shí)間系數(shù))相互正交,其中空間模態(tài)(EOF)在一定程度上可反映要素場(chǎng)的空間分布特點(diǎn),而時(shí)間系數(shù)(PC)則反映相應(yīng)空間模態(tài)隨時(shí)間的權(quán)重變化。除此之外,EOF方法還具有展開收斂速度快,容易將原始要素場(chǎng)的變化信息濃縮在前幾個(gè)模態(tài)上的特點(diǎn)。因此,該方法常用于氣象要素場(chǎng)時(shí)空變化特征規(guī)律研究。

    EOF方法是氣候研究中的主要方法之一,在過去幾十年,中國(guó)學(xué)者利用該方法對(duì)中國(guó)降水時(shí)空變化特征做了大量研究(例如,王紹武和趙宗慈,1979;鄧愛軍等,1989;李棟梁等,1997;宋正山和楊輝,2001;魏鳳英,2003;孫林海等,2005;Wang et al.,2013;封國(guó)林等,2013 ;Zhang et al., 2013)。其中,王紹武和趙宗慈(1979)利用此方法分析了我國(guó)旱澇分布特點(diǎn),并在此基礎(chǔ)上提出了中國(guó)6種主要雨型的觀點(diǎn)。宋正山和楊輝(2001)利用此方法分析了500 hPa環(huán)流場(chǎng)與我國(guó)夏季降水(雨型)相關(guān)分布,指出該環(huán)流場(chǎng)前兩個(gè)模態(tài)與中國(guó)夏季降水主要分布型有很好的對(duì)應(yīng)關(guān)系。除此之外,學(xué)者們還將EOF方法應(yīng)用于短期氣候和區(qū)域統(tǒng)計(jì)降尺度預(yù)測(cè)中,通過EOF各模態(tài)的PC系數(shù)預(yù)測(cè)和重建來預(yù)測(cè)未來要素場(chǎng)的變化(章基嘉等,1979;Harnack and Lanzante, 1985;袁景鳳和吳曉曦,2000;何慧等,2004)。例如,何慧等(2004)通過對(duì)廣西前汛期降水距平百分率作 EOF 分解,選取累積方差貢獻(xiàn)接近70%的前五個(gè)模態(tài)為預(yù)報(bào)量,利用統(tǒng)計(jì)方法建立降水預(yù)報(bào)模型。眾所周知,在EOF重建過程中,參與模態(tài)越多,重建場(chǎng)與原始場(chǎng)越接近。然而,在利用EOF方法預(yù)測(cè)時(shí),參與模態(tài)的選取則有嚴(yán)格的條件限制。章基嘉等(1981)指出,只有EOF模態(tài)在預(yù)測(cè)時(shí)效內(nèi)穩(wěn)定,該模態(tài)的信息才可用于未來的氣候預(yù)測(cè)。Charles et al.(1999)、Wilby and Wigley(2000)指出預(yù)報(bào)因子和預(yù)測(cè)函數(shù)隨時(shí)間穩(wěn)定是建立統(tǒng)計(jì)預(yù)測(cè)模型的基本條件。因此,在利用 EOF方法建立氣候預(yù)測(cè)模型時(shí),需首先討論各EOF模態(tài)是否隨時(shí)間發(fā)生變化。為方便討論,本文將 EOF模態(tài)不隨時(shí)間發(fā)生變化稱作時(shí)間穩(wěn)定。

    夏季降水一直是我國(guó)短期氣候預(yù)測(cè)的重點(diǎn)內(nèi)容,近年來許多氣象工作者在前人研究的基礎(chǔ)上提出了新的預(yù)測(cè)方法及思路并取得了不錯(cuò)的預(yù)測(cè)效果(郭玲,2011;郭彥和李建平,2012;Liu and Fan,2012),但就現(xiàn)階段而言,全國(guó)范圍的降水預(yù)測(cè)水平還有待進(jìn)一步提高。而研究中國(guó)夏季降水 EOF模態(tài)的時(shí)間穩(wěn)定性對(duì)于認(rèn)識(shí)中國(guó)夏季降水主要變化特征,并從降水全局變化角度探索可能的預(yù)測(cè)方法具有重要的意義。因此,本文分析和討論了1980~2012年中國(guó)夏季降水距平和距平百分率前四個(gè)EOF模態(tài)的特征和時(shí)間穩(wěn)定性,在此基礎(chǔ)上探討 EOF方法在整個(gè)中國(guó)地區(qū)夏季降水預(yù)測(cè)中的潛在能力和應(yīng)用條件。

    2 資料和分析方法

    本文采用的資料為中國(guó)氣象局信息中心提供的1951~2012年全國(guó)160臺(tái)站夏季(6~8月)逐月降水資料,該套資料常用于國(guó)家氣候中心的短期氣候預(yù)測(cè)工作。研究過程中,本文首先對(duì) 1980~2012年降水距平和距平百分率進(jìn)行EOF分解,取前四個(gè)模態(tài)進(jìn)行分析。在此基礎(chǔ)上,利用滑動(dòng)交叉檢驗(yàn)、獨(dú)立樣本檢驗(yàn)、距平相關(guān)系數(shù)(Anomaly correlation coefficient,ACC)、時(shí)間相關(guān)系數(shù)(Temporal correlation coefficient,TCC)等方法討論EOF各模態(tài)的時(shí)間穩(wěn)定性。為定量評(píng)估EOF各模態(tài)的時(shí)間穩(wěn)定性,本文定義了時(shí)間穩(wěn)定性指數(shù)(Stability Index, SI),即資料時(shí)間長(zhǎng)度處理前后對(duì)應(yīng)空間模態(tài)ACC和時(shí)間系數(shù)TCC的乘積,并將穩(wěn)定標(biāo)準(zhǔn)設(shè)為 0.5。其中,資料時(shí)間長(zhǎng)度處理分別通過滑動(dòng)交叉檢驗(yàn)[見公式(1)]和獨(dú)立樣本檢驗(yàn)[見公式(2)]兩種方法實(shí)現(xiàn)。

    其中,Ri,t為時(shí)間長(zhǎng)度處理后的資料(i=1980,1981,…,2012;t為剔除時(shí)間長(zhǎng)度:1,2,3,4,5);R0為原始資料(1980~2012年夏季降水場(chǎng));r為樣本元素。

    其中,Ri為時(shí)間長(zhǎng)度處理后的資料(i=2003,2004,…,2012);R0為原始資料(1980~2002年夏季降水場(chǎng));r為樣本元素。

    3 中國(guó)夏季降水異常主要模態(tài)的特征分析

    過去很多研究認(rèn)為中國(guó)東西部降水氣候特點(diǎn)不同(徐國(guó)昌,1997;孫國(guó)武,1997;朱乾根等,2007),基于這種觀點(diǎn),對(duì)整個(gè)中國(guó)地區(qū)夏季降水進(jìn)行EOF分析很可能無法區(qū)分東西部降水的特征。為此,我們對(duì)比了 1980~2012年中國(guó)和中國(guó)東部地區(qū)(100°E以東)夏季降水距平、距平百分率的前四個(gè)EOF模態(tài)。對(duì)比結(jié)果表明,兩者夏季降水距平各模態(tài)的解釋方差差異較小,對(duì)應(yīng)空間分布形態(tài)相似,PC系數(shù)變化一致,相關(guān)系數(shù)接近1.0(圖略)。而兩者降水距平百分率各模態(tài)的差異則主要表現(xiàn)在前兩個(gè)模態(tài)的順序發(fā)生了對(duì)調(diào),其對(duì)應(yīng)PC系數(shù)的時(shí)間相關(guān)關(guān)系仍然通過顯著性水平為 0.01的檢驗(yàn)。由此可見,盡管中國(guó)東西部地區(qū)降水特性有所差異,但由于EOF分析結(jié)果突出的是降水方差貢獻(xiàn)較大的區(qū)域,所以從整體角度進(jìn)行EOF分析不僅能很好地反映方差貢獻(xiàn)較大的東部降水特征,同時(shí)也能反映西部地區(qū)降水的一些特性,適用于全國(guó)夏季降水預(yù)測(cè)。

    基于上述討論,本文重點(diǎn)分析中國(guó)夏季降水距平和距平百分率的前四個(gè)EOF模態(tài)。圖1表示降水距平前四個(gè)EOF模態(tài)的空間分布特征,其中,各模態(tài)的解釋方差分別為 18.45%、10.93%、8.14%和7.44%,累計(jì)方差貢獻(xiàn)約 45%。如圖所示,第一模態(tài)在100°E以東地區(qū)表現(xiàn)出南北向“+、-、+”的三極型分布,是東亞夏季風(fēng)降水的主要空間分布型;第二模態(tài)表現(xiàn)出以長(zhǎng)江為界,南北方降水反相變化的特征;第三模態(tài)表現(xiàn)為江南地區(qū)、東北部分地區(qū)和西北地區(qū)降水偏少,其余地區(qū)降水偏多的特征;第四模態(tài)則呈現(xiàn)華南、黃淮區(qū)域降水偏少,其他區(qū)域降水偏多的特征。而從降水距平百分率前四個(gè)EOF模態(tài)分析結(jié)果(圖2)可以看出,其各模態(tài)的解釋方差依次是 14.46%、10.64%、7.95%和6.39%,累計(jì)方差貢獻(xiàn)為39.44%。其中第一模態(tài)的空間分布表現(xiàn)為黃河以北地區(qū)與其余地區(qū)降水反向變化的特征;第二模態(tài)則呈現(xiàn)出南北向“-、+、-”三極型分布特征;第三模態(tài)表現(xiàn)為河套地區(qū)、華南沿海地區(qū)和其余地區(qū)降水反向變化特征;第四模態(tài)則表現(xiàn)為全國(guó)大部分地區(qū)降水偏多,零星地區(qū)降水偏少的特征。

    從上述分析可以發(fā)現(xiàn),作為指示降水異常的兩種物理量,降水距平和距平百分率前四個(gè)模態(tài)的空間分布表現(xiàn)出多樣性和復(fù)雜性,且各模態(tài)所占的方差貢獻(xiàn)率最大不超過20%,這說明中國(guó)夏季降水的年際變化特征多樣而復(fù)雜,單個(gè)EOF模態(tài)反映的信息十分有限,因此在利用EOF方法預(yù)測(cè)中國(guó)夏季降水時(shí),必須盡可能考慮多個(gè)模態(tài)的信息。另外,由于降水距平和距平百分率側(cè)重不同,前者反映降水偏離氣候值的大小,后者反映降水偏離氣候值的程度,所以兩者各模態(tài)的空間分布差異明顯。這說明在利用 EOF方法預(yù)測(cè)中國(guó)夏季降水時(shí),這兩種物理量不能相互替代。

    4 中國(guó)夏季降水異常主要模態(tài)的時(shí)間穩(wěn)定性分析

    如前所述,時(shí)間穩(wěn)定是各EOF模態(tài)應(yīng)用到氣候預(yù)測(cè)中的重要前提條件。為此,本文首先通過滑動(dòng)交叉檢驗(yàn)方法討論中國(guó)夏季降水距平和距平百分率前四個(gè)模態(tài)的時(shí)間穩(wěn)定性,以探究各EOF模態(tài)是否可參與中國(guó)夏季降水預(yù)測(cè)。

    圖3表示通過滑動(dòng)交叉檢驗(yàn)得到的中國(guó)夏季降水距平各EOF模態(tài)SI隨時(shí)間和剔除時(shí)間長(zhǎng)度(t)的變化情況,如圖所示,當(dāng)t≤5時(shí),第一個(gè)模態(tài)的SI恒大于 0.9,高度穩(wěn)定。第二模態(tài)在t=1時(shí),SI恒大于0.8,高度穩(wěn)定;在t=2時(shí),SI在個(gè)別年份(如1998年)略低于穩(wěn)定度標(biāo)準(zhǔn)0.5,其余年份高于0.8,基本穩(wěn)定;在t>2時(shí),SI的時(shí)間波動(dòng)幅度增大,低于穩(wěn)定標(biāo)準(zhǔn)的年份增多,時(shí)間穩(wěn)定性下降。第三、四模態(tài)在t=1時(shí),除個(gè)別年份(2008年)外,SI總體高于穩(wěn)定性標(biāo)準(zhǔn),基本穩(wěn)定;在t>1時(shí),SI波動(dòng)幅度隨t增加而增大,數(shù)值明顯下降。由此可見,在利用EOF方法預(yù)測(cè)中國(guó)夏季降水距平時(shí),隨著預(yù)測(cè)時(shí)效增加,可參與預(yù)測(cè)的模態(tài)減少。當(dāng)預(yù)測(cè)時(shí)效為一年時(shí),前四個(gè)模態(tài)穩(wěn)定,可參與未來氣候預(yù)測(cè);當(dāng)預(yù)測(cè)時(shí)效為兩年時(shí),前兩個(gè)模態(tài)基本穩(wěn)定,可參與預(yù)測(cè);當(dāng)預(yù)測(cè)時(shí)效大于兩年,僅第一模態(tài)的穩(wěn)定性滿足參與預(yù)測(cè)條件。

    與降水距平相比,2012年在降水距平百分率各模態(tài) SI隨時(shí)間和剔除時(shí)間長(zhǎng)度(t)的變化中表現(xiàn)出顯著異常(見圖4)。當(dāng)被剔除資料包含2012年時(shí),各模態(tài) SI急劇下降。如果忽略 2012年的影響,則前兩個(gè)EOF模態(tài)的SI始終高于0.5;第三模態(tài)的SI在t≤3時(shí)高于0.5,在t=4、5時(shí)除1994年

    略低于0.5,其余高于0.5。而第四模態(tài)的SI波動(dòng)較大,低于0.5的年份較多,穩(wěn)定性較差。由此可見,在利用 EOF方法預(yù)測(cè)中國(guó)夏季降水距平百分率時(shí),第四模態(tài)的穩(wěn)定性較低,無法參與預(yù)測(cè)。另外,基于2012年之前的EOF模態(tài)特征無法準(zhǔn)確預(yù)測(cè)該年的降水距平百分率。

    圖1 1980~2012年中國(guó)夏季降水距平場(chǎng)前四個(gè)EOF模態(tài)空間分布:(a)EOF-1;(b)EOF-2;(c)EOF-3;(d)EOF-4Fig. 1 The first four EOF modes of China summer precipitation anomalies during 1980–2012: (a) EOF-1; (b) EOF-2; (c) EOF-3; (d) EOF-4

    圖2 同圖1,但表示的是距平百分率場(chǎng)Fig. 2 Same as Fig. 1, but for the field of precipitation anomaly percentage

    圖3 1980~2012年中國(guó)夏季降水距平場(chǎng)前四個(gè)EOF模態(tài)穩(wěn)定度(SI)隨時(shí)間和剔除時(shí)間長(zhǎng)度的變化:(a)EOF-1;(b)EOF-2;(c)EOF-3;(d)EOF-4Fig. 3 Stability indices of the first four EOF modes of China summer precipitation anomaly during 1980–2012, and their change with time and length of excluded time: (a) EOF-1; (b) EOF-2; (c) EOF-3; (d) EOF-4

    圖4 同圖3,但表示的是距平百分率場(chǎng)Fig. 4 Same as Fig. 3, but for the field of precipitation anomaly percentage

    雖然滑動(dòng)交叉檢驗(yàn)結(jié)果在一定程度上反映了各模態(tài)的時(shí)間穩(wěn)定性,但在實(shí)際預(yù)測(cè)時(shí)不知道未來的情況,預(yù)測(cè)結(jié)果仍然存在較大的不確定性。為此,我們通過獨(dú)立樣本檢驗(yàn)方法進(jìn)一步檢驗(yàn)了降水距平和距平百分率前四個(gè)EOF模態(tài)在近10年的時(shí)間穩(wěn)定性。結(jié)果表明,近 10年內(nèi)中國(guó)夏季降水距平前三個(gè)模態(tài)的SI均大于0.9,表現(xiàn)出高度的穩(wěn)定性;第四模態(tài)的SI在近10年內(nèi)波動(dòng)相對(duì)較大,盡管如此,其SI僅有一年(2007年)低于穩(wěn)定度標(biāo)準(zhǔn)(圖略)。而中國(guó)夏季降水距平百分率第一、四模態(tài)的SI在近10年中均大于或等于0.9,第二、三模態(tài)的SI除在2012年低于0.5,其余均大于0.8??傮w而言,中國(guó)夏季降水距平和距平百分率前四個(gè)模態(tài)在近 10年中基本穩(wěn)定。然而與滑動(dòng)交叉檢驗(yàn)結(jié)果類似,降水距平百分率前四個(gè)模態(tài)的時(shí)間穩(wěn)定性在2012年相對(duì)較低。

    5 極端降水事件對(duì) EOF模態(tài)時(shí)間穩(wěn)定性的影響

    由于 EOF特征向量取決于樣本資料的協(xié)方差矩陣,因此某一年的極端降水氣候事件可能會(huì)對(duì)多年降水(距平和距平百分率)的EOF模態(tài)時(shí)間穩(wěn)定性產(chǎn)生影響。2012年中國(guó)北方一些地區(qū)出現(xiàn)極端降水事件,以北京“7.21”大暴雨為例,全市日平均降水量達(dá) 170 mm,其中房山區(qū)高達(dá) 460 mm(http://www.cma.gov.cn/2011xwzx/2011xmtjj/201207/t20120724_179464.html)。因此,降水距平百分率各模態(tài)時(shí)間穩(wěn)定性在 2012年的特殊性很可能由極端降水引起。

    為驗(yàn)證極端降水氣候事件對(duì)中國(guó)夏季降水EOF模態(tài)穩(wěn)定性的影響,本文將1951~2012年中國(guó)160臺(tái)站夏季降水(R)進(jìn)行閾值修訂(見公式3)以削弱極端降水的強(qiáng)度。經(jīng)檢驗(yàn),該修訂方法不會(huì)破壞原始降水的變化趨勢(shì)和其降水距平(距平百分率)各 EOF模態(tài)的基本特征(圖略)。在此基礎(chǔ)上,對(duì)比閾值修訂前后降水距平和距平百分率前四個(gè)模態(tài)SI的滑動(dòng)交叉檢驗(yàn)結(jié)果。

    其中,Rm為修訂后的降水資料,RMin和RMax則被定義為1951~2012年各站點(diǎn)降水序列(由低到高)5%和95%分位數(shù)的值。

    通過對(duì)比修訂前后降水距平各模態(tài)的 SI(圖5)可以發(fā)現(xiàn),當(dāng)t≤5時(shí),前四個(gè)模態(tài)的 SI在大多數(shù)年份中減小,在少數(shù)年份中增加,差異值小于 0.3,改變幅度不大。但比較特殊的是,當(dāng)t=2時(shí),第三、四模態(tài)的SI在一些年份中下降程度大于0.5。

    與降水距平不同,降水距平百分率第一、二、四模態(tài)的SI在閾值修訂前后總體提高(圖6),尤其當(dāng)剔除資料中含2012年時(shí)第一模態(tài)(t≤5)和第二模態(tài)(t≤2)SI的增加值超過0.5。而其第三模態(tài)SI的下降年份數(shù)偏多,但改變幅度不大。閾值修訂后,距平百分率各模態(tài) SI隨時(shí)間和剔除時(shí)間長(zhǎng)度(t)的變化如圖7所示,當(dāng)t≤5時(shí),其第一模態(tài)的SI恒大于0.8,高度穩(wěn)定。而第二、三模態(tài)的SI在t=1時(shí)大于0.6,基本穩(wěn)定;在t>1時(shí),低于0.5的年份數(shù)增多,穩(wěn)定性較差。而第四模態(tài)SI在t≤5時(shí)波動(dòng)明顯,低于0.5的年份較多,穩(wěn)定性較差。

    從以上分析可以看出,閾值修訂方法對(duì)降水距平百分率各模態(tài) SI有明顯改善作用,消除了其在2012年的特殊性,而對(duì)降水距平各模態(tài)SI的改善作用則不明顯。這說明,與降水距平相比,降水距平百分率各模態(tài)的時(shí)間穩(wěn)定性更易受到極端降水的影響。然而極端降水的產(chǎn)生不僅受到氣候變化的影響,天氣尺度和中小尺度氣象條件變化的影響也十分重要,但在短期氣候預(yù)測(cè)中,天氣尺度和中小尺度的氣象要素變化無法準(zhǔn)確預(yù)測(cè)。因此在利用EOF方法預(yù)測(cè)中國(guó)夏季降水距平百分率時(shí),可首先對(duì)降水資料進(jìn)行閾值修訂。修訂后,當(dāng)預(yù)測(cè)時(shí)效為1年時(shí),前三個(gè)模態(tài)可參加預(yù)測(cè);當(dāng)預(yù)測(cè)時(shí)效為 1年以上時(shí),僅第一模態(tài)可參加預(yù)測(cè)。

    綜合第4、5節(jié)研究結(jié)果可以發(fā)現(xiàn),當(dāng)利用EOF方法預(yù)測(cè)未來一年中國(guó)夏季降水時(shí),可參與預(yù)測(cè)的模態(tài)最多,即降水距平前四個(gè)模態(tài)或修訂后降水距平百分率前三個(gè)模態(tài)。假設(shè)各模態(tài)PC系數(shù)預(yù)測(cè)完全準(zhǔn)確,則可以通過它們的EOF重建場(chǎng)與原始場(chǎng)的TCC和ACC來討論EOF方法的潛在預(yù)測(cè)能力。研究結(jié)果表明(圖略),降水距平前四個(gè)模態(tài)能解釋原始降水 45%的變化方差,其重建場(chǎng)與原始場(chǎng)的TCC顯著相關(guān)(α=0.05)站點(diǎn)集中分布在黃河以南大部分地區(qū),平均ACC約為0.6。而修訂后降水距平百分率前三個(gè)模態(tài)能解釋原始降水 30.74%的變化方差,其重建場(chǎng)與原始場(chǎng)的 TCC顯著相關(guān)(α=0.05)站點(diǎn)分布在中國(guó)大部分地區(qū),平均 ACC為0.48。由此可見,EOF方法在中國(guó)夏季降水年際變化預(yù)測(cè)中有一定潛力,當(dāng)以降水距平百分率為預(yù)測(cè)對(duì)象時(shí),可預(yù)測(cè)范圍相對(duì)較大,當(dāng)以降水距平為預(yù)測(cè)對(duì)象時(shí),預(yù)測(cè)結(jié)果相對(duì)準(zhǔn)確。

    圖5 1980~2012年閾值修訂前后中國(guó)夏季降水距平場(chǎng)EOF模態(tài)的穩(wěn)定度(SI)差異(后者減前者):(a)EOF-1;(b)EOF-2;(c)EOF-3;(d)EOF-4Fig. 5 Differences of stability indices between the original and modified Chinasummer precipitation anomalies for the first four EOF modes during 1980–2012: (a) EOF-1; (b) EOF-2; (c) EOF-3; (d) EOF-4

    圖6 同圖5,表示的是距平百分率場(chǎng)Fig. 6 Same as Fig. 5, but for the precipitation anomaly percentage

    圖7 閾值修訂后1980~2012年中國(guó)夏季降水距平百分率場(chǎng)前四個(gè)EOF模態(tài)穩(wěn)定度(SI)隨時(shí)間和剔除時(shí)間長(zhǎng)度的變化:(a)EOF-1;(b)EOF-2;(c)EOF-3;(d)EOF-4Fig. 7 Stability indices of the first four EOF modes of the modified China summer precipitation anomaly percentage during 1980–2012, and their change with time and length of excluded time: (a) EOF-1; (b) EOF-2; (c) EOF-3; (d) EOF-4

    6 結(jié)論和討論

    本文基于1980~2012年中國(guó)160個(gè)臺(tái)站夏季降水資料,分析和研究了降水距平與距平百分率前四個(gè)EOF模態(tài)的特征和時(shí)間穩(wěn)定性,討論了極端降水事件對(duì)各模態(tài)時(shí)間穩(wěn)定性的影響,并提出了一種通過削弱極端降水強(qiáng)度來提高降水距平百分率各模態(tài)時(shí)間穩(wěn)定性的方法。在此基礎(chǔ)上,探討了EOF方法在中國(guó)夏季降水短期氣候預(yù)測(cè)中的潛在能力。

    研究發(fā)現(xiàn),中國(guó)夏季降水距平第一模態(tài)在預(yù)測(cè)時(shí)效為5年內(nèi)時(shí)間穩(wěn)定,第二模態(tài)在2年內(nèi)時(shí)間穩(wěn)定,第三、四模態(tài)在1年內(nèi)時(shí)間穩(wěn)定。而受極端降水事件影響,中國(guó)夏季降水距平百分率各模態(tài)時(shí)間穩(wěn)定性較低,無法參與預(yù)測(cè)。當(dāng)對(duì)資料進(jìn)行閾值修訂,人為削弱極端降水影響后,距平百分率前三個(gè)模態(tài)的時(shí)間穩(wěn)定性顯著提高,其中第一模態(tài)在預(yù)測(cè)時(shí)效為5年內(nèi)時(shí)間穩(wěn)定,第二、三模態(tài)在1年內(nèi)時(shí)間穩(wěn)定。因此,可以利用降水距平前四個(gè)EOF模態(tài)或修訂后降水距平百分率前三個(gè) EOF模態(tài)的信息預(yù)測(cè)未來一年中國(guó)夏季降水的變化。假設(shè)PC系數(shù)預(yù)測(cè)完全準(zhǔn)確,降水距平前四個(gè)模態(tài)可以解釋原始降水 45%的變化方差,其重建場(chǎng)與原始場(chǎng)的 TCC顯著相關(guān)站點(diǎn)集中分布在黃河以南地區(qū),平均ACC約為 0.6。而修訂后降水距平百分率前三個(gè)模態(tài)能解釋原始降水 30.74%的變化方差,其重建場(chǎng)與原始場(chǎng)的 TCC顯著相關(guān)站點(diǎn)在全國(guó)均勻分布,平均ACC為0.48。相較而言,在利用EOF方法預(yù)測(cè)中國(guó)夏季降水時(shí),距平比距平百分率更適合作預(yù)測(cè)對(duì)象。

    本文主要討論了 EOF方法在中國(guó)夏季降水短期氣候預(yù)測(cè)中的潛在應(yīng)用價(jià)值。但在實(shí)際預(yù)測(cè)之前,我們還需要進(jìn)一步探索各模態(tài)PC系數(shù)的影響因子以及多因子配置與各模態(tài)PC之間的關(guān)系,以建立中國(guó)夏季降水預(yù)測(cè)模型。此外,新研究表明在降水預(yù)測(cè)中以年際增量作為預(yù)測(cè)對(duì)象能取得不錯(cuò)的預(yù)測(cè)結(jié)果(范可等,2007),因而通過年際增量討論 EOF方法的潛在預(yù)測(cè)能力亦十分重要。這兩個(gè)問題將是作者今后研究的重要內(nèi)容。

    References)

    Charles S P, Bates B C, Whetton P H, et al. 1999. Validation of downscaling models for changed climate conditions: Case study of southwestern Australia [J]. Climate Research, 12: 1–14.

    鄧愛軍, 陶詩言, 陳烈庭. 1989. 我國(guó)汛期降水的EOF分析 [J]. 大氣科學(xué), 13 (3): 289–294. Deng Aijun, Tao Shiyan, Chen Lieting. 1989. The EOF analysis of rainfall in China during monsoon season [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 13 (3): 289–294.

    范可, 王會(huì)軍, Young-Jean C. 2007. 一個(gè)長(zhǎng)江中下游夏季降水的物理統(tǒng)計(jì)預(yù)測(cè)模型 [J]. 科學(xué)通報(bào), 52 (24): 2900–2905. Fan Ke, Wang Huijun, Young-Jean C. 2007. A physical-statistical prediction model for summer rainfall in the middle and lower reaches of the Yangtze River [J].Chinese Science Bulletin (in Chinese), 52 (24): 2900–2905.

    封國(guó)林, 孫樹鵬, 趙俊虎, 等. 2013. 基于 2009年初長(zhǎng)江中下游地區(qū)持續(xù)陰雨過程的10~30天延伸期穩(wěn)定分量的提取及配置分析 [J]. 中國(guó)科學(xué): 地球科學(xué), 43 (5): 836–847. Feng Guolin, Sun Shupeng, Zhao Junhu, et al. 2013. Analysis of stable components for extended-range(10–30 days) weather forecast: A case study of continuous overcast-rainy process in early 2009 over the mid-lower reaches of the Yangtze River [J].Science China: Earth Sciences. 56 (9): 1576–1587.

    郭玲. 2011. 中國(guó)夏季降水影響因子評(píng)估和統(tǒng)計(jì)預(yù)測(cè)研究 [D]. 南京信息工程大學(xué)碩士學(xué)位論文. Guo Ling. 2011. Potential predictors evaluations and statistical prediction of summer rainfall in China [D]. M.S. thesis (in Chinese), Nanjing University of Information Science and Technology.

    郭彥, 李建平. 2012. 一種分離時(shí)間尺度的統(tǒng)計(jì)降尺度模型的建立和應(yīng)用——以華北汛期降水為例 [J]. 大氣科學(xué), 36 (2): 385–396. Guo Yan, Li Jianping. 2012. A time-scale decomposition statistical downscaling model: Case study of North China rainfall in rainy season[J]. Chinese Journal of Atmospheric Sciences (in Chinese), 36 (2): 385–396.

    Harnack R P, Lanzante J R. 1985. Specification of United States seasonal precipitation [J]. Mon. Wea. Rev., 113: 319–325.

    李棟梁, 謝金南, 王文. 1997. 中國(guó)西北夏季降水特征及其異常研究 [J].大氣科學(xué), 21 (3): 331–340. Li Dongliang, Xie Jinnan, Wang Wen.1997. A study of summer precipitation features and anomaly in northwest China [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 21 (3):331–340.

    Lorenz E N. 1956. Empirical orthogonal functions and statistical weather prediction [R]. Science Report 1, Statistical Forecasting Project,Department of Meteorology, MIT (NTIS AD 110268), 49 PP.

    Liu Y, Fan K. 2012. A new statistical downscaling model for autumn precipitation in China [J]. International Journal of Climatology, 33 (6):1321–1336.

    何慧, 金龍, 覃志年. 2004. 降維方法在廣西前汛期降水預(yù)測(cè)中的應(yīng)用[J]. 災(zāi)害學(xué), 19 (4): 7–12. He Hui, Jin Long, Qin Zhinian. 2004. The application of EOF method to the precipitation forecast during the first raining season in Guangxi [J]. Journal of Catastrophology (in Chinese),19 (4): 7–12.

    Pearson K. 1902. On lines and planes of closest fit to systems of points in space [J]. Philos. Mag., 2 (6): 559–572.

    宋正山, 楊輝. 2001. 夏季東亞季風(fēng)區(qū)500 hPa月環(huán)流異常及與我國(guó)降水關(guān)系的向量 EOF分析 [J]. 大氣科學(xué), 25 (3): 401–410. Song Zhengshan, Yang Hui. 2001. Vector EOF expansion of summer monthly anomalous wind fields at 500 hPa over East Asia relation to rainfall over China [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 25 (3):401–410.

    孫國(guó)武. 1997. 中國(guó)西北干旱氣候研究 [M]. 北京: 氣象出版社, 384.Sun Guowu. 1997. Research of Arid Climate in Northwest China (in Chinese) [M]. Beijing: China Meteorological Press, 384.

    孫林海, 趙振國(guó), 許力, 等. 2005. 中國(guó)東部季風(fēng)區(qū)夏季雨型的劃分及其環(huán)流成因分析 [J]. 應(yīng)用氣象學(xué)報(bào), 16 (增刊): 56–62. Sun Linhai,Zhao Zhenguo, Xu Li, et al. 2005. Study of summer rain pattern in monsoon region of East China and its circulation cause [J]. Journal of Applied Meteorological Science (in Chinese), 16 (S1): 56–62.

    王紹武, 趙宗慈. 1979. 近五百年我國(guó)旱澇史料的分析 [J]. 地理學(xué)報(bào),34 (4): 329–341. Wang Shaowu, Zhao Zongci. 1979. An analyses of historical data of droughts and floods in last 500 years in China [J]. Acta Geographica Sinica (in Chinese), 34 (4): 329–341.

    魏鳳英. 2003. 現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù) [M]. 北京: 氣象出版社,105–113. Wei Fengying. 2003. Statistical Diagnosis and Prediction Techniques of Modern Climate (in Chinese) [M]. Beijing: China Meteorological Press, 105–113.

    Wang B, Xiang B Q, Lee J Y. 2013. Subtropical High predictability establishes a promising way for monsoon and tropical storm predictions[J]. Proc. Natl. Acad. Sci. USA, 110: 2718–2722.

    Wilby R L, Wigley T M L. 2000. Precipitation predictors for downscaling:Observed and general circulation model relationships [J]. International Journal of Climatology, 20: 641–661.

    徐國(guó)昌. 1997. 中國(guó)干旱半干旱區(qū)氣候變化 [M]. 北京: 氣象出版社,101. Xu Guochang. 1997. Climate Change in China Arid and Semi-arid Areas (in Chinese) [M]. Beijing: China Meteorological Press, 101.

    袁景鳳, 吳曉曦. 2000. 一種新的多因子綜合預(yù)報(bào)方法及其預(yù)測(cè)能力評(píng)估 [J]. 應(yīng)用氣象學(xué)報(bào), 11 (S1): 98–102. Yuan Jingfeng, Wu Xiaoxi.2000. A new multifactor synthesize prediction method and its forecast performance assessment [J]. Journal of Applied Meteorological Science(in Chinese), 11(S1): 98–102.

    章基嘉, 孫照渤, 陳松軍. 1981. 對(duì)自然正交函數(shù)穩(wěn)定性條件的討論 [J].氣象學(xué)報(bào), 39 (1): 82–88. Zhang Jijia, Sun Zhaobo, Chen Songjun.1981. On the stability of empirical orthogonal function (EOF) [J]. Acta Meteorologica Sinica (in Chinese), 39 (1): 82–88.

    章基嘉, 孫照渤, 蘭國(guó)明. 1979. 亞洲500毫巴候平均環(huán)流的天氣統(tǒng)計(jì)學(xué)分析 [J]. 南京氣象學(xué)院院報(bào), (1): 33–41. Zhang Jijia, Sun Zhaobo,Lan Guoming. 1979. Weather statistical analysis of Asian 500 mb mean circulations [J]. Journal of Nanjing Institute of Meteorology (in Chinese),(1): 33–41.

    朱乾根, 林錦瑞, 壽紹文, 等. 2007. 天氣學(xué)原理和方法 [M]. 北京: 氣象出版社, 343. Zhu Qian’gen, Lin Jinrui, Shou Shaowen, et al. 2007.Principles and Methods of Meteorology (in Chinese) [M]. Beijing: China Meteorological Press, 343.

    Zhang X S, Wu K J, Wang B, et al. 2013. The low-frequency variance of the ocean surface wave field in the area of the Antarctic Circumpolar Current [J]. Acta Oceanol. Sin., 32 (5): 15–21.

    猜你喜歡
    距平百分率降水
    簡(jiǎn)圖法理解標(biāo)記有絲分裂百分率法測(cè)定細(xì)胞周期時(shí)間
    颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對(duì)比分析
    黑龍江省玉米生長(zhǎng)季自然降水與有效降水對(duì)比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    近40年阿里地區(qū)云量和氣溫的年際變化
    西藏科技(2018年9期)2018-10-17 05:51:30
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    甘肅省降水和冰雹天氣氣候分析
    降水現(xiàn)象儀模擬軟件設(shè)計(jì)與實(shí)現(xiàn)
    基于距平百分率的萬源市氣象干旱灰色預(yù)測(cè)
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    美國(guó)“富人稅”最重的十個(gè)州
    海外星云 (2012年13期)2012-04-29 00:44:03
    观看免费一级毛片| 精品欧美国产一区二区三| 18禁在线播放成人免费| 日韩强制内射视频| 亚洲av成人av| 久久人人精品亚洲av| 日韩欧美 国产精品| 人妻夜夜爽99麻豆av| 小说图片视频综合网站| 黄色一级大片看看| 久久草成人影院| 丝袜美腿在线中文| 国产伦人伦偷精品视频| 色吧在线观看| 精品久久久久久,| 亚洲美女视频黄频| 久久精品国产清高在天天线| 国产精品国产三级国产av玫瑰| 少妇人妻一区二区三区视频| 人妻夜夜爽99麻豆av| 亚洲人成网站在线播放欧美日韩| 精品人妻一区二区三区麻豆 | 看免费成人av毛片| 精品人妻一区二区三区麻豆 | av天堂中文字幕网| 国产毛片a区久久久久| 国产黄片美女视频| 白带黄色成豆腐渣| 中国美白少妇内射xxxbb| 国语自产精品视频在线第100页| 日韩国内少妇激情av| 午夜福利在线观看免费完整高清在 | 午夜激情欧美在线| 男插女下体视频免费在线播放| 三级毛片av免费| 老司机福利观看| 成人欧美大片| 18+在线观看网站| 欧美高清成人免费视频www| 小蜜桃在线观看免费完整版高清| 欧美一级a爱片免费观看看| 久久久精品大字幕| 又黄又爽又免费观看的视频| 免费在线观看影片大全网站| 久久人人精品亚洲av| 美女黄网站色视频| 久久久精品欧美日韩精品| 我要搜黄色片| 色播亚洲综合网| 亚洲精品影视一区二区三区av| 国国产精品蜜臀av免费| 亚洲av一区综合| 国产欧美日韩一区二区精品| 免费看av在线观看网站| 国产高清三级在线| 色综合色国产| 村上凉子中文字幕在线| 日韩欧美三级三区| 少妇裸体淫交视频免费看高清| 欧美色视频一区免费| 中文字幕免费在线视频6| 真人做人爱边吃奶动态| 国产白丝娇喘喷水9色精品| 1000部很黄的大片| 亚洲美女视频黄频| 俄罗斯特黄特色一大片| 一区二区三区免费毛片| 国产69精品久久久久777片| 亚洲在线自拍视频| 国产精品乱码一区二三区的特点| 亚洲中文字幕一区二区三区有码在线看| 日本免费一区二区三区高清不卡| 97人妻精品一区二区三区麻豆| 免费高清视频大片| 深爱激情五月婷婷| 嫩草影院入口| 亚洲七黄色美女视频| 美女cb高潮喷水在线观看| 久久久久性生活片| 91麻豆av在线| 国产在线精品亚洲第一网站| 观看免费一级毛片| 啦啦啦观看免费观看视频高清| 国产激情偷乱视频一区二区| 观看免费一级毛片| 级片在线观看| 欧美日韩亚洲国产一区二区在线观看| 久久草成人影院| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产鲁丝片午夜精品 | 日日夜夜操网爽| 精品一区二区三区av网在线观看| 国产aⅴ精品一区二区三区波| 国产 一区 欧美 日韩| 十八禁国产超污无遮挡网站| 亚洲在线自拍视频| 窝窝影院91人妻| 亚洲av中文字字幕乱码综合| 亚洲经典国产精华液单| 亚洲电影在线观看av| 日韩欧美 国产精品| 悠悠久久av| 免费观看在线日韩| 一区二区三区高清视频在线| 免费看a级黄色片| 五月伊人婷婷丁香| 国内精品美女久久久久久| 亚洲成人久久性| 精品人妻1区二区| 少妇丰满av| 88av欧美| 我要搜黄色片| 麻豆成人午夜福利视频| 夜夜夜夜夜久久久久| 国产精品一区www在线观看 | 国产av在哪里看| 露出奶头的视频| 亚洲成a人片在线一区二区| 欧美日韩国产亚洲二区| 男女下面进入的视频免费午夜| av视频在线观看入口| 中文字幕免费在线视频6| 国产黄片美女视频| 亚洲av.av天堂| 高清毛片免费观看视频网站| 成人欧美大片| 97人妻精品一区二区三区麻豆| 天美传媒精品一区二区| 久久久久精品国产欧美久久久| 三级男女做爰猛烈吃奶摸视频| 99久久精品热视频| 麻豆国产av国片精品| 亚洲欧美日韩高清在线视频| 国产精品综合久久久久久久免费| 我要看日韩黄色一级片| 老司机福利观看| 婷婷精品国产亚洲av| 99久久精品热视频| bbb黄色大片| 亚洲国产日韩欧美精品在线观看| 美女免费视频网站| 成人无遮挡网站| 国产又黄又爽又无遮挡在线| 日本免费一区二区三区高清不卡| 我要搜黄色片| 天堂影院成人在线观看| 日韩一本色道免费dvd| 91久久精品国产一区二区三区| 99热6这里只有精品| 91狼人影院| 精华霜和精华液先用哪个| 亚洲七黄色美女视频| 久久精品国产亚洲av天美| 久久香蕉精品热| 亚洲精品一区av在线观看| 亚洲色图av天堂| 欧美中文日本在线观看视频| 麻豆成人午夜福利视频| 亚洲精品456在线播放app | 午夜免费成人在线视频| 香蕉av资源在线| 欧美性猛交╳xxx乱大交人| 嫁个100分男人电影在线观看| 国产69精品久久久久777片| 免费人成在线观看视频色| h日本视频在线播放| av国产免费在线观看| 我要搜黄色片| 桃色一区二区三区在线观看| 男人狂女人下面高潮的视频| 女同久久另类99精品国产91| 一级毛片久久久久久久久女| 日韩大尺度精品在线看网址| 精华霜和精华液先用哪个| 性色avwww在线观看| a级一级毛片免费在线观看| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区精品| 亚洲欧美日韩东京热| 97超视频在线观看视频| 蜜桃亚洲精品一区二区三区| 亚洲专区中文字幕在线| 俺也久久电影网| 国产精品不卡视频一区二区| 国产精品一区二区三区四区久久| 蜜桃久久精品国产亚洲av| 成年女人看的毛片在线观看| 一级毛片久久久久久久久女| 99精品在免费线老司机午夜| 久久精品国产亚洲av香蕉五月| 可以在线观看的亚洲视频| 日韩,欧美,国产一区二区三区 | 韩国av一区二区三区四区| 久久99热这里只有精品18| 欧美色视频一区免费| 两性午夜刺激爽爽歪歪视频在线观看| 麻豆国产97在线/欧美| 老司机午夜福利在线观看视频| 日本-黄色视频高清免费观看| 亚洲自偷自拍三级| 免费看av在线观看网站| 亚洲第一电影网av| 在线天堂最新版资源| 国产高清激情床上av| 亚洲欧美日韩高清专用| 国产高清不卡午夜福利| 美女黄网站色视频| 蜜桃久久精品国产亚洲av| 亚洲色图av天堂| 久久精品综合一区二区三区| 欧美3d第一页| 久久久久国产精品人妻aⅴ院| АⅤ资源中文在线天堂| 国产69精品久久久久777片| 日本黄色视频三级网站网址| 亚洲无线在线观看| 99热这里只有是精品在线观看| 九九在线视频观看精品| 少妇的逼好多水| 亚洲男人的天堂狠狠| 国产黄片美女视频| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 国产毛片a区久久久久| eeuss影院久久| 亚洲五月天丁香| 少妇人妻精品综合一区二区 | 搡老妇女老女人老熟妇| 亚洲av不卡在线观看| 国产精品乱码一区二三区的特点| 女生性感内裤真人,穿戴方法视频| 久久久国产成人免费| 午夜精品一区二区三区免费看| 人妻夜夜爽99麻豆av| 哪里可以看免费的av片| 亚洲男人的天堂狠狠| 婷婷精品国产亚洲av| 日韩欧美一区二区三区在线观看| 久久久色成人| 国产成年人精品一区二区| 午夜福利视频1000在线观看| 亚洲在线自拍视频| 日韩精品有码人妻一区| 99在线视频只有这里精品首页| 人妻丰满熟妇av一区二区三区| 999久久久精品免费观看国产| 亚洲av第一区精品v没综合| 国产精品久久久久久久久免| 婷婷精品国产亚洲av| 国产伦精品一区二区三区视频9| 久久天躁狠狠躁夜夜2o2o| 熟女人妻精品中文字幕| 国产亚洲精品av在线| 国产真实伦视频高清在线观看 | 91精品国产九色| 一区二区三区激情视频| 日韩强制内射视频| 欧美中文日本在线观看视频| 黄色一级大片看看| 国国产精品蜜臀av免费| 成年女人毛片免费观看观看9| 国产毛片a区久久久久| 国产免费av片在线观看野外av| 在线看三级毛片| 欧美激情久久久久久爽电影| 亚洲18禁久久av| 亚洲国产色片| 看黄色毛片网站| 色在线成人网| 99riav亚洲国产免费| 午夜日韩欧美国产| 精品久久久久久久久亚洲 | 日韩精品中文字幕看吧| 黄色视频,在线免费观看| 日本黄色片子视频| 国产精品一区二区三区四区免费观看 | 看十八女毛片水多多多| 久久久久久久久大av| 亚洲欧美精品综合久久99| .国产精品久久| 国产精品女同一区二区软件 | 三级男女做爰猛烈吃奶摸视频| 波多野结衣高清作品| 女的被弄到高潮叫床怎么办 | 深夜精品福利| 午夜老司机福利剧场| 免费观看精品视频网站| 亚洲国产精品sss在线观看| 简卡轻食公司| 伦理电影大哥的女人| 国产精品自产拍在线观看55亚洲| 特大巨黑吊av在线直播| 黄色配什么色好看| 最近视频中文字幕2019在线8| 国产v大片淫在线免费观看| 亚洲专区国产一区二区| 国产主播在线观看一区二区| 成年免费大片在线观看| 日韩中文字幕欧美一区二区| 最近视频中文字幕2019在线8| 国产精品电影一区二区三区| 亚洲成人精品中文字幕电影| 成年人黄色毛片网站| 啦啦啦观看免费观看视频高清| 久久久久久久久大av| 麻豆成人午夜福利视频| 又黄又爽又免费观看的视频| 免费高清视频大片| 国产色爽女视频免费观看| 欧美又色又爽又黄视频| 久久午夜福利片| 亚洲欧美日韩东京热| 亚州av有码| 麻豆av噜噜一区二区三区| 精品福利观看| 男女之事视频高清在线观看| 天堂影院成人在线观看| 国产精品无大码| 国产淫片久久久久久久久| 亚洲最大成人中文| 哪里可以看免费的av片| 久久久国产成人免费| 动漫黄色视频在线观看| 校园人妻丝袜中文字幕| 99热这里只有是精品50| 国产私拍福利视频在线观看| 亚洲欧美激情综合另类| 精品不卡国产一区二区三区| 亚洲va在线va天堂va国产| 亚洲色图av天堂| 精品久久国产蜜桃| 赤兔流量卡办理| 久久精品国产自在天天线| 国产单亲对白刺激| 九九在线视频观看精品| 成人无遮挡网站| 在线国产一区二区在线| 亚洲av日韩精品久久久久久密| 亚洲美女搞黄在线观看 | 九九久久精品国产亚洲av麻豆| 亚洲av日韩精品久久久久久密| 又黄又爽又刺激的免费视频.| 国产三级在线视频| 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 免费av不卡在线播放| av视频在线观看入口| 中文字幕久久专区| 亚州av有码| 国产单亲对白刺激| 亚洲性久久影院| 日韩强制内射视频| 精品久久久久久久久久免费视频| 波多野结衣高清作品| 中出人妻视频一区二区| 国产三级在线视频| 亚洲av免费高清在线观看| 精品国内亚洲2022精品成人| 亚洲成人中文字幕在线播放| 在线免费十八禁| 亚洲人成网站在线播放欧美日韩| 日本免费a在线| 久久人妻av系列| 悠悠久久av| 午夜精品久久久久久毛片777| 欧美激情国产日韩精品一区| 美女高潮喷水抽搐中文字幕| 日本欧美国产在线视频| 免费观看在线日韩| 国内毛片毛片毛片毛片毛片| 九九爱精品视频在线观看| 国产一区二区在线观看日韩| 最近在线观看免费完整版| 国产伦精品一区二区三区四那| 色哟哟·www| 精品99又大又爽又粗少妇毛片 | 麻豆国产av国片精品| 欧美+日韩+精品| 精品一区二区免费观看| 久久精品国产亚洲网站| 欧美激情久久久久久爽电影| 亚洲精品成人久久久久久| 超碰av人人做人人爽久久| 欧美中文日本在线观看视频| 22中文网久久字幕| 免费不卡的大黄色大毛片视频在线观看 | 最好的美女福利视频网| 一本久久中文字幕| 成熟少妇高潮喷水视频| 看免费成人av毛片| 一进一出抽搐gif免费好疼| 久久精品影院6| 一进一出抽搐动态| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻视频免费看| а√天堂www在线а√下载| 国产69精品久久久久777片| 国产亚洲精品综合一区在线观看| eeuss影院久久| 国产淫片久久久久久久久| 伦理电影大哥的女人| 琪琪午夜伦伦电影理论片6080| 日日夜夜操网爽| 丰满的人妻完整版| 亚洲一区高清亚洲精品| 亚洲电影在线观看av| 久9热在线精品视频| 成人高潮视频无遮挡免费网站| 精品乱码久久久久久99久播| 日本一本二区三区精品| 3wmmmm亚洲av在线观看| 午夜精品在线福利| 在线a可以看的网站| 黄色视频,在线免费观看| 久久精品国产亚洲av涩爱 | 老师上课跳d突然被开到最大视频| 黄色女人牲交| 国产69精品久久久久777片| 国产精品久久久久久久电影| 欧美日韩中文字幕国产精品一区二区三区| 欧美+日韩+精品| 国产精品爽爽va在线观看网站| 国产成人一区二区在线| 国产av不卡久久| 欧美日韩亚洲国产一区二区在线观看| 久久久午夜欧美精品| 一级av片app| 久久热精品热| 亚洲四区av| 亚洲第一区二区三区不卡| 亚洲黑人精品在线| 成人永久免费在线观看视频| 91在线精品国自产拍蜜月| 国产精品野战在线观看| 美女黄网站色视频| 国产av麻豆久久久久久久| 九九久久精品国产亚洲av麻豆| 午夜激情福利司机影院| 少妇人妻一区二区三区视频| x7x7x7水蜜桃| 国产免费一级a男人的天堂| 日本黄色片子视频| 俄罗斯特黄特色一大片| 91麻豆av在线| 国产极品精品免费视频能看的| 久久久久久大精品| 欧洲精品卡2卡3卡4卡5卡区| 色哟哟哟哟哟哟| 我要搜黄色片| 久久精品国产鲁丝片午夜精品 | 动漫黄色视频在线观看| 国产黄a三级三级三级人| 在线看三级毛片| 亚洲av二区三区四区| 在线免费观看的www视频| 精品欧美国产一区二区三| 尾随美女入室| 亚洲国产精品sss在线观看| 欧美黑人巨大hd| x7x7x7水蜜桃| 联通29元200g的流量卡| 亚洲经典国产精华液单| 五月伊人婷婷丁香| 99久久精品国产国产毛片| 国产在线精品亚洲第一网站| 熟女电影av网| 亚洲内射少妇av| 中文字幕精品亚洲无线码一区| 国产精品99久久久久久久久| 久久久成人免费电影| 久久精品国产亚洲网站| 女同久久另类99精品国产91| а√天堂www在线а√下载| 在线观看舔阴道视频| 精品不卡国产一区二区三区| 最近中文字幕高清免费大全6 | 日韩精品中文字幕看吧| 亚洲专区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品影视一区二区三区av| 国产v大片淫在线免费观看| av福利片在线观看| 91麻豆精品激情在线观看国产| 午夜福利成人在线免费观看| 赤兔流量卡办理| 久久国产乱子免费精品| 日韩欧美精品免费久久| 成年版毛片免费区| 中亚洲国语对白在线视频| 国产精品亚洲一级av第二区| 久久午夜亚洲精品久久| 不卡一级毛片| 欧美一区二区国产精品久久精品| 亚洲三级黄色毛片| 一进一出抽搐gif免费好疼| 狠狠狠狠99中文字幕| 淫妇啪啪啪对白视频| 联通29元200g的流量卡| 99久久久亚洲精品蜜臀av| 丰满乱子伦码专区| 欧美中文日本在线观看视频| 亚洲美女黄片视频| 亚洲精华国产精华精| 十八禁国产超污无遮挡网站| 亚洲人成网站在线播放欧美日韩| 色噜噜av男人的天堂激情| 国内精品宾馆在线| 国产精品电影一区二区三区| 免费电影在线观看免费观看| 很黄的视频免费| 黄色一级大片看看| 国产不卡一卡二| 啦啦啦啦在线视频资源| 亚洲四区av| 欧美精品啪啪一区二区三区| 亚洲精品色激情综合| 男插女下体视频免费在线播放| 国产欧美日韩精品一区二区| 一a级毛片在线观看| 亚洲avbb在线观看| 免费在线观看影片大全网站| 成人特级av手机在线观看| 国产毛片a区久久久久| 狠狠狠狠99中文字幕| 少妇的逼水好多| 亚洲黑人精品在线| 欧美最新免费一区二区三区| 国产精品美女特级片免费视频播放器| 女同久久另类99精品国产91| 欧美色欧美亚洲另类二区| 亚洲真实伦在线观看| 久久精品国产鲁丝片午夜精品 | 啦啦啦啦在线视频资源| 午夜福利18| 老女人水多毛片| 我的女老师完整版在线观看| 色5月婷婷丁香| 99riav亚洲国产免费| 亚洲中文字幕日韩| 欧美中文日本在线观看视频| www日本黄色视频网| 日本一本二区三区精品| 美女免费视频网站| 成人综合一区亚洲| 日本爱情动作片www.在线观看 | 欧美日韩综合久久久久久 | 大又大粗又爽又黄少妇毛片口| 3wmmmm亚洲av在线观看| 国产激情偷乱视频一区二区| 色哟哟·www| 黄色日韩在线| 欧美成人一区二区免费高清观看| av福利片在线观看| 夜夜夜夜夜久久久久| 久久国产乱子免费精品| 欧美中文日本在线观看视频| www日本黄色视频网| 精品国内亚洲2022精品成人| 毛片一级片免费看久久久久 | 亚洲va在线va天堂va国产| 国产男靠女视频免费网站| 中文在线观看免费www的网站| a级毛片a级免费在线| 国产熟女欧美一区二区| 一级av片app| 国产欧美日韩精品亚洲av| 黄色配什么色好看| 熟女人妻精品中文字幕| 国产爱豆传媒在线观看| 亚洲不卡免费看| 精华霜和精华液先用哪个| 午夜福利视频1000在线观看| 精品人妻视频免费看| 婷婷丁香在线五月| 大型黄色视频在线免费观看| 99久久成人亚洲精品观看| 国产91精品成人一区二区三区| 国产午夜精品论理片| 我要搜黄色片| 国产精品综合久久久久久久免费| 一进一出抽搐动态| 日韩人妻高清精品专区| 日韩欧美在线二视频| 男女下面进入的视频免费午夜| 国产精品永久免费网站| 成人国产麻豆网| 99久久成人亚洲精品观看| 99热只有精品国产| 国产探花在线观看一区二区| 久久6这里有精品| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲欧美98| 国产三级在线视频| 午夜激情欧美在线| 久久中文看片网| 一本精品99久久精品77| 欧美性感艳星| 久久久久国产精品人妻aⅴ院| 久久亚洲真实| av在线观看视频网站免费| 国产视频内射| 波多野结衣高清作品| 国产探花在线观看一区二区| videossex国产| 亚洲第一区二区三区不卡| 给我免费播放毛片高清在线观看| 在线观看一区二区三区| av黄色大香蕉| av视频在线观看入口| 国产亚洲精品久久久久久毛片| 欧美丝袜亚洲另类 | 91久久精品国产一区二区成人| 波多野结衣高清无吗|