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

    北極中央?yún)^(qū)海冰密集度與云量相關(guān)性分析

    2015-01-05 06:45:20紀旭鵬趙進平
    海洋學(xué)報 2015年11期
    關(guān)鍵詞:低云高云云量

    紀旭鵬,趙進平

    (1 中國海洋大學(xué) 物理海洋教育部重點實驗室,山東 青島 266100)

    北極中央?yún)^(qū)海冰密集度與云量相關(guān)性分析

    紀旭鵬1,趙進平1

    (1 中國海洋大學(xué) 物理海洋教育部重點實驗室,山東 青島 266100)

    本文使用海冰密集度以及低云、中云、高云的日平均數(shù)據(jù),借助滑動相關(guān)分析方法,研究了北極中央?yún)^(qū)海冰密集度與云量之間的相關(guān)性,分析了海冰與云的相互作用機制。研究表明,在春季海冰融化季節(jié)(4、5月)、秋季海冰凍結(jié)季節(jié)(10、11月),低云與海冰密集度之間表現(xiàn)為較好的負相關(guān),表明在這段時間內(nèi)冰區(qū)海面蒸發(fā)強烈,對低云的形成有重要貢獻。在10月和11月,中云與海冰密集度也有很好的負相關(guān),表明秋季低云可以通過抬升形成中云。高云與海冰密集度之間并沒有明顯的相關(guān)性,可能原因:一方面海冰的空間分布對高云無影響,另一方面,高云主要影響到達的短波輻射,從而影響海冰的融化和凍結(jié)速度,與海冰厚度有直接顯著的關(guān)系,而與海冰密集度的關(guān)系不明顯。此外,在海冰密集度與低云存在較好負相關(guān)的情況下會出現(xiàn)某些年份相關(guān)性不好的情況,我們的研究發(fā)現(xiàn)這是北極中央?yún)^(qū)與周邊海區(qū)發(fā)生了海冰交換或云交換的結(jié)果。

    北極中央?yún)^(qū);海冰密集度;云量;滑動相關(guān)

    1 引言

    海冰作為北極重要的自然現(xiàn)象,強烈制約著海氣熱量交換,對北極氣候乃至北半球氣候起到了重要的調(diào)控作用,北極海冰通過影響大洋表面的輻射平衡、物質(zhì)平衡、能量平衡以及大洋溫、鹽流的形成和循環(huán)而影響全球氣候變化[1—3]。自20世紀70年代末開始,北極夏季海冰覆蓋范圍始終在不斷減少,厚度不斷變薄,多年冰的數(shù)量不斷變少,且最近幾年減少的速度明顯加快,整個北極地區(qū)一年冰及多年冰的總面積從1979-1996年的每10年減少2.2%和3.0%變?yōu)?998-2007年的每10年減少10.1%和10.7%[4]。在2012年9月16日,海冰覆蓋范圍出現(xiàn)了有歷史記錄以來的最低值[5]。

    云作為一種自然現(xiàn)象,對輻射有兩種近乎相反的重要影響:一方面,云的存在阻礙了太陽短波輻射到達北極地區(qū)表面,海冰吸收熱量減少,有利于海冰的維持;另一方面,云的存在阻礙了地球表面長波輻射的射出,造成氣溫升高,促進海冰融化[6]。大氣中云和霧的大幅度減少會導(dǎo)致海冰吸收熱量的大幅度增加,是導(dǎo)致北冰洋密集冰區(qū)的海冰快速融化的重要因素[7]。因此,北極地區(qū)的云在北極海冰變化中起到非常重要的作用。另外,北極地區(qū)海冰的增加或減少同樣可能對云有重要的影響。所以探究北極地區(qū)海冰與云之間的關(guān)系十分重要。

    云的高度實際上指示了云來源的差異。低云主要來自下墊面的蒸發(fā)過程,早期形態(tài)可能是海霧。在太陽輻射加熱條件下,表面蒸發(fā)的水汽發(fā)生抬升,形成低云。由于在北極有強烈的逆溫層,有很好的靜力穩(wěn)定度,不易產(chǎn)生強對流天氣,因而低云難以成為高云。高云的形成主要是在強對流區(qū)域發(fā)生的蒸發(fā)形成的,在北極周邊形成強對流的區(qū)域通??拷鶏u低壓和阿留申低壓,形成的高云隨大氣環(huán)流擴展到北極。

    前人已經(jīng)對北極地區(qū)海冰與云之間的相關(guān)性進行了研究。Schweiger等研究發(fā)現(xiàn),在海冰邊緣區(qū),云量變化與海冰變化有很強的關(guān)系,海冰減少,低云減少,同時中云增加[8]。Liu等的研究表明,海冰密集度每減少1%,云量會增加0.36%~0.47%[9]。海冰的厚度受到云量保溫作用的顯著影響,在春季,云量的改變影響到達冰面的短波輻射,是造成海冰厚度變化的主要驅(qū)動因子;但是,云量的變化對海冰密集度幾乎沒有影響[10—11]。在夏季,低云對海冰減少幾乎無響應(yīng),云量與海冰沒有表現(xiàn)出顯著的關(guān)系[12—13]。然而,2007年夏季北極海冰面積出現(xiàn)的極端低值,與云量減少造成的地面接收短波輻射增加有關(guān)[14]。在秋季,海冰邊緣區(qū)的開闊水域上方有云形成,但在海冰上方并沒有云形成[15]。在靠近太平洋的北極邊緣海中,新形成的開闊水域上方確實有低云量的增加[16],這可能促成一個初秋云-冰反饋機制的產(chǎn)生[13]。在冬季,云量的減少有助于海冰的增加[10]。Eastman指出,相比于中云和高云,低云是總云量的重要貢獻因子[16]。前人的研究表明,云量和海冰面積在春季和秋季有明顯的負相關(guān)性,夏季和冬季的相關(guān)性較弱。

    此外,數(shù)值模擬也被用于研究北極海冰與云之間的相關(guān)性。Beesley通過一個理想化能量平衡模式和單柱模式(SCM)研究北冰洋海冰厚度與云量的相關(guān)性時發(fā)現(xiàn),低云的增加有利于冰變厚,高云的增加將會導(dǎo)致冰變薄[17]。Vavrus等通過國際氣候系統(tǒng)模式CCSM3進行模擬,結(jié)果表明,海冰與中云及高云在各個季節(jié)都沒有明顯的相關(guān)關(guān)系,但是海冰減少與低云增加是相對應(yīng)的,只是兩者的超前滯后關(guān)系并不明確,并且夏季低云的減少可能會促進秋季海冰減少及云量增加[18]。Koenigk等利用全球耦合氣候模式EC-Earth2.3模擬的結(jié)果表明,在巴倫支海,海冰減少最為明顯,并且伴隨著海冰減少,低云大量減少,與此同時,中云與高云有少量增加[19]。

    前人的研究表明了海冰與云之間存在著季節(jié)性的相關(guān)性。不同研究者的結(jié)果表現(xiàn)出了較大的差異,可能與區(qū)域性差異有關(guān)。海冰與云量存在顯著的多年變化,人們對這種變化的認識不多。另外,以往的研究主要分析北極海冰邊緣區(qū)海冰密集度與云量的相關(guān)性,而北極海冰邊緣區(qū)的海冰變化會受到多種因素的影響。北極中央?yún)^(qū)主體由密集冰區(qū)覆蓋,與云的相關(guān)性與海冰邊緣區(qū)有很大不同,但隨著海冰的減少,北極中央?yún)^(qū)海冰逐漸也呈現(xiàn)出海冰邊緣區(qū)的特征。因此,就目前北極海冰的物理狀態(tài)來說,研究北極中央?yún)^(qū)海冰密集度和云量的相關(guān)關(guān)系顯得重要起來。為此,本文采用云量數(shù)據(jù)和海冰密集度數(shù)據(jù),研究北極中央?yún)^(qū)不同季節(jié)云與海冰密集度的相關(guān)性,研究云與海冰之間關(guān)系的多年變化,較全面地認識北極云量與北極中央?yún)^(qū)海冰密集度的關(guān)系。

    2 數(shù)據(jù)和分析方法

    2.1 數(shù)據(jù)

    本文使用的海冰密集度數(shù)據(jù)是85 GHz的SSM/I與/或SSM/IS通道測量的亮溫數(shù)據(jù)經(jīng)過ARTIST Sea Ice(ASI)算法得到的日平均數(shù)據(jù),并且在數(shù)據(jù)發(fā)布之前,整個時間序列都使用了中值濾波算法,空間分辨率為12.5 km×12.5 km。本文使用的云量數(shù)據(jù)是NOAA-CIRES 20th Century Reanalysis version 2的日平均數(shù)據(jù),分為3種,分別是低云數(shù)據(jù)、中云數(shù)據(jù)、高云數(shù)據(jù),該數(shù)據(jù)是空間分辨率為1.875°(經(jīng)度)×~1.904°(緯度)的經(jīng)緯度網(wǎng)格數(shù)據(jù),其中低云是地面至650 hPa之間的云,中云為400~650 hPa之間的云,高云是400 hPa至大氣層頂之間的云。海冰密集度數(shù)據(jù)與云量數(shù)據(jù)的時間跨度都為1992年1月1日至2010年12月31日。此外,為了剔除極端天氣現(xiàn)象的影響,我們對3種云量數(shù)據(jù)都進行了7點時間平均的平滑處理。

    同時,由于本文使用的云量數(shù)據(jù)是經(jīng)緯度網(wǎng)格數(shù)據(jù),在北半球,隨著緯度的增加,1°×1°范圍的面積在不斷變小,而且越靠近極地區(qū)域,面積變小的速率越大。所以我們需要引進一個面積系數(shù)Ac,將經(jīng)緯度網(wǎng)格數(shù)據(jù)轉(zhuǎn)化為面積數(shù)據(jù)。

    Ac=R2cosφΔλΔφ,

    (1)

    式中,R為地球半徑,λ為經(jīng)度,φ為緯度。每一個云量數(shù)據(jù)都乘以其對應(yīng)的面積系數(shù)便可以將云量的經(jīng)緯度網(wǎng)格數(shù)據(jù)轉(zhuǎn)化為面積數(shù)據(jù)。

    2.2 分析區(qū)域

    由于云具有很強的流動性,局地的云與海冰之間很難有顯著的相關(guān)。海冰密集度與云量如果有聯(lián)系,也只能在大范圍云和大范圍海冰之間存在相關(guān)性。因此,本文將北極地區(qū)分為北極太平洋扇區(qū)、北極大西洋扇區(qū)以及北極中央?yún)^(qū)(圖1)。確定北極中央?yún)^(qū)的形狀是比較任意的,主要是反映密集冰區(qū)的大致范圍,我們的試驗表明,形狀的改變對結(jié)果影響很小。本文重點分析北極中央?yún)^(qū)云與海冰的關(guān)系。

    圖1 北極地區(qū)分區(qū)圖Fig.1 Regional divisions of Arctic region

    2.3 滑動相關(guān)分析法

    本文借鑒Zhao等在研究北極濤動指數(shù)(AOI)與海表面壓強(SLP)之間相關(guān)性分析時使用的滑動相關(guān)方法[20],具體原理如下:

    i=1+n,…,N-n

    (2)

    式中,N為總的數(shù)據(jù)點數(shù),X代表海冰密集度數(shù)據(jù),Y代表云量數(shù)據(jù),Ri是第i個滑動窗口(-n,n)中心點對應(yīng)的海冰密集度與云量之間的滑動相關(guān)系數(shù)。經(jīng)過順序移動計算窗口,可以得到海冰密集度與云量數(shù)據(jù)在2n時間段相關(guān)系數(shù)的時間序列。通過分析滑動相關(guān)系數(shù)在各月以及各季節(jié)的高低,就能夠分析海冰密集度與云量之間的相關(guān)性隨時間的變化。

    滑動窗口n的選擇是滑動相關(guān)分析的關(guān)鍵,需要根據(jù)研究的物理量變化的時間尺度來確定。通常情況下,隨著n值的增加,得到的滑動相關(guān)系數(shù)會不斷變小。如果n值選擇太小,就會顯現(xiàn)更多高頻成分的影響;反之,如果n值選擇太大,就會突出低頻的成分,一些重要的較高頻率的現(xiàn)象將被消除。

    圖2顯示的是取不同滑動窗口得到的時間范圍為1993-1994年的北極中央?yún)^(qū)域海冰密集度與云量滑動相關(guān)系數(shù)繪制的曲線。n取15、30、45時分別代表1、2、3個月的滑動窗口。從圖中可以看出,n=15時,滑動相關(guān)系數(shù)的波動范圍太大,這是由天氣尺度的高頻擾動造成的。而n=45時,滑動相關(guān)系數(shù)與n=30時接近。因此,本文中進行北極中央?yún)^(qū)域海冰密集度與云量之間滑動相關(guān)分析時使用n=30,即兩個月的滑動窗口。需要說明的是,滑動窗口取兩個月,那么1月份的滑動相關(guān)系數(shù)應(yīng)該涉及12月至2月這一時間范圍內(nèi)的所有數(shù)據(jù),只是滑動窗口中心在1月份。

    圖2 不同時間窗口滑動相關(guān)系數(shù)對比圖Fig.2 Comparison of running correlation coefficients with different window lengths

    3 海冰密集度變化與云量的相關(guān)性

    本文在使用滑動相關(guān)方法進行北極中央?yún)^(qū)海冰密集度與云量之間相關(guān)性分析時,主要是分析兩者在各月以至各季節(jié)的相關(guān)性。以1月份為例,具體的相關(guān)性分析是將1992-2010年總共19年的海冰密集度與云量的日平均數(shù)據(jù)進行總體滑動相關(guān),得到滑動相關(guān)系數(shù)的19年時間序列,然后從這19年的滑動相關(guān)系數(shù)中,取出1月份的滑動相關(guān)系數(shù),繪制圖像,從圖像中就可以看出在這19年中,海冰密集度與云量在1月份所表現(xiàn)出的相關(guān)性。

    下面,我們對北極中央?yún)^(qū)海冰密集度與低云、高云、中云3種云在各季節(jié)的相關(guān)性進行分析。

    3.1 海冰密集度與低云的相關(guān)性

    春季的時間范圍為3-5月,在北極,海冰面積在3月達到最大值,4、5月份是北極海冰開始融化的季節(jié),但是在北極中央?yún)^(qū),這個時期海冰變化較小。從圖3中可以看到,北極中央?yún)^(qū)的海冰密集度與低云在不同年份的3月表現(xiàn)出的相關(guān)性差異很大,所以兩者并沒有較好的相關(guān)性;4月除了1994、2001、2004、2008年之外,海冰密集度與低云的滑動相關(guān)系數(shù)都為負值,并且大部分都通過了95%置信度檢驗;5月,除了在1995、1997、1999、2007、2008、2010年的部分時間以及2004年,海冰密集度與云量表現(xiàn)為較好的負相關(guān)。以上分析表明,4月和5月北極中央?yún)^(qū)的海冰密集度與低云表現(xiàn)為較好的負相關(guān)性。但是,有些年份相關(guān)性不好。

    6-8月是北冰洋的夏季,在北極中央?yún)^(qū),低云與海冰密集度在整個夏季都沒有表現(xiàn)出顯著的相關(guān)性。

    秋季的情況與夏季不同。9月在北極中央?yún)^(qū),海冰密集度與低云在不同年份表現(xiàn)出的相關(guān)性差異很大,表明兩者并沒有較好的相關(guān)性。10月除了在1996、1998、2000、2004年的小部分時間,以及2008-2010年海冰密集度與低云在10月表現(xiàn)為較好的負相關(guān)性。而在11月海冰密集度與低云的相關(guān)性最好,只有在1994、2008年的小部分時間以及2009-2010年相關(guān)性較差。10月和11月是北極地區(qū)的秋季海冰凍結(jié)季節(jié),上面的結(jié)果表明,在這段時間內(nèi),北極中央?yún)^(qū)的海冰密集度與低云表現(xiàn)為較好的負相關(guān)性。

    12月至翌年2月是北冰洋的冬季,雖然從圖中看出冬季海冰密集度與低云在整體上表現(xiàn)為負相關(guān)性,但是冬季海冰密集度很高,只在小范圍內(nèi)變化,冰間水道的比例很低,蒸發(fā)量小不足以形成大范圍的低云,不能較好地研究海冰密集度與低云的相關(guān)性,因此本文對冬季的情況不進行深入研究。

    低云與海面蒸發(fā)有明顯的關(guān)系,本文的結(jié)果表明,春季4、5月份和秋季10、11月份,低云與海冰密集度有較好的相關(guān),這與Eastman[16]的結(jié)果一致。這兩個季節(jié)的氣溫較低,開闊水道可以形成強烈的蒸發(fā),產(chǎn)生低云。負相關(guān)表明,海冰密集度越低云量越高,體現(xiàn)了蒸發(fā)對低云的貢獻。夏季海冰大幅度減少,而海冰密集度與低云的相關(guān)性卻很差,原因可能是夏季海冰大幅度減少,有更多的開闊水域,因而有更多的蒸發(fā)[21],但是由于夏季溫度高,飽和水汽壓高,蒸發(fā)的水汽并不會立即形成較多的低云,到了秋季,當(dāng)溫度降低時,更多的低云才會形成[22]。此外,也可能是由于夏季海冰密集度時間變化率比較大,云不存在同步的時間變化率,兩者即使有關(guān)聯(lián),也可能被海冰自身的時間變換率所掩蓋。冬季由于海冰密集度很高,氣溫低,開闊水道稀少,蒸發(fā)率小,不足以形成大范圍低云。

    3.2 海冰密集度與高云的相關(guān)性

    由于北極下表面溫度較低、逆溫層較強,沒有很強的對流形成高云,因而,北極的高云都不是當(dāng)?shù)睾Q笳舭l(fā)的結(jié)果,而是在對流區(qū)形成并被輸送到北冰洋。高云極大地影響到達海冰的太陽輻射,對于海冰的融化產(chǎn)生強烈的影響,云量與海冰之間應(yīng)該有顯著的正相關(guān)。然而,從圖4中我們可以看出,在北極中央?yún)^(qū),海冰密集度與高云之間在全年都沒有表現(xiàn)出較好的相關(guān)性,這與Vavrus等[18]的研究結(jié)果有很好的一致性。

    研究結(jié)果表明,海冰的空間分布不會影響高云,并且高云雖然通過影響到達的太陽輻射影響海冰融化過程,但其影響的主要是海冰厚度的變化,對海冰密集度的影響很小。

    3.3 海冰密集度與中云的相關(guān)性

    中云的高度介于低云與高云之間。中云的來源比較復(fù)雜,低云垂直對流上升可以轉(zhuǎn)化為中云,高云下沉可以轉(zhuǎn)化為中云,鄰近區(qū)域中云的平流輸入也是重要來源。在北極中央?yún)^(qū),4月和5月的海冰密集度與低云有很好的負相關(guān),但圖5表明,海冰密集度與中云之間并沒有表現(xiàn)出較好的相關(guān)性。

    10月,除了在1998、2000、2004、2008年的小部分時間以及1996年、2010年,海冰密集度與中云表現(xiàn)為較好的負相關(guān)性。同樣,在11月除了在2000、2004年的小部分時間以及1994、1996、2010年,海冰密集度與中云表現(xiàn)為較好的負相關(guān)性(圖5)。10月和11月的負相關(guān)特性與低云的情形相似,即海冰密集度越低,云量越高。因而,在秋季,北極中央?yún)^(qū)的中云來自于低云的對流上升或抬升。

    北極地區(qū)大氣普遍存在著逆溫層。逆溫層限制了垂向混合,增加了大氣穩(wěn)定性,在很大程度上限制了低云的對流上升過程,而逆溫層主要發(fā)生在1 km高度以下,2~3 km以上的大氣穩(wěn)定性并不是很強,并且隨著季節(jié)并不會發(fā)生太大變化[23—25]。Schweiger等[8]將秋季海冰與中云的負相關(guān)關(guān)系歸因于近地表溫度增加造成的大氣靜力穩(wěn)定度降低,進一步導(dǎo)致低云經(jīng)垂直對流轉(zhuǎn)化成為中云,而1997-1998年的SHEBA(Surface Heat Budget of the Arctic Ocean Experiment)計劃實測數(shù)據(jù)分析結(jié)果表明,秋冬季節(jié)的逆溫層強度要強于春夏季節(jié)[26],僅僅依據(jù)逆溫層強度無法解釋春季海冰密集度與中云相關(guān)性差。

    圖3 海冰密集度與低云的滑動相關(guān)圖Fig.3 Running correlation between sea ice concentration and low cloud圖中的紅線以及藍線表示95%置信度檢驗線,由于滑動窗口為2個月(61天),95%置信度檢驗的相關(guān)系數(shù)值為±0.254The red line and blue line indicate the lines of 95% confidence level,a 2-month(61-day)window is adopted,the critical value of the correlation coefficient for a 95% confidence level is ±0.254

    圖4 海冰密集度與高云的滑動相關(guān)圖Fig.4 Running correlation between sea ice concentration and high cloud圖中的紅線以及藍線表示95%置信度檢驗線,由于滑動窗口為2個月(61天),95%置信度檢驗的相關(guān)系數(shù)值為±0.254The red line and blue line indicate the lines of 95% confidence level,a 2-month(61-day)window is adopted,the critical value of the correlation coefficient for a 95% confidence level is ±0.254

    圖5 海冰密集度與中云的滑動相關(guān)圖Fig.5 Running correlation between sea ice concentration and medium cloud圖中的紅線以及藍線表示95%置信度檢驗線,由于滑動窗口為2個月(61天),95%置信度檢驗的相關(guān)系數(shù)值為±0.254The red line and blue line indicate the lines of 95% confidence level,a 2-month(61-day)window is adopted,the critical value of the correlation coefficient for a 95% confidence level is ±0.254

    我們認為,低云是否抬升與逆溫層的結(jié)構(gòu)有關(guān)。春夏以懸浮逆溫層為主,在低云的范圍內(nèi)存在很強的靜力穩(wěn)定度,不容易抬升到中云的高度,這是4、5月份的中云與海冰密集度相關(guān)性差的原因。秋冬季節(jié)以近地表逆溫層為主,低云通常發(fā)生在逆溫層之外,比較容易抬升[23—25]。此外,風(fēng)會破壞逆溫層、降低大氣穩(wěn)定度,而北極秋季平均風(fēng)速要比春季大,所以北極低云的垂直對流上升過程在秋季要強于春季[23]。

    為進一步明確北極中央?yún)^(qū)中云的來源及季節(jié)差異,我們分別對北極中央?yún)^(qū)的中云與低云,中云與高云以及中央?yún)^(qū)的中云與大西洋扇區(qū)和太平洋扇區(qū)的中云進行滑動相關(guān)分析(n=30),得到的平均相關(guān)系數(shù)如表1所示。

    表1 各區(qū)域云量與北極中央?yún)^(qū)中云的平均相關(guān)系數(shù)Tab.1 The average correlation coefficients between each regional cloud amount and the central Arctic medium cloud

    從表1中可以看出,北極中央?yún)^(qū)的中云與低云有非常好的相關(guān)性。北極中央?yún)^(qū)的中云與高云相關(guān)性很差,表明中云不是高云下沉的結(jié)果。同樣,北極中央?yún)^(qū)的中云與太平洋扇區(qū)和大西洋扇區(qū)的中云沒有顯著的相關(guān)性,并非是水平輸運的結(jié)果。因此,北極中央?yún)^(qū)的中云主要是北極中央?yún)^(qū)的低云經(jīng)對流上升轉(zhuǎn)化而來的,并且秋季要強于春季。這也解釋了北極中央?yún)^(qū)的中云與海冰密集度只有在秋季的10月和11月表現(xiàn)出較好的負相關(guān)關(guān)系,與低云具有相似的性質(zhì),而在春季沒有。

    3.4 個別年份相關(guān)性差異的原因分析

    從滑動相關(guān)曲線可以看出,海冰密集度與云量的滑動相關(guān)系數(shù)很好地體現(xiàn)了云與海冰之間的相關(guān)性,基本可以確定相關(guān)性較好的云類和月份。例如:低云在4、5、10、11月與海冰的滑動相關(guān)最好,反映了云與冰之間的密切物理聯(lián)系。

    然而,在相關(guān)性較好的時期也有相關(guān)度很差的時間段,包括:在1994、2001、2004、2008年的4月份;在1995、1997、1999、2004、2007、2008、2010年的5月份;在1996、1998、2000、2004、2008、2009、2010年的10月份;在1994、2008、2009、2010年的11月份。我們很希望知道,在這些年份究竟發(fā)生了什么而導(dǎo)致二者的相關(guān)性降低。只有找出了相關(guān)性差的原因,才能真正認識云量與海冰密集度的關(guān)系。

    如前所述,在相關(guān)性較好的年份,海冰密集度與云量呈現(xiàn)較高的負相關(guān)。從系統(tǒng)的角度看,北極中央?yún)^(qū)更像是一個封閉系統(tǒng),海冰與低云云量的相關(guān)性表現(xiàn)為系統(tǒng)內(nèi)部的耦合變化,即:海冰密集度降低,云量增大。某些年份兩者相關(guān)性降低,甚至呈現(xiàn)正相關(guān),表明系統(tǒng)的封閉性變差,在風(fēng)或流的作用下北極中央?yún)^(qū)云或者海冰與鄰近區(qū)域發(fā)生了交換。與鄰近區(qū)域的交換有3種類型,云交換、海冰交換、海冰和云同時交換。

    為了深入探究4、5、10、11月相關(guān)性較差的年份中,北極中央?yún)^(qū)究竟與哪一鄰近區(qū)域發(fā)生了海冰和云的輸入輸出交換,我們將北極中央?yún)^(qū)以南5°的范圍劃分了7個環(huán)區(qū),各個環(huán)區(qū)以所對應(yīng)的邊緣海命名,如圖6所示。

    另外,如果北極中央?yún)^(qū)與某一環(huán)區(qū)發(fā)生了海冰(或云)的交換,那么這兩個區(qū)域之間的海冰密集度(或云量)有可能呈反向變化,相關(guān)性表現(xiàn)為負相關(guān),因此我們計算了北極中央?yún)^(qū)與各個環(huán)區(qū)之間海冰密集度(或云量)的滑動相關(guān)系數(shù)(滑動窗口為2個月),并且給出了兩個區(qū)域之間的海冰密集度(或云量)變化趨勢,從而探究了4、5、10、11月相關(guān)性較差的年份中,北極中央?yún)^(qū)究竟與哪一環(huán)區(qū)發(fā)生了海冰和云的交換。在此以2009年(圖7)和2010年(圖8)的11月份為例說明。

    圖6 北極地區(qū)細化分區(qū)圖Fig.6 Detailed regional divisions of Arctic region(1)格陵蘭海環(huán)區(qū),(2)巴倫支海環(huán)區(qū),(3)喀拉海環(huán)區(qū),(4)拉普捷夫海環(huán)區(qū),(5)東西伯利亞海環(huán)區(qū),(6)楚科奇海環(huán)區(qū),(7)波弗特海環(huán)區(qū)(1)Greenland Sea region,(2)Barents Sea region,(3)Kara Sea region,(4)Laptev Sea region,(5)East Siberian Sea region,(6)Chukchi Sea region,(7)Beaufort Sea region

    圖7 2009年11月北極中央?yún)^(qū)與巴倫支海環(huán)區(qū)之間海冰交換分析圖Fig.7 Sea ice exchange between the central Arctic and Barents Sea region during November 2009a.2009年11月北極中央?yún)^(qū)海冰密集度與低云云量滑動相關(guān)圖;b.2009年11月北極中央?yún)^(qū)與巴倫支海環(huán)區(qū)海冰密集度滑動相關(guān)圖;c.2009年11月份北極中央?yún)^(qū)與巴倫支海環(huán)區(qū)海冰密集度變化趨勢對比圖a.Running correlation between sea ice concentration and low cloud amount in the central Arctic during November 2009;b.running correlation of sea ice concentration between the central Arctic and Barents Sea region during November 2009;c.comparison of sea ice concentration variation trend between the central Arctic and Barents Sea region during November 2009

    圖8 2010年11月北極中央?yún)^(qū)與拉普捷夫海環(huán)區(qū)之間低云交換分析圖Fig.8 Low cloud exchange between the central Arctic and Barents Sea region during November 2010a.2010年11月北極中央?yún)^(qū)海冰密集度與低云云量滑動相關(guān)圖;b.2010年11月北極中央?yún)^(qū)與拉普捷夫海環(huán)區(qū)低云云量滑動相關(guān)圖;c.2010年11月北極中央?yún)^(qū)與拉普捷夫海環(huán)區(qū)低云云量變化趨勢對比圖a.Running correlation between sea ice concentration and low cloud amount in the central Arctic during November 2010;b.running correlation of low cloud amount between the central Arctic and Laptev Sea region during November 2010;c.comparison of low cloud amount variation trend between the central Arctic and Laptev Sea region during November 2010

    由圖7可見,2009年11月北極中央?yún)^(qū)與巴倫支海環(huán)區(qū)的海冰密集度呈負相關(guān),由兩個區(qū)域之間海冰密集度變化趨勢的對比看出,兩個區(qū)域的海冰密集度發(fā)生了反向變化,表明兩個區(qū)域可能發(fā)生了海冰交換。圖8表明,2010年11月份,與北極中央?yún)^(qū)進行交換的區(qū)域是拉普捷夫海,交換類型是云交換。需要說明的是:第一,北極中央?yún)^(qū)與相鄰區(qū)域的反向變化僅僅是可能發(fā)生交換的條件,其他非交換因素也可以導(dǎo)致反向變化;第二,如果發(fā)生交換導(dǎo)致二者反向變化可以用這種方法加以判別,有些情況下發(fā)生交換卻并不能產(chǎn)生反向變化,用這種方法就無法判別。雖然這種方法具有這些局限性,但其結(jié)果為我們理解云量與海冰密集度的相關(guān)性提供了方法。

    使用該方法對所有相關(guān)性差的時間段進行分析,匯總結(jié)果如表2所示。

    由表2~5可以看出,4月份,全部是海冰交換,5、10、11月份,海冰交換和云交換都有,整體上以海冰交換為主,未見同一環(huán)區(qū)發(fā)生海冰和云同時交換的情況。與北極中央?yún)^(qū)發(fā)生交換的區(qū)域以格陵蘭海環(huán)區(qū)和東西伯利亞海環(huán)區(qū)為主。此外,2007和2008年的5月份、2008年的10月份以及1994年的11月份,這4個時段較差的相關(guān)性無法用北極中央?yún)^(qū)與環(huán)區(qū)之間發(fā)生海冰交換或云交換進行解釋。我們的研究結(jié)果表明,低云與海冰密集度相關(guān)性較差的時段中,有81.8%的時段可以用北極中央?yún)^(qū)與環(huán)區(qū)之間發(fā)生海冰交換或云交換進行解釋,表明在風(fēng)或流的作用下北極中央?yún)^(qū)與鄰近區(qū)域之間云及海冰的輸入輸出交換是最主要的原因。

    表2 4月份與北極中央?yún)^(qū)進行海冰和云交換的各環(huán)區(qū)匯總表Tab.2 The summary of each region with which the central Arctic has sea ice and cloud exchange during April

    表3 5月份與北極中央?yún)^(qū)進行海冰和云交換的各環(huán)區(qū)匯總表Tab.3 The summary of each region with which the central Arctic has sea ice and cloud exchange during May

    表4 10月份與北極中央?yún)^(qū)進行海冰和云交換的各環(huán)區(qū)匯總表Tab.4 The summary of each region with which the central Arctic has sea ice and cloud exchange during October

    表5 11月份與北極中央?yún)^(qū)進行海冰和云交換的各環(huán)區(qū)匯總表Tab.5 The summary of each region with which the central Arctic has sea ice and cloud exchange during November

    4 結(jié)果與討論

    本文將北極中央?yún)^(qū)的海冰密集度和低云、中云和高云云量作整體平均,得到二者的時間序列。為了體現(xiàn)二者相關(guān)程度的變化,我們計算了滑動相關(guān)系數(shù)?;瑒酉嚓P(guān)的主要價值是能夠體現(xiàn)那些有時相關(guān)好,有時相關(guān)不好的現(xiàn)象。由于海冰的密集度不僅取決于云,還與多種因素有關(guān),因此,預(yù)期海冰密集度與云量的相關(guān)系數(shù)不會很高。但是,本文的結(jié)果表明,海冰密集度與云量之間還是有比較顯著的聯(lián)系。

    結(jié)果表明,春季4、5月份海冰密集度與低云相關(guān)度很高,主要是春季冰間水道強烈的蒸發(fā)形成低云;由于春季低云的垂直對流上升過程較弱,海冰密集度對中云的貢獻不大。秋季10月和11月海冰密集度與低云有很高的相關(guān)度,也是由于氣溫逐漸降低過程中冰間水域的強烈蒸發(fā)形成低云。與春季不同的是,10月和11月的海冰密集度與中云也有很高的相關(guān)度,意味著秋季的對流上升過程強烈,低云可以上升到中云的高度。冬季海冰密集度接近1,與云的相關(guān)性應(yīng)該較差。而夏季雖然蒸發(fā)強烈,但是飽和水汽壓高,難以形成低云,此外,海冰密集度時間變化率較大,會將它與云量可能存在的關(guān)聯(lián)信號掩蓋。春季和秋季的相關(guān)性主要體現(xiàn)了海冰對云形成的影響。而高云對海冰密集度的影響從二者的相關(guān)度來看并不明顯,表明海冰的空間分布不會影響高云,高云主要體現(xiàn)為對海冰厚度變化的影響,而不是影響海冰密集度的變化。

    在海冰密集度與低云有較好相關(guān)性的4、5、10、11月份中,某些年份出現(xiàn)相關(guān)性不好的情況。我們對此進行了詳細的分析。在正常情況下,北極中央?yún)^(qū)在大多數(shù)年份體現(xiàn)為相對封閉的特性,海冰密集度與低云呈現(xiàn)較好的負相關(guān)。而在某些相關(guān)性較差或呈正相關(guān)的年份,北極中央?yún)^(qū)與周邊的環(huán)區(qū)進行了海冰或云的交換,交換的區(qū)域主要以格陵蘭海環(huán)區(qū)以及東西伯利亞海環(huán)區(qū)為主。北極中央?yún)^(qū)與周邊海域發(fā)生冰或云的交換可以對81.8%的情形進行解釋,表明北極中央?yún)^(qū)之外的海冰和低云會通過海流和風(fēng)場與北極中央?yún)^(qū)進行交換。

    我們認為,北極中央?yún)^(qū)海冰與云的關(guān)系有兩種。一種是云削弱了太陽輻射,減緩了海冰的融化;另一種是海冰減少,蒸發(fā)增大,形成低云。如果云是冰間水道產(chǎn)生的,海冰密集度與云量之間應(yīng)該是負相關(guān),即海冰密集度越低云量應(yīng)該越大;如果是云影響了海冰,二者應(yīng)該是正相關(guān),即云量越大,海冰密集度應(yīng)該越高。在我們的研究結(jié)果中,海冰密集度與云量之間沒有出現(xiàn)顯著的正相關(guān),表明在北極中央?yún)^(qū),云雖然改變了太陽輻射,但海冰的變化主要體現(xiàn)在海冰融化加劇導(dǎo)致的厚度的變化,海冰密集度與云量聯(lián)系不大。而在海冰開裂、密集度減少的春季和秋季,由于海洋蒸發(fā)強烈,容易形成低云,因此,海冰密集度與低云有密切的關(guān)系。夏季的情況不同,夏季雖然海冰密集度最低,但由于夏季海面氣溫高,飽和水汽壓高,不利于低云形成,并且海冰密集度的時間變化率較大,所以海冰密集度與低云相關(guān)性差。

    [1] Morison J H,McPhee M G,Curtin T B,et al. The oceanography of winter leads[J]. Journal of Geophysical Research: Oceans (1978-2012),1992,97(C7): 11199-11218.

    [2] Alam A,Curry J. Lead-induced atmospheric circulations[J]. Journal of Geophysical Research: Oceans (1978-2012),1995,100(C3): 4643-4651.

    [3] 張樹剛. 夏季北冰洋海冰融化過程中能量分配的研究[D]. 青島: 中國海洋大學(xué),2009.

    Zhang Shugang. Study on the distribution of solar energy in marginal ice zone during ice melting season[D]. Qingdao: Ocean University of China,2009.

    [4] Comiso J C,Parkinson C L,Gersten R,et al. Accelerated decline in the Arctic sea ice cover[J]. Geophysical Research Letters,2008,35(1): L01703.

    [5] Zhang J,Lindsay R,Schweiger A,et al. The impact of an intense summer cyclone on 2012 Arctic sea ice retreat[J]. Geophysical Research Letters,2013,40(4): 720-726.

    [6] Gorodetskaya I V,Tremblay L. Arctic cloud properties and radiative forcing from observations and their role in sea ice decline predicted by the NCAR CCSM3 model during the 21st century[J]. Arctic Sea Ice Decline: Observations,Projections,Mechanisms,and Implications,2008,180: 47-62.

    [7] 趙進平,李濤,張樹剛,等. 北冰洋中央密集冰區(qū)海冰對太陽短波輻射能吸收的觀測研究[J]. 地球科學(xué)進展,2009,24(1): 33-41.

    Zhao Jinping,Li Tao,Zhang Shugang,et al. The shortwave solar radiation energy absorbed by packed sea ice in the central Arctic[J]. Advances in Earth Science,2009,24(1): 33-41.

    [8] Schweiger A J,Lindsay R W,Vavrus S,et al. Relationships between Arctic sea ice and clouds during autumn[J]. Journal of Climate,2008,21(18): 4799-4810.

    [9] Liu Y,Key J R,Liu Z,et al. A cloudier Arctic expected with diminishing sea ice[J]. Geophysical Research Letters,2012,39(5): L05705.

    [10] Schweiger A J. Changes in seasonal cloud cover over the Arctic seas from satellite and surface observations[J]. Geophysical Research Letters,2004,31(12): L12207.

    [11] Francis J A,Hunter E. New insight into the disappearing Arctic sea ice[J]. Eos,Transactions American Geophysical Union,2006,87(46): 509-511.

    [12] Herman G,Goody R. Formation and persistence of summertime Arctic stratus clouds[J]. Journal of the Atmospheric Sciences,1976,33(8): 1537-1553.

    [13] Kay J E,Gettelman A. Cloud influence on and response to seasonal Arctic sea ice loss[J]. Journal of Geophysical Research: Atmospheres (1984-2012),2009,114(D18):D18204.

    [14] Kay J E,L’Ecuyer T,Gettelman A,et al. The contribution of cloud and radiation anomalies to the 2007 Arctic sea ice extent minimum[J]. Geophysical Research Letters,2008,35(8): L08503.

    [15] Paluch I R,Lenschow D H,Wang Q. Arctic boundary layer in the fall season over open and frozen sea[J]. Journal of Geophysical Research: Atmospheres (1984-2012),1997,102(D22): 25955-25971.

    [16] Eastman R. Interannual variations of Arctic cloud types in relation to sea ice[D]. Seattle: University of Washington,2009.

    [17] Beesley J A. Estimating the effect of clouds on the arctic surface energy budget[J]. Journal of Geophysical Research: Atmospheres (1984-2012),2000,105(D8): 10103-10117.

    [18] Vavrus S,Holland M M,Bailey D A. Changes in Arctic clouds during intervals of rapid sea ice loss[J]. Climate Dynamics,2011,36(7/8): 1475-1489.

    [19] Koenigk T,Brodeau L,Graversen R G,et al. Arctic climate change in 21st century CMIP5 simulations with EC-Earth[J]. Climate dynamics,2013,40(11/12): 2719-2743.

    [20] Zhao J,Cao Y,Shi J. Core region of Arctic Oscillation and the main atmospheric events impact on the Arctic[J]. Geophysical research letters,2006,33(22): L22708.

    [21] Perovich D K,Light B,Eicken H,et al. Increasing solar heating of the Arctic Ocean and adjacent seas,1979-2005: Attribution and role in the ice-albedo feedback[J]. Geophysical Research Letters,2007,34(19):L19505.

    [22] Wu D L,Lee J N. Arctic low cloud changes as observed by MISR and CALIOP: Implication for the enhanced autumnal warming and sea ice loss[J]. Journal of Geophysical Research: Atmospheres (1984-2012),2012,117(D7):D07107.

    [23] Serreze M C,Barrett A P,Stroeve J C,et al. The emergence of surface-based Arctic amplification[J]. The Cryosphere,2009,3(1): 11-19.

    [24] Pavelsky T M,Boé J,Hall A,et al. Atmospheric inversion strength over polar oceans in winter regulated by sea ice[J]. Climate Dynamics,2011,36(5/6): 945-955.

    [25] 程剛,高志球,鄭友飛,等.北極低空急流和低層逆溫特征觀測分析[J].氣候與環(huán)境研究,2013,18(1): 23-31.

    Cheng Gang,Gao Zhiqiu,Zheng Youfei,et al. A study on low-level jets and temperature inversion over the Arctic ocean by using SHEBA data[J]. Climatic and Environmental Research,2013,18(1): 23-31.

    [26] Tjernstr?m M,Graversen R G. The vertical structure of the lower Arctic troposphere analysed from observations and the ERA-40 reanalysis[J]. Quarterly Journal of the Royal Meteorological Society,2009,135(639): 431-443.

    Analysis of correlation between sea ice concentration and cloudiness in the central Arctic

    Ji Xupeng1,Zhao Jinping1

    (1.KeyLaboratoryofPhysicalOceanography,MinistryofEducation,OceanUniversityofChina,Qingdao266100,China)

    Correlation between sea ice concentration and cloudiness in the central Arctic is studied. The running correlation coefficients of daily sea ice concentration averaged for whole central Arctic with the daily averaged low cloud,medium cloud and high cloud are calculated. During the melting period in spring (April and May) and the freezing period in autumn (October and November),sea ice concentration and low cloud present significant negative correlation. This could be explained that the low clouds were formed by strong evaporation from the open water in gaps of sea ice during these periods. In the same period,sea ice concentration is negatively correlated with medium cloud only during the freezing period in autumn,which suggests that in autumn the low cloud could upwell to form medium cloud,but in spring the upwell is much week because of its stable stratification. There is no a significant correlation between sea ice concentration and high cloud,which is speculated that sea ice can’t influence the high cloud and the high cloud is not formed by local eveporation and the influence of high cloud on sea ice is not ice concentration,but is ice thickness. In the periods of significant correlation between ice concentration and low cloud,there still exist some occasional inconsistent betweeen them. The exchanges of sea ice or cloud between the central Arctic and adjacent regions is suggested to cause the inconsistent.

    central Arctic Ocean; sea ice concentration; cloudiness; running correlation

    2015-04-16;

    2015-06-16。

    國家重大科學(xué)研究計劃項目(2015CB953903);國家自然科學(xué)委員會重點基金項目(41330960)。

    紀旭鵬(1990—),男,山東省萊陽市人,從事極地海洋學(xué)與全球海洋變化研究。E-mail:jixupeng90@126.com

    10.3969/j.issn.0253-4193.2015.11.009

    P732.6

    A

    0253-4193(2015)11-0092-13

    紀旭鵬,趙進平. 北極中央?yún)^(qū)海冰密集度與云量相關(guān)性分析[J]. 海洋學(xué)報,2015,37(11): 92-104,

    Ji Xupeng,Zhao Jinping. Analysis of correlation between sea ice concentration and cloudiness in the central Arctic[J]. Haiyang Xuebao,2015,37(11): 92-104,doi:10.3969/j.issn.0253-4193.2015.11.009

    猜你喜歡
    低云高云云量
    2007—2020年威海機場4—9月平流低云的變化特征及影響因子分析
    晉江機場低云天氣統(tǒng)計分析
    科技資訊(2023年21期)2023-11-22 08:35:46
    贛州地區(qū)云量變化特征及其與降水的關(guān)系
    ASC200型地基雙波段全天空云量自動觀測儀云量比對分析
    氣象科技(2022年2期)2022-04-28 09:35:46
    走進新疆
    國畫家(2022年2期)2022-04-13 09:07:46
    遺愛湖畔看雨
    人生之秋也當(dāng)天高云闊
    1971—2010年虎林市云量與氣溫、降水的年際變化特征分析
    迷宮闖一闖
    Our Lovely School
    免费不卡的大黄色大毛片视频在线观看| 超碰97精品在线观看| 少妇被粗大的猛进出69影院 | 色哟哟·www| 国产永久视频网站| 亚洲精华国产精华液的使用体验| 99热6这里只有精品| 丰满迷人的少妇在线观看| 3wmmmm亚洲av在线观看| 少妇裸体淫交视频免费看高清| 国产淫语在线视频| 另类精品久久| 精品卡一卡二卡四卡免费| 日韩强制内射视频| 3wmmmm亚洲av在线观看| 男女啪啪激烈高潮av片| 亚洲人成网站在线观看播放| 一级黄片播放器| 人妻人人澡人人爽人人| 国产视频首页在线观看| 国产成人91sexporn| 国产亚洲5aaaaa淫片| 亚洲电影在线观看av| 老女人水多毛片| 丝袜在线中文字幕| 国产精品99久久99久久久不卡 | 天天操日日干夜夜撸| 国产又色又爽无遮挡免| 国产高清三级在线| a 毛片基地| av卡一久久| 国产亚洲午夜精品一区二区久久| 国产午夜精品久久久久久一区二区三区| 国产欧美日韩精品一区二区| 美女cb高潮喷水在线观看| 五月玫瑰六月丁香| 91久久精品电影网| 精品少妇内射三级| 又大又黄又爽视频免费| 成年美女黄网站色视频大全免费 | 在线观看国产h片| 婷婷色综合大香蕉| 永久网站在线| 亚洲国产日韩一区二区| 国产一级毛片在线| 好男人视频免费观看在线| 9色porny在线观看| 欧美日韩精品成人综合77777| 久久久久久久久久成人| 午夜福利影视在线免费观看| 亚洲av福利一区| 久久久午夜欧美精品| 成人亚洲欧美一区二区av| 亚洲av成人精品一区久久| 亚洲人与动物交配视频| 中文字幕精品免费在线观看视频 | 伊人久久国产一区二区| 国产成人91sexporn| 日韩,欧美,国产一区二区三区| 老司机影院成人| 777米奇影视久久| 亚洲欧美精品专区久久| 九九在线视频观看精品| 亚洲性久久影院| 观看免费一级毛片| 女的被弄到高潮叫床怎么办| 一本大道久久a久久精品| 91久久精品国产一区二区三区| 国产成人精品久久久久久| 久久毛片免费看一区二区三区| 一级av片app| 免费大片18禁| 最近最新中文字幕免费大全7| 老司机影院毛片| 国产免费又黄又爽又色| 亚洲av免费高清在线观看| 久久97久久精品| 狂野欧美激情性xxxx在线观看| 一级,二级,三级黄色视频| 少妇高潮的动态图| av黄色大香蕉| 日本黄色片子视频| 男男h啪啪无遮挡| 国产综合精华液| 日韩大片免费观看网站| 日本午夜av视频| 欧美xxⅹ黑人| 亚洲成人av在线免费| 欧美成人精品欧美一级黄| 亚洲欧美清纯卡通| 欧美少妇被猛烈插入视频| 菩萨蛮人人尽说江南好唐韦庄| 伊人久久精品亚洲午夜| 色哟哟·www| 亚洲欧美成人综合另类久久久| 欧美日韩一区二区视频在线观看视频在线| 一级毛片aaaaaa免费看小| 久热久热在线精品观看| a 毛片基地| 免费观看性生交大片5| 成人国产麻豆网| 两个人的视频大全免费| 看十八女毛片水多多多| 精品一区二区三卡| 亚洲精品第二区| 国模一区二区三区四区视频| 丰满乱子伦码专区| 丰满人妻一区二区三区视频av| 国产伦在线观看视频一区| av视频免费观看在线观看| 精品亚洲成国产av| 香蕉精品网在线| 久久国产亚洲av麻豆专区| 国产一区亚洲一区在线观看| freevideosex欧美| 久久精品久久久久久噜噜老黄| 日日摸夜夜添夜夜添av毛片| 精品亚洲乱码少妇综合久久| 成人美女网站在线观看视频| 亚洲三级黄色毛片| 久久久久国产网址| 夫妻性生交免费视频一级片| 三级国产精品欧美在线观看| 久久人妻熟女aⅴ| 亚洲久久久国产精品| 日韩欧美一区视频在线观看 | 男女啪啪激烈高潮av片| 亚洲成色77777| 亚洲,欧美,日韩| 欧美变态另类bdsm刘玥| 乱人伦中国视频| 精品人妻熟女av久视频| 美女cb高潮喷水在线观看| 99久久综合免费| 亚洲,一卡二卡三卡| 高清不卡的av网站| 免费av中文字幕在线| 91久久精品国产一区二区三区| 看非洲黑人一级黄片| 99热国产这里只有精品6| 国产欧美日韩综合在线一区二区 | 成人午夜精彩视频在线观看| 国产淫语在线视频| 男女无遮挡免费网站观看| 久久精品国产亚洲av天美| 日韩免费高清中文字幕av| 热re99久久国产66热| 精品卡一卡二卡四卡免费| 午夜91福利影院| 美女内射精品一级片tv| 久久韩国三级中文字幕| 亚洲欧洲精品一区二区精品久久久 | 国产av码专区亚洲av| 成人综合一区亚洲| 日本黄色片子视频| 国产精品一区www在线观看| 啦啦啦视频在线资源免费观看| 精品一区二区免费观看| 欧美日韩精品成人综合77777| 国内少妇人妻偷人精品xxx网站| 成年av动漫网址| 国产美女午夜福利| 人妻人人澡人人爽人人| 久久韩国三级中文字幕| 色视频在线一区二区三区| 成人影院久久| av免费在线看不卡| 午夜av观看不卡| freevideosex欧美| 中文天堂在线官网| 国产av码专区亚洲av| 免费看光身美女| 久久毛片免费看一区二区三区| 韩国av在线不卡| 男人爽女人下面视频在线观看| 在线观看美女被高潮喷水网站| 大码成人一级视频| 美女大奶头黄色视频| 最新的欧美精品一区二区| 国产精品成人在线| 久久久亚洲精品成人影院| 久久免费观看电影| 亚洲真实伦在线观看| 欧美日韩在线观看h| 亚洲精品日韩av片在线观看| 免费高清在线观看视频在线观看| 免费观看在线日韩| 免费观看在线日韩| 国产免费福利视频在线观看| 边亲边吃奶的免费视频| 午夜日本视频在线| 国产真实伦视频高清在线观看| 亚洲国产精品一区二区三区在线| 日本黄色日本黄色录像| h日本视频在线播放| 伊人久久精品亚洲午夜| 国产黄片美女视频| 女人久久www免费人成看片| 中文字幕人妻丝袜制服| av专区在线播放| 少妇的逼水好多| 亚洲精华国产精华液的使用体验| 三级国产精品片| 9色porny在线观看| av在线播放精品| 黄色视频在线播放观看不卡| 国产亚洲一区二区精品| 亚洲精品国产av成人精品| 色婷婷久久久亚洲欧美| 青春草国产在线视频| 国产免费又黄又爽又色| 亚洲欧美日韩东京热| 在线观看av片永久免费下载| 菩萨蛮人人尽说江南好唐韦庄| 一级片'在线观看视频| 91精品国产国语对白视频| 欧美日韩av久久| 国产精品麻豆人妻色哟哟久久| 国产在线男女| 成人免费观看视频高清| 国产 一区精品| 久久人人爽人人爽人人片va| 一二三四中文在线观看免费高清| 亚洲综合精品二区| 亚洲欧美日韩另类电影网站| 欧美一级a爱片免费观看看| 少妇高潮的动态图| 午夜av观看不卡| 免费av中文字幕在线| 大香蕉97超碰在线| 黑人巨大精品欧美一区二区蜜桃 | av网站免费在线观看视频| 亚洲国产av新网站| 视频中文字幕在线观看| 午夜老司机福利剧场| 国产视频内射| 精品亚洲成a人片在线观看| 久久精品国产亚洲av涩爱| 欧美精品人与动牲交sv欧美| 爱豆传媒免费全集在线观看| 日日爽夜夜爽网站| 国产69精品久久久久777片| 成年女人在线观看亚洲视频| 肉色欧美久久久久久久蜜桃| 亚洲欧美成人综合另类久久久| 丰满人妻一区二区三区视频av| 国国产精品蜜臀av免费| 99久国产av精品国产电影| 香蕉精品网在线| 欧美人与善性xxx| 汤姆久久久久久久影院中文字幕| 日本wwww免费看| 午夜免费鲁丝| 欧美丝袜亚洲另类| 国产国拍精品亚洲av在线观看| 中国国产av一级| 婷婷色麻豆天堂久久| 在线观看av片永久免费下载| 啦啦啦中文免费视频观看日本| 日本黄色片子视频| 国产一区二区三区综合在线观看 | 国产免费一级a男人的天堂| 少妇的逼水好多| 最新中文字幕久久久久| 国产69精品久久久久777片| 成年美女黄网站色视频大全免费 | 国产精品国产三级国产专区5o| 日韩电影二区| 成人国产麻豆网| 国产精品99久久99久久久不卡 | 又大又黄又爽视频免费| 国产精品一二三区在线看| 九草在线视频观看| 成人漫画全彩无遮挡| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆| 热99国产精品久久久久久7| 亚洲,欧美,日韩| 自拍欧美九色日韩亚洲蝌蚪91 | 五月玫瑰六月丁香| 老熟女久久久| 久久精品国产自在天天线| 欧美xxxx性猛交bbbb| 九九在线视频观看精品| 91精品伊人久久大香线蕉| 成年女人在线观看亚洲视频| 日韩大片免费观看网站| 蜜桃久久精品国产亚洲av| 日日啪夜夜撸| 下体分泌物呈黄色| 永久网站在线| 午夜福利影视在线免费观看| 性高湖久久久久久久久免费观看| 亚洲精品国产色婷婷电影| 如日韩欧美国产精品一区二区三区 | 男的添女的下面高潮视频| 99久久综合免费| 另类亚洲欧美激情| 少妇的逼好多水| 熟女人妻精品中文字幕| 欧美激情极品国产一区二区三区 | 人体艺术视频欧美日本| 久久 成人 亚洲| 久久婷婷青草| 亚洲自偷自拍三级| 亚洲av在线观看美女高潮| 在线观看免费日韩欧美大片 | 中文字幕制服av| 国产永久视频网站| 最近手机中文字幕大全| 亚洲成人av在线免费| 校园人妻丝袜中文字幕| 成人午夜精彩视频在线观看| 国产精品人妻久久久影院| 亚洲精品456在线播放app| 国产成人精品无人区| 黄色怎么调成土黄色| 日本欧美视频一区| av又黄又爽大尺度在线免费看| 性色av一级| 男的添女的下面高潮视频| 大香蕉97超碰在线| 国产白丝娇喘喷水9色精品| 亚洲国产精品一区三区| 97超碰精品成人国产| 免费久久久久久久精品成人欧美视频 | 日韩伦理黄色片| 久久99热6这里只有精品| 两个人的视频大全免费| 在线 av 中文字幕| 国产成人精品久久久久久| 王馨瑶露胸无遮挡在线观看| 人人妻人人澡人人爽人人夜夜| 九色成人免费人妻av| 乱码一卡2卡4卡精品| 成人国产麻豆网| 欧美成人精品欧美一级黄| 日日啪夜夜撸| 建设人人有责人人尽责人人享有的| 亚洲欧美精品自产自拍| 日本欧美国产在线视频| 午夜福利在线观看免费完整高清在| 久久韩国三级中文字幕| 女人久久www免费人成看片| 综合色丁香网| 久久精品国产亚洲av涩爱| 伦理电影免费视频| 久久毛片免费看一区二区三区| 成人国产av品久久久| 久久久a久久爽久久v久久| 国产成人精品无人区| 一级片'在线观看视频| 大片免费播放器 马上看| 人妻系列 视频| 中文资源天堂在线| 最黄视频免费看| 一级毛片我不卡| 久热这里只有精品99| 中文欧美无线码| 午夜91福利影院| 日韩精品免费视频一区二区三区 | 亚洲av电影在线观看一区二区三区| 欧美日韩av久久| 色5月婷婷丁香| 久久久久精品久久久久真实原创| 精品卡一卡二卡四卡免费| 成人美女网站在线观看视频| 男女啪啪激烈高潮av片| 久久久久久久精品精品| 日本91视频免费播放| freevideosex欧美| 日本av手机在线免费观看| 成人二区视频| 亚洲成色77777| 国产亚洲91精品色在线| 午夜老司机福利剧场| 99热这里只有是精品在线观看| 狠狠精品人妻久久久久久综合| 18+在线观看网站| 亚洲av中文av极速乱| 黄色毛片三级朝国网站 | 伊人亚洲综合成人网| 九色成人免费人妻av| 国产淫片久久久久久久久| 日韩电影二区| 街头女战士在线观看网站| 成人18禁高潮啪啪吃奶动态图 | 欧美3d第一页| 美女主播在线视频| 国产69精品久久久久777片| 99视频精品全部免费 在线| 大码成人一级视频| 国产视频内射| 亚洲成人av在线免费| 97在线视频观看| 99热这里只有精品一区| 久久女婷五月综合色啪小说| 亚洲精品色激情综合| 欧美 亚洲 国产 日韩一| 乱人伦中国视频| 中文欧美无线码| 丝袜在线中文字幕| 日韩精品有码人妻一区| 老熟女久久久| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| 国产视频首页在线观看| 高清视频免费观看一区二区| 少妇裸体淫交视频免费看高清| 一级毛片aaaaaa免费看小| 九色成人免费人妻av| 人妻 亚洲 视频| 久久久精品94久久精品| 国产精品99久久久久久久久| 国产成人a∨麻豆精品| 欧美日韩综合久久久久久| 一区二区三区乱码不卡18| 插逼视频在线观看| 女性被躁到高潮视频| 久久免费观看电影| 婷婷色av中文字幕| 又大又黄又爽视频免费| tube8黄色片| 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 亚洲欧洲精品一区二区精品久久久 | 九九爱精品视频在线观看| 亚洲怡红院男人天堂| 丝袜喷水一区| 少妇的逼好多水| 久久久欧美国产精品| 观看美女的网站| 国产欧美亚洲国产| 午夜福利影视在线免费观看| 22中文网久久字幕| 国产中年淑女户外野战色| 国产黄频视频在线观看| 纯流量卡能插随身wifi吗| 一区二区av电影网| 男人添女人高潮全过程视频| 国产国拍精品亚洲av在线观看| 日本91视频免费播放| 2018国产大陆天天弄谢| 色5月婷婷丁香| 久久97久久精品| 色网站视频免费| 国产精品熟女久久久久浪| 一级,二级,三级黄色视频| 欧美3d第一页| 亚洲国产精品专区欧美| 精品一品国产午夜福利视频| 在线观看美女被高潮喷水网站| 高清视频免费观看一区二区| 美女内射精品一级片tv| 天堂中文最新版在线下载| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 自线自在国产av| 国产乱来视频区| 国产精品偷伦视频观看了| 亚洲国产欧美日韩在线播放 | 一级爰片在线观看| 中国三级夫妇交换| 国产av一区二区精品久久| 亚洲精品久久午夜乱码| 欧美区成人在线视频| 性色avwww在线观看| 纵有疾风起免费观看全集完整版| 少妇熟女欧美另类| 最黄视频免费看| 少妇精品久久久久久久| 亚洲av电影在线观看一区二区三区| 久久久久久久久久成人| 大片免费播放器 马上看| 成人毛片60女人毛片免费| 韩国av在线不卡| 精华霜和精华液先用哪个| 国产亚洲av片在线观看秒播厂| xxx大片免费视频| 午夜免费男女啪啪视频观看| 亚洲欧美日韩卡通动漫| 日本免费在线观看一区| 免费在线观看成人毛片| 99热国产这里只有精品6| 免费黄频网站在线观看国产| 激情五月婷婷亚洲| 九草在线视频观看| 日韩亚洲欧美综合| 国产精品偷伦视频观看了| 乱人伦中国视频| 国产精品三级大全| 伊人久久精品亚洲午夜| 国内精品宾馆在线| 国产精品无大码| videos熟女内射| 国产精品欧美亚洲77777| 欧美日本中文国产一区发布| 美女cb高潮喷水在线观看| 日韩一本色道免费dvd| 中文字幕亚洲精品专区| 大片电影免费在线观看免费| 欧美+日韩+精品| 99九九在线精品视频 | 日本黄色片子视频| 99热网站在线观看| 少妇被粗大的猛进出69影院 | 亚洲天堂av无毛| 中文字幕人妻熟人妻熟丝袜美| 大码成人一级视频| 一个人看视频在线观看www免费| 黄色毛片三级朝国网站 | 人妻系列 视频| 菩萨蛮人人尽说江南好唐韦庄| 久久国产精品男人的天堂亚洲 | av免费在线看不卡| a级片在线免费高清观看视频| 精品少妇久久久久久888优播| 老司机影院毛片| 欧美精品亚洲一区二区| 亚洲欧美成人综合另类久久久| 国产亚洲av片在线观看秒播厂| 国产亚洲午夜精品一区二区久久| 久久av网站| 亚洲欧美清纯卡通| 国产精品国产三级专区第一集| 欧美成人精品欧美一级黄| 国产成人91sexporn| 一级爰片在线观看| 午夜激情久久久久久久| 青春草视频在线免费观看| 亚洲精品色激情综合| 又黄又爽又刺激的免费视频.| 一区在线观看完整版| 欧美丝袜亚洲另类| 少妇人妻一区二区三区视频| 热re99久久精品国产66热6| 成人黄色视频免费在线看| 中文字幕制服av| 亚洲美女视频黄频| 国产在线一区二区三区精| 女人精品久久久久毛片| 免费黄频网站在线观看国产| 国产亚洲91精品色在线| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级| 在线精品无人区一区二区三| 免费观看无遮挡的男女| 国产精品.久久久| 插阴视频在线观看视频| 成年美女黄网站色视频大全免费 | 亚洲人与动物交配视频| 亚洲欧美精品专区久久| 在线观看一区二区三区激情| 视频中文字幕在线观看| 午夜福利在线观看免费完整高清在| a级毛色黄片| 久久久久久久久大av| 中国美白少妇内射xxxbb| 人妻人人澡人人爽人人| 只有这里有精品99| 在线观看一区二区三区激情| 久久午夜综合久久蜜桃| 伊人久久精品亚洲午夜| 好男人视频免费观看在线| 亚洲av日韩在线播放| 日韩 亚洲 欧美在线| 日韩中字成人| videos熟女内射| 午夜久久久在线观看| 久久精品国产亚洲av涩爱| av国产精品久久久久影院| 国产精品一区二区在线观看99| 国产淫片久久久久久久久| 最新中文字幕久久久久| 亚洲人成网站在线观看播放| videossex国产| 国产午夜精品久久久久久一区二区三区| 乱人伦中国视频| 9色porny在线观看| 久久人妻熟女aⅴ| 中文字幕亚洲精品专区| 伦理电影免费视频| 欧美日韩一区二区视频在线观看视频在线| 三级国产精品片| 一本—道久久a久久精品蜜桃钙片| 日韩在线高清观看一区二区三区| 成人无遮挡网站| 91aial.com中文字幕在线观看| 中文在线观看免费www的网站| 男女国产视频网站| 最近中文字幕2019免费版| 亚洲精品中文字幕在线视频 | 久久精品国产亚洲av天美| 亚洲精品国产成人久久av| 人体艺术视频欧美日本| 亚洲国产欧美日韩在线播放 | 亚洲,欧美,日韩| videossex国产| 日本av免费视频播放| 久久这里有精品视频免费| 六月丁香七月| av在线观看视频网站免费| 一级爰片在线观看| 国产成人a∨麻豆精品| 自拍偷自拍亚洲精品老妇| 青春草视频在线免费观看| 精品国产国语对白av| 国产黄频视频在线观看| 免费人成在线观看视频色| av国产久精品久网站免费入址| 亚洲精品乱码久久久久久按摩| 亚洲精品乱码久久久v下载方式| √禁漫天堂资源中文www| 99精国产麻豆久久婷婷| 亚洲精品中文字幕在线视频 |