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

    華南前汛期降水的氣候特征研究

    2015-03-01 02:20:43李一鳴范華良顧志夏王迪菲
    安徽農(nóng)業(yè)科學(xué) 2015年4期
    關(guān)鍵詞:散度華南高值

    劉 輝,李一鳴,范華良,顧志夏,王迪菲

    (蘇州市防雷中心,江蘇蘇州 215000)

    ?

    華南前汛期降水的氣候特征研究

    劉 輝,李一鳴,范華良,顧志夏,王迪菲

    (蘇州市防雷中心,江蘇蘇州 215000)

    利用1968~2007年華南94站逐日降水資料,通過EOF、REOF分析表明,EOF分析得到的前4個典型場可以很好地反映華南前汛期降水空間分布的異常結(jié)構(gòu),分別為整體一致的空間結(jié)構(gòu)、東西相反的空間結(jié)構(gòu)、南北相反的空間結(jié)構(gòu)、中間與周圍地區(qū)相反變化的空間結(jié)構(gòu),且通過分析NCEP/NCAR資料場發(fā)現(xiàn)不同的環(huán)流形勢是引起模態(tài)不同的原因。REOF旋轉(zhuǎn)后的前5個空間模態(tài)可較好地為華南前汛期降水進(jìn)行氣候區(qū)劃,分為華南東南部、華南中部、華南東北部、華南西南部和華南西北部這5個區(qū)域,具體分析了與EOF第4模態(tài)類似的REOF第2模態(tài)的環(huán)流成因,并發(fā)現(xiàn)此模態(tài)與2月南方濤動指數(shù)和1月Nino3.4指數(shù)存在正相關(guān)性。

    華南;前汛期降水;氣候特征;EOF;REOF

    華南地區(qū)位于我國最南部,是指武夷山-南嶺以南的廣西、廣東、福建和海南等省區(qū)區(qū)域,主要受熱帶季風(fēng)和副熱帶季風(fēng)的共同影響,是我國汛期最長、降水最多的地區(qū),也是我國旱澇災(zāi)害發(fā)生頻繁的地區(qū)之一,其中以華南前汛期的降水最為典型。嚴(yán)重的自然災(zāi)害對人民生活以及國家經(jīng)濟(jì)、社會的發(fā)展等均會帶來不容忽視的影響。因此,關(guān)于華南前汛期降水的研究一直受到氣象工作者的重視[1-2]。

    華南前汛期的降水以5、6月份為主,研究華南地區(qū)的降水氣候特征對旱澇的發(fā)生規(guī)律有一定的指示作用,同時對預(yù)報華南地區(qū)的氣候也起著至關(guān)重要的作用。因此,近百年以來國內(nèi)外都有著大量研究[3-7],如池艷珍等利用1957~2001年華南地區(qū)74個測站逐日降水資料和同期NCEP/NCAR逐日再分析格點資料,對華南前汛期(4~6月)不同的降水時段特征進(jìn)行了比較,發(fā)現(xiàn)華南前汛期(4~6月)存在2個降水比較集中的階段,即存在有雙峰,鋒面降水大致集中在4月,夏季風(fēng)降水盛期集中于6月[6];林愛蘭等分析華南前汛期旱澇變化及趨勢表明,華南經(jīng)歷了多個干-濕轉(zhuǎn)換的時期,形成的華南月、季尺度干濕的環(huán)流系統(tǒng)有明顯差異[7-8];馬慧等利用國家氣候中心整編的1951~2000年中國160個站的月降水資料,選取了16個代表站研究年降水和前汛期降水,發(fā)現(xiàn)華南地區(qū)降水量的空間分布不太均勻,各站與周圍站點的降水量差值也均比較大,前汛期和年降水極大值均在河源,極小值均在百色,華南地區(qū)降水量東部比西部大[9]。盡管大家對華南前汛期降水的空間和時間分布的規(guī)律已經(jīng)做過不少研究,但大多都是從氣候統(tǒng)計學(xué)的角度來進(jìn)行,筆者利用EOF和REOF 2種分析方法,探究華南前汛期降水的氣候特征,重點結(jié)合環(huán)流形勢、相關(guān)物理量以及相關(guān)性研究分析空間分布規(guī)律的原因,給出華南地區(qū)前汛期降水的分布型,為預(yù)報其降水提供依據(jù)。

    1 資料和方法

    1.1 降水資料采用的資料包括中國氣象局提供的1968~2007年每年5月10日~6月10日的華南94站(圖1)的逐日降水資料;1968~2007年的NCEP/NCAR月平均分析資料,空間分辨率為2.5°×2.5°。美國氣候預(yù)測中心提供的1968~2007年的南方濤動指數(shù),該指數(shù)反映了太平洋東西兩側(cè)氣壓的協(xié)同變化;美國氣候預(yù)測中心提供的1968~2007年的Nino3.4指數(shù),即Nino3區(qū)和Nino 4區(qū)的平均海溫距平。

    圖1 華南94站分布示意圖

    1.2 分析方法

    1.2.1EOF方法。也稱為自然正交函數(shù)[10]的分解,用EOF方法可將m個站點、n次觀測的氣象要素Xmn分解為包含空間系數(shù)矩陣和時間系數(shù)矩陣的形式。用EOF展開后的系數(shù)T是時間的函數(shù),所以T的時間變程中就反映場隨時間的變化,而特征向量V只是空間坐標(biāo)函數(shù),一經(jīng)求出就不隨時間變化。由于EOF方法的收斂速度很快,所以取前幾個特征值較大的特征向量就可以充分描述展開要素場的主要特征[11]。

    1.2.2REOF方法。它是旋轉(zhuǎn)經(jīng)驗正交展開方法。氣象上經(jīng)常采用REOF方法對氣象要素場進(jìn)行分析,它使原來的特征向量結(jié)構(gòu)簡化,反映的氣候特征更明顯。它具體是在EOF的經(jīng)驗基礎(chǔ)上,選擇一個正交旋轉(zhuǎn)矩陣,使原始矩陣旋轉(zhuǎn)以后的列向量元素平方的方差達(dá)最大,從而使原要素場的信息特征集中映射到荷載場所表示的優(yōu)勢空間上;旋轉(zhuǎn)后每一個空間點僅與一個主成分存在高相關(guān),旋轉(zhuǎn)后的特征場比旋轉(zhuǎn)前在時間上更穩(wěn)定[12]。根據(jù)REOF理論,可以證明各主分量的最高荷載中心實際上可作為聚類分區(qū)的中心站點[13]。

    1.2.4相關(guān)系數(shù)r的計算。相關(guān)系數(shù)有-1≤r≤1,當(dāng)r為正時,表示有正相關(guān);r為負(fù)時,有負(fù)相關(guān);當(dāng)r=±1時,表示有一一對應(yīng)的函數(shù)關(guān)系。所以︱r︱越大,表示兩者關(guān)系越密切[15]。當(dāng)然,計算出的相關(guān)系數(shù)是否顯著,需要經(jīng)過顯著性檢驗[14]。

    2 全域空間型及成因分析

    通過對1968~2007年處理為月均的日均降水量資料采取標(biāo)準(zhǔn)化距平之后,進(jìn)行EOF分析,得出10個模態(tài)相應(yīng)的空間變化和時間變化特征。同時,利用North檢驗,發(fā)現(xiàn)EOF分析方法中前4個模態(tài)通過了檢驗,說明前4個模態(tài)對華南前汛期降水的貢獻(xiàn)是較大的,方差貢獻(xiàn)率分別為26.658%、10.783%、7.452%和5.119%,且它們的累積方差貢獻(xiàn)率達(dá)50.012%,能很好地代表華南前汛期降水的空間分布特征,因此這里給出前4個模態(tài)對應(yīng)的特征場。

    2.1 第1模態(tài)從第1模態(tài)空間典型場和其時間序列(圖2a、b)可以看出,華南整體呈現(xiàn)出一致的態(tài)勢,其最大值中心位于福建省中南部和廣東東部地區(qū),第1模態(tài)的方差貢獻(xiàn)率為26.658%,方差貢獻(xiàn)較大,是華南前汛期降水的較為重要的一種分布型;第1模態(tài)對應(yīng)的時間序列中高值年與圖2a中的正值相乘的結(jié)果表示華南整體降水多,反之亦然(圖2b)。相對應(yīng)高低值年的500 hPa高度距平差值圖(圖2c)顯示,中高緯度的環(huán)流形勢大體為兩高一低型,低值中心位于巴爾喀什湖以西,高值中心分別位于中西伯利亞和阿留申地區(qū),低緯度環(huán)流形勢不明顯;相對應(yīng)高低值年的850 hPa假相當(dāng)位溫[16]差值圖(圖2d),可以很明顯地看到在整個華南區(qū)域內(nèi)均為正能量區(qū),有利于降水,與在高值年是全域一致型的降水偏多很好地吻合;從相對應(yīng)高低值年的850 hPa風(fēng)場差值圖(圖2e)可看到,副高中心位于海上,同時在華南北部有一切變線,在副高西北部、切變線東南部有利于降水,而華南正好處于這個區(qū)域,因此在高值年間降水偏多;從相對應(yīng)高低值年的850 hPa水汽通量及其散度差值合成圖(圖2f)可以看出,華南地區(qū)有從南海輸送來的水汽,且大部分區(qū)域處在水汽通量散度<0的陰影區(qū)中及水汽通量輻合,水汽條件有利于降水,也能說明華南此模態(tài)高值年整個區(qū)域降水多的情況,低值年的情況正好相反。

    2.2 第2模態(tài)從第2模態(tài)空間典型場(圖3a)可以看出,這一模態(tài)華南地區(qū)降水呈現(xiàn)的是西北和東南相反的空間結(jié)構(gòu)特征,其負(fù)值的最大值中心位于華南西北部即湖南西部,正值的最大值中心位于華南東南部即廣東東部福建南部,另外,第2模態(tài)的方差貢獻(xiàn)率為10.783%;其時間序列圖中高值年與空間模態(tài)正值區(qū)的乘積對應(yīng)的是降水多的年份、高值年與負(fù)值區(qū)的乘積對應(yīng)的是降水少的年份,反之亦然(圖3b)。從相對應(yīng)高低值年中的500 hPa高度距平差值(圖3c)可以看出主要的環(huán)流異常形勢,在中高緯度地區(qū),有一異常低壓帶,中心分別位于貝加爾湖以北和鄂霍次克海附近,低緯度環(huán)流形勢并不明顯。相對應(yīng)高低值年中的850 hPa風(fēng)場和散度差值合成圖(圖3d)顯示,在華南東南部有小片的負(fù)值區(qū)及輻合區(qū),易產(chǎn)生上升運動而引起降水,與高值年華南東南部降水偏多的情況相吻合,同時在華南的西北部為正值的最大中心區(qū)及輻散最明顯的區(qū)域,與高值年華南西北部降水偏少的結(jié)果也一致,因此綜合起來,可以很好地反映這一模態(tài)東西反向的情況;低值年的情況正好相反。從相對應(yīng)高低值年中的850 hPa水汽通量和水汽通量散度差值合成圖(圖3e)可以明顯看到,在華南東部有大片的負(fù)值區(qū),這里有較明顯的水汽輻合,最大輻合中心位于華南的東南部,同時從水汽通量可以看出,有水汽從南海輸送到華南,水汽條件較好,可見,東部較西部更利于形成降水,也在一定程度上說明了高值年的西北和東南相反且東部降水多的空間結(jié)構(gòu)特征,低值年的情況正好相反。

    2.3 第3模態(tài)從第3模態(tài)空間典型場(圖4a)可以清楚看出這一模態(tài)與前2個模態(tài)均有很大的差異,華南地區(qū)降水呈現(xiàn)的是南北相反的空間結(jié)構(gòu)特征,其正值的最大值中心位于江西省,負(fù)值的最大值中心位于廣東和廣西省南部,另外,第3模態(tài)的方差貢獻(xiàn)率為7.452%,通過了檢驗,說明這也是華南前汛期降水較為重要的一種模態(tài)。其時間序列圖中高值年與空間模態(tài)正值區(qū)的乘積對應(yīng)的是降水多的年份、高值年與負(fù)值區(qū)的乘積對應(yīng)的是降水少的年份,反之亦然(圖4b)。從相對應(yīng)高低值年中的500 hPa高度距平差值(圖4c)可以看出,在中高緯度,主要為兩低一高型,低值中心分別位于巴爾喀什湖以北、貝加爾湖以西地區(qū)和阿留申地區(qū),高值中心位于西西伯利亞地區(qū),低緯度環(huán)流形勢不太明顯。從相對應(yīng)高低值年中的850 hPa假相當(dāng)位溫差值圖(圖4d)可以看到,華南北部為正值區(qū),負(fù)值區(qū)位于南部,即北部為高能量區(qū),容易引起對流不穩(wěn)定,有利于形成降水,而南部的情況則剛好

    注:a.為空間典型場;b為時間序列;c為相對應(yīng)高低值年的500 hPa高度距平差值,單位為gpm;d為相對應(yīng)高低值年的850 hPa位溫差值場,單位為℃;e為相對應(yīng)高低值年的850 hPa風(fēng)場差值,單位為m/s;f為相對應(yīng)高低值年850 hPa的水汽通量及其散度差值合成圖,單位分別為g/(cm·hPa·s)、10-6 g/(cm2·hPa·s),陰影區(qū)代表負(fù)值區(qū)。圖2 EOF的第1模態(tài)對應(yīng)的特征場

    注:a為空間典型場;b為時間序列;c為相對應(yīng)高低值年的500 hPa高度距平差值,單位為gpm;d為相對應(yīng)高低值年的850 hPa風(fēng)場和散度差值合成圖,單位分別為m/s、10-6s-1;e為相對應(yīng)高低值年850 hPa的水汽通量及其散度差值合成圖,單位分別為g/(cm·hPa·s)、10-6 g/(cm2·hPa·s)。圖d、e中陰影區(qū)代表負(fù)值區(qū)。圖3 EOF的第2模態(tài)對應(yīng)的特征場

    注:a.為空間典型場;b為時間序列;c為相對應(yīng)高低值年的500 hPa高度距平差值,單位為gpm;d為相對應(yīng)高低值年的850 hPa位溫差值場,單位為℃;e為相對應(yīng)高低值年的850 hPa風(fēng)場和散度差值合成圖,單位分別為m/s、10-6s-1;f為相對應(yīng)高低值年850 hPa的水汽通量及其散度差值合成圖,單位分別為g/(cm·hPa·s)、10-6 g/(cm2·hPa·s),陰影區(qū)代表負(fù)值區(qū)。圖4 EOF的第3模態(tài)對應(yīng)的特征場

    相反,很好地解釋了高值年間這一模態(tài)北部降水多、南部偏少的情況。由相對應(yīng)高低值年中的850 hPa風(fēng)場和散度差值(圖4e)可見,北部多為負(fù)值區(qū),即輻合區(qū),有上升運動,利于降水形成,而南部正好相反,為輻散區(qū),有下沉運動,與在高值年北多南少的情況正好相符。從相對應(yīng)高低值年中的850 hPa水汽通量與水汽通量散度差值合成圖(圖4f)可以看到,在華南有從南海輸送而來的水汽,且在北部有負(fù)的水汽通量散度,即代表水汽通量輻合,水汽條件較好,利于降水,能很好地說明這一模態(tài)高值年的分布特征;低值年的情況則正好相反。

    2.4 第4模態(tài)從第4模態(tài)空間典型場(圖5a)可以看出,華南地區(qū)呈現(xiàn)的是中間區(qū)域變化趨勢一致,其余地方變化趨勢相反的空間結(jié)構(gòu)特征,其正值的最大值中心位于貴州地區(qū)中北部附近,負(fù)值的最大值中心位于湖南南部、廣東西北部地區(qū)附近,且第4模態(tài)的方差貢獻(xiàn)率為5.119%;其時間序列中高值年與空間模態(tài)正值區(qū)的乘積對應(yīng)的是降水多的年份、高值年與空間模態(tài)負(fù)值區(qū)的乘積對應(yīng)的是降水少的年份,反之亦然(圖5b)。從相對應(yīng)高低值年中的500 hPa高度距平差值(圖5c)可以看出,在中高緯度,主要為兩低一高型,低值中心分別位于貝加爾湖以西地區(qū)和白令海以南地區(qū),高值中心位于我國東北及其以北地區(qū),低緯度環(huán)流形勢不太明顯。850 hPa多年的風(fēng)場平均狀況(圖5d)顯示,低緯度盛行東風(fēng),風(fēng)速較大,且副高脊線穩(wěn)定在21°N左右,華南地區(qū)處于西南氣流中即太平洋副熱帶高壓的西北部,同時華南中部為較強(qiáng)的偏西氣流,有南海和孟加拉灣的水汽輸送,水汽條件好,綜合起來有利于降水。相對應(yīng)高低值年的850 hPa的風(fēng)場差值圖(圖5e)與圖5d作對比后發(fā)現(xiàn)有明顯的不同,在華南中部為異常的偏東氣流,且風(fēng)速較平均年小,即較平均年西南氣流減小,水汽通量輸送明顯減少,水汽條件不利于降水,說明高值年華南中部降水偏少的空間分布特征。此外,其他變量所反映的情況不具有很好的代表性,故在這一模態(tài)不再進(jìn)一步討論。

    注:a為空間典型場;b為時間序列;c為相對應(yīng)高低值年的500 hPa高度距平差值,單位為gpm;d為850 hPa風(fēng)場平均,單位為m/s;e為相對應(yīng)高低值年850 hPa風(fēng)場差值,單位為m/s。圖5 EOF的第4模態(tài)對應(yīng)的特征場

    以上分別敘述了EOF前4個模態(tài)的空間和時間結(jié)構(gòu)特征以及成因分析,與郭其蘊(yùn)等的時空結(jié)構(gòu)特征分析結(jié)果相似,前3個模態(tài)的分布型具有很好的一致性,但其僅給出了前3個模態(tài),總方差貢獻(xiàn)率為64.2%[17]。

    在上面討論的基礎(chǔ)上,對比通過檢驗的前4個模態(tài)的EOF方差貢獻(xiàn)率及其對應(yīng)的降水情況(表1)和通過檢驗的前4個模態(tài)的EOF方差貢獻(xiàn)率及相應(yīng)的去除不同模態(tài)相同年份后的降水情況(表2),可以很直觀清晰地看到各個空間模態(tài)的貢獻(xiàn)率大小和旱澇年的分布情況,且發(fā)現(xiàn)EOF中有許多其他的因子會影響旱澇年的選取,因此這里剔除了4個時間序列中存在2個以上的年份,從而使結(jié)果更加準(zhǔn)確。

    3 局域空間型分析

    3.1 REOF各個模態(tài)分析通過EOF方法的討論可以看出,華南前汛期降水的空間格局既有全局一致的多雨或少雨現(xiàn)象,也有南北、東西和中間存在的內(nèi)部的較大差異,但其仍以緯向分布型為主,不能更為精細(xì)地描述不同地理區(qū)域的特征,因此在EOF分析的基礎(chǔ)上,在此通過旋轉(zhuǎn)EOF將華南地區(qū)通過載荷向量值最大的地區(qū)來分區(qū),這樣可以得出比較精細(xì)的地理分區(qū)特征。

    利用選取的10個模態(tài)的REOF分析進(jìn)行檢驗后,發(fā)現(xiàn)前5個模態(tài)均通過了檢驗。因為REOF方法可以更進(jìn)一步地反映區(qū)域性的空間結(jié)構(gòu)特征,所以在每個模態(tài)的空間典型場中僅選取絕對值最大的高載荷區(qū),這樣可以更加直觀地反映出華南前汛期的區(qū)域性降水特征。從華南地區(qū)前5個模態(tài)的分區(qū)示意圖(圖6)可以很明顯地看出,華南地區(qū)被很好地分為華南東南部、華南中部、華南東北部、華南西南部和華南西北部5個區(qū)域,對應(yīng)的具體模態(tài)為第1模態(tài)的高載荷區(qū)位于福建廣東地區(qū),其方差貢獻(xiàn)率為20.281%;第2模態(tài)的高載荷區(qū)位于湖南南部和廣東、廣西北部地區(qū),其方差貢獻(xiàn)率為9.300%;第3模態(tài)的高載荷區(qū)位于江西和湖南地區(qū),其方差貢獻(xiàn)率為8.185%;第4模態(tài)的高載荷區(qū)位于廣西南部地區(qū),其方差貢獻(xiàn)率為7.974%;第5模態(tài)的高載荷區(qū)位于貴州地區(qū),其方差貢獻(xiàn)率為5.187%;前5個模態(tài)的累積方差貢獻(xiàn)率為50.927%,累積方差貢獻(xiàn)較大。因此,在研究華南前汛期降水氣候特征及其異常變化時,上述5個類型區(qū)域的降水特征有一定的重要性。對比通過檢驗的前5個模態(tài)的REOF的方差貢獻(xiàn)率及對應(yīng)的降水情況(表3)和通過檢驗的前5個模態(tài)的REOF的方差貢獻(xiàn)率及相應(yīng)的去除不同模態(tài)相同年份后的降水情況(表4),可以很直觀清晰地看出各個空間模態(tài)的貢獻(xiàn)大小,且發(fā)現(xiàn)REOF中有很多其他的因子會影響旱澇年的選取,因此這里也做了剔除5個時間序列中存在2個以上的年份,從而使結(jié)果更加準(zhǔn)確。

    3.2 子域的環(huán)流對比和相關(guān)性分析對比EOF所得到的華南分型模態(tài)發(fā)現(xiàn),REOF方法反映的空間典型場中的第2模態(tài)為華南中部型,且其方差貢獻(xiàn)率為9.300%,是較為重要且有一定代表性的分型,而在EOF中,中間與周圍反向型為第4模態(tài),其方差貢獻(xiàn)為5.119%,較其他模態(tài)并不是有很大的代表性,因此為了進(jìn)一步研究兩者的區(qū)別,對REOF的第2模態(tài)做了進(jìn)一步的研究。

    結(jié)合(圖6)華南地區(qū)前5個模態(tài)的分區(qū)示意圖發(fā)現(xiàn),高載荷區(qū)位于廣西東北、湖南南部地區(qū),對應(yīng)時間序列圖中(圖7b),高值年與空間模態(tài)正值區(qū)的乘積對應(yīng)的是降水多的年份、高值年與空間模態(tài)負(fù)值區(qū)的乘積對應(yīng)的是降水少的年份;反之亦然。此模態(tài)反映的是華南中部地區(qū)40年中的降水情況。但實際上此模態(tài)的中心為負(fù)值(圖7a),圖6已做了絕對值處理,因此在高值年中華南基本上降水偏少,處于干旱年。配合REOF這一模態(tài)時間序列所對應(yīng)的高低值年的850 hPa風(fēng)場與散度差值合成圖(圖7c)可以看到,在華南除了東南、西南角有小范圍的負(fù)值區(qū)外,大部分區(qū)域均處于正值區(qū),正值為輻散區(qū),且華南中部的數(shù)值最大最為明顯,輻散有下沉運動不利于降水,降水易偏少,很好地反映了這一模態(tài)高值年的情況;低值年的情況正好相反。從相對應(yīng)的高低值年的850 hPa水汽通量及其散度差值合成(圖7d)也可以清楚看到,除了華南東西兩側(cè)有小范圍的負(fù)值區(qū)外,其余均為正值區(qū)且中部數(shù)值最大,即華南大部分區(qū)域水汽通量散度輻散,中部最為明顯,水汽條件差,不利于降水,也能很好地說明高值年的情況;低值年的情況也正好相反。通過分析發(fā)現(xiàn),這個模態(tài)與EOF的第4模態(tài)存在較大的環(huán)流差別。

    為了進(jìn)一步分析REOF的這一模態(tài)的降水特征,用逐月的南方濤動指數(shù)(SOI)[18]與時間序列做相關(guān),發(fā)現(xiàn)2月份的南方濤動指數(shù)與這一模態(tài)有較好的正相關(guān)性,相關(guān)系數(shù)為0.284 5,通過了90%的信度檢驗,證明這兩者有一定的相關(guān)性,且這種相關(guān)是3個月的滯后性相關(guān);但此模態(tài)為高值年降水偏少年,即南方濤動指數(shù)為正值時,華南中部降水偏少,反之亦然。此外,還研究了Nino3.4指數(shù)與華南前汛期的降水相關(guān)性(圖8),在實際工作中,人們發(fā)現(xiàn)用Nino 3區(qū)和Nino 4區(qū)的平均海溫距平來進(jìn)行研究比較有代表性,因此常采用后來增加的新區(qū)Nino 3.4區(qū)(5°N~5°S、170°~120°W)來進(jìn)行分析[18]。結(jié)果表明,華南前汛期降水與1月份的Nino 3.4指數(shù)有很好的相關(guān)性,相關(guān)性最大的區(qū)域大致位于華南中部區(qū)域,與這一模態(tài)的分布情況類似,基本上均通過了90%的信度檢驗,甚至是95%的信度檢驗,證明這兩者有很好的正相關(guān)關(guān)系,且是4個月的滯后性相關(guān),即Nino 3.4指數(shù)為正值時,華南中部降水增多。研究發(fā)現(xiàn),當(dāng)Nino 3.4指數(shù)為正值時,代表中東太平洋海溫偏高,為厄爾尼諾年,形成在西太平洋下沉、東太平洋上升的異常環(huán)流,并引發(fā)在赤道附近下沉、華南區(qū)域上升的經(jīng)向環(huán)流,容易引起降水,即正相關(guān)關(guān)系[18]。

    表1 前4個模態(tài)EOF的方差貢獻(xiàn)率及對應(yīng)的降水情況

    表2 前4個模態(tài)EOF的方差貢獻(xiàn)率及相應(yīng)的去除不同模態(tài)相同年份后的降水情況

    表3 前5個模態(tài)REOF的方差貢獻(xiàn)率及對應(yīng)的降水情況

    表4 前5個模態(tài)REOF的方差貢獻(xiàn)率及相應(yīng)的去除不同模態(tài)的相同年份后的降水情況

    圖6 華南地區(qū)前5個模態(tài)的分區(qū)示意圖

    注:a為空間模態(tài);b為時間序列;c為時間序列相對應(yīng)的高低值年的850 hPa風(fēng)場和散度差值合成圖,單位分別為m/s 、10-6s-1;d為850 hPa水汽通量及其散度差值合成圖,單位分別為g/(cm·hPa·s) 、10-6g/(cm2·hPa·s)。圖c、d中陰影區(qū)表示負(fù)值區(qū)。圖7 REOF的第2模態(tài)對應(yīng)的特征場

    注:空心標(biāo)識為負(fù)相關(guān),實心為正相關(guān),且從大到小分別為通過95%、90%信度檢驗,圈內(nèi)為大片通過檢驗的明顯區(qū)域。圖8 1月份Nino3.4指數(shù)與華南前汛期降水的相關(guān)分布

    4 結(jié)論

    (1)通過EOF方法分析了華南前汛期的降水典型分布,前4個特征向量累積方差貢獻(xiàn)率高達(dá)50.012%,其中第1典型場說明華南前汛期全區(qū)降水是一致的,即在不同年份,華南前汛期降水一致偏多或偏少;第2典型場說明是東西反向型的,即在高值年,華南西部少東部多,反之亦然;第3典型場說明是南北反向型的,即在高值年,華南北部多南部少,反之亦然;第4典型場說明是中間型的,高值年華南中部雨少,反之亦然。

    (2)通過對位勢場、風(fēng)場、散度場、假相當(dāng)位溫場、水汽通量場和水汽通量散度場等方面的討論,取時間序列對應(yīng)的高低值年的差值進(jìn)行分析,發(fā)現(xiàn)環(huán)流形勢均能很好地說明EOF前4個模態(tài)的分布特征。

    (3)通過REOF方法研究了華南前汛期地區(qū)內(nèi)部的分型情況,華南降水區(qū)大致可以分為5個高載荷區(qū)域,即華南東南、華南中部、華南東北部、華南西南部和華南西北部,且第2模態(tài)與EOF第4模態(tài)類似,也能從環(huán)流形勢很好地解釋這一模態(tài)的分布情況。

    (4)從相關(guān)性研究發(fā)現(xiàn),2月份的南方濤動指數(shù)與REOF第2模態(tài)時間序列存在3個月的滯后性正相關(guān)關(guān)系,但此模態(tài)在高值年為華南中部區(qū)域降水偏少,即表明當(dāng)2月南方濤動指數(shù)為正值時,前汛期華南中部降水偏少,低值年相反。此外,1月份的Nino 3.4指數(shù)與華南前汛期降水也有正相關(guān)關(guān)系,為4個月的滯后性相關(guān),相關(guān)性較大的區(qū)域在華南中部,也與這一模態(tài)的分布類似,即當(dāng)Nino 3.4指數(shù)為正值時,華南中部降水偏多,反之亦然。

    [1] 李江南,王安宇.廣東省前汛期和后汛期降水的氣候特征[J].中山大學(xué)學(xué)報:自然科學(xué)版,2002,41(3):91-96.

    [2] FISCHER T,GEMMER M,LIU L L,et al.Trends in Monthly Temperature and Precipitation Extremes in the Zhujiang River Basin,South China(1961-2007)[J].Advances in Climate Change Research,2010,1(2):63-70.

    [3] CHEN H,LI H Y,LIN Z H.Interdecadal variability of spring precipitation over South China and its associated atmospheric water vapot transport[J].Armospheric and Oceanic Science Letters,2009,2(2):113-118.

    [4] 鄭彬,谷德軍,李春暉,等.華南前汛期的鋒面降水和夏季風(fēng)降水II空間分布特征[J].大氣科學(xué),2007,31(3):495-504.

    [5] ZHENG B,LIANG J Y,LIN A L,et al.Frontal rainfall and summer monsoon rainfall during pre-rainy season in South China.Part I:Determination of the division dates[J].ChineseJ Atmos Sci,2006,30:1207-1216.

    [6] 池艷珍,何金海,吳志偉.華南前汛期不同降水時段的特征分析[J].南京氣象學(xué)院學(xué)報,2005,28(2):163-171.

    [7] 林愛蘭,吳尚森.近40 年華南汛期旱澇變化及趨勢預(yù)測[J].熱帶氣象學(xué)報,1996,12(2):160-166.

    [8] 鄧立平,王謙謙.華南前汛期(4~6月)降水異常特征及其與我國近海海溫的關(guān)系[J].熱帶氣象學(xué)報,2002,18(1):45-55.

    [9] 馬慧,王謙謙,陳楨華.華南前汛期降水異常的時空變化特征[J].高原氣象,2006,25(2):325-329.

    [10] 屠其璞,王俊德,丁裕國,等.氣象應(yīng)用概率統(tǒng)計學(xué)[M].北京:氣象出版社,1984:445.

    [11] 章基嘉.中長期天氣預(yù)報基礎(chǔ)[M].北京:氣象出版社,1994:66-108.

    [12] 黃嘉佑.氣象統(tǒng)計分析與預(yù)報方法[M].北京:氣象出版社,2000:135-139.

    [13] 么枕生.載荷相關(guān)模式用于氣候分類與天氣氣候描述[M].北京:氣象出版社,1998:1-9.

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

    [15] 施能.氣象科研與預(yù)報中的多元分析方法[M].北京:氣象出版社,2008:15.

    [16] 陳創(chuàng)買,郭英瓊,劉麗英.論假相當(dāng)位溫的數(shù)值計算(一)[J].中山大學(xué)學(xué)報,1989,28(2):111-116.

    [17] 郭其蘊(yùn),沙萬英.華南前汛期降水變率的分析[J].應(yīng)用氣象學(xué)報,1998(9):9-15.

    [18] 孫照渤,陳海山,譚桂榮,等.短期氣候預(yù)測基礎(chǔ)[M].北京:氣象出版社,2010:94-99.

    Study on Climatic Features of Pre-Summer Flood Period Precipitation in South China

    LIU Hui, LI Yi-ming, FAN Hua-liang et al

    (Suzhou Lightning Protection Center, Suzhou, Jiangsu 215000)

    The rainfall data of 94 stations in south China from 1968-2007 are analyzed by the way of EOF and REOF. According to the analysis results, the former four typical fields from the EOF analysis can successfully reflect the anomalous spatial distribution of rainfall in South China, to be specific, consistent as a whole, disparate in the south and the north, disparate in the east and the west, disparate in the central and the peripheral areas. Through the analysis of the NCEP/NCAR data, it can be found that different modals are caused by different circulation situation. The former five by the REOF rotated spatial patterns can make climate zoning in pre-summer flood period precipitation of south China: divided into southeast, central, northeast, southwest and southnorth. And the circulation situation of the second mode which is similar from the forth mode of EOF was analyzed. This mode has a positive correlation with soi index in February and Nino3.4 index in January.

    South China; Precipitation in pre-summer flood period; Climatic characteristics; EOF; REOF

    劉輝(1985- ),男,江蘇蘇州人,助理工程師,從事雷電科學(xué)與預(yù)警研究。

    2014-12-15

    S 161.6

    A

    0517-6611(2015)04-217-08

    猜你喜歡
    散度華南高值
    帶勢加權(quán)散度形式的Grushin型退化橢圓算子的Dirichlet特征值的上下界
    養(yǎng)殖廢棄物快速發(fā)酵及高值轉(zhuǎn)化土壤修復(fù)生物肥料關(guān)鍵技術(shù)
    麻文化發(fā)展與高值利用前景展望
    具有部分BMO系數(shù)的非散度型拋物方程的Lorentz估計
    華南風(fēng)采
    海峽姐妹(2019年8期)2019-09-03 01:00:46
    記華南女院前三任校長
    海峽姐妹(2018年10期)2018-12-26 01:20:56
    華南掠影
    海峽姐妹(2018年10期)2018-12-26 01:20:54
    蘇萌娜 初心不渝 情牽華南
    海峽姐妹(2018年10期)2018-12-26 01:20:52
    H型群上一類散度形算子的特征值估計
    H?rmander 向量場上散度型拋物方程弱解的Orlicz估計
    床上黄色一级片| 九九久久精品国产亚洲av麻豆| 哪里可以看免费的av片| 亚洲五月天丁香| 成人二区视频| 亚洲无线观看免费| 午夜福利18| 哪里可以看免费的av片| 少妇人妻精品综合一区二区 | 一卡2卡三卡四卡精品乱码亚洲| 大又大粗又爽又黄少妇毛片口| 日韩成人伦理影院| 日韩欧美国产在线观看| 中文字幕精品亚洲无线码一区| 日本撒尿小便嘘嘘汇集6| 麻豆av噜噜一区二区三区| 一进一出抽搐动态| 老司机影院成人| 亚洲av免费高清在线观看| 51国产日韩欧美| 夜夜夜夜夜久久久久| 亚洲欧美精品综合久久99| 九九爱精品视频在线观看| 又爽又黄a免费视频| 免费人成在线观看视频色| 欧美激情国产日韩精品一区| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲av香蕉五月| 乱人视频在线观看| 成人鲁丝片一二三区免费| 久久久久久久久中文| 91久久精品电影网| 成人特级av手机在线观看| 简卡轻食公司| 亚洲人成网站在线观看播放| 老熟妇仑乱视频hdxx| 免费高清视频大片| 嫩草影院新地址| 给我免费播放毛片高清在线观看| 三级经典国产精品| 蜜臀久久99精品久久宅男| 亚洲综合色惰| a级毛片免费高清观看在线播放| 舔av片在线| 免费人成视频x8x8入口观看| 欧美成人一区二区免费高清观看| 久久久色成人| 亚洲欧美清纯卡通| 欧美丝袜亚洲另类| 毛片一级片免费看久久久久| 欧美色视频一区免费| 成年女人毛片免费观看观看9| 亚洲美女搞黄在线观看 | 日韩欧美 国产精品| 91av网一区二区| 亚洲精品国产成人久久av| 亚洲五月天丁香| 国产精品一区www在线观看| 成熟少妇高潮喷水视频| 联通29元200g的流量卡| 午夜福利18| avwww免费| 国产精品爽爽va在线观看网站| а√天堂www在线а√下载| 午夜激情欧美在线| 日韩大尺度精品在线看网址| 亚洲成人av在线免费| 深夜a级毛片| a级毛片a级免费在线| 亚洲va在线va天堂va国产| 久久草成人影院| 欧美另类亚洲清纯唯美| 天堂网av新在线| 内射极品少妇av片p| 你懂的网址亚洲精品在线观看 | 久久精品国产亚洲网站| 久久精品国产鲁丝片午夜精品| 国产亚洲精品av在线| 成人综合一区亚洲| 日本熟妇午夜| 成年女人看的毛片在线观看| 久久国产乱子免费精品| avwww免费| 美女xxoo啪啪120秒动态图| 一边摸一边抽搐一进一小说| 亚洲成a人片在线一区二区| 老司机午夜福利在线观看视频| 欧美色欧美亚洲另类二区| aaaaa片日本免费| 欧美成人精品欧美一级黄| 久久久久久久亚洲中文字幕| 国产免费男女视频| 国产爱豆传媒在线观看| www日本黄色视频网| 久久久久久久久久黄片| 男插女下体视频免费在线播放| 国产aⅴ精品一区二区三区波| 日本黄色视频三级网站网址| 国产一区二区三区在线臀色熟女| 蜜桃久久精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 精品乱码久久久久久99久播| 国产日本99.免费观看| 亚洲性久久影院| 久久国产乱子免费精品| 久久九九热精品免费| 免费人成在线观看视频色| 波多野结衣高清无吗| 欧美性猛交╳xxx乱大交人| 欧美中文日本在线观看视频| 国产高清不卡午夜福利| 亚洲精品国产成人久久av| 我要搜黄色片| 午夜福利成人在线免费观看| 免费在线观看影片大全网站| 国产黄色小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| ponron亚洲| 黄色视频,在线免费观看| 欧美中文日本在线观看视频| .国产精品久久| 亚洲av美国av| 最近在线观看免费完整版| 国产色婷婷99| .国产精品久久| 亚洲中文字幕一区二区三区有码在线看| 18禁在线无遮挡免费观看视频 | 色尼玛亚洲综合影院| 国产欧美日韩精品一区二区| 特级一级黄色大片| 成人毛片a级毛片在线播放| 日本 av在线| 中国国产av一级| 亚洲精品久久国产高清桃花| 久久草成人影院| 国产在线男女| 亚洲中文字幕一区二区三区有码在线看| 3wmmmm亚洲av在线观看| a级毛片免费高清观看在线播放| 麻豆成人午夜福利视频| aaaaa片日本免费| 一级黄片播放器| 亚洲内射少妇av| 成人漫画全彩无遮挡| 午夜激情福利司机影院| 免费观看精品视频网站| av福利片在线观看| 国产高清激情床上av| 欧美性猛交╳xxx乱大交人| 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 97人妻精品一区二区三区麻豆| 欧美一区二区亚洲| 成年女人毛片免费观看观看9| 国产探花极品一区二区| 欧美又色又爽又黄视频| 99在线人妻在线中文字幕| 国产淫片久久久久久久久| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 日本精品一区二区三区蜜桃| 老司机影院成人| 日本在线视频免费播放| 看免费成人av毛片| 少妇人妻一区二区三区视频| 日韩欧美免费精品| 亚洲成人久久爱视频| 高清毛片免费看| 国产成人一区二区在线| 国产高清不卡午夜福利| 深爱激情五月婷婷| 国产在线精品亚洲第一网站| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 精品少妇黑人巨大在线播放 | 日本黄大片高清| 俄罗斯特黄特色一大片| 一区二区三区高清视频在线| 少妇裸体淫交视频免费看高清| 高清毛片免费看| 九九在线视频观看精品| 最近中文字幕高清免费大全6| 国产高清三级在线| 精品久久久久久久久亚洲| 国内精品宾馆在线| 国产精品亚洲美女久久久| 国产精品久久久久久精品电影| 国产美女午夜福利| 婷婷精品国产亚洲av| 有码 亚洲区| 欧美成人精品欧美一级黄| h日本视频在线播放| 久久精品影院6| 欧美成人精品欧美一级黄| 国内精品美女久久久久久| 国产精品精品国产色婷婷| 久久久久久久久久久丰满| 日本免费一区二区三区高清不卡| 成人无遮挡网站| 亚洲中文日韩欧美视频| 日韩欧美三级三区| 老司机福利观看| 日本黄大片高清| 日本爱情动作片www.在线观看 | 激情 狠狠 欧美| 亚洲图色成人| 亚洲乱码一区二区免费版| 在线观看午夜福利视频| 亚洲欧美日韩东京热| 国产成年人精品一区二区| 成人二区视频| 日韩精品青青久久久久久| 亚洲婷婷狠狠爱综合网| 伊人久久精品亚洲午夜| 亚洲无线在线观看| 啦啦啦观看免费观看视频高清| 少妇人妻一区二区三区视频| 最近的中文字幕免费完整| 国语自产精品视频在线第100页| 99在线人妻在线中文字幕| 午夜精品国产一区二区电影 | 搡女人真爽免费视频火全软件 | 麻豆久久精品国产亚洲av| 日本a在线网址| 精品国内亚洲2022精品成人| 久久婷婷人人爽人人干人人爱| 国产麻豆成人av免费视频| 人人妻人人澡人人爽人人夜夜 | 国产欧美日韩精品一区二区| 精品乱码久久久久久99久播| 久久精品国产亚洲av天美| 91午夜精品亚洲一区二区三区| 午夜福利在线观看吧| 亚洲无线观看免费| 十八禁网站免费在线| 如何舔出高潮| 最好的美女福利视频网| 成人性生交大片免费视频hd| 在线观看66精品国产| 乱人视频在线观看| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片免费观看直播| 国产片特级美女逼逼视频| 久久久久国产精品人妻aⅴ院| 久久久久久国产a免费观看| 别揉我奶头 嗯啊视频| av天堂在线播放| 亚洲精品久久国产高清桃花| 伦理电影大哥的女人| 中文字幕熟女人妻在线| 99久久精品一区二区三区| 男女做爰动态图高潮gif福利片| 中文字幕精品亚洲无线码一区| 全区人妻精品视频| 久久韩国三级中文字幕| 免费看光身美女| 国产一区二区激情短视频| 天堂影院成人在线观看| 日韩精品有码人妻一区| 麻豆av噜噜一区二区三区| 不卡视频在线观看欧美| 久久久久国产精品人妻aⅴ院| 俺也久久电影网| 精品99又大又爽又粗少妇毛片| 亚洲专区国产一区二区| 日日干狠狠操夜夜爽| 亚洲中文日韩欧美视频| 色哟哟·www| 国产又黄又爽又无遮挡在线| 一a级毛片在线观看| 一级黄片播放器| av在线蜜桃| 性插视频无遮挡在线免费观看| 国产午夜精品论理片| 有码 亚洲区| 国产综合懂色| 国产成人91sexporn| 中出人妻视频一区二区| 欧美日韩在线观看h| 国产美女午夜福利| 在线天堂最新版资源| 久久6这里有精品| 午夜免费男女啪啪视频观看 | 真实男女啪啪啪动态图| 日韩中字成人| 性欧美人与动物交配| 22中文网久久字幕| 人妻少妇偷人精品九色| 夜夜看夜夜爽夜夜摸| 99久久无色码亚洲精品果冻| 国产毛片a区久久久久| 欧美性感艳星| 国产精品不卡视频一区二区| 亚洲无线在线观看| 国产精品野战在线观看| 国产精品三级大全| 真实男女啪啪啪动态图| 精品乱码久久久久久99久播| 国产 一区 欧美 日韩| 少妇丰满av| 桃色一区二区三区在线观看| 午夜福利在线在线| 非洲黑人性xxxx精品又粗又长| 免费在线观看影片大全网站| 99在线视频只有这里精品首页| 亚洲成人av在线免费| 丝袜美腿在线中文| 国产精品国产高清国产av| 免费看美女性在线毛片视频| 国产伦一二天堂av在线观看| 免费观看人在逋| 亚洲人成网站在线观看播放| 欧美性猛交黑人性爽| 天堂网av新在线| 99久久中文字幕三级久久日本| 女的被弄到高潮叫床怎么办| 五月伊人婷婷丁香| 老司机福利观看| 欧美日韩在线观看h| 麻豆乱淫一区二区| 97碰自拍视频| 一级黄片播放器| 国产爱豆传媒在线观看| 精品人妻熟女av久视频| 亚洲激情五月婷婷啪啪| 免费在线观看成人毛片| 蜜臀久久99精品久久宅男| 国产麻豆成人av免费视频| 1000部很黄的大片| 成人性生交大片免费视频hd| 天天躁夜夜躁狠狠久久av| 99在线视频只有这里精品首页| 久久精品91蜜桃| videossex国产| 三级经典国产精品| 一区二区三区免费毛片| 欧美日韩综合久久久久久| 久久韩国三级中文字幕| 亚洲专区国产一区二区| 天天躁日日操中文字幕| 亚洲专区国产一区二区| 久久久久九九精品影院| 黄色视频,在线免费观看| 国产av麻豆久久久久久久| 亚洲av成人av| a级一级毛片免费在线观看| 精品乱码久久久久久99久播| 亚洲国产精品国产精品| 久久精品国产清高在天天线| 丝袜喷水一区| 国产精品一区二区三区四区久久| 成人一区二区视频在线观看| 精品久久久久久久久av| 国产av麻豆久久久久久久| 国产国拍精品亚洲av在线观看| 成人综合一区亚洲| 成年女人永久免费观看视频| 日韩欧美在线乱码| 国产精品精品国产色婷婷| 赤兔流量卡办理| a级毛片免费高清观看在线播放| 看片在线看免费视频| 国国产精品蜜臀av免费| 久久亚洲国产成人精品v| 乱人视频在线观看| 日韩国内少妇激情av| 淫秽高清视频在线观看| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 久久鲁丝午夜福利片| 嫩草影院新地址| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人| 白带黄色成豆腐渣| 日韩高清综合在线| 国产免费一级a男人的天堂| 老司机影院成人| 国产高潮美女av| 成人av在线播放网站| 久久热精品热| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 国产单亲对白刺激| 亚洲欧美成人综合另类久久久 | 国产v大片淫在线免费观看| 成人三级黄色视频| 亚洲精品乱码久久久v下载方式| 久久精品国产亚洲av香蕉五月| 国产乱人偷精品视频| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 国产欧美日韩精品亚洲av| 啦啦啦韩国在线观看视频| 哪里可以看免费的av片| 91在线观看av| av视频在线观看入口| 网址你懂的国产日韩在线| 国产精品久久视频播放| 国产午夜福利久久久久久| 麻豆成人午夜福利视频| 99国产极品粉嫩在线观看| 亚洲一区高清亚洲精品| 精品国产三级普通话版| av在线老鸭窝| 国产国拍精品亚洲av在线观看| 日日摸夜夜添夜夜添小说| 亚洲无线在线观看| 国产探花极品一区二区| 午夜福利在线观看吧| 亚洲国产高清在线一区二区三| 老司机影院成人| 22中文网久久字幕| 啦啦啦啦在线视频资源| 最新中文字幕久久久久| 十八禁国产超污无遮挡网站| 麻豆国产97在线/欧美| 国产精品人妻久久久久久| 超碰av人人做人人爽久久| 欧美激情国产日韩精品一区| 嫩草影院入口| 三级国产精品欧美在线观看| 51国产日韩欧美| 自拍偷自拍亚洲精品老妇| 老司机午夜福利在线观看视频| 一区福利在线观看| 国产 一区 欧美 日韩| 国产av不卡久久| 看黄色毛片网站| 久久综合国产亚洲精品| 亚洲av二区三区四区| 日韩欧美精品免费久久| 国产免费一级a男人的天堂| 美女 人体艺术 gogo| 国产精品福利在线免费观看| 亚洲不卡免费看| 国产精品,欧美在线| 青春草视频在线免费观看| or卡值多少钱| 白带黄色成豆腐渣| 久久午夜亚洲精品久久| 日本五十路高清| 男人狂女人下面高潮的视频| 国产精品av视频在线免费观看| 真实男女啪啪啪动态图| 91久久精品国产一区二区三区| 特级一级黄色大片| 一边摸一边抽搐一进一小说| 欧美一区二区精品小视频在线| 99热这里只有是精品在线观看| 国产一区二区在线观看日韩| 观看免费一级毛片| a级毛片a级免费在线| 亚洲最大成人中文| 亚洲精品成人久久久久久| 欧美在线一区亚洲| 久久亚洲国产成人精品v| 国产黄a三级三级三级人| 一本精品99久久精品77| 色噜噜av男人的天堂激情| 女的被弄到高潮叫床怎么办| 国产真实伦视频高清在线观看| 韩国av在线不卡| 国内少妇人妻偷人精品xxx网站| 99久久精品国产国产毛片| av天堂中文字幕网| 免费搜索国产男女视频| 伊人久久精品亚洲午夜| 乱人视频在线观看| 久久久久久久久中文| 麻豆国产97在线/欧美| 男女那种视频在线观看| 日日撸夜夜添| 亚洲第一电影网av| 久久婷婷人人爽人人干人人爱| 亚洲久久久久久中文字幕| 日韩中字成人| 亚洲久久久久久中文字幕| 少妇丰满av| 成年av动漫网址| 国产精品国产三级国产av玫瑰| 天堂av国产一区二区熟女人妻| 少妇裸体淫交视频免费看高清| 在线a可以看的网站| 人妻少妇偷人精品九色| www日本黄色视频网| 欧美高清性xxxxhd video| 日韩欧美精品v在线| 成年av动漫网址| 国产男人的电影天堂91| 国产伦一二天堂av在线观看| 插阴视频在线观看视频| 99国产极品粉嫩在线观看| 九九爱精品视频在线观看| 久久精品影院6| 欧美激情国产日韩精品一区| 亚洲欧美清纯卡通| 国产男靠女视频免费网站| 亚洲不卡免费看| 亚洲内射少妇av| 日韩,欧美,国产一区二区三区 | 日日干狠狠操夜夜爽| 久久精品91蜜桃| 亚洲成a人片在线一区二区| 久久99热这里只有精品18| 国产精品野战在线观看| 91久久精品国产一区二区成人| 国产免费一级a男人的天堂| 国产单亲对白刺激| 国产一区亚洲一区在线观看| 特级一级黄色大片| 99久国产av精品国产电影| 亚洲美女黄片视频| 十八禁国产超污无遮挡网站| 免费看美女性在线毛片视频| 国产精品久久久久久久电影| 国产亚洲欧美98| 国产精品一及| 国产精品不卡视频一区二区| 最后的刺客免费高清国语| 国产精品久久久久久亚洲av鲁大| 日本三级黄在线观看| 精品久久国产蜜桃| 成年版毛片免费区| 亚洲人与动物交配视频| 国内揄拍国产精品人妻在线| 看片在线看免费视频| 卡戴珊不雅视频在线播放| 国产女主播在线喷水免费视频网站 | 国内久久婷婷六月综合欲色啪| 国产成年人精品一区二区| 久久久久久久久久成人| 免费观看的影片在线观看| 亚洲美女黄片视频| 国产成人a区在线观看| 午夜亚洲福利在线播放| 人妻夜夜爽99麻豆av| 老司机福利观看| 国产成人福利小说| 三级经典国产精品| 在线免费观看不下载黄p国产| 国产精品av视频在线免费观看| 性插视频无遮挡在线免费观看| 22中文网久久字幕| 91精品国产九色| 欧美精品国产亚洲| 亚洲国产精品合色在线| 久久久国产成人免费| 男人舔女人下体高潮全视频| 亚洲丝袜综合中文字幕| 欧美成人a在线观看| 男人舔奶头视频| 国内久久婷婷六月综合欲色啪| or卡值多少钱| 亚洲国产高清在线一区二区三| 免费看av在线观看网站| 亚洲国产欧洲综合997久久,| 黄色一级大片看看| 一进一出抽搐动态| 久久久久久久久大av| 亚洲性久久影院| 一本一本综合久久| 少妇丰满av| 天堂av国产一区二区熟女人妻| 六月丁香七月| 最好的美女福利视频网| 亚洲在线观看片| 99国产精品一区二区蜜桃av| 又粗又爽又猛毛片免费看| 国产在视频线在精品| 欧美一区二区亚洲| av在线蜜桃| 亚洲美女搞黄在线观看 | 网址你懂的国产日韩在线| 精品无人区乱码1区二区| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av涩爱 | 成人精品一区二区免费| videossex国产| 欧美日韩在线观看h| 成人一区二区视频在线观看| 亚洲高清免费不卡视频| 日产精品乱码卡一卡2卡三| 亚洲av成人精品一区久久| 欧美zozozo另类| 亚洲精品亚洲一区二区| 丝袜喷水一区| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品亚洲一区二区| av在线老鸭窝| 99国产精品一区二区蜜桃av| 欧美又色又爽又黄视频| 成人美女网站在线观看视频| 久久韩国三级中文字幕| 亚洲精品一区av在线观看| 淫妇啪啪啪对白视频| 丝袜美腿在线中文| 午夜福利在线在线| 午夜精品一区二区三区免费看| 亚洲国产精品成人综合色| 婷婷精品国产亚洲av在线| 久久人妻av系列| 亚洲aⅴ乱码一区二区在线播放| 欧美又色又爽又黄视频| 精品久久久噜噜| 99久久无色码亚洲精品果冻| 亚洲精品一区av在线观看| 精品国内亚洲2022精品成人| 欧美xxxx黑人xx丫x性爽| 国产 一区精品| 看片在线看免费视频| 日本一二三区视频观看| 在线观看一区二区三区| av在线亚洲专区| 在线免费十八禁| 99久久中文字幕三级久久日本|