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

    熱帶大西洋對(duì)全球變暖的響應(yīng)

    2020-06-28 00:21:26董文靜劉福凱羅義勇
    海洋科學(xué) 2020年6期
    關(guān)鍵詞:正位赤道大西洋

    董文靜, 劉福凱, 羅義勇

    熱帶大西洋對(duì)全球變暖的響應(yīng)

    董文靜, 劉福凱, 羅義勇

    (中國(guó)海洋大學(xué) 海洋與大氣學(xué)院, 山東 青島 266100)

    與太平洋和印度洋不同, 全球變暖下熱帶大西洋變化的研究較少。本文使用地球系統(tǒng)模型CESM(Community Earth System Model), 發(fā)現(xiàn)全球變暖后熱帶大西洋在秋季的升溫類(lèi)似大西洋尼諾(Atlantic Ni?o)的正位相, 即大西洋西部增暖幅度小于東部; 在夏季類(lèi)似大西洋尼諾的負(fù)位相, 即大西洋西部增暖幅度大于東部。利用覆蓋(overriding)技術(shù), 分離了風(fēng)應(yīng)力、風(fēng)速和CO2的直接熱效應(yīng)對(duì)海洋升溫的作用, 探討了大西洋尼諾本身和全球變暖作用下類(lèi)似大西洋尼諾正位相(下文簡(jiǎn)稱(chēng)“類(lèi)大西洋尼諾升溫”)的形成機(jī)制。結(jié)果表明, 這兩種情況下的形成機(jī)制基本相同, 風(fēng)應(yīng)力的變化是導(dǎo)致大西洋東部暖異常的主要機(jī)制。但兩者之間也存在區(qū)別: 1) 全球變暖下海表溫度的季節(jié)變化振幅減小, 而大西洋尼諾時(shí)變化不大; 2) 全球變暖下西風(fēng)異常主要集中在大西洋東部, 而大西洋尼諾時(shí)主要集中在大西洋中部; 3) 除風(fēng)應(yīng)力外, CO2的熱效應(yīng)對(duì)類(lèi)尼諾升溫的變化也有一定影響。

    熱帶大西洋; 全球變暖; 大西洋尼諾; 類(lèi)大西洋尼諾升溫

    大西洋尼諾事件, 也稱(chēng)為大西洋赤道變化模態(tài)或赤道大西洋緯向變化模態(tài), 是赤道大西洋年際變化的主要模態(tài)。它能夠引起赤道大西洋海表溫度的冷暖交替變化, 對(duì)區(qū)域降水和生態(tài)系統(tǒng)也有深遠(yuǎn)影響[1-2]。將大西洋尼諾現(xiàn)象與太平洋厄爾尼諾現(xiàn)象對(duì)比, 峰值期時(shí)它們有許多相似之處: 最高值發(fā)生在冷舌區(qū)、沿赤道的緯向溫度梯度減小、沿赤道的溫躍層傾斜減弱[3-4]。大西洋尼諾和太平洋厄爾尼諾兩個(gè)事件的形成均與海氣Bjerknes正反饋密切相關(guān)[5]。關(guān)于大西洋尼諾正位相向負(fù)位相的過(guò)渡, Foltz等[6]提出了一種與赤道羅斯比波和西邊界反射相關(guān)的延遲負(fù)反饋, 類(lèi)似于解釋厄爾尼諾-南方濤動(dòng)(El Ni?o - southern oscillation , ENSO)現(xiàn)象循環(huán)的延遲振蕩反饋[7-8]。然而, 與太平洋厄爾尼諾相比, 大西洋尼諾的強(qiáng)度更小, 峰值的季節(jié)不同(大西洋尼諾的峰值期是北半球夏季), 峰值期的持續(xù)時(shí)間更短(3個(gè)月), 可預(yù)測(cè)性較低[3, 9-10]。而且, 大西洋的Bjerknes正反饋強(qiáng)度也比太平洋弱, 這說(shuō)明大西洋尼諾的形成機(jī)制與厄爾尼諾存在差異[1, 3]。

    全球變暖下, 熱帶太平洋和熱帶印度洋的響應(yīng)已經(jīng)有大量的研究。通常在模型預(yù)測(cè)中, 溫室氣體排放會(huì)導(dǎo)致熱帶太平洋的類(lèi)厄爾尼諾升溫(El Ni?o-like)和熱帶印度洋的類(lèi)印度洋偶極子正位相(pIOD-like)[11-13]。而類(lèi)厄爾尼諾(El Ni?o-like)響應(yīng)與厄爾尼諾本身有較大區(qū)別。例如, 前者情景下上層海洋層結(jié)增強(qiáng)、溫躍層變淺, 而后者層結(jié)減弱、溫躍層加深[18]。此外, 研究已經(jīng)發(fā)現(xiàn), 熱帶太平洋氣候變化由動(dòng)力學(xué)機(jī)制驅(qū)動(dòng)[14], 厄爾尼諾發(fā)展的主要原因是東風(fēng)減弱。而全球變暖下類(lèi)厄爾尼諾發(fā)展的主要原因則是東風(fēng)減弱和CO2熱效應(yīng)減小。前者對(duì)季節(jié)的變化貢獻(xiàn)很大, 后者則對(duì)熱帶太平洋年平均溫度升高的貢獻(xiàn)更大[13]。在熱帶印度洋, 全球變暖引起的類(lèi)印度洋偶極子正位相(pIOD-like)與印度洋偶極子正位相(pIOD)具有幾乎相同的形成過(guò)程和相關(guān)季節(jié)特性, 兩種情況下產(chǎn)生東部熱帶冷卻異常的主要機(jī)制都是風(fēng)應(yīng)力(WStr)反饋機(jī)制[13]。

    與熱帶太平洋和印度洋相比, 目前為止關(guān)于熱帶大西洋對(duì)全球變暖響應(yīng)的研究很少, 最近Tokinaga等[16]通過(guò)分析1950—2009年的觀測(cè)數(shù)據(jù), 發(fā)現(xiàn)全球變暖后熱帶大西洋在北半球秋季有一個(gè)類(lèi)大西洋尼諾正位相的響應(yīng), 特點(diǎn)為赤道東大西洋(Equatorial Eastern Atlantic, EEA)海表面溫度(sea surface temperature)SST變暖增強(qiáng), 貿(mào)易風(fēng)減弱, 溫躍層沿赤道變平。在此基礎(chǔ)上, 本文擬利用地球系統(tǒng)模型(Community Earth System Model), 分析全球變暖導(dǎo)致的類(lèi)大西洋尼諾升溫現(xiàn)象的變化特征, 并將這種變化特征與大西洋尼諾現(xiàn)象本身對(duì)比, 探討兩者形成機(jī)制的異同。通過(guò)分析, 我們發(fā)現(xiàn)在全球變暖下熱帶大西洋呈現(xiàn)類(lèi)大西洋尼諾升溫, 這證實(shí)了Tokinaga等[16]的發(fā)現(xiàn)。在季節(jié)變化尺度上, 全球變暖引起的類(lèi)大西洋尼諾升溫在北半球夏季(5—7月)較弱,而在北半球秋季(8—10月)最明顯。此外, 我們還進(jìn)行了一系列覆蓋實(shí)驗(yàn), 以分離風(fēng)應(yīng)力、風(fēng)速和CO2的直接熱效應(yīng)對(duì)類(lèi)大西洋尼諾升溫及大西洋尼諾的作用, 結(jié)果表明這兩種情況下, 風(fēng)應(yīng)力反饋對(duì)赤道東大西洋溫度異常的形成均起主要作用。本文還進(jìn)一步分析了熱收支平衡, 以進(jìn)一步驗(yàn)證類(lèi)尼諾升溫和大西洋尼諾的形成機(jī)制。

    本文的其余部分結(jié)構(gòu)如下: 下一節(jié)將介紹模型和數(shù)值實(shí)驗(yàn)。在第3節(jié)介紹合成大西洋尼諾事件和計(jì)算熱收支平衡的方法。第4節(jié)中, 比較了類(lèi)大西洋尼諾升溫和大西洋尼諾之間的海洋與氣象的變化特征, 其中包括平均態(tài)變化特征、季節(jié)變化特征、次表層變化特征和北向熱輸運(yùn)變化特征。在第5節(jié)中進(jìn)行了熱收支平衡分析。第6節(jié)則對(duì)研究結(jié)果進(jìn)行了總結(jié)和討論。

    1 數(shù)據(jù)與方法

    1.1 CESM模式介紹

    本文使用CESM模式, 具體版本是1.1.2(http:// www.cesm.ucar.edu/models/cesm1.1/index.html)。它由大氣、海洋、陸面等分量模式組成, 并由耦合器進(jìn)行各個(gè)分量模式之間的數(shù)據(jù)信息交換和模式運(yùn)行。其中大氣模式為CAM5, 陸面模式為CLM4, 海洋模式為POP2。CAM5和CLM4的空間分辨率為1.9°×1.9°, 垂向分為30層。POP2的空間分辨率為1°, 在赤道附近南北向加密為約0.3°。垂直方向有60個(gè)非均勻分層, 上層海洋分辨率高, 從上到下分辨率依次遞減, 表層10 m, 底層接近250 m。

    從國(guó)家大氣研究中心 (National Center for Atmospheric Research, NCAR)控制實(shí)驗(yàn)(1861—2005年)的結(jié)束時(shí)間起始, 首先, 我們使用CESM模型在RCP8.5場(chǎng)景下運(yùn)行了94 a(2006—2099年), 得到海洋和大氣變量日平均的輸出, 此實(shí)驗(yàn)記為“CPL85”(表1)。需要注意的是, RCP8.5是IPCC第五次報(bào)告中提出的一種碳排放場(chǎng)景假設(shè), 即到2100年時(shí), 空氣中的二氧化碳濃度達(dá)到工業(yè)革命前的濃度高3~4倍。在覆蓋實(shí)驗(yàn)中, 我們首先用CPL85實(shí)驗(yàn)中2006年的大氣強(qiáng)迫數(shù)據(jù)(包括風(fēng)、氣溫、氣壓、比濕、降水率、空氣密度、凈短波輻射和向下的長(zhǎng)波輻射)反復(fù)驅(qū)動(dòng)POP2, 從而得到94 a的“CTRL”實(shí)驗(yàn)(表1)。接下來(lái)我們使用CPL85實(shí)驗(yàn)94 a的大氣強(qiáng)迫場(chǎng)驅(qū)動(dòng)POP2, 得到“FULL”(表1)實(shí)驗(yàn)。另外, 為了隔離風(fēng)應(yīng)力(風(fēng)速)變化的影響, 我們用2006年的風(fēng)應(yīng)力(風(fēng)速)反復(fù)驅(qū)動(dòng)海洋, 所有其他場(chǎng)與FULL一樣使用94 a數(shù)據(jù)的循環(huán), 從而得到STRS(“SPED”)實(shí)驗(yàn)。需要注意的是, STRS實(shí)驗(yàn)中的風(fēng)應(yīng)力和SPED中的風(fēng)速對(duì)SST的作用機(jī)制有很大不同, 風(fēng)應(yīng)力的貢獻(xiàn)主要通過(guò)風(fēng)應(yīng)力-垂向流速-SST機(jī)制(即Bjerknes正反饋機(jī)制)實(shí)現(xiàn), 而風(fēng)速對(duì)海表面溫度的貢獻(xiàn)則主要是通過(guò)風(fēng)速-蒸發(fā)-SST機(jī)制(即WES機(jī)制)實(shí)現(xiàn)。另外, 我們同時(shí)固定2006年的風(fēng)應(yīng)力和風(fēng)速?gòu)亩玫絎IND實(shí)驗(yàn), 以驗(yàn)證STRS實(shí)驗(yàn)和SPED實(shí)驗(yàn)的可加性, 結(jié)果也驗(yàn)證了兩個(gè)實(shí)驗(yàn)有較好的線(xiàn)性可加性。

    表1 CESM1.1和POP2實(shí)驗(yàn)介紹

    由此, FULL-STRS實(shí)驗(yàn)可以計(jì)算出風(fēng)應(yīng)力對(duì)海洋變化的影響(WStr反饋), FULL-SPED實(shí)驗(yàn)可以計(jì)算出風(fēng)速的貢獻(xiàn)(WES反饋), WIND-CTRL實(shí)驗(yàn)可以計(jì)算出CO2的直接熱效應(yīng)(CO2)。需要強(qiáng)調(diào)的是, 這里WES反饋僅考慮了風(fēng)速變化對(duì)海洋的直接熱效應(yīng),不包括完全耦合模式中WES反饋對(duì)大氣過(guò)程的間接影響[11]。FULL-CTRL實(shí)驗(yàn)?zāi)M了耦合CESM1.1模型中的全部效應(yīng), 包含了上面的所有反饋。我們后續(xù)的研究結(jié)果都基于月平均場(chǎng), 大西洋尼諾下各物理量的異常值來(lái)源于各實(shí)驗(yàn)的數(shù)值場(chǎng)減去本身的氣候態(tài)平均, 而全球變暖下的類(lèi)尼諾升溫則是計(jì)算了2006—2099年的線(xiàn)性趨勢(shì)場(chǎng)。

    這一實(shí)驗(yàn)方法已經(jīng)被成功運(yùn)用于太平洋類(lèi)厄爾尼諾升溫等現(xiàn)象的研究中[13, 18], 充分說(shuō)明了該技術(shù)的可靠性。從物理性質(zhì)上講, 前人運(yùn)用該實(shí)驗(yàn)得出的結(jié)論和理論結(jié)果較吻合, 如: 風(fēng)應(yīng)力效應(yīng)造成的海溫變化主要集中在次表層, 而CO2的直接效應(yīng)所造成的溫度變化大多從表層向下遞減。前人得到的結(jié)論也發(fā)表在一系列期刊上, 得到了廣泛認(rèn)可。

    1.2 分析方法

    1.2.1 熱收支平衡方程

    熱收支平衡分析是分析海洋混合層響應(yīng)機(jī)制的重要方法, 本文選擇55 m作為固定底部計(jì)算赤道東大西洋的熱收支平衡。

    熱收支平衡方程為:

    T=Q+Q+Q+Q+Q,(1)

    式中,T為混合層溫度(mixing layer temprature, MLT)隨時(shí)間的變化;Q=(0-Q)/0p其中0為海表面熱通量(surface heat flux),Q為短波輻射在混合層底的量值,0=1 025 kg/m3和p=2 850 J/(kg?K)分別是海水密度和海水絕熱系數(shù);QQ分別為緯向(zonal advection)和經(jīng)向(meridional advection)的溫度平流;Q是垂直對(duì)流(vertical advection)過(guò)程;r是方程的余項(xiàng), 代表了次網(wǎng)格尺度和時(shí)間尺度為月以下的運(yùn)動(dòng)過(guò)程, 以及在方程計(jì)算過(guò)程中可能出現(xiàn)的誤差。雖然由于缺乏更高時(shí)空分辨率的數(shù)據(jù), 次網(wǎng)格尺度和月以下尺度的海洋動(dòng)力過(guò)程不能被精確計(jì)算, 但是混合層熱收支方程中的余項(xiàng)主要是由于垂直擴(kuò)散作用導(dǎo)致, 因而可以被看作垂直擴(kuò)散(vertical diffusion)項(xiàng)。在方程中, 右邊四項(xiàng)的和代表三個(gè)維度上海洋的運(yùn)動(dòng)過(guò)程和擴(kuò)散作用導(dǎo)致的海洋熱輸運(yùn)(ocean heat transport)。

    混合層熱收支方程的一個(gè)優(yōu)點(diǎn)是可以定量地比較導(dǎo)致溫度變化的各種過(guò)程的貢獻(xiàn), 從而確定各項(xiàng)的相對(duì)重要性。在熱收支方程中, 正的熱收支項(xiàng)代表海洋從大氣吸收熱量, 而負(fù)的熱收支項(xiàng)代表海洋向大氣釋放熱量。

    本文使用公式(1), 計(jì)算了MLT的年際變化及其相關(guān)過(guò)程。

    1.2.2 合成大西洋尼諾(Atlantic Ni?o)

    與太平洋厄爾尼諾指數(shù)區(qū)域不同, 大西洋沒(méi)有普遍共識(shí)的類(lèi)似指數(shù)區(qū)域[2]。在本研究中, MLT變化最大的區(qū)域位于赤道東大西洋區(qū)域(對(duì)應(yīng)Zebiak[9]在1993年定義的ATL3區(qū)域)。定義MLT變化最大的區(qū)域ATL0(3°N—3°S, 15°W—10°E)為本文的研究區(qū)域。為了篩選大西洋尼諾事件, 我們首先去除CPL85數(shù)據(jù)時(shí)間序列(2006—2099年)的趨勢(shì)項(xiàng), 得到赤道東大西洋區(qū)域的MLT異常(圖1)。本文將兩倍標(biāo)準(zhǔn)差(0.98℃)作為選擇標(biāo)準(zhǔn), 大于0.98℃的時(shí)間定義為一個(gè)大西洋尼諾事件。通過(guò)這種方法, 我們從94 a的模擬時(shí)間序列中識(shí)別出13個(gè)大西洋尼諾事件。圖2展示了這13個(gè)事件合成的季節(jié)變化: 溫度異常在北半球夏季(6月)達(dá)到峰值(圖2中虛線(xiàn)黑線(xiàn)), 這個(gè)結(jié)果與前人[1-2]的研究相近。為了便于討論, 我們將大西洋尼諾的變化分為三個(gè)不同的階段: (a)1—4月的形成期; (b)5—7月的峰值期; (c)8—12月的衰減期。在下一節(jié)中, 可以看到合成大西洋尼諾事件能夠很好地表現(xiàn)出大西洋尼諾現(xiàn)象的主要特征。

    圖1 赤道東大西洋混合層溫度(mixing layer temprature, MLT)異常的時(shí)間序列

    注: 藍(lán)色虛線(xiàn)表示此時(shí)間序列的兩倍標(biāo)準(zhǔn)差(約0.98 ℃)

    由于FULL-CTRL實(shí)驗(yàn)結(jié)果準(zhǔn)確地再現(xiàn)了CPL85實(shí)驗(yàn)結(jié)果的年際變化, 我們使用CPL85實(shí)驗(yàn)中確定的13個(gè)事件的發(fā)生時(shí)間, 用于覆蓋實(shí)驗(yàn)中大西洋尼諾事件的合成(圖2黑色實(shí)線(xiàn))。發(fā)現(xiàn)CPL85實(shí)驗(yàn)合成大西洋尼諾事件的變化特征與FULL-CTRL實(shí)驗(yàn)(圖2黑色虛線(xiàn))基本相似, 但振幅略大, 這可能是由于ocean-only實(shí)驗(yàn)中缺乏高頻海氣通量的原因。如圖所示, WStr效應(yīng)在大西洋尼諾變化中起主導(dǎo)作用(圖2藍(lán)色線(xiàn)), 而WES效應(yīng)的貢獻(xiàn)可以忽略不計(jì)(圖2紅色線(xiàn))。有趣的是, 我們發(fā)現(xiàn)在沒(méi)有風(fēng)應(yīng)力和風(fēng)速影響時(shí), CO2的直接熱效應(yīng)能夠在大西洋尼諾期間引起暖異常(圖2綠色線(xiàn))。

    圖2 合成大西洋尼諾事件在各實(shí)驗(yàn)中的季節(jié)變化

    本文后面的分析主要基于CPL85模擬實(shí)驗(yàn), 以及通過(guò)覆蓋實(shí)驗(yàn)進(jìn)一步提取各反饋過(guò)程在大西洋尼諾和類(lèi)尼諾升溫中的作用。由于CPL85實(shí)驗(yàn)結(jié)果和FULL-CTRL實(shí)驗(yàn)幾乎相同, 所以下文中只展示FULL-CTRL實(shí)驗(yàn)的結(jié)果。此外, 由于FULL-SPED (WES反饋)對(duì)熱帶大西洋變化的影響微不足道, 因此在本文的剩余部分中, 我們并沒(méi)有給出WES反饋的結(jié)果。為了便于與大西洋尼諾MLT特征相比較, 本文去掉了全球變暖下熱帶大西洋海域(20°S—20°N, 3°S—3°N的平均)的平均升溫。在這種情況下, 如果MLT升溫小于海域變暖平均, 即表現(xiàn)為圖3b、圖4b和圖5b中的冷信號(hào)。

    2 結(jié)果與分析

    2.1 熱帶大西洋的海洋和大氣變化

    2.1.1 季節(jié)變化

    圖3分別展示了大西洋尼諾和類(lèi)尼諾升溫情況下沿赤道MLT和環(huán)流的季節(jié)變化。與之前的研究[1, 3]一致, 在大西洋尼諾正位相期間, 赤道大西洋中、西部的西風(fēng)異常導(dǎo)致赤道東大西洋MLT異常, 并且延遲約1個(gè)月(圖3a和圖3c)。赤道東大西洋的MLT異常在峰值期最高達(dá)到約2℃。此外, 赤道大西洋表層流減弱(圖3e), 上升流減弱(圖3g), 上層海洋層結(jié)減弱(圖3i)。赤道大西洋西風(fēng)異常的逐漸減弱, 也導(dǎo)致了大西洋尼諾現(xiàn)象的衰減。赤道東大西洋的大西洋尼諾異常最為突出, 覆蓋實(shí)驗(yàn)的結(jié)果顯示風(fēng)應(yīng)力在大西洋尼諾事件中起到至關(guān)重要的作用(圖4左側(cè))。也就是說(shuō), WStr效應(yīng)對(duì)赤道東大西洋MLT異常變暖以及季節(jié)變化起主導(dǎo)作用, 而CO2的熱效應(yīng)對(duì)其異常變暖以及季節(jié)變化幾乎沒(méi)有貢獻(xiàn)。

    在熱帶大西洋, 類(lèi)尼諾升溫和大西洋尼諾引起的MLT異常季節(jié)變化及風(fēng)場(chǎng)異常的季節(jié)變化均有明顯差異。北半球夏季, 類(lèi)大西洋尼諾升溫時(shí), 東風(fēng)正異常最大, 而在大西洋尼諾正位相期間, 東風(fēng)負(fù)異常最大(圖3d)。類(lèi)似地, 熱帶大西洋對(duì)類(lèi)大西洋尼諾升溫和大西洋尼諾正位相的反應(yīng)基本是相反的: 赤道東大西洋MLT暖異常達(dá)到最小(圖3b), 緯向流速增大(圖3f), 上升流加強(qiáng)(圖3h), 上層海洋層結(jié)增大(圖3j)。在北半球秋季, 緯向風(fēng)應(yīng)力異常方向快速轉(zhuǎn)換, 而在大西洋尼諾事件發(fā)生時(shí), 風(fēng)異常隨時(shí)間逐漸減弱。因此, 在北半球夏季出現(xiàn)類(lèi)大西洋尼諾降溫之后, 在秋季又出現(xiàn)了類(lèi)大西洋尼諾升溫。覆蓋實(shí)驗(yàn)進(jìn)一步表明, 類(lèi)尼諾升溫和大西洋尼諾引起的MLT的季節(jié)變化均由WStr效應(yīng)主導(dǎo)(圖4a和圖4b); 與大西洋尼諾不同(圖4c), 類(lèi)尼諾升溫引發(fā)的層結(jié)變化是WStr效應(yīng)和CO2的熱效應(yīng)正變化的疊加(圖4d)。

    注: 類(lèi)尼諾升溫情景下表現(xiàn)的是94 a的趨勢(shì)值; 等值線(xiàn)分別是各量CPL85的原始值

    圖4 大西洋尼諾(左)和全球變暖類(lèi)尼諾升溫(右)情景下海洋特征在赤道東大西洋的變化

    注: a、b: MLT; c、d: 55 m深的層結(jié); e、f: 表面到55 m深垂直流速平均

    有趣的是, 我們的研究結(jié)果表明, 在類(lèi)尼諾升溫下, MLT的季節(jié)循環(huán)位相提前, 幅度減小(圖3b)。從氣候?qū)W的角度看, 赤道大西洋MLT呈現(xiàn)出明顯的季節(jié)周期。季節(jié)周期的特點(diǎn)是MLT在4—6月(7—9月)變暖(冷), 歷經(jīng)1~2個(gè)月的相位提前和振幅下降, 暖異常在5—7月(8—10月)達(dá)到最小(最大)。季節(jié)循環(huán)信號(hào)減弱這一現(xiàn)象在其他海洋特征中也有體現(xiàn)(圖3d、圖3f、圖3h和圖3j)。這一結(jié)果與Tokinaga等[16]的觀察研究一致, 他們也發(fā)現(xiàn)在過(guò)去60 a里赤道大西洋的季節(jié)循環(huán)減弱。

    2.1.2 大西洋尼諾正相位和類(lèi)尼諾升溫的空間特征

    我們之后重點(diǎn)研究了大西洋尼諾正位相和類(lèi)尼諾升溫峰值期的空間特征, 以及二者之間的異同。圖5左側(cè)各圖顯示的是大西洋尼諾在其峰值期的各種特征: 圖5a中的大西洋尼諾正位相與太平洋厄爾尼諾正位相非常相似, 但僅限于赤道區(qū)域; 赤道東大西洋暖異常最大, 其最大值區(qū)西部的信風(fēng)減弱(圖5a和圖5b); 信風(fēng)減弱導(dǎo)致表面西向流減弱, 赤道上升流減弱(圖5e和圖5g), 溫躍層變平(圖6c), 以及進(jìn)一步的升溫, 即Bjerknes正反饋機(jī)制。此外, 海洋的變化特征還包括整個(gè)赤道的凈熱損失(圖5k), 東(西)部層結(jié)(圖5i)的顯著減弱(增強(qiáng))等。

    圖5 大西洋尼諾(左)和全球變暖類(lèi)尼諾升溫(右)情景下在其峰值期各海洋特征的空間分布

    注: 類(lèi)尼諾升溫情境下表現(xiàn)的是94 a的趨勢(shì)值; 等值線(xiàn)分別是各量CPL85的原始值

    圖5右側(cè)各圖展示了北半球秋季類(lèi)尼諾升溫情況下各海洋特征量趨勢(shì)的空間分布。在赤道大西洋, MLT呈現(xiàn)出明顯的類(lèi)似大西洋尼諾正位相的模態(tài): 中部和東部赤道變暖加劇, 西部變暖減弱(圖5b)。對(duì)比圖5的左右兩側(cè), 盡管類(lèi)尼諾升溫與大西洋尼諾正位相之間有很多相似之處, 但差異也很明顯: 大西洋尼諾比類(lèi)尼諾升溫關(guān)于赤道的對(duì)稱(chēng)性更大; 后者在熱帶大西洋東南部暖異常最小, 這可能與信風(fēng)變化(見(jiàn)圖5d)的非對(duì)稱(chēng)性有關(guān)[11, 13, 17]。類(lèi)尼諾升溫和大西洋尼諾的另一個(gè)顯著區(qū)別在于赤道東大西洋上層海洋的層結(jié)特征不同: 在大西洋尼諾期間, 西部層結(jié)增強(qiáng), 東部層結(jié)減少; 但在類(lèi)尼諾升溫情況下, 這種特征不存在(圖5j和圖5i比較)。造成這種差異的原因在于風(fēng)應(yīng)力異常不同: 類(lèi)尼諾升溫情況下赤道東大西洋西風(fēng)異常而大西洋尼諾期間赤道東大西洋則為東風(fēng)異常(圖5d和圖5c, 圖3c和圖3d)。風(fēng)應(yīng)力的緯向位移則使赤道東大西洋層結(jié)加強(qiáng), 赤道西大西洋層結(jié)減弱。

    2.1.3 類(lèi)尼諾升溫和大西洋尼諾的次表層變化

    圖6所示為大西洋尼諾正位相和類(lèi)尼諾升溫兩種情況下沿赤道的溫度和溫躍層深度的次表層變化。模式結(jié)果再現(xiàn)了大西洋尼諾正位相的溫躍層變化的主要觀測(cè)特征: 夏季西部加深, 東部變淺(圖6c)。符合Foltz等[6]提出的負(fù)反饋機(jī)制。

    在類(lèi)大西洋尼諾負(fù)位相時(shí)期, 因?yàn)轭?lèi)尼諾升溫的風(fēng)應(yīng)力異常與大西洋尼諾相反, 所以類(lèi)尼諾升溫的次表層溫度異常也與大西洋尼諾完全相反(圖6d與圖6c)。赤道大西洋溫躍層的溫度異常在東西部正負(fù)不同, 在類(lèi)大西洋尼諾正位相時(shí)期, 東西部全部變成正值(圖6f)。

    圖6 大西洋尼諾(左)和全球變暖類(lèi)尼諾升溫(右)情境下沿赤道(2.5°S—2.5°N的平均)溫度異常的季節(jié)變化

    注: 綠線(xiàn)和藍(lán)線(xiàn): 混合層底深度(綠色是2006—2025年的平均, 左側(cè)藍(lán)色是合成大西洋尼諾的平均, 右側(cè)藍(lán)色則是2080—2099年的平均); 等值線(xiàn): CPL85的94 a平均值; 溫躍層深度是取的溫度梯度的最大值

    在類(lèi)尼諾升溫和大西洋尼諾現(xiàn)象下, 次表層溫度異常都要比表面溫度異常大得多且東部溫度異常最大值均在200 m左右。說(shuō)明在這兩種情況下海洋動(dòng)力都很重要。

    2.2 混合層熱收支分析

    為了進(jìn)一步了解大西洋變暖機(jī)制, 我們利用公式(1)分析了赤道大西洋熱量收支平衡, 并分別計(jì)算了它們對(duì)大西洋尼諾和類(lèi)尼諾升溫的響應(yīng)。圖7展示了大西洋尼諾和類(lèi)尼諾升溫下MLT趨勢(shì)的季節(jié)變化。根據(jù)結(jié)果我們發(fā)現(xiàn), 兩種情況下垂向?qū)α?升溫)和擴(kuò)散(降溫)對(duì)赤道東大西洋溫度異常的貢獻(xiàn)都非常顯著, 而經(jīng)緯向水平平流和凈熱通量等其他項(xiàng)的貢獻(xiàn)幾乎可以忽略不計(jì)。

    2.2.1 垂直對(duì)流

    在氣候態(tài)平衡中, 垂直對(duì)流是赤道東大西洋區(qū)域的主要冷卻項(xiàng)(圖7e等值線(xiàn)), 是由風(fēng)應(yīng)力引起的強(qiáng)上升流(圖3g等值線(xiàn))和垂向溫度梯度增大(圖3i等值線(xiàn))造成的。根據(jù)圖8, 垂直對(duì)流冷異常變?nèi)鯇?duì)赤道東大西洋的溫度暖異常貢獻(xiàn)最大, 從而促進(jìn)大西洋尼諾現(xiàn)象的發(fā)展和衰減, 這與Richter[1]的實(shí)驗(yàn)結(jié)果一致。垂向熱輸運(yùn)異常是赤道東大西洋垂向速度減小(圖4g)和層結(jié)變?nèi)?圖3i)共同的結(jié)果。而垂向速度和層結(jié)的變化都與風(fēng)應(yīng)力有關(guān), 前者是由于赤道東風(fēng)減弱, 艾克曼抽吸減弱, 上升流速減弱; 后者則是因?yàn)槌嗟罇|風(fēng)減弱, 赤道大西洋西部水團(tuán)堆積減少, 溫躍層坡度變緩, 赤道大西洋西部的層結(jié)變?nèi)?。因此WStr效應(yīng)的主導(dǎo)作用可以被證實(shí)。并且覆蓋實(shí)驗(yàn)也可以驗(yàn)證WStr效應(yīng)的主導(dǎo)作用(圖8)。

    圖7 大西洋尼諾(左)和全球變暖類(lèi)尼諾升溫(右)情景下熱収支平衡方程各項(xiàng)的季節(jié)分布

    注: 等值線(xiàn)分別是各量CPL85的原始值

    類(lèi)尼諾升溫情況下的垂直對(duì)流導(dǎo)致北半球夏季溫度負(fù)異常, 秋季溫度正異常, 即決定了類(lèi)尼諾升溫下MLT的季節(jié)變化(圖7d和7b比較)。覆蓋實(shí)驗(yàn)進(jìn)一步表明, 與大西洋尼諾現(xiàn)象不同, CO2的熱效應(yīng)對(duì)赤道東大西洋的類(lèi)尼諾升溫的季節(jié)變化也有一定貢獻(xiàn)。然而, 就類(lèi)尼諾升溫相位變化而言, WStr效應(yīng)的貢獻(xiàn)要大很多。

    2.2.2 擴(kuò)散

    氣候態(tài)平衡中, 擴(kuò)散項(xiàng)是赤道東大西洋主要的加熱項(xiàng)。在類(lèi)尼諾升溫和大西洋尼諾情況下, 在擴(kuò)散和垂直熱輸運(yùn)異常之間都存在補(bǔ)償關(guān)系(圖7c和圖7e, 圖7d和圖7f)。例如, 在大西洋尼諾正位相和類(lèi)尼諾升溫的峰值期, 赤道東大西洋的特征均為擴(kuò)散項(xiàng)減小, 加熱減弱, 垂向熱輸運(yùn)項(xiàng)增大, 冷卻加強(qiáng)。覆蓋實(shí)驗(yàn)進(jìn)一步表明, 與大西洋尼諾現(xiàn)象(圖8)不同, WStr和CO2的熱效應(yīng)對(duì)類(lèi)尼諾升溫的季節(jié)變化都有貢獻(xiàn), 但是前者絕對(duì)占主導(dǎo)地位(圖8)。

    2.2.3 凈熱通量

    凈表面熱通量對(duì)赤道大西洋MLT的正異常有次要貢獻(xiàn)。在大西洋尼諾正位相時(shí)期, 海洋的失熱阻礙了類(lèi)大西洋尼諾升溫現(xiàn)象的發(fā)生(圖7g)。當(dāng)把變化響應(yīng)分解為WStr和CO2的熱效應(yīng)時(shí), 我們發(fā)現(xiàn), 如果只在WStr效應(yīng)的作用下, 海洋的凈熱通量會(huì)減少更多(圖8)。這證明了Bjerknes正反饋在大西洋尼諾的生成中起到正面作用。

    圖8 覆蓋實(shí)驗(yàn)各元素在峰值期對(duì)熱収支平衡方程各項(xiàng)貢獻(xiàn)的平均

    2.2.4 緯向和經(jīng)向平流

    緯向和經(jīng)向熱輸運(yùn)在熱帶大西洋與其他收支項(xiàng)相比貢獻(xiàn)很少。覆蓋的實(shí)驗(yàn)進(jìn)一步證實(shí), 無(wú)論是在類(lèi)尼諾升溫還是大西洋尼諾中, 赤道東大西洋緯向熱輸運(yùn)異常變化的原因都以WStr效應(yīng)為主(圖8)。

    3 總結(jié)和討論

    研究結(jié)果表明, 熱帶大西洋對(duì)全球變暖的響應(yīng)在北半球秋季和夏季相反。前者與大西洋尼諾正位相特征非常相似, 包括赤道信風(fēng)減弱, 赤道東大西洋變暖加劇, 同時(shí)溫躍層加深; 與之相反, 后者的特征是與北半球秋季相反的大西洋尼諾負(fù)位相對(duì)應(yīng)。

    大西洋尼諾和全球變暖下的類(lèi)尼諾升溫的特征及相關(guān)的形成機(jī)制都非常相似。特別是, 在這兩種情況下, Bjerknes正反饋對(duì)赤道東大西洋溫度異常都有非常大貢獻(xiàn)。然而, 它們之間的區(qū)別也很明顯:

    1) 類(lèi)尼諾升溫下西風(fēng)異常主要集中在大西洋東部, 而大西洋尼諾時(shí)主要集中在大西洋中部;

    2) 除風(fēng)應(yīng)力外, CO2的熱效應(yīng)對(duì)類(lèi)尼諾升溫的變化也有一定影響。

    與之前對(duì)熱帶太平洋和印度洋的研究相比[13], 我們發(fā)現(xiàn)了三個(gè)熱帶海洋之間有趣的相似和不同之處。WStr效應(yīng)在三個(gè)熱帶海洋表面溫度變化的季節(jié)變化中都占絕對(duì)主導(dǎo)地位, 這表明Bjerknes正反饋在調(diào)節(jié)熱帶氣候系統(tǒng)中的重要性。盡管有這些相似之處, 它們之間的區(qū)別也很明顯:

    1) 太平洋和印度洋的升溫分別類(lèi)厄爾尼諾和印度洋偶極子的正位相, 而大西洋在不同季節(jié)的升溫呈現(xiàn)大西洋尼諾的不同位相。

    2) 太平洋, CO2的熱效應(yīng)對(duì)赤道海溫變化的貢獻(xiàn)很大(47%)[13], 而在印度洋和大西洋, CO2的熱效應(yīng)對(duì)赤道海溫變化的貢獻(xiàn)非常小。

    [1] Richter I, Behera S K, Masumoto Y, et al. Multiple causes of interannual sea surface temperature variability in the equatorial Atlantic Ocean[J]. Nat Geosci, 2012, 6: 43-47.

    [2] Lutz K, Rathmann J, Jacobeit J. Classification of warm and cold water events in the eastern tropical Atlantic Ocean[J]. Atmos Sci Lett, 2013, 14: 102-106.

    [3] Keenlyside N S, Latif M. Understanding equatorial atlantic interannual variability[J]. J Clim, 2007, 20: 131-142.

    [4] Philander S G H. Unusual conditions in the tropical Atlantic Ocean in 1984[J]. Nature, 1986, 322: 236-238.

    [5] Bjerknes J. Atmospheric teleconnections from the equatorial Pacific[J]. Mon Wea Rev, 1969, 97: 163-172.

    [6] Folland C K, Colman a W, Rowell D P, et al. Predictability of northeast Brazil rainfall and real-time forecast skill, 1987-98[J]. J Clim, 2001, 14: 1937-1958.

    [7] Suarez M J, Schopf P S. A delayed action oscillator for ENSO[J]. J Atmos Sci, 1998, 45: 3283-3287.

    [8] Battisti D S. Dynamics and thermodynamics of a warming event in a coupled tropical atmosphere-ocean model[J]. J Atmos Sci, 1988, 45: 2889-2919.

    [9] Zebiak S E. Air-Sea Interaction in the Equatorial Atlantic Region[J]. J Clim, 1993, (6): 1567-1586.

    [10] Carton J A, Huang B. Warm Events in the Tropical Atlantic[J]. J Phys Oceanogr, 1994, 24: 888-903.

    [11] Lu J, Zhao B. The role of oceanic feedback in the climate response to doubling CO2[J]. J Clim, 2012, 25: 7544-7563.

    [12] Xie S P. On the genesis of the equatorial annual cycle[J]. J Clim, 1994, 7: 2008-2013.

    [13] Luo Y, Lu J, Liu F, et al. Understanding the El Ni?o-like oceanic response in the tropical Pacific to global warming[J]. Clim Dyn, 2014, 45: 1945-1964.

    [14] DiNezio P N, Clement A C, Vecchi G A, et al. Climate response of the equatorial Pacific to global warming[J]. J Clim, 2009, 22: 4873-4892.

    [15] Collins M, An S I, Cai W. et al. The impact of global warming on the tropical Pacific Ocean and El Ni?o[J]. Nature Geosci, 2010, 3: 391-397.

    [16] Tokinaga H, Xie S P. Weakening of the equatorial Atlantic cold tongue over the past six decades[J]. Nature Geoscience, 2010, 4(4): 222-226.

    [17] Xie S P, Deser C, Vecchi G A, et al. Global warming pattern formation: sea surface temperature and rainfall[J]. J Clim, 2010, 23: 966-986.

    [18] Liu F, Luo Y, Lu J, et al. Response of the tropical Pacific Ocean to El Ni?o versus global warming[J]. Climate Dynamics, Springer Berlin Heidelberg, 2017, 48(3-4): 935-956.

    Response of the tropical Atlantic Ocean to global warming in CESM

    DONG Wen-jing, LIU Fu-kai, LUO Yi-yong

    (College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China)

    There has been little discussion of the changes in the tropical Atlantic attributable to global warming (GW). Here, we used the Community Earth System Model (CESM), and found that GW induces an Atlantic Nino-like warming pattern in boreal autumn and a cold Atlantic Nino-like warming pattern in boreal summer in the tropical Atlantic Ocean. Applying an overriding technique to the CESM’s ocean component (version 2 of the Parallel Ocean Program), this study isolated the effects of wind stress, wind speed, and thermal warming on the response of the tropical Atlantic Ocean to global warming, and, further, investigated the similarities and differences in formation mechanisms underlying the changes in the tropical Atlantic Ocean between the Atlantic Ni?o and GW. Results showed that the underlying formation mechanism was quite similar in both cases, with wind-thermocline-SST feedback being the leading mechanism producing the anomalous warming over the eastern tropics in both cases. There were also some significant distinctions between them: 1) GW featured an overall reduction in amplitude and phase advance of the seasonal cycle of the mesosphere-lower troposphere (MLT) while the Atlantic Nino caused less disruption of the seasonal cycle; 2) westerly wind anomalies were most prominent in the eastern equatorial Atlantic during the warm phase of GW while in the central equatorial Atlantic during the warm Atlantic Nino; 3) air-sea interactions made some contribution to GW, however to a lesser degree than the wind stress effect.

    Tropical Atlantic Ocean; global warming; Atlantic Ni?o; Atlantic Ni?o-like response

    Jan. 26, 2019

    P47

    A

    1000-3096(2020)06-0012-10

    10.11759/hykx20190126001

    2019-01-26;

    2019-03-18

    太平洋環(huán)流對(duì)全球變暖的響應(yīng)的形成機(jī)制及其氣候效應(yīng)(41376009)

    [The formation mechanism of Pacific circulation in response to global warming and its climate effects, No. 41376009]

    董文靜(1991-), 女, 山東濱州人,漢族,學(xué)生,碩士, 主要從事海氣相互作用研究, 電話(huà):18765951916, E-mail: mirror_ouc@163.com

    (本文編輯: 李曉燕)

    猜你喜歡
    正位赤道大西洋
    “雪龍2”號(hào)過(guò)赤道,救下一船人!
    軍事文摘(2024年2期)2024-01-10 01:58:44
    滌綸正位/換位紡復(fù)合結(jié)構(gòu)紗的試制及性能評(píng)價(jià)
    關(guān)于鐵路貨車(chē)制動(dòng)管系漏泄原因分析及對(duì)策建議
    距赤道越遠(yuǎn)越喜愛(ài)黃色
    最熱的地方不在赤道
    大西洋海雀,你真倔
    飛越大西洋
    金融委的意義在于為金融“正位”
    暢游于大西洋彼岸
    需要給合作學(xué)習(xí)“正位”
    欧美久久黑人一区二区| 国产欧美日韩综合在线一区二区| 久久精品亚洲av国产电影网| 国产国语露脸激情在线看| 色婷婷av一区二区三区视频| 中出人妻视频一区二区| 久久久久国内视频| 欧美精品一区二区免费开放| e午夜精品久久久久久久| 操出白浆在线播放| 熟女少妇亚洲综合色aaa.| 两个人看的免费小视频| 中文字幕av电影在线播放| 搡老乐熟女国产| 亚洲国产欧美网| 一边摸一边抽搐一进一出视频| 制服人妻中文乱码| 黄色怎么调成土黄色| 成人精品一区二区免费| 精品久久久久久电影网| 欧美日韩乱码在线| 亚洲av美国av| 91麻豆精品激情在线观看国产 | 久久久国产成人免费| 天堂动漫精品| 另类亚洲欧美激情| 一区二区日韩欧美中文字幕| 国产在线观看jvid| 91成人精品电影| av天堂在线播放| 日本五十路高清| 97人妻天天添夜夜摸| 欧美人与性动交α欧美软件| 国产免费男女视频| 午夜日韩欧美国产| 黑人欧美特级aaaaaa片| 午夜精品在线福利| 人人妻人人爽人人添夜夜欢视频| 91九色精品人成在线观看| 极品人妻少妇av视频| 叶爱在线成人免费视频播放| 欧美黑人精品巨大| 一进一出抽搐gif免费好疼 | 天天躁日日躁夜夜躁夜夜| 黑人欧美特级aaaaaa片| 99国产精品99久久久久| 一边摸一边做爽爽视频免费| 两个人免费观看高清视频| 亚洲成人国产一区在线观看| 久久99一区二区三区| 成年人黄色毛片网站| 国产麻豆69| √禁漫天堂资源中文www| 中出人妻视频一区二区| 成人永久免费在线观看视频| 亚洲精品成人av观看孕妇| 新久久久久国产一级毛片| 国产一卡二卡三卡精品| 久久精品国产亚洲av香蕉五月 | 午夜免费成人在线视频| 18禁国产床啪视频网站| 成在线人永久免费视频| 久久久久久久午夜电影 | 男女床上黄色一级片免费看| 国产亚洲精品久久久久5区| 大香蕉久久成人网| 国产单亲对白刺激| 一本大道久久a久久精品| 建设人人有责人人尽责人人享有的| 女人高潮潮喷娇喘18禁视频| 亚洲男人天堂网一区| 香蕉丝袜av| 中文欧美无线码| 99久久精品国产亚洲精品| 亚洲情色 制服丝袜| 高清欧美精品videossex| 999久久久精品免费观看国产| 三上悠亚av全集在线观看| 午夜福利一区二区在线看| 成人亚洲精品一区在线观看| 免费看十八禁软件| 久久99一区二区三区| 亚洲视频免费观看视频| 精品一品国产午夜福利视频| 一个人免费在线观看的高清视频| 欧美色视频一区免费| 多毛熟女@视频| 咕卡用的链子| 一区二区三区激情视频| 久久久国产欧美日韩av| 90打野战视频偷拍视频| 十八禁人妻一区二区| 99久久国产精品久久久| 久久人人爽av亚洲精品天堂| 成年版毛片免费区| 人人妻,人人澡人人爽秒播| 亚洲第一青青草原| 91av网站免费观看| 69精品国产乱码久久久| svipshipincom国产片| 黑人操中国人逼视频| 身体一侧抽搐| 久久精品aⅴ一区二区三区四区| 日韩 欧美 亚洲 中文字幕| 黄色成人免费大全| 久久国产精品人妻蜜桃| 国产亚洲欧美精品永久| 国产精品亚洲av一区麻豆| 大片电影免费在线观看免费| 精品午夜福利视频在线观看一区| 大片电影免费在线观看免费| 操美女的视频在线观看| 亚洲欧美激情综合另类| 首页视频小说图片口味搜索| 两人在一起打扑克的视频| 色婷婷av一区二区三区视频| 后天国语完整版免费观看| 18在线观看网站| 夜夜爽天天搞| 成人亚洲精品一区在线观看| 国产精品久久久av美女十八| 天堂俺去俺来也www色官网| 久久精品国产a三级三级三级| 人人妻人人爽人人添夜夜欢视频| av网站免费在线观看视频| 日韩一卡2卡3卡4卡2021年| 精品国产乱子伦一区二区三区| 在线十欧美十亚洲十日本专区| 欧美成人午夜精品| 他把我摸到了高潮在线观看| 午夜福利在线观看吧| 久久久国产成人精品二区 | 国产成人欧美在线观看 | 精品一区二区三卡| 大型av网站在线播放| e午夜精品久久久久久久| 一级毛片高清免费大全| 精品少妇久久久久久888优播| 无限看片的www在线观看| 国产视频一区二区在线看| 另类亚洲欧美激情| 色精品久久人妻99蜜桃| 中文字幕最新亚洲高清| 精品一区二区三区视频在线观看免费 | 国产男女超爽视频在线观看| 免费在线观看亚洲国产| 黄片播放在线免费| 丝袜在线中文字幕| 亚洲精品国产精品久久久不卡| 大香蕉久久成人网| av国产精品久久久久影院| 99国产精品一区二区三区| av线在线观看网站| 叶爱在线成人免费视频播放| 久久久国产一区二区| 国精品久久久久久国模美| 国产成人欧美| 午夜福利在线免费观看网站| 午夜老司机福利片| av电影中文网址| 老熟女久久久| 久久精品国产亚洲av高清一级| 国产单亲对白刺激| 热99国产精品久久久久久7| 好男人电影高清在线观看| 久久久水蜜桃国产精品网| 亚洲一区高清亚洲精品| 精品卡一卡二卡四卡免费| 中文字幕av电影在线播放| 欧美 亚洲 国产 日韩一| 怎么达到女性高潮| 男女午夜视频在线观看| 精品国产国语对白av| 亚洲av片天天在线观看| 91麻豆av在线| 国产黄色免费在线视频| 最近最新中文字幕大全免费视频| 国产精品自产拍在线观看55亚洲 | 最近最新中文字幕大全免费视频| 久久香蕉国产精品| 亚洲国产看品久久| 亚洲精品在线观看二区| 天堂俺去俺来也www色官网| 成人特级黄色片久久久久久久| 久久中文看片网| 亚洲性夜色夜夜综合| 国产单亲对白刺激| 欧美日韩瑟瑟在线播放| 亚洲国产看品久久| 国产黄色免费在线视频| xxx96com| 日韩免费av在线播放| 国产欧美亚洲国产| 国产精品亚洲av一区麻豆| 久久香蕉精品热| 精品久久久久久久毛片微露脸| 大片电影免费在线观看免费| 乱人伦中国视频| 亚洲av欧美aⅴ国产| 国产成人av激情在线播放| av天堂在线播放| 麻豆av在线久日| 国产1区2区3区精品| 日本欧美视频一区| 国产精品美女特级片免费视频播放器 | 两性夫妻黄色片| 欧美乱色亚洲激情| 人人妻人人爽人人添夜夜欢视频| 精品午夜福利视频在线观看一区| 午夜两性在线视频| 国产午夜精品久久久久久| 日韩有码中文字幕| av不卡在线播放| 亚洲国产毛片av蜜桃av| 久久精品91无色码中文字幕| 91老司机精品| 新久久久久国产一级毛片| 狠狠婷婷综合久久久久久88av| 国产精品亚洲av一区麻豆| 一区二区日韩欧美中文字幕| 欧美精品高潮呻吟av久久| 国产精品自产拍在线观看55亚洲 | 亚洲情色 制服丝袜| 中文亚洲av片在线观看爽 | 国产一卡二卡三卡精品| 久久草成人影院| 国产精品1区2区在线观看. | 亚洲中文av在线| 波多野结衣一区麻豆| 人妻丰满熟妇av一区二区三区 | 亚洲精品粉嫩美女一区| 亚洲人成伊人成综合网2020| 少妇猛男粗大的猛烈进出视频| 三上悠亚av全集在线观看| 一区福利在线观看| 在线观看免费日韩欧美大片| 国精品久久久久久国模美| 亚洲va日本ⅴa欧美va伊人久久| 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 成人国语在线视频| 国产精品99久久99久久久不卡| 久热这里只有精品99| 美女高潮到喷水免费观看| 国产成人一区二区三区免费视频网站| 国精品久久久久久国模美| 建设人人有责人人尽责人人享有的| 国产91精品成人一区二区三区| 久久性视频一级片| 亚洲成人手机| 国产精品久久久久久精品古装| 精品亚洲成国产av| 久久久久国内视频| 黄色毛片三级朝国网站| 精品一区二区三卡| 十八禁人妻一区二区| 18禁裸乳无遮挡动漫免费视频| 精品熟女少妇八av免费久了| 18禁美女被吸乳视频| 亚洲av熟女| 三级毛片av免费| 久99久视频精品免费| 成人手机av| 国产精品 国内视频| 在线永久观看黄色视频| 欧美在线黄色| 两个人免费观看高清视频| 久久久精品免费免费高清| 国产又爽黄色视频| 久久精品aⅴ一区二区三区四区| 欧美精品高潮呻吟av久久| 国产亚洲av高清不卡| 国产精品98久久久久久宅男小说| 亚洲伊人色综图| 中文字幕另类日韩欧美亚洲嫩草| 熟女少妇亚洲综合色aaa.| 丰满饥渴人妻一区二区三| 精品人妻熟女毛片av久久网站| 免费女性裸体啪啪无遮挡网站| 老司机午夜十八禁免费视频| www日本在线高清视频| 久久精品国产亚洲av香蕉五月 | 亚洲精品在线美女| 精品免费久久久久久久清纯 | 美女高潮到喷水免费观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩制服丝袜自拍偷拍| 91精品三级在线观看| 99久久综合精品五月天人人| 无遮挡黄片免费观看| 女警被强在线播放| 中文字幕制服av| 欧美日韩成人在线一区二区| 国产精品国产高清国产av | 91精品国产国语对白视频| 亚洲片人在线观看| 国产成人一区二区三区免费视频网站| 午夜成年电影在线免费观看| av线在线观看网站| a在线观看视频网站| 看黄色毛片网站| 精品国产一区二区三区四区第35| 亚洲人成电影免费在线| 免费一级毛片在线播放高清视频 | 欧美日韩福利视频一区二区| 黄色丝袜av网址大全| 亚洲 国产 在线| 久久久久久久精品吃奶| 黄色a级毛片大全视频| 伊人久久大香线蕉亚洲五| 国产精品二区激情视频| 免费不卡黄色视频| 国产99白浆流出| 18禁美女被吸乳视频| 国产日韩一区二区三区精品不卡| 精品电影一区二区在线| 亚洲欧美一区二区三区黑人| 18禁黄网站禁片午夜丰满| 国产精品乱码一区二三区的特点 | 中出人妻视频一区二区| 老司机亚洲免费影院| 日韩欧美三级三区| 免费看十八禁软件| 搡老熟女国产l中国老女人| 国产精品国产av在线观看| 久久久国产成人免费| 老司机深夜福利视频在线观看| 黄色女人牲交| 操美女的视频在线观看| 麻豆乱淫一区二区| 狠狠婷婷综合久久久久久88av| 两个人看的免费小视频| 女人高潮潮喷娇喘18禁视频| 成年人黄色毛片网站| 好看av亚洲va欧美ⅴa在| 国产真人三级小视频在线观看| 久久精品亚洲av国产电影网| 久久ye,这里只有精品| 中文字幕精品免费在线观看视频| 久久精品亚洲精品国产色婷小说| 免费观看精品视频网站| 老司机在亚洲福利影院| 水蜜桃什么品种好| 久久久久国产精品人妻aⅴ院 | av国产精品久久久久影院| 9191精品国产免费久久| 久久精品亚洲熟妇少妇任你| 久久热在线av| 国产免费男女视频| 成人特级黄色片久久久久久久| videos熟女内射| 亚洲欧美一区二区三区黑人| 男男h啪啪无遮挡| 久久精品aⅴ一区二区三区四区| 一进一出抽搐动态| 精品国产亚洲在线| 老司机影院毛片| 国产人伦9x9x在线观看| 韩国精品一区二区三区| 99riav亚洲国产免费| 精品福利观看| 免费高清在线观看日韩| 精品一区二区三区av网在线观看| 巨乳人妻的诱惑在线观看| 久久午夜亚洲精品久久| 亚洲欧美精品综合一区二区三区| 免费观看人在逋| 亚洲在线自拍视频| 啦啦啦免费观看视频1| 大香蕉久久成人网| 欧美日韩瑟瑟在线播放| 99国产精品免费福利视频| 日本撒尿小便嘘嘘汇集6| 老熟妇乱子伦视频在线观看| 女人久久www免费人成看片| 最新在线观看一区二区三区| 欧美成人午夜精品| 精品福利永久在线观看| 国产高清videossex| 婷婷精品国产亚洲av在线 | 国产精品乱码一区二三区的特点 | 精品电影一区二区在线| 一区二区三区精品91| 一区二区三区激情视频| 制服诱惑二区| 一进一出抽搐gif免费好疼 | 成人永久免费在线观看视频| 高清黄色对白视频在线免费看| 热99国产精品久久久久久7| 高清黄色对白视频在线免费看| 色老头精品视频在线观看| 国产精品久久久久久人妻精品电影| www日本在线高清视频| 国产成人精品无人区| 欧美日韩瑟瑟在线播放| 亚洲一区二区三区欧美精品| 中文亚洲av片在线观看爽 | 女人高潮潮喷娇喘18禁视频| 精品欧美一区二区三区在线| 精品视频人人做人人爽| 色尼玛亚洲综合影院| 久久精品国产清高在天天线| 国产精品一区二区免费欧美| 国产一区二区三区综合在线观看| 色播在线永久视频| а√天堂www在线а√下载 | 香蕉国产在线看| 国产99白浆流出| 国产欧美日韩一区二区精品| 久久精品国产a三级三级三级| 99热网站在线观看| 久久狼人影院| 久久国产精品人妻蜜桃| 国产一区二区激情短视频| 看片在线看免费视频| 亚洲性夜色夜夜综合| 欧美性长视频在线观看| 亚洲国产欧美一区二区综合| 久久久久久免费高清国产稀缺| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 精品久久蜜臀av无| 国产高清激情床上av| av天堂久久9| 99国产精品免费福利视频| 日日摸夜夜添夜夜添小说| 免费高清在线观看日韩| 夜夜夜夜夜久久久久| 国产精品 欧美亚洲| 亚洲av成人不卡在线观看播放网| 黄片大片在线免费观看| 久久精品亚洲精品国产色婷小说| www.999成人在线观看| 久久草成人影院| av天堂在线播放| 99精品久久久久人妻精品| 欧美精品亚洲一区二区| 露出奶头的视频| 在线观看舔阴道视频| 久久国产精品大桥未久av| 日韩免费高清中文字幕av| cao死你这个sao货| 侵犯人妻中文字幕一二三四区| 欧美亚洲 丝袜 人妻 在线| 久久久精品区二区三区| 午夜精品久久久久久毛片777| 亚洲全国av大片| bbb黄色大片| 国产日韩欧美亚洲二区| 久久久久国产一级毛片高清牌| 国产av一区二区精品久久| 国产高清videossex| 一二三四社区在线视频社区8| 国产成人精品久久二区二区免费| 久久中文字幕一级| 欧美激情久久久久久爽电影 | 亚洲成av片中文字幕在线观看| 成人精品一区二区免费| 高清黄色对白视频在线免费看| 国产精品一区二区精品视频观看| 国产欧美日韩一区二区精品| 91成人精品电影| 午夜成年电影在线免费观看| 国产欧美日韩精品亚洲av| 亚洲中文日韩欧美视频| 亚洲欧美一区二区三区黑人| 亚洲一区中文字幕在线| 制服诱惑二区| 免费黄频网站在线观看国产| 国产精品欧美亚洲77777| 国产精品国产av在线观看| 黑人欧美特级aaaaaa片| 在线永久观看黄色视频| 在线观看日韩欧美| 欧美激情久久久久久爽电影 | 好男人电影高清在线观看| 后天国语完整版免费观看| 一级黄色大片毛片| 狠狠婷婷综合久久久久久88av| 丁香六月欧美| 亚洲av美国av| 正在播放国产对白刺激| 免费少妇av软件| 操美女的视频在线观看| 老熟女久久久| 女警被强在线播放| 久久午夜综合久久蜜桃| 色在线成人网| 亚洲欧美激情综合另类| 亚洲全国av大片| 免费在线观看黄色视频的| 国产99白浆流出| 国产色视频综合| 久久久精品区二区三区| 在线观看免费日韩欧美大片| 国产亚洲一区二区精品| 嫁个100分男人电影在线观看| 变态另类成人亚洲欧美熟女 | 丝袜在线中文字幕| 如日韩欧美国产精品一区二区三区| 97人妻天天添夜夜摸| 久久久久久久精品吃奶| 亚洲va日本ⅴa欧美va伊人久久| 国产免费现黄频在线看| 黄色怎么调成土黄色| 亚洲第一青青草原| 91大片在线观看| 国产精品国产av在线观看| 丰满人妻熟妇乱又伦精品不卡| 精品一品国产午夜福利视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久久久久精品古装| 大香蕉久久成人网| 久久久精品区二区三区| 制服诱惑二区| 老汉色∧v一级毛片| 99riav亚洲国产免费| 久久香蕉精品热| 亚洲五月婷婷丁香| 欧美日韩黄片免| 欧美成狂野欧美在线观看| 12—13女人毛片做爰片一| 18禁黄网站禁片午夜丰满| av在线播放免费不卡| 夜夜躁狠狠躁天天躁| 性色av乱码一区二区三区2| 一二三四在线观看免费中文在| 人人妻人人澡人人看| 一区二区日韩欧美中文字幕| 真人做人爱边吃奶动态| 电影成人av| 伊人久久大香线蕉亚洲五| 日韩欧美在线二视频 | 99国产精品一区二区三区| 成年人黄色毛片网站| 亚洲国产毛片av蜜桃av| 亚洲欧美日韩高清在线视频| 国产亚洲精品久久久久5区| 亚洲精品在线美女| 精品欧美一区二区三区在线| 国产一区二区三区综合在线观看| 欧美激情久久久久久爽电影 | 桃红色精品国产亚洲av| 91成人精品电影| 欧美久久黑人一区二区| 性色av乱码一区二区三区2| 婷婷成人精品国产| 久久久久国产一级毛片高清牌| 一二三四社区在线视频社区8| www日本在线高清视频| 午夜精品久久久久久毛片777| 精品视频人人做人人爽| 成人精品一区二区免费| 一边摸一边抽搐一进一小说 | x7x7x7水蜜桃| 亚洲国产毛片av蜜桃av| 丝袜人妻中文字幕| 一级毛片女人18水好多| 亚洲精品一卡2卡三卡4卡5卡| 成年人黄色毛片网站| 丝瓜视频免费看黄片| 高清av免费在线| 国产精品免费一区二区三区在线 | 久久久久久久国产电影| 精品欧美一区二区三区在线| 他把我摸到了高潮在线观看| 人妻丰满熟妇av一区二区三区 | 国产精品亚洲av一区麻豆| 亚洲国产看品久久| 亚洲成国产人片在线观看| 国产淫语在线视频| 久久热在线av| 午夜福利免费观看在线| 欧美激情久久久久久爽电影 | 国产91精品成人一区二区三区| 少妇猛男粗大的猛烈进出视频| 成在线人永久免费视频| 变态另类成人亚洲欧美熟女 | 亚洲美女黄片视频| 国产精品免费大片| 国产不卡av网站在线观看| 看免费av毛片| 国产男女超爽视频在线观看| 国产免费男女视频| 亚洲精华国产精华精| 国产一区二区三区视频了| 18禁黄网站禁片午夜丰满| 精品久久久久久,| 国产精品永久免费网站| 高清视频免费观看一区二区| 中文字幕高清在线视频| 亚洲欧美一区二区三区黑人| 一级黄色大片毛片| 国产熟女午夜一区二区三区| 新久久久久国产一级毛片| 亚洲国产中文字幕在线视频| 国产精品二区激情视频| 王馨瑶露胸无遮挡在线观看| 亚洲av日韩在线播放| 在线观看免费视频日本深夜| 制服诱惑二区| 男女午夜视频在线观看| 国产精品 欧美亚洲| 欧美日本中文国产一区发布| 色综合婷婷激情| 男女免费视频国产| 丝瓜视频免费看黄片| 精品久久久久久久久久免费视频 | 一级作爱视频免费观看| 人妻 亚洲 视频| 中文字幕精品免费在线观看视频| 日韩有码中文字幕| 亚洲,欧美精品.| 老熟妇乱子伦视频在线观看| 国产成人免费观看mmmm|