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

    CMIP5模式對(duì)拉尼娜生命史模擬能力的評(píng)估

    2017-05-30 10:48:04黃玉蓉張福穎陳明誠
    大氣科學(xué)學(xué)報(bào) 2017年3期

    黃玉蓉 張福穎 陳明誠

    摘要利用19個(gè)CMIP5模式輸出資料,評(píng)估模式對(duì)于拉尼娜事件特殊生命史發(fā)展過程的模擬能力。評(píng)估結(jié)果顯示,僅有少數(shù)模式可以很好地再現(xiàn)拉尼娜事件緩慢衰減并再次增強(qiáng)的生命史發(fā)展過程,而多數(shù)模式中拉尼娜事件持續(xù)衰減直至消亡。觀測分析結(jié)果表明,一個(gè)可能導(dǎo)致拉尼娜再次增強(qiáng)的原因是風(fēng)場強(qiáng)迫作用下的海洋赤道波動(dòng)過程。模擬能力較好的模式可以建立起“SST一對(duì)流一風(fēng)場”正反饋過程,使得拉尼娜事件再次發(fā)展。而模擬能力較弱的模式中正反饋過程無法建立,因此拉尼娜事件最終消亡。另一個(gè)可能導(dǎo)致拉尼娜事件再次增強(qiáng)的原因是海洋平均經(jīng)圈環(huán)流的作用。模擬能力較好的模式可以很好地模擬出氣候態(tài)海洋經(jīng)圈環(huán)流強(qiáng)度,因此海洋平均經(jīng)向冷平流會(huì)幫助赤道地區(qū)負(fù)海溫距平再次增強(qiáng)。而模擬能力較弱的模式中海洋經(jīng)圈環(huán)流強(qiáng)度較弱,因此赤道地區(qū)負(fù)海溫距平持續(xù)衰減,最終回歸到氣候態(tài)。

    關(guān)鍵詞海氣間相互作用;模式評(píng)估;ENSO動(dòng)力學(xué);生命史不對(duì)稱性;熱量收支診斷分析

    厄爾尼諾一南方濤動(dòng)(ENSO)事件是地球氣候系統(tǒng)中最顯著的年際變率(Rasmusson andCarpenter,1982;Philander,1990)。作為熱帶地區(qū)年際變率的主導(dǎo)模態(tài),ENSO對(duì)全球氣候變化都有重要的影響作用(Latif et a1.,1998;Trenberth et a1.,1998;Wallace et a1.,1998;Alexander et a1.,2002;張人禾和巢紀(jì)平,2002;Jin et a1.,2016)。在厄爾尼諾事件發(fā)展階段,不穩(wěn)定的大氣一海洋間相互作用過程(Bjerknes,1969;Philander et a1.,1984;Cane andZebiak,1985;Hirst,1986,1988),以及海洋赤道波動(dòng)過程(Suarez and Schopf,1988;Battisti and Hirst,1989;Jin,1997;Li,1997)將使得海溫正距平不斷增長,發(fā)展年底時(shí)達(dá)到成熟時(shí)期。在此之后,厄爾尼諾事件在西北太平洋異常反氣旋性環(huán)流強(qiáng)迫作用下迅速衰減,發(fā)展年次年底時(shí)形成拉尼娜事件(Wang eta1.,2000,2003;Wu et a1.,2010;Chen et a1.,2016)。

    作為厄爾尼諾事件的負(fù)位相性事件,拉尼娜事件的生命史演變過程具有顯著地不對(duì)稱性特征(Kang and Kug,2002;Kessler,2002;Larkin and Har-rison,2002;McPhaden et a1.,2009;Okumura and De-ser,2010;Dommenget et a1.,2013)。拉尼娜事件在發(fā)展年底達(dá)到最強(qiáng)值之后往往緩慢的衰減,次年底時(shí)又再次增強(qiáng)形成第二次拉尼娜事件。Chen et a1.(2016)的診斷分析結(jié)果表明,西北太平洋上大氣低層風(fēng)場空間分布的不對(duì)稱性以及赤道東太平洋地區(qū)短波輻射和潛熱釋放等熱力通量場空間分布的不對(duì)稱性,都對(duì)厄爾尼諾和拉尼娜事件生命史不對(duì)稱性的形成具有重要作用。

    耦合模式比較計(jì)劃(CMIP)是在大氣模式比較計(jì)劃(AMIP)基礎(chǔ)之上發(fā)展而來。作為目前包含內(nèi)容最為全面的氣候變化模式資料集,CMIP5模式結(jié)果是我們研究過往氣候變化特征和預(yù)估未來氣候變化的重要依據(jù)之一(張學(xué)洪等,2003;張雅樂等,2012;杜美芳等,2015;智海等,2015)。然而受限于氣候變化特征的復(fù)雜性以及模式自身的局限性,CMIP5模式對(duì)于ENSO事件生命史特征的模擬仍存在很多不足之處,特別是對(duì)于拉尼娜事件特殊生命史發(fā)展過程的模擬情況具有顯著的不確定性,并呈現(xiàn)出兩種不同的發(fā)展特征(圖1)。例如以CESMl-CAM5模式(圖1f)為代表的第一類CMIP5模式可以較好地還原觀測中拉尼娜事件真實(shí)的生命史發(fā)展過程(圖1a)。在發(fā)展年次年秋季時(shí),拉尼娜事件具有再次增長過程。然而在以ACCESS1.0模式(圖1b)為代表的大部分CMIP5模式模擬結(jié)果中,拉尼娜事件與觀測結(jié)果表現(xiàn)出較強(qiáng)的不對(duì)稱性發(fā)展特征,在發(fā)展年次年底時(shí)并沒有再次增長為第二次拉尼娜事件,而是持續(xù)衰減直至回復(fù)到氣候態(tài)。

    在本文研究工作中,將通過海洋次表層海溫?zé)崃渴罩г\斷分析方法,定量分析兩類CMIP5模式中拉尼娜事件生命史不對(duì)稱性的形成原因。首先提出如下幾個(gè)科學(xué)問題:在兩類CMIP5模式中,拉尼娜事件生命史不對(duì)稱性形成的前期信號(hào)是什么?導(dǎo)致拉尼娜事件生命史不對(duì)稱性形成的主要?jiǎng)恿σ蜃佑心男扛鶕?jù)上述提出的科學(xué)問題,設(shè)計(jì)論文結(jié)構(gòu)如下:第1節(jié)簡要介紹文章中使用的CMIP5多模式輸出資料和觀測的海洋一大氣數(shù)據(jù)集,以及海洋次表層海溫?zé)崃渴罩г\斷分析方法;第2節(jié)對(duì)兩類CMIP5模式中拉尼娜事件生命史不對(duì)稱性的形成機(jī)理進(jìn)行研究;最后是本文的結(jié)論部分。

    1資料與方法

    文章中為了評(píng)估模式對(duì)于拉尼娜事件生命史發(fā)展過程的模擬能力,選取了19個(gè)CMIP5氣候模式的歷史模擬試驗(yàn)數(shù)據(jù),并選取相應(yīng)的拉尼娜事件進(jìn)行合成分析。文章中拉尼娜事件的選取標(biāo)準(zhǔn)與NO.AA ONI(Oceanic NINO Index)指數(shù)選取標(biāo)準(zhǔn)一致,當(dāng)3個(gè)月滑動(dòng)平均的Nifio-3.4區(qū)海洋表面溫度(sea Surface Temperature,SST)距平連續(xù)5個(gè)月小于-0.5 K時(shí),將其定義為一次拉尼娜事件。CMIP5模式資料研究時(shí)段統(tǒng)一選取為1950年1月至1999年12月,共計(jì)50 a。由于各模式資料間水平分辨率存在較大差異,將所有模式數(shù)據(jù)結(jié)果統(tǒng)一插值為l?!?°。相關(guān)模式基本信息請見表1,更多詳細(xì)內(nèi)容請參閱http://cmip-pcmdi.llnl.gov/cmip5/。

    研究中使用的觀測的海洋再分析資料集包括美國馬里蘭大學(xué)提供的簡單海洋同化資料(SO-DAv2.1.6;Carton and Giese,2008),以及歐洲中期天氣預(yù)報(bào)中心提供的海洋再分析資料系統(tǒng)(ORAS4;Balmaseda et a1.,2013)。緯向風(fēng)應(yīng)力數(shù)據(jù)同樣來自于SODAv2.1.6和ORAS4海洋再分析資料。SST資料來自于美國國家海洋和大氣局提供的延伸重建海表面溫度資料(ERSSTv3b;Smith et a1.,2008)。向外長波輻射場資料同樣由美國國家海洋和大氣局提供(Outgoing Longwave Radiation,OLR;Liebmann and Smith,1996)。

    為了深入理解導(dǎo)致兩類CMIP5模式中拉尼娜事件生命史不對(duì)稱性形成的動(dòng)力學(xué)機(jī)理,我們利用上述海洋再分析資料集對(duì)海洋次表層海溫?zé)崃渴罩н^程進(jìn)行診斷分析。海溫變化傾向診斷方程可表示為:上述方程中各項(xiàng)計(jì)算結(jié)果均為混合層內(nèi)垂直平均計(jì)算得到。由于文章中所使用的觀測數(shù)據(jù)資料分別具有不同的時(shí)問長度,為了使合成結(jié)果具有一致性,統(tǒng)一選取了自1980年起的時(shí)間段進(jìn)行合成分析。診斷方程中海洋混合層海溫傾向項(xiàng)和三維動(dòng)力平流項(xiàng)結(jié)果均由SODAv2.1.6和ORAS4資料集合平均得到。

    2兩類CMIP5模式中拉尼娜事件生命史不對(duì)稱性形成機(jī)制

    2.1拉尼娜事件中生命史不對(duì)稱性發(fā)展特征

    相關(guān)觀測分析結(jié)果表明,相較于厄爾尼諾事件而言,拉尼娜事件通常具有更長的生命史發(fā)展周期,往往可以持續(xù)2 a或3 a以上的時(shí)間(Chen et a1.,2016)。在拉尼娜事件衰減過程中存在著一個(gè)十分有趣的現(xiàn)象,當(dāng)負(fù)海溫距平發(fā)展至冬季到達(dá)最強(qiáng)值之后通常會(huì)進(jìn)入緩慢衰減階段。然而在發(fā)展年次年后半年中,負(fù)海溫距平將會(huì)再次增強(qiáng)形成第二次拉尼娜事件(圖1a)。因此在拉尼娜事件發(fā)展年次年,負(fù)海溫距平場經(jīng)歷了緩慢衰減和再次增強(qiáng)兩種生命史發(fā)展過程。但是在CMIP5多模式模擬結(jié)果中,拉尼娜事件生命史發(fā)展過程卻具有顯著的不確定性,并呈現(xiàn)出兩種截然不同的發(fā)展特征(圖1)。

    在第一類模式中,模式模擬結(jié)果可以較好地再現(xiàn)出觀測中拉尼娜事件的生命史發(fā)展特征。以CESM1-CAM5模式模擬結(jié)果為例(圖1f),拉尼娜事件自發(fā)展年1月起由較強(qiáng)的正海溫距平開始衰減,6月時(shí)回到氣候平均態(tài),之后海溫繼續(xù)降低,在發(fā)展年12月時(shí)負(fù)海溫距平達(dá)到最強(qiáng),最大值約為-1.3℃,與觀測中拉尼娜事件成熟時(shí)期振幅相近。自發(fā)展年次年1月起,負(fù)海溫距平開始緩慢衰減,7月時(shí)依然維持著負(fù)海溫距平,振幅約為-0.6℃,振幅強(qiáng)度約為最強(qiáng)值時(shí)期1/2。在后續(xù)秋季時(shí)期,負(fù)海溫距平再次增長,年底12月時(shí)負(fù)海溫距平超過-1℃,與前期拉尼娜事件成熟時(shí)期振幅強(qiáng)度相當(dāng)。

    然而在第二類模式中,模式模擬結(jié)果相比于觀測結(jié)果呈現(xiàn)出不對(duì)稱性生命史發(fā)展特征。以AC-CESS1.0模式模擬結(jié)果為例(圖lb),拉尼娜事件同樣自發(fā)展年1月起由正海溫距平開始衰減,表明拉尼娜事件同樣由厄爾尼諾事件轉(zhuǎn)變而來,但是正海溫距平強(qiáng)度明顯弱于觀測結(jié)果中。正海溫距平在4月時(shí)衰減至氣候平均態(tài),之后轉(zhuǎn)變?yōu)樨?fù)海溫距平開始增長,在發(fā)展年底12月時(shí)到達(dá)成熟期。成熟期最強(qiáng)振幅強(qiáng)度為-1℃,與觀測結(jié)果相近??傮w而言,在發(fā)展年時(shí)期CMIP5多模式模擬的拉尼娜事件發(fā)展過程與觀測結(jié)果保持了較好的一致性,不對(duì)稱性生命史特征主要出現(xiàn)在發(fā)展年次年衰減階段。在以ACCESS1.0模式為例的第二類CMIP5模式中,拉尼娜事件自發(fā)展年次年1月起開始緩慢衰減,7月時(shí)同樣維持著負(fù)海溫距平,振幅強(qiáng)度約為峰值時(shí)期1/2。然而在后續(xù)秋季時(shí)期,拉尼娜事件并沒有再次增強(qiáng)形成第二次拉尼娜事件,而是持續(xù)著緩慢的衰減過程,最終在發(fā)展年次年底時(shí)海溫回歸到氣候平均態(tài)。

    進(jìn)一步分析拉尼娜事件中赤道東太平洋區(qū)域海溫距平變化傾向隨時(shí)問的演變特征(圖2)可以發(fā)現(xiàn),在觀測結(jié)果和CMIP5多模式模擬結(jié)果中,拉尼娜事件發(fā)展年時(shí)期海溫變率均以負(fù)距平為主。表明在這一時(shí)期赤道東太平洋區(qū)域海溫持續(xù)衰減,拉尼娜事件為發(fā)展時(shí)期。在發(fā)展年次年衰減時(shí)期,CMIP5多模式模擬結(jié)果中海溫變化的不對(duì)稱性特征開始顯現(xiàn)。在以CESM1-CAM5模式(圖2f)為代表的第一類模式模擬結(jié)果中,發(fā)展年次年前期(1-6月)海溫距平變率以正距平為主,表明拉尼娜正處于衰減階段。發(fā)展年次年后期夏秋季節(jié)中,海溫距平變率迅速由正距平轉(zhuǎn)變?yōu)樨?fù)距平并得以維持,表明拉尼娜事件生命史發(fā)展過程再次發(fā)生轉(zhuǎn)變,進(jìn)入再次增強(qiáng)時(shí)期。這一類CMIP5模式模擬的海溫距平變率發(fā)展過程與觀測結(jié)果基本保持一致(圖2a)。

    相反的,在以ACCESS1.0模式(圖2b)為代表的第二類模式模擬結(jié)果中,在發(fā)展年次年全年中海溫距平變率均以正距平為主,前期強(qiáng)度較強(qiáng),后期強(qiáng)度減弱。這表明在此類CMIP5模式模擬的拉尼娜事件中,負(fù)海溫距平始終處于衰減階段,衰減速率先急后緩。這一類CMIP5模式模擬的海溫距平變率隨時(shí)間的演變特征與觀測結(jié)果相比呈現(xiàn)出明顯的不對(duì)稱性特征(圖2a)。

    綜上所述,根據(jù)拉尼娜事件合成結(jié)果中生命史發(fā)展過程的不同,可以將上述19個(gè)CMIP5模式模擬結(jié)果分為兩類。在第一類CMIP5模式中(包括CESM1-BGC,CESMl-CAM5和MPI-ESM-MR三個(gè)模式),拉尼娜事件的生命史發(fā)展過程與觀測結(jié)果(圖3a)相似,在發(fā)展年底時(shí)達(dá)到最強(qiáng),最強(qiáng)負(fù)海溫距平中心位于150°w的赤道中太平洋地區(qū),之后負(fù)海溫距平緩慢衰減,并在發(fā)展年次年7月后開始再次增強(qiáng),最終在次年底時(shí)發(fā)展形成第二次拉尼娜事件(圖3b)。將這一類CMIP5模式歸類為對(duì)拉尼娜事件生命史發(fā)展過程模擬較好的GM(GoodModel)模式。在第二類CMIP5模式中(包含剩余的16個(gè)CMIP5模式),拉尼娜事件的生命史發(fā)展過程與觀測結(jié)果表現(xiàn)出明顯的不對(duì)稱性(圖3c)。在發(fā)展年底成熟時(shí)期,負(fù)海溫距平中心相比于觀測結(jié)果向東偏移至120°w的赤道東太平洋地區(qū)。在此之后負(fù)海溫距平開始緩慢衰減,但是并沒有出現(xiàn)再次增長過程。在發(fā)展年次年年底時(shí),海溫回歸到氣候平均態(tài)附近。將這一類CMIP5模式歸類為對(duì)拉尼娜事件生命史發(fā)展過程模擬能力較弱的PM(Pool Model)模式。

    2.2海洋赤道波動(dòng)過程對(duì)風(fēng)場強(qiáng)迫的響應(yīng)

    在拉尼娜事件發(fā)展年次年時(shí)期,一個(gè)可能導(dǎo)致冷海溫再次增強(qiáng)的原因是風(fēng)場強(qiáng)迫作用下的海洋赤道波動(dòng)過程(Chen et a1.,2016)。根據(jù)經(jīng)典的海洋延遲振子理論(Suarez and Schopf,1988;Battisti andHirst,1989),拉尼娜事件發(fā)展年時(shí)期赤道外地區(qū)由于Rossby波活動(dòng)產(chǎn)生的向西運(yùn)動(dòng)的次表層暖海溫信號(hào)將會(huì)在西北太平洋區(qū)域集聚,并通過熱帶太平洋西邊界的反射作用轉(zhuǎn)變成沿赤道向東傳播的下沉Kelvin波從而帶來緯向異常暖平流。因此在拉尼娜事件衰減初期(發(fā)展年次年1-6月),赤道東太平洋地區(qū)負(fù)海溫距平不斷衰減(圖3a)。與此同時(shí),在拉尼娜事件發(fā)展年冬季成熟時(shí)期和后續(xù)春季衰減時(shí)期中,西北太平洋上空以反氣旋性環(huán)流為主,反氣旋性環(huán)流南側(cè)赤道西太平洋為東風(fēng)距平。東風(fēng)距平將會(huì)激發(fā)出海洋中的上翻Kelvin波向東傳播至赤道東太平洋地區(qū),從而阻礙拉尼娜事件衰減。這兩種相反的赤道海洋波動(dòng)過程將使得拉尼娜事件衰減過程中衰減速率緩慢,因此在發(fā)展年次年夏季時(shí)期赤道東太平洋地區(qū)海溫依舊維持負(fù)距平(Chen et a1.,2016)。由于北半球秋季是海氣相互作用最劇烈的時(shí)期(Li and Philander,1996),因此負(fù)海溫距平將會(huì)激發(fā)出更強(qiáng)的赤道東風(fēng)距平(圖4a)。更強(qiáng)的異常東風(fēng)環(huán)流可以加強(qiáng)赤道東太平洋地區(qū)海洋上翻運(yùn)動(dòng)以及海洋表層蒸發(fā)作用,使得負(fù)海溫距平進(jìn)一步增強(qiáng),形成“SST一對(duì)流一風(fēng)場”正反饋過程(Lj,1997;Chen et a1.,2016),因此拉尼娜事件再次發(fā)展增強(qiáng)。

    對(duì)拉尼娜事件生命史發(fā)展過程模擬能力較強(qiáng)的GM模式可以較好地再現(xiàn)出這一正反饋過程(圖4b)。在拉尼娜事件發(fā)展年次年初期(1-6月),赤道西太平洋地區(qū)有較強(qiáng)的東風(fēng)距平出現(xiàn),因此對(duì)拉尼娜事件的衰減過程起到延緩作用,在這一時(shí)期拉尼娜事件衰減緩慢。當(dāng)發(fā)展年次年后期(7-12月)熱帶東太平洋地區(qū)海氣相互作用增強(qiáng)時(shí),赤道中太平洋區(qū)域的東風(fēng)距平將會(huì)再次增強(qiáng),并通過“SST-對(duì)流一風(fēng)場”正反饋?zhàn)饔脦椭崮仁录俅卧鰪?qiáng)。

    相反的,在模擬能力較弱的PM模式模擬結(jié)果中,衰減年初期赤道西太平洋區(qū)域東風(fēng)距平強(qiáng)度相比于觀測結(jié)果明顯偏弱(圖4c),這將無法抑制初期時(shí)拉尼娜事件衰減過程。在發(fā)展年次年后期,雖然赤道東太平洋地區(qū)依然存在較弱的負(fù)海溫距平(圖3c),但是在赤道中太平洋區(qū)域東風(fēng)距平并沒有再次增強(qiáng),而是逐漸減弱直至消亡(圖4c)。因此與觀測結(jié)果相比,PM模式中并未能建立起有效的海洋一大氣問正反饋過程幫助拉尼娜事件再次增強(qiáng),最終赤道東太平洋區(qū)域海溫回歸至氣候態(tài)。

    作為“SST一對(duì)流一風(fēng)場”正反饋過程中至關(guān)重要的環(huán)節(jié),對(duì)流活動(dòng)的生成和發(fā)展過程構(gòu)建起了海溫場與大氣風(fēng)場間的聯(lián)系。因此,赤道太平洋區(qū)域?qū)α骰顒?dòng)隨時(shí)問的演變情況也可以很好地反映出拉尼娜事件的生命史發(fā)展過程(圖5)。合成分析結(jié)果表明(圖5b),模擬能力較好的GM模式與觀測結(jié)果保持了較好的一致性。在拉尼娜事件發(fā)展年次年初期,赤道中太平洋區(qū)域的負(fù)海溫距平使得局地對(duì)流活動(dòng)受到抑制,因此OLR場表現(xiàn)為較強(qiáng)的正距平值(圖5a、5b),對(duì)流活動(dòng)西側(cè)存在較強(qiáng)的東風(fēng)距平(圖4a、4b)。此后伴隨著拉尼娜事件的不斷衰減,赤道東太平洋地區(qū)負(fù)海溫距平強(qiáng)度不斷減弱,對(duì)流活動(dòng)逐漸增強(qiáng),因此OLR正距平值逐漸減弱。在發(fā)展年次年后期,由于氣候背景場下海氣相互作用增強(qiáng),“SST一對(duì)流一風(fēng)場”正反饋過程建立。因此對(duì)流活動(dòng)再次被削弱,具體表現(xiàn)為OLR正距平再次增強(qiáng)(圖5a、5b)。相應(yīng)的,對(duì)流活動(dòng)西側(cè)東風(fēng)距平再次增強(qiáng)(圖4a、4b)。

    與此不同的是,在模擬能力較弱的PM模式中,拉尼娜事件發(fā)展年次年初期赤道中太平洋地區(qū)OLR場正距平強(qiáng)度較弱(圖5c),意味著在這一時(shí)期抑制相對(duì)流較弱。考慮到在拉尼娜事件發(fā)展成熟時(shí)期,GM模式和PM模式與觀測結(jié)果相比負(fù)海溫距平強(qiáng)度相當(dāng)(圖3),因此PM模式對(duì)赤道中太平洋地區(qū)對(duì)流活動(dòng)的模擬能力明顯較弱。造成PM模式中對(duì)流活動(dòng)模擬效果較差的可能原因是在拉尼娜事件成熟時(shí)期,對(duì)負(fù)海溫距平中心位置的模擬存在偏差(圖3c)。在PM模式模擬的拉尼娜事件發(fā)展年底成熟時(shí)期,負(fù)海溫距平中心相比于觀測結(jié)果向東偏移至120°w的赤道東太平洋地區(qū)。在氣候態(tài)冷舌區(qū)背景場作用下,赤道東太平洋區(qū)域?qū)α骰顒?dòng)相對(duì)較少。因此即使給定相同強(qiáng)度的冷海溫異常值,冷舌區(qū)抑制相對(duì)流活動(dòng)相比于赤道中太平洋地區(qū)也會(huì)更弱。

    PM模式中對(duì)流活動(dòng)模擬的偏差導(dǎo)致對(duì)流西側(cè)東風(fēng)距平強(qiáng)度較弱(圖4c),此后伴隨著拉尼娜事件不斷衰減,正OLR距平(圖5c)和東風(fēng)距平(圖4c)也不斷減弱。在發(fā)展年次年后期,PM模式中“SST一對(duì)流一風(fēng)場”正反饋過程并未建立,具體表現(xiàn)為赤道太平洋地區(qū)OLR場正距平值(圖5c)和東風(fēng)距平值(圖4c)并沒有再次增強(qiáng),而是持續(xù)衰減直至氣候平均態(tài),相應(yīng)的海表溫度也回歸到平衡態(tài)(圖3c)。由此可見,PM模式對(duì)于北半球秋季時(shí)期冷舌區(qū)海氣相互作用過程的模擬能力同樣較弱。

    2.3氣候態(tài)STC環(huán)流圈經(jīng)向平流作用

    與海洋表層海溫信號(hào)相類似的,海洋次表層信號(hào)同樣可以清晰的表征出拉尼娜事件發(fā)展年次年中先“緩慢衰減”后“再次增強(qiáng)”的特殊生命史發(fā)展特征。圖6給出了合成的觀測結(jié)果和19個(gè)CMIP5模式中,拉尼娜事件發(fā)展年次年赤道太平洋地區(qū)溫躍層厚度距平隨時(shí)間的變化情況。從結(jié)果中可以清晰地看到拉尼娜事件發(fā)展年次年初期,在熱帶太平洋西邊界反射而來的暖性Kelvin波東傳的作用之下,赤道東太平洋地區(qū)原本淺薄的溫躍層厚度逐漸加深,拉尼娜事件進(jìn)入衰減階段(圖6a)。有趣的是在這一時(shí)期,原本位于赤道地區(qū)淺薄的溫躍層信號(hào)并沒有隨之消失,而是向兩極方向傳播至赤道外地區(qū),具體表現(xiàn)為在發(fā)展年次年5-7月期間,熱帶東太平洋南北緯5度附近形成了兩個(gè)負(fù)溫躍層厚度距平中心。在拉尼娜事件發(fā)展年次年后期,赤道東太平洋地區(qū)不斷加強(qiáng)的異常東風(fēng)通過加強(qiáng)海洋上翻運(yùn)動(dòng)從而抬升溫躍層高度,使得赤道地區(qū)溫躍層厚度再次變得淺薄。與此同時(shí),還能看到位于赤道外的負(fù)溫躍層厚度距平中心逐漸向赤道地區(qū)移動(dòng),最終回歸到赤道東太平洋地區(qū),幫助拉尼娜事件再次發(fā)展(圖6a)。

    GM模式可以較好的再現(xiàn)海洋次表層負(fù)溫躍層厚度距平“分裂一再回歸”的發(fā)展過程,模擬結(jié)果與觀測結(jié)果基本保持一致(圖6b),稍有不同的是分裂而來的南北兩支負(fù)溫躍層厚度距平中心向兩極方向傳播至更遠(yuǎn)的區(qū)域(南北緯10°附近)。這也表明,即使在模擬拉尼娜事件效果較好的CMIP5模式中,依然存在著許多可以改進(jìn)的地方。相對(duì)而言,在PM模式中海洋次表層負(fù)溫躍層厚度距平信號(hào)并不能再現(xiàn)“分裂一再回歸”的生命史發(fā)展過程(圖6c)。在拉尼娜事件發(fā)展年次年初期,赤道東太平洋地區(qū)負(fù)溫躍層厚度距平逐漸減弱,并向赤道外地區(qū)傳播,拉尼娜事件進(jìn)入衰減階段。在發(fā)展年次年5月時(shí),赤道地區(qū)轉(zhuǎn)變?yōu)檎郎剀S層厚度距平。位于赤道外的兩支負(fù)溫躍層厚度距平始終未能回歸至赤道地區(qū),因此赤道地區(qū)正溫躍層距平信號(hào)持續(xù)至發(fā)展年次年底,相應(yīng)的拉尼娜事件最終消亡(圖3c)。

    Chen et a1.(2016)提出氣候平均態(tài)subtropicalcell(STC;McCreary and Lu,1994)環(huán)流圈的經(jīng)向平流作用可以導(dǎo)致拉尼娜事件衰減時(shí)期次表層海溫信號(hào)“分裂一再回歸”現(xiàn)象的產(chǎn)生。通常將熱帶太平洋地區(qū)氣候平均態(tài)的經(jīng)向環(huán)流過程稱之為STC環(huán)流圈,具體表現(xiàn)為海洋次表層區(qū)向赤道的經(jīng)向洋流,赤道區(qū)平均垂直上翻運(yùn)動(dòng)以及海洋表層向極方向的經(jīng)向洋流運(yùn)動(dòng)(圖7a)。在拉尼娜事件發(fā)展年次年后期,氣候態(tài)STC環(huán)流圈中次表層區(qū)向赤道的平均經(jīng)向洋流會(huì)將赤道外區(qū)更強(qiáng)的冷信號(hào)平流至赤道地區(qū)(圖6a),并通過垂直上翻運(yùn)動(dòng)輸送至海洋表層,幫助拉尼娜事件在這一時(shí)期再次增強(qiáng)。

    在GM模式中氣候態(tài)海洋STC環(huán)流圈的模擬效果出色(圖7b),次表層區(qū)向赤道的平均經(jīng)向洋流運(yùn)動(dòng)速度可達(dá)到0.04 m/s(約1(。)/mon),強(qiáng)度與觀測結(jié)果中相當(dāng)(圖7a)。在次表層平均經(jīng)向冷平流的作用下,赤道外的負(fù)溫躍層厚度距平逐漸向赤道區(qū)移動(dòng)(圖6b)。然而在PM模式的模擬結(jié)果中,次表層區(qū)向赤道的平均經(jīng)向洋流運(yùn)動(dòng)速度約為0.02 m/s,僅為觀測結(jié)果中強(qiáng)度的1/2(圖7c)。由于PM模式中對(duì)氣候態(tài)STC環(huán)流圈模擬能力較弱,因此在溫躍層變化過程中并沒有出現(xiàn)“分裂一再回歸”的發(fā)展過程(圖6c)。

    為了進(jìn)一步檢驗(yàn)CMIP5模式對(duì)“STC環(huán)流圈經(jīng)向平流”機(jī)制的表現(xiàn)能力,對(duì)熱帶東太平洋區(qū)域次表層海溫進(jìn)行熱量收支診斷分析。圖8給出的是合成的觀測結(jié)果和19個(gè)CMIP5模式模擬結(jié)果中,拉尼娜事件發(fā)展年次年熱帶太平洋地區(qū)次表層海溫距平隨時(shí)間的演變特征。與溫躍層距平演變過程相類似的(圖6),在觀測結(jié)果和GM模式模擬結(jié)果中都具有明顯的赤道外負(fù)海溫距平逐漸向赤道地區(qū)移動(dòng)的過程(圖8a、8b)。但是在PM模式模擬結(jié)果中,赤道外區(qū)負(fù)海溫距平始終停留在原地,并沒有再次向赤道地區(qū)傳播。因此在發(fā)展年次年5月時(shí),PM模式中赤道地區(qū)次表層海溫距平轉(zhuǎn)變?yōu)檎盘?hào)并得以維持(圖8c)。

    相應(yīng)的,觀測結(jié)果和GM模式模擬結(jié)果中同樣顯示出赤道地區(qū)次表層海溫變化傾向經(jīng)歷了由正轉(zhuǎn)負(fù)的過程(圖9a、9b)。在拉尼娜事件發(fā)展年次年6-7月間,次表層海溫由不斷增暖轉(zhuǎn)變?yōu)樵俅巫兝?。相?duì)而言,GM模式中赤道地區(qū)負(fù)海溫變率的經(jīng)向尺度更大。這也許是由于赤道外區(qū)次表層冷信號(hào)距離赤道更遠(yuǎn)(圖6b、8b),經(jīng)向平均冷平流輸送距離同樣更遠(yuǎn)所導(dǎo)致。在模擬能力較弱的PM模式模擬結(jié)果中,次表層海溫變率也體現(xiàn)出由正轉(zhuǎn)負(fù)的生命史發(fā)展過程(圖9c)。但此時(shí)負(fù)海溫變率遠(yuǎn)小于觀測結(jié)果和GM模式中,并不足以使拉尼娜事件再次增長。

    海洋次表層海溫?zé)崃渴罩г\斷分析結(jié)果進(jìn)一步表明,經(jīng)向平均平流項(xiàng)是導(dǎo)致拉尼娜事件在發(fā)展年次年中,次表層海溫變率由正轉(zhuǎn)負(fù)的主要?jiǎng)恿σ蜃?。在觀測結(jié)果中,赤道外南北緯3度附近(特別是南支)有較強(qiáng)的經(jīng)向平均冷平流作用(圖10a)。GM模式模擬結(jié)果與觀測結(jié)果保持了較好的一致性,在拉尼娜事件發(fā)展年次年3-9月期間,赤道外區(qū)經(jīng)向平均冷平流作用顯著(圖10b)。而PM模式對(duì)“STC環(huán)流圈經(jīng)向平流”過程的模擬能力較弱,赤道外區(qū)經(jīng)向平均冷平流在發(fā)展年次年后期迅速衰減,無法幫助赤道地區(qū)拉尼娜事件再次發(fā)展

    3結(jié)論

    利用19個(gè)CMIP5模式輸出資料,針對(duì)拉尼娜事件中特殊生命史發(fā)展過程模擬能力進(jìn)行評(píng)估。評(píng)估結(jié)果顯示,僅有少數(shù)CMIP5模式(包括CESMl-BGC,CESM1-CAM5和MPI-ESM-MR三個(gè)模式)可以很好地再現(xiàn)拉尼娜事件發(fā)展年次年中負(fù)海溫距平場緩慢衰減和再次增強(qiáng)兩種不同的生命史發(fā)展特征,而多數(shù)CMIP5模式(包含剩余的16個(gè)CMIP5模式)中負(fù)海溫距平場持續(xù)衰減過程直至回歸到氣候平均態(tài)。因此,兩類CMIP5模式模擬結(jié)果中拉尼娜事件生命史不對(duì)稱性特征主要出現(xiàn)在發(fā)展年次年后期,模擬能力較好的GM模式中拉尼娜事件再次增強(qiáng)形成第二次拉尼娜事件,而模擬能力較弱的PM模式中拉尼娜事件持續(xù)衰減最終消亡。

    一個(gè)可能導(dǎo)致冷海溫再次增強(qiáng)的原因是風(fēng)場強(qiáng)迫作用下的海洋赤道波動(dòng)過程。在北半球秋季時(shí)期逐漸增強(qiáng)的海氣相互作用背景場作用下,GM模式中初始時(shí)刻較弱的負(fù)海溫距平可以激發(fā)出更強(qiáng)的抑制相對(duì)流活動(dòng)。抑制相對(duì)流西側(cè)東風(fēng)距平可以加強(qiáng)氣候態(tài)東風(fēng)信風(fēng),使得赤道東太平洋地區(qū)海洋上翻運(yùn)動(dòng)以及海洋表層蒸發(fā)作用加強(qiáng)。因此負(fù)海溫距平進(jìn)一步增強(qiáng),形成“SST一對(duì)流一風(fēng)場”正反饋過程,拉尼娜事件再次發(fā)展。相反的,由于PM模式對(duì)北半球秋季時(shí)期冷舌區(qū)海氣相互作用過程的模擬能力較弱,因此“SST一對(duì)流一風(fēng)場”正反饋過程并未建立,拉尼娜事件不斷減弱直至回歸到氣候平均態(tài)。

    另一個(gè)可能導(dǎo)致拉尼娜事件再次增強(qiáng)的原因是“STC環(huán)流圈經(jīng)向平流”的作用。在拉尼娜事件發(fā)展年次年中,赤道東太平洋地區(qū)次表層海溫經(jīng)歷了“分裂一再回歸”的生命史發(fā)展過程。在發(fā)展年次年后期,氣候平均態(tài)STC環(huán)流圈次表層區(qū)向赤道的平均經(jīng)向洋流運(yùn)動(dòng)會(huì)將赤道外地區(qū)更強(qiáng)的冷信號(hào)平流至赤道地區(qū),并通過赤道地區(qū)垂直上翻運(yùn)動(dòng)輸送至海洋表層,幫助拉尼娜事件在這一時(shí)期再次增強(qiáng)。GM模式對(duì)氣候態(tài)STC環(huán)流圈的模擬效果出色,強(qiáng)度與觀測結(jié)果中相當(dāng)。次表層海溫診斷分析結(jié)果同樣表明,經(jīng)向平均平流項(xiàng)是導(dǎo)致次表層海溫變化傾向由正向負(fù)轉(zhuǎn)變的主要?jiǎng)恿σ蜃?。相反的,由于PM模式對(duì)氣候態(tài)STC環(huán)流圈模擬能力較弱,強(qiáng)度僅為觀測結(jié)果中的1/2,因此次表層海溫并沒有形成“分裂一再回歸”的發(fā)展過程。赤道東太平洋地區(qū)次表層海溫診斷分析結(jié)果進(jìn)一步表明,這一時(shí)期PM模式中次表層海溫變率和經(jīng)向平均平流項(xiàng)強(qiáng)度遠(yuǎn)小于觀測結(jié)果。因此PM模式中拉尼娜事件并沒有再次增強(qiáng)形成第二次拉尼娜事件,而是不斷衰減直至消亡。

    国产精品久久久久久亚洲av鲁大| 久久精品国产清高在天天线| 国产探花在线观看一区二区| 中文在线观看免费www的网站 | 夜夜躁狠狠躁天天躁| 亚洲av中文字字幕乱码综合| 中文字幕高清在线视频| 99热这里只有是精品50| 精品午夜福利视频在线观看一区| 午夜免费激情av| 大型黄色视频在线免费观看| 午夜福利高清视频| 巨乳人妻的诱惑在线观看| 日韩欧美精品v在线| 欧美黑人精品巨大| 在线永久观看黄色视频| 色噜噜av男人的天堂激情| 国产精品免费一区二区三区在线| 国产成年人精品一区二区| 淫秽高清视频在线观看| 中文字幕熟女人妻在线| 欧美日韩乱码在线| 成人手机av| 男人的好看免费观看在线视频 | 国产亚洲av高清不卡| 国产又色又爽无遮挡免费看| 舔av片在线| 欧美日韩一级在线毛片| 成人欧美大片| 亚洲专区中文字幕在线| 听说在线观看完整版免费高清| 欧美日韩一级在线毛片| 欧美乱码精品一区二区三区| 久久国产乱子伦精品免费另类| 久久久国产精品麻豆| 成人av在线播放网站| 欧美一级a爱片免费观看看 | 国产69精品久久久久777片 | av欧美777| 男女之事视频高清在线观看| 欧美中文综合在线视频| 成人高潮视频无遮挡免费网站| 久久久久亚洲av毛片大全| 亚洲国产欧美人成| 日本一区二区免费在线视频| 搡老熟女国产l中国老女人| 欧美+亚洲+日韩+国产| 床上黄色一级片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 小说图片视频综合网站| 此物有八面人人有两片| 久久久国产精品麻豆| 18禁黄网站禁片午夜丰满| 国产1区2区3区精品| 亚洲全国av大片| 国产亚洲欧美在线一区二区| 免费在线观看黄色视频的| 精品国产乱子伦一区二区三区| 久久久久久免费高清国产稀缺| 久久久久久久久久黄片| 高清毛片免费观看视频网站| 人妻丰满熟妇av一区二区三区| 妹子高潮喷水视频| 久久久久久国产a免费观看| 国产精品美女特级片免费视频播放器 | 精品欧美一区二区三区在线| 欧美日韩国产亚洲二区| 一级a爱片免费观看的视频| 黄频高清免费视频| 亚洲欧洲精品一区二区精品久久久| 欧美日韩瑟瑟在线播放| 国产精品一及| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 此物有八面人人有两片| 两个人视频免费观看高清| 色综合欧美亚洲国产小说| 日本熟妇午夜| 国产蜜桃级精品一区二区三区| 色精品久久人妻99蜜桃| 成人18禁高潮啪啪吃奶动态图| 黑人巨大精品欧美一区二区mp4| 欧美午夜高清在线| 禁无遮挡网站| 中文字幕高清在线视频| 日韩精品免费视频一区二区三区| 男人的好看免费观看在线视频 | 制服丝袜大香蕉在线| 欧美中文日本在线观看视频| 婷婷亚洲欧美| 国产91精品成人一区二区三区| 男插女下体视频免费在线播放| 男女那种视频在线观看| 亚洲 国产 在线| 90打野战视频偷拍视频| 国产精品亚洲av一区麻豆| 亚洲自偷自拍图片 自拍| 啦啦啦免费观看视频1| 日韩av在线大香蕉| 国产精品 欧美亚洲| 在线观看免费视频日本深夜| 岛国视频午夜一区免费看| 久久久国产成人免费| 真人做人爱边吃奶动态| 美女大奶头视频| 欧美在线一区亚洲| 精品久久久久久久久久免费视频| 神马国产精品三级电影在线观看 | 精品不卡国产一区二区三区| 岛国视频午夜一区免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机深夜福利视频在线观看| 18禁黄网站禁片免费观看直播| 亚洲国产精品成人综合色| 欧美另类亚洲清纯唯美| 国产一级毛片七仙女欲春2| 久久精品影院6| 18禁黄网站禁片午夜丰满| 国产亚洲精品一区二区www| 亚洲国产精品999在线| 成人午夜高清在线视频| 丰满人妻熟妇乱又伦精品不卡| 国产视频一区二区在线看| 国产精品一区二区三区四区久久| 午夜老司机福利片| 欧美性长视频在线观看| 亚洲男人天堂网一区| 久9热在线精品视频| 中文字幕高清在线视频| 成人三级黄色视频| 午夜福利视频1000在线观看| 亚洲成av人片在线播放无| 97人妻精品一区二区三区麻豆| 久久亚洲真实| 淫妇啪啪啪对白视频| 97碰自拍视频| 精品久久久久久久末码| 啦啦啦韩国在线观看视频| 黄色 视频免费看| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲人成网站高清观看| 又爽又黄无遮挡网站| 夜夜夜夜夜久久久久| 国产片内射在线| 亚洲午夜理论影院| 一级黄色大片毛片| 一级a爱片免费观看的视频| 国产精品九九99| 国产精品久久视频播放| 亚洲成a人片在线一区二区| 啦啦啦观看免费观看视频高清| 亚洲成av人片在线播放无| 在线观看www视频免费| 搡老岳熟女国产| 久久久久久久久久黄片| 亚洲av熟女| 操出白浆在线播放| 熟女少妇亚洲综合色aaa.| 久久国产精品人妻蜜桃| www日本在线高清视频| 午夜福利高清视频| 美女午夜性视频免费| 99在线人妻在线中文字幕| 老鸭窝网址在线观看| 免费在线观看成人毛片| 久久人妻av系列| 成人高潮视频无遮挡免费网站| 国产三级中文精品| 久久精品国产亚洲av香蕉五月| 亚洲18禁久久av| 国内毛片毛片毛片毛片毛片| 18禁裸乳无遮挡免费网站照片| 99久久综合精品五月天人人| 亚洲国产欧美网| 午夜视频精品福利| 国产精品久久久久久精品电影| www日本黄色视频网| 久久久精品欧美日韩精品| 日韩高清综合在线| 国产日本99.免费观看| 国产1区2区3区精品| 欧美日韩黄片免| 精品国产美女av久久久久小说| 一个人观看的视频www高清免费观看 | 中文亚洲av片在线观看爽| 国产av一区二区精品久久| 男女之事视频高清在线观看| 免费在线观看亚洲国产| 日本熟妇午夜| 老熟妇仑乱视频hdxx| 香蕉丝袜av| 91大片在线观看| 午夜福利在线在线| tocl精华| av免费在线观看网站| 国产一区二区在线av高清观看| 久久草成人影院| 欧美成人午夜精品| 国产黄色小视频在线观看| 高清在线国产一区| 色在线成人网| 99精品久久久久人妻精品| 亚洲va日本ⅴa欧美va伊人久久| 黄片小视频在线播放| 无人区码免费观看不卡| 中文亚洲av片在线观看爽| 黄色视频,在线免费观看| 日本在线视频免费播放| 天天躁夜夜躁狠狠躁躁| 国模一区二区三区四区视频 | 999久久久国产精品视频| av视频在线观看入口| 高潮久久久久久久久久久不卡| 丁香欧美五月| 无人区码免费观看不卡| 999久久久精品免费观看国产| 黄色视频,在线免费观看| 最近最新中文字幕大全电影3| 午夜福利高清视频| 青草久久国产| 亚洲专区国产一区二区| 欧美性猛交黑人性爽| 日本三级黄在线观看| 日韩精品青青久久久久久| a级毛片在线看网站| 看黄色毛片网站| 黄色女人牲交| 中出人妻视频一区二区| 日韩免费av在线播放| 久久亚洲真实| 好看av亚洲va欧美ⅴa在| 美女高潮喷水抽搐中文字幕| 日本三级黄在线观看| 久久久久亚洲av毛片大全| www.www免费av| 欧美乱色亚洲激情| 久久久久久国产a免费观看| 丰满人妻熟妇乱又伦精品不卡| 久久久久亚洲av毛片大全| 18禁裸乳无遮挡免费网站照片| 成人永久免费在线观看视频| 亚洲国产欧美人成| 精品国产乱子伦一区二区三区| 中文字幕人成人乱码亚洲影| 在线观看日韩欧美| 午夜两性在线视频| 中亚洲国语对白在线视频| 91九色精品人成在线观看| 国产成年人精品一区二区| 午夜福利成人在线免费观看| 日韩成人在线观看一区二区三区| 热99re8久久精品国产| 久久国产精品影院| 亚洲av熟女| 在线国产一区二区在线| 长腿黑丝高跟| 国产片内射在线| 亚洲欧美日韩无卡精品| 午夜激情福利司机影院| 亚洲精品在线美女| 国产成人影院久久av| 亚洲一区二区三区不卡视频| 久久久国产精品麻豆| 精品第一国产精品| 他把我摸到了高潮在线观看| 校园春色视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 精品无人区乱码1区二区| 久久久久久免费高清国产稀缺| 久久久久久久精品吃奶| 特级一级黄色大片| 国内毛片毛片毛片毛片毛片| 中文字幕av在线有码专区| 亚洲激情在线av| 狂野欧美激情性xxxx| ponron亚洲| 国产成人影院久久av| 岛国在线观看网站| 国产亚洲精品久久久久5区| 怎么达到女性高潮| 99精品欧美一区二区三区四区| 久久久久国产一级毛片高清牌| 成人av一区二区三区在线看| 1024视频免费在线观看| 欧美日韩精品网址| 亚洲一码二码三码区别大吗| 特级一级黄色大片| 国产精品av久久久久免费| 黄色a级毛片大全视频| 国产高清激情床上av| x7x7x7水蜜桃| 欧美又色又爽又黄视频| 黄色片一级片一级黄色片| 在线a可以看的网站| 国产高清有码在线观看视频 | 国产av不卡久久| www.熟女人妻精品国产| 男人舔女人下体高潮全视频| 婷婷丁香在线五月| 啦啦啦免费观看视频1| 97超级碰碰碰精品色视频在线观看| 色哟哟哟哟哟哟| 午夜福利在线观看吧| 在线十欧美十亚洲十日本专区| 国产精品爽爽va在线观看网站| 男女那种视频在线观看| 国产私拍福利视频在线观看| 国产高清videossex| av在线播放免费不卡| xxxwww97欧美| 亚洲精品在线观看二区| 免费搜索国产男女视频| 日本三级黄在线观看| 亚洲人成伊人成综合网2020| 成人一区二区视频在线观看| 亚洲av成人不卡在线观看播放网| 色综合站精品国产| 国产成人欧美在线观看| 亚洲av成人一区二区三| 在线视频色国产色| 高清毛片免费观看视频网站| 国产av麻豆久久久久久久| 日韩欧美国产一区二区入口| 51午夜福利影视在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产99白浆流出| 听说在线观看完整版免费高清| 亚洲avbb在线观看| 亚洲一区二区三区色噜噜| 久久精品影院6| 女人被狂操c到高潮| 国产亚洲精品久久久久久毛片| 中文字幕熟女人妻在线| 国产精品久久电影中文字幕| 日本成人三级电影网站| 亚洲五月天丁香| 少妇熟女aⅴ在线视频| 最近视频中文字幕2019在线8| 黄色毛片三级朝国网站| 久久香蕉国产精品| 波多野结衣巨乳人妻| 国产精品九九99| 国产精品电影一区二区三区| 女生性感内裤真人,穿戴方法视频| 亚洲av成人精品一区久久| 欧美午夜高清在线| 黑人欧美特级aaaaaa片| 亚洲男人天堂网一区| 欧美大码av| 两个人免费观看高清视频| 色综合欧美亚洲国产小说| 亚洲欧美日韩无卡精品| 特大巨黑吊av在线直播| 国产高清激情床上av| 精品久久久久久成人av| 日韩欧美一区二区三区在线观看| 亚洲第一电影网av| 久久天堂一区二区三区四区| xxx96com| 床上黄色一级片| 窝窝影院91人妻| 麻豆国产97在线/欧美 | 亚洲av成人不卡在线观看播放网| 97超级碰碰碰精品色视频在线观看| 亚洲av成人不卡在线观看播放网| 两人在一起打扑克的视频| 国产精品影院久久| e午夜精品久久久久久久| 亚洲精品在线观看二区| 欧美3d第一页| 欧美在线黄色| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久精品电影| 精品欧美一区二区三区在线| 成人av一区二区三区在线看| 国产午夜精品久久久久久| 亚洲美女黄片视频| 亚洲无线在线观看| 久久久国产成人免费| 亚洲自偷自拍图片 自拍| 宅男免费午夜| 两人在一起打扑克的视频| 国产99白浆流出| 欧美日韩亚洲综合一区二区三区_| 一卡2卡三卡四卡精品乱码亚洲| 午夜两性在线视频| 一个人免费在线观看电影 | 特大巨黑吊av在线直播| 久久精品人妻少妇| 在线观看午夜福利视频| 国产精品亚洲美女久久久| 国产精品免费一区二区三区在线| 成人国语在线视频| 我的老师免费观看完整版| 天天添夜夜摸| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 好男人在线观看高清免费视频| 精品福利观看| 天天添夜夜摸| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区在线av高清观看| 波多野结衣巨乳人妻| 成在线人永久免费视频| 久久精品亚洲精品国产色婷小说| 亚洲国产精品成人综合色| 亚洲专区国产一区二区| 好男人在线观看高清免费视频| 久久久久久九九精品二区国产 | av在线播放免费不卡| √禁漫天堂资源中文www| 午夜免费观看网址| 窝窝影院91人妻| 国产一区二区三区在线臀色熟女| 亚洲最大成人中文| 国产精品久久视频播放| 我的老师免费观看完整版| bbb黄色大片| 波多野结衣高清无吗| 亚洲成a人片在线一区二区| 丁香欧美五月| 一本精品99久久精品77| 一二三四社区在线视频社区8| 国产1区2区3区精品| 欧美性猛交╳xxx乱大交人| 不卡av一区二区三区| 啪啪无遮挡十八禁网站| 妹子高潮喷水视频| 国内精品久久久久久久电影| 美女免费视频网站| 黄色丝袜av网址大全| 99热只有精品国产| 国产成人精品久久二区二区免费| 精品免费久久久久久久清纯| 精品少妇一区二区三区视频日本电影| 香蕉国产在线看| 欧美+亚洲+日韩+国产| 国产视频内射| 日韩精品免费视频一区二区三区| 男女那种视频在线观看| 91成年电影在线观看| 国产伦一二天堂av在线观看| 欧美最黄视频在线播放免费| 麻豆一二三区av精品| 国产视频一区二区在线看| 叶爱在线成人免费视频播放| 国产精品电影一区二区三区| 久久精品aⅴ一区二区三区四区| 久久香蕉精品热| 亚洲精品美女久久久久99蜜臀| 一边摸一边抽搐一进一小说| 男女那种视频在线观看| 男人舔女人下体高潮全视频| 国模一区二区三区四区视频 | 久久久久国产精品人妻aⅴ院| 一夜夜www| 欧美日韩瑟瑟在线播放| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩亚洲综合一区二区三区_| 一级毛片精品| 宅男免费午夜| 成人18禁高潮啪啪吃奶动态图| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 午夜免费成人在线视频| 90打野战视频偷拍视频| 亚洲精品粉嫩美女一区| 高清毛片免费观看视频网站| 好看av亚洲va欧美ⅴa在| 88av欧美| 色av中文字幕| 国内少妇人妻偷人精品xxx网站 | 国产精品久久久人人做人人爽| 给我免费播放毛片高清在线观看| 久久热在线av| 亚洲自拍偷在线| 久久精品成人免费网站| 亚洲人成网站高清观看| 日本五十路高清| 国产欧美日韩一区二区精品| 搡老熟女国产l中国老女人| 欧美黑人欧美精品刺激| 国产视频内射| 99久久久亚洲精品蜜臀av| 国产亚洲精品久久久久5区| 中文在线观看免费www的网站 | 免费看美女性在线毛片视频| 色av中文字幕| av天堂在线播放| 每晚都被弄得嗷嗷叫到高潮| 久久天堂一区二区三区四区| 亚洲欧美日韩高清专用| 不卡av一区二区三区| 国产野战对白在线观看| 在线观看免费日韩欧美大片| 一本大道久久a久久精品| 国产欧美日韩精品亚洲av| 国产成人啪精品午夜网站| 人人妻人人澡欧美一区二区| 亚洲精华国产精华精| 超碰成人久久| 久久久精品欧美日韩精品| 免费观看精品视频网站| 在线永久观看黄色视频| xxxwww97欧美| 老熟妇仑乱视频hdxx| 国产成人影院久久av| 老熟妇乱子伦视频在线观看| 久久中文字幕人妻熟女| 天堂av国产一区二区熟女人妻 | 免费av毛片视频| 欧美3d第一页| 999久久久精品免费观看国产| 亚洲精华国产精华精| 又黄又爽又免费观看的视频| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 日韩欧美一区二区三区在线观看| 国产区一区二久久| 嫩草影视91久久| 国产欧美日韩一区二区精品| 国产精品久久久久久久电影 | 少妇粗大呻吟视频| 亚洲五月婷婷丁香| 国产精品98久久久久久宅男小说| 免费观看人在逋| 久久精品91无色码中文字幕| 欧美乱色亚洲激情| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| 在线免费观看的www视频| 国内久久婷婷六月综合欲色啪| 18禁黄网站禁片免费观看直播| 狂野欧美白嫩少妇大欣赏| 999精品在线视频| 午夜免费激情av| 又紧又爽又黄一区二区| 亚洲av熟女| 亚洲欧美日韩东京热| 在线观看美女被高潮喷水网站 | 99精品在免费线老司机午夜| 欧美午夜高清在线| 中文在线观看免费www的网站 | 正在播放国产对白刺激| 草草在线视频免费看| 色av中文字幕| 精华霜和精华液先用哪个| 国产精品亚洲av一区麻豆| 久久国产乱子伦精品免费另类| 精品无人区乱码1区二区| 天天躁夜夜躁狠狠躁躁| 制服诱惑二区| 亚洲中文字幕一区二区三区有码在线看 | 欧美精品啪啪一区二区三区| 免费看美女性在线毛片视频| 欧美日韩国产亚洲二区| 嫩草影院精品99| 真人做人爱边吃奶动态| 九色成人免费人妻av| 日韩大尺度精品在线看网址| 亚洲av成人av| 免费在线观看日本一区| 亚洲国产欧洲综合997久久,| 久久久精品大字幕| АⅤ资源中文在线天堂| 91老司机精品| 五月玫瑰六月丁香| 久久久久精品国产欧美久久久| 久久婷婷人人爽人人干人人爱| 亚洲欧美精品综合久久99| 成人三级黄色视频| 成年免费大片在线观看| 久久欧美精品欧美久久欧美| 国产亚洲欧美98| av中文乱码字幕在线| 国产成人av教育| 色哟哟哟哟哟哟| 高清毛片免费观看视频网站| 岛国在线观看网站| 91av网站免费观看| 国产探花在线观看一区二区| 国产日本99.免费观看| 亚洲中文日韩欧美视频| 日韩成人在线观看一区二区三区| 叶爱在线成人免费视频播放| 一本大道久久a久久精品| 18禁美女被吸乳视频| 不卡av一区二区三区| 欧美成人性av电影在线观看| 久久久久久久久中文| 精品乱码久久久久久99久播| 欧美久久黑人一区二区| 十八禁网站免费在线| 久久久精品大字幕| 一二三四社区在线视频社区8| 免费在线观看日本一区| 91成年电影在线观看| 亚洲av中文字字幕乱码综合| 九色成人免费人妻av| 亚洲av片天天在线观看| 亚洲精品美女久久av网站| 免费一级毛片在线播放高清视频| 国产午夜精品久久久久久| 黄色毛片三级朝国网站| 国产av不卡久久| 男女那种视频在线观看| 久久久久久久久久黄片| 夜夜躁狠狠躁天天躁| 露出奶头的视频| 日韩中文字幕欧美一区二区| 久久久久国产一级毛片高清牌| 国产黄a三级三级三级人|