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

    伴隨IOD型和獨(dú)立型ENSO對山東氣候年際變化的影響

    2017-12-04 08:44:27郭飛燕左衛(wèi)強(qiáng)郭飛龍曹艷察
    海洋氣象學(xué)報(bào) 2017年4期
    關(guān)鍵詞:海溫印度洋年際

    郭飛燕,左衛(wèi)強(qiáng),郭飛龍,曹艷察

    (1. 青島市氣象局,山東 青島 266003; 2. 青島市氣象災(zāi)害防御工程技術(shù)研究中心,山東 青島 266003; 3. 膠州市氣象局,山東 膠州 266318; 4. 中南林業(yè)科技大學(xué),湖南 長沙 410004; 5. 國家氣象中心,北京 100081)

    伴隨IOD型和獨(dú)立型ENSO對山東氣候年際變化的影響

    郭飛燕1,2,左衛(wèi)強(qiáng)3,郭飛龍4,曹艷察5

    (1. 青島市氣象局,山東 青島 266003; 2. 青島市氣象災(zāi)害防御工程技術(shù)研究中心,山東 青島 266003; 3. 膠州市氣象局,山東 膠州 266318; 4. 中南林業(yè)科技大學(xué),湖南 長沙 410004; 5. 國家氣象中心,北京 100081)

    使用1951年以來66 a的觀測和再分析資料,通過合成分析的方法對比分析了厄爾尼諾/拉尼娜(El Nio /La Nia)伴隨正/負(fù)印度洋偶極子(positive/negative Indian Ocean Dipole,pIOD/nIOD)發(fā)生年或獨(dú)立發(fā)生年山東夏、秋季氣溫和降水的年際變化特征,結(jié)果表明,伴隨IOD型和獨(dú)立型El Nio/La Nia對山東夏、秋季氣溫和降水的影響在強(qiáng)度、范圍、正負(fù)位相、空間型態(tài)上存在很大的差異。在氣溫方面,El Nio在pIOD的調(diào)制作用下對山東南部地區(qū)夏季氣溫年際變化的影響加強(qiáng);El Nio與pIOD伴隨發(fā)生時(shí),山東秋季氣溫較常年偏高,而獨(dú)立發(fā)生時(shí)氣溫則偏低,呈反位相變化;La Nia與nIOD伴隨發(fā)生年夏季魯西北氣溫較常年偏低,La Nia獨(dú)立發(fā)生年夏季半島東部氣溫較常年偏高,氣溫異常呈反位相變化;nIOD對La Nia的調(diào)制促進(jìn)作用有利于山東秋季氣溫較常年異常偏高;850 hPa氣溫異常與山東表面氣溫異常有很強(qiáng)的正相關(guān)關(guān)系。在降水方面,El Nio在pIOD的調(diào)制作用下容易引起山東北部地區(qū)夏季降水偏少,但會削弱其對山東中部地區(qū)秋季降水負(fù)異常的影響;La Nia在nIOD的調(diào)制作用下山東境內(nèi)降水都較常年偏多,但降水異常地域分布非常不均,魯西北降水較常年顯著偏多;獨(dú)立型La Nia更易引起魯西北西部、魯中、魯南大部分地區(qū)夏季降水偏少。850 hPa環(huán)流異常配合溫度場異常對山東夏、秋季降水異常分布有一定的影響。

    伴隨型; 獨(dú)立型; ENSO; IOD; 山東年際氣候變化

    引言

    ENSO對全球氣候有極大的影響,在ENSO發(fā)生時(shí)可以通過大氣遙相關(guān)波列即太平洋北美型(Pacific North America,PNA)影響赤道外氣候[7]。ENSO發(fā)生年夏季,熱帶西太平洋和亞洲上空的東亞-太平洋(EAP)[8]遙相關(guān)型會影響東亞阻塞高壓的年際變化和西北太平洋副熱帶高壓的位置變化[9-10],從而引起我國江淮流域的降水[11]和氣溫異常,如1954年和1998年夏季長江流域的洪災(zāi)。孫舒悅和任榮彩[12]、楚純潔等[13]均研究指出ENSO冷暖位相對中國的不同區(qū)域、不同季節(jié)的影響程度、方式均有較大的差異。

    同樣,印太暖池區(qū)的夏季降水異常也可以激發(fā)遙相關(guān)Pacific-Japan(PJ)波列,從而影響赤道外氣候異常[14]。此外,印度洋海表面溫度異常還可以影響ENSO的變化,通過ENSO間接地影響天氣氣候[15]。肖子牛等[16]通過相關(guān)分析發(fā)現(xiàn)印度洋地區(qū)東西海溫的偶極振蕩與中國 6—8 月汛期的降水有較好的相關(guān)關(guān)系。郝立生等[17]研究指出近60 a華北夏季降水出現(xiàn)減少趨勢可能與熱帶印度洋海溫升高和IOD指數(shù)升高有關(guān),其中熱帶印度洋海溫整體升高的影響更為重要。劉宣飛和袁慧珍[18]使用相關(guān)分析研究指出,僅有IOD發(fā)生時(shí),其正位相年使得中國西南地區(qū)和黃河流域的秋季降水出現(xiàn)正異常,而當(dāng)IOD與ENSO伴隨出現(xiàn)時(shí),IOD正位相和El Nio使得中國西南地區(qū)秋季降水正異常區(qū)域維持并向東擴(kuò)展,還使得黃河流域秋季降水轉(zhuǎn)為負(fù)異常。

    熱帶太平洋的ENSO事件和熱帶印度洋的IOD事件均具有較穩(wěn)定的季節(jié)鎖相特征:ENSO的冷暖位相El Nio/La Nia通常在夏季開始生成,秋季發(fā)展,冬季達(dá)到峰值[2,6,19-20],次年春季開始衰亡;IOD往往在ENSO年的夏季生成,秋季發(fā)展并達(dá)到峰值[4,6,20-22],冬季消亡。這表明在某些年份ENSO與IOD會同時(shí)出現(xiàn),兩者會在同時(shí)出現(xiàn)的夏、秋季對氣候產(chǎn)生一定的影響。2017年Fan et al.[23]研究指出,在早春爆發(fā)的ENSO往往伴隨IOD的發(fā)生?!皢渭儭盓NSO事件、“單純”IOD事件對中國夏季降水和氣溫均有顯著影響,但二者產(chǎn)生的影響作用不同[24]。山東地處中國東部,屬于中國南北氣候過渡帶,受ENSO的影響很明顯[25-26]。馬露等[27]研究分析指出,ENSO強(qiáng)度對山東省氣溫變化影響不顯著,而與降水量呈顯著負(fù)相關(guān)[28]。郭飛燕等[22]通過分類分析研究發(fā)現(xiàn),ENSO對山東地區(qū)不同季節(jié)氣溫和降水異常強(qiáng)度的影響并非嚴(yán)格反對稱,且存在顯著的空間分布不均勻性:El Nio主要影響山東南部地區(qū)夏季和冬季氣溫,La Nia幾乎可以影響山東全省夏季和秋季的氣溫,即La Nia對山東氣溫年際變化的影響要比El Nio強(qiáng),且影響范圍要廣;El Nio年夏季魯西北降水較常年顯著減少,秋季半島南部降水較常年顯著減少,La Nia年僅半島地區(qū)夏季降水較常年偏多,也就是說El Nio對山東降水年際變化的影響更明顯。但目前ENSO獨(dú)立發(fā)生年或與IOD共同發(fā)生年對山東年際氣候變化的影響研究幾乎沒有。目前在做氣候預(yù)測時(shí)僅考慮了熱帶太平洋ENSO的作用,而忽略了熱帶印度洋IOD的影響,這樣得到的結(jié)果將可能是兩類現(xiàn)象疊加后或者抵消后的結(jié)果,這樣的結(jié)果有欠合理。因此非常有必要對ENSO發(fā)生的同時(shí)伴隨或者不伴隨印度洋IOD進(jìn)行分類研究,得到更為精確和完善的結(jié)論,從而為氣候預(yù)測提供新的參考依據(jù)。

    1 資料與分析方法

    為了研究熱帶太平洋和熱帶印度洋在夏季和秋季的海表面溫度年際變化情況,本文采用了美國海洋大氣管理局(National Oceanic and Atmospheric Administration,NOAA)擴(kuò)充重建海表面溫度 (Extended Reconstructed Sea Surface Temperature, Version 4, ERSST V4)[32-33]。這一再分析資料融合了國際海洋大氣資料集(International Comprehensive Ocean-Atmosphere Data Set,release 2.5,ICOADS 2.5)、Hadley中心夜間海洋大氣氣溫?cái)?shù)據(jù)集 (Hadley Centre Nighttime Marine Air Temperature Dataset Version 2,HadNMATA V2)、浮標(biāo)以及船測資料,并改進(jìn)了再分析技術(shù)和提高了質(zhì)量控制標(biāo)準(zhǔn),因此該資料更能反映真實(shí)的海洋變化情況。此外,該資料從1854年1月至今,長時(shí)間尺度的海溫變化數(shù)據(jù)能夠較準(zhǔn)確地反映大尺度氣候變化特征。文中為了和氣象資料長度保持一致,截取了1951年1月—2017年5月的月平均海表面資料數(shù)據(jù)。同時(shí)采用了美國國家環(huán)境預(yù)報(bào)中心/國家大氣研究中心 (National Centers for Environmental Prediction/National Center for Atmospheric Research, NCEP/NCAR) 的再分析850 hPa位勢高度和氣溫資料。文中氣象觀測資料采用國家氣候中心提供的中國160個(gè)站點(diǎn)的月平均氣溫和降水資料,該資料的時(shí)間長度為1951年1月—2016年12月。

    本文研究的ENSO事件和IOD事件是發(fā)生在熱帶太平洋和印度洋的年際變化現(xiàn)象,因此所有物理量在去掉年循環(huán)之后,又經(jīng)過了帶通濾波處理,也就是說將年際變化信號之外即季節(jié)變化以下和9 a以上的其他噪音過濾掉,僅保留了4~108個(gè)月之間的信號。ENSO有東部和中部型模態(tài),為了包含所有ENSO事件,本文采用了國際上通用的Nio3.4區(qū)(170°W~120°W,5°S~5°N)指數(shù)作為判別ENSO事件的關(guān)鍵區(qū)。ENSO事件的挑選標(biāo)準(zhǔn)參考Lyon et al.[34],且與NOAA官網(wǎng)采用的標(biāo)準(zhǔn)一致,即Nio3.4區(qū)海溫距平3個(gè)月滑動平均超過+/-0.5 ℃并持續(xù)5個(gè)月以上,就認(rèn)為該年為El Nio/La Nia年。IOD指數(shù)為赤道西印度洋(50°E~70°E,10°S~10°N)和東印度洋(90~110°E,10°S~0°)區(qū)域平均的海表面溫度異常之差[4],挑選標(biāo)準(zhǔn)參考Guo et al.[6]IOD指數(shù)5個(gè)月滑動平均超過+/-0.8 ℃并持續(xù)5個(gè)月以上,就認(rèn)為該年為正/負(fù)IOD年。

    本文首先采用相關(guān)分析法說明ENSO與IOD的關(guān)系,以及ENSO、IOD指數(shù)與山東夏、秋季氣溫和降水年際變化之間的相關(guān)性。使用分類合成分析法研究ENSO獨(dú)立發(fā)生或伴隨IOD發(fā)生對山東夏、秋季氣溫和降水年際變化的影響,同時(shí)從大氣環(huán)流角度出發(fā)探究其影響成因。

    2 熱帶太平洋、印度洋海表面溫度異常時(shí)空分布特征

    基于Lyon et al.[34]定義的ENSO標(biāo)準(zhǔn),1951—2016年熱帶太平洋共發(fā)生16次El Nio事件(表1),其中有9次El Nio(1957、1963、1972、1982、1987、1991、1997、2006、2015年)發(fā)生時(shí)伴隨發(fā)生正(positive)IOD(pIOD)事件(簡稱伴隨型El Nio),7次El Nio(1953、1965、1969、1976、1986、2002、2009年)事件獨(dú)立發(fā)生(簡稱獨(dú)立型El Nio);同時(shí)在過去66 a內(nèi)共發(fā)生17次La Nia (表1)事件,其中8次La Nia(1964、1970、1984、1996、1998、2005、2010、2016年)發(fā)生時(shí)伴隨發(fā)生負(fù)(negative)IOD(nIOD)事件(簡稱伴隨型La Nia),9次La Nia(1955、1967、1973、1975、1983、1988、1995、1999、2007年)事件獨(dú)立發(fā)生(簡稱獨(dú)立型La Nia)。各類型事件樣本數(shù)相當(dāng),具有可比性。

    表11951—2016年期間ElNio/LaNia與正/負(fù)IOD伴隨發(fā)生年及其獨(dú)立發(fā)生年

    Table 1 The El Nio/ La Nia events with or without pIOD/nIOD events during 1951—2016

    類型年份ElNi?o+pIOD1957,1963,1972,1982,1987,1991,1997,2006,2015ElNi?o1953,1965,1969,1976,1986,2002,2009LaNi?a+nIOD1964,1970,1984,1996,1998,2005,2010,2016LaNi?a1955,1967,1973,1975,1983,1988,1995,1999,2007

    圖1 1951年1月—2016年12月的ENSO(a)和IOD(b)指數(shù)時(shí)間序列Fig.1 Time Series of ENSO (a) and IOD (b) indexes from January 1951 to December 2016

    為了驗(yàn)證所選個(gè)例的合理性,對熱帶太平洋和熱帶印度洋海溫異常進(jìn)行了合成分析。圖3是伴隨型和獨(dú)立型El Nio(左圖)/La Nia(右圖)的海表面溫度異常合成圖。由圖3a-b和圖3e-f可以看到,太平洋的El Nio/La Nia和印度洋的pIOD/nIOD在夏季就開始發(fā)展,秋季增強(qiáng)。由圖3c-d和圖3g-h可以看到,當(dāng)El Nio或La Nia在夏季發(fā)展,秋季增強(qiáng)的過程中,印度洋海溫異常很弱,不存在印度洋偶極子現(xiàn)象,表明本文分類結(jié)果可靠。此外,與伴隨pIOD發(fā)生的El Nio年相比,El Nio獨(dú)立發(fā)生年(圖3c-d),El Nio強(qiáng)度偏弱,且海表面溫度異常大值中心偏西,同時(shí)熱帶印度洋的海表面溫度異常也較弱;然而La Nia與nIOD伴隨發(fā)生年熱帶太平洋海溫異常較其獨(dú)立發(fā)生年較弱,且海表面溫度異常大值中心位置偏西,體現(xiàn)了ENSO正負(fù)位相和IOD正負(fù)位相的不對稱性。無論是伴隨型El Nio還是獨(dú)立型El Nio,熱帶太平洋海溫異常強(qiáng)度要比相應(yīng)的La Nia型強(qiáng),這與前人的結(jié)果一致。pIOD的強(qiáng)度也要比nIOD強(qiáng)度偏強(qiáng),范圍也較大。

    圖2 夏季(a)和秋季(b)的ENSO(實(shí)線)和IOD(虛線)指數(shù)時(shí)間序列Fig.2 Time Series of ENSO (solid line) and IOD (dash line) indexes in summer (a) and autumn (b)

    圖3 厄爾尼諾(左圖)/拉尼娜(右圖)與正/負(fù)IOD伴隨發(fā)生年(a、e為夏季, b、f為秋季) 及其獨(dú)立發(fā)生年(c、g為夏季,d、h為秋季)的季節(jié)平均海表面溫度異常合成圖Fig.3 Seasonal averaged SSTA composites for the years of El Nio(left panel)/La Nia(right panel) with pIOD/nIOD concurrence (a, e for summer; b, f for autumn) and El Nio only (c, g for summer; d, h for autumn)

    既然ENSO和IOD的正負(fù)位相存在很強(qiáng)的振幅和空間不對稱性,因此有必要對熱帶太平洋的El Nio/La Nia事件以及熱帶印度洋的pIOD/nIOD事件進(jìn)行分類研究。ENSO和IOD均具有很強(qiáng)的季節(jié)鎖相特征,ENSO在夏、秋季發(fā)展,冬季達(dá)到峰值,IOD在夏季發(fā)展,秋季達(dá)到峰值,且夏、秋季是山東省降水較多的季節(jié),因此將對山東省夏季和秋季的年際變化展開研究。

    3 ENSO、IOD指數(shù)與山東夏、秋季氣溫和降水的相關(guān)性分布特征

    對比夏季ENSO、IOD與山東氣溫(圖4a、c)相關(guān)性分析可知,兩類熱帶年際變化模態(tài)與山東氣溫的相關(guān)性并不高,這可能與兩者的正負(fù)位相對山東氣溫影響高度不對稱有關(guān)。但從有限的弱相關(guān)分布型來看,ENSO、IOD與魯西北部分地區(qū)氣溫變化具有一定的正相關(guān)性,兩者對魯西北的影響是同位相的。ENSO與山東南部地區(qū)夏季氣溫呈負(fù)相關(guān)關(guān)系(圖4a),而IOD只與膠州灣西部部分地區(qū)呈負(fù)相關(guān)關(guān)系(圖4c),可見ENSO對山東南部地區(qū)的影響較IOD強(qiáng)。秋季,ENSO以及IOD與威海和煙臺部分地區(qū)氣溫變化具有弱負(fù)相關(guān)關(guān)系,也就是說兩者對半島東部秋季氣溫的影響是一致的。但是IOD與半島南部沿岸地區(qū)氣溫有弱正相關(guān)關(guān)系。

    圖4 1951—2016年ENSO (上圖)、IOD(下圖) 指數(shù)與山東夏(a、c)、秋(b、d)季氣溫異常相關(guān)系數(shù)分布圖 Fig.4 Correlation coefficients of the ENSO (top panel) and IOD (bottom panel) indexes with the simultaneous air temperature anomalies at each point of Shandong in summer (a, c) and autumn (b, d) for the period 1951-2016

    與氣溫的相關(guān)性相比(圖4),ENSO、IOD與降水的相關(guān)性要明顯偏強(qiáng)(圖5),且ENSO與山東夏季降水相關(guān)性最強(qiáng)(圖5a)。在山東主降水期夏季,除魯西南地區(qū)外,山東大部分地區(qū)降水與ENSO指數(shù)有負(fù)相關(guān)性,尤其以魯西北地區(qū)最為顯著,相關(guān)系數(shù)可達(dá)-0.4(超過99%信度檢驗(yàn))。IOD與魯西北東部、半島西部以及魯西部分地區(qū)夏季降水有弱的負(fù)相關(guān)性,這與ENSO的影響是一致的。但I(xiàn)OD與魯中、魯西南半島東部部分地區(qū)夏季降水相關(guān)性較差(圖5c)。ENSO與山東南部秋季降水有較強(qiáng)的負(fù)相關(guān)性(膠州灣西部高達(dá)-0.3,超過95%信度檢驗(yàn)),而與魯西北相關(guān)性較差(圖5b),這樣的相關(guān)分布型結(jié)果正好與夏季相反。IOD與魯西北部分地區(qū)秋季降水呈弱正相關(guān),而與山東南部地區(qū)呈弱負(fù)相關(guān)(圖5d),IOD與ENSO對山東南部秋季降水的影響基本一致。

    圖5 1951—2016年ENSO(上圖)、IOD(下圖)指數(shù)與山東夏(a,c)、秋(b,d)季降水異常相關(guān)系數(shù)分布圖 Fig.5 Same as Fig.4, but for precipitation

    從以上相關(guān)分析可知,ENSO與IOD均對山東不同地區(qū)夏、秋季氣溫和降水有不同程度的影響,且ENSO的影響要更強(qiáng)。在不同季節(jié),ENSO與IOD對山東某些地區(qū)的影響是同位相的,但對其他一些地區(qū)的影響則不一致,甚至是相反的。前人的研究結(jié)果[43-45]表明,IOD不僅會與ENSO同時(shí)發(fā)生,還可以通過對流層準(zhǔn)兩年振蕩相互影響[46]。Annamalai et al.[47]使用觀測和數(shù)值模式資料指出,pIOD會引起赤道西太平洋出現(xiàn)西風(fēng)異常,從而有利于El Nio的維持加強(qiáng)。但是ENSO與IOD對不同地區(qū)大氣的影響可能并不是簡單線性疊加的結(jié)果。因此當(dāng)研究ENSO對山東氣候的影響時(shí)非常有必要將獨(dú)立型ENSO和伴隨IOD型ENSO進(jìn)行分類。鑒于ENSO以及IOD均具有很高的正負(fù)位相不對稱性,將分別對ENSO和IOD的正負(fù)冷暖位相進(jìn)行分類研究。

    4 ENSO獨(dú)立發(fā)生和伴隨IOD發(fā)生對山東氣溫年際變化的影響

    圖6 山東省El Nio與pIOD伴隨發(fā)生年(a為夏季,b為秋季)以及獨(dú)立發(fā)生年(c為夏季,d為秋季)的季節(jié)平均氣溫異常合成圖(等值線表示通過80%信度的顯著性檢驗(yàn))Fig.6 Seasonal averaged air temperature anomaly composites for the years of El Nio with pIOD concurrence (a for summer, b for autumn) and El Nio only (c for summer, d for autumn) in Shandong Province (contour areas are significant at the 80% confidence level)

    圖7 同圖6,但為La Nia事件的合成圖Fig.7 Same as Fig. 6, but for La Nia

    5 ENSO獨(dú)立發(fā)生和伴隨IOD發(fā)生對山東降水年際變化的影響

    圖8 山東省El Nio與pIOD伴隨發(fā)生年(a為夏季,b為秋季)以及獨(dú)立發(fā)生年(c為夏季,d為秋季)的季節(jié)平均降水異常合成圖(等值線表示通過80%信度的顯著性檢驗(yàn))Fig.8 Seasonal averaged precipitation anomaly composites for the years of El Nio with pIOD concurrence (a for summer, b for autumn) and El Nio only (c for summer, d for autumn) in Shandong Province (contour areas are significant at the 80% confidence level)

    圖9 同圖8,但為La Nia事件的合成圖Fig.9 Same as Fig. 8, but for La Nia

    6 ENSO獨(dú)立發(fā)生和伴隨IOD發(fā)生大氣環(huán)流型分析

    圖10 El Nio與pIOD伴隨發(fā)生年(a為夏季,b為秋季)以及獨(dú)立發(fā)生年(c為夏季,d為秋季)的季節(jié)平均850 hPa氣溫(彩色填圖,單位:℃)和位勢高度(等值線,單位:dagpm)異常合成圖Fig.10 Seasonal averaged 850 hPa air temperature (shaded in ℃) and geopotential height (contours in dagpm) anomaly composites for the years of El Nio with pIOD concurrence (a for summer, b for autumn) and El Nio only (c for summer, d for autumn)

    圖11 同圖10,但為La Nia事件的合成圖Fig.11 Same as Fig. 10, but for La Nia

    7 結(jié)論與討論

    本文采用1951年1月—2016年12月NOAA ERSST V4月平均海表面溫度數(shù)據(jù)、NCEP/NCAR的位勢高度和氣溫?cái)?shù)據(jù)以及國家氣候中心提供的中國160站氣溫和降水站點(diǎn)觀測資料,首先對熱帶太平洋和印度洋的年際變化現(xiàn)象進(jìn)行分析,同時(shí)使用ENSO、IOD指數(shù)與山東夏、秋季氣溫和降水做相關(guān)分析,得到ENSO、IOD各自對山東氣候的影響。并從大氣環(huán)流角度分析了造成兩類ENSO事件對山東氣候影響差異的原因。通過對66 a間發(fā)生的El Nio/La Nia與pIOD/nIOD的伴隨發(fā)生事件和El Nio/La Nia獨(dú)立發(fā)生事件進(jìn)行分類合成分析,揭示了ENSO背景下IOD發(fā)生與否對山東夏季和秋季氣溫和降水年際變化的影響。文中將ENSO和IOD事件進(jìn)行了細(xì)致分類,導(dǎo)致每一類事件樣本數(shù)較少,因此只進(jìn)行了80%的信度檢驗(yàn),但對于文中通過80%信度檢驗(yàn)的大部分地區(qū)同時(shí)通過90%甚至95%信度檢驗(yàn)(圖略),表明結(jié)果顯著可信。本文將El Nio/La Nia與pIOD/nIOD伴隨發(fā)生定義為伴隨型El Nio/La Nia,將El Nio/La Nia獨(dú)立發(fā)生定義為獨(dú)立型El Nio/La Nia。ENSO背景下IOD發(fā)生與否對山東不同季節(jié)氣溫和降水的影響在強(qiáng)度、范圍、正負(fù)位相、空間型態(tài)上存在很大的差異,結(jié)果如下:

    [1] Philander S G. El Nio, La Nia, and the Southern Oscillation[J].San Diego: Academic Press, 1990.

    [2] 劉秦玉,謝尚平,鄭小童.熱帶海洋-大氣相互作用[M].北京:高等教育出版社,2013.

    [3] 董宇佳,孟祥鳳.兩類ENSO事件分類的研究進(jìn)展[J].海洋環(huán)境科學(xué)學(xué)報(bào),2015,34(3):473-480.

    [4] Saji N H, Goswami B N, Vinayachandran P N, et al. A dipole mode in the tropical Indian Ocean [J].Nature,1999, 401(6751):360-363.

    [5] Webster P J, Moore A M, Loschnigg J P, et al. Coupled ocean-atmosphere dynamics in the Indian Ocean during 1997-98 [J].Nature,1999, 401(6751):356-360.

    [6] Guo F, Liu Q-Y, Sun S, et al. Three types of Indian ocean dipoles[J].J Climate, 2015, 28 (8): 3073-3092.

    [7] Alexander M A, Blade I, Newman M, et al. The atmospheric bridge: the influence of ENSO teleconnections on air-sea interaction over the global oceans[J].J Climate,2002,15(16): 2205-2231.

    [8] Huang R-H, Li W-J. Influence of heat source anomaly over the tropical western Pacific on the subtropical western Pacific over the East Asia[C]. Chendu: Proc International Conference on the General Circulation of East Asia,1987:40-45.

    [9] 陸日宇,黃榮輝.東亞-太平洋遙相關(guān)波列對夏季東北亞阻塞高壓年際變化的影響[J].大氣科學(xué),1998,22(5):727-734.

    [10] 陸日宇.與華北地區(qū)春季降水量異常關(guān)聯(lián)的大氣環(huán)流異常[J].氣候與環(huán)境研究,2001,6(4):400-408.

    [11] 朱益民,楊修群,陳小穎,等. ENSO與中國夏季年際氣候異常關(guān)系的年代際變化[J].熱帶氣象學(xué)報(bào),2007,23(2):105-116.

    [12] 孫舒悅,任榮彩.ENSO與印度洋海盆海溫多尺度相互作用及其對氣候影響的研究進(jìn)展 [J].氣象科技,2015,43(5):858-865.

    [13] 楚純潔,趙景波,安春華.1957—2011年中國中部不同氣候帶氣候變化及其與ENSO的關(guān)系[J].地域研究與開發(fā),34(5):121-127.

    [14] Nitta T. Convective activities in the tropical western Pacific and their impact on the northern hemisphere summer circulation[J].J Meteor Soc Japan,1987,65(3):373-390.

    [15] Yang J, Liu Q Y, Xie S-P, et al. Impact of the Indian Ocean SST basin mode on the Asian summer monsoon[J]. Geophys Res Lett,2007,34(2):155-164. DOI:10.1029/2006GL028571.

    [16] 肖子牛,晏紅明,李崇銀. 印度洋地區(qū)異常海溫的偶極振蕩與中國降水及溫度的關(guān)系[J]. 熱帶氣象學(xué)報(bào),2002,18(4):335-344.

    [17] 郝立生,丁一匯,康文英,等.印度洋海溫變化與華北夏季降水減少的關(guān)系[J].氣候變化研究快報(bào),2012,1(1):13-21.

    [18] 劉宣飛,袁慧珍. ENSO對印度洋偶極子與中國秋季降水關(guān)系的影響[J].大氣科學(xué)學(xué)報(bào), 2006,29(6):762-768.

    [19] Neelin J D, Jin F-F, Syu H-H. Variations in ENSO phase locking[J].J Climate,2000,13(14): 2570-2590.

    [20] Liu Q, Guo F, Zheng X-T. Relationships of interannual variability between the equatorial Pacific and tropical Indian Ocean in 17 CMIP5 models[J].J Ocean Univ China, 2013, 12(2): 237-244.DOI:10.1007/s11802-013-2195-8.

    [21] Zheng X-T, Xie S-P, Vecchi G A, et al. Indian Ocean dipole response to global warming: Analysis of ocean-atmospheric feedbacks in a coupled model[J].J Climate,2010,23(5): 1240-1253.

    [22] 郭飛燕,畢瑋,郭飛龍,等.山東氣候年際變化特征及其與ENSO的關(guān)系[J].海洋與湖沼,2017,48(3):465-474.

    [23] Fan L, Liu Q, Wang C, et al. Indian Ocean dipole modes associated with different types of ENSO development[J].J Climate,2017,30(6):2233-2249.DOI:10.1175/JCLI-D-16-0426.1.

    [24] 楊霞.IOD對ENSO影響亞澳地區(qū)降水和氣溫的干擾作用[D].南京:南京信息工程大學(xué),2006.

    [25] 孟翠玲,徐宗學(xué). ENSO與山東省降水的關(guān)系[J].人民黃河,2007,29(1):33-35.

    [26] 王常山,王慶,夏艷玲,等. 1961—2011年山東地區(qū)氣溫變化的區(qū)域差異特征及突變分析[J].魯東大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,29(4):353-360.

    [27] 馬露,楊東,曾婷. 1961—2012年山東省氣候變化特征及其與ENSO的關(guān)系[J].中國農(nóng)學(xué)通報(bào),2015,31(16):241-249.

    [28] 馬露,楊東,錢大文. ENSO事件對山東省區(qū)域降水量及干旱指數(shù)的影響[J].中國農(nóng)業(yè)氣象,2015(6):666-673.

    [29] Cai W, Borlace S, Lengaigne M, et al. Increasing frequency of extreme El Nino events due to greenhouse warming[J].Nature Climate Change,2014, 4(2):111-116.

    [30] Cai W, Wang G, Agus S, et al. Increased frequency of extreme La Nina events under greenhouse warming[J].Nature Climate Change, 2015, 5(2):132-137.

    [31] Cai W, Santoso A, Wang G, et al. Increased frequency of extreme Indian Ocean dipole events due to greenhouse warming[J].Nature,2014,510(7504):254-258.

    [32] Huang B, Banzon V F, Freeman E, et al. Extended Reconstructed Sea Surface Temperature version 4 (ERSST.v4): Part I. Upgrades and intercomparisons [J]. J Climate, 2015, 28(3): 911-930. DOI:10.1175/JCLI-D-14-00006.1

    [33] Liu W, Huang B, Thorne P W, et al. Extended Reconstructed Sea Surface Temperature version 4 (ERSST.v4):Part II. Parametric and structural uncertainty estimations[J].J Climate,2015,28(3):931-951.DOI:10.1175/JCLI-D-14-00007.1.

    [34] Lyon B, Cristi H, Verceles E R, et al. Seasonal reversal of the ENSO rainfall signal in the Philippines[J]. Geophys Res Lett,2006,33(24):194-199.DOI:10.1029/2006GL028182.

    [35] Meinen C S, McPhaden M J. Observations of warm water volume changes in the equatorial Pacific and their relationship to El Nio and La Nia[J].J Climate,2000,13(20):3551-3559.

    [36] Larkin N K, Harrison D E.ENSO warm (El Nio) and cold (La Nia) event life cycles: Ocean surface anomaly patterns, their symmetries, asymmetries, and implications[J].J Climate, 2002, 15(10):1118-1140.

    [37] Jin F-F, An S-I, Timmermann A, et al. Strong El Nio events and nonlinear dynamical heating [J]. Geophys Res Lett,2003,30(3):1120. DOI:10.1029/2002GL016356.

    [38] McPhaden M J, Zhang X. Asymmetry in zonal phase propagation of ENSO sea surface temperature anomalies[J].Geophys Res Lett,2009,36(13):88-97.

    [39] 張人禾.東太平洋赤道地區(qū)厄爾尼諾和拉尼娜幅度不對稱的原因[J].科學(xué)中國人,2010 (8):61.

    [40] Saji N H, Yamagata T. Structure of SST and surface wind variability during Indian Ocean dipole mode events: COADS observations[J].J Climate,2003,16(16):2735-2751.

    [41] Hong C, Li T, Lin H. Asymmetry of the Indian Ocean dipole.Part I:Observational analysis[J].J Climate,2008, 21(18):4834-4848.

    [42] Cai W, Van Rensch P, Cowan T. An asymmetry in the IOD and ENSO teleconnection pathway and its impact on Australian climate[J].J Climate,2012,25(18):6318-6329.

    [43] Shinoda T, Hendon H H, Alexander M A. Surface and subsurface dipole variability in the Indian Ocean and its relation with ENSO[J]. Deep Sea Res I, 2004, 51(5):619-635.

    [44] Izumo T, Vialard J, Lengaigne M, et al. Influence of the state of the Indian Ocean dipole on the following year’s El Nio[J].Nat Geosci,2010,3(3):168-172.

    [45] 郭飛燕.熱帶印度洋海溫年際變化主模態(tài)的分類及其與ENSO的聯(lián)系[D].青島:中國海洋大學(xué),2015.

    [46] Meehl G A, Arblaster J M, Loschnigg J. Coupled ocean-atmosphere dynamical processes in the tropical Indian and Pacific Oceans and the TBO[J].J Climate,2003,16(13):2138-2158.

    [47] Annamalai H, Xie S P, McCreary J P, et al. Impact of Indian Ocean sea surface temperature on developing El Nio[J].J Climate,2005,18(2):302-319.

    TheimpactsofENSOforaccompaniedwithIODtypeandisolatedtypeontheinterannualclimatevariabilityinShandong

    GUO Feiyan1,2, ZUO Weiqiang3, GUO Feilong4, CAO Yancha5

    (1.QingdaoMeteorologicalBureau,Qingdao266003,China; 2.QingdaoEngineeringTechnologyResearchCenterforMeteorologicalDisasterPrevention,Qingdao266003,China; 3.JiaozhouMeteorologicalBureau,Jiaozhou266318,China; 4.CentralSouthUniversityofForestryandTechnology,Changsha410004,China; 5.NationalMeteorologicalCenter,Beijing100081,China)

    Based on the observation and reanalysis datasets since 1951, the interannual variability of air temperature and precipitation in Shandong under the influence of El Nio/La Nia-Southern Oscillation (ENSO) with (accompanied type) or without (isolated type) Indian Ocean Dipole (IOD) are conducted through the composite analysis. It is discovered that there are distinct differences on the effect intensity, extent, phase and spatial pattern of interannual variability between the accompanied and isolated type of ENSO: for the air temperature, with the reinforce of positive IOD(pIOD), El Nio has more notable effect on the air temperature anomalies over the south of Shandong in summer; when El Nio occurs with pIOD, the air temperature in Shandong is higher than normal in autumn, while it is lower when El Nio occurs only, which varies out of phase; when La Nia appears with negative IOD (nIOD), the air temperature over the northwest of Shandong is lower than normal, while it is higher over the east of Shandong Peninsula when La Nia appears only, which varies out of phase too; as a joint influence of La Nia and nIOD, the air temperature over Shandong is higher than normal in autumn; a close correlation exists between the air temperature anomalies at surface and 850 hPa over Shandong during the summer and autumn of El Nio/La Nia. As for the precipitation, under the influence of pIOD, El Nio tends to reduce the precipitation of the north Shandong in summer, but weakens its influence on the negative precipitation anomalies in the middle of Shandong. La Nia occurring with nIOD has a stronger effect on the interannual variability of summer precipitation in Shandong; the precipitation displays significant inhomogenous distribution accompanying with significant positive anomalies in the northwest of Shandong; for the isolated La Nia, the summer precipitation is less than the climatology in the northwest, middle and south of Shandong; the circulation anomalies accompanied with the air temperature anomalies at 850 hPa play a role in the precipitation anomalies over Shandong in summer and autumn.

    accompanied type; isolated type; ENSO; IOD; interannual climate variability of Shandong

    郭飛燕,左衛(wèi)強(qiáng),郭飛龍,等.伴隨IOD型和獨(dú)立型ENSO對山東氣候年際變化的影響[J].海洋氣象學(xué)報(bào),2017,37(4):34-48.

    Guo Feiyan,Zuo Weiqiang,Guo Feilong,et al. The impacts of ENSO for accompanied with IOD type and isolated type on the interannual climate variability in Shandong[J]. Journal of Marine Meteorology,2017,37(4):34-48.

    10.19513/j.cnki.issn2096-3599.2017.04.005. (in Chinese)

    P46

    A

    2096-3599(2017)04-0034-15

    10.19513/j.cnki.issn2096-3599.2017.04.005

    2017-08-25;

    2017-10-15

    青島市局校合作科研基金項(xiàng)目(2015qdqxh02);淮河流域氣象開放研究基金項(xiàng)目(HRM201601);黃渤海區(qū)域科技協(xié)同創(chuàng)新基金項(xiàng)目(QYXM201614)

    郭飛燕(1986—),女,博士,工程師,主要從事海洋與大氣相互作用、氣候變化、強(qiáng)對流天氣等研究,guofeiyan01@163.com。

    猜你喜歡
    海溫印度洋年際
    基于深度學(xué)習(xí)的海溫觀測數(shù)據(jù)質(zhì)量控制應(yīng)用研究
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    與南亞高壓相聯(lián)的歐亞大陸-印度洋經(jīng)向環(huán)流
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    南印度洋偶極型海溫與中國西南地區(qū)初秋降水的關(guān)系
    去印度洋
    滇池(2018年1期)2018-01-17 21:34:30
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    波浪輸運(yùn)通量剩余量對上層海溫的影響分析
    印度洋上接“嫦娥”回家:遠(yuǎn)望3號船精確測控探月三期試驗(yàn)返回器順利再入返回側(cè)記
    太空探索(2014年12期)2014-07-12 15:17:10
    一级黄片播放器| 亚洲第一av免费看| 国产国语露脸激情在线看| 午夜福利影视在线免费观看| 国产精品一区二区在线观看99| 亚洲国产最新在线播放| 国产 一区精品| 丁香六月欧美| 成年人免费黄色播放视频| 久热爱精品视频在线9| 高清黄色对白视频在线免费看| 亚洲五月色婷婷综合| 亚洲国产精品999| 大片免费播放器 马上看| 欧美成人精品欧美一级黄| 午夜av观看不卡| 亚洲国产精品国产精品| 人妻人人澡人人爽人人| 国产精品av久久久久免费| 最近2019中文字幕mv第一页| 国产一区二区在线观看av| 欧美日本中文国产一区发布| 一级黄片播放器| 日本猛色少妇xxxxx猛交久久| 欧美日韩综合久久久久久| 亚洲欧美激情在线| 成年女人毛片免费观看观看9 | 老汉色av国产亚洲站长工具| 亚洲,一卡二卡三卡| 亚洲国产av影院在线观看| 十八禁网站网址无遮挡| 伊人久久国产一区二区| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 激情视频va一区二区三区| 亚洲av男天堂| 国产乱人偷精品视频| 十分钟在线观看高清视频www| 女人久久www免费人成看片| 水蜜桃什么品种好| 久久人人爽av亚洲精品天堂| 国产成人精品久久久久久| 国产成人系列免费观看| 青春草亚洲视频在线观看| 亚洲天堂av无毛| 亚洲av中文av极速乱| 五月开心婷婷网| 精品亚洲乱码少妇综合久久| a级毛片黄视频| 男女国产视频网站| 国产欧美日韩综合在线一区二区| 看免费成人av毛片| 无遮挡黄片免费观看| 中文乱码字字幕精品一区二区三区| 国产视频首页在线观看| 国产黄频视频在线观看| 久久女婷五月综合色啪小说| 久久久欧美国产精品| 黄色视频在线播放观看不卡| 久久久久视频综合| 日本wwww免费看| 欧美人与性动交α欧美精品济南到| 国产av码专区亚洲av| 看免费成人av毛片| 国产欧美日韩一区二区三区在线| 卡戴珊不雅视频在线播放| 天堂8中文在线网| 国产男女内射视频| 国产成人精品久久二区二区91 | 波多野结衣av一区二区av| 我的亚洲天堂| 成人毛片60女人毛片免费| 在线精品无人区一区二区三| 啦啦啦啦在线视频资源| 久久精品久久精品一区二区三区| 如日韩欧美国产精品一区二区三区| 久久狼人影院| 久久久久久久国产电影| 久久99热这里只频精品6学生| 日本av免费视频播放| 精品酒店卫生间| 老熟女久久久| 精品一区二区三区av网在线观看 | 日韩不卡一区二区三区视频在线| av视频免费观看在线观看| 黄网站色视频无遮挡免费观看| 精品一区二区三区av网在线观看 | 日韩人妻精品一区2区三区| 国产在线一区二区三区精| 久久综合国产亚洲精品| 桃花免费在线播放| 黄频高清免费视频| 欧美激情 高清一区二区三区| 午夜精品国产一区二区电影| 久久久久人妻精品一区果冻| 欧美日韩视频精品一区| 一级毛片黄色毛片免费观看视频| 久热这里只有精品99| 成年人午夜在线观看视频| 国产有黄有色有爽视频| 啦啦啦啦在线视频资源| 日本wwww免费看| 一区二区日韩欧美中文字幕| 男女国产视频网站| 日韩熟女老妇一区二区性免费视频| 日韩免费高清中文字幕av| 欧美激情高清一区二区三区 | 久久久久网色| 丁香六月欧美| 日韩制服骚丝袜av| 在线天堂最新版资源| 精品人妻在线不人妻| 水蜜桃什么品种好| 啦啦啦中文免费视频观看日本| 免费不卡黄色视频| 肉色欧美久久久久久久蜜桃| 国产av一区二区精品久久| 国产精品蜜桃在线观看| 啦啦啦在线免费观看视频4| 午夜老司机福利片| 99精国产麻豆久久婷婷| 精品人妻在线不人妻| 精品少妇一区二区三区视频日本电影 | 十分钟在线观看高清视频www| 十八禁网站网址无遮挡| 纯流量卡能插随身wifi吗| 日本欧美国产在线视频| 亚洲欧洲日产国产| 哪个播放器可以免费观看大片| 香蕉丝袜av| 在线观看人妻少妇| 99热全是精品| 伦理电影免费视频| 永久免费av网站大全| 亚洲av电影在线进入| 熟妇人妻不卡中文字幕| 国产无遮挡羞羞视频在线观看| av女优亚洲男人天堂| 丰满乱子伦码专区| 大码成人一级视频| 91aial.com中文字幕在线观看| 欧美激情 高清一区二区三区| 97人妻天天添夜夜摸| 看免费av毛片| xxx大片免费视频| 免费观看a级毛片全部| 免费人妻精品一区二区三区视频| 国产熟女欧美一区二区| 日韩精品免费视频一区二区三区| www.熟女人妻精品国产| 在线观看免费午夜福利视频| 亚洲人成网站在线观看播放| 国产精品99久久99久久久不卡 | 国产成人啪精品午夜网站| 男女下面插进去视频免费观看| 两个人免费观看高清视频| 性色av一级| 我要看黄色一级片免费的| 夫妻午夜视频| 日韩不卡一区二区三区视频在线| 欧美日韩成人在线一区二区| 久久国产精品大桥未久av| 三上悠亚av全集在线观看| 啦啦啦在线免费观看视频4| 亚洲av欧美aⅴ国产| 国产在线视频一区二区| www.自偷自拍.com| 丝袜美足系列| 国产精品亚洲av一区麻豆 | 青草久久国产| 亚洲欧美一区二区三区黑人| 久久婷婷青草| 亚洲av综合色区一区| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区久久| 国产精品秋霞免费鲁丝片| 在线天堂最新版资源| 亚洲成av片中文字幕在线观看| 亚洲第一av免费看| 男人舔女人的私密视频| 嫩草影院入口| 国产1区2区3区精品| 99热国产这里只有精品6| 熟女av电影| 婷婷成人精品国产| 亚洲欧美清纯卡通| 精品久久蜜臀av无| 国产日韩欧美在线精品| 悠悠久久av| 久久ye,这里只有精品| 久久久精品国产亚洲av高清涩受| 国产免费福利视频在线观看| 久久鲁丝午夜福利片| 男女无遮挡免费网站观看| 国产日韩欧美视频二区| 午夜福利影视在线免费观看| 久久精品久久精品一区二区三区| 最近最新中文字幕免费大全7| 99香蕉大伊视频| 热99久久久久精品小说推荐| 国产精品久久久久久久久免| 久久久久久久大尺度免费视频| 日韩制服丝袜自拍偷拍| 一边摸一边做爽爽视频免费| 国产精品免费视频内射| av片东京热男人的天堂| 纵有疾风起免费观看全集完整版| 侵犯人妻中文字幕一二三四区| 黄片播放在线免费| 亚洲精品自拍成人| 欧美激情高清一区二区三区 | 午夜激情av网站| 黄色怎么调成土黄色| 人人澡人人妻人| 毛片一级片免费看久久久久| 欧美乱码精品一区二区三区| 9191精品国产免费久久| 精品国产国语对白av| 在现免费观看毛片| 欧美日韩视频高清一区二区三区二| 各种免费的搞黄视频| 亚洲精品aⅴ在线观看| 午夜福利,免费看| 亚洲五月色婷婷综合| 我的亚洲天堂| 女人精品久久久久毛片| 这个男人来自地球电影免费观看 | 色精品久久人妻99蜜桃| 国产成人精品在线电影| 国产精品国产av在线观看| 国产xxxxx性猛交| 建设人人有责人人尽责人人享有的| 亚洲自偷自拍图片 自拍| xxx大片免费视频| 一区二区三区精品91| 黄片无遮挡物在线观看| 久久精品久久久久久噜噜老黄| 中文欧美无线码| 99香蕉大伊视频| 久久鲁丝午夜福利片| 赤兔流量卡办理| 国产毛片在线视频| 中文字幕人妻丝袜一区二区 | 男女午夜视频在线观看| av女优亚洲男人天堂| 欧美精品亚洲一区二区| 亚洲 欧美一区二区三区| 精品一区二区三区av网在线观看 | 欧美激情高清一区二区三区 | 久久久久精品性色| 老熟女久久久| 人妻一区二区av| 波野结衣二区三区在线| 人成视频在线观看免费观看| 黄频高清免费视频| 国产成人精品在线电影| 99re6热这里在线精品视频| 欧美黄色片欧美黄色片| 国产淫语在线视频| 午夜影院在线不卡| 91成人精品电影| 另类精品久久| 午夜免费鲁丝| 欧美人与性动交α欧美精品济南到| 免费黄频网站在线观看国产| 嫩草影院入口| 亚洲欧美精品自产自拍| 男女午夜视频在线观看| 激情五月婷婷亚洲| 最新的欧美精品一区二区| 亚洲视频免费观看视频| 色视频在线一区二区三区| 黑人猛操日本美女一级片| 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 国产无遮挡羞羞视频在线观看| 在线免费观看不下载黄p国产| 久久久亚洲精品成人影院| 电影成人av| 亚洲第一青青草原| 亚洲国产欧美日韩在线播放| 日韩 亚洲 欧美在线| 久久精品亚洲熟妇少妇任你| 国产精品亚洲av一区麻豆 | 如日韩欧美国产精品一区二区三区| a级毛片在线看网站| 又大又爽又粗| 视频在线观看一区二区三区| 国产黄频视频在线观看| 久久久精品国产亚洲av高清涩受| 欧美日韩福利视频一区二区| 视频区图区小说| 国产成人欧美在线观看 | 亚洲国产看品久久| www.自偷自拍.com| 自线自在国产av| 久久久久久久大尺度免费视频| 天堂俺去俺来也www色官网| av女优亚洲男人天堂| 我要看黄色一级片免费的| 一区二区三区激情视频| 国产黄色视频一区二区在线观看| 如何舔出高潮| 免费看不卡的av| 久久久久精品人妻al黑| 男人添女人高潮全过程视频| 十八禁高潮呻吟视频| 纯流量卡能插随身wifi吗| 老汉色∧v一级毛片| 一区二区三区四区激情视频| 青草久久国产| 最近最新中文字幕大全免费视频 | 亚洲国产中文字幕在线视频| 中文字幕人妻熟女乱码| 男人操女人黄网站| 制服诱惑二区| 大香蕉久久成人网| 亚洲一区二区三区欧美精品| 男女下面插进去视频免费观看| 国产成人精品在线电影| 欧美激情 高清一区二区三区| 欧美黑人精品巨大| 高清在线视频一区二区三区| 高清在线视频一区二区三区| 欧美在线黄色| 精品亚洲成a人片在线观看| 亚洲美女搞黄在线观看| 亚洲欧美色中文字幕在线| 不卡视频在线观看欧美| 亚洲激情五月婷婷啪啪| av又黄又爽大尺度在线免费看| 人人妻人人添人人爽欧美一区卜| 国产免费视频播放在线视频| 一边摸一边做爽爽视频免费| 伊人久久大香线蕉亚洲五| avwww免费| 亚洲成色77777| 波多野结衣一区麻豆| 久久久国产一区二区| 9191精品国产免费久久| 国产精品国产三级国产专区5o| 国产免费又黄又爽又色| 夫妻性生交免费视频一级片| 国产 一区精品| 中国国产av一级| 狂野欧美激情性bbbbbb| tube8黄色片| 天天添夜夜摸| 中文天堂在线官网| 爱豆传媒免费全集在线观看| 精品卡一卡二卡四卡免费| 亚洲国产欧美日韩在线播放| 51午夜福利影视在线观看| 免费观看a级毛片全部| av天堂久久9| 免费日韩欧美在线观看| 亚洲专区中文字幕在线 | 日韩中文字幕视频在线看片| 9191精品国产免费久久| 亚洲精品中文字幕在线视频| 久久久国产欧美日韩av| 成人亚洲欧美一区二区av| 成人国产麻豆网| 中文字幕av电影在线播放| 亚洲精品国产区一区二| 人妻 亚洲 视频| 欧美日韩视频精品一区| 黑人欧美特级aaaaaa片| 下体分泌物呈黄色| 80岁老熟妇乱子伦牲交| 欧美日韩一区二区视频在线观看视频在线| 亚洲成国产人片在线观看| 午夜91福利影院| 无遮挡黄片免费观看| 丝袜美腿诱惑在线| 国产精品人妻久久久影院| 亚洲欧美精品综合一区二区三区| 女人高潮潮喷娇喘18禁视频| www.av在线官网国产| 婷婷色综合大香蕉| 99久久综合免费| 国产亚洲av片在线观看秒播厂| 国产福利在线免费观看视频| 国产精品无大码| 国产精品三级大全| 少妇 在线观看| 免费看不卡的av| 国产极品天堂在线| 综合色丁香网| 一边摸一边做爽爽视频免费| 国产深夜福利视频在线观看| 亚洲,欧美精品.| av在线观看视频网站免费| xxx大片免费视频| 国产精品免费视频内射| 久久狼人影院| 十八禁人妻一区二区| 午夜久久久在线观看| 韩国av在线不卡| 久久精品国产亚洲av涩爱| 日韩不卡一区二区三区视频在线| av国产久精品久网站免费入址| 国产野战对白在线观看| 国产精品一区二区在线观看99| 国产麻豆69| 美女福利国产在线| 老熟女久久久| 一级毛片黄色毛片免费观看视频| 久久久久视频综合| 亚洲精品第二区| 久久av网站| 伦理电影大哥的女人| 日韩不卡一区二区三区视频在线| 一边亲一边摸免费视频| 免费黄色在线免费观看| 精品少妇黑人巨大在线播放| 99久久精品国产亚洲精品| 亚洲欧美一区二区三区久久| 久久狼人影院| 日韩av不卡免费在线播放| 久久久久精品国产欧美久久久 | 久久久久久久国产电影| 成年美女黄网站色视频大全免费| 美女主播在线视频| 欧美日韩视频高清一区二区三区二| 国产欧美日韩一区二区三区在线| 久久久久久久久久久免费av| 大码成人一级视频| 国产精品免费视频内射| 中文字幕人妻丝袜一区二区 | 亚洲精品成人av观看孕妇| www.av在线官网国产| 久久久久国产一级毛片高清牌| 高清不卡的av网站| 国产精品久久久av美女十八| 亚洲一码二码三码区别大吗| 欧美国产精品一级二级三级| 欧美 亚洲 国产 日韩一| 国产不卡av网站在线观看| 亚洲美女黄色视频免费看| 美女高潮到喷水免费观看| 在线观看免费视频网站a站| 两性夫妻黄色片| 9191精品国产免费久久| 久久精品aⅴ一区二区三区四区| 满18在线观看网站| 校园人妻丝袜中文字幕| 男女午夜视频在线观看| 天天操日日干夜夜撸| 啦啦啦中文免费视频观看日本| 国产精品麻豆人妻色哟哟久久| 国产精品久久久av美女十八| 精品福利永久在线观看| 青春草亚洲视频在线观看| 亚洲图色成人| 色精品久久人妻99蜜桃| 青春草亚洲视频在线观看| 波多野结衣av一区二区av| 18在线观看网站| 国产成人a∨麻豆精品| 精品一品国产午夜福利视频| 婷婷色综合大香蕉| av女优亚洲男人天堂| av有码第一页| 国产精品欧美亚洲77777| 伊人久久国产一区二区| 亚洲精品视频女| 中文字幕高清在线视频| 在线观看人妻少妇| 18在线观看网站| 国产成人免费无遮挡视频| 精品一区二区三卡| 91国产中文字幕| 毛片一级片免费看久久久久| 日韩欧美一区视频在线观看| 又大又爽又粗| 亚洲,一卡二卡三卡| 欧美在线黄色| 汤姆久久久久久久影院中文字幕| 中文字幕高清在线视频| 一级爰片在线观看| 国产精品成人在线| 亚洲综合精品二区| 国产又色又爽无遮挡免| 视频区图区小说| 亚洲欧美一区二区三区国产| 女性生殖器流出的白浆| 日本猛色少妇xxxxx猛交久久| 老司机靠b影院| 人人妻人人澡人人爽人人夜夜| 美女扒开内裤让男人捅视频| 久久久精品区二区三区| videosex国产| 高清av免费在线| 亚洲精品国产av蜜桃| 亚洲av欧美aⅴ国产| 美女福利国产在线| 另类亚洲欧美激情| 性少妇av在线| 狂野欧美激情性bbbbbb| 亚洲人成77777在线视频| 国产亚洲最大av| 午夜影院在线不卡| 最黄视频免费看| 韩国av在线不卡| 成年美女黄网站色视频大全免费| 亚洲av电影在线观看一区二区三区| 一个人免费看片子| 别揉我奶头~嗯~啊~动态视频 | 人人妻人人爽人人添夜夜欢视频| 99久久精品国产亚洲精品| 一二三四在线观看免费中文在| 亚洲av电影在线进入| 一级毛片电影观看| videos熟女内射| 99香蕉大伊视频| 国产精品99久久99久久久不卡 | 人人澡人人妻人| 日韩一区二区视频免费看| 亚洲av男天堂| 在线亚洲精品国产二区图片欧美| 男女下面插进去视频免费观看| xxxhd国产人妻xxx| 精品一品国产午夜福利视频| 国产成人免费观看mmmm| 中文字幕亚洲精品专区| 丝袜人妻中文字幕| 日韩免费高清中文字幕av| 国产欧美亚洲国产| 久久久久精品国产欧美久久久 | 美女高潮到喷水免费观看| 亚洲精品乱久久久久久| 久久久久久久大尺度免费视频| 超色免费av| 精品久久久精品久久久| 欧美激情极品国产一区二区三区| √禁漫天堂资源中文www| 国产亚洲av片在线观看秒播厂| bbb黄色大片| 日韩 亚洲 欧美在线| 免费观看av网站的网址| 美女福利国产在线| 亚洲,欧美,日韩| 欧美国产精品一级二级三级| 韩国精品一区二区三区| 午夜日本视频在线| 色婷婷久久久亚洲欧美| 国产高清不卡午夜福利| 看免费av毛片| 成人午夜精彩视频在线观看| 桃花免费在线播放| 18禁观看日本| 国产黄色免费在线视频| av女优亚洲男人天堂| 香蕉丝袜av| av国产精品久久久久影院| 日韩欧美一区视频在线观看| 少妇人妻 视频| 欧美成人精品欧美一级黄| 日韩大码丰满熟妇| 啦啦啦视频在线资源免费观看| 熟女av电影| 男女高潮啪啪啪动态图| 日本91视频免费播放| 国产精品国产三级专区第一集| 丝袜脚勾引网站| 免费女性裸体啪啪无遮挡网站| 人人澡人人妻人| 精品一品国产午夜福利视频| 中文字幕色久视频| 久久人人97超碰香蕉20202| 日韩大码丰满熟妇| 在现免费观看毛片| 黄频高清免费视频| 热re99久久国产66热| 国产一区有黄有色的免费视频| av片东京热男人的天堂| 老汉色∧v一级毛片| 视频在线观看一区二区三区| av电影中文网址| 国产精品.久久久| 日韩欧美精品免费久久| 国产成人精品久久久久久| 日韩 亚洲 欧美在线| 肉色欧美久久久久久久蜜桃| 欧美乱码精品一区二区三区| 啦啦啦 在线观看视频| 成人影院久久| 国产成人午夜福利电影在线观看| √禁漫天堂资源中文www| 国产精品香港三级国产av潘金莲 | 美女主播在线视频| 悠悠久久av| 99精品久久久久人妻精品| 免费看av在线观看网站| 欧美精品一区二区免费开放| 亚洲精品在线美女| 亚洲精品自拍成人| 亚洲欧美日韩另类电影网站| 午夜日本视频在线| 亚洲三区欧美一区| 久久天躁狠狠躁夜夜2o2o | 亚洲综合精品二区| 免费观看av网站的网址| 亚洲精品日本国产第一区| 国产精品 国内视频| 国产黄色视频一区二区在线观看| 电影成人av| 日韩 亚洲 欧美在线| 日韩大码丰满熟妇| 国产黄色视频一区二区在线观看| 国产男女超爽视频在线观看| 91精品伊人久久大香线蕉| 午夜福利免费观看在线| 卡戴珊不雅视频在线播放|