• <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í)“正位”
    变态另类丝袜制服| 在线a可以看的网站| 黄色毛片三级朝国网站| 好看av亚洲va欧美ⅴa在| 99国产精品一区二区三区| 天天一区二区日本电影三级| 成熟少妇高潮喷水视频| 法律面前人人平等表现在哪些方面| 精品日产1卡2卡| 91麻豆av在线| 国内久久婷婷六月综合欲色啪| 午夜精品一区二区三区免费看| 成人一区二区视频在线观看| 亚洲国产欧美网| 日韩国内少妇激情av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产高清在线一区二区三| 久久婷婷成人综合色麻豆| 日本五十路高清| 亚洲精品久久成人aⅴ小说| 亚洲男人天堂网一区| 99国产综合亚洲精品| 老鸭窝网址在线观看| 99热6这里只有精品| 丁香六月欧美| 一区二区三区高清视频在线| 亚洲精品粉嫩美女一区| 久久久国产成人免费| 久久久精品大字幕| 免费在线观看影片大全网站| 在线免费观看的www视频| 中亚洲国语对白在线视频| 亚洲男人的天堂狠狠| 亚洲最大成人中文| 草草在线视频免费看| 国产成人精品久久二区二区91| 国产精品亚洲av一区麻豆| 国产精品野战在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲国产欧美一区二区综合| 丰满人妻熟妇乱又伦精品不卡| 国产精品国产高清国产av| 亚洲国产精品999在线| 听说在线观看完整版免费高清| 午夜激情福利司机影院| 99久久99久久久精品蜜桃| 国产三级中文精品| av在线天堂中文字幕| 亚洲欧美激情综合另类| 久久久久免费精品人妻一区二区| 亚洲国产精品999在线| 日韩欧美免费精品| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 国产蜜桃级精品一区二区三区| 亚洲av熟女| 精品久久蜜臀av无| 免费无遮挡裸体视频| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 久久久国产成人精品二区| 久久精品亚洲精品国产色婷小说| www日本黄色视频网| 男女之事视频高清在线观看| 成人三级做爰电影| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| av超薄肉色丝袜交足视频| 亚洲精品国产一区二区精华液| 欧美日韩瑟瑟在线播放| av欧美777| 一级片免费观看大全| 欧洲精品卡2卡3卡4卡5卡区| 欧美激情久久久久久爽电影| 国产真人三级小视频在线观看| 18禁美女被吸乳视频| 亚洲第一欧美日韩一区二区三区| 国产精品一区二区三区四区免费观看 | 亚洲一区高清亚洲精品| 国产伦在线观看视频一区| 色在线成人网| 亚洲一码二码三码区别大吗| 国产精品久久电影中文字幕| 999久久久精品免费观看国产| 亚洲熟妇中文字幕五十中出| 国产成人aa在线观看| 日本在线视频免费播放| 小说图片视频综合网站| 亚洲精品国产精品久久久不卡| 天堂√8在线中文| 禁无遮挡网站| 日本一本二区三区精品| 亚洲专区中文字幕在线| 亚洲七黄色美女视频| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| 久久精品aⅴ一区二区三区四区| 99国产精品一区二区蜜桃av| 亚洲精品国产一区二区精华液| 亚洲精品一卡2卡三卡4卡5卡| 91麻豆精品激情在线观看国产| 91国产中文字幕| 精品久久久久久久久久免费视频| 欧美一区二区精品小视频在线| 久久久久久久久久黄片| bbb黄色大片| 人人妻人人看人人澡| 国产男靠女视频免费网站| 久久中文字幕一级| 亚洲av电影不卡..在线观看| 国产黄片美女视频| 九色成人免费人妻av| 亚洲人与动物交配视频| 欧美日韩国产亚洲二区| 色综合亚洲欧美另类图片| 欧美乱妇无乱码| 99热这里只有是精品50| 人人妻,人人澡人人爽秒播| 亚洲九九香蕉| 国产人伦9x9x在线观看| 18禁观看日本| 大型av网站在线播放| 在线播放国产精品三级| 大型黄色视频在线免费观看| а√天堂www在线а√下载| 欧洲精品卡2卡3卡4卡5卡区| 老司机靠b影院| 琪琪午夜伦伦电影理论片6080| 亚洲欧美激情综合另类| 在线观看66精品国产| 国产高清视频在线观看网站| 欧美黄色淫秽网站| 国产亚洲av高清不卡| 国产欧美日韩精品亚洲av| 很黄的视频免费| 一本大道久久a久久精品| 久久精品夜夜夜夜夜久久蜜豆 | 搡老熟女国产l中国老女人| 啦啦啦免费观看视频1| 色精品久久人妻99蜜桃| 在线国产一区二区在线| 国产不卡一卡二| 嫁个100分男人电影在线观看| 欧洲精品卡2卡3卡4卡5卡区| av免费在线观看网站| ponron亚洲| svipshipincom国产片| 亚洲av成人一区二区三| 精品一区二区三区视频在线观看免费| 日韩欧美三级三区| 欧美zozozo另类| 欧美日韩一级在线毛片| 日本一本二区三区精品| 亚洲中文字幕一区二区三区有码在线看 | 人人妻人人看人人澡| 伊人久久大香线蕉亚洲五| 19禁男女啪啪无遮挡网站| 国产av又大| 国产探花在线观看一区二区| 不卡一级毛片| 亚洲熟妇中文字幕五十中出| 毛片女人毛片| 国产精品永久免费网站| 757午夜福利合集在线观看| 精华霜和精华液先用哪个| 欧美日韩中文字幕国产精品一区二区三区| 1024手机看黄色片| 欧美极品一区二区三区四区| 熟妇人妻久久中文字幕3abv| 伊人久久大香线蕉亚洲五| 亚洲中文av在线| 国产成人影院久久av| 可以免费在线观看a视频的电影网站| 欧美最黄视频在线播放免费| 岛国在线观看网站| 妹子高潮喷水视频| av欧美777| 男女之事视频高清在线观看| 国内精品久久久久久久电影| 亚洲 欧美一区二区三区| 精品国产超薄肉色丝袜足j| 欧美黑人精品巨大| 黄频高清免费视频| 中文字幕精品亚洲无线码一区| 一级毛片女人18水好多| 全区人妻精品视频| 国产免费av片在线观看野外av| 99国产精品一区二区蜜桃av| 久久草成人影院| 不卡av一区二区三区| av国产免费在线观看| 国产亚洲精品久久久久5区| 亚洲欧美激情综合另类| 在线观看www视频免费| 黄色视频,在线免费观看| 国产精品永久免费网站| 波多野结衣高清无吗| 99国产精品一区二区蜜桃av| 久久久久久久久中文| 国产主播在线观看一区二区| xxxwww97欧美| 悠悠久久av| 51午夜福利影视在线观看| 久久国产精品人妻蜜桃| 色噜噜av男人的天堂激情| 1024手机看黄色片| 黄色成人免费大全| 性色av乱码一区二区三区2| 日本一本二区三区精品| 男人舔女人下体高潮全视频| 欧美色欧美亚洲另类二区| 国产精品免费视频内射| 天天添夜夜摸| 一级毛片女人18水好多| 美女午夜性视频免费| 亚洲成av人片免费观看| 中文资源天堂在线| 日韩三级视频一区二区三区| 亚洲天堂国产精品一区在线| 国产一区在线观看成人免费| 亚洲人成伊人成综合网2020| 妹子高潮喷水视频| 精品免费久久久久久久清纯| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月| 日本一区二区免费在线视频| 熟女电影av网| 欧美日韩亚洲国产一区二区在线观看| 人人妻,人人澡人人爽秒播| 午夜福利在线观看吧| 日韩欧美在线乱码| 亚洲精品美女久久久久99蜜臀| 91麻豆av在线| 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 亚洲av中文字字幕乱码综合| 精品乱码久久久久久99久播| 国产精品久久电影中文字幕| 亚洲美女黄片视频| 精品欧美国产一区二区三| 日韩国内少妇激情av| 亚洲熟妇熟女久久| 久久人妻av系列| 午夜影院日韩av| 日本一二三区视频观看| 成人高潮视频无遮挡免费网站| 免费高清视频大片| 国产又黄又爽又无遮挡在线| 国内揄拍国产精品人妻在线| 少妇的丰满在线观看| 老汉色∧v一级毛片| 欧美zozozo另类| 一区二区三区国产精品乱码| 国产精品一区二区三区四区久久| 在线播放国产精品三级| av在线播放免费不卡| 午夜老司机福利片| 我的老师免费观看完整版| 夜夜爽天天搞| 成人高潮视频无遮挡免费网站| 国产午夜精品久久久久久| 亚洲国产中文字幕在线视频| 免费看a级黄色片| 香蕉国产在线看| 中文资源天堂在线| 国内揄拍国产精品人妻在线| 欧美乱妇无乱码| 一本大道久久a久久精品| 深夜精品福利| 特大巨黑吊av在线直播| 日本黄大片高清| 欧美日韩乱码在线| 国产成人av激情在线播放| 亚洲成人精品中文字幕电影| 午夜精品久久久久久毛片777| 日韩国内少妇激情av| 深夜精品福利| 蜜桃久久精品国产亚洲av| 午夜福利高清视频| 一本久久中文字幕| 免费在线观看视频国产中文字幕亚洲| 免费高清视频大片| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久亚洲av鲁大| 男插女下体视频免费在线播放| 91九色精品人成在线观看| 久久中文看片网| 99热6这里只有精品| 麻豆成人av在线观看| 舔av片在线| 国产aⅴ精品一区二区三区波| 欧美+亚洲+日韩+国产| 亚洲五月天丁香| 午夜日韩欧美国产| 亚洲激情在线av| 欧美3d第一页| 老熟妇仑乱视频hdxx| 免费人成视频x8x8入口观看| 久久精品亚洲精品国产色婷小说| 女生性感内裤真人,穿戴方法视频| 一卡2卡三卡四卡精品乱码亚洲| 动漫黄色视频在线观看| 制服丝袜大香蕉在线| 日日爽夜夜爽网站| 亚洲av美国av| 国产精品 欧美亚洲| 美女大奶头视频| 日本一二三区视频观看| 可以在线观看毛片的网站| 亚洲成人国产一区在线观看| 亚洲avbb在线观看| 午夜a级毛片| 久久久久久亚洲精品国产蜜桃av| 午夜免费成人在线视频| 午夜激情福利司机影院| 嫩草影视91久久| 欧美成狂野欧美在线观看| 欧美日韩精品网址| 免费搜索国产男女视频| 黄色a级毛片大全视频| 国产欧美日韩一区二区三| 校园春色视频在线观看| 午夜福利18| 午夜激情福利司机影院| 国产三级中文精品| 美女 人体艺术 gogo| 免费在线观看影片大全网站| 毛片女人毛片| 999精品在线视频| 丝袜人妻中文字幕| 日韩中文字幕欧美一区二区| 男人的好看免费观看在线视频 | 最近最新中文字幕大全免费视频| 老司机深夜福利视频在线观看| 亚洲狠狠婷婷综合久久图片| 波多野结衣巨乳人妻| АⅤ资源中文在线天堂| 成人欧美大片| 亚洲国产欧洲综合997久久,| 国产久久久一区二区三区| 啦啦啦免费观看视频1| 成人国语在线视频| 美女大奶头视频| 美女 人体艺术 gogo| 亚洲国产高清在线一区二区三| 搡老妇女老女人老熟妇| 夜夜躁狠狠躁天天躁| 全区人妻精品视频| 又紧又爽又黄一区二区| 日本精品一区二区三区蜜桃| 无人区码免费观看不卡| 国产精品一区二区精品视频观看| 国产一区二区三区在线臀色熟女| 国产在线精品亚洲第一网站| 村上凉子中文字幕在线| 18禁黄网站禁片午夜丰满| 两性午夜刺激爽爽歪歪视频在线观看 | 18禁黄网站禁片免费观看直播| 露出奶头的视频| 欧美中文综合在线视频| av福利片在线| 亚洲精品av麻豆狂野| 国产高清激情床上av| 亚洲人成电影免费在线| av免费在线观看网站| or卡值多少钱| 国产高清videossex| 国产精品影院久久| 嫁个100分男人电影在线观看| 又黄又粗又硬又大视频| 九色成人免费人妻av| 亚洲精品国产一区二区精华液| 每晚都被弄得嗷嗷叫到高潮| 亚洲av第一区精品v没综合| netflix在线观看网站| 美女黄网站色视频| 欧美另类亚洲清纯唯美| 国产精品影院久久| 麻豆av在线久日| 又黄又粗又硬又大视频| 丰满人妻一区二区三区视频av | 白带黄色成豆腐渣| 极品教师在线免费播放| 丁香六月欧美| 久久天堂一区二区三区四区| 精品日产1卡2卡| 亚洲精品中文字幕在线视频| 婷婷精品国产亚洲av在线| 亚洲国产欧美一区二区综合| 国产成人av激情在线播放| 18禁黄网站禁片免费观看直播| 精品久久久久久久末码| 在线国产一区二区在线| 国产免费av片在线观看野外av| 欧美色欧美亚洲另类二区| 久久婷婷成人综合色麻豆| 精品久久久久久久末码| 特大巨黑吊av在线直播| 亚洲专区中文字幕在线| 日本 欧美在线| 免费在线观看视频国产中文字幕亚洲| 亚洲片人在线观看| 很黄的视频免费| 悠悠久久av| 免费看十八禁软件| 色老头精品视频在线观看| 欧美日韩乱码在线| 男男h啪啪无遮挡| 婷婷六月久久综合丁香| 免费在线观看成人毛片| 精品少妇一区二区三区视频日本电影| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 久久中文看片网| 国产激情久久老熟女| 亚洲精品一区av在线观看| 在线观看日韩欧美| 亚洲欧美精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 国产在线精品亚洲第一网站| 国产精品,欧美在线| 国产激情久久老熟女| 美女免费视频网站| 好男人在线观看高清免费视频| 亚洲一区中文字幕在线| 国产精品日韩av在线免费观看| 观看免费一级毛片| 男女午夜视频在线观看| 久9热在线精品视频| 亚洲精品在线美女| 成人午夜高清在线视频| 9191精品国产免费久久| avwww免费| 99热6这里只有精品| 我的老师免费观看完整版| 国产精品久久视频播放| 国产精品野战在线观看| 中文字幕人妻丝袜一区二区| 精品午夜福利视频在线观看一区| 亚洲精品中文字幕在线视频| 国产精品 国内视频| 欧美精品啪啪一区二区三区| 国内毛片毛片毛片毛片毛片| 欧美3d第一页| 久久久久久人人人人人| 久久久久国内视频| 激情在线观看视频在线高清| 中文字幕熟女人妻在线| 波多野结衣巨乳人妻| 哪里可以看免费的av片| 亚洲天堂国产精品一区在线| 丁香欧美五月| 三级毛片av免费| 精品一区二区三区视频在线观看免费| 一边摸一边抽搐一进一小说| 三级男女做爰猛烈吃奶摸视频| 夜夜爽天天搞| 婷婷亚洲欧美| 亚洲国产欧美人成| xxx96com| 免费av毛片视频| 欧美中文综合在线视频| 女人高潮潮喷娇喘18禁视频| 日韩有码中文字幕| 亚洲精品中文字幕在线视频| 一本久久中文字幕| 欧美日韩精品网址| 九九热线精品视视频播放| 亚洲 欧美一区二区三区| 国产私拍福利视频在线观看| 级片在线观看| 亚洲成人国产一区在线观看| 亚洲天堂国产精品一区在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美激情综合另类| 动漫黄色视频在线观看| 俺也久久电影网| 又紧又爽又黄一区二区| 麻豆av在线久日| 亚洲欧美精品综合一区二区三区| 亚洲五月天丁香| 亚洲av五月六月丁香网| 亚洲精品中文字幕一二三四区| 人人妻人人看人人澡| 91成年电影在线观看| 色综合婷婷激情| 人人妻,人人澡人人爽秒播| 国产成人aa在线观看| 91字幕亚洲| 中文字幕精品亚洲无线码一区| www.自偷自拍.com| 久9热在线精品视频| av免费在线观看网站| 母亲3免费完整高清在线观看| 午夜激情av网站| 中亚洲国语对白在线视频| 三级男女做爰猛烈吃奶摸视频| 欧美黑人巨大hd| 欧美成人午夜精品| 色精品久久人妻99蜜桃| 亚洲色图av天堂| av超薄肉色丝袜交足视频| 亚洲欧美精品综合一区二区三区| 日韩欧美免费精品| www日本在线高清视频| 亚洲国产精品成人综合色| 国产精品一区二区三区四区免费观看 | 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 伦理电影免费视频| 国产av一区二区精品久久| 狂野欧美激情性xxxx| 国产真人三级小视频在线观看| 精品久久蜜臀av无| 精品人妻1区二区| 亚洲成a人片在线一区二区| 亚洲美女黄片视频| 亚洲欧美激情综合另类| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 97碰自拍视频| 午夜激情福利司机影院| 色在线成人网| 成人三级做爰电影| 美女大奶头视频| 国产精品亚洲av一区麻豆| 日本免费一区二区三区高清不卡| 免费在线观看黄色视频的| 嫩草影视91久久| 午夜免费激情av| 精品电影一区二区在线| 欧美色欧美亚洲另类二区| 无遮挡黄片免费观看| 老司机福利观看| 色在线成人网| 亚洲欧美精品综合久久99| 日韩精品青青久久久久久| 两人在一起打扑克的视频| 人人妻人人澡欧美一区二区| 久久草成人影院| 男女做爰动态图高潮gif福利片| 观看免费一级毛片| 一二三四社区在线视频社区8| 黑人欧美特级aaaaaa片| 伦理电影免费视频| 久久婷婷人人爽人人干人人爱| 亚洲第一欧美日韩一区二区三区| avwww免费| 中文在线观看免费www的网站 | 可以免费在线观看a视频的电影网站| 色在线成人网| 99热这里只有精品一区 | 亚洲成人精品中文字幕电影| 成人亚洲精品av一区二区| 国产精品一及| 免费观看精品视频网站| 日本精品一区二区三区蜜桃| 免费观看精品视频网站| 久久久久久久久久黄片| 看免费av毛片| 国产精品一区二区三区四区久久| 午夜日韩欧美国产| 欧美日韩乱码在线| 国产精品1区2区在线观看.| 高清毛片免费观看视频网站| 精品福利观看| 啪啪无遮挡十八禁网站| 在线观看66精品国产| 国产欧美日韩一区二区精品| 亚洲精品一区av在线观看| 成人亚洲精品av一区二区| 亚洲男人天堂网一区| 黄色 视频免费看| 久久这里只有精品中国| 老司机深夜福利视频在线观看| 久久久精品大字幕| 婷婷精品国产亚洲av在线| 免费高清视频大片| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 黄片小视频在线播放| 欧美中文日本在线观看视频| 动漫黄色视频在线观看| 国内毛片毛片毛片毛片毛片| 夜夜躁狠狠躁天天躁| 午夜a级毛片| 免费在线观看成人毛片| 久久精品91蜜桃| 日韩av在线大香蕉| 91九色精品人成在线观看| 亚洲激情在线av| 国产精品久久电影中文字幕| 91字幕亚洲| 欧美人与性动交α欧美精品济南到| 亚洲国产高清在线一区二区三| 搡老妇女老女人老熟妇| 精品久久久久久久久久免费视频| 他把我摸到了高潮在线观看| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av在线| 免费在线观看亚洲国产| 亚洲精品在线观看二区| 老司机午夜十八禁免费视频| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 97碰自拍视频| av福利片在线| 久久久久久九九精品二区国产 | 可以免费在线观看a视频的电影网站| 麻豆久久精品国产亚洲av| 日韩av在线大香蕉| 超碰成人久久| 露出奶头的视频| 中文字幕久久专区| 亚洲,欧美精品.|