亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        北太平洋濤動(dòng)和東亞冬季風(fēng)對(duì)ENSO影響差異分析

        2021-08-17 05:16:34馮立成陳權(quán)亮丁瑞強(qiáng)
        關(guān)鍵詞:大氣區(qū)域影響

        溫 濤 ,馮立成 ,陳權(quán)亮 ,丁瑞強(qiáng) ,史 良

        (1.成都信息工程大學(xué)、高原大氣與環(huán)境四川省重點(diǎn)實(shí)驗(yàn)室,四川 成都 610225; 2.國(guó)家海洋環(huán)境預(yù)報(bào)中心, 北京 100081; 3.北京師范大學(xué)、地表過(guò)程與資源生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100089; 4.蘭州理工大學(xué)理學(xué)院, 甘肅 蘭州 730050)

        早在1969年,Bjerknes就歷史性地提出厄爾尼諾-南方濤動(dòng)(EI Nio-Southern Oscillation, ENSO)是熱帶太平洋大尺度海氣相互作用的產(chǎn)物,是全球海洋-大氣耦合系統(tǒng)在年際尺度上最主要的模態(tài)[1]。前人的系列研究表明,ENSO不僅對(duì)熱帶地區(qū)天氣、氣候異常和大氣環(huán)流有直接影響,還能通過(guò)大氣遙相關(guān)的方式影響到熱帶外區(qū)域甚至全球范圍內(nèi)的氣候[2-3]。因此對(duì)ENSO的認(rèn)識(shí)一直是氣候?qū)W上的熱點(diǎn)問(wèn)題。一些研究發(fā)現(xiàn),熱帶區(qū)域的海氣相互作用和海洋動(dòng)力學(xué)作用對(duì)ENSO的發(fā)生、發(fā)展過(guò)程具有重要意義[4-5]。其中Bjerknes[1]應(yīng)用海溫異常和信風(fēng)理論解釋了ENSO的發(fā)生、發(fā)展問(wèn)題,Jin(1997)結(jié)合赤道溫躍層海溫變化和赤道風(fēng)應(yīng)力間的Sverdrup平衡關(guān)系提出充放電理論[5],再結(jié)合Bjerknes正反饋機(jī)制,從動(dòng)力學(xué)上解釋了ENSO發(fā)生、發(fā)展的物理機(jī)制以及位相轉(zhuǎn)換的過(guò)程。

        除了熱帶系統(tǒng)影響ENSO外,近些年來(lái)一些研究也表明熱帶外的氣候系統(tǒng)對(duì)ENSO的發(fā)生、發(fā)展同樣具有重要的調(diào)節(jié)作用[6-9]。作為熱帶外區(qū)域重要的氣候系統(tǒng),東亞冬季風(fēng)(East Asian Winter Monsoon, EAWM)[10]和冬季北太平洋濤動(dòng)(North Pacific Oscillation, NPO)[11]已經(jīng)被證明對(duì)ENSO的發(fā)生、發(fā)展同樣具有積極的意義[12]。EAWM是空間上呈三維結(jié)構(gòu),南北跨越熱帶、亞熱帶和溫帶,主要包含西伯利亞高壓、阿留申低壓等系統(tǒng)成員,對(duì)東亞冬季氣候起關(guān)鍵作用的環(huán)流系統(tǒng)[10,13]。Li(1989)的研究中就指出強(qiáng)EAWM可以激發(fā)赤道西太平洋的西風(fēng)異常和對(duì)流活動(dòng)異常[13];隨后穆明權(quán)等(2001)證實(shí),受EAWM影響的赤道西風(fēng)異常會(huì)向東傳播,致使赤道東太平洋出現(xiàn)正的海表溫度(SST)異常,因此伴隨著海氣相互作用進(jìn)而觸發(fā)了厄爾尼諾(El Nio)事件[14]。另一方面,NPO作為北太平洋區(qū)域內(nèi)海平面氣壓在高、低緯度間的蹺蹺板式變化的大氣遙相關(guān)模態(tài),它的變化對(duì)東亞和北美地區(qū)的大氣環(huán)流系統(tǒng)均有重要影響[15]。近年來(lái),有學(xué)者指出NPO可以通過(guò)季節(jié)足印機(jī)制來(lái)影響次年冬季ENSO的發(fā)生、發(fā)展[6-7]。即在當(dāng)年冬季時(shí),與NPO相關(guān)的大氣環(huán)流異常使得北太平洋熱帶地區(qū)在次年春季出現(xiàn)明顯的SST異常(SSTa)暖海溫,同時(shí)又反作用于大氣環(huán)流,使得赤道西太平洋區(qū)域出現(xiàn)西風(fēng)異常,異常西風(fēng)又通過(guò)觸發(fā)向東傳播的赤道暖性開(kāi)爾文波,進(jìn)一步影響了次年冬季ENSO的發(fā)生。

        通過(guò)對(duì)前人研究成果的分析,可以發(fā)現(xiàn)EAWM和NPO在北太平洋區(qū)域都呈經(jīng)向的高低壓中心分布,都屬于北半球重要的氣候系統(tǒng),且都會(huì)對(duì)ENSO的發(fā)生、發(fā)展有一定的作用。然而,EAWM和NPO對(duì)ENSO的影響是否存在差異?如果有,這些差異產(chǎn)生的原因又是什么呢?本研究將嘗試回答這些科學(xué)問(wèn)題。對(duì)于這些問(wèn)題的研究有助于加深對(duì)熱帶外氣候系統(tǒng)與ENSO之間相互關(guān)系的認(rèn)識(shí),具有重要的科學(xué)意義。

        1 材料與方法

        本研究采用了相關(guān)分析、線性回歸等統(tǒng)計(jì)方法。使用了HadISST海溫資料[16],美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心和大氣研究中心(NECP-NCAR)提供的風(fēng)場(chǎng)資料、海表面氣壓(SLP)數(shù)據(jù)、500 hPa位勢(shì)高度場(chǎng)數(shù)據(jù)[17]和降水資料[18]。時(shí)間范圍選取1948年1月至2017年12月,除了海溫資料的分辨率是1°×1°之外,其余資料分辨率均是2.5°×2.5°。

        參照Furtado等(2012)對(duì)NPO指數(shù)的定義[19],本研究將范圍在(55.0°N—72.5°N, 180.0°—140.0°W)和(15.0°N—27.5°N, 175.0°E—147.5°W)的SLP距平(SLPa)場(chǎng)的區(qū)域平均值之差定義為NPO指數(shù)。本研究使用的EAWM指數(shù)定義為(25.0°N—45.0°N, 110.0°E—145.0° E)范圍內(nèi)500 hPa位勢(shì)高度場(chǎng)的平均值[20]。本研究得到的結(jié)果對(duì)定義NPO、EAWM指數(shù)所選取區(qū)域的合理變化并不敏感。此外,使用正交經(jīng)驗(yàn)函數(shù)(EOF)方法來(lái)定義NPO和EAWM指數(shù)也可以得到類(lèi)似的結(jié)果。本研究中冬季定義為12月至次年2月。另外,我們選取了Nio3.4(5.0°S—5.0°N, 170.0°E—120.0°W)指數(shù)來(lái)表征ENSO。本研究采用T檢驗(yàn)估計(jì)相關(guān)系數(shù)和回歸系數(shù)的顯著性水平。

        2 結(jié)果與討論

        2.1 NPO和EAWM的時(shí)空特征對(duì)比分析

        首先我們對(duì)NPO和EAWM對(duì)北半球特別是東亞大氣環(huán)流的線性影響進(jìn)行了回顧。圖1分別給出了冬季NPO指數(shù)和EAWM指數(shù)回歸的東亞地區(qū)SLPa場(chǎng)的空間分布。由圖1(a)可以看出,SLPa正負(fù)異常中心呈南部為異常低壓,北部為異常高壓的南北偶極態(tài)分布,就是人們熟知的NPO模態(tài)[21]。EAWM可分為正負(fù)位相,當(dāng)EAWM正位相時(shí),在北太平洋中緯度上空是異常低壓,其北部為異常高壓[圖1(b)]。對(duì)比圖1(a)發(fā)現(xiàn)EAWM和NPO異常低壓中心在北太平洋中緯度基本重合,但EAWM的異常低壓中心位置比NPO的更靠北,位于35°N附近;而NPO和EAWM的北部異常高壓空間分布則存在較大差異,即NPO正位相時(shí)北部異常高壓范圍集中,中心位于北太平洋北部邊緣;而EAWM高壓區(qū)域橫跨歐亞大陸,中心主要位于西伯利亞區(qū)域??梢?jiàn)NPO和EAWM對(duì)北半球特別是東亞大氣環(huán)流的線性影響既有相似之處,又各有不同。圖2分別給出了冬季NPO指數(shù)、EAWM指數(shù)和Nio3.4指數(shù)的時(shí)間演變曲線,由圖可見(jiàn)3個(gè)指數(shù)均有明顯的年際變化,但EAWM還具有明顯的年代際變化特征(下文詳細(xì)做介紹),在20世紀(jì)80年代和20世紀(jì)末均出現(xiàn)了明顯的變化。這與前人的結(jié)果一致[20]。

        圖1 1948—2017年冬季NPO指數(shù)和EAWM指數(shù)回歸的SLP異常場(chǎng)的空間分布

        圖2 冬季NPO指數(shù)、EAWM指數(shù)和Nio3.4指數(shù)的年際變化(1948—2017年)

        綜上所述,冬季NPO和EAWM在時(shí)空分布上既有一定的相似性,又有明顯區(qū)別。那么他們對(duì)ENSO的影響是否也具有類(lèi)似的特征呢?

        2.2 冬季NPO和EAWM對(duì)ENSO影響的異同

        為了進(jìn)一步揭示冬季NPO和EAWM與次年冬季ENSO的關(guān)系,我們首先分別用Lanczos濾波器對(duì)3種指數(shù)進(jìn)行了7 a的高通濾波,以提取1948—2017年間3種指數(shù)的年際變化,隨后分別計(jì)算了冬季NPO指數(shù)和EAWM指數(shù)與次年冬季Nio3.4指數(shù)的相關(guān)系數(shù),分別為0.49和0.24,均超過(guò)了置信度為95%的相關(guān)性檢驗(yàn)。這表明冬季NPO、EAWM與次年冬季ENSO有著密切的聯(lián)系,這與前人的研究結(jié)果一致[6-7,20]。

        隨后我們比較了冬季NPO和EAWM指數(shù)回歸到SST異常場(chǎng)的空間季節(jié)演化特征(圖3),同時(shí)也比較了與冬季NPO和EAWM相關(guān)的SST(圖略)、降水和850 hPa風(fēng)場(chǎng)異常場(chǎng)的空間季節(jié)演化特征(圖4)。需要指出的是,由于已有研究表明ENSO具有很強(qiáng)的自相關(guān)性以及與ENSO自身相關(guān)的SSTa可以對(duì)熱帶外大氣異常產(chǎn)生影響[22],因此本部分的冬季NPO指數(shù)和EAWM指數(shù)是通過(guò)線性回歸的方式去除了同期冬季ENSO信號(hào)(Nio3.4指數(shù)表征)。由圖3可見(jiàn),冬季NPO(圖3左邊兩列)、EAWM(圖3右邊兩列)回歸得到的SSTa場(chǎng)的季節(jié)演化與Vimont(2001、2003)[6-7]研究提出的季節(jié)足印機(jī)制相似。

        2.3 NPO和EAWM影響ENSO差異性的機(jī)制分析

        由圖3(b)、(d)可見(jiàn),春季在北太平洋區(qū)域NPO和EAWM引起的SSTa都呈現(xiàn)出正-負(fù)-正的三極子海溫模態(tài)。特別地,在北太平洋中部NPO和EAWM引起的SSTa強(qiáng)度基本相當(dāng)。然而,春季與NPO和EAWM相關(guān)的北太平洋區(qū)域降水正異常以及赤道西太平洋區(qū)域異常氣旋式風(fēng)場(chǎng)和西風(fēng)異常卻都有一定差異,即NPO強(qiáng)迫出的SSTa引起的大氣響應(yīng)強(qiáng)于EAWM[圖4(b)、(d)]。相應(yīng)地,當(dāng)強(qiáng)迫出的SSTa持續(xù)到次年夏季時(shí),在赤道西太平洋區(qū)域與NPO有關(guān)的SSTa場(chǎng)也是強(qiáng)于EAWM強(qiáng)迫出來(lái)的[圖3(e)、(g)]。這種差異產(chǎn)生的影響一直隨著季節(jié)變化,傳播到冬季就以回歸場(chǎng)不同強(qiáng)度的SSTa來(lái)呈現(xiàn)了[圖3(m)、(o)]。另外,對(duì)比圖4(e)、(f)和圖4(g)、(h)的異常西風(fēng)位置可知,NPO引起的異常西風(fēng)更靠近赤道,而EAWM引起的異常西風(fēng)相對(duì)偏北,這就導(dǎo)致了NPO引起的異常西風(fēng)更容易持續(xù)和加強(qiáng),從而能更有效的影響到次年冬季ENSO發(fā)生、發(fā)展。因此EAWM在熱帶太平洋區(qū)域激發(fā)并建立起的海氣正反饋機(jī)制并不能夠像NPO一樣強(qiáng),這可能就直接導(dǎo)致了EAWM對(duì)次年冬季ENSO的影響并不如NPO那么強(qiáng)。

        圖3 1948—2017年冬季NPO和EAWM回歸的SST異常場(chǎng)的季節(jié)演化的空間分布

        接下來(lái)我們通過(guò)計(jì)算西太平洋指定區(qū)域(2.5°N—17.5°N,120.0°E—170.0°E,EWP區(qū)域),即圖4紅色矩形框內(nèi)的區(qū)域平均降水和緯向風(fēng)異常分別與去除了前冬ENSO信號(hào)的冬季NPO指數(shù)和EAWM指數(shù)的逐月交叉相關(guān)系數(shù)(圖5),量化考察NPO和EAWM對(duì)EWP區(qū)域海氣相互作用的影響。圖5實(shí)線分別表示的是降水異常與冬季NPO和EAWM指數(shù)的交叉相關(guān)系數(shù)逐月演化情況。由圖可見(jiàn),冬季NPO、EAWM指數(shù)和當(dāng)年春季EWP區(qū)域平均降水異常相關(guān)系數(shù)均較好,通過(guò)了95%的顯著性檢驗(yàn)。

        圖4 1948—2017年與冬季NPO和EAWM指數(shù)相關(guān)的降水和850 hPa風(fēng)場(chǎng)異常場(chǎng)的季節(jié)演化圖

        圖5 西北太平洋區(qū)域平均降水和緯向風(fēng)異常與去除了前冬ENSO信號(hào)的冬季NPO和EAWM指數(shù)的交叉相關(guān)圖

        隨著季節(jié)更替,相關(guān)系數(shù)均有所降低,但在當(dāng)年夏季之前與冬季NPO指數(shù)有關(guān)的相關(guān)系數(shù)均大于和EAWM有關(guān)的,且差值在春季最大。此外,虛線表示的是緯向風(fēng)場(chǎng)異常與冬季NPO和EAWM指數(shù)的交叉相關(guān)逐月演化情況。由圖可見(jiàn),冬季NPO、EAWM指數(shù)和EWP區(qū)域平均緯向風(fēng)場(chǎng)整體滿足與NPO有關(guān)的交叉相關(guān)系數(shù)大于與EAWM有關(guān)的,同時(shí)相關(guān)系數(shù)在同期時(shí)最大,隨后雖有降低但一直持續(xù)到夏季也通過(guò)了95%顯著性檢驗(yàn)。另外兩條虛線間的差值隨著季節(jié)變化而縮小,即區(qū)域平均緯向風(fēng)場(chǎng)與冬季NPO和EAWM指數(shù)的交叉相關(guān)系數(shù)差值隨著季節(jié)變化而減小。這也體現(xiàn)了隨著春季EWP西風(fēng)異常被激發(fā)起之后,該區(qū)域異常西風(fēng)持續(xù),并通過(guò)觸發(fā)赤道暖性開(kāi)爾文波的傳播使得接下來(lái)夏季熱帶中東太平洋海溫變暖[23]??梢?jiàn),圖5顯示出的結(jié)論基本與圖4(b)至(d)和(j)至(l)得到的結(jié)果一致,即EAWM對(duì)ENSO發(fā)生、發(fā)展的影響不如NPO強(qiáng)的可能原因是:前一年冬季EAWM引起的熱帶中北太平洋春季降水異常增多不及前冬NPO引起的強(qiáng),從而導(dǎo)致EWP春季西風(fēng)異常增強(qiáng)的強(qiáng)度也不及NPO引起的強(qiáng),進(jìn)而對(duì)次年冬季ENSO發(fā)生和發(fā)展的影響就弱于前冬NPO[24]。

        綜上,對(duì)比圖6(a)、(b)可知:前一年冬季NPO對(duì)次年冬季ENSO影響強(qiáng)于前冬EAWM的主要原因是前冬NPO強(qiáng)迫出的正-負(fù)-正三極子海溫模態(tài)引起的大氣響應(yīng)強(qiáng)于前冬EAWM,直觀表現(xiàn)為夏季赤道西太平洋區(qū)域降水和風(fēng)場(chǎng)強(qiáng)度及位置上的差異,這種差異又導(dǎo)致二者在熱帶太平洋區(qū)域建立的海氣正反饋機(jī)制的強(qiáng)度不同,遂直接導(dǎo)致次年冬季ENSO的發(fā)展程度不同。

        圖6 NPO和EAWM影響ENSO機(jī)制圖

        此外,由圖2(b)可知,EAWM還具有明顯的年代際變化特征。那么這種年代際變化特征是否是造成冬季EAWM對(duì)次年冬季ENSO影響弱于冬季NPO的另一個(gè)原因呢?因此我們通過(guò)滑動(dòng)相關(guān)分析對(duì)比檢測(cè)冬季NPO、EAWM和次年ENSO之間的關(guān)系。由于在這部分我們最關(guān)心的是年代際變化,所以采取的滑動(dòng)窗口寬度為23 a[20]。圖7就顯示了冬季NPO、EAWM指數(shù)與次年冬季Nio3.4指數(shù)的23 a滑動(dòng)相關(guān)系數(shù)。由圖可見(jiàn):與NPO有關(guān)的滑動(dòng)相關(guān)系數(shù)雖在20世紀(jì)70年代和20世紀(jì)末有小幅度的波動(dòng),但整體均通過(guò)90%顯著性檢驗(yàn);而與EAWM有關(guān)的滑動(dòng)相關(guān)系數(shù)僅在20世紀(jì)50~70年代和20世紀(jì)90年代通過(guò)90%的顯著性檢驗(yàn),且有大幅波動(dòng),這與前人關(guān)于EAWM在20世紀(jì)80年代和20世紀(jì)末分別有一次明顯的減弱和增強(qiáng)過(guò)程的結(jié)果相一致[20,25]。上述結(jié)果表明:從年代際尺度上來(lái)說(shuō)冬季NPO指數(shù)與次年冬季ENSO的關(guān)系明顯好于EAWM;且與NPO有關(guān)的滑動(dòng)相關(guān)系數(shù)波動(dòng)遠(yuǎn)小于EAWM,即冬季NPO與次年冬季ENSO的關(guān)系明顯比EAWM與次年冬季ENSO關(guān)系穩(wěn)定。因此在年代際尺度上,冬季NPO和EAWM的穩(wěn)定性差異也可能是造成對(duì)次年冬季ENSO的發(fā)生、發(fā)展影響不同的原因之一。

        圖7 冬季NPO指數(shù)和EAWM指數(shù)分別與次年冬季Nio3.4指數(shù)在1948—2017年間的23 a滑動(dòng)相關(guān)系數(shù)

        需要注意的是,受NPO強(qiáng)迫產(chǎn)生的維多利亞模態(tài)(VM)已被證明會(huì)對(duì)次年的ENSO產(chǎn)生影響[9],同時(shí)VM在春季北太平洋強(qiáng)迫出的海溫模態(tài)與圖3(b)、(j)相似,亦會(huì)對(duì)西北太平洋的環(huán)流場(chǎng)產(chǎn)生一定影響,同時(shí)VM對(duì)熱帶大氣氣旋等熱帶系統(tǒng)也有一定的影響[9],那么NPO和EAWM影響ENSO的差異是否與VM有關(guān)呢?這種關(guān)系又如何去驗(yàn)證或者排除?其次,本研究只從強(qiáng)度上考查了EAWM與ENSO的關(guān)系,而EAWM除了強(qiáng)弱變化外,還有路徑的變化[26]。這些問(wèn)題都值得進(jìn)一步深入研究。

        3 結(jié)論

        基于HadISST海溫資料以及NCEP月平均氣象場(chǎng)再分析資料,本研究采用回歸、相關(guān)分析等方法研究了冬季NPO和EAWM對(duì)次年冬季ENSO發(fā)生、發(fā)展影響的差異性,并分析了產(chǎn)生這種差異的可能原因。結(jié)果表明,冬季NPO和EAWM對(duì)次年冬季ENSO的影響無(wú)論是從相關(guān)性大小還是觸發(fā)的ENSO強(qiáng)度上來(lái)說(shuō)都存在明顯的差異。進(jìn)一步利用SST場(chǎng)、降水場(chǎng)以及風(fēng)場(chǎng)的季節(jié)演化分析發(fā)現(xiàn),冬季NPO引起的異常西風(fēng)更靠近赤道,而EAWM引起的異常西風(fēng)相對(duì)靠北,這就導(dǎo)致了NPO引起的異常西風(fēng)更容易持續(xù)以及加強(qiáng);另外前一年冬季EAWM引起的熱帶中北太平洋春季降水異常增多不及前冬NPO引起的強(qiáng),從而導(dǎo)致EWP春季西風(fēng)異常增強(qiáng)的強(qiáng)度也不及NPO引起的強(qiáng),進(jìn)而對(duì)次年冬季ENSO發(fā)生、發(fā)展的影響就弱于前冬NPO。也就是說(shuō)冬季NPO和EAWM在熱帶太平洋區(qū)域通過(guò)引起赤道西太平洋區(qū)域西風(fēng)異常和降水異常建立起的海氣正反饋機(jī)制在強(qiáng)度上的不同,亦是造成上述差異的主要原因之一。

        同時(shí)通過(guò)對(duì)冬季NPO、EAWM指數(shù)與次年冬季ENSO年代際關(guān)系變化的考察發(fā)現(xiàn),冬季NPO指數(shù)與次年冬季ENSO的關(guān)系明顯好于EAWM;且與NPO有關(guān)的滑動(dòng)相關(guān)系數(shù)波動(dòng)遠(yuǎn)小于EAWM,即冬季NPO與次年冬季ENSO的關(guān)系明顯比EAWM與次年冬季ENSO關(guān)系穩(wěn)定。因此在年代際尺度上,冬季NPO和EAWM的穩(wěn)定性差異也可能是造成對(duì)次年冬季ENSO的發(fā)生、發(fā)展不同的又一原因。

        致謝:感謝中國(guó)科學(xué)院大氣物理研究所鐘權(quán)加博士提出的寶貴意見(jiàn)。

        猜你喜歡
        大氣區(qū)域影響
        大氣的呵護(hù)
        軍事文摘(2023年10期)2023-06-09 09:15:06
        是什么影響了滑動(dòng)摩擦力的大小
        哪些顧慮影響擔(dān)當(dāng)?
        大氣古樸揮灑自如
        大氣、水之后,土十條來(lái)了
        擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
        關(guān)于四色猜想
        分區(qū)域
        基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
        區(qū)域
        民生周刊(2012年10期)2012-10-14 09:06:46
        久久99免费精品国产| 亚洲精品国产老熟女久久| 日本高清一区二区三区视频| 一区二区视频在线国产| 国产做无码视频在线观看| 亚洲永久无码7777kkk| 国产精品三级在线专区1| 国产精品一品二区三区| 日韩av无码一区二区三区| 无码精品a∨在线观看| 精品中文字幕久久久久久| 男女动态视频99精品| 久青草影院在线观看国产| 国产成人亚洲精品无码mp4| 香蕉视频免费在线| 亚洲精品中文字幕熟女| 丰满熟妇乱又伦精品| 国产精品黄在线观看免费软件| 国产九九在线观看播放| 国产一区资源在线播放| 97精品国产一区二区三区| 国产剧情av麻豆香蕉精品| 午夜无码亚| 久久综合久久综合久久| 国产亚洲一本大道中文在线| 国产欧美精品在线一区二区三区| 日韩精品夜色二区91久久久| 亚洲av午夜一区二区三| 永久免费av无码入口国语片| 免费国产一级特黄aa大片在线| 中文字幕一区二区三区综合网| 国产成人精品日本亚洲专区61| 亚洲丁香五月激情综合| 国产男女做爰猛烈视频网站| 一区二区三区国产高清视频| 色偷偷噜噜噜亚洲男人| 日本手机在线| 国产亚洲精品一品二品| 在线观看视频播放| 一本大道香蕉视频在线观看| 青青草在线成人免费视频|