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

    MJO與北半球高緯地區(qū)冬季地表氣溫異常的聯(lián)系*

    2021-02-02 07:26:24楊雙艷LITim
    氣象 2021年1期
    關(guān)鍵詞:比濕位勢北半球

    崔 靜 楊雙艷 LI Tim,2

    1 南京信息工程大學氣象災(zāi)害教育部重點實驗室/氣候與環(huán)境變化國際合作聯(lián)合實驗室/氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心,南京 210044 2 夏威夷大學國際太平洋研究中心,夏威夷 美國

    提 要: 基于1979—2016年NCEP-NCAR逐日再分析資料研究了熱帶季節(jié)內(nèi)振蕩(MJO)和北半球冬季高緯地區(qū)地表氣溫(SAT)之間的聯(lián)系。利用實時多變量MJO(RMM)指數(shù),將MJO分為8個位相,其中位相2(位相6)對應(yīng)于位于印度洋地區(qū)的正(負)對流。不同MJO位相下的SAT合成結(jié)果顯示MJO第二位相后的5~15 d,北半球高緯地區(qū)(60°~90°N、180°~60°W)有明顯的負SAT異常;由于熱帶異常加熱信號的改變,在MJO第六位相后的5~15 d該地區(qū)則對應(yīng)于顯著的正SAT異常。該地區(qū)溫度的垂直結(jié)構(gòu)在各個位相下也表現(xiàn)出類似的分布特征。合成的500 hPa位勢高度異常場顯示,在溫度負(正)異常的位相對應(yīng)有明顯的位勢高度負(正)異常,這種環(huán)流異常主要是由與熱帶對流異常相聯(lián)系的向東北方向傳播的羅斯貝波列所引起的。通過對波活動通量的計算,推斷該東北方向傳播的羅斯貝波列很可能是羅斯貝波能量頻散的結(jié)果。合成的700 hPa比濕異常場和SAT之間在存在著較好的對應(yīng)關(guān)系,考慮到對流層中層的比濕與向下長波輻射之間存在著正相關(guān)關(guān)系,說明該溫度異常也可能與輻射過程相關(guān)。上述分析表明與MJO對流相關(guān)的大尺度環(huán)流異常對高緯地區(qū)季節(jié)內(nèi)SAT變率有重要影響,該異常SAT信號可能來自平流輸送和輻射過程等。準確把握MJO位相與SAT異常信號的聯(lián)系也可以為北半球高緯地區(qū)SAT的延伸期預(yù)報提供一些可靠線索。

    引 言

    熱帶大氣季節(jié)內(nèi)振蕩(Madden-Julian oscillation,MJO)是熱帶地區(qū)大氣環(huán)流異常的一個最重要的模態(tài),主要的時間尺度為30~60 d,主要特征表現(xiàn)為沿全球熱帶地區(qū)東傳的對流活動,具有緯向1波為主的行星尺度空間結(jié)構(gòu)(Madden and Julian,1971;1972)。MJO活動表現(xiàn)出明顯的季節(jié)演變特征。在北半球冬季,MJO以沿赤道向東的傳播為主;在北半球夏季,向東傳播會減弱,在印度夏季風區(qū)以及東亞季風區(qū)以北傳的對流活動為主導,在西北太平洋區(qū)以西北向的傳播為主導。Li(2014)指出北半球夏季熱赤道的北移是造成北印度洋和中國南海地區(qū)MJO對流北向傳播的重要原因之一,同時背景場中的東風切變和水汽的經(jīng)向分布也起到了重要貢獻(Jiang et al,2004)。關(guān)于MJO動力機制(包括其觸發(fā)機制和東傳機制)的研究工作已取得了顯著的進展(Wang and Rui,1990;Wang and Li,1994;Hsu and Li,2012;Sobel and Maloney,2013;Zhao et al,2013)。

    MJO作為熱帶最主要的環(huán)流演變模態(tài)可以直接調(diào)節(jié)熱帶地區(qū)的天氣氣候(Madden and Julian,1994)。比如,它會對熱帶氣旋和颶風的生成和加強(Maloney and Hartmann,2000;Hall et al,2001;祝從文等,2004;李崇銀等,2012)產(chǎn)生影響。MJO還會影響印度洋偶極子事件(Rao and Yamagata,2004)、北美環(huán)流異常(Higgins and Mo,1997)、南海夏季風(林愛蘭等,2016;李春暉等,2016)、孟加拉灣西南夏季風季節(jié)內(nèi)振蕩(李汀和琚建華,2013)以及ENSO事件的發(fā)展(Zhang and Gottschalck,2002)。

    前人的許多研究表明,MJO對環(huán)流異常的影響不僅僅局限于熱帶,由于對流加熱異常可以激發(fā)羅斯貝波列(Ferranti et al,1990;Kiladis and Weickmann,1992),這種影響可以延伸至中高緯地區(qū)(Lau and Phillips,1986)。比如,Bond and Vecchi(2003)指出了MJO對于俄勒岡州和華盛頓早冬和晚冬降水的重要影響;Lin and Brunet(2009)表明MJO位相與加拿大冬季地表氣溫(SAT)之間存在顯著的相關(guān)關(guān)系。除熱帶季節(jié)內(nèi)振蕩以外,中高緯的季節(jié)內(nèi)振蕩(ISO)同樣會對熱帶和中高緯的天氣氣候產(chǎn)生重要影響(Wen et al,2011;許樂心等,2017;Yang and Li,2016;2017;Yang et al,2019;李健穎和毛江玉,2019)。近年來,針對MJO與我國氣候異常間的聯(lián)系也開展了一系列研究,研究發(fā)現(xiàn)MJO對我國降水有重要的調(diào)制作用(吳捷等,2018),對中國東部地區(qū)夏季(白旭旭等,2011)和冬季(劉冬晴和楊修群,2010)、華南前汛期(李文鎧等,2014;章麗娜等,2011)、夏季西南(李永華等,2016)、云南地區(qū)(牛法寶等,2013;李汀等,2012)以及夏季廣東(林愛蘭等,2013)的降水異常均存在一定程度的影響。同時,MJO對我國氣溫異常也存在顯著影響,比如2008年初南方持續(xù)性冰凍雨雪天氣(吳俊杰等,2009;馬寧等,2011)以及2009年11月我國東部大范圍雨雪天氣(賈小龍和梁瀟云,2011)等極端天氣都與MJO存在著密切的聯(lián)系。因此,MJO的發(fā)展和演變,對于把握這種氣候異常有很大的幫助。MJO作為次季節(jié)時間尺度上最具可預(yù)報性的模態(tài)(Gottschalck et al,2010),是延伸期(10~30 d)預(yù)報主要的可預(yù)報來源(丁一匯和梁萍,2010;賈小龍等,2012)。MJO和大氣低頻季節(jié)內(nèi)振蕩的準確預(yù)報對于填補短期天氣預(yù)報和長期氣候預(yù)測之間的空隙至關(guān)重要。

    在Vecchi and Bond(2004)的工作中分析了MJO對高緯地區(qū)冬季SAT的影響,他們的合成分析均為基于MJO各個位相同時次的合成。前人的研究表明中高緯地區(qū)對于熱帶強迫的響應(yīng)在時間上存在一定的滯后,完全發(fā)展建立需要兩周左右的時間(Jin and Hoskins,1995)。為使這種異常SAT信號更加顯著,我們的工作主要關(guān)注了MJO與北半球冬季高緯地區(qū)SAT之間的滯后相關(guān)關(guān)系,從環(huán)流分析等角度進一步揭示兩者之間的聯(lián)系。旨在將MJO和高緯地區(qū)的大氣變率相聯(lián)系,從而為該地區(qū)季節(jié)內(nèi)SAT變率的延伸期預(yù)報提供參考。

    1 資料和方法

    1.1 資 料

    本文使用了美國國家海洋與大氣中心(NOAA)的逐日向外長波輻射(OLR)資料以表征熱帶大尺度對流特征;同時使用了美國國家環(huán)境預(yù)測中心-國家大氣研究中心(NCEP-NCAR)發(fā)布的全球逐日再分析資料(水平分辨率為2.5°×2.5°)。本文的分析時段為1979/1980年到2015/2016年共37個冬季,冬季定義為11月到次年3月共5個月。再分析資料中使用到的氣象要素包括2 m處的地表氣溫(T2 m)以及在各個高度層次的緯向風(u)、經(jīng)向風(v)、溫度(T)、p坐標系下的垂直速度(ω)等。

    采用由Wheeler and Hendon(2004)所定義的實時多變量MJO指數(shù)(real-time multivariate MJO indices,簡稱RMM指數(shù))來描述MJO的活動特征,RMM指數(shù)最常用于MJO的監(jiān)測和預(yù)報。本文使用到的逐日RMM指數(shù)可以在澳大利亞氣象局官網(wǎng)下載(http:∥www.bom.gov.au/climate/mjo/)。該指數(shù)是將逐日資料投影到熱帶地區(qū)(15°S~15°N)緯向平均后的逐日的850 hPa緯向風(u850)、200 hPa 緯向風(u200)和OLR的前兩個多變量EOF模態(tài)上,得到的前兩個主成分時間序列,將該時間序列標準化即可得到RMM指數(shù),分別記為RMM1和RMM2。上述資料均去除了年循環(huán)和年際變率。根據(jù)RMM1和RMM2所確定的二維位相空間,可以將強的MJO 沿全球熱帶地區(qū)的緯向傳播分為8個空間位相。MJO的地理位置可以由一個以RMM1和RMM2指數(shù)所定義的二維位相空間中的一個點來表示,而該點與整個空間中心的距離則表示MJO的振幅,代表MJO事件的強度,這8個位相分別對應(yīng)于MJO生命史的特定階段,表示MJO對流從印度洋到太平洋的向東傳播。

    1.2 方 法

    為提取與MJO相關(guān)的季節(jié)內(nèi)信號,對原始數(shù)據(jù)進行了相關(guān)處理。首先去除了1979—2016年多年氣候態(tài)的0~3波以去除季節(jié)循環(huán),然后通過使用Lanczos帶通濾波方法(Duchon,1979)進行30~60 d 的帶通濾波,濾波后選擇了11月至次年3月的資料用于分析,但1979年1月至次年3月以及2016年11—12月的資料未包括在內(nèi)。

    為檢驗每個MJO位相合成結(jié)果是否顯著,使用到的t檢驗公式為:

    式中:ρ1為滯后一天的相關(guān)系數(shù),m2是yi的樣本數(shù)。

    2 MJO位相合成的高緯地表氣溫異常

    2.1 MJO位相

    為了證實赤道對流活動是和MJO的位相相聯(lián)系的,圖1給出了各個位相下合成的OLR場及850 hPa的經(jīng)向風和緯向風異常場。在位相1,MJO對流中心位于非洲和西印度洋,此時太平洋以西風異常為主,而印度洋存在東風異常;在位相2和位相3,對流中心主要位于印度洋,繼續(xù)發(fā)展并逐漸東移至海洋性大陸(位相4、位相5);位相6~位相8,對流中心由西太平洋繼續(xù)東移至中太平洋。

    使用到的RMM指數(shù)時間段包括從1979/1980—2015/2016年的37個冬季,因此可以得到每個位相所對應(yīng)的天數(shù)。圖2給出了每個MJO位相所對應(yīng)的天數(shù)及平均振幅的大小??梢钥吹轿幌?~位相4以及位相6~位相8 MJO事件的發(fā)生頻率較高,其中平均振幅的最大值出現(xiàn)位相3。

    圖1 MJO 8個位相(a~h)下合成的30~60 d濾波后的冬季OLR異常場(填色,單位:W·m-2)和850 hPa異常風場(箭頭,單位:m·s-1)(陰影區(qū)為通過0.05顯著性水平檢驗區(qū)域,t檢驗;右下角數(shù)字為每個位相對應(yīng)的天數(shù))

    圖2 二維MJO指數(shù)(RMM1和RMM2)定義的MJO空間位相以及各位相對應(yīng)的天數(shù)和平均振幅

    2.2 與MJO相聯(lián)系的高緯地區(qū)地表氣溫異常

    為探究北半球冬季SAT和MJO之間的聯(lián)系,對不同MJO位相下的SAT進行了合成分析。在北半球高緯地區(qū)這種SAT異常信號較為顯著。Jinand Hoskins(1995)研究表明:給定的熱帶加熱響應(yīng)下的羅斯貝波列會在15 d以內(nèi)完全建立。因此在下文的合成分析中將分別采取0,5,10,15 d的滯后。

    MJO的位相1~位相4代表其一半的生命周期,而位相5~位相8則代表相反的另一半周期,這種生命史的對稱特征也表現(xiàn)在SAT的合成場上,隨著滯后時間的變化,每個位相下高緯SAT異常的強度也會有一定的變化,總體來說在滯后5~15 d時位相2~位相7下高緯地區(qū)SAT信號比較顯著,第0天時(即同時次的合成圖中)這種異常SAT信號相對較弱,隨著滯后時間的增加,該異常信號得到加強。這里給出了滯后5~15 d的MJO的8個位相下的異常SAT合成圖(圖3),可以看到位相2和位相6以及位相3和位相7下這種SAT信號比較顯著,并且該異常SAT信號存在著明顯的反向?qū)?yīng)關(guān)系(即位相2、6表現(xiàn)為明顯的負異常而位相3、位相7則表現(xiàn)為明顯的正異常),后續(xù)為了簡潔起見選擇位相2和位相6進行研究(對位相3和位相7也做了同樣分析,得到相似結(jié)論)。從圖4中可以看到在MJO的第2(第6)位相后5 d左右北半球高緯地區(qū)開始出現(xiàn)負(正)的SAT異常,在滯后10~15 d時這種異常信號表現(xiàn)得更加明顯。

    此外,值得關(guān)注的是在高緯地區(qū)(60°N以北),其SAT的異常情況并不是完全一致的,部分區(qū)域會呈現(xiàn)出相反的變化,在滯后5 d以后其SAT異常信號均較為顯著和一致的主要為60°~90°N、180°~60°W這部分區(qū)域。為了探究該區(qū)域溫度的垂直結(jié)構(gòu)是否也有類似的變化特征,圖5給出了MJO位相2和位相6滯后0,5,10,15 d后的溫度異常及位勢高度異常沿緯向平均(180°~60°W)后的經(jīng)向-垂直分布。和地表氣溫類似,整個對流層的溫度異常也呈現(xiàn)出相似的分布特征,即在滯后5~10 d時,位相2(位相6)呈現(xiàn)出明顯的負(正)溫度異常。另一個顯著的特征是:例如在位相2后的5 d,一個位于300 hPa左右的負位勢高度異常中心出現(xiàn)在70°~80°N,在該中心以上(以下)溫度異常表現(xiàn)為暖(冷)異常,較好地滿足了靜力平衡關(guān)系。在第6位相也能看到這種較好的對應(yīng)關(guān)系,靜力平衡關(guān)系使得低層的與MJO相關(guān)的溫度異常信號與對流層高層的位勢高度異常之間建立了密切的關(guān)系。最后,由于在北半球高緯地區(qū)冬季有相似幅度的增暖(冷卻)發(fā)生在MJO第1,3,4(5,7,8)位相后的1~2周,所以該地區(qū)冬季SAT與MJO位相2和位相6之間的聯(lián)系同樣應(yīng)該可適用于MJO的其他位相。

    圖3 MJO 8個位相(a~h)下冬季SAT滯后5~15 d的合成(打點代表通過0.05的顯著性水平檢驗)

    圖4 MJO位相2(a)和位相6(b)下冬季SAT滯后0 d(a1,b1),5 d(a2,b2),10 d(a3,b3),15 d(a4,b4)的合成(打點代表通過0.05的顯著性水平檢驗)

    圖5 MJO位相2(a)和位相6(b)下滯后0 d(a1,b1),5 d(a2,b2),10 d(a3,b3),15 d(a4,b4)的冬季溫度異常場(填色)以及位勢高度異常場(等值線,單位:gpm)沿180°~60°W平均后的經(jīng)向-垂直分布(打點代表通過0.05顯著性水平檢驗)

    3 大尺度環(huán)流異常

    為了進一步理解上述北半球高緯地區(qū)與MJO相聯(lián)系的SAT信號,探究了與MJO位相相聯(lián)系的大尺度異常環(huán)流場的分布特征。Ferranti et al(1990),Matthews et al(2004)研究表明:與MJO相關(guān)的熱帶非絕熱加熱異常會激發(fā)羅斯貝波列并影響中高緯的大尺度環(huán)流異常。

    圖6給出了在冬季MJO位相2和位相6下合成的北半球500 hPa位勢高度異常場。對于MJO位相2,在滯后5 d時,一個顯著的特征是有一個位于北太平洋,中心位于40°N、180°附近的正位勢高度異常中心。在阿拉斯加和加拿大西岸有一個負位勢高度異常中心,加拿大東部則為正位勢高度異常,可以看到波的模態(tài)主要表現(xiàn)出向下游頻散的特征,與PNA的模態(tài)很類似,但是相對于PNA有略微的西移;到滯后10 d時,北大西洋分支有明顯的環(huán)流異常在發(fā)展,形成了近似于正的NAO的模態(tài)。位相6的合成圖表現(xiàn)出和位相2相似的演變特征,但由于此時的熱帶強迫發(fā)生了變化,因此呈現(xiàn)出的波列模態(tài)也相反。

    由圖6可以看到北半球高緯表現(xiàn)出負(正)SAT異常的地區(qū)主要對應(yīng)于負(正)位勢高度異常,這表明MJO對北半球高緯地區(qū)SAT信號的影響可能與水平溫度平流的輸送相關(guān)。這種位勢高度的異常在對流層整層基本是一致的,呈現(xiàn)出正壓的結(jié)構(gòu)。參考前人的工作,可以推測與MJO相關(guān)的熱帶對流所引起的加熱異??赡軙|發(fā)羅斯貝波列,由此持續(xù)性地影響了北半球高緯地區(qū)的環(huán)流異常,從而對季節(jié)內(nèi)SAT信號產(chǎn)生重要影響。

    圖6 同圖5,但為合成的冬季500 hPa位勢高度異常場(等值線,間隔為2 gpm)(陰影為通過0.05的顯著性水平檢驗)

    為探究熱力因素與異常SAT信號之間是否也存在著一定的對應(yīng)關(guān)系,對MJO位相下合成的700 hPa比濕場異常場(圖8)進行了簡單分析。可以看到MJO位相2(位相6)下較冷(暖)的SAT分別對應(yīng)于700 hPa較干(濕)的空氣,考慮到對流層中層比濕和向下長波輻射之間又存在正相關(guān)關(guān)系(Overland et al,1997; Adams et al,2000),可以推測不同MJO位相下的比濕正(負)異常會導致向下長波輻射的正(負)異常,從而會對該地區(qū)的正(負)SAT信號有一定的正貢獻。

    由上述分析發(fā)現(xiàn)與MJO對流相聯(lián)系熱帶加熱異常會觸發(fā)東北向傳播的羅斯貝波列,該波列會影響到研究所關(guān)注的高緯地區(qū)的大尺度環(huán)流異常,從而對該地區(qū)的季節(jié)內(nèi)SAT異常信號產(chǎn)生重要影響。此外,除異常環(huán)流場以外,700 hPa比濕場和高緯SAT信號間也有較好的對應(yīng)關(guān)系,表明與輻射過程相關(guān)的熱力因素可能對該地區(qū)的SAT異常有一定的正貢獻。而具體哪一過程在溫度異常的變化中起主導作用還需借助于溫度方程進行診斷分析,在之后的工作中我們會作進一步的研究。

    圖7 MJO位相2(a)和位相6(b)合成的200 hPa波活動通量(箭頭,單位: m2·s-2)及羅斯貝波源(填色,單位:10-10 s-2)

    圖8 同圖5,但為合成的冬季SAT(填色,僅顯示通過0.05顯著性水平檢驗的部分)和700 hPa比濕異常場(等值線,間隔為2×10-5 kg·kg-1)

    4 結(jié)論和討論

    本文研究了MJO和北半球冬季高緯地區(qū)SAT之間的聯(lián)系,基于38年(1979—2016年)的逐日2 m地表氣溫資料以及RMM指數(shù),得到了MJO熱帶對流和北半球冬季高緯地區(qū)SAT異常之間顯著的相關(guān)關(guān)系。全文主要關(guān)注了MJO生命周期中的兩個位相:位相2和位相6,分別對應(yīng)于位于東印度洋和西太平洋地區(qū)的MJO正對流??紤]到中高緯對于熱帶對流的響應(yīng)存在1~2周的延遲,為使與MJO相聯(lián)系的異常信號更顯著,在進行MJO位相的合成分析時考慮了5~15 d的滯后。可以明顯看到:在北半球高緯地區(qū),尤其是60°~90°N、180°~60°W這部分區(qū)域,在MJO位于位相2后的5~15 d有明顯的負SAT異常,此時的正對流主要位于印度洋和海洋性大陸地區(qū),負對流位于太平洋西部及中部;同樣的,在MJO位于位相6后的5~15 d,該區(qū)域有顯著的正SAT異常,對應(yīng)位于印度洋和海洋性大陸的負對流以及位于太平洋西部及中部的正對流。該地區(qū)的溫度異常的垂直分布特征與SAT類似,即在MJO位于位相2后的5~15 d表現(xiàn)為負溫度異常,而在MJO位于位相6后的5~15 d則表現(xiàn)為正溫度異常。滯后合成的大尺度異常環(huán)流場表明熱帶對流活動會對中高緯的SAT產(chǎn)生影響,通常是以羅斯貝波列的形式。500 hPa位勢高度異常的合成場上可以明顯看到在太平洋以及北美地區(qū)形成了類似于PNA模態(tài)的波列。羅斯貝波源及波活動通量的分析表明該波列應(yīng)該是由與MJO對流相聯(lián)系的熱帶加熱異常所觸發(fā)的,同時在加熱場東北側(cè)的正波源區(qū)有東北指向的波活動通量,表明觀測的羅斯貝波列可能是東北方向頻散的羅斯貝能量的結(jié)果。熱帶對流通過這種遙相關(guān)的方式對中高緯的環(huán)流形勢產(chǎn)生影響,從而影響季節(jié)內(nèi)SAT變率。此外700 hPa的比濕異常場的合成分析圖也呈現(xiàn)出與SAT異常信號間較好的對應(yīng)關(guān)系,在高緯地區(qū)溫度負異常時有明顯的比濕負異常,相反的,在溫度正異常時則有明顯的比濕正異常,而對流層中層的比濕和向下長波輻射之間存在正相關(guān)關(guān)系,表明輻射過程等熱力因素也可能會對該地區(qū)的SAT異常有一定的貢獻。

    本文中MJO位相和SAT異常間顯著的相關(guān)關(guān)系可能會給北半球高緯地區(qū)SAT的預(yù)報提供一些幫助,因為熱帶地區(qū)緩慢變化的對流活動可以為預(yù)報提供季節(jié)內(nèi)時間尺度上的信息。但是也認識到了本次工作中的幾個局限,首先,上文所討論的滯后合成分析僅僅是多個MJO事件的平均結(jié)果,而這些MJO事件通常具有不同的傳播特征和強度。此外,北極的SAT異常更多是受到高緯地區(qū)自身的變化的影響,比如西伯利亞高壓和北極濤動(AO)等,所以MJO的影響只能被視作是一個附加部分并且它的影響也是有限的。最后,還應(yīng)該注意到,研究不可避免地會受到和厄爾尼諾等相關(guān)的更長時間尺度變率的影響。熱帶MJO是進行10~30 d延伸期預(yù)報的重要可預(yù)報來源,如果能準確預(yù)測MJO的發(fā)展和演變,對于氣候異常的預(yù)測可以提供一些依據(jù)。

    猜你喜歡
    比濕位勢北半球
    北半球最強“星空攝影師”開工啦
    軍事文摘(2023年24期)2023-12-19 06:50:06
    雅安地區(qū)近50 年濕度變化特征分析
    含Hardy位勢的非線性Schr?dinger-Poisson方程正規(guī)化解的多重性
    一類帶強制位勢的p-Laplace特征值問題
    清涼一夏
    基于探空資料的1961—2018年新疆高空大氣比濕氣候特征分析
    南方比濕特征及其與暴雨的關(guān)系
    南北半球天象
    軍事文摘(2019年18期)2019-09-25 08:09:22
    2017年朝陽市地面比濕特征分析
    含變號位勢的ρ-Kirchhoff型方程組無窮多個高能量解的存在性
    91精品伊人久久大香线蕉| 欧美xxxx性猛交bbbb| 亚洲国产欧美日韩在线播放| 建设人人有责人人尽责人人享有的| 自线自在国产av| 老女人水多毛片| 久久精品国产鲁丝片午夜精品| 婷婷色综合大香蕉| 亚洲av日韩在线播放| www.av在线官网国产| 亚洲精品乱码久久久久久按摩| 国产精品无大码| 国产国语露脸激情在线看| 超色免费av| 18在线观看网站| 亚洲国产av新网站| 国产一区二区在线观看av| 欧美bdsm另类| 一区二区三区精品91| 国产黄频视频在线观看| 最黄视频免费看| 男女边摸边吃奶| 精品久久久久久久久亚洲| 久久影院123| 狂野欧美激情性bbbbbb| a级片在线免费高清观看视频| 天天躁夜夜躁狠狠久久av| 视频在线观看一区二区三区| 久久99精品国语久久久| 婷婷色麻豆天堂久久| 国产精品国产三级专区第一集| 国产男女超爽视频在线观看| 欧美精品一区二区大全| 亚洲激情五月婷婷啪啪| 日本色播在线视频| 九九久久精品国产亚洲av麻豆| 美女xxoo啪啪120秒动态图| 水蜜桃什么品种好| 熟女av电影| 婷婷色麻豆天堂久久| 一级毛片电影观看| 欧美丝袜亚洲另类| 国产探花极品一区二区| 免费播放大片免费观看视频在线观看| 美女福利国产在线| 国产日韩欧美亚洲二区| 男女免费视频国产| 亚洲内射少妇av| 男人爽女人下面视频在线观看| 日本欧美视频一区| 一级,二级,三级黄色视频| xxx大片免费视频| 能在线免费看毛片的网站| 乱人伦中国视频| 久久精品人人爽人人爽视色| 亚洲欧洲国产日韩| 男男h啪啪无遮挡| 熟妇人妻不卡中文字幕| 久久国产精品大桥未久av| 免费播放大片免费观看视频在线观看| videossex国产| 亚洲图色成人| 男女免费视频国产| 久久久久人妻精品一区果冻| 晚上一个人看的免费电影| 日本黄色日本黄色录像| 五月伊人婷婷丁香| 中文字幕av电影在线播放| 中文字幕亚洲精品专区| 亚洲人成77777在线视频| 精品一区二区免费观看| av不卡在线播放| 久久国产精品男人的天堂亚洲 | 黑人巨大精品欧美一区二区蜜桃 | 如日韩欧美国产精品一区二区三区 | 亚洲美女黄色视频免费看| 亚洲国产精品一区三区| 另类精品久久| 免费大片黄手机在线观看| 国产欧美日韩综合在线一区二区| 91精品国产国语对白视频| 精品99又大又爽又粗少妇毛片| 成人免费观看视频高清| 18禁裸乳无遮挡动漫免费视频| 欧美 亚洲 国产 日韩一| 免费少妇av软件| 男女免费视频国产| 最黄视频免费看| 亚洲美女黄色视频免费看| 亚洲人与动物交配视频| 精品少妇久久久久久888优播| 久久精品久久久久久噜噜老黄| 国产av精品麻豆| 久久av网站| 乱人伦中国视频| 久久久久久久久久久丰满| 在线观看免费视频网站a站| 国产av一区二区精品久久| 99精国产麻豆久久婷婷| 欧美精品一区二区大全| 亚洲精品一区蜜桃| 欧美亚洲日本最大视频资源| 亚洲成人一二三区av| 精品久久久噜噜| 国产高清不卡午夜福利| 亚洲精品色激情综合| 91精品伊人久久大香线蕉| 中文字幕制服av| 男女无遮挡免费网站观看| 在线观看国产h片| 各种免费的搞黄视频| 欧美性感艳星| videossex国产| 在线亚洲精品国产二区图片欧美 | 日本免费在线观看一区| 日本wwww免费看| av卡一久久| 99热这里只有精品一区| 麻豆成人av视频| freevideosex欧美| 一本色道久久久久久精品综合| 插阴视频在线观看视频| 国产日韩欧美亚洲二区| 女性生殖器流出的白浆| 成人国产麻豆网| 视频中文字幕在线观看| 中文字幕精品免费在线观看视频 | 亚洲人成网站在线播| 国产成人一区二区在线| 久久精品国产a三级三级三级| 日日爽夜夜爽网站| 午夜精品国产一区二区电影| 精品视频人人做人人爽| 国产免费现黄频在线看| 视频在线观看一区二区三区| 男人操女人黄网站| 免费日韩欧美在线观看| 人妻 亚洲 视频| 国产黄频视频在线观看| 我要看黄色一级片免费的| 午夜免费观看性视频| 欧美精品人与动牲交sv欧美| 美女福利国产在线| 午夜视频国产福利| 成人手机av| av国产精品久久久久影院| 中文字幕人妻熟人妻熟丝袜美| 永久免费av网站大全| 18禁在线播放成人免费| 一个人免费看片子| 亚洲精品视频女| 精品少妇久久久久久888优播| 国产乱人偷精品视频| 亚洲av日韩在线播放| 在线观看国产h片| 国产爽快片一区二区三区| 中文字幕人妻丝袜制服| 男女国产视频网站| 久久久久久久大尺度免费视频| 国产av精品麻豆| 夫妻午夜视频| 亚洲av在线观看美女高潮| 久久99一区二区三区| 成人影院久久| 亚洲精品乱码久久久v下载方式| 高清黄色对白视频在线免费看| 久久影院123| 自线自在国产av| 纵有疾风起免费观看全集完整版| 日日摸夜夜添夜夜爱| 日本欧美视频一区| 精品国产乱码久久久久久小说| 婷婷色综合www| 18+在线观看网站| 少妇丰满av| 婷婷色综合www| 国产又色又爽无遮挡免| 国产成人91sexporn| 丝袜脚勾引网站| 伊人久久精品亚洲午夜| 少妇的逼好多水| 蜜桃国产av成人99| 夜夜看夜夜爽夜夜摸| 日韩av免费高清视频| 99久久精品国产国产毛片| 国产精品人妻久久久久久| 2022亚洲国产成人精品| videossex国产| 两个人免费观看高清视频| 看非洲黑人一级黄片| 热99国产精品久久久久久7| 在线播放无遮挡| 国产伦精品一区二区三区视频9| 王馨瑶露胸无遮挡在线观看| 男女无遮挡免费网站观看| 一区在线观看完整版| 一本大道久久a久久精品| 在线播放无遮挡| 久久久久久伊人网av| 日韩 亚洲 欧美在线| 在线播放无遮挡| 七月丁香在线播放| 欧美精品亚洲一区二区| 80岁老熟妇乱子伦牲交| 欧美人与性动交α欧美精品济南到 | 夫妻性生交免费视频一级片| 亚洲欧洲日产国产| 精品99又大又爽又粗少妇毛片| 人成视频在线观看免费观看| 国产黄片视频在线免费观看| 亚洲欧美色中文字幕在线| 久久99蜜桃精品久久| 人妻人人澡人人爽人人| 97超碰精品成人国产| 国产成人免费观看mmmm| 韩国高清视频一区二区三区| 在线天堂最新版资源| av天堂久久9| 久久久国产欧美日韩av| 欧美日韩成人在线一区二区| 亚洲国产最新在线播放| 国产又色又爽无遮挡免| 日本vs欧美在线观看视频| 久久人人爽人人片av| 人体艺术视频欧美日本| 国产亚洲最大av| 十八禁高潮呻吟视频| 成人毛片a级毛片在线播放| 天堂8中文在线网| 亚洲第一av免费看| 亚洲天堂av无毛| 中文字幕精品免费在线观看视频 | 少妇被粗大猛烈的视频| 曰老女人黄片| 亚洲久久久国产精品| 国产成人精品婷婷| 国产精品蜜桃在线观看| 纵有疾风起免费观看全集完整版| 91久久精品国产一区二区成人| 人妻 亚洲 视频| www.色视频.com| 青春草国产在线视频| 欧美xxⅹ黑人| 一区二区三区四区激情视频| 美女中出高潮动态图| 久久精品国产自在天天线| 欧美精品亚洲一区二区| 久久久久久久久久久免费av| av天堂久久9| 国产成人精品一,二区| 亚洲色图 男人天堂 中文字幕 | 久久国内精品自在自线图片| 中国美白少妇内射xxxbb| 亚洲av.av天堂| 另类精品久久| 亚洲av欧美aⅴ国产| 国产成人a∨麻豆精品| 午夜福利视频在线观看免费| 亚洲精品亚洲一区二区| 一级毛片我不卡| 伊人久久精品亚洲午夜| 一级毛片电影观看| 亚洲美女搞黄在线观看| 亚洲精品,欧美精品| 最后的刺客免费高清国语| 中文字幕久久专区| 亚洲欧洲国产日韩| 天堂俺去俺来也www色官网| 人妻夜夜爽99麻豆av| 国精品久久久久久国模美| 亚洲精华国产精华液的使用体验| 久久久久视频综合| 国产欧美日韩一区二区三区在线 | 久久久久久人妻| 最新的欧美精品一区二区| 大香蕉97超碰在线| 色94色欧美一区二区| 日韩av免费高清视频| 飞空精品影院首页| 狠狠婷婷综合久久久久久88av| 黑人猛操日本美女一级片| 一级黄片播放器| 免费av不卡在线播放| 各种免费的搞黄视频| 婷婷色av中文字幕| 国产成人精品久久久久久| 在线观看国产h片| 精品亚洲乱码少妇综合久久| freevideosex欧美| 18禁在线无遮挡免费观看视频| 99久久精品一区二区三区| 午夜av观看不卡| 日韩免费高清中文字幕av| kizo精华| 国产一区亚洲一区在线观看| 国产乱来视频区| 中文字幕人妻丝袜制服| 人人妻人人添人人爽欧美一区卜| 香蕉精品网在线| 一级毛片黄色毛片免费观看视频| 亚洲情色 制服丝袜| 99热这里只有精品一区| av.在线天堂| 黑人高潮一二区| 曰老女人黄片| 免费看不卡的av| 看非洲黑人一级黄片| 亚洲成人手机| 精品亚洲成国产av| 黑丝袜美女国产一区| 亚洲欧洲日产国产| 欧美 日韩 精品 国产| 色94色欧美一区二区| 97在线视频观看| 九九爱精品视频在线观看| 麻豆乱淫一区二区| 日本黄色日本黄色录像| 国产成人免费观看mmmm| 国产日韩欧美视频二区| 久久久久久久久久人人人人人人| 国产亚洲最大av| 成年人午夜在线观看视频| 国产欧美亚洲国产| 精品久久蜜臀av无| 91精品国产九色| 三级国产精品欧美在线观看| 91久久精品国产一区二区成人| 精品一区二区三卡| 亚洲精品日韩av片在线观看| 国产爽快片一区二区三区| 免费高清在线观看日韩| 99国产综合亚洲精品| 国产一区二区三区av在线| 熟女电影av网| 亚洲av男天堂| 黑人欧美特级aaaaaa片| 少妇 在线观看| 赤兔流量卡办理| 美女大奶头黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久精品精品| 啦啦啦啦在线视频资源| 国产在线视频一区二区| 男女啪啪激烈高潮av片| 少妇的逼好多水| 国产老妇伦熟女老妇高清| 又黄又爽又刺激的免费视频.| 日本-黄色视频高清免费观看| 一本色道久久久久久精品综合| 午夜福利视频精品| 纯流量卡能插随身wifi吗| 黄色配什么色好看| 日本av免费视频播放| 观看美女的网站| 久热久热在线精品观看| 免费看不卡的av| 亚洲激情五月婷婷啪啪| √禁漫天堂资源中文www| 老司机影院成人| 99久久精品一区二区三区| 国产精品 国内视频| av天堂久久9| 免费观看的影片在线观看| 涩涩av久久男人的天堂| 寂寞人妻少妇视频99o| av一本久久久久| 三上悠亚av全集在线观看| 你懂的网址亚洲精品在线观看| 欧美精品国产亚洲| av在线播放精品| 女的被弄到高潮叫床怎么办| 五月开心婷婷网| 女人精品久久久久毛片| av不卡在线播放| 国产午夜精品一二区理论片| 汤姆久久久久久久影院中文字幕| 色视频在线一区二区三区| 亚洲图色成人| 亚洲国产精品一区二区三区在线| 七月丁香在线播放| av不卡在线播放| 考比视频在线观看| 如日韩欧美国产精品一区二区三区 | 国产乱人偷精品视频| 欧美xxⅹ黑人| 免费久久久久久久精品成人欧美视频 | 黄色一级大片看看| 少妇丰满av| 亚洲综合精品二区| 免费少妇av软件| 边亲边吃奶的免费视频| 国国产精品蜜臀av免费| 麻豆乱淫一区二区| 亚洲成色77777| 久久久精品94久久精品| 久久精品久久久久久久性| 免费看不卡的av| 亚洲av.av天堂| 99久国产av精品国产电影| 精品久久久噜噜| 欧美精品一区二区免费开放| 一级二级三级毛片免费看| 久久精品久久久久久久性| 久久久欧美国产精品| 国产成人精品久久久久久| 九色成人免费人妻av| 中文乱码字字幕精品一区二区三区| 久久久欧美国产精品| 亚洲av成人精品一二三区| 日日爽夜夜爽网站| 国产在线一区二区三区精| 91午夜精品亚洲一区二区三区| 中文字幕av电影在线播放| 精品亚洲成a人片在线观看| 日韩av不卡免费在线播放| 观看美女的网站| 国产一区有黄有色的免费视频| 午夜日本视频在线| 亚洲av男天堂| 久久99热这里只频精品6学生| 久久精品久久久久久久性| 久久久久国产网址| 五月玫瑰六月丁香| 建设人人有责人人尽责人人享有的| 亚洲经典国产精华液单| 日韩av免费高清视频| 如何舔出高潮| 亚洲激情五月婷婷啪啪| 久久国产亚洲av麻豆专区| 久热这里只有精品99| 亚洲精品亚洲一区二区| 中文字幕av电影在线播放| 久久婷婷青草| 亚洲熟女精品中文字幕| 一本大道久久a久久精品| 欧美变态另类bdsm刘玥| 国产精品成人在线| 免费观看无遮挡的男女| 国产在线一区二区三区精| 国产女主播在线喷水免费视频网站| 精品少妇久久久久久888优播| 国产在线一区二区三区精| 高清黄色对白视频在线免费看| 日韩av免费高清视频| 免费看av在线观看网站| 在线天堂最新版资源| 男女高潮啪啪啪动态图| 自线自在国产av| 成人毛片60女人毛片免费| 久久人人爽av亚洲精品天堂| 久久精品国产a三级三级三级| 国产永久视频网站| 久久综合国产亚洲精品| 蜜桃久久精品国产亚洲av| 久久久久久久久久人人人人人人| 国产亚洲欧美精品永久| 午夜激情福利司机影院| 熟女av电影| 色吧在线观看| 亚洲欧美日韩另类电影网站| 亚洲在久久综合| 一本色道久久久久久精品综合| 九色亚洲精品在线播放| 亚洲无线观看免费| 搡女人真爽免费视频火全软件| 国产片特级美女逼逼视频| 久久久久久久久久成人| 人人妻人人澡人人看| 国产精品嫩草影院av在线观看| 日韩精品免费视频一区二区三区 | 美女xxoo啪啪120秒动态图| 精品99又大又爽又粗少妇毛片| 国产有黄有色有爽视频| 伦理电影免费视频| 丝袜在线中文字幕| 日本免费在线观看一区| 国语对白做爰xxxⅹ性视频网站| 久久精品国产鲁丝片午夜精品| 国产成人精品久久久久久| 国产高清国产精品国产三级| 亚洲av综合色区一区| 老司机影院毛片| 日韩av免费高清视频| 国产精品久久久久久久久免| 欧美日韩一区二区视频在线观看视频在线| 97在线视频观看| 国产乱来视频区| 亚洲人成77777在线视频| 免费久久久久久久精品成人欧美视频 | 丝袜在线中文字幕| 成年美女黄网站色视频大全免费 | 高清毛片免费看| 国产免费一级a男人的天堂| 国产综合精华液| 一级爰片在线观看| 免费高清在线观看视频在线观看| 国产精品久久久久久av不卡| 亚洲av成人精品一区久久| 日韩av不卡免费在线播放| 日本av手机在线免费观看| 国产高清三级在线| 国产 一区精品| 水蜜桃什么品种好| 人体艺术视频欧美日本| 人人妻人人添人人爽欧美一区卜| 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 亚洲综合精品二区| 精品人妻熟女毛片av久久网站| 日韩av免费高清视频| 伊人久久国产一区二区| 只有这里有精品99| 女性被躁到高潮视频| 欧美+日韩+精品| 日本猛色少妇xxxxx猛交久久| 99久久人妻综合| 91国产中文字幕| 欧美xxxx性猛交bbbb| 在线看a的网站| 一二三四中文在线观看免费高清| 又黄又爽又刺激的免费视频.| 久久精品久久久久久噜噜老黄| 亚洲精品第二区| 亚洲av日韩在线播放| 精品国产国语对白av| 亚洲国产日韩一区二区| 搡老乐熟女国产| 成人综合一区亚洲| 久久久久久伊人网av| 熟妇人妻不卡中文字幕| 欧美亚洲 丝袜 人妻 在线| 免费日韩欧美在线观看| 久久国产精品大桥未久av| av在线播放精品| 亚洲,欧美,日韩| av一本久久久久| 亚洲精品久久成人aⅴ小说 | 久久 成人 亚洲| 亚洲无线观看免费| 日韩av免费高清视频| 乱码一卡2卡4卡精品| 精品国产一区二区久久| 久热久热在线精品观看| 欧美精品一区二区大全| 人妻人人澡人人爽人人| 成年人免费黄色播放视频| 丰满迷人的少妇在线观看| 免费人成在线观看视频色| 99热网站在线观看| 18禁观看日本| 亚洲精品日韩av片在线观看| 青春草国产在线视频| 国产熟女欧美一区二区| 精品少妇久久久久久888优播| 午夜免费男女啪啪视频观看| 色哟哟·www| av在线app专区| 亚洲丝袜综合中文字幕| 十分钟在线观看高清视频www| 熟妇人妻不卡中文字幕| 亚洲欧美成人综合另类久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品乱久久久久久| 欧美成人午夜免费资源| 国产精品一国产av| 一边摸一边做爽爽视频免费| 久久久亚洲精品成人影院| a级毛片黄视频| 亚洲av日韩在线播放| 2021少妇久久久久久久久久久| 成人亚洲欧美一区二区av| 丝袜美足系列| 在线天堂最新版资源| 99久国产av精品国产电影| 日本黄大片高清| 在现免费观看毛片| 亚洲欧美一区二区三区黑人 | 亚洲第一区二区三区不卡| 精品久久久精品久久久| 国产熟女欧美一区二区| 欧美 亚洲 国产 日韩一| 免费观看av网站的网址| 久久久久精品性色| 国产免费一级a男人的天堂| 国产在视频线精品| 亚洲国产精品一区二区三区在线| 精品人妻熟女av久视频| 男女啪啪激烈高潮av片| 秋霞在线观看毛片| 91国产中文字幕| 欧美成人精品欧美一级黄| 欧美日韩精品成人综合77777| 国产爽快片一区二区三区| 最新中文字幕久久久久| 久久精品国产a三级三级三级| 我的女老师完整版在线观看| 精品亚洲乱码少妇综合久久| 日本wwww免费看| 国产日韩欧美亚洲二区| 国产成人免费观看mmmm| 91成人精品电影| xxx大片免费视频| 黄色欧美视频在线观看| 欧美精品国产亚洲| h视频一区二区三区| a级毛片免费高清观看在线播放| 日韩一区二区视频免费看| 日本猛色少妇xxxxx猛交久久| 日韩欧美一区视频在线观看| √禁漫天堂资源中文www| 午夜免费男女啪啪视频观看| 中文字幕最新亚洲高清| 国产精品不卡视频一区二区| 天美传媒精品一区二区| videosex国产|