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

    北歐海主要海盆海面熱通量的多年變化?

    2014-06-24 14:04:19趙進(jìn)平KenDrinkwater
    關(guān)鍵詞:海盆格陵蘭長(zhǎng)波

    趙進(jìn)平,Ken Drinkwater

    (1.中國(guó)海洋大學(xué),山東青島266100;2.挪威海洋研究所Bergen,N-5024)

    北歐海主要海盆海面熱通量的多年變化?

    趙進(jìn)平1,Ken Drinkwater2

    (1.中國(guó)海洋大學(xué),山東青島266100;2.挪威海洋研究所Bergen,N-5024)

    北歐海有暖流和寒流注入,又發(fā)生大量回流,水團(tuán)特性異常復(fù)雜。由于北歐海的環(huán)流受地形控制,其水團(tuán)的分布與4個(gè)海盆的分布有密切的關(guān)系。本文研究各個(gè)海盆熱通量變化的差異,以研究獲取對(duì)北歐海海氣相互作用區(qū)域差異的認(rèn)識(shí)。北歐海的熱量夏季以太陽短波輻射為主,冬季以來自海洋的長(zhǎng)波輻射、感熱和潛熱通量為主。海盆間的差異主要體現(xiàn)在格陵蘭海,其變化幅度短波輻射高出50%,長(zhǎng)波輻射高出大約40%,潛熱高出大約60%,感熱高出近4倍??赡艿脑蚴?,格陵蘭海強(qiáng)烈的感熱和潛熱釋放導(dǎo)致海溫降低,氣溫也受北極冷空氣的影響,形成與暖流區(qū)迥異的自然環(huán)境。過去30年發(fā)生了2次顯著的熱量減少事件。其中,1987年的事件很可能與冰島的火山噴發(fā)有關(guān),火山噴發(fā)對(duì)短波輻射的影響長(zhǎng)達(dá)一年之久,導(dǎo)致感熱和潛熱也同步減少。1998年格陵蘭海的潛熱和感熱明顯減少,與北極海冰輸出導(dǎo)致的海溫偏低有顯著關(guān)系。文章分析了4個(gè)主要海盆熱通量的變化與北極濤動(dòng)(AO)指數(shù)的關(guān)系。結(jié)果表明,發(fā)生在冰島海的向下短波輻射和發(fā)生在格陵蘭海的感熱和潛熱與AO相關(guān)度較高,體現(xiàn)了與AO的密切關(guān)系。這些熱通量與AO指數(shù)的滑動(dòng)相關(guān)系數(shù)表明,1992年以前冰島海的短波輻射與AO的相關(guān)性非常高,而格陵蘭海的感熱和潛熱在1993年后與AO高度相關(guān),是值得深入研究的現(xiàn)象。本文的結(jié)果支持以下觀點(diǎn):北歐海對(duì)北極濤動(dòng)的貢獻(xiàn)主要是格陵蘭海的感熱和潛熱釋放通過冰島低壓區(qū)的上升氣流影響冰島低壓的云量,從而影響到達(dá)的太陽輻射而導(dǎo)致大氣環(huán)流的變化。

    北歐海;熱通量;海氣耦合;北極濤動(dòng);北極濤動(dòng)核心區(qū)

    1 導(dǎo)論

    北歐海(Nordic Sea)是連接北大西洋與北冰洋的海域,由格陵蘭海、冰島海和挪威海組成,也稱GIN海(見圖1)。來自北大西洋的暖流——挪威流沿北歐海東部進(jìn)入,從巴倫支海和Fram海峽流出北歐海。另一支來自北冰洋的寒流——東格陵蘭流沿北歐海西部進(jìn)入,并通過Denmark海峽流出北歐海。受海底地貌的影響和地形的約束,暖流多次發(fā)生從主流分離的現(xiàn)象,稱為回流,在北歐海中部形成大范圍回流區(qū)[1]。回流區(qū)的水體與周邊水體交換,形成特有的水團(tuán),稱為大西洋回流水(Recirculating Atlantic Water)[2]。由于回流區(qū)的水體要比暖流區(qū)和寒流區(qū)水體的滯留時(shí)間更長(zhǎng),有更多的時(shí)間發(fā)生海氣相互作用,成為北歐海海氣熱交換的主體。在相對(duì)較冷的冬季,海洋向大氣大量釋放熱量,導(dǎo)致了深對(duì)流和高密度水的形成[3]。

    據(jù)此,北歐海的上層分為3個(gè)主要區(qū)域:在西部受到極地冷水直接影響的極區(qū),在東部受到大西洋暖水直接影響的大西洋區(qū),以及在北歐海中部占據(jù)了格陵蘭海的中部及冰島海的中部和東北部的過渡區(qū)域[4-5]。各個(gè)區(qū)域之間由海洋鋒來分隔[6-8]。

    北歐海發(fā)生2種海洋對(duì)流,形成巨大的下沉流量。最顯著的對(duì)流是表面冷卻和鹵水釋放導(dǎo)致的對(duì)流,即開放海洋對(duì)流(Open-Ocean Convection)[9]。開放海洋對(duì)流過程引發(fā)海洋中垂直熱交換,加強(qiáng)海洋中的熱量釋放,產(chǎn)生強(qiáng)海氣相互作用,同時(shí)影響海洋和大氣的運(yùn)動(dòng)[10-11]。此外,在格陵蘭海的極鋒附近,發(fā)生大范圍的等密度對(duì)流,主要是回流水與來自北冰洋水體的混合增密導(dǎo)致的對(duì)流[12-13]。

    由于環(huán)流和回流都與海底地形有關(guān),北歐海的水團(tuán)與各個(gè)主要海盆緊密聯(lián)系,北歐海海氣界面熱通量的差異主要是各個(gè)海盆之間熱通量的差異。因此,本文研究各個(gè)主要海盆海面熱通量的差異。在北歐海存在4個(gè)大的海盆,分別為格陵蘭海、冰島海、挪威海盆、羅弗敦海盆。其中,冰島海的深度較淺,大約1000m,其余3個(gè)都是深海盆,深度在2800m以上。由于水體來源不同,海水的溫度不同,各個(gè)海盆的海氣交換有很大的差別。其中,挪威海盆和羅弗敦海盆主要是挪威暖流的暖水回流區(qū),水溫較高,是主要的感熱釋放區(qū)。而格陵蘭海盆是東格陵蘭流的回流水和暖流的回流水匯合區(qū),水溫較低,感熱和潛熱的釋放都很強(qiáng)烈[14]。各個(gè)海盆的海氣交換對(duì)北歐海的大氣環(huán)流有重要的貢獻(xiàn),其中,挪威海盆和羅弗敦海盆的暖水回流在不間斷的循環(huán)過程中不斷向大氣釋放熱量,是歐洲熱量的主要來源。

    北歐海是北大西洋濤動(dòng)/北極濤動(dòng)(NAO/AO)的關(guān)鍵海域[15-16],在氣候系統(tǒng)中起到重要作用。NAO明顯與系統(tǒng)性的風(fēng)場(chǎng)距平變化、潛熱和感熱通量的變化以及海面溫度場(chǎng)變化有關(guān)[17-18],亞極區(qū)渦旋的海表溫度、海冰范圍與氣壓十年尺度變化幾乎是同步的[19],來自格陵蘭島的下降風(fēng)減弱有助于抑制蒸發(fā),改變了海冰的分布,導(dǎo)致海冰與氣壓場(chǎng)變化呈正相關(guān)關(guān)系[20]。文獻(xiàn)[21]通過研究北極濤動(dòng)(AO)的變化與北半球海面氣壓場(chǎng)(SLP)變化的關(guān)系,指出存在1個(gè)以北歐海為主體的北極濤動(dòng)核心區(qū)(AO核心區(qū)),AO核心區(qū)的平均SLP與北極濤動(dòng)指數(shù)(AOI)高度負(fù)相關(guān),相關(guān)系數(shù)達(dá)到0.945,二者幾乎可以相互替代。然而,北歐海東西兩側(cè)最大溫差達(dá)到10℃以上,暖流區(qū)域寒流區(qū)的海氣通量也有很大的差別。文獻(xiàn)[21]的結(jié)果提出1個(gè)重要的命題,是什么因素導(dǎo)致AO核心區(qū)的氣壓與AOI一致變化?空間不均勻的海洋如何導(dǎo)致大氣參數(shù)的均勻特性?為此,需要研究北歐海表面熱通量如何影響海洋和大氣的變化。

    迄今為止,北歐海內(nèi)部的長(zhǎng)期連續(xù)數(shù)據(jù)主要有:各國(guó)沿岸氣象站,位于揚(yáng)馬延島的氣象站,以及位于挪威南部近海的長(zhǎng)期監(jiān)測(cè)船(M站,1949—2009)。然而,由于北歐海各項(xiàng)參數(shù)的空間變化范圍大,這樣幾個(gè)站的數(shù)據(jù)不足以了解北歐海海氣相互作用的全貌。

    本文的目的是了解各個(gè)熱通量的時(shí)間變化及其與AO的相關(guān)性,采用NCEP-DOE Reanalysis-2在北歐海的格點(diǎn)化數(shù)據(jù),研究北歐海海氣通量與北極濤動(dòng)指數(shù)的關(guān)系,揭示北歐海海洋對(duì)北極濤動(dòng)核心區(qū)SLP變化的作用以及對(duì)北極濤動(dòng)的貢獻(xiàn)形式。本文的結(jié)果深入揭示了北歐海熱通量各個(gè)分量對(duì)海氣耦合的響應(yīng)和反饋,指出北歐海是1個(gè)潛熱主導(dǎo)的海氣耦合海域。

    2 熱通量的季節(jié)變化

    2.1 熱通量多年月平均分布特征

    熱通量和氣象數(shù)據(jù)取自NCEP-DOE Reanalysis-2數(shù)據(jù)集(http://www.esrl.noaa.gov/psd/dat a/gridded/data.ncep.reanalysis2.html#references),該數(shù)據(jù)集使用了高級(jí)分析與預(yù)報(bào)系統(tǒng)進(jìn)行數(shù)據(jù)同化,是NCEP Reanalysis-1的改進(jìn)版本。所有數(shù)據(jù)的網(wǎng)格如圖1所示,經(jīng)度間隔1.875°,緯度間隔1.903°。各個(gè)熱通量以進(jìn)入大氣的為正,向下輻射以到達(dá)海面為正,向上輻射以離開海面為正,潛熱和感熱以進(jìn)入大氣為正。凈短波輻射是指向下的短波輻射減去向上的短波輻射,凈長(zhǎng)波輻射是向上的長(zhǎng)波輻射減去向下的長(zhǎng)波輻射。同期的北極濤動(dòng)指數(shù)數(shù)據(jù)來自(http://www.cpc.ncep. noaa.gov/products/precip/CWlink/daily_ao_index/ao_index.html)。所使用的數(shù)據(jù)長(zhǎng)度為30年,從1979年1月~2010年12月。

    北歐海的海氣熱通量有非常顯著的季節(jié)變化。作者在圖2中給出了凈短波輻射、凈長(zhǎng)波輻射、感熱和潛熱在1月和7月的空間分布圖。

    北歐海大部分海域位于北極圈之內(nèi)。冬季(11、12、1月)大部分海域沒有短波輻射,夏季短波輻射很好的體現(xiàn)了緯度差異。在東格陵蘭流海域,每年大部分時(shí)間有海冰覆蓋[22],海冰較強(qiáng)的反照率導(dǎo)致凈短波輻射很低,大部分太陽輻射能被反射回太空。

    長(zhǎng)波輻射的冬夏差異并不大,冬季的凈長(zhǎng)波輻射略大一些,表明冬季的云量要顯著低于夏季,海洋失熱更為嚴(yán)重。在東格陵蘭流域,冬季凈長(zhǎng)波輻射略小一些,夏季凈長(zhǎng)波輻射略大一些。

    感熱通量在夏季比較弱,主要是海氣溫差較小的緣故。在冬季,發(fā)生非常強(qiáng)的感熱通量,尤其在格陵蘭海盆,感熱通量達(dá)到極大值,可達(dá)200 W/m2以上。在其他海域,感熱釋放量也達(dá)到80 W/m2以上。這樣大的感熱通量體現(xiàn)了該海域的海洋熱量通過感熱釋放給大氣,形成冬季大氣運(yùn)動(dòng)的主要熱源。這種強(qiáng)大的感熱釋放從10月持續(xù)到翌年4月。

    潛熱通量在夏季較弱,是由于夏季大氣濕度較大的緣故。而在冬季潛熱釋放量非常可觀。在挪威海流流域,潛熱釋放超過100 W/m2,而在回流區(qū),潛熱也達(dá)到60 W/m2以上,只有在東格陵蘭流海域,潛熱通量低于30 W/m2。顯然,北歐海冬季感熱和潛熱的釋放對(duì)歐洲冬季溫暖的氣候有重要貢獻(xiàn)。

    可見,在夏季,太陽短波輻射是北歐海大氣運(yùn)動(dòng)的主要熱源,凈長(zhǎng)波輻射通量與冬季差別不大,感熱和潛熱通量都很小。而在冬季,長(zhǎng)波輻射、感熱和潛熱提供給大氣的熱量接近300 W/m2,接近夏季太陽輻射提供的熱量,是該海域冬季大氣運(yùn)動(dòng)的主要熱源。

    2.2 4個(gè)海盆熱通量的季節(jié)變化

    圖3給出了4個(gè)海盆各個(gè)熱通量的季節(jié)變化情況。凈短波輻射夏季最大,其他3種熱通量都是夏季最小。各個(gè)海盆的凈短波輻射相差不大,只有格陵蘭海略高。凈長(zhǎng)波輻射通量也是格陵蘭海最高,其他3個(gè)海域相差不大。其實(shí),格陵蘭海的海表溫度要低于挪威海盆和羅弗敦海盆。感熱通量格陵蘭海是最高的,其他3個(gè)海盆相差不大。潛熱通量當(dāng)屬暖水區(qū)的挪威海盆和羅弗敦海盆最大,其他2個(gè)海盆較小,其中,冰島海的潛熱通量要遠(yuǎn)小于其他海域。值得注意的是,格陵蘭海夏季的潛熱接近0,表明該海域夏季氣溫低,濕度大,無法形成蒸發(fā)。

    圖1 北歐海地形和數(shù)據(jù)示意圖Fig.1 Topography of Nordic Seas and the grid of dada

    圖2 多年月平均1和7月凈短波輻射(a),凈長(zhǎng)波輻射(b),感熱(c)和潛熱(d)通量的空間分布Fig.2 Multiyear monthly average of net shortwave radiation(a),Net longwave radiation(b),Sensible heat(c)and latent heat(d)at January(left)and July(right),respectively

    圖3 4個(gè)海盆各個(gè)熱通量分量的季節(jié)變化(多年月平均)Fig.3 Seasonal variation(multiyear monthly average)of each heat flux component at four basins

    2.3各個(gè)海盆熱通量的方差

    從數(shù)據(jù)中去掉圖3中的多年月平均數(shù)據(jù),得到各個(gè)熱通量的距平量。統(tǒng)計(jì)各個(gè)海盆代表點(diǎn)距平量的方差列于表1??梢?,短波輻射的方差在各個(gè)海盆相差不多,只有在格陵蘭海盆高出50%,與該海域的云量相對(duì)變化較大有關(guān)。各個(gè)海域的長(zhǎng)波輻射方差相差不多,格陵蘭海高出大約40%,也是由于云量變幅較大的緣故。感熱的方差區(qū)域差別最大,暖水區(qū)的挪威海盆和羅弗敦海盆感熱的方差都略高于5 W·m-2,而在冰島海,感熱方差增大到8 W·m-2以上,在格陵蘭海盆,感熱方差驟升到18 W·m-2以上。這2個(gè)海域感熱方差變幅增大的原因是近海面氣溫變幅較大的緣故。潛熱方差在格陵蘭海變幅較大,達(dá)到8 W·m-2以上,而在其他3個(gè)海盆都在5 W·m-2左右,格陵蘭海潛熱變幅較大的原因是其濕度變化較大造成的。因此,格陵蘭海的感熱和潛熱通量不僅是海洋失去熱量的最主要方式,而且也是變幅最大的熱通量,是海洋和大氣發(fā)生較大幅度變化的主要原因。

    表1 各個(gè)海盆各熱通量的方差Table 1 Variances of heat fluxes in four basins/W·m2

    3 各個(gè)熱通量的空間差異

    3.1 各個(gè)熱通量的空間差異和時(shí)間變化

    為了體現(xiàn)各個(gè)海盆的季節(jié)變化,本文將圖1給出的4個(gè)海盆的中心點(diǎn)為核心,計(jì)算中心點(diǎn)與周邊海盆各個(gè)熱通量變化的相關(guān)性,以體現(xiàn)各個(gè)海盆之間的差異。圖4(a)~7(a)給出了各個(gè)海區(qū)相關(guān)系數(shù)大于0.6的海域,并計(jì)算各個(gè)海盆中心點(diǎn)熱通量的時(shí)間變化的差異(見圖4(b)~7(b))。高相關(guān)海域的范圍大,意味著其變化對(duì)更大的海域具有代表性。

    圖4(a)表明,挪威海盆和羅弗敦海盆凈短波輻射的相關(guān)系數(shù)場(chǎng)相互重疊,表明二者的空間一致性較高,表明這2個(gè)海盆的云量相差不大。而格陵蘭海和冰島海近乎自成體系。圖4(b)為凈短波輻射變化的距平,表現(xiàn)出格陵蘭海與眾不同的較大變幅,冰島海也與其他2個(gè)海盆有較明顯的差異。凈長(zhǎng)波輻射的空間一致性最好,4個(gè)海盆的相關(guān)系數(shù)場(chǎng)高度重疊(見圖5(a)),其時(shí)間變化的位相和幅度又有較好的一致性(見圖5(b))。

    圖4 凈短波輻射的區(qū)域差異和時(shí)間變化Fig.4 Spatial difference and temporal variation of net shortwave radiation

    圖5 凈長(zhǎng)波輻射的區(qū)域差異和時(shí)間變化Fig.5 Spatial difference and temporal variation of net longwave radiation

    從圖6(a)和圖7(a)可見,冰島海、挪威海盆和羅弗敦海盆的感熱和潛熱相關(guān)性都很好,相關(guān)系數(shù)場(chǎng)相互重疊的區(qū)域很大,表明暖流影響區(qū)域的感熱和潛熱特性比較接近。唯有格陵蘭海,其感熱和潛熱都與周邊不一樣,相關(guān)區(qū)域的范圍也很小,說明在格陵蘭海的感熱和潛熱特性主要集中在不大的范圍內(nèi)。圖6(b),7(b)也表明,格陵蘭海的感熱和潛熱都異常高,其中感熱產(chǎn)生的距平非常顯著,潛熱也是高于其他海盆。

    圖6 感熱通量的區(qū)域差異和時(shí)間變化Fig.6 Spatial difference and temporal variation of sensible heat flux

    圖7 潛熱通量的區(qū)域差異和時(shí)間變化Fig.7 Spatial difference and temporal variation of latent heat flux

    3.2各個(gè)熱通量變化的一致性

    表2給出了潛熱與其它各個(gè)輻射通量距平之間的相關(guān)系數(shù),以查看這些熱通量變化的一致性。結(jié)果表明,短波輻射與潛熱通量變化的相關(guān)性最小,主要原因是短波輻射主要發(fā)生在夏季,而潛熱主要發(fā)生在冬季。而潛熱與長(zhǎng)波輻射和感熱都有很好的正相關(guān),表明各種與海洋因素密切相關(guān)的熱通量都有很好的一致性。但是,潛熱與感熱的相關(guān)性更高一些,在格陵蘭海,二者的相關(guān)系數(shù)達(dá)到0.95;在挪威海盆,相關(guān)系數(shù)也達(dá)到0.82。而潛熱與長(zhǎng)波輻射的相關(guān)系數(shù)在0.51~0.66之間。

    表2 各個(gè)海盆潛熱與其他熱通量的相關(guān)系數(shù)Table 2 Correlation coefficients of latent heat with other heat fluxes at each basin

    總體來看,挪威海盆與羅弗敦海盆的熱通量變化一致,表明暖流區(qū)熱通量有比較一致的變化。冰島海雖然有時(shí)有比較大的偏差,但更趨向于暖流區(qū)的特征,與進(jìn)入冰島的伊爾明哥流有關(guān)。北冰島伊爾明哥流沿冰島北部沿岸流動(dòng),將大西洋水帶入冰島海,占據(jù)上層200~300 m的深度[23]。只有格陵蘭海非常特殊,其短波輻射、感熱和潛熱都有很大的變幅,與相鄰海域完全不同。這與格陵蘭海的特殊物理特征有關(guān)。在格陵蘭海,沒有暖流直接流入,強(qiáng)烈的感熱和潛熱釋放,導(dǎo)致海溫降低,氣溫也受北極冷空氣入侵的影響,形成與暖流區(qū)迥異的自然環(huán)境。此外,格陵蘭海的上升氣流最弱,云量少,造成凈短波輻射高,凈長(zhǎng)波輻射也高的特點(diǎn)。

    3.3各個(gè)熱通量變化的特殊事件

    從圖4(b)~7(b)可以看出,在過去30年發(fā)生了幾次顯著的熱量減少事件。其中,短波輻射的減少事件發(fā)生在1986—1987年和1999—2000年。這2次強(qiáng)負(fù)距平事件很可能與冰島的火山噴發(fā)有關(guān)。冰島在1986年7月發(fā)生了強(qiáng)火山噴發(fā),噴發(fā)期持續(xù)了大約15 d;冰島2000年的火山噴發(fā)也很強(qiáng)?;鹕絿姲l(fā)產(chǎn)生的火山灰和氣溶膠對(duì)太陽短波輻射的影響長(zhǎng)達(dá)一年以上[24]。

    感熱和潛熱通量的負(fù)距平事件與短波輻射有所不同。1986年冬季發(fā)生的感熱和潛熱的負(fù)距平事件,可以認(rèn)為也受到火山噴發(fā)的影響;而2000年發(fā)生火山噴發(fā),卻沒有發(fā)生感熱和潛熱的負(fù)距平。1998—1999年發(fā)生顯著的感熱和潛熱的負(fù)距平事件。云量的變化也表明,1987年的短波輻射負(fù)距平事件與云量增大有直接的關(guān)系。而在1998—1999年間,云量以低值為主。低云量通常對(duì)應(yīng)于高凈長(zhǎng)波輻射,但那年的長(zhǎng)波輻射也處于低值。因此,只有一種可能,就是海溫異常低。這種情形是合理的,因?yàn)楹氐蜁?huì)導(dǎo)致感熱和潛熱的下降。1998—1999年的異常低溫不會(huì)影響后續(xù)年份的短波輻射。

    從圖6和7可以看到,感熱和潛熱通量在1990—1996年之間一直保持高值,這段時(shí)間恰逢北極濤動(dòng)的正位相,表明感熱通量在這段時(shí)間與北極濤動(dòng)有很好的關(guān)系。除了特殊的事件之外,各個(gè)熱通量的變化更多地與北極濤動(dòng)相聯(lián)系。下文將仔細(xì)研究各個(gè)熱通量的變化與北極濤動(dòng)的關(guān)系。

    4 熱通量變化與北極濤動(dòng)指數(shù)的關(guān)系

    由于北歐海位于北極濤動(dòng)核心區(qū),北歐海輻射量的變化與北極濤動(dòng)應(yīng)該有密切的聯(lián)系。但是,由于北歐海的流系和水團(tuán)復(fù)雜,各個(gè)海盆的輻射特性相差很大。為此,本文在此討論各個(gè)海盆的輻射差異,以了解各個(gè)海盆對(duì)北極濤動(dòng)的貢獻(xiàn)。為了使結(jié)果對(duì)各個(gè)海盆有更好的代表性,本文選取各個(gè)海盆中所有點(diǎn)(見圖1中綠色三角所示)對(duì)各個(gè)熱通量進(jìn)行平均,將平均的結(jié)果與北極濤動(dòng)指數(shù)進(jìn)行相關(guān)分析。這里,采用了[25]北極濤動(dòng)核心區(qū)的概念,即北極濤動(dòng)指數(shù)與北歐海平均氣壓的高度一致性,用北極濤動(dòng)指數(shù)代表北歐海的平均氣壓場(chǎng)的變化。其中,由于考慮輻射通量與大氣的關(guān)系,本文采用向下的短波輻射代替凈短波輻射,因?yàn)橄蛳碌亩滩ㄝ椛渑c大氣吸收的關(guān)系更密切。分析結(jié)果見表3。

    表3 各海盆各熱通量與AO指數(shù)的相關(guān)性Table 3 Correlation coefficients of each heat flux with AO index at each basin

    從表3可見,向下短波輻射與AO指數(shù)相關(guān)性較高的在冰島海,其他海域相關(guān)性都比較差,表明發(fā)生在冰島海的短波輻射與AO的關(guān)系比較密切。在冰島海和其他2個(gè)暖水海盆凈長(zhǎng)波輻射與AO有較高的相關(guān)性。格陵蘭海的感熱和潛熱與北極濤動(dòng)指數(shù)正相關(guān)程度較高,在其他海盆,感熱和潛熱的變化與北極濤動(dòng)關(guān)系不明顯。表3的結(jié)果表明,向下短波輻射和凈長(zhǎng)波輻射與AO相關(guān)性最好的都在冰島海,CC均為-0.35。感熱和潛熱與AO相關(guān)最好的在格陵蘭海,CC均為0.30(見圖8)。

    上述結(jié)果表明,以下2個(gè)過程是發(fā)生在北歐海的重要海氣相互作用過程:第一,冰島海的短波輻射與AO指數(shù)相關(guān)度較高,是與北極濤動(dòng)關(guān)系密切的過程。短波輻射的變化通常與云量高度相關(guān),4個(gè)海盆向下短波輻射變化與云量的相關(guān)系數(shù)分別為-0.74,-0.74,-0.65和-0.71。因此,云量是影響北極濤動(dòng)最重要因子之一。第二,發(fā)生在暖水海盆中的長(zhǎng)波輻射與AO有密切關(guān)系,保持大范圍的較高相關(guān)度。第三,格陵蘭海的潛熱和感熱通量的變化與AO正相關(guān),有可能是北極濤動(dòng)的驅(qū)動(dòng)因素,至少與AO耦合變化,值得深入研究。

    本文通過對(duì)2個(gè)時(shí)間序列進(jìn)行滑動(dòng)相關(guān)來表現(xiàn)二者的一致程度。設(shè)Xk和Yk2個(gè)時(shí)間序列,Ri是2種信號(hào)以第i時(shí)刻為中心的相關(guān)系數(shù),

    計(jì)算相關(guān)系數(shù)的時(shí)間窗口為±n。對(duì)各點(diǎn)滑動(dòng)求取相關(guān)系數(shù),就可得到相關(guān)系數(shù)的滑動(dòng)變化曲線。在本文中,取n=12個(gè)月,相當(dāng)于對(duì)2年的數(shù)據(jù)取相關(guān)。圖8給出了與北極濤動(dòng)相關(guān)系數(shù)較高的這4個(gè)熱通量的滑動(dòng)相關(guān)系數(shù)。

    在30年的時(shí)間里,相關(guān)系數(shù)0.3~0.35超過了99.9%的置信度,已經(jīng)體現(xiàn)了顯著的相關(guān)性。但從圖8和表3可見,如果將這30年分為兩段,1989—1992和1993—2009,兩段的相關(guān)系數(shù)相差很大。冰島海前段的短波和長(zhǎng)波輻射與AO的相關(guān)系數(shù)達(dá)到-0.59和-0.63,而后段相關(guān)性很差;而格陵蘭海正好相反,感熱和潛熱前段的相關(guān)性很差,而后段的相關(guān)系數(shù)竟高達(dá)0.72和0.73。這個(gè)結(jié)果表明,在1992年以前,冰島海的短波輻射和長(zhǎng)波輻射對(duì)AO的貢獻(xiàn)可能更大一些,而1993年以后,格陵蘭海的感熱和潛熱對(duì)AO的貢獻(xiàn)更大一些。作者認(rèn)為,這可能與較長(zhǎng)時(shí)間尺度內(nèi)北極濤動(dòng)的位相變化有關(guān)。1989—1995是北極濤動(dòng)的正位相期,1992年以前AO總體上在升高,而1993年以后AO總體上在降低,各個(gè)熱通量的在相關(guān)性上的差別可能與此有關(guān)。對(duì)這個(gè)現(xiàn)象還無法在物理上給出合理的解釋,需要在此基礎(chǔ)上深入研究。

    圖8 與北極濤動(dòng)相關(guān)性較好的4個(gè)熱通量Fig.8 Heat fluxes well correlated with AO index

    5 結(jié)果與討論

    北歐海包含有4個(gè)海盆,格陵蘭海、冰島海、挪威海盆和羅弗敦海盆。北歐海有暖流和寒流注入,又發(fā)生大量回流,其環(huán)流和水團(tuán)特性異常復(fù)雜。由于北歐海的環(huán)流受地形控制,其水團(tuán)的分布與海盆的分布有密切的關(guān)系。因此,本文研究各個(gè)海盆熱通量變化的差異,以研究獲取對(duì)北歐海海氣相互作用的全面認(rèn)識(shí),并通過研究這些熱通量與北極濤動(dòng)指數(shù)的相關(guān)性,研究北歐??赡艿暮怦詈线^程。

    由于特殊的流場(chǎng)特點(diǎn)和地域特征,北歐海短波輻射主要發(fā)生在夏季,而來自海洋的長(zhǎng)波輻射、感熱和潛熱通量主要發(fā)生在冬季。海洋在冬季強(qiáng)烈地向大氣輸送熱量,與夏季太陽輻射的熱量相當(dāng)。海洋輸送的熱量在各個(gè)海盆很不相同,其中,格陵蘭海的月平均長(zhǎng)波輻射和感熱最大。用方差表現(xiàn)各個(gè)熱通量的變化幅度得出,最為特殊的是格陵蘭海,其短波輻射高出50%,長(zhǎng)波輻射高出大約40%,潛熱高出大約60%,感熱驟升到18 W·m-2以上,高出其他海盆近4倍。

    各海盆熱通量變化的相關(guān)性分析表明,挪威海盆與羅弗敦海盆的熱通量變化比較一致,表明暖流區(qū)有一致的變化。冰島海雖然有時(shí)有比較大的偏差,但更趨向于暖流區(qū)的特征。只有格陵蘭海非常特殊,其短波輻射、感熱和潛熱都有很大的變幅,與相鄰海域完全不同。在格陵蘭海,強(qiáng)烈的感熱和潛熱釋放導(dǎo)致海溫降低,氣溫也深受北極冷空氣入侵的影響,形成與暖流區(qū)迥異的自然環(huán)境。此外,格陵蘭海云量少,造成凈短波輻射高,凈長(zhǎng)波輻射也高的特點(diǎn)。

    在過去30年發(fā)生了幾次顯著的熱量減少事件,其中,短波輻射的減少事件發(fā)生在1986—1987年和1999—2000年,可能與1987年和2000年冰島的火山噴發(fā)有關(guān)。由于到達(dá)的太陽輻射減弱,感熱和潛熱也同步減少。感熱和潛熱的變化也與北極海冰的輸出量增大,海面溫度偏低有顯著關(guān)系。

    文章分析了4個(gè)主要海盆熱通量的變化與北極濤動(dòng)指數(shù)的關(guān)系,結(jié)果表明,發(fā)生在冰島海的向下短波輻射和發(fā)生在格陵蘭海的感熱和潛熱與AO相關(guān)度較高,體現(xiàn)了與AO的耦合變化特性。由于到達(dá)海面的短波輻射與云量變化高度相關(guān),短波輻射與AO的高相關(guān)是云量變化的結(jié)果,因此,AO的驅(qū)動(dòng)因素是云量的變化。而云量主要是海洋的蒸發(fā)和大氣的上升氣流生成的,北歐海位于冰島低壓區(qū),發(fā)生在格陵蘭海的潛熱通量過程可能與云的形成有密切關(guān)系,成為北極濤動(dòng)的驅(qū)動(dòng)因素。本文的統(tǒng)計(jì)分析只能指出這種可能性,證實(shí)這個(gè)機(jī)理需要更多的動(dòng)力學(xué)研究。

    這些熱通量與AO指數(shù)的滑動(dòng)相關(guān)系數(shù)揭示了一個(gè)現(xiàn)象,1992年以前冰島海的短波輻射與AO的相關(guān)性非常高,1993年后相關(guān)性較差;反之,格陵蘭海的感熱和潛熱在1992年以前與AO的相關(guān)性并不好,而1993年后與AO高度相關(guān),表明在AO增強(qiáng)過程中主要與輻射特性有關(guān),而AO減弱過程主要與海洋的感熱和潛熱有關(guān)。

    [1] Mauritzen C.Production of dense overflow waters feeding the North Atlantic across the Greenland-Scotland Ridge.Part 2:An inverse model[J].Deep Sea Research I,1996b,43(6):807-835.

    [2] Rudels B,F(xiàn)ahrbach E,Meincke J,et al.The East Greenland Current and its contribution to the Denmark Strait overflow[J].ICES Journal of Marine Science,2002,59(6):1133-1154.

    [3] Meincke J,Rudels B,F(xiàn)riedrich H.The Arctic Ocean-Nordic Seas thermohaline system[J].ICES Journal of Marine Science,1997,54(3):283-299.

    [4] Swift J H,Aagaard K.Seasonal transitions and water mass formation in the Iceland and Greenland seas[J].Deep Sea Research Part A.Oceanographic Research Papers,1981,28(10):1107-1129.

    [5] Rossby T,Prater M D,S?iland H.Pathways of inflow and dispersion of warm waters in the Nordic seas[J].J Geophys Res,2009,114,C04011.doi:10.1029/2008JC005073.

    [6] Van Aken H M,Budéus G,Héhnel M.The anatomy of the Arctic Frontal Zone in the Greenland Sea[J].J Geophys Res,1995,100(C8):15999-16014.doi:10.1029/95JC01176.

    [7] Kostianoy A G,Nihoul J C J.Frontal Zones in the Norwegian,Greenland,Barents and Bering Seas[M].∥Nihoul J C J,Kostianoy A G,eds.Influence of Climate Change on the Changing Arctic and Sub-Arctic Conditions.Netherlands:Springer,2009.

    [8] 何琰,趙進(jìn)平.北歐海的鋒面分布特征及其季節(jié)變化[J].地球科學(xué)進(jìn)展,2011,26(10):1079-1091.

    [9] Marshall J,Schott F.Open-ocean convection:Observations,theory,and models[J].Rev Geophys,1999,37(1):1-64.

    [10] Cayan D R.Latent and sensible heat flux anomalies over the northern oceans:the connection to monthly atmospheric circulation[J].J of Clim,1992b,5(16):354-369.

    [11] 邵秋麗,趙進(jìn)平.北歐海深層水形成的研究進(jìn)展[J].地球科學(xué)進(jìn)展,2014,29(1):42-45.

    [12] Cottier F R,Venables E J.On the double-diffusive and cabbeling environment of the Arctic Front,West Spitsbergen[J].Polar Research,2007,26(2):152-159.

    [13] 史文奇、趙進(jìn)平.北冰洋水體對(duì)格陵蘭海混合增密對(duì)流的可能影響分析[J].海洋學(xué)報(bào),2012,34(6):19-29.

    [14] Rudels B,Bj?rk G,Nilsson J,et al.Theinteraction between waters from the Arctic Ocean and the Nordic Seas north of Fram Strait and along the East Greenland Current:results from the Arctic Ocean-02 Oden expedition[J].Journal of Marine Systems,2005,55(1):1-30.

    [15] Hurrell JW.Decadal trends in the North Atlantic Oscillation:regional temperatures and precipitation[J].Science,1995,269:676-679.

    [16] Thompson D W J,Wallace J M.The Arctic oscillation signature in wintertime geopotential height and temperature fields[J].Geophys Res Lett,1998,25:1297-1300.

    [17] Cayan D R.Latent and sensible heat flux anomalies over the northern oceans:driving the sea surface temperature[J].J of Phys Oceanogr,1992a,22(17):859-881.

    [18] Kawamura R.A rotated EOF analysis of global sea surface temperature variability with interannual and interdecadal scales[J].J Phys Oceanogr,1994,24:707-715.

    [19] Deser C,Blackmon M L.Surface climate variability over the North Atlantic Ocean during winter:1900-1989[J].J Climate,1993,6:1743-1753.

    [20] Slonosky V C,Mysak L A,Derome J.Linking arctic sea-ice and atmospheric circulation anomalies on interannual and decadal time scales[J].Atmosphere-Ocean,1997,35(3):333-366.

    [21] Zhao Jinping,Cao Yong,Shi Jiuxin.Core region of Arctic Oscillation and the main atmospheric events impact on the Arctic[J].Geophys Res Lett,2006,33:L22708.doi:10.1029/2006GL027590.

    [22] 牟龍江,趙進(jìn)平.格陵蘭海海冰外緣線變化特征分析[J].地球科學(xué)進(jìn)展,2013,28(6):709-717.

    [23] Blindheim J,?sterhus S.The Nordic Seas,Main Oceanographic Features[M].//Drange H,ed.The Nordic seas:an integrated perspective oceanoraphy,climatology,biogeochemistry,and modeling.Washington D.C:AGU,2005:11-37.

    [24] Moune S,Sigmarsson O,Thordarson T,et al.Recent volatile evolution in the magmatic system of Hekla volcano[J].Earth and Planetary Science Letters,2007,255:373-389.

    [25] Zhao Jinping,Cao Yong,Shi Jiuxin.Spatial variation of the Arctic Oscillation and its long-term change[J].Tellus,2010,62A:661-672.

    Multiyear Variation of the Main Heat Flux Components in the Four Basins of Nordic Seas

    ZHAO Jin-Ping1,Ken Drinkwater2
    (1.Ocean University of China,Qingdao 266100,China;2.Institute of Marine Research,Bergen N-5024,Norway)

    Water masses in the Nordic Seas are complex as there are imports of warm and cold ocean current,and a number of recycles.The spatial distribution of water masses is closely related with the four basins,as the water transport in Nordic Seas are carried by topography-followed currents.The differences of heat flux components in the four basins are studied to understand the spatial difference of the air-sea interaction in the Nordic Seas.The heat flux components of Greenland Sea are quite different,whose variance of shortwave radiation is 50%higher,longwave radiation is 40%higher,latent heat is 60%higher,and the sensible heat is four times than the other basins.The possible reason is its low sea surface temperature reduced by enhanced release of sensible and latent,and its low air temperature influenced by the cold air mass from the Arctic.As the highly consistent of the averaged SLP of Nordic Seas with AO index,the correlation of the heat fluxes of the four basins with AO index is analyzed to show the influence of oceanic heat on atmospheric movement.The high correlation with AOindex appeared in both regions:the shortwave radiation in Icelandic Sea,and the sensible and latent heat flux in the Greenland Sea.Running correlation coefficient indicates that the high correlation processes could be divided into two periods:before 1992,the shortwave radiation in Icelandic Sea is highly correlated with AO index;After 1993,the sensible and latent heat in the Greenland Sea are closely correlated with AO index.The result of this study supports the following possibility:the sensible and latent heat release in the Greenland Sea is dominant to the atmospheric processes in the Icelandic area through the upwelled air flow to influence the cloud cover,and hence,the arriving solar radiation.

    Nordic Seas;heat flux;air-sea coupling;Arctic Oscillation(AO);AO core region

    P731.21

    A

    1672-5174(2014)10-009-11

    責(zé)任編輯 龐 旻

    國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目“北極海冰與上層海洋環(huán)流耦合變化及其氣候效應(yīng)(41330960)”;南北極環(huán)境綜合考察與評(píng)估專項(xiàng)(CHINARE0403)資助

    2014-08-20;

    2014-09-08

    趙進(jìn)平(1954-),男,教授,博導(dǎo)。E-mail:jpzhao@ouc.edu.cn

    猜你喜歡
    海盆格陵蘭長(zhǎng)波
    廣義對(duì)稱正則長(zhǎng)波方程的孤波解和周期波解及它們與Hamilton能量的關(guān)系
    鶯歌海盆地超高溫高壓井?dāng)D水泥承壓堵漏技術(shù)
    法爾維海盆構(gòu)造特征及演化
    冰雪王國(guó)格陵蘭
    公民與法治(2020年2期)2020-05-30 12:28:54
    格陵蘭也曾“脫歐”
    特朗普為何執(zhí)意買格陵蘭
    東太平洋秘魯海盆海域夏季聲速特點(diǎn)分析
    1982年,格陵蘭公投“脫歐”
    基于構(gòu)架點(diǎn)頭角速度的軌道垂向長(zhǎng)波不平順在線檢測(cè)
    掃描型長(zhǎng)波紅外連續(xù)變焦光學(xué)系統(tǒng)
    精品国产美女av久久久久小说| 国产高清视频在线播放一区| 国产成年人精品一区二区| 波多野结衣巨乳人妻| 欧美黄色淫秽网站| 亚洲一卡2卡3卡4卡5卡精品中文| 岛国在线观看网站| 免费在线观看亚洲国产| 中文字幕熟女人妻在线| 久久久久久久精品吃奶| 久久精品成人免费网站| 中出人妻视频一区二区| 淫妇啪啪啪对白视频| 啦啦啦韩国在线观看视频| 激情在线观看视频在线高清| 日韩 欧美 亚洲 中文字幕| a级毛片a级免费在线| 久久九九热精品免费| 欧美中文综合在线视频| 婷婷精品国产亚洲av| 亚洲美女黄片视频| 日本免费a在线| xxxwww97欧美| 又黄又粗又硬又大视频| 国产又色又爽无遮挡免费看| а√天堂www在线а√下载| 久久久久性生活片| 国产精品亚洲美女久久久| 亚洲专区中文字幕在线| 99在线人妻在线中文字幕| 成人手机av| 国产人伦9x9x在线观看| 亚洲精品在线美女| 亚洲,欧美精品.| 69av精品久久久久久| 中出人妻视频一区二区| av福利片在线观看| 精品乱码久久久久久99久播| 久久久国产成人免费| 久久久久久久久免费视频了| 欧美日韩黄片免| 日本 av在线| 免费搜索国产男女视频| 中文字幕最新亚洲高清| 999精品在线视频| 亚洲精品色激情综合| 一本久久中文字幕| 国产激情欧美一区二区| 午夜福利在线在线| 一夜夜www| 成年免费大片在线观看| 欧美激情久久久久久爽电影| 中文字幕人成人乱码亚洲影| 1024视频免费在线观看| 高潮久久久久久久久久久不卡| 欧美三级亚洲精品| 国产亚洲欧美在线一区二区| 成年免费大片在线观看| 少妇熟女aⅴ在线视频| 国产在线观看jvid| 国产亚洲精品久久久久5区| 国产97色在线日韩免费| 亚洲av五月六月丁香网| 国产精品久久久久久精品电影| 色精品久久人妻99蜜桃| 国产成人啪精品午夜网站| 亚洲av五月六月丁香网| 五月玫瑰六月丁香| 在线看三级毛片| 韩国av一区二区三区四区| 人人妻,人人澡人人爽秒播| 国产精品爽爽va在线观看网站| 亚洲人与动物交配视频| 91九色精品人成在线观看| 婷婷精品国产亚洲av在线| 国产精品自产拍在线观看55亚洲| 亚洲欧美精品综合久久99| 亚洲欧美日韩东京热| 国产真人三级小视频在线观看| 亚洲七黄色美女视频| 欧美成人免费av一区二区三区| 欧美乱色亚洲激情| 国产av一区二区精品久久| 久久精品91无色码中文字幕| 怎么达到女性高潮| 亚洲精品一区av在线观看| 天天添夜夜摸| 熟女少妇亚洲综合色aaa.| 曰老女人黄片| 岛国视频午夜一区免费看| 青草久久国产| 在线国产一区二区在线| 男人舔奶头视频| 午夜福利免费观看在线| 三级毛片av免费| 亚洲av成人一区二区三| 日韩 欧美 亚洲 中文字幕| 人成视频在线观看免费观看| 精品第一国产精品| 精品电影一区二区在线| 国产一区二区三区在线臀色熟女| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| av在线播放免费不卡| 啪啪无遮挡十八禁网站| 欧美久久黑人一区二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲性夜色夜夜综合| 成年免费大片在线观看| 麻豆成人av在线观看| 黄片小视频在线播放| 午夜精品在线福利| 日韩欧美在线乱码| 国产精品香港三级国产av潘金莲| 俺也久久电影网| 国产成人av激情在线播放| 欧美一区二区精品小视频在线| a级毛片a级免费在线| 国产成人欧美在线观看| 中文亚洲av片在线观看爽| 亚洲男人天堂网一区| 99精品久久久久人妻精品| 中文字幕高清在线视频| 国产精品香港三级国产av潘金莲| 国产精华一区二区三区| 黄色毛片三级朝国网站| 国产三级在线视频| 亚洲国产精品成人综合色| √禁漫天堂资源中文www| 亚洲熟女毛片儿| www.999成人在线观看| 亚洲成人久久爱视频| 国产伦一二天堂av在线观看| 亚洲男人天堂网一区| 亚洲人与动物交配视频| 国产黄片美女视频| 日韩免费av在线播放| 草草在线视频免费看| 美女黄网站色视频| 国产亚洲欧美在线一区二区| 亚洲国产精品成人综合色| 99久久无色码亚洲精品果冻| 毛片女人毛片| 校园春色视频在线观看| 熟妇人妻久久中文字幕3abv| 又紧又爽又黄一区二区| 人妻夜夜爽99麻豆av| 少妇人妻一区二区三区视频| a级毛片在线看网站| 欧美大码av| 老熟妇仑乱视频hdxx| 听说在线观看完整版免费高清| 99久久综合精品五月天人人| 99久久久亚洲精品蜜臀av| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看| 每晚都被弄得嗷嗷叫到高潮| 色播亚洲综合网| 亚洲人成网站高清观看| 在线视频色国产色| 亚洲欧美日韩高清专用| 少妇熟女aⅴ在线视频| 亚洲成人国产一区在线观看| 午夜激情福利司机影院| 亚洲人成网站在线播放欧美日韩| 欧美一区二区国产精品久久精品 | 成人国语在线视频| 午夜激情av网站| 波多野结衣高清作品| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 国产av不卡久久| 亚洲成a人片在线一区二区| 亚洲欧美日韩东京热| 在线免费观看的www视频| 18禁黄网站禁片午夜丰满| 老司机深夜福利视频在线观看| 两个人免费观看高清视频| 国产欧美日韩一区二区精品| 又黄又爽又免费观看的视频| 亚洲国产欧美人成| 色av中文字幕| 999久久久国产精品视频| av视频在线观看入口| 最近最新中文字幕大全电影3| 大型黄色视频在线免费观看| 18美女黄网站色大片免费观看| 日日夜夜操网爽| 欧美av亚洲av综合av国产av| 99精品久久久久人妻精品| 欧美性长视频在线观看| 国产午夜精品久久久久久| 中出人妻视频一区二区| 亚洲精品中文字幕在线视频| 黄频高清免费视频| 真人做人爱边吃奶动态| 久久中文看片网| x7x7x7水蜜桃| 久久精品91蜜桃| 成在线人永久免费视频| 一级毛片女人18水好多| 国产精品av视频在线免费观看| 国产欧美日韩一区二区精品| 国产精品一区二区精品视频观看| 精品久久久久久久毛片微露脸| 蜜桃久久精品国产亚洲av| 少妇粗大呻吟视频| 亚洲一码二码三码区别大吗| 精品高清国产在线一区| 亚洲18禁久久av| 国产精品亚洲av一区麻豆| 国产亚洲精品久久久久久毛片| a级毛片a级免费在线| 国产精品永久免费网站| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 亚洲午夜精品一区,二区,三区| 国产黄a三级三级三级人| 国产单亲对白刺激| 大型黄色视频在线免费观看| 91国产中文字幕| 国产成+人综合+亚洲专区| 宅男免费午夜| 亚洲av电影不卡..在线观看| 波多野结衣高清无吗| 日本 欧美在线| 9191精品国产免费久久| 黄色视频不卡| 久久午夜综合久久蜜桃| 国产精品香港三级国产av潘金莲| 国产97色在线日韩免费| 成人国产一区最新在线观看| 好男人在线观看高清免费视频| 特级一级黄色大片| 麻豆国产97在线/欧美 | 国产真人三级小视频在线观看| 麻豆成人午夜福利视频| 女同久久另类99精品国产91| 狠狠狠狠99中文字幕| 欧美最黄视频在线播放免费| 18禁国产床啪视频网站| 亚洲男人的天堂狠狠| 亚洲 国产 在线| 亚洲黑人精品在线| 搞女人的毛片| 极品教师在线免费播放| 一级毛片高清免费大全| 国产成年人精品一区二区| 人妻久久中文字幕网| 国产高清视频在线播放一区| 给我免费播放毛片高清在线观看| 曰老女人黄片| 性色av乱码一区二区三区2| 91大片在线观看| 一级a爱片免费观看的视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区视频在线观看免费| 国产精品久久久久久人妻精品电影| 黄色毛片三级朝国网站| 免费高清视频大片| 成人18禁高潮啪啪吃奶动态图| 五月伊人婷婷丁香| 国产区一区二久久| 男人舔奶头视频| 91国产中文字幕| 日韩精品中文字幕看吧| 91av网站免费观看| 夜夜躁狠狠躁天天躁| 亚洲片人在线观看| 免费看日本二区| 国内少妇人妻偷人精品xxx网站 | 一级毛片女人18水好多| 日韩欧美三级三区| svipshipincom国产片| 亚洲成人精品中文字幕电影| 中文字幕精品亚洲无线码一区| 欧美成人午夜精品| 国产精品久久视频播放| 亚洲va日本ⅴa欧美va伊人久久| 变态另类丝袜制服| 日本一本二区三区精品| 日韩欧美一区二区三区在线观看| 亚洲国产中文字幕在线视频| 欧美在线黄色| 非洲黑人性xxxx精品又粗又长| 成人特级黄色片久久久久久久| 99国产精品一区二区蜜桃av| 日韩欧美在线乱码| 90打野战视频偷拍视频| 一级毛片高清免费大全| 久久久久久久午夜电影| 国产成人av激情在线播放| 久久婷婷人人爽人人干人人爱| 久久久久性生活片| 婷婷精品国产亚洲av在线| 啦啦啦观看免费观看视频高清| av中文乱码字幕在线| 俺也久久电影网| 亚洲色图av天堂| 国产精品av久久久久免费| 窝窝影院91人妻| 国产精品一及| 91大片在线观看| 亚洲精品久久成人aⅴ小说| 日本撒尿小便嘘嘘汇集6| 色哟哟哟哟哟哟| 欧美成人一区二区免费高清观看 | 国产在线观看jvid| 两人在一起打扑克的视频| 99久久久亚洲精品蜜臀av| 桃红色精品国产亚洲av| 欧美日韩国产亚洲二区| 精品乱码久久久久久99久播| 精品久久久久久久久久久久久| 日本免费一区二区三区高清不卡| 国产成人精品久久二区二区91| 亚洲av中文字字幕乱码综合| 欧美一区二区国产精品久久精品 | 久久婷婷人人爽人人干人人爱| 99久久综合精品五月天人人| 国产一区在线观看成人免费| 老司机午夜福利在线观看视频| 免费看a级黄色片| 国产成人av教育| 麻豆久久精品国产亚洲av| 在线观看午夜福利视频| 熟女电影av网| 嫁个100分男人电影在线观看| 激情在线观看视频在线高清| 香蕉av资源在线| 午夜激情福利司机影院| 脱女人内裤的视频| 香蕉久久夜色| 日韩欧美国产一区二区入口| or卡值多少钱| 久久久久久亚洲精品国产蜜桃av| 日韩精品中文字幕看吧| 丰满人妻一区二区三区视频av | 老司机在亚洲福利影院| 99精品欧美一区二区三区四区| 无人区码免费观看不卡| 激情在线观看视频在线高清| 长腿黑丝高跟| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 制服丝袜大香蕉在线| 级片在线观看| 欧美又色又爽又黄视频| 亚洲黑人精品在线| 真人做人爱边吃奶动态| 少妇被粗大的猛进出69影院| 欧美色欧美亚洲另类二区| 此物有八面人人有两片| 亚洲欧美激情综合另类| 久久国产乱子伦精品免费另类| 99热6这里只有精品| 999久久久国产精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 大型黄色视频在线免费观看| 在线看三级毛片| 国产精品影院久久| 欧美黄色淫秽网站| 五月伊人婷婷丁香| 日韩欧美在线二视频| avwww免费| 午夜福利成人在线免费观看| 亚洲国产看品久久| 国产久久久一区二区三区| 午夜a级毛片| 精品久久久久久成人av| 草草在线视频免费看| 日本在线视频免费播放| 欧美午夜高清在线| 中文在线观看免费www的网站 | 啦啦啦韩国在线观看视频| 国产av不卡久久| 日日夜夜操网爽| 国产一区二区在线av高清观看| 午夜精品一区二区三区免费看| www.999成人在线观看| 午夜精品一区二区三区免费看| 国产精品一区二区三区四区久久| 99精品在免费线老司机午夜| 久久久久久大精品| 特级一级黄色大片| 日韩av在线大香蕉| 久久香蕉精品热| 成人手机av| 黄色 视频免费看| 好男人在线观看高清免费视频| 两个人的视频大全免费| 99久久国产精品久久久| 亚洲乱码一区二区免费版| 又爽又黄无遮挡网站| 欧美av亚洲av综合av国产av| 亚洲国产中文字幕在线视频| 免费电影在线观看免费观看| 欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 亚洲一区二区三区不卡视频| 成人午夜高清在线视频| 女同久久另类99精品国产91| 最近在线观看免费完整版| 999久久久国产精品视频| 中文字幕精品亚洲无线码一区| 99国产极品粉嫩在线观看| 可以在线观看的亚洲视频| 欧美+亚洲+日韩+国产| 很黄的视频免费| 日本五十路高清| 成人国语在线视频| 国产精品av视频在线免费观看| 亚洲欧洲精品一区二区精品久久久| 国产欧美日韩一区二区三| 欧美大码av| 一二三四社区在线视频社区8| 国产精品电影一区二区三区| 俄罗斯特黄特色一大片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品成人综合色| 日本成人三级电影网站| 久久久国产成人免费| 又粗又爽又猛毛片免费看| 夜夜夜夜夜久久久久| 久久国产乱子伦精品免费另类| 成熟少妇高潮喷水视频| 国产精品免费一区二区三区在线| 免费观看精品视频网站| av在线天堂中文字幕| 757午夜福利合集在线观看| 国产精品久久久久久人妻精品电影| 一级黄色大片毛片| 久久香蕉激情| 人人妻人人澡欧美一区二区| 亚洲男人的天堂狠狠| 国产91精品成人一区二区三区| 91老司机精品| 国产午夜精品久久久久久| 久久精品国产综合久久久| 国产伦人伦偷精品视频| 在线观看66精品国产| 香蕉国产在线看| ponron亚洲| 国产69精品久久久久777片 | 久久精品成人免费网站| 国产精品久久久人人做人人爽| 国产成人啪精品午夜网站| 亚洲,欧美精品.| 成人欧美大片| 国产精品av久久久久免费| 亚洲午夜精品一区,二区,三区| 日韩大码丰满熟妇| 日本免费一区二区三区高清不卡| www国产在线视频色| 国产av一区二区精品久久| 久久久精品大字幕| or卡值多少钱| 免费无遮挡裸体视频| 国产亚洲av嫩草精品影院| 天天躁夜夜躁狠狠躁躁| 欧美中文日本在线观看视频| 草草在线视频免费看| 亚洲在线自拍视频| 欧美性猛交╳xxx乱大交人| 天天一区二区日本电影三级| 黄片大片在线免费观看| 黑人操中国人逼视频| 国产高清有码在线观看视频 | 国产高清有码在线观看视频 | 人成视频在线观看免费观看| 亚洲中文字幕一区二区三区有码在线看 | 91麻豆精品激情在线观看国产| 免费在线观看日本一区| 久久人妻福利社区极品人妻图片| 三级男女做爰猛烈吃奶摸视频| 亚洲色图 男人天堂 中文字幕| 亚洲av成人不卡在线观看播放网| 亚洲精品国产一区二区精华液| 久久久久久久久久黄片| 亚洲国产欧美网| 国产高清有码在线观看视频 | 亚洲国产精品成人综合色| av福利片在线观看| 久久精品91无色码中文字幕| av免费在线观看网站| 亚洲五月婷婷丁香| 91字幕亚洲| 久久香蕉激情| 天天躁夜夜躁狠狠躁躁| 亚洲欧美精品综合久久99| 久久中文看片网| 午夜日韩欧美国产| 亚洲精华国产精华精| 老熟妇仑乱视频hdxx| 久久久久国产一级毛片高清牌| 黄片小视频在线播放| 亚洲av成人一区二区三| 中亚洲国语对白在线视频| 麻豆久久精品国产亚洲av| 欧美丝袜亚洲另类 | 久久天躁狠狠躁夜夜2o2o| 一进一出抽搐gif免费好疼| 欧美中文日本在线观看视频| 黄片小视频在线播放| 亚洲成av人片在线播放无| 欧美黑人巨大hd| 成人av一区二区三区在线看| 免费看日本二区| 亚洲av美国av| 搡老熟女国产l中国老女人| 夜夜爽天天搞| 悠悠久久av| 久久久久国产一级毛片高清牌| 久久久久久久久久黄片| 国产成+人综合+亚洲专区| a在线观看视频网站| 亚洲av成人不卡在线观看播放网| 色尼玛亚洲综合影院| 色综合站精品国产| 国语自产精品视频在线第100页| 麻豆成人午夜福利视频| 美女午夜性视频免费| 男女那种视频在线观看| 国产一区二区在线av高清观看| 国产1区2区3区精品| 亚洲无线在线观看| 99国产精品99久久久久| 女生性感内裤真人,穿戴方法视频| 久久午夜综合久久蜜桃| 一a级毛片在线观看| www.熟女人妻精品国产| 日韩欧美三级三区| 高潮久久久久久久久久久不卡| 国产高清视频在线观看网站| 亚洲天堂国产精品一区在线| 久久精品国产清高在天天线| 最近最新中文字幕大全电影3| 久久久久久久久久黄片| 夜夜躁狠狠躁天天躁| 变态另类丝袜制服| 99国产精品99久久久久| 久久久水蜜桃国产精品网| 一级毛片高清免费大全| 亚洲狠狠婷婷综合久久图片| 韩国av一区二区三区四区| 久久人妻福利社区极品人妻图片| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成人久久爱视频| 精品人妻1区二区| 国产视频内射| 久久国产乱子伦精品免费另类| 色综合欧美亚洲国产小说| 久久精品aⅴ一区二区三区四区| 级片在线观看| 十八禁网站免费在线| 免费一级毛片在线播放高清视频| videosex国产| 国产伦一二天堂av在线观看| 成年人黄色毛片网站| 五月伊人婷婷丁香| 免费一级毛片在线播放高清视频| 色精品久久人妻99蜜桃| 人人妻人人澡欧美一区二区| 91老司机精品| 大型av网站在线播放| 亚洲 国产 在线| 丰满人妻熟妇乱又伦精品不卡| 国产精品自产拍在线观看55亚洲| 日本黄色视频三级网站网址| 成年女人毛片免费观看观看9| 哪里可以看免费的av片| 手机成人av网站| bbb黄色大片| 久久精品国产亚洲av香蕉五月| 黑人操中国人逼视频| 成人手机av| 国产不卡一卡二| 亚洲av中文字字幕乱码综合| 又大又爽又粗| 午夜影院日韩av| 怎么达到女性高潮| 母亲3免费完整高清在线观看| 国产精品影院久久| 亚洲av成人精品一区久久| 一个人免费在线观看电影 | 亚洲精品粉嫩美女一区| 国产精品98久久久久久宅男小说| 两人在一起打扑克的视频| 波多野结衣高清作品| 亚洲av五月六月丁香网| 免费在线观看成人毛片| 亚洲18禁久久av| 成年版毛片免费区| 亚洲成a人片在线一区二区| 亚洲国产欧洲综合997久久,| 亚洲色图av天堂| 亚洲成a人片在线一区二区| 亚洲精品美女久久久久99蜜臀| 曰老女人黄片| 99久久精品国产亚洲精品| 亚洲av美国av| 成年版毛片免费区| 桃色一区二区三区在线观看| 天堂√8在线中文| 国产乱人伦免费视频| e午夜精品久久久久久久| 亚洲熟妇中文字幕五十中出| 国产精品99久久99久久久不卡| 国产亚洲欧美98| 亚洲精品粉嫩美女一区| 国产视频一区二区在线看| 又爽又黄无遮挡网站| av有码第一页| 午夜两性在线视频| 精品日产1卡2卡|