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

    基于BP-CCA的貴州夏季降水與太平洋海溫耦合關(guān)系

    2015-05-04 08:41:59唐浩鵬
    中低緯山地氣象 2015年1期
    關(guān)鍵詞:海溫赤道太平洋

    唐浩鵬,王 芬,2

    (1.貴州省黔西南自治州氣象局,貴州 興義 562400 2;2.貴州省山地氣候與資源重點(diǎn)實(shí)驗(yàn)室,貴州 貴陽(yáng) 550002)

    基于BP-CCA的貴州夏季降水與太平洋海溫耦合關(guān)系

    唐浩鵬1,王 芬1,2

    (1.貴州省黔西南自治州氣象局,貴州 興義 562400 2;2.貴州省山地氣候與資源重點(diǎn)實(shí)驗(yàn)室,貴州 貴陽(yáng) 550002)

    利用貴州83個(gè)臺(tái)站1979—2011年的夏季月降水資料、海溫、氣壓場(chǎng)、風(fēng)場(chǎng)等再分析格點(diǎn)資料及副高指數(shù)資料,利用點(diǎn)相關(guān)找出影響貴州夏季降水的海溫關(guān)鍵區(qū)及關(guān)鍵時(shí)段,利用變形的典型相關(guān)分析(BP-CCA)方法對(duì)貴州夏季降水與關(guān)鍵區(qū)海溫的耦合關(guān)系進(jìn)行了分析,并對(duì)海溫影響貴州夏季降水的可能機(jī)制進(jìn)行了探討。結(jié)果表明:①影響貴州夏季降水的海溫關(guān)鍵區(qū)域?yàn)楸碧窖蠹爸袞|赤道太平洋,其中尤以加利福尼亞冷流區(qū)、北太平洋暖流區(qū)對(duì)貴州夏季降水的影響更為重要,關(guān)鍵時(shí)段為上一年7-9月;②BP-CCA第一對(duì)典型相關(guān)場(chǎng)分析結(jié)果表明:加利福尼亞冷流區(qū)海溫與貴州夏季降水呈同位相變化,而北太平洋暖流區(qū)海溫與貴州夏季降水呈反相位變化,關(guān)鍵區(qū)海溫對(duì)貴州北部降水的耦合好于南部,而對(duì)西南部的耦合最差;③前期7-9月海溫與副高強(qiáng)度指數(shù)的相關(guān)系數(shù)分布呈EI Nino 型,與副高西伸脊線指數(shù)的相關(guān)系數(shù)分布在中低緯度西太平洋呈正相關(guān),與中東赤道太平洋呈負(fù)相關(guān),與副高脊點(diǎn)的相關(guān)系數(shù)與海溫和副高的強(qiáng)度相關(guān)分布相反,在赤道中東太平洋為負(fù)相關(guān)區(qū),西北太平洋為正相關(guān)區(qū),海溫分布型為L(zhǎng)A Nina型。東北太平洋及中東赤道太平洋海溫異常與貴州500 hPa氣壓場(chǎng)及風(fēng)場(chǎng)顯著正相關(guān),而西北太平洋海溫與貴州地區(qū)500 hPa氣壓場(chǎng)及風(fēng)場(chǎng)成負(fù)相關(guān),為典型的EI Nino 型 。

    太平洋海溫;貴州;夏季降水;BP-CCA

    1 引言

    在認(rèn)識(shí)氣候變化過(guò)程中最主要的挑戰(zhàn)是對(duì)各種反饋過(guò)程的定量表述。當(dāng)一個(gè)氣候變量變化時(shí),通過(guò)一定的途徑也改變了另一個(gè)變量,從而對(duì)引起變化的初始變量產(chǎn)生影響。美國(guó)哥倫比亞大學(xué)拉蒙特—多爾蒂地質(zhì)觀象臺(tái)的布魯克爾第一次指出,世界洋流的變化可以解釋過(guò)去一些氣候變化[1]。20世紀(jì)60年代,Bjerkness[2]開(kāi)始對(duì)熱帶東太平洋海氣相互作用進(jìn)行了研究,發(fā)現(xiàn)赤道東太平洋海溫異常會(huì)對(duì)中高緯度地區(qū)甚至全球的大氣環(huán)流和天氣氣候產(chǎn)生影響。近年來(lái),我國(guó)氣象工作者針對(duì)海溫和降水開(kāi)展了一系列的研究工作,一方面研究集中在對(duì)區(qū)域性場(chǎng)用海溫建立預(yù)報(bào)模型及數(shù)據(jù)模擬上,另一方面研究則集中在海溫對(duì)降水的影響規(guī)律研究上。賈小龍等[3]使用變形的典型相關(guān)分析(BP-CCA)建立了東亞冬季500 hPa大尺度環(huán)流和中國(guó)冬季溫度及降水的降尺度預(yù)測(cè)模型。姚愚等[4]利用OSR(最優(yōu)子集回歸)及BP-CCA兩種方法建立了對(duì)降水的預(yù)報(bào)模型。王冰[5]利用典型相關(guān)理論分析了云南5月降水與上年5月至當(dāng)年4月各月太平洋熱帶海溫距平場(chǎng)之間的關(guān)系。胡桂芳[6]建立了位勢(shì)高度場(chǎng)前13個(gè)主分量場(chǎng)與季降水場(chǎng)之間的EOF-CCA模型,并對(duì)模型的擬合、預(yù)測(cè)效果進(jìn)行了評(píng)估。金建德等[7]利用3因子最佳子集回歸把前期海溫場(chǎng)作為預(yù)報(bào)因子與貴州夏季降水求相關(guān),許丹等[8]用合成分析、奇異值分解SVD方法對(duì)冬季、春季北太平洋海溫與貴州夏季降水的關(guān)系進(jìn)行了分析計(jì)算,許可等[9]利用經(jīng)驗(yàn)正交函數(shù)(EOF)對(duì)印度洋海表溫度及其與貴州夏季降水的關(guān)系進(jìn)行了分析。所有這些研究為預(yù)報(bào)貴州夏季降水提供了一定的思路。BP-CCA 來(lái)源于經(jīng)典多元統(tǒng)計(jì)理論,具有嚴(yán)密的數(shù)理基礎(chǔ),而應(yīng)用BP-CCA模型來(lái)研究貴州夏季降水場(chǎng)與前期海溫的關(guān)系的研究尚無(wú),本文將在已有的研究基礎(chǔ)上,選用貴州83站1979—2011年的夏季降水月資料,使用點(diǎn)相關(guān)的分析方法找出影響貴州夏季降水的海溫關(guān)鍵區(qū)及關(guān)鍵時(shí)段,使用BP-CCA模型分析貴州夏季降水與關(guān)鍵區(qū)海溫的關(guān)系,并對(duì)前期海溫異常與貴州夏季降水異常進(jìn)行了合成分析,以揭示前期海溫異常影響貴州夏季降水的可能機(jī)制。

    2 資料與方法

    2.1 資料來(lái)源及處理

    貴州省氣象局提供的83個(gè)氣象臺(tái)站的夏季月降水觀測(cè)資料。

    美國(guó)國(guó)家大氣中心和環(huán)境預(yù)報(bào)中心(NCEP/NCAR)發(fā)布的第二套格點(diǎn)資料。

    美國(guó)國(guó)家海洋大氣局(NOAA)提供的全球逐月平均海表溫度資料。

    國(guó)家氣候中心提供的1951—2011年74項(xiàng)環(huán)流指數(shù)。

    2.2 方法

    本文使用點(diǎn)相關(guān)的分析方法找出影響貴州夏季降水的海溫關(guān)鍵區(qū)及關(guān)鍵時(shí)段,對(duì)前期海溫異常與貴州夏季降水異常進(jìn)行了合成分析,運(yùn)用t檢驗(yàn)方法對(duì)海溫距平等變量場(chǎng)合成平均進(jìn)行顯著性檢驗(yàn),并對(duì)海溫影響降水的可能機(jī)制進(jìn)行了探討。

    BP-CCA方法的主要功能是找出關(guān)鍵區(qū)海溫與貴州夏季降水關(guān)聯(lián)最密切的空間型態(tài)。本文選取的研究對(duì)象為1977—2011年海溫場(chǎng)資料及1979—2011年夏季(6—8月平均)貴州區(qū)域內(nèi)(24.5~29.5°N,104~109°E)的降水量。CCA是氣象上常用的方法之一,它將較多的變量轉(zhuǎn)化為少數(shù)幾個(gè)典型變量,通過(guò)研究典型變量之間的相關(guān)系數(shù),分析兩組變量場(chǎng)之間的關(guān)系。使用典型相關(guān)分析時(shí),要求樣本量大于兩組變量的個(gè)數(shù),如果變量的個(gè)數(shù)大于樣本量時(shí),要減少變量個(gè)數(shù),才能使用CCA方法,這種做法有可能會(huì)損失有價(jià)值的信息,Barnett等[10-11]提出了一種新方法—BP典型相關(guān)分析(BP-CCA):首先對(duì)預(yù)報(bào)因子場(chǎng)和對(duì)像場(chǎng)進(jìn)行EOF分解,截取相應(yīng)的特征向量、時(shí)間系數(shù)及特征根,對(duì)新構(gòu)造的主分量變量場(chǎng)進(jìn)行CCA分析,對(duì)得到的典型相關(guān)系數(shù)進(jìn)行檢驗(yàn)[12-13]。

    計(jì)算步驟概括如下:

    ①對(duì)預(yù)報(bào)因子場(chǎng)A和預(yù)報(bào)對(duì)象場(chǎng)B作標(biāo)準(zhǔn)化或距平處理。

    ②對(duì)①處理后的變量場(chǎng)進(jìn)行EOF分析分解,分別截取前幾個(gè)模態(tài)的時(shí)間系數(shù),構(gòu)造主分量矩陣X和Y。

    ③分別計(jì)算兩主分量X和Y矩陣的協(xié)方差陣syy-1SyxSxx-1Sxy,及交叉協(xié)方差矩陣Sxy及Syx。

    ④用奇異值分解計(jì)算方法求出矩陣Syy-1SyxSxx-1Sxy的特征值λ及對(duì)應(yīng)的載荷特征向量d。

    ⑧對(duì)典型相關(guān)系數(shù)采用大樣本的x2進(jìn)行顯著性檢驗(yàn)。

    3 結(jié)果

    3.1 海溫關(guān)鍵區(qū)和關(guān)鍵時(shí)段的選取

    為找出影響貴州夏季降水的關(guān)鍵海溫區(qū)域及關(guān)鍵時(shí)段,將貴州夏季降水?dāng)?shù)據(jù)進(jìn)行了標(biāo)準(zhǔn)化處理,對(duì)此標(biāo)準(zhǔn)化后的降水場(chǎng)進(jìn)行EOF分解,提取第一時(shí)間序列,與1977年6月—2011年5月超前2 a共24個(gè)月海溫逐月滑動(dòng)求相關(guān),共得到逐月滑動(dòng)的24張相關(guān)系數(shù)空間分布圖(圖略),經(jīng)過(guò)對(duì)比分析, 綜合考慮相關(guān)系數(shù)的大小、穩(wěn)定性、持續(xù)時(shí)間等因素, 可以得到與貴州夏季降水顯著相關(guān)的海域及關(guān)鍵時(shí)段,圖1給出了24張相關(guān)圖中與貴州夏季降水的相關(guān)系數(shù)較顯著的4個(gè)月的空間分布圖。

    綜合24張相關(guān)系數(shù)空間分布圖可知(圖略),影響貴州翌年夏季降水的最佳海溫主要集中在重要洋流區(qū)(黑潮、親潮、印度洋暖池區(qū)、北太平洋暖流區(qū)、加里福尼亞冷流區(qū)、西風(fēng)漂流區(qū)等),這與過(guò)去的研究一致[7-9],但是這些對(duì)貴州夏季降水產(chǎn)生重要影響的洋流區(qū)面積和持續(xù)時(shí)間不盡相同,超前3-6個(gè)月的日本海以東洋面海溫與貴州夏季降水相關(guān)通過(guò)了顯著性檢驗(yàn),但面積小,持續(xù)時(shí)間短,即從上年11月開(kāi)始對(duì)貴州夏季降水影響開(kāi)始明顯,1月與相關(guān)系數(shù)通過(guò)顯著性區(qū)域面積達(dá)到最大,2月其面積又迅速減小,此后對(duì)貴州夏季降水的影響宣告結(jié)束(圖略)。而上年9月印度洋暖池海溫雖與降水有較大面積的顯著負(fù)相關(guān),但是持續(xù)時(shí)間短,僅持續(xù)1個(gè)月,到10月的時(shí)候幾乎消失。加利福尼亞冷流區(qū)及北太平洋暖流區(qū)與降水的相關(guān)系數(shù)通過(guò)顯著性檢驗(yàn)的面積較之日本海附近更大,而且持續(xù)時(shí)間長(zhǎng),達(dá)到4個(gè)月(圖1),而貴州夏季降水與其它海域相關(guān)不夠穩(wěn)定、持續(xù)時(shí)間短, 相關(guān)系數(shù)小。因此, 將影響貴州夏季降水的海溫關(guān)鍵區(qū)初選為位于東北太平洋的加利福尼亞冷流區(qū)及北太平洋暖流區(qū),180° ~100°W,10° ~ 50°N,31*21個(gè)格點(diǎn),共計(jì)651個(gè)格點(diǎn),選定影響貴州夏季降水的顯著相關(guān)時(shí)段為上年的7-9月。

    圖1 海溫與貴州夏季降水的相關(guān)系數(shù)空間分布

    (a圖上年7月,b圖為上年8月,c圖為上年9月,d圖為上年10月,0.34是0.05信度閾值,陰影區(qū)是通過(guò)0. 05 顯著性檢驗(yàn), 方框區(qū)為顯著關(guān)鍵區(qū))

    對(duì)關(guān)鍵區(qū)海溫1978—2010 年共33 a 7-9月標(biāo)準(zhǔn)化場(chǎng)進(jìn)行EOF分析,表2給出了EOF統(tǒng)計(jì)結(jié)果,圖2給出了方差貢獻(xiàn)率較大的前3個(gè)模態(tài)的空間分布,它反映出關(guān)鍵區(qū)海溫的幾種主要的大范圍空間分布型,其累積方差貢獻(xiàn)率為63%。第1 模態(tài)( 圖2a) 占總方差的38%,表明關(guān)鍵區(qū)海溫的空間分布型為全區(qū)基本一致,僅在西北部即北太平洋冷流區(qū)出現(xiàn)小片負(fù)值,高值區(qū)主要位于加利福尼亞冷流區(qū)。第2 模態(tài)( 圖2b) 占總方差的14%,其特征向量場(chǎng)由西北向東南表現(xiàn)為負(fù)—正—負(fù)分布,反映了海溫關(guān)鍵區(qū)西北—東南之間的空間局地差異。第3特征向量由南到北表現(xiàn)為正—負(fù)—正分布。North檢驗(yàn)[14]表明前10個(gè)模態(tài)中僅有4個(gè)模態(tài)通過(guò)了North檢驗(yàn)(其中第2模態(tài)并未通過(guò)檢驗(yàn))。

    表1 海溫關(guān)鍵區(qū)7-9月數(shù)據(jù)EOF分析統(tǒng)計(jì)結(jié)果

    圖2 海溫關(guān)鍵區(qū)的前3模態(tài)空間型態(tài)

    3.2 貴州夏季降水與海溫關(guān)鍵區(qū)的BP-CCA分析

    采用1978—2010年(即樣本數(shù)為33)的7-9月關(guān)鍵區(qū)海溫資料(651個(gè)格點(diǎn))和1979—2011貴州夏季83個(gè)站點(diǎn)降水資料。使用經(jīng)典典型相關(guān)分析要求樣本量大于兩組變量的個(gè)數(shù)或者說(shuō)兩個(gè)變量場(chǎng)的空間格點(diǎn)、站點(diǎn)數(shù),才能保證方差陣可逆,使計(jì)算得以進(jìn)行,因而在做典型相關(guān)分析之前分別對(duì)海溫場(chǎng)和降水場(chǎng)進(jìn)行EOF分解,并截取前15個(gè)主分量作為新的因子場(chǎng),降水場(chǎng)前15主分量累積方差貢獻(xiàn)率為90.97%。為了找出影響貴州夏季降水的最佳時(shí)段,我們將上年7-9月的海溫場(chǎng)采用單月滑動(dòng)、2個(gè)月滑動(dòng)和3個(gè)月滑動(dòng)3種不同時(shí)間步長(zhǎng)與貴州夏季降水作BP-CCA分析,結(jié)果見(jiàn)表2。

    表2 不同時(shí)段海溫與貴州夏季降水的BP-CCA分析結(jié)果

    據(jù)表2可看出,海溫取不同時(shí)段所提供給降水的信息存在差異,從相關(guān)系數(shù)來(lái)看,7月的最高,達(dá)到0.997 2,通過(guò)了0.05的顯著性檢驗(yàn),9月最低為0.97。降水對(duì)9月海溫的解釋方差最高,達(dá)到60%,其次為7-9月,為54%,7月的最低為51%。7-9月海溫對(duì)降水場(chǎng)解釋最高,為62%,其次為8-9月,為61%,9月最低為54%,綜合其相關(guān)系數(shù)及解釋方差,因子的時(shí)間尺度取連續(xù)3個(gè)月(7-9月)的整體效果最好。

    由于用BP-CCA方法建立兩場(chǎng)關(guān)系時(shí)與兩場(chǎng)截取的EOF模態(tài)的個(gè)數(shù)有關(guān),為使建立的二者關(guān)系具有最佳的效果,這里對(duì)1978—2010年7-9月關(guān)鍵區(qū)海溫及1979—2011年貴州6-8月降水分別截取1~32個(gè)EOF模態(tài),重復(fù)進(jìn)行CCA分析,圖3給出了截取不同EOF模態(tài)個(gè)數(shù)時(shí)所對(duì)應(yīng)的第一模態(tài)相關(guān)系數(shù)??梢钥吹剑?dāng)截取的EOF模態(tài)個(gè)數(shù)增加時(shí),模態(tài)相關(guān)系數(shù)是迅速增加的,當(dāng)取14個(gè)模態(tài)后其相關(guān)系數(shù)都在0.99以上,隨后隨著截?cái)嗄B(tài)個(gè)數(shù)的增加,模態(tài)相關(guān)系數(shù)趨于相對(duì)穩(wěn)定。

    圖3 取不同EOF模態(tài)個(gè)數(shù)進(jìn)行CCA重建

    根據(jù)表2及圖3的比較結(jié)果,最終選取1978-2010年7-9月的關(guān)鍵區(qū)海溫資料與1979-2011年貴州夏季降水資料進(jìn)行EOF分解,截取其前15個(gè)模態(tài)構(gòu)建新的主分量場(chǎng),對(duì)此主分量場(chǎng)作CCA分析,由于只有第1模態(tài)通過(guò)了顯著性檢驗(yàn),所以本文只分析第1對(duì)模態(tài),對(duì)其它模態(tài)不予分析。

    圖4 前期7-9月海溫與貴州夏季降水BP-CC的第1對(duì)典型載荷特征向量

    圖4a載荷特征向量反映了海溫場(chǎng)東南為負(fù)值、西北為正值的分布,正值區(qū)位于北太平洋的西風(fēng)漂流區(qū)(北太平洋暖流區(qū)),負(fù)值區(qū)位于加利福尼亞冷流區(qū),且此區(qū)域與沃克環(huán)流下降區(qū)域一致,可認(rèn)為上年7-9月海溫場(chǎng)作用于沃克環(huán)流,通過(guò)這種大氣環(huán)流模式影響貴州翌年夏季的降水[5]。而對(duì)應(yīng)的降水場(chǎng)(圖4b)則表現(xiàn)為全場(chǎng)一致的負(fù)值,即當(dāng)上年7-9月加利福尼亞冷流區(qū)的海溫偏高(低)時(shí),翌年貴州夏季的降水將會(huì)偏多(少),而當(dāng)同期北太平洋暖流區(qū)的海溫偏高(低)時(shí),翌年貴州夏季的降水將會(huì)發(fā)生與之相反的變化,將會(huì)偏少(多)。關(guān)鍵區(qū)海溫對(duì)貴州北部降水的解釋效果好于南部,對(duì)黔西南地區(qū)的解釋效果最差,這與金建德等研究結(jié)果[7]一致。這種空間分布可能由于貴州降水不僅受東南季風(fēng)的影響,還受西南季風(fēng)、南支槽、印度洋海溫等其它因素的影響,因此,要詳盡分析貴州西南部降水的影響因子,還要考慮到其它因素。

    3.3 前期海溫異常影響貴州翌年夏季旱澇的天氣學(xué)成因

    由上述分析可知,前期關(guān)鍵區(qū)海溫異常與貴州夏季降水有很好的對(duì)應(yīng)關(guān)系。這種對(duì)應(yīng)關(guān)系從天氣學(xué)意義上看,是否有較合理解釋?下面對(duì)關(guān)鍵區(qū)海溫與大氣環(huán)流進(jìn)行探討,以揭示海溫影響貴州夏季水的可能機(jī)制。

    圖5 前期7-9月海溫與翌年夏季副高3個(gè)指數(shù)(副高強(qiáng)度、脊線、脊點(diǎn))的相關(guān)系數(shù)

    圖5給出了前期7-9月海溫與夏季西太平洋副高3個(gè)指數(shù)(副高強(qiáng)度、脊線、脊點(diǎn))平均值的相關(guān)系數(shù),可以看到前期7-9月海溫副高強(qiáng)度指數(shù)的相關(guān)系數(shù)分布呈EI Nino 型,在赤道中東太平洋及東北太平洋為大范圍內(nèi)的正相關(guān)區(qū),且通過(guò)了0.05的顯著性檢驗(yàn),即當(dāng)前期7-9月中東赤道太平洋及東北太平洋海溫異常偏高時(shí),次年夏季西太平洋副高強(qiáng)度加強(qiáng)。前期7-9月海溫與副高西伸脊線指數(shù)的相關(guān)系數(shù)分布在中低緯度西太平洋呈正相關(guān),與中東赤道太平洋呈負(fù)相關(guān)。前期7-9月海溫與副高脊點(diǎn)的相關(guān)系數(shù)與海溫和副高的強(qiáng)度相關(guān)分布相反,在赤道中東太平洋為負(fù)相關(guān)區(qū),西北太平洋為正相關(guān)區(qū),海溫分布型為L(zhǎng)A Nina型。綜合西太平洋副高3個(gè)指數(shù),前期7-9月海溫對(duì)夏季副高的影響較大,強(qiáng)度及脊點(diǎn)的顯著相關(guān)區(qū)域在赤道中東太平洋及東北太平洋海區(qū),即當(dāng)赤道中東太平洋及東北太平洋海溫偏高時(shí),西太平洋副高強(qiáng)度偏強(qiáng),位置偏西、偏南,反之,西太平洋副高強(qiáng)度偏弱,位置偏東、偏北。

    對(duì)前期5-9月海溫與貴州夏季同期500 hPa高度場(chǎng)及緯向風(fēng)場(chǎng)求相關(guān)系數(shù) (圖6),由圖6可知,東北太平洋及中東赤道太平洋海溫異常與貴州500 hPa氣壓場(chǎng)及風(fēng)場(chǎng)顯著正相關(guān),而西北太平洋海溫與貴州地區(qū)500 hPa氣壓場(chǎng)及風(fēng)場(chǎng)成負(fù)相關(guān),為典型的EI Nino 型,說(shuō)明前期海溫偏高時(shí),中東赤道太平洋及東北太平洋副高偏強(qiáng),偏西,有利于貴州夏季降水。

    圖6 夏季500 hPa平均氣壓場(chǎng)及風(fēng)場(chǎng)與上年7-9月

    (a:500 hPa氣壓場(chǎng)與海溫相關(guān),b:500 hPa高度上風(fēng)場(chǎng)與海溫相關(guān))

    以1979—2011年33 a貴州夏季降水標(biāo)準(zhǔn)化距平數(shù)據(jù)大于1.0倍標(biāo)準(zhǔn)差為降水偏多年, 小于-1.0倍標(biāo)準(zhǔn)差為降水偏少年的標(biāo)準(zhǔn), 選出降水旱澇年(表3),與選取貴州夏季旱澇年的標(biāo)準(zhǔn)一樣,用此方法選出1978—2010年關(guān)鍵區(qū)海溫的冷(暖)水年(表3)。分析表3可以發(fā)現(xiàn),貴州夏季旱澇年與上年7-9月關(guān)鍵區(qū)海溫冷暖水年有4 a是對(duì)應(yīng)的,分別為1990年的關(guān)鍵區(qū)海溫暖水年對(duì)應(yīng)1991年貴州夏季的澇年,1992年的關(guān)鍵區(qū)海溫暖水年對(duì)應(yīng)1993年貴州夏季的澇年,1988年的關(guān)鍵區(qū)海溫冷水年對(duì)應(yīng)1989年貴州夏季的旱年,2010年的關(guān)鍵區(qū)海溫冷水年對(duì)應(yīng)2011年貴州夏季的旱年,即4個(gè)暖水年對(duì)應(yīng)2個(gè)澇年,占50%,4個(gè)冷水年對(duì)應(yīng)2個(gè)旱年,也占50%。綜上所述,上年7-9月關(guān)鍵區(qū)海溫冷水年預(yù)示翌年貴州夏季偏少,而暖水年則預(yù)示著翌年夏季降水偏多。

    表3 貴州夏季旱澇年及海溫冷暖水年的對(duì)應(yīng)

    根據(jù)前面降水和海溫場(chǎng)的相關(guān)分析發(fā)現(xiàn),貴州夏季降水與前期7-9月海溫場(chǎng)相關(guān)較好,所以選擇選擇旱澇年前期7-9月的海溫做距平合成,并用兩組樣本平均值差異的方法進(jìn)行顯著性檢驗(yàn)(黃嘉佑,2000)。圖7給出了貴州夏季旱澇年前期7-9月海溫距平合成圖及合成差的t檢驗(yàn)。分析圖7可看到旱澇年前期7-9月海溫距平分布有明顯不同。在澇年上年7-9月,東北太平洋、印度洋及赤道中太平洋為負(fù)距平,負(fù)值中心位于東北太平洋,中心海溫偏低1.2℃,加利福尼亞冷流區(qū)、赤道東太平洋及黑潮、親潮交匯處為海溫正距平,中心區(qū)海溫偏高1.2℃(圖7a)。而旱年前期7-9月海溫距平分布與澇年有明顯不同,在赤道中東太平洋及加利福尼亞冷流區(qū)大范圍海溫為負(fù)距平,中心偏低1.6℃,而在西北太平洋、黑潮為正距平區(qū)(圖7b)。分析圖7c旱澇差檢驗(yàn)圖可以看出,在赤道中東太平洋及加利福尼亞冷流區(qū)旱澇年的海溫差值較大,中心達(dá)到了3℃,且部分區(qū)域都通過(guò)0.05的顯著性t檢驗(yàn),這與前面做的點(diǎn)相關(guān)的結(jié)果較一致,說(shuō)明影響貴州夏季旱澇的前期7-9月海溫關(guān)鍵區(qū)主要在赤道中東太平洋及北太平洋海域,其中尤以加利福尼亞冷流區(qū)及赤道中東太平洋的海溫對(duì)貴州夏季旱澇的影響較明顯。澇年前期7-9月北太平洋的海溫分布從西到東為“+-+”分布,中東赤道太平洋從西到東也呈現(xiàn)“+-+”分布,但此分布趨勢(shì)不如北太平洋明顯,而旱年前期7-9月北太平洋的海溫分布從西到東為“-+-”分布,中東赤道太平洋呈現(xiàn)全場(chǎng)較一致的負(fù)值分布,為典型的La Nina型分布。

    4 結(jié)語(yǔ)

    ①點(diǎn)相關(guān)分析結(jié)果顯示:影響貴州夏季降水的海溫關(guān)鍵區(qū)域?yàn)槲挥跂|北太平洋的加利福尼亞冷流區(qū)及北太平洋暖流區(qū),關(guān)鍵時(shí)段為上一年的7-9月。

    圖7 貴州夏季旱澇年前期7-9月海溫距平合成圖及合成差的t檢驗(yàn)

    (a:澇年SST,b:旱年SST,c:澇-旱SST,陰影區(qū)為通過(guò)0.05信度的t檢驗(yàn)區(qū))

    ②BP-CCA第一對(duì)典型相關(guān)場(chǎng)顯示:當(dāng)上年7-9月加利福尼亞冷流區(qū)的海溫偏高(低)時(shí),翌年貴州夏季的降水將會(huì)偏多(少),而當(dāng)上年7-9月北太平洋暖流區(qū)的海溫偏高(低)時(shí),翌年貴州夏季的降水將會(huì)發(fā)生與之相反的變化,將會(huì)偏少(多)。關(guān)鍵區(qū)海溫對(duì)貴州北部降水的耦合好于南部,尤其對(duì)黔西南地區(qū)耦合最差。

    ③前期7-9月海溫與副高強(qiáng)度指數(shù)的相關(guān)系數(shù)分布呈EI Nino 型,與副高西伸脊線指數(shù)的相關(guān)系數(shù)分布在中低緯度西太平洋呈正相關(guān),與中東赤道太平洋呈負(fù)相關(guān),與副高脊點(diǎn)的相關(guān)系數(shù)與海溫和副高的強(qiáng)度相關(guān)分布相反,在赤道中東太平洋為負(fù)相關(guān)區(qū),西北太平洋為正相關(guān)區(qū),海溫分布型為L(zhǎng)A Nina型。即當(dāng)赤道中東太平洋及東北太平洋海溫偏高時(shí),西太平洋副高強(qiáng)度偏強(qiáng),位置偏西、偏南,反之,西太平洋副高強(qiáng)度偏弱,位置偏東、偏北。東北太平洋及中東赤道太平洋海溫異常與貴州500 hPa氣壓場(chǎng)及風(fēng)場(chǎng)顯著正相關(guān),而西北太平洋海溫與貴州地區(qū)500 hPa氣壓場(chǎng)及風(fēng)場(chǎng)成負(fù)相關(guān),為典型的EI Nino 型,說(shuō)明前期海溫偏高時(shí),中東赤道太平洋及東北太平洋副高偏強(qiáng),偏西,有利于貴州夏季降水。

    ④前期海溫異常年與貴州夏季降水異常合成分析結(jié)果:影響貴州夏季旱澇的前期7-9月海溫關(guān)鍵區(qū)主要在北太平洋及赤道中東太平洋,澇年前期7-9月北太平洋的海溫分布從西到東為“+-+”分布,中東赤道太平洋從西到東呈現(xiàn)“+-+分布,而旱年前期7-9月北太平洋的海溫分布從西到東為“-+-”分布,中東赤道太平洋呈現(xiàn)全場(chǎng)較一致的負(fù)分布,為典型的La Nina型分布。

    [1] 威廉·伯勒斯.21世紀(jì)的氣候[M].北京:氣象出版社,2007.

    [2] Bjerkness J, Atmospheric teleconnections from the equatorial Pacific,Mon. Wea. Rev., 1969,97:163-172.

    [3] 賈小龍,陳麗娟,李維京,等.BP-CCA方法用于中國(guó)冬季溫度和降水的可預(yù)報(bào)性研究和降尺度季節(jié)預(yù)測(cè)[J].氣象學(xué)報(bào),2010,68(3):398-410.

    [4] 姚愚,晏紅明,李衛(wèi)紅,等.云南5月降水量可預(yù)報(bào)性的時(shí)空分布[J].廣東氣象,2010,33(1):28-31.

    [5] 王冰, 張健挺.云南降水場(chǎng)與太平洋海溫場(chǎng)典型相關(guān)分析[J].云南大學(xué)學(xué)報(bào),2002,24(6):428:430.

    [6] 胡桂芳,張?zhí)K平,謝考憲.EOF-CCA模型在山東季降水預(yù)測(cè)中的應(yīng)用[J].氣象,2000,26(7):12-16.

    [7] 金建德,嚴(yán)小冬,雷云,等.西北太平洋海溫變化對(duì)貴州夏季降水的影響[J].熱帶氣象學(xué)報(bào),2006,22(2):192-197.

    [8] 許丹,王瑾. 貴州夏季降水場(chǎng)與北太平洋海溫場(chǎng)的非同步相關(guān)研究[J].貴州氣象,2000,24(1):3-7.

    [9] 許可,嚴(yán)小冬.印度洋海溫異常對(duì)貴州夏季降水的影響[J].貴州氣象,2005,29(增刊):12-14.

    [10] Barnett T P.1983.Interaction of the monsoon and Pactific trade wind systems at interannual time scales,PartⅠ:The equatorial zone.Monthly Weather Review,111:756-773.

    [11] Barnett T P,Preisendorfer R.1987.Origins and levels of monthly and seasonal forecast skill for united states surface air temperature determined by canonical correlation analysis.Monthly Weather Review,115:1825-1850.

    [12] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京:氣象出版社,1999.

    [13] 黃嘉佑.氣象統(tǒng)計(jì)分析與預(yù)報(bào)方法[M].北京:氣象出版社,2004.

    [14] North G R,Bell T,Cahalan R,et al.Sampling errors on the estimation of empirical orthogonal function.Mon.Wea.Rev.,1982,110:669-706.

    [15] 李志方.太平洋海溫與我國(guó)大陸降水的EOF分析[J].貴州氣象,2012,36(1):18-20.

    [16] 舒國(guó)勇,高紅梅,嚴(yán)小冬.貴州省銅仁市雨季開(kāi)始期演變特征及環(huán)流分析[J]. 貴州氣象,2012,36(5):9-13.

    [17] 楊洋,李光強(qiáng). 厄爾尼諾與貴州夏季降水的關(guān)系[J]. 貴州氣象,1996,19(2):12-15.

    Coupling relationship between the Guizhou summer precipitation and the surface temperature in Pacific Ocean based on BP-CCA

    TANG Haopeng1,WANG Fen1,2

    (1.Meteorological Office of Southwestern Guizhou,Xingyi 562400;2.The mountainous climate and resource key laboratories of Guizhou Province,Guiyang 550002)

    Based on the precipitation of summer at 83 stations from 1979 to 2011 over Guizhou and the mean monthly reanalysis SST ,pressure, wind data of NCEP/ NCAR , the relationship between summer precipitation and the SST,in the Pacific Ocean are investigated in terms of empirical orthog onal function( EOF) , composite analysis, BP-CCA ,and so on. The results are as follows: ①the key SST region which influences lies in the north Pacific Ocean and middle-east Pacific Ocean, the key period is from July of the former year to September;②the results of BP-CCA :the summer precipitation is increasing ( decreasing ) in the corresponding period when SST of the california cold curren from Jul to Sep were warmer ( colder) than average,but when the SST of north pacific ocean warm current were warmer ( colder) than average, the summer precipitation over Guizhou is decreasing ( increasing ) ,the influence of the key SST region to Guizhou is better in the middle-north of Guizhou, but,best in the west-south of Guizhou;③From Early July to September, the correlation coefficient distribution of SST and subtropical high strength index Was EL Nino, the correlation coefficient distribution of SST and subtropical high of west ridge line index in the low latitudes western pacific was positive, but which was negative in middle east equatorial pacific,the correlation coefficient distribution of SST and subtropical high of west ridge point index in the middle-east equatorial pacific was negative, but which was positive in west-north pacific. The distribution of SST was LA Nina. the correlation coefficient between northeast pacific and the middle east equatorial pacific SST and 500hPa wind and pressure in Guizhou was significantly positive, the correlation coefficient between northwest pacific SST and 500hPa wind and pressure in Guizhou was negative, which was for typical EL Nino.

    Pacific SST field;Guizhou;Precipitation in flood season;BP-CCA

    2014-07-27

    唐浩鵬(1985—),男,助工,主要從事短期氣象預(yù)報(bào)預(yù)測(cè)工作。

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41375097);黔西南州科技局社會(huì)發(fā)展攻關(guān)項(xiàng)目(2014-32);黔西南州氣象局項(xiàng)目(2014-08,2014-09)共同資助。

    1003-6598(2015)01-0014-07

    P466

    B

    猜你喜歡
    海溫赤道太平洋
    “雪龍2”號(hào)過(guò)赤道,救下一船人!
    軍事文摘(2024年2期)2024-01-10 01:58:44
    愛(ài)管閑事的“太平洋警察”
    基于深度學(xué)習(xí)的海溫觀測(cè)數(shù)據(jù)質(zhì)量控制應(yīng)用研究
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    決勝太平洋
    南印度洋偶極型海溫與中國(guó)西南地區(qū)初秋降水的關(guān)系
    距赤道越遠(yuǎn)越喜愛(ài)黃色
    最熱的地方不在赤道
    跨越太平洋的愛(ài)戀
    太平洋還是北冰洋
    netflix在线观看网站| 成年人免费黄色播放视频| 村上凉子中文字幕在线| 国产精品98久久久久久宅男小说| 人人妻人人添人人爽欧美一区卜| 日韩一卡2卡3卡4卡2021年| 香蕉丝袜av| 在线观看免费视频网站a站| 欧美亚洲日本最大视频资源| 搡老熟女国产l中国老女人| 国产精品 欧美亚洲| 一进一出抽搐动态| 色精品久久人妻99蜜桃| 美女高潮到喷水免费观看| 国产又爽黄色视频| 飞空精品影院首页| 高清欧美精品videossex| 久久久久精品国产欧美久久久| 久久久国产成人免费| 一夜夜www| 国产欧美亚洲国产| 亚洲精品美女久久久久99蜜臀| 女人被躁到高潮嗷嗷叫费观| 老司机亚洲免费影院| 亚洲情色 制服丝袜| 国产区一区二久久| 香蕉国产在线看| 国产亚洲精品久久久久5区| 国产不卡av网站在线观看| 黄色 视频免费看| 无限看片的www在线观看| 久久热在线av| 一级片'在线观看视频| 在线观看免费日韩欧美大片| 99在线人妻在线中文字幕 | 国产成人精品在线电影| 一级a爱片免费观看的视频| 建设人人有责人人尽责人人享有的| 中文字幕高清在线视频| 老司机午夜福利在线观看视频| av超薄肉色丝袜交足视频| 精品国产一区二区三区久久久樱花| 黄片小视频在线播放| 久久狼人影院| 嫩草影视91久久| 妹子高潮喷水视频| 午夜91福利影院| 9色porny在线观看| 国产亚洲欧美精品永久| 中文欧美无线码| 中文字幕另类日韩欧美亚洲嫩草| 人人妻人人澡人人爽人人夜夜| 精品人妻在线不人妻| 极品少妇高潮喷水抽搐| 在线观看免费视频日本深夜| 好看av亚洲va欧美ⅴa在| 黄片播放在线免费| 91大片在线观看| 最近最新中文字幕大全电影3 | 热99国产精品久久久久久7| 国产日韩一区二区三区精品不卡| 大码成人一级视频| 成人18禁在线播放| 亚洲国产中文字幕在线视频| 亚洲成a人片在线一区二区| 午夜老司机福利片| 桃红色精品国产亚洲av| 免费看十八禁软件| 一级,二级,三级黄色视频| 日本a在线网址| 日韩三级视频一区二区三区| 久久久精品免费免费高清| 亚洲伊人色综图| 成人av一区二区三区在线看| 天堂俺去俺来也www色官网| 亚洲国产精品合色在线| 老司机深夜福利视频在线观看| 天天躁日日躁夜夜躁夜夜| 国产成人免费无遮挡视频| 久久久久久久久免费视频了| 国产在线精品亚洲第一网站| svipshipincom国产片| 多毛熟女@视频| 久久草成人影院| 欧美黄色淫秽网站| 国产精华一区二区三区| 亚洲精品自拍成人| 最近最新免费中文字幕在线| 亚洲黑人精品在线| 大型av网站在线播放| 国产精品自产拍在线观看55亚洲 | 成年女人毛片免费观看观看9 | 欧美久久黑人一区二区| 日韩有码中文字幕| 桃红色精品国产亚洲av| 交换朋友夫妻互换小说| 少妇裸体淫交视频免费看高清 | 久久香蕉国产精品| 人妻 亚洲 视频| 国产高清激情床上av| 天天躁狠狠躁夜夜躁狠狠躁| 涩涩av久久男人的天堂| 自线自在国产av| 欧美 日韩 精品 国产| 丝袜在线中文字幕| 午夜视频精品福利| 国产欧美日韩一区二区三区在线| 国产精品久久久久久精品古装| 不卡一级毛片| 女人被躁到高潮嗷嗷叫费观| 免费在线观看影片大全网站| 又黄又爽又免费观看的视频| 久久中文字幕一级| 80岁老熟妇乱子伦牲交| 欧美日本中文国产一区发布| 老司机午夜福利在线观看视频| 老鸭窝网址在线观看| 日韩欧美免费精品| 日本a在线网址| 人成视频在线观看免费观看| av网站免费在线观看视频| 久久久久久久精品吃奶| 交换朋友夫妻互换小说| 交换朋友夫妻互换小说| 一级片免费观看大全| 男女下面插进去视频免费观看| av超薄肉色丝袜交足视频| 这个男人来自地球电影免费观看| 18禁裸乳无遮挡免费网站照片 | 亚洲精品在线观看二区| 男人舔女人的私密视频| 久久精品国产a三级三级三级| 亚洲成人免费av在线播放| 欧美在线黄色| 一区福利在线观看| av视频免费观看在线观看| 宅男免费午夜| 捣出白浆h1v1| 久久中文看片网| 人妻一区二区av| 在线天堂中文资源库| 亚洲人成77777在线视频| 亚洲av第一区精品v没综合| 久久99一区二区三区| 国产精品久久视频播放| 18在线观看网站| 国产色视频综合| 午夜福利在线观看吧| av网站免费在线观看视频| 欧美人与性动交α欧美精品济南到| 国产极品粉嫩免费观看在线| 女人精品久久久久毛片| 国产精品 欧美亚洲| 亚洲va日本ⅴa欧美va伊人久久| 久久久精品区二区三区| 成人特级黄色片久久久久久久| 老司机靠b影院| 精品电影一区二区在线| 涩涩av久久男人的天堂| 国产又爽黄色视频| 两个人免费观看高清视频| 免费在线观看完整版高清| 啦啦啦 在线观看视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲中文字幕日韩| 午夜两性在线视频| 中文字幕制服av| 亚洲伊人色综图| 巨乳人妻的诱惑在线观看| 少妇被粗大的猛进出69影院| 他把我摸到了高潮在线观看| a级片在线免费高清观看视频| 午夜福利,免费看| 岛国毛片在线播放| 亚洲精品国产精品久久久不卡| videosex国产| 91av网站免费观看| 久久精品国产亚洲av高清一级| 国产黄色免费在线视频| tocl精华| 免费在线观看亚洲国产| 亚洲综合色网址| 国产高清videossex| 50天的宝宝边吃奶边哭怎么回事| 一区二区三区国产精品乱码| 超碰成人久久| 麻豆成人av在线观看| 男人舔女人的私密视频| 两个人免费观看高清视频| 人人妻人人爽人人添夜夜欢视频| 大香蕉久久网| av网站免费在线观看视频| 国产熟女午夜一区二区三区| 中文字幕色久视频| 91大片在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲五月色婷婷综合| 午夜精品国产一区二区电影| 免费少妇av软件| 人人妻人人爽人人添夜夜欢视频| 国产一区二区三区在线臀色熟女 | 欧美日韩亚洲综合一区二区三区_| 电影成人av| 别揉我奶头~嗯~啊~动态视频| 欧美一级毛片孕妇| 我的亚洲天堂| 一本大道久久a久久精品| 丰满饥渴人妻一区二区三| 91字幕亚洲| 日韩欧美一区二区三区在线观看 | 超色免费av| 老司机午夜十八禁免费视频| 99久久综合精品五月天人人| 中文字幕av电影在线播放| 老熟妇乱子伦视频在线观看| 身体一侧抽搐| 国产又爽黄色视频| 欧美日本中文国产一区发布| 狠狠狠狠99中文字幕| 亚洲色图 男人天堂 中文字幕| 精品欧美一区二区三区在线| 夜夜躁狠狠躁天天躁| 不卡av一区二区三区| 身体一侧抽搐| 国产不卡一卡二| 女人久久www免费人成看片| 老鸭窝网址在线观看| 亚洲成a人片在线一区二区| www.熟女人妻精品国产| 色播在线永久视频| 黄网站色视频无遮挡免费观看| 国产日韩欧美亚洲二区| 一进一出好大好爽视频| 99精国产麻豆久久婷婷| 一级,二级,三级黄色视频| 精品久久久久久电影网| 黑人猛操日本美女一级片| 亚洲欧美一区二区三区久久| 精品一区二区三区视频在线观看免费 | 亚洲七黄色美女视频| av福利片在线| www.精华液| 精品久久久精品久久久| 久久久久久免费高清国产稀缺| 一级a爱片免费观看的视频| 啦啦啦 在线观看视频| 色精品久久人妻99蜜桃| 操美女的视频在线观看| 午夜免费鲁丝| 一边摸一边抽搐一进一出视频| 亚洲av成人不卡在线观看播放网| 18禁国产床啪视频网站| 在线看a的网站| 王馨瑶露胸无遮挡在线观看| 超碰成人久久| 999久久久国产精品视频| 久久精品国产亚洲av香蕉五月 | 精品少妇久久久久久888优播| 日韩免费av在线播放| 久热这里只有精品99| 国产一区二区激情短视频| 亚洲精华国产精华精| 99热网站在线观看| www.熟女人妻精品国产| 久久久久精品国产欧美久久久| 黄色怎么调成土黄色| 在线观看免费高清a一片| 国产精品99久久99久久久不卡| videos熟女内射| 19禁男女啪啪无遮挡网站| 久久久久国内视频| av网站在线播放免费| 欧美日韩国产mv在线观看视频| 超碰97精品在线观看| 一进一出好大好爽视频| 亚洲欧美一区二区三区黑人| 热99re8久久精品国产| 最近最新中文字幕大全电影3 | 国产精品 欧美亚洲| 91精品国产国语对白视频| 麻豆成人av在线观看| 国产又色又爽无遮挡免费看| 国产亚洲欧美在线一区二区| 丝袜人妻中文字幕| 午夜久久久在线观看| 国产亚洲精品一区二区www | 丝袜人妻中文字幕| 老司机午夜福利在线观看视频| 在线观看66精品国产| 国产xxxxx性猛交| 国产av精品麻豆| 欧美日韩福利视频一区二区| 国产日韩一区二区三区精品不卡| 亚洲一区二区三区不卡视频| 国产高清视频在线播放一区| 国产精品免费一区二区三区在线 | 老司机午夜十八禁免费视频| 人人妻人人添人人爽欧美一区卜| 黄色成人免费大全| 欧美在线黄色| 久久人妻熟女aⅴ| 欧美黄色淫秽网站| 女人爽到高潮嗷嗷叫在线视频| 男人舔女人的私密视频| 老汉色av国产亚洲站长工具| 99riav亚洲国产免费| 久久久国产一区二区| 亚洲自偷自拍图片 自拍| 十八禁高潮呻吟视频| 成人亚洲精品一区在线观看| 精品无人区乱码1区二区| 精品国产乱子伦一区二区三区| 9热在线视频观看99| tube8黄色片| 黄色视频,在线免费观看| 欧美日韩福利视频一区二区| 脱女人内裤的视频| 建设人人有责人人尽责人人享有的| 午夜福利视频在线观看免费| 国产精品一区二区在线观看99| 午夜精品久久久久久毛片777| 国产日韩一区二区三区精品不卡| 欧美不卡视频在线免费观看 | 中文字幕人妻丝袜制服| 久9热在线精品视频| 淫妇啪啪啪对白视频| 亚洲欧美精品综合一区二区三区| 亚洲国产看品久久| 18禁美女被吸乳视频| 首页视频小说图片口味搜索| 看片在线看免费视频| 超碰97精品在线观看| 日韩欧美一区二区三区在线观看 | 久99久视频精品免费| 夜夜躁狠狠躁天天躁| 日韩中文字幕欧美一区二区| 精品国产美女av久久久久小说| 亚洲精华国产精华精| 国产精品久久电影中文字幕 | 亚洲熟女毛片儿| 亚洲成人国产一区在线观看| 一级片免费观看大全| 欧洲精品卡2卡3卡4卡5卡区| 老司机亚洲免费影院| 久久国产乱子伦精品免费另类| 久久中文字幕一级| 国产单亲对白刺激| 老司机深夜福利视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩中文字幕国产精品一区二区三区 | 国产91精品成人一区二区三区| 久久精品91无色码中文字幕| av一本久久久久| 欧美成人午夜精品| 国产精品二区激情视频| 波多野结衣一区麻豆| 日日摸夜夜添夜夜添小说| 久久人人97超碰香蕉20202| 99精品欧美一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 99国产精品一区二区蜜桃av | 国产又爽黄色视频| 亚洲欧美色中文字幕在线| 国产亚洲欧美在线一区二区| 免费在线观看黄色视频的| 中文字幕高清在线视频| 极品少妇高潮喷水抽搐| 精品免费久久久久久久清纯 | 国产高清激情床上av| 伦理电影免费视频| 亚洲欧美日韩高清在线视频| 免费在线观看影片大全网站| 久久精品人人爽人人爽视色| av福利片在线| 在线观看舔阴道视频| 午夜免费成人在线视频| 少妇裸体淫交视频免费看高清 | 人妻丰满熟妇av一区二区三区 | 日韩欧美一区视频在线观看| 亚洲黑人精品在线| 搡老乐熟女国产| 欧美亚洲日本最大视频资源| 国产欧美亚洲国产| 日本一区二区免费在线视频| 一区二区三区国产精品乱码| 亚洲全国av大片| 久久中文字幕一级| 老司机午夜福利在线观看视频| 国产精品98久久久久久宅男小说| 丰满的人妻完整版| 成人手机av| 91精品国产国语对白视频| tube8黄色片| 母亲3免费完整高清在线观看| 成人影院久久| 国产亚洲一区二区精品| 岛国毛片在线播放| 精品久久久久久电影网| 亚洲精品在线美女| 18禁美女被吸乳视频| 中文字幕av电影在线播放| 亚洲成人国产一区在线观看| 亚洲欧美激情综合另类| www日本在线高清视频| 人妻 亚洲 视频| 男人操女人黄网站| 国产成人精品久久二区二区免费| 免费看a级黄色片| 男女之事视频高清在线观看| 欧美中文综合在线视频| 999精品在线视频| 亚洲人成77777在线视频| 中出人妻视频一区二区| 变态另类成人亚洲欧美熟女 | 久久久国产一区二区| 欧美成人午夜精品| 日韩熟女老妇一区二区性免费视频| 在线av久久热| 久久香蕉精品热| 午夜福利在线免费观看网站| 亚洲综合色网址| 一边摸一边做爽爽视频免费| 亚洲中文av在线| 狂野欧美激情性xxxx| 久久香蕉精品热| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看 | 窝窝影院91人妻| 国产日韩欧美亚洲二区| 大型av网站在线播放| 久久九九热精品免费| 成在线人永久免费视频| 啦啦啦免费观看视频1| 黄网站色视频无遮挡免费观看| 中文字幕av电影在线播放| 久久 成人 亚洲| 色在线成人网| 一个人免费在线观看的高清视频| 婷婷丁香在线五月| 午夜两性在线视频| 亚洲av成人av| 高清av免费在线| 超碰97精品在线观看| 久久亚洲精品不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久精品久久久| 在线免费观看的www视频| 国产麻豆69| 亚洲熟女毛片儿| 亚洲中文av在线| 校园春色视频在线观看| 丝袜在线中文字幕| 最近最新中文字幕大全免费视频| 怎么达到女性高潮| 久久久国产欧美日韩av| 成年人黄色毛片网站| 色播在线永久视频| 午夜福利一区二区在线看| 视频区图区小说| 一级a爱视频在线免费观看| 美女高潮到喷水免费观看| 欧洲精品卡2卡3卡4卡5卡区| 老熟妇仑乱视频hdxx| cao死你这个sao货| 久久人妻av系列| 99国产极品粉嫩在线观看| 久久精品亚洲av国产电影网| 成人精品一区二区免费| 丰满的人妻完整版| 日本撒尿小便嘘嘘汇集6| 国产视频一区二区在线看| 岛国毛片在线播放| 国内久久婷婷六月综合欲色啪| ponron亚洲| 日韩人妻精品一区2区三区| 搡老熟女国产l中国老女人| 日韩欧美一区二区三区在线观看 | 成熟少妇高潮喷水视频| 国产高清视频在线播放一区| 天天操日日干夜夜撸| 一级毛片精品| 久久久精品免费免费高清| 国产成人免费观看mmmm| 19禁男女啪啪无遮挡网站| 精品第一国产精品| 亚洲精品一二三| 国产高清videossex| 久久狼人影院| 三级毛片av免费| 欧美成人午夜精品| 91av网站免费观看| 在线观看舔阴道视频| 大香蕉久久网| 久久影院123| 欧美在线黄色| 三级毛片av免费| 亚洲中文字幕日韩| 久久99一区二区三区| 美女福利国产在线| 欧美午夜高清在线| 亚洲成国产人片在线观看| 在线观看www视频免费| 18禁美女被吸乳视频| 淫妇啪啪啪对白视频| 国产高清激情床上av| 99热网站在线观看| 欧美 日韩 精品 国产| 欧美亚洲 丝袜 人妻 在线| 欧美精品啪啪一区二区三区| 国产精品九九99| 中文字幕最新亚洲高清| 日本欧美视频一区| 国产麻豆69| 婷婷成人精品国产| 最近最新免费中文字幕在线| 久久精品亚洲精品国产色婷小说| 欧美日本中文国产一区发布| www.熟女人妻精品国产| 丁香欧美五月| 亚洲欧美激情综合另类| 18禁观看日本| 成年人黄色毛片网站| 777久久人妻少妇嫩草av网站| 丝袜美腿诱惑在线| 成人精品一区二区免费| 国产精品二区激情视频| 伦理电影免费视频| 国产淫语在线视频| 18禁裸乳无遮挡动漫免费视频| 日韩免费av在线播放| 欧美成狂野欧美在线观看| 变态另类成人亚洲欧美熟女 | 午夜免费鲁丝| 1024香蕉在线观看| 亚洲国产中文字幕在线视频| 波多野结衣av一区二区av| 午夜老司机福利片| 欧美成人午夜精品| 一级,二级,三级黄色视频| 热99国产精品久久久久久7| 国产av精品麻豆| 色综合婷婷激情| 日韩熟女老妇一区二区性免费视频| 一级毛片女人18水好多| 99re在线观看精品视频| 狂野欧美激情性xxxx| 村上凉子中文字幕在线| 精品一区二区三区四区五区乱码| 超碰97精品在线观看| 嫁个100分男人电影在线观看| 亚洲一区中文字幕在线| 日本a在线网址| 在线观看舔阴道视频| 色94色欧美一区二区| 亚洲 欧美一区二区三区| 美女高潮喷水抽搐中文字幕| 亚洲欧洲精品一区二区精品久久久| 精品视频人人做人人爽| 亚洲第一欧美日韩一区二区三区| 久久精品国产a三级三级三级| 亚洲色图 男人天堂 中文字幕| 麻豆乱淫一区二区| 久久精品成人免费网站| 九色亚洲精品在线播放| 夜夜爽天天搞| 老司机亚洲免费影院| 国产成人精品在线电影| 黄色毛片三级朝国网站| 岛国毛片在线播放| 久久青草综合色| 亚洲久久久国产精品| av一本久久久久| 在线永久观看黄色视频| 国内久久婷婷六月综合欲色啪| 国产主播在线观看一区二区| 大码成人一级视频| 无遮挡黄片免费观看| 久久久久视频综合| 无遮挡黄片免费观看| 中出人妻视频一区二区| 国产成人系列免费观看| 两人在一起打扑克的视频| 乱人伦中国视频| 午夜老司机福利片| 99国产极品粉嫩在线观看| 性少妇av在线| 亚洲精品在线观看二区| 亚洲免费av在线视频| 欧美av亚洲av综合av国产av| 久久草成人影院| 欧美av亚洲av综合av国产av| 少妇的丰满在线观看| 免费少妇av软件| 欧美成人免费av一区二区三区 | 黄色视频不卡| 一级毛片高清免费大全| 国产免费现黄频在线看| 亚洲成人手机| 黑丝袜美女国产一区| 成人永久免费在线观看视频| 狠狠婷婷综合久久久久久88av| 咕卡用的链子| 久久人妻福利社区极品人妻图片| 亚洲自偷自拍图片 自拍| 交换朋友夫妻互换小说| 国产精品99久久99久久久不卡| 亚洲av欧美aⅴ国产| 又紧又爽又黄一区二区| 久久性视频一级片| 国产高清视频在线播放一区| 无人区码免费观看不卡| 新久久久久国产一级毛片| 老司机午夜福利在线观看视频| av在线播放免费不卡| 最新美女视频免费是黄的|