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

    熱帶太平洋海溫異常對北極海冰的可能影響

    2017-07-07 13:19:36羅蕊祁莉張文君何金海
    海洋學報 2017年7期
    關鍵詞:北極海海溫海冰

    羅蕊,祁莉*,張文君,何金海

    (1.南京信息工程大學 氣象災害教育部重點實驗室/氣候與環(huán)境變化國際合作聯(lián)合實驗室/氣象災害預報預警與評估協(xié)同創(chuàng)新中心,江蘇 南京 210044)

    ?

    熱帶太平洋海溫異常對北極海冰的可能影響

    羅蕊1,祁莉1*,張文君1,何金海1

    (1.南京信息工程大學 氣象災害教育部重點實驗室/氣候與環(huán)境變化國際合作聯(lián)合實驗室/氣象災害預報預警與評估協(xié)同創(chuàng)新中心,江蘇 南京 210044)

    本文利用1950-2015年間Hadley環(huán)流中心海冰和海溫資料及NCEP/NCAR再分析資料,研究了熱帶太平洋海溫異常對北極海冰的可能影響,并從大氣環(huán)流和凈表面熱通量兩個角度探討了可能的物理機制。結果表明,在ENSO事件發(fā)展年的夏、秋季節(jié),EP型與CP型El Nio事件與北極海冰異常的聯(lián)系無明顯信號。而La Nia事件期間北極海冰出現(xiàn)顯著異常,并且EP型與CP型La Nia之間存在明顯差異。EP型La Nia發(fā)生時,北極地區(qū)巴倫支海、喀拉海關鍵區(qū)海冰異常減少,CP型La Nia事件則對應著東西伯利亞海、楚科奇海地區(qū)海冰異常增加。在EP型La Nia發(fā)展年的夏、秋季節(jié),熱帶太平洋海溫異常通過遙相關波列,使得巴倫支海、喀拉海海平面氣壓為負異常并與中緯度氣壓正異常共同構成類似AO正位相的結構,形成的風場異常有利于北大西洋暖水的輸入,同時造成暖平流,偏高的水汽含量進一步加強了凈表面熱通量收入,使得巴倫支海、喀拉海海冰異常減少。而在CP型La Nia發(fā)展年的夏季,東西伯利亞海、楚科奇海關鍵區(qū)受其東側氣旋式環(huán)流的影響,以異常北風分量占主導,將海冰從極點附近由北向南輸送到關鍵區(qū),海冰異常增加,而凈表面熱通量的作用較小。

    北極海冰;EP/CP La Nia;大氣環(huán)流;凈表面熱通量

    1 引言

    北極海冰作為氣候系統(tǒng)五大圈層之一冰雪圈的重要組成部分,是一個強大的冷源,它通過隔絕海洋和大氣的直接熱量輸送,改變了地表輻射平衡和熱量平衡,進而影響大氣環(huán)流和全球氣候的變化。IPCC第五次評估報告(AR5)[1]指出,在全球變暖的背景下20世紀后期北極海冰出現(xiàn)了年代際尺度的急劇消融,這個現(xiàn)象引起了全球廣泛的關注,因為北極海冰的消融使得北冰洋相對海冰顏色較深的開闊水域面積大大增加,改變了海洋表層的能量收支,這對北極甚至全球氣候起著非常重要的調制作用[2]。裸露的洋面增加了夏季對太陽輻射的吸收,同時減緩了秋冬季節(jié)北極海冰的結冰速度,因此,北極海冰的這種溫度—反照率的正反饋機制將全球變暖效應放大,稱為“北極放大效應”(Arctic amplification)[3]。正是由于北極海冰對氣候變化的高度敏感性,它成為了氣候變化重要的指示器和記憶器[4]。ENSO作為全球氣候系統(tǒng)中海氣相互作用年際變化的最強信號,它引起的赤道太平洋海溫異常,不僅在緯向上對Walker環(huán)流產生影響,直接造成熱帶太平洋、印度洋地區(qū)的天氣、氣候異常,同時還在經(jīng)向上通過三圈環(huán)流和遙相關等方式,間接影響中高緯地區(qū)[5-9]。

    那么,北極海冰和ENSO作為全球氣候系統(tǒng)中至關重要的兩個外強迫源,它們之間是否存在一定的相互聯(lián)系呢?已有的研究主要分為兩個部分。一部分以北極海冰為立足點,分析這一強大的冷性外強迫源對熱帶海洋的影響。從海冰和ENSO的周期相關來看,北極海冰面積變化具有準兩年和準四年的變化周期,這與ENSO事件的周期有較好的一致性[10],當Nio3區(qū)海表面溫度(SST)滯后喀拉海、拉普捷夫海和東西伯利亞海海冰面積指數(shù)16個月時,兩者達到最大負相關[11]。北極海冰的變化存在明顯的區(qū)域性和季節(jié)性[12],冬春季節(jié)喀拉海、巴倫支海海冰面積均與春季白令海海冰面積呈反向變化關系,與巴芬灣戴維斯海峽海冰面積也存在相反的變化趨勢[13-15]。不同區(qū)域不同季節(jié)的海冰變化對熱帶海洋的影響存在很大的差異,冬季挪威海區(qū)域海冰對厄爾尼諾的產生具有強迫作用,而夏季巴芬灣區(qū)域海冰的變化則是對厄爾尼諾響應的結果[16]。格陵蘭海、挪威海和巴倫支海冬季海冰面積偏大能夠使赤道海溫持續(xù)偏高,拉普帖夫海、東西伯利亞海和鄂霍次克海夏秋季海冰面積偏大能夠使翌年夏秋季赤道海溫偏低,這種持續(xù)影響可能與極渦和東亞冷空氣有關[17],海冰可以通過冷空氣活動和低渦西移,引起副熱帶高壓位置和強度的變化,影響了熱帶太平洋地區(qū)的信風強弱,從而與熱帶太平洋海溫的異常聯(lián)系在一起[13-15,18-19];也會通過影響中亞氣壓來改變赤道環(huán)流,影響印度洋和太平洋赤道海溫[19-20];還有研究把海冰與遙相關(比如AO、NAO)波列聯(lián)系起來,通過影響中緯度地區(qū)的風暴軸事件或者西風異常來影響海溫[21-27]。另一部分的研究工作從ENSO出發(fā),把熱帶太平洋海溫異常看做外強迫源,研究它對北極海冰的影響[28-31]。ENSO事件海溫異常會激發(fā)大氣的PNA(Pacific North American pattern)遙相關,使得中緯度的西風偏強,阿留申低壓和北美高壓都加強,造成太平洋風暴軸偏南偏東[17],或者與NAO聯(lián)系在一起[30],正位相的ENSO事件會改變局地Ferrel環(huán)流,調控平均經(jīng)向熱通量,造成楚科奇和波弗特海的海冰減少[31]。有研究表明,中部型El Nio(La Nia)對北極海冰產生的影響,主要以AO作為聯(lián)系熱帶和北極的紐帶[32],或者通過“大氣橋”來影響極地渦旋,進而改變海冰狀況[33]。印太暖池區(qū)域海溫異常與北極海冰變化的聯(lián)系,則強調印太暖池區(qū)域海溫異常對北極海冰快速減少的作用,同樣也是通過AO作為中間橋梁,來解釋熱帶海洋對北極海冰產生影響的機制[34]。

    從國內外的研究進展可知,現(xiàn)有的研究工作大多側重分析北極海冰對熱帶太平洋的影響,且較為成熟;但是對熱帶太平洋影響北極海冰的研究工作較少。此外,近20年來一種新的ENSO事件頻繁發(fā)生,它的異常海溫中心位于赤道中太平洋[9],本文中稱之為CP型El Nio/La Nia,傳統(tǒng)型ENSO(異常海溫中心位于熱帶東太平洋)稱之為EP型El Nio/La Nia。大量的研究指出CP型ENSO與EP型ENSO對熱帶、中高緯地區(qū)的大氣環(huán)流和氣候異常的影響存在顯著的差異[35-39]。那么,它們對北極海冰的影響是否也不同?其中的物理機制是什么?本文以熱帶太平洋海溫異常的不同模態(tài)入手,分析不同海溫模態(tài)對北極海冰的可能影響,以及其中的可能物理機制,為更好理解北極海冰銳減及其氣候效應奠定一定的基礎。

    2 資料與方法

    本文所用到的資料主要有Hadley環(huán)流中心月平均的海冰密集度和海溫資料HadISST(Hadley Centre Sea Ice and Sea Surface Temperature Data Set),空間分辨率1.0°×1.0°,海冰密集度單位為%,表示單位網(wǎng)格內海冰覆蓋的百分比[40],需要進行說明的是,雖然該海冰數(shù)據(jù)在1979年之前的可靠性受到爭議,但本文發(fā)現(xiàn)僅使用1979年后的資料進行分析,不影響本文的定性結論。為了去除全球變暖的信號以及北極放大效應,本文中所有與海冰相關的分析,均進行了去趨勢處理。大氣環(huán)流資料來自NCEP/NCAR(National Centers for Environmental Prediction /National Center for Atmospheric Research)月平均再分析資料[41],水平分辨率2.5°×2.5°。包括海平面氣壓、位勢高度、風場、相對濕度、2 m氣溫、向上向下的短波和長波輻射、潛熱通量和感熱通量。本文主要采用回歸、合成分析的方法,并用t檢驗對合成的結果進行顯著性檢驗。本文的研究時段為1950年1月至2015年12月。

    北極海冰具有明顯季節(jié)變化的特點,每年夏季,隨著太陽直射點的向北移動,北極地區(qū)進入極晝,接收到的太陽輻射全年最多,海冰進入融冰期,直到9月海冰密集度達到極小值,秋季隨著氣溫的降低,北極地區(qū)慢慢進入結冰期,直到翌年的3月,海冰密集度又恢復到極大值,因此每年3-9月,海冰處于逐漸減少的過程,夏季減少得尤其快;9月到翌年3月海冰逐漸增加,秋季增加最快,所以夏、秋季節(jié)是海冰年際變率最大的季節(jié)[42],同時也是年際變化尺度上ENSO事件的發(fā)展期。如果海冰在其年際變率最大的季節(jié),同時有ENSO事件開始發(fā)展,它是否會對海冰的變化產生影響?因此,為了分析ENSO對海冰的可能影響,本文重點關注ENSO事件發(fā)展年夏季和秋季北極海冰的變化,北極地區(qū)地理分布如圖1所示。

    圖1 北極地區(qū)地理分布示意圖Fig.1 Sketch map of Arctic geographic distribution

    關于ENSO事件的分類,不同的研究選用不同的指數(shù)得到不同的分類結果[9,35-36]。大部分工作較多地關注EP和CP型El Nio的有效區(qū)分以及它們不同的氣候效應。王磊等[37]和Zhang[39]發(fā)現(xiàn)La Nia事件也存在截然不同的EP型和CP型,它們不同的海氣耦合特征在大西洋地區(qū)引起的NAO響應幾乎相反,這可能意味著EP型和CP型La Nia事件將對北極海冰的不同影響。因此,本文參考何珊珊等[38]和王磊等[37]的分類方法,依據(jù)海溫異常的空間分布型對ENSO事件進行分類得到如表1所示的分類結果。

    圖2給出了根據(jù)表1中4類ENSO事件分類合成得到的發(fā)展年夏、秋季節(jié)海溫異常的空間分布。EP型El Nio發(fā)展年夏季(圖2a),海溫正距平已經(jīng)出現(xiàn),主要位于赤道東太平洋90°W附近,最大值超過1.2℃,且海溫異常自東向西逐漸擴展減弱;到了秋季(圖2b),海溫異常大值中心強度繼續(xù)加強,范圍顯著西伸,緯向跨度增大,但其主體仍然位于赤道東太平洋。而CP型El Nio在發(fā)展年的夏季(圖2c)海溫異常主體位于赤道中太平洋160°W附近,且赤道東、西太平洋海溫距平在零值附近,秋季(圖2d)海溫正異常進一步加強,達到0.8℃以上,但大值中心位置幾乎沒有變化。兩類La Nia的情況與El Nio類似,EP型La Nia發(fā)展年夏季(圖2e),海溫最大負距平為-0.8℃,主體位于赤道東太平洋90°W附近,向西逐漸減弱,秋季(圖2f)海溫負異常加強且緯向范圍明顯擴大,大值中心向西延伸到110°W附近。CP型La Nia發(fā)展年夏季(圖2g),海溫異常強度雖然與EP型一致,但在空間分布上主體位于赤道中太平洋150°W偏西,秋季時(圖2h)其強度迅速加強達到-1.2℃以上,位置幾乎沒有變化??梢姡还苁荅l Nio還是La Nia,EP型事件的海溫異常主要分布在赤道東太平洋地區(qū),而CP型事件異常海溫偏西,表明了本文采取的ENSO事件的分類方式是合理有效的。并且,EP型和CP型海溫異常中心的緯向差異在發(fā)展年的夏季和秋季均十分明顯,這將對局地以及熱帶外的大氣環(huán)流造成不同的影響。

    表1 1950-2015年間4類不同ENSO事件的分類

    圖2 4類ENSO事件發(fā)展年夏季(a, c, e, g)和秋季(b, d, f, h)赤道太平洋地區(qū)海表面溫度異常(SSTA, 單位:℃)的空間分布Fig.2 Spatial distribution of summer and autumn mean SST anomalies (SSTA,unit:℃) over the tropical Pacific during the development phase of four types of ENSO eventsa、b為EP型El Nio,c、d為CP型El Nio,e、f為EP型La Nia,g、h為CP型La Nia;打點區(qū)域表示通過90%顯著性檢驗;兩類La Nia合成圖中紅色框為圖7和圖8中海溫回歸位勢高度場時的關鍵區(qū)范圍EP El Nio(a, b),CP El Nio(c, d),EP La Nia(e, f),CP La Nia(g, h); doting indicates anomalies exceeding 90% confidence level; red boxes in two types of La Nia events indicate the key regions of SSTA regressing to geopotential height in Fig.7 and Fig.8

    3 4類ENSO事件與北極海冰的聯(lián)系

    4類ENSO事件與北極海冰是否存在聯(lián)系?它們截然不同的海溫異常分布型能否引起差異顯著的北極海冰響應?

    圖3給出了北極區(qū)域(60°~90°N,0°~360°)平均的海冰標準化指數(shù)與同期4類ENSO事件強度的散點圖,以美國氣候預測中心(Climate Prediction Center,CPC)Nio3和Nio4指數(shù)分別表征EP和CP型事件強度。在ENSO事件發(fā)展年的夏季(圖3a),個別La Nia事件發(fā)展較晚,海溫異常中心有所偏移,因此關鍵區(qū)的SSTA偏弱,甚至符號相反(例如2005/2006和1983/1984年),El Nio事件時海冰以正距平居多,占總數(shù)的比例為10/17,但EP型El Nio(紅色)和CP型El Nio(橙色)對應的海冰變化沒有表現(xiàn)出明顯的區(qū)別,它們中各有2次和3次事件中海冰異常偏少。La Nia事件中,EP型La Nia(藍色)發(fā)展年夏季,8次事件有7次北極海冰均為負距平,而CP型La Nia時海冰有正有負。然而,到發(fā)展年秋季(圖3b),6/7的CP型La Nia對應北極海冰正距平,而8次EP型La Nia事件中有5次海冰明顯偏少,3次雖然數(shù)值較小,但均為負距平,可以認為兩類La Nia對應北極海冰的變化出現(xiàn)了顯著相反的異常特征,即EP型La Nia時北極海冰偏少而CP型La Nia時海冰異常偏多;但北極海冰對El Nio事件的響應并不穩(wěn)定,無論是EP型還是CP型El Nio海冰異常均有正有負。

    圖3 4類ENSO事件發(fā)展年夏季(a)和秋季(b)北極海冰區(qū)域(60°~90°N,0°~360°)平均標準化指數(shù)與同期CPC Nio指數(shù)散點圖分布。Nio3指數(shù)代表EP型事件,Nio4指數(shù)代表CP型事件Fig.3 Scatter of standardization index of regional average of Arctic sea ice (60°-90°N,0°-360°) and CPC Nio index in summer (a) and autumn (b) during the development phase of four types of ENSO e-vents. Nio3 and Nio 4 index indicate EP and CP events respectively

    圖4 EP/CP型El Nio發(fā)展年夏季(a,c)和秋季(b,d)北極海冰密集度異常(%)合成結果Fig.4 Composites of Arctic sea ice concentration anomalies (%) in summer (a,c) and autumn (b,d) of EP/CP El Nio during their development phasea、b為EP型El Nio,c、d為CP型El Nio;圖中填色區(qū)域通過了90%顯著性檢驗EP El Nio (a,b),CP El Nio (c,d);shading indicates anomalies exceeding 90% confidence level

    圖5 EP/CP型La Nia發(fā)展年夏季(a,c)和秋季(b,d)北極海冰密集度異常(%)合成結果Fig.5 Composites of Arctic sea ice concentration anomalies (%) in summer (a,c) and autumn (b,d) of EP/CP La Nia during their development phasea、b為EP型El Nio,c、d為CP型El Nio;圖中填色區(qū)域通過了90%顯著性檢驗EP El Nio(a,b),CP El Nio(c,d);shading indicates anomalies exceeding 90% confidence level

    圖6 EP型La Nia事件發(fā)展年夏季(a,c)和秋季(b,d)海平面氣壓異常(a,b,單位:hPa)以及2 m氣溫(單位:℃)和風場異常(c,d,單位:m/s)的合成結果Fig.6 Composites of sea level pressure (SLP) anomalies (a,b, unit:hPa) , 2 m air temperature (unit:℃) and wind anomalies (c,d, unit:m/s) in summer (a,c) and autumn (b,d) during EP La Nia development phase紅色扇形框表示海冰異常偏少的關鍵區(qū)范圍,網(wǎng)格區(qū)域(a,b)和粗實線內區(qū)域(c,d)表示通過90%顯著性檢驗Red fans indicate key areas that sea ice is less than normal, griddings(a,b)and areas within the thick line(c,d)indicate anomalies exceeding 90% confidence level

    圖7 EP La Nia事件發(fā)展年夏季(a,c)和秋季(b,d)500 hPa位勢高度異常合成(a,b)和回歸(c,d)結果對比Fig.7 Contrast of 500 hPa geopotential height composite (a,b) and regression (c,d) in summer (a,c) and autumn (b,d) during EP La Nia development phase回歸序列為圖2e、f紅色框內海溫異常的區(qū)域平均;紅色扇形框表示海冰異常偏少的關鍵區(qū)范圍;網(wǎng)格區(qū)域表示通過90%顯著性檢驗Regression sequence is the SSTA regional average of red boxes in Fig.2e and 2f;red fans indicate key areas that sea ice is less than normal, griddings indicate anomalies exceeding 90% confidence level

    圖8 EP型La Nia事件發(fā)展年夏季(a)和秋季(b)凈表面熱通量(Qnet)的合成結果Fig.8 Composites of net surface heat flux in summer (a) and autumn (b) during EP La Nia development phase紅色扇形框表示海冰異常偏少的關鍵區(qū)范圍,網(wǎng)格區(qū)域表示通過90%顯著性檢驗Red fans indicate key areas that sea ice is less than normal, griddings indicate anomalies exceeding 90% confidence level

    圖9 EP型La Nia發(fā)展年夏季(a)和秋季(b)海冰關鍵區(qū)(70°~90°N,30°~90°E)區(qū)域平均凈表面熱通量各項貢獻Fig.9 Contributions of each items of net surface heat flux averaged in the sea ice key area (70°-90°N,30°-90°E) in summer (a) and autumn (b) during EP La Nia development phase

    圖10 EP型La Nia事件發(fā)展年夏季(a)和秋季(b)850 hPa比濕的合成結果Fig.10 Composites of 850 hPa specific humidity in summer (a) and autumn (b) during EP La Nia development phase紅色扇形框表示海冰異常偏少的關鍵區(qū)范圍,網(wǎng)格區(qū)域表示通過90%顯著性檢驗Red fans indicate key areas that sea ice is less than normal, griddings indicate anomalies exceeding 90% confidence level

    圖12 CP La Nia事件發(fā)展年夏季(a,c)和秋季(b,d)500 hPa位勢高度異常合成(a,b)和回歸(c,d)結果對比Fig.12 Contrast of 500 hPa geopotential height composite (a,b) and regression (c,d) in summer (a,c) and autumn (b,d) during CP La Nia development phase回歸序列為圖2g、h紅色框內海溫異常的區(qū)域平均;紅色扇形框表示海冰異常偏多的關鍵區(qū)范圍,網(wǎng)格區(qū)域表示通過90%顯著性檢驗Regression sequence is the SSTA regional average of red boxes in Fig.2g,h; red fans indicate key areas that sea ice is more than normal,griddings indicate anomalies exceeding 90% confidence level

    圖14 CP型La Nia發(fā)展年夏季(a)和秋季(b)海冰關鍵區(qū)(70°~90°N,150°E~160°W)區(qū)域平均凈表面熱通量各項貢獻Fig.14 Contributions of each items of net surface heat flux averaged in the sea ice key area (70°-90°N,150°E-160°W) in summer (a) and autumn (b) during CP La Nia development phase

    4 兩類La Nia事件與北極海冰聯(lián)系的可能物理機制

    北極海冰的變化與局地大氣環(huán)流和凈表面熱通量異常關系密切,下面將主要從這兩個角度對兩類La Nia事件對應的北極海冰變化進行對比分析。

    圖6的分析表明關鍵區(qū)附近的氣旋式環(huán)流異常對海冰的減少有重要作用,這種環(huán)流異常在對流層中層也表現(xiàn)明顯(圖7a,7b)。圖7c和7d(圖2e、2f紅框范圍)給出了夏秋季節(jié)海溫異常對500 hPa位勢高度的回歸(為方便分析,回歸結果均已進行反號處理),對比合成和回歸結果可以發(fā)現(xiàn),無論是夏季還是秋季,在副熱帶中東太平洋至北美東部地區(qū)均呈現(xiàn)出了PNA位勢高度距平波列,海冰關鍵區(qū)位勢高度以負異常為主,秋季,該負異常與中緯度出現(xiàn)的位勢高度正異常共同形成了類似AO正位相氣壓北負南正的結構,使海冰關鍵區(qū)及其周邊地區(qū)出現(xiàn)氣旋式環(huán)流,氣溫增加,有利于海冰的減少。需要進行說明的是,當用逐月的海溫異常對海冰密集度進行超前滯后相關分析時,發(fā)現(xiàn)6-11月北極海冰關鍵區(qū)的異常變化滯后海溫異常約1-2個月(圖略),這種滯后性同樣體現(xiàn)在大氣環(huán)流的變化上,但由于本文研究的是夏秋季節(jié)平均海冰的異常變化,因此1-2月的滯后性不明顯。以上分析表明,熱帶太平洋海溫異??梢酝ㄟ^遙相關波列影響北極高緯度地區(qū)的大氣環(huán)流異常,在巴倫支海、喀拉海形成氣旋式環(huán)流,促進該地區(qū)海冰的異常偏少,進一步研究還發(fā)現(xiàn),該地區(qū)的局地凈表面熱通量也出現(xiàn)了異常變化。局地凈表面熱通量由以下公式給出:

    Qnet=Qdsw-Qusw+Qdlw-Qulw-Qlh-Qsh,

    (1)

    式中,方程右邊各項分別為向下(上)短波輻射Qdsw(Qusw),向下(上)長波輻射Qdlw(Qulw),潛熱通量Qlh和感熱通量Qsh。根據(jù)該公式計算得到EP型La Nia發(fā)展年夏季和秋季凈表面熱通量的異常分布,如圖8所示??梢姡贓P型La Nia發(fā)展年夏季(圖8a)巴倫支海、喀拉海關鍵區(qū)為顯著的熱通量凈收入,且通過90%顯著性檢驗,異常增加的熱通量造成了海冰減少。然而秋季(圖8b)凈表面熱通量反向變?yōu)樨摼嗥剑瑸榻忉屵@一現(xiàn)象,給出海冰關鍵區(qū)在EP型La Nia發(fā)展年期間公式(1)中右邊每一項的貢獻(圖9)。

    夏季關鍵區(qū)凈熱通量(Qnet)異常偏高(圖9a),海冰異常偏少,產生正貢獻的項為向上短波輻射(Qusw)、向下長波輻射(Qdlw)、潛熱通量(Qlh)和感熱通量(Qsh)。向上短波輻射(Qusw)取決于下表面反照率,但海冰-反照率強烈的正反饋機制,使得該項主要表現(xiàn)為海冰變化的結果,因此在進行海冰變化原因分析時暫不考慮;向下長波輻射(Qdlw)、潛熱通量(Qlh)和感熱通量(Qsh)主要與空氣中水汽含量、風場和溫度有關,水汽含量越高(低)、風場與溫度場形成暖(冷)平流,對凈熱通量貢獻越大(小),關鍵區(qū)海冰則異常減少(增加),這與上文分析的2 m溫度場與10 m風場的合成結果一致(圖6c,6d)。此外,夏季關鍵區(qū)比濕也表現(xiàn)為顯著的正異常(圖10),這與熱通量的分析是吻合的。秋季(圖9b),關鍵區(qū)凈熱通量(Qnet)為異常負距平,主要貢獻項為向上長波輻射(Qulw),該項取決于地表溫度,反映結冰量的多少,結冰量越大(小),釋放熱量使得下表面溫度越高(低),向上長波輻射越強(弱),則凈熱通量增加(減少),這說明秋季凈熱通量異常主要是海冰異常作用的結果,而非原因。Walsh和Johnson[46]認為在海冰結冰的季節(jié)(8月至翌1月),海冰對大氣環(huán)流的影響稍大于大氣對海冰的影響,而其他季節(jié)主要是大氣影響海冰,以上的分析驗證了這一說法,同時也證明了雖然秋季凈熱通量不能作為影響海冰變化的原因進行分析,但對夏季凈熱通量的分析仍然是必要的。

    在500 hPa位勢高度上(圖12a,12b),關鍵區(qū)氣壓異常與低層相似并有所加強,形成的正壓結構對北風異常的維持具有重要作用,北半球中高緯出現(xiàn)沿緯圈的位勢高度異常波列,與同期海溫回歸的結果十分相似(圖12c,12d),關鍵區(qū)都表現(xiàn)為夏季西高東低而秋季全區(qū)一致偏低的分布形式,這表明該區(qū)域的夏秋季節(jié)海冰異常與熱帶太平洋海溫異常是緊密聯(lián)系的,熱帶太平洋的海溫異常通過遙相關波列使得關鍵區(qū)出現(xiàn)氣壓異常,北風的維持與凈熱通量的顯著減少(圖13a),共同促進了關鍵區(qū)海冰的增加。秋季(圖13b),極區(qū)通過檢驗的區(qū)域較少,部分區(qū)域還出現(xiàn)了凈通量的正距平,為分析其原因,同樣計算了公式(1)右邊各項對該地區(qū)凈表面熱通量的貢獻,結果如圖14。夏季(圖14a)關鍵區(qū)凈表面熱通量為顯著的負異常,主要貢獻項分別為向上短波輻射(Qusw)和感熱通量(Qsh),上文分析已經(jīng)指出暫不考慮向上短波輻射項,但感熱通量數(shù)值較小,因此凈表面熱通量的作用較弱,同時冷暖平流合成結果中(圖11c,11d)通過檢驗的區(qū)域也較少,可見該區(qū)域海冰的異常增加可能主要來源于風的驅動作用,從圖11c、11d可以看到,關鍵區(qū)北風異常將極點附近的海冰由北向南輸送,因此引起海冰密集度的異常增加。秋季(圖14b)情況與EP型La Nia相似,凈熱通量變化是海冰作用的結果而非原因,故不再進行分析。

    5 結論

    本文利用Hadley環(huán)流中心月平均海冰密集度、海溫資料及NCEP/NCAR再分析資料,分析了4類不同ENSO事件與北極海冰變化之間的可能關系,得到結論如下:

    (1)通過合成分析發(fā)現(xiàn),在ENSO事件發(fā)展年的夏、秋季節(jié),EP與CP型El Nio事件與北極海冰異常沒有明顯聯(lián)系。而La Nia事件期間北極海冰出現(xiàn)顯著的異常,并且EP型與CP型La Nia之間存在明顯差異,在EP型La Nia事件中北極地區(qū)的巴倫支海、喀拉海關鍵區(qū)海冰異常減少,而CP型La Nia事件對應東西伯利亞海、楚科奇海關鍵區(qū)海冰異常增加。

    [1] Ofipcc W G I. Climate Change 2013: The Physical Science Basis.[J]. Contribution of Working, 2013, 43(22):866-871.

    [2] Bader J, Mesquita M D S, Hodges K I, et al. A review on Northern Hemisphere sea-ice, storminess and the North Atlantic Oscillation: Observations and projected changes[J]. Atmospheric Research, 2011, 101(4):809-834.

    [3] Screen J A, Ian S. The central role of diminishing sea ice in recent Arctic temperature amplification[J]. Nature, 2010, 464(7293):1334-1337.

    [4] 俞永強. 海-氣相互作用對我國氣候變化的影響[M]. 北京:氣象出版社, 2005.

    Yu Yongqiang. The Effect of Air-sea Interaction on the Climate Change in China[M]. Beijing:China Meteorological Press,2005.

    [5] Wallace J M, Gutzler D S. Teleconnections in the geopotential height field during the Northern Hemisphere winter[J]. Monthly Weather Review, 1981, 109(4):784-812

    [6] King J C. Currents of change: El Nio’s impact on climate and society[J]. Weather, 1997, 52(5):159-160.

    [7] Trenberth K E, Branstator G W, Karoly D, et al. Progress during TOGA in understanding and modeling global teleconnections associated with tropical sea surface temperatures[J]. Journal of Geophysical Research Oceans, 1998, 103(C7):14291-14324.

    [8] Ropelewski C F, Halpert M S. Global and regional scale precipitation patterns associated with the El Nio/Southern Oscillation[J]. Monthly Weather Review, 1987, 115(8):1606-1626.

    [9] Ashok K, Behera S K, Rao S A, et al. El Nio Modoki and its possible teleconnection[J]. Journal of Geophysical Research Oceans, 2007, 112(C11):C11007.

    [10] Gloersen P. Modulation of hemispheric sea-ice cover by ENSO events[J]. Nature International Weekly Journal of Science, 1995, 373(6514):503-506.

    [11] 朱艷峰, 陳隆勛. 北極海冰與ENSO事件在準四年時間尺度上的可能聯(lián)系[J]. 大氣科學, 2003, 27(5):834-846.

    Zhu Yanfeng, Chen Longxun. Potential relationship between Arctic sea-ice and ENSO events on timescales of quasi-quadrennial-year[J].Chinese Journal of Atmosphere Sciences, 2003, 27(5):834-846.

    [12] 汪代維, 楊修群. 北極海冰變化的時間和空間型[J]. 氣象學報, 2002, 60(2):129-138.

    Wang Daiwei,Yang Xiuqun. Temporal and spatial patterns of Arctic sea ice variations[J]. Acta Meteorologica Sinica, 2002, 60(2):129-138.

    [13] 武炳義, 黃榮輝, 高登義. 冬季喀拉海、巴倫支海海冰面積變化對后期北太平洋海溫的影響[J]. 氣候與環(huán)境研究, 1999(2):38-48.

    Wu Bingyi, Huang Ronghui, Gao Dengyi. The effect of variation of sea-ice extent in the Kara and Barents Seas in winter on SST in the Northern Pacific in the later period[J]. Climate and Environmental Research, 1999(2):38-48.

    [14] 武炳義, 高登義, 黃榮輝. 冬春季節(jié)北極海冰的年際和年代際變化[J]. 氣候與環(huán)境研究, 2000(3):249-258.

    Wu Bingyi, Gao Dengyi, Huang Ronghui. Interannual and interdecadal variations in Arctic sea-ice in spring and winter[J]. Climate and Environmental Research, 2000(3):249-258.

    [15] 武炳義, 高登義. 冬季格陵蘭、喀拉海和巴倫支海海冰年際變化與ENSO事件[J]. 科學通報, 1997(18):1979-1981.

    Wu Bingyi, Gao Dengyi. The relationship between interannual changes of winter sea ice in Greenland, Kara Sea, Barents Sea and ENSO[J]. Chinese Science Bulletin, 1997(18):1979-1981.

    [16] 胡增臻, 莊麗. 北極海冰與埃爾—尼諾[J]. 海洋預報, 1989(3):12-16.

    Hu Zengzhen, Zhuang Li. The sea ice of north pole (sinp) and El-Nio[J]. Marine Forecasts, 1989(3):12-16.

    [17] 錢步東, 范鐘秀, 彭公炳,等. 北極海冰與赤道東太平洋海溫的相互影響及其與El Nio的聯(lián)系[J]. 熱帶氣象學報, 1994(4):325-334.

    Qian Budong, Fan Zhongxiu, Peng Gongbing, et al. Interaction between Arctic sea ice and the equatorial eastern Pacific sea surface temperature and the relation to El Nio[J]. Journal of Tropical Meteorology, 1994(4):325-334.

    [18] 方之芳. 北半球副熱帶高壓與極地海冰的相互作用[J]. 科學通報, 1986(4):286-289.

    Fang Zhifang. Interaction between northern hemisphere sub-tropical high and polar sea ice[J]. Chinese Science Bulletin, 1986(4):286-289.

    [19] 李崇銀. 頻繁的強東亞大槽活動與El Nio發(fā)生[J]. 中國科學:B輯, 1988, 18(6):667-674.

    Li Chongyin. Frequent activities of strong east Asian trough and the occurrence of El Nio[J]. Science in China,Series B, 1988, 18(6):667-674.

    [20] Yasunari T. Zonally propagating modes of the global east-west circulation associated with the Southern Oscillat[J]. Journal of the Meteorological Society of Japan, 1985, 63(6):1013-1029.

    [21] Kvamsto Ng S P S D. Impact of Labrador sea-ice extent on the North Atlantic Oscillation[J]. International Journal of Climatology, 2004, 24(5): 603-612.

    [22] Seierstad I A, Bader J. Impact of projected future Arctic Sea Ice reduction on extratropical storminess and the NAO[J]. Climate Dynamics, 2009, 33(7): 937-943.

    [23] Honda M, Inoue J, Yamane S. Influence of low Arctic sea-ice minima on anomalously cold Eurasian winters[J]. Geophysical Research Letters, 2009, 36(8): L8707.

    [24] Outten S D, Esau I. A link between Arctic sea ice and recent cooling trends over Eurasia[J]. Climatic Change, 2012, 110(3): 1069-1075.

    [25] Peings Y, Magnusdottir G. Response of the wintertime Northern Hemisphere atmospheric circulation to current and projected Arctic sea ice decline: a numerical study with CAM5[J]. Journal of Climate, 2014, 27(1):244-264.

    [26] Francis J A, Chan W, Leathers D J, et al. Winter Northern Hemisphere weather patterns remember summer Arctic sea-ice extent[J]. Geophysical Research Letters, 2009, 36(7):157-163.

    [27] Wu Qigang, Zhang Xiangdong. Observed forcing-feedback processes between Northern Hemisphere atmospheric circulation and Arctic sea ice coverage[J]. Journal of Geophysical Research Atmospheres, 2010, 115(D14), doi:10.1029/2009JD013574.

    [28] Yu J Y, Zou Y. The enhanced drying effect of Central-Pacific El Nio on US winter[J]. Environmental Research Letters, 2013, 8(1):3865-3879.

    [29] Karori M A, Li J, Jin F F. The asymmetric influence of the two types of El Nio and La Nia on summer rainfall over southeast China[J]. Journal of Climate, 2013, 26(13):4567-4582.

    [30] Mysak L A, Ingram R G, Wang J,et al. The anomalous sea-ice extent in Hudson Bay, Baffin Bay and the Labrador Sea during three simultaneous NAO and ENSO episodes[J]. Atmosphere-Ocean, 1996, 34(2): 313-343.

    [31] Liu Jiping, Curry J A. Recent Arctic sea ice variability connections to the Arctic Oscillation and ENSO[J]. Geophysical Research Letters, 2004, 31: L9211.

    [32] 左濤, 陳錦年, 王凡. 中部型El Nio與北極海冰變化的聯(lián)系[J]. 海洋湖沼通報, 2015(3):1-13.

    Zuo Tao, Chen Jinnian, Wang Fan. The contribution of Central Pacific El Nio(La Nia) to the Arctic Sea Ice variation[J]. Transactions of Oceanology and Limnology, 2015(3):1-13.

    [33] Hu Chundi, Yang Song, Wu Qigang, et al. Shifting El Nio inhibits summer Arctic warming and Arctic sea-ice melting over the Canada Basin[J]. Nature Communications, 2016, 11721, doi:10.1038/ncomms11721.

    [34] 陳迪, 高山紅, 陳錦年. 印太暖池區(qū)域海溫異常與北極海冰變化的聯(lián)系[J]. 極地研究, 2016, 28(1):49-57.

    Chen Di, Gao Shanhong, Chen Jinnian. Impact of the Indo-Pacific warm pool sst anomaly on Arctic sea ice variation[J]. Chinese Journal of Polar Research, 2016, 28(1):49-57.

    [35] Kug J S, Jin F F, An S I. Two-types of El Nio events: Cold tongue El Nio and warm pool El Nio[J]. Journal of Climate, 2009, 22(6):1499-1515.

    [36] Ren H L, Jin F F. Nio indices for two types of ENSO[J]. Geophysical Research Letters, 2011, 38(4):L04704.

    [37] 王磊, 張文君, 祁莉,等. 兩類La Nia季節(jié)演變過程的海氣耦合特征對比[J]. 海洋學報, 2014, 36(1):72-85.

    Wang Lei, Zhang Wenjun, Qi Li, et al. Contrasting air-sea features associated with two types of La Nia during the seasonal evolution[J]. Haiyang Xuebao, 2014, 36(1):72-85.

    [38] 何珊珊, 張文君, 祁莉,等. 兩類厄爾尼諾事件發(fā)展年秋季印度洋海溫異常特征對比[J]. 氣象學報, 2015, 73(3):515-528.

    He Shanshan, Zhang Wenjun, Qi Li, et al. Contrasting SST anomalies over the Indian Ocean between the two types of El Nio events during boreal autumn[J]. Acta Meteorologica Sinica, 2015, 73(3):515-528.

    [39] Zhang Wenjun, Wang Lei, Xiang Baoqiang, et al. Impacts of two types of La Nia on the NAO during boreal winter[J]. Climate Dynamics, 2014, 44(5):1351-1366.

    [40] Rayner N A, Parker D E, Horton E B, et al. Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century[J]. Journal of Geophysical Research Atmospheres, 2003, 108(D14):1063-1082.

    [41] Kalnay E, Kanamitsu M, Kistler R, et al. The NCEP/NCAR 40-year reanalysis project[J]. Bulletin of the American Meteorological Society, 1995, 77(3):437-472.

    [42] Serreze M C, Holland M M, Stroeve J. Perspectives on the Arctic’s shrinking sea-ice cover[J]. Science, 2007, 315(5818):1533-1536.

    [43] Yeh S W, Kug J S, Dewitte B, et al. El Nio in a changing climate[J]. Nature, 2009, 461(7263):511-514.

    [44] 朱大勇, 趙進平, 史久新. 北極楚科奇海海冰面積多年變化的研究[J]. 海洋學報, 2007, 29(2):25-33.

    Zhu Dayong, Zhao Jinping, Shi Jiuxin. Study on the multi-year variations of sea ice cover of Chukchi Sea in Arctic Ocean[J]. Haiyang Xuebao, 2007, 29(2):25-33.

    [45] 武炳義, 黃榮輝, 高登義. 與北大西洋接壤的北極海冰和年際氣候變化[J]. 科學通報, 2000, 45(18):1993-1998.

    Wu Bingyi, Huang Ronghui, Gao Dengyi. Arctic sea ice and interannual climate change bordering north Atlantic Ocean[J]. Chinese Science Bulletin, 2000, 45(18):1993-1998.

    [46] Walsh J E, Johnson C M. Interannual atmospheric variability and associated fluctuations in Arctic Sea ice extent[J]. Journal of Geophysical Research Atmospheres, 1979, 84(C11): 6915-6928.

    [47] Maslanik J A, Fowler C, Stroeve J, et al. A younger, thinner Arctic ice cover: Increased potential for rapid, extensive sea-ice loss[J]. Geophysical Research Letters, 2007, 34(24):497-507.

    [48] Hudson S R, Granskog M A, Arild S, et al. Energy budget of first-year Arctic sea ice in advanced stages of melt[J]. Geophysical Research Letters, 2013, 40(11):2679-2683.

    [49] 程彥杰, 卞林根, 陸龍驊. 南極海冰濤動與ENSO的關系[J]. 應用氣象學報, 2002, 13(6):711-717.

    Cheng Yanjie, Bian Lingen, Lu Longhua. The relationship between Antarctic sea ice oscillation and ENSO[J]. Journal of Applied Meteorological Science, 2002, 13(6):711-717.

    [50] 何金海, 武豐民, 祁莉. 秋季北極海冰與歐亞冬季氣溫在年代際和年際尺度上的不同聯(lián)系[J]. 地球物理學報, 2015(4):1089-1102.

    He Jinhai, Wu Fengmin, Qi Li, et al. 2015. Decadal/interdecadal linking between autumn Arctic sea ice and following winter Eurasian air temperature[J]. Chinese Journal of Geophysics, 2015(4):1089-1102.

    [51] 武豐民, 何金海, 祁莉,等. 海冰消融背景下北極增溫的季節(jié)差異及其原因探討[J]. 海洋學報, 2014, 36(3):39-47.

    Wu Fengmin, He Jinhai, Qi Li,et al. The seasonal difference of Arctic warming and it’s mechanism under sea ice cover diminishing[J]. Haiyang Xuebao, 2014, 36(3):39-47.

    [52] 劉喜迎, 劉海龍. 大氣變率對北極地區(qū)近期海冰變化趨勢影響數(shù)值模擬研究[J]. 地球物理學報, 2012, 55(9):2867-2875.

    Liu Xiying, Liu Hailong. Investigation of influence of atmospheric variability on sea variation trend in recent years in the Arctic with numerical sea ice-ocean coupled model[J]. Chinese Journal of Geophysics, 2012, 55(9):2867-2875.

    [53] 鄧娟. 北半球海冰變化及其與氣候要素的關系[D]. 南京:南京大學, 2014.

    Deng Juan. Northern Hemisphere sea ice variability and its relationship with climate factors[D]. Nanjing: Nanjing University, 2014.

    The possible influence of sea surface temperature anomalies over the tropical Pacific on the Arctic sea ice

    Luo Rui1, Qi Li1, Zhang Wenjun1, He Jinhai1

    (1.KeyLaboratoryofMeteorologicalDisaster,MinistryofEducation/InternationalJointLaboratoryonClimateandEnvironmentChange(ILCEC)/CollaborativeInnovationCenteronForecastandEvaluationofMeteorologicalDisasters,NanjingUniversityofInformationScience&Technology,Nanjing210044,China)

    The possible influence of sea surface temperature (SST) anomalies over the tropical Pacific on the Arctic sea ice and their possible mechanisms are researched from atmospheric circulation and net surface heat flux by using monthly data of HadISST and NCEP/NCAR spanning from 1950 to 2015. The main results obtained show that during summer and autumn of ENSO development phase, EP and CP El Nio events have no apparent connection with the Arctic sea ice anomalies. However, strong signals are detected when it happens to La Nia and there are distinctive differences between EP and CP events.Sea ice reduce in Barents/Kara Sea in EP La Nia but increase in East Siberia/Chukchi Sea in CP La Nia. During summer and autumn of EP La Nia’s development, SST anomalies over the tropical Pacific leads to negative SLP in Barents/Kara Sea by teleconnection wave trains and forms AO-like pattern together with positive SLP in middle latitude. It changes the wind which brings warm advection and southerly anomalies helps to transfer more warm water from North Atlantic, high water vapor content at the same time enhance the income of net surface heat flux and result in the decrease of sea ice in Barents/Kara Sea. In the situation of CP La Nia, East Siberia/Chukchi Sea are dominated by north wind in summer because of a cyclonic circulation in the eastern side, sea ice are transported here from North Pole and thus sea ice concentration increased, effect of net surface heat flux is small.

    Arctic sea ice; EP/CP La Nia; atmospheric circulation; net surface heat flux

    2016-08-30;

    2016-12-26。

    國家自然科學基金資助項目(41475086);長江學者和創(chuàng)新團隊發(fā)展計劃資助(PCSIRT);江蘇高校優(yōu)勢學科建設工程資助項目(PAPD)。

    羅蕊(1991-),女,云南省文山壯族苗族自治州人,主要從事海-冰-氣相互作用研究。E-mail:lorry_17723@163.com

    *通信作者:祁莉,女,教授,主要從事季風和海-陸-氣相互作用研究。E-mail:qili@nuist.edu.cn

    10.3969/j.issn.0253-4193.2017.07.006

    P732.6

    A

    0253-4193(2017)07-0053-17

    羅蕊,祁莉,張文君,等. 熱帶太平洋海溫異常對北極海冰的可能影響[J]. 海洋學報, 2017, 39(7): 53-69,

    Luo Rui, Qi Li, Zhang Wenjun, et al. The possible influence of sea surface temperature anomalies over the tropical Pacific on the Arctic sea ice[J]. Haiyang Xuebao, 2017, 39(7): 53-69, doi:10.3969/j.issn.0253-4193.2017.07.006

    猜你喜歡
    北極海海溫海冰
    基于深度學習的海溫觀測數(shù)據(jù)質量控制應用研究
    海洋通報(2021年3期)2021-08-14 02:20:48
    末次盛冰期以來巴倫支海-喀拉海古海洋環(huán)境及海冰研究進展
    海洋通報(2021年3期)2021-08-14 02:20:38
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    海洋通報(2021年1期)2021-07-23 01:55:22
    南印度洋偶極型海溫與中國西南地區(qū)初秋降水的關系
    基于SIFT-SVM的北冰洋海冰識別研究
    南、北極海冰的長期變化趨勢及其與大氣環(huán)流的聯(lián)系
    中部型El Nino與北極海冰變化的聯(lián)系
    累積海冰密集度及其在認識北極海冰快速變化的作用
    波浪輸運通量剩余量對上層海溫的影響分析
    應用MODIS數(shù)據(jù)監(jiān)測河北省近海海域海冰
    河北遙感(2014年4期)2014-07-10 13:54:59
    人妻人人澡人人爽人人| 国产高清videossex| 男女高潮啪啪啪动态图| 久热爱精品视频在线9| 嫩草影视91久久| tube8黄色片| 亚洲成人免费av在线播放| 国产一区有黄有色的免费视频| 十八禁网站免费在线| 欧美日韩亚洲综合一区二区三区_| 亚洲五月婷婷丁香| 正在播放国产对白刺激| 亚洲国产精品成人久久小说| 18禁观看日本| 99国产精品免费福利视频| 亚洲欧美一区二区三区黑人| 亚洲国产av新网站| 中国国产av一级| 亚洲自偷自拍图片 自拍| 老司机靠b影院| 在线观看免费视频网站a站| 视频区欧美日本亚洲| 亚洲av国产av综合av卡| 久久久国产成人免费| 亚洲国产中文字幕在线视频| 一边摸一边抽搐一进一出视频| 涩涩av久久男人的天堂| 免费高清在线观看日韩| 亚洲伊人久久精品综合| 天天躁日日躁夜夜躁夜夜| 亚洲专区中文字幕在线| 欧美日本中文国产一区发布| 久久香蕉激情| 窝窝影院91人妻| 超碰97精品在线观看| 亚洲视频免费观看视频| a在线观看视频网站| 久久久久视频综合| 亚洲国产欧美网| 免费在线观看影片大全网站| 成人免费观看视频高清| 免费在线观看黄色视频的| 高潮久久久久久久久久久不卡| 97人妻天天添夜夜摸| 啦啦啦免费观看视频1| 国产男女超爽视频在线观看| 国产精品亚洲av一区麻豆| 黑人操中国人逼视频| 久久久久久人人人人人| 人人妻人人澡人人爽人人夜夜| 黄网站色视频无遮挡免费观看| 久久精品国产亚洲av香蕉五月 | 岛国毛片在线播放| 在线观看免费午夜福利视频| 大片免费播放器 马上看| 老汉色∧v一级毛片| 国产精品久久久人人做人人爽| 人成视频在线观看免费观看| 黑人猛操日本美女一级片| 成年人黄色毛片网站| 青春草亚洲视频在线观看| 欧美精品亚洲一区二区| 考比视频在线观看| 一区二区三区精品91| 制服人妻中文乱码| 12—13女人毛片做爰片一| 精品乱码久久久久久99久播| 日本一区二区免费在线视频| av电影中文网址| 日本91视频免费播放| 美女脱内裤让男人舔精品视频| 国产极品粉嫩免费观看在线| 丰满少妇做爰视频| 啦啦啦免费观看视频1| 69av精品久久久久久 | 亚洲欧美清纯卡通| 在线观看免费视频网站a站| a级毛片在线看网站| 免费久久久久久久精品成人欧美视频| 国产视频一区二区在线看| 精品视频人人做人人爽| 一区二区三区精品91| 91老司机精品| 亚洲熟女毛片儿| 国产不卡av网站在线观看| 又大又爽又粗| 亚洲 欧美一区二区三区| 12—13女人毛片做爰片一| 性高湖久久久久久久久免费观看| 国产成人免费观看mmmm| 成年人午夜在线观看视频| 美女中出高潮动态图| 高清黄色对白视频在线免费看| 国产精品免费大片| 在线观看www视频免费| 国产精品国产av在线观看| www.av在线官网国产| 大型av网站在线播放| 日韩一卡2卡3卡4卡2021年| 午夜日韩欧美国产| 欧美在线黄色| 免费高清在线观看视频在线观看| 窝窝影院91人妻| 亚洲伊人色综图| 黄色视频,在线免费观看| 国产区一区二久久| 老司机深夜福利视频在线观看 | 精品一区在线观看国产| 免费在线观看影片大全网站| 黄色视频,在线免费观看| 可以免费在线观看a视频的电影网站| 男人爽女人下面视频在线观看| 一区福利在线观看| 桃红色精品国产亚洲av| 亚洲 欧美一区二区三区| www.自偷自拍.com| 久久99一区二区三区| 午夜成年电影在线免费观看| 啦啦啦啦在线视频资源| 天天添夜夜摸| 久久久精品94久久精品| 久久久久久人人人人人| 在线看a的网站| 国产男女内射视频| 不卡一级毛片| 久久精品人人爽人人爽视色| av网站在线播放免费| 欧美激情久久久久久爽电影 | 亚洲精华国产精华精| 精品一区二区三区av网在线观看 | avwww免费| 久久人妻熟女aⅴ| 多毛熟女@视频| 亚洲三区欧美一区| 国产一区二区三区综合在线观看| 在线av久久热| 曰老女人黄片| 国产又色又爽无遮挡免| 如日韩欧美国产精品一区二区三区| 看免费av毛片| 欧美在线黄色| 国产熟女午夜一区二区三区| 国产亚洲精品一区二区www | 国产日韩欧美亚洲二区| 视频在线观看一区二区三区| 亚洲精品av麻豆狂野| 国产在线一区二区三区精| 在线看a的网站| 国产日韩一区二区三区精品不卡| 成人av一区二区三区在线看 | 欧美日韩精品网址| 视频区欧美日本亚洲| 国产男人的电影天堂91| 久久人人爽人人片av| 精品少妇黑人巨大在线播放| 操出白浆在线播放| 曰老女人黄片| 国产成人精品在线电影| 久久午夜综合久久蜜桃| 日本91视频免费播放| 亚洲激情五月婷婷啪啪| 国产成人精品久久二区二区免费| 欧美黑人精品巨大| 操出白浆在线播放| 久久国产精品影院| 成人免费观看视频高清| 中文字幕最新亚洲高清| 日韩欧美国产一区二区入口| 伊人亚洲综合成人网| 久久久久精品国产欧美久久久 | 亚洲激情五月婷婷啪啪| 91九色精品人成在线观看| 亚洲第一av免费看| 中文欧美无线码| 男男h啪啪无遮挡| 亚洲欧美色中文字幕在线| av超薄肉色丝袜交足视频| av国产精品久久久久影院| 精品亚洲成a人片在线观看| 亚洲黑人精品在线| 一级毛片精品| 日韩电影二区| 亚洲欧美色中文字幕在线| 最新的欧美精品一区二区| 啦啦啦在线免费观看视频4| 国产亚洲av片在线观看秒播厂| 久久精品成人免费网站| 国产av国产精品国产| 日本91视频免费播放| 国产黄色免费在线视频| 大片免费播放器 马上看| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 99国产精品一区二区三区| 另类亚洲欧美激情| 最新在线观看一区二区三区| 蜜桃国产av成人99| 69av精品久久久久久 | 午夜精品国产一区二区电影| 电影成人av| av国产精品久久久久影院| 日本欧美视频一区| 别揉我奶头~嗯~啊~动态视频 | 超色免费av| 丝袜美足系列| 人妻一区二区av| 国产区一区二久久| 久久久久久久大尺度免费视频| 免费女性裸体啪啪无遮挡网站| 男女国产视频网站| 无遮挡黄片免费观看| 丝袜人妻中文字幕| 美女主播在线视频| 日本猛色少妇xxxxx猛交久久| 国产精品一二三区在线看| 国产一区二区三区在线臀色熟女 | 国产欧美日韩精品亚洲av| 亚洲免费av在线视频| 下体分泌物呈黄色| 一个人免费在线观看的高清视频 | 国产亚洲精品第一综合不卡| 日本av免费视频播放| 一区二区三区乱码不卡18| 国产免费av片在线观看野外av| 欧美日韩一级在线毛片| 99国产精品一区二区三区| 精品人妻熟女毛片av久久网站| 国产成人欧美在线观看 | 色婷婷久久久亚洲欧美| av有码第一页| 我的亚洲天堂| 涩涩av久久男人的天堂| 一进一出抽搐动态| 亚洲三区欧美一区| 久久免费观看电影| 在线观看舔阴道视频| 1024香蕉在线观看| 女人精品久久久久毛片| 妹子高潮喷水视频| 伊人亚洲综合成人网| 国产成人精品久久二区二区免费| 国产欧美日韩精品亚洲av| 欧美日韩亚洲国产一区二区在线观看 | 天天影视国产精品| 日本av免费视频播放| 亚洲一区中文字幕在线| 免费观看人在逋| 两个人看的免费小视频| 美女中出高潮动态图| 午夜老司机福利片| 九色亚洲精品在线播放| 天堂中文最新版在线下载| 精品少妇久久久久久888优播| 国产日韩欧美亚洲二区| 午夜老司机福利片| 在线亚洲精品国产二区图片欧美| 亚洲人成电影观看| 美女国产高潮福利片在线看| 俄罗斯特黄特色一大片| 久久青草综合色| 80岁老熟妇乱子伦牲交| 2018国产大陆天天弄谢| www.自偷自拍.com| 操出白浆在线播放| 欧美 日韩 精品 国产| 亚洲激情五月婷婷啪啪| 欧美日韩精品网址| 满18在线观看网站| 99国产精品99久久久久| 黄片小视频在线播放| 一本—道久久a久久精品蜜桃钙片| 在线观看www视频免费| 人成视频在线观看免费观看| 在线十欧美十亚洲十日本专区| 欧美黄色片欧美黄色片| 久久午夜综合久久蜜桃| 热99久久久久精品小说推荐| 午夜成年电影在线免费观看| 精品视频人人做人人爽| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av成人不卡在线观看播放网 | 国产日韩欧美亚洲二区| 国产熟女午夜一区二区三区| 不卡一级毛片| 如日韩欧美国产精品一区二区三区| 欧美精品亚洲一区二区| 久久精品国产亚洲av高清一级| 91老司机精品| 国产免费av片在线观看野外av| 成年女人毛片免费观看观看9 | 王馨瑶露胸无遮挡在线观看| 欧美激情 高清一区二区三区| 精品久久蜜臀av无| 91成人精品电影| 国产熟女午夜一区二区三区| 国产成人欧美| 大香蕉久久成人网| 少妇 在线观看| 性色av乱码一区二区三区2| 国精品久久久久久国模美| 日韩欧美一区二区三区在线观看 | 亚洲av欧美aⅴ国产| 午夜91福利影院| av视频免费观看在线观看| 日韩一卡2卡3卡4卡2021年| 美女视频免费永久观看网站| 久久精品亚洲av国产电影网| 丰满饥渴人妻一区二区三| 日本a在线网址| 91九色精品人成在线观看| www.av在线官网国产| 黄片播放在线免费| 久久亚洲国产成人精品v| svipshipincom国产片| 午夜福利视频在线观看免费| 女警被强在线播放| 天天躁夜夜躁狠狠躁躁| 欧美人与性动交α欧美精品济南到| 亚洲专区国产一区二区| 视频在线观看一区二区三区| 99热全是精品| 黄片大片在线免费观看| 国产精品一二三区在线看| 免费看十八禁软件| 国产精品国产三级国产专区5o| 最近最新中文字幕大全免费视频| 王馨瑶露胸无遮挡在线观看| 久久国产精品大桥未久av| 国产成人免费无遮挡视频| 美女高潮喷水抽搐中文字幕| 精品一区二区三区av网在线观看 | av欧美777| 精品久久蜜臀av无| 日韩 欧美 亚洲 中文字幕| 免费不卡黄色视频| 一个人免费看片子| 一本—道久久a久久精品蜜桃钙片| 99久久人妻综合| 18禁黄网站禁片午夜丰满| 婷婷色av中文字幕| 亚洲色图 男人天堂 中文字幕| 日韩中文字幕欧美一区二区| 成年人黄色毛片网站| 韩国高清视频一区二区三区| 免费看十八禁软件| 国产亚洲一区二区精品| 亚洲五月婷婷丁香| 美女视频免费永久观看网站| 最黄视频免费看| 国产亚洲一区二区精品| av超薄肉色丝袜交足视频| 一边摸一边抽搐一进一出视频| 国产免费福利视频在线观看| 国产精品久久久久久精品电影小说| 满18在线观看网站| 777久久人妻少妇嫩草av网站| 纵有疾风起免费观看全集完整版| av免费在线观看网站| 国产精品久久久久久人妻精品电影 | 丝袜脚勾引网站| 国产免费福利视频在线观看| 丝袜脚勾引网站| 啪啪无遮挡十八禁网站| 国产精品熟女久久久久浪| 精品免费久久久久久久清纯 | 日日摸夜夜添夜夜添小说| 美女高潮喷水抽搐中文字幕| 国产99久久九九免费精品| 日本精品一区二区三区蜜桃| 午夜免费鲁丝| 午夜福利,免费看| 91老司机精品| 国产一区二区激情短视频 | 亚洲一码二码三码区别大吗| 日韩 亚洲 欧美在线| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜精品一区,二区,三区| 少妇裸体淫交视频免费看高清 | 欧美日韩中文字幕国产精品一区二区三区 | 午夜成年电影在线免费观看| 久久天堂一区二区三区四区| av一本久久久久| 日韩 欧美 亚洲 中文字幕| 久久久国产一区二区| 巨乳人妻的诱惑在线观看| 一本色道久久久久久精品综合| 欧美精品高潮呻吟av久久| 成人国语在线视频| av欧美777| 久久久精品94久久精品| 人妻一区二区av| 久久99热这里只频精品6学生| 久久亚洲国产成人精品v| 久久久久久久久免费视频了| 国产欧美日韩精品亚洲av| 国精品久久久久久国模美| 视频区欧美日本亚洲| 亚洲色图综合在线观看| 日韩视频一区二区在线观看| 老司机亚洲免费影院| 亚洲国产精品一区二区三区在线| 一级a爱视频在线免费观看| 一级片免费观看大全| 亚洲精品成人av观看孕妇| 91国产中文字幕| 动漫黄色视频在线观看| 一进一出抽搐动态| 亚洲少妇的诱惑av| 大片电影免费在线观看免费| 午夜福利一区二区在线看| 中文精品一卡2卡3卡4更新| 欧美精品高潮呻吟av久久| 久久久国产一区二区| 欧美变态另类bdsm刘玥| 国产激情久久老熟女| 日韩视频在线欧美| 黑人猛操日本美女一级片| 免费在线观看完整版高清| 免费在线观看日本一区| 免费高清在线观看日韩| 欧美激情久久久久久爽电影 | 久久影院123| 午夜福利,免费看| 亚洲精品av麻豆狂野| 久久国产精品大桥未久av| 超碰成人久久| 久久久久国内视频| 婷婷成人精品国产| 热99re8久久精品国产| 亚洲成人免费av在线播放| av网站在线播放免费| 国产精品国产av在线观看| 淫妇啪啪啪对白视频 | av一本久久久久| 欧美黑人欧美精品刺激| 亚洲国产欧美网| 亚洲欧美一区二区三区黑人| 一边摸一边做爽爽视频免费| 国产欧美日韩一区二区精品| 日韩制服丝袜自拍偷拍| 久久99热这里只频精品6学生| 国产三级黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 一级毛片精品| 免费观看人在逋| 老司机在亚洲福利影院| 久久久久久久久久久久大奶| 欧美亚洲日本最大视频资源| 两人在一起打扑克的视频| 精品一区二区三卡| 国产成人啪精品午夜网站| 久久av网站| 国产成人a∨麻豆精品| 成人免费观看视频高清| 老司机影院成人| 亚洲av欧美aⅴ国产| 久久久精品区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 日本欧美视频一区| 国产一区二区 视频在线| 亚洲中文av在线| 亚洲欧美一区二区三区久久| 老鸭窝网址在线观看| 国产亚洲av高清不卡| av网站免费在线观看视频| 啦啦啦啦在线视频资源| 99久久精品国产亚洲精品| 搡老岳熟女国产| 中文字幕最新亚洲高清| 色94色欧美一区二区| 女性被躁到高潮视频| 电影成人av| 国产一级毛片在线| 在线av久久热| 亚洲成人国产一区在线观看| 亚洲 欧美一区二区三区| 精品国产乱码久久久久久男人| 国产精品99久久99久久久不卡| 性色av乱码一区二区三区2| 又紧又爽又黄一区二区| 超色免费av| 国产精品久久久久久精品古装| 人成视频在线观看免费观看| 午夜激情av网站| av有码第一页| 丝袜美足系列| 母亲3免费完整高清在线观看| 婷婷色av中文字幕| 日韩大片免费观看网站| a级片在线免费高清观看视频| 久久久久精品国产欧美久久久 | 亚洲欧美一区二区三区黑人| 成年人午夜在线观看视频| 免费观看av网站的网址| 在线观看舔阴道视频| 久久人妻熟女aⅴ| 老司机深夜福利视频在线观看 | 久久精品人人爽人人爽视色| 久久这里只有精品19| 久久久久久人人人人人| 人妻 亚洲 视频| 又紧又爽又黄一区二区| 久久精品人人爽人人爽视色| 99久久精品国产亚洲精品| 啦啦啦啦在线视频资源| 国产精品一区二区在线观看99| 亚洲精品国产色婷婷电影| 欧美黑人欧美精品刺激| 少妇猛男粗大的猛烈进出视频| 大陆偷拍与自拍| 国产1区2区3区精品| 国产精品.久久久| 亚洲 欧美一区二区三区| 女性生殖器流出的白浆| 中文字幕av电影在线播放| 国产三级黄色录像| www.精华液| 久久精品aⅴ一区二区三区四区| 日本猛色少妇xxxxx猛交久久| 日韩制服丝袜自拍偷拍| 欧美日韩福利视频一区二区| 妹子高潮喷水视频| 国产成人av教育| 精品人妻1区二区| 波多野结衣一区麻豆| 国产一级毛片在线| 精品国产乱码久久久久久小说| 99久久人妻综合| 日韩一卡2卡3卡4卡2021年| 一个人免费看片子| 91麻豆av在线| 久久毛片免费看一区二区三区| 亚洲综合色网址| 亚洲av成人不卡在线观看播放网 | 亚洲 欧美一区二区三区| 亚洲成av片中文字幕在线观看| 国产高清videossex| 91麻豆精品激情在线观看国产 | 精品国产乱子伦一区二区三区 | 日日爽夜夜爽网站| 亚洲av电影在线进入| xxxhd国产人妻xxx| 18禁国产床啪视频网站| 精品久久蜜臀av无| 亚洲视频免费观看视频| 国产1区2区3区精品| 人人妻人人澡人人爽人人夜夜| 亚洲欧美日韩高清在线视频 | 人人妻,人人澡人人爽秒播| 别揉我奶头~嗯~啊~动态视频 | 欧美亚洲 丝袜 人妻 在线| 欧美成人午夜精品| 国产精品成人在线| 纯流量卡能插随身wifi吗| 精品国产一区二区久久| 乱人伦中国视频| 大码成人一级视频| 国产精品香港三级国产av潘金莲| 久久ye,这里只有精品| 少妇被粗大的猛进出69影院| 国产精品一区二区免费欧美 | 啦啦啦在线免费观看视频4| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 国内毛片毛片毛片毛片毛片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品麻豆人妻色哟哟久久| 成人三级做爰电影| a级片在线免费高清观看视频| 国产高清国产精品国产三级| 成年人免费黄色播放视频| 777久久人妻少妇嫩草av网站| 国产麻豆69| 午夜福利,免费看| 欧美精品啪啪一区二区三区 | 久久久久久免费高清国产稀缺| 欧美一级毛片孕妇| 欧美日本中文国产一区发布| a级毛片黄视频| 精品国产一区二区久久| 黄片小视频在线播放| 欧美性长视频在线观看| 国产熟女午夜一区二区三区| 80岁老熟妇乱子伦牲交| 国产在视频线精品| 制服诱惑二区| 美女中出高潮动态图| 热re99久久国产66热| 一区二区三区精品91| 欧美在线黄色| 亚洲av片天天在线观看| 青草久久国产| 久久精品熟女亚洲av麻豆精品| 真人做人爱边吃奶动态| 超色免费av| 男人添女人高潮全过程视频| 大片免费播放器 马上看| 国产黄频视频在线观看| 午夜影院在线不卡| 国产成人a∨麻豆精品| 91国产中文字幕| 一区福利在线观看| 99久久综合免费| 黄片大片在线免费观看| 男人添女人高潮全过程视频| 99久久99久久久精品蜜桃| 欧美成狂野欧美在线观看| h视频一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 91国产中文字幕| 亚洲精品久久成人aⅴ小说| 王馨瑶露胸无遮挡在线观看| 后天国语完整版免费观看| 国产xxxxx性猛交| 9191精品国产免费久久|