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

    北極水汽輸送的時(shí)空分布特征研究

    2022-10-15 04:15:30溫世強(qiáng)常亮
    極地研究 2022年3期
    關(guān)鍵詞:北大西洋水汽北極

    溫世強(qiáng) 常亮,2,3

    (1上海海洋大學(xué), 海洋科學(xué)學(xué)院, 上海 201306;2自然資源部第二海洋研究所, 衛(wèi)星海洋動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 浙江 杭州 310012;3中國(guó)極地研究中心, 自然資源部極地科學(xué)重點(diǎn)實(shí)驗(yàn)室, 上海 200136)

    0 引言

    過(guò)去的半個(gè)世紀(jì), 北極是全球地表氣溫增暖最劇烈的地區(qū), 增溫幅度高達(dá)1.2 ℃ · (10 a)-1, 平均增溫幅度是全球的2倍以上, 在冬季尤為顯著達(dá)到4倍以上, 這種現(xiàn)象被稱(chēng)為“北極放大”[1-4]。2015年12月至2016年2月期間, 北極地表氣溫創(chuàng)造了有觀測(cè)記錄以來(lái)的最暖記錄, 比2011—2012年冬季高出0.7℃, 隨后在2016—2017年和2017—2018年冬季又連續(xù)打破最暖記錄[5]。北極的快速升溫不但引起了本區(qū)域的海冰、冰蓋和凍土等的快速消融, 同時(shí)還通過(guò)改變大氣環(huán)流特征而影響北半球中低緯度的天氣和氣候[6-8]。研究發(fā)現(xiàn)北半球大氣的向極能量傳輸能夠?qū)⒋罅颗瘽駳饬鬏斔偷奖睒O地區(qū), 使得北極大氣中水汽集聚、云量增加, 促使大氣的下行長(zhǎng)波輻射增強(qiáng), 進(jìn)而影響地表輻射收支, 導(dǎo)致地表升溫[9-12]。

    水汽輸送是北半球大氣向極能量傳輸?shù)闹匾绞? 主要由濕度梯度驅(qū)動(dòng), 并受大尺度環(huán)流的影響, 包括平均經(jīng)向環(huán)流、準(zhǔn)靜止渦旋和瞬態(tài)渦旋三種基本輸送形式[13]。北極水汽輸送深刻影響了局地淡水循環(huán)、能量循環(huán)以及海洋冰凍過(guò)程。從熱力學(xué)角度而言, 水汽作為最顯著的溫室氣體,通過(guò)直接影響輻射和間接影響云的形成[14], 加強(qiáng)了云-反照率-地表氣溫之間的反饋過(guò)程[15], 增強(qiáng)了北極放大效應(yīng), 導(dǎo)致北極升溫和海冰減少[16]。從動(dòng)力學(xué)角度而言, 水汽輸送會(huì)破壞冬季北極地區(qū)穩(wěn)定的“逆溫”結(jié)構(gòu), 加快大氣活動(dòng), 影響海冰漂流的過(guò)程[17]?;诖髿饧t外探測(cè)儀(Atmospheric Infrared Sounder, AIRS)、先進(jìn)微波探測(cè)器(Advanced Microwave Sounding Unit-A,AMSU-A)遙感數(shù)據(jù)和ERA-Interim再分析資料,Boisvert等人[18]發(fā)現(xiàn)北極地區(qū)水汽輸送具有巨大的年際變化, 2003年、2004年和2010年均是水汽輸送的較小值年, 而2005年、2007年和2008年則是較大值年。Naakka等人[19]使用ERA-Interim再分析資料, 發(fā)現(xiàn)北極地區(qū)絕大多數(shù)向北水汽輸送被向南輸送抵消, 且在不同高度水汽輸送的方向不同。Dufour等人[20]認(rèn)為再分析資料會(huì)過(guò)高估算水汽輸送, 在1979—2013年水汽輸送呈現(xiàn)下降趨勢(shì), 這與他人的研究結(jié)論相悖[9,21-22], 此外還發(fā)現(xiàn)北大西洋、北太平洋和拉布拉多海是中緯度地區(qū)水汽進(jìn)入北極的主要通道[20]。Gimeno等人[21]使用拉格朗日方法對(duì)ERA-Interim再分析資料估算的北極水汽輸送進(jìn)行后向軌跡分析后認(rèn)為, 水汽輸送的源頭位于北大西洋、北太平洋和地中海。Johansson等人[17]使用AIRS和中分辨率成像光譜儀(Moderate-resolution Imaging Spectroradiometer,MODIS)遙感數(shù)據(jù)發(fā)現(xiàn)水汽輸送不僅造成北極地表溫度的升高, 同時(shí)破壞了近地表大氣的穩(wěn)定性。Liu等人[23]根據(jù)ERA-Interim再分析資料和MODIS云產(chǎn)品發(fā)現(xiàn), 水汽入侵會(huì)導(dǎo)致北極地區(qū)出現(xiàn)更多的云, 以及更高的云水路徑。王岑等人[24]使用NCEP/NCAR再分析資料發(fā)現(xiàn)2015年12月29日北極地區(qū)增加的向北水汽輸送通過(guò)云輻射強(qiáng)迫和非絕熱作用共同導(dǎo)致了當(dāng)?shù)氐谋l(fā)性增溫。Boisvert和Stroeve[25]使用AIRS遙感數(shù)據(jù)分析發(fā)現(xiàn), 5月之前的水汽輸送使得北冰洋海冰提前融化, 而5月之后的水汽輸送促使云層變厚, 進(jìn)而導(dǎo)致低層大氣溫度和輻射正異常, 推遲了北極海冰的凍結(jié)時(shí)間。

    已有的相關(guān)研究針對(duì)北極水汽輸送的氣候特征、來(lái)源和異常變化, 以及其對(duì)局地氣候的影響等開(kāi)展了一系列的分析, 但是在北極水汽輸送通道的分布和增長(zhǎng)趨勢(shì)等方面還存在一些分歧。為了進(jìn)一步認(rèn)識(shí)全球變暖背景下北極地區(qū)水汽輸送的時(shí)空變化特征, 本文基于ERA5、JRA-55和MERRA-2三種再分析資料探討了北極水汽輸送的時(shí)空分布、變化趨勢(shì)等方面的特征, 并分析比較了三者間的差異性。本文在1.1節(jié)和1.2節(jié)介紹所使用的數(shù)據(jù)和處理方法。為了評(píng)價(jià)不同再分析資料進(jìn)行水汽輸送估計(jì)的精度, 在1.3節(jié)將再分析資料與第二代高質(zhì)量全球無(wú)線(xiàn)電探空數(shù)據(jù)集(Integrated Global Radiosonde Archive 2,IGRA2)估計(jì)的水汽通量進(jìn)行了對(duì)比驗(yàn)證。在第2節(jié)基于1980—2018年的再分析資料對(duì)沿70°N的水汽通量的水平分布、垂直分布和月平均分布進(jìn)行了對(duì)比研究, 并系統(tǒng)分析了北極地區(qū)整層水汽通量在不同季節(jié)的變化特征和線(xiàn)性趨勢(shì);同時(shí)還探討了沿70°N水汽入侵的時(shí)間變化和水平分布特征。最后在第3節(jié)歸納了本文的主要結(jié)論。

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

    1.1 再分析資料

    再分析資料的時(shí)空分辨率高, 可以用于在整個(gè)北極地區(qū)開(kāi)展長(zhǎng)時(shí)間尺度的水汽輸送變化研究。本文主要引入了ERA5、JRA-55和MERRA-2三種再分析資料進(jìn)行北極水汽輸送的時(shí)空分布特征研究。歐洲中期數(shù)值天氣預(yù)報(bào)中心(ECMWF)的第五代再分析資料產(chǎn)品ERA5重新處理了ECMWF衛(wèi)星氣候數(shù)據(jù)集, 并運(yùn)用了RTTOV11輻射傳輸模式, 使得數(shù)據(jù)具有更高的時(shí)間分辨率。另外, ERA5同化了更多的觀測(cè)數(shù)據(jù)和衛(wèi)星數(shù)據(jù),能更為準(zhǔn)確地估計(jì)大氣狀況[26]。JRA-55是日本氣象廳(Japan Meteorological Agency, JMA)在2013年發(fā)布的再分析資料, 與上一代再分析資料JRA-25相比, JRA-55采用了更先進(jìn)的四維同化變分方案, 并對(duì)衛(wèi)星觀測(cè)數(shù)據(jù)進(jìn)行了修正[27]。MERRA-2是美國(guó)宇航局全球建模與同化辦公室(Global Modeling And Assimilation Office, GMAO)于2017年發(fā)布的一項(xiàng)再分析資料, 與上一代版本MERRA相比, MERRA-2同化了微波探測(cè)儀和高光譜紅外輻射儀等新的衛(wèi)星觀測(cè)數(shù)據(jù), 并使用了新的地球觀測(cè)系統(tǒng)模式(GEOS-5)[28]。本文使用的ERA5、JRA-55和MERRA-2的基本參數(shù)如表1所示。

    表1 ERA5、JRA-55和MERRA-2基本參數(shù)Table 1.Basic parameters of ERA5, JRA-55 and MERRA-2

    本文選取了1980年1月1日至2018年12月31日ERA5、JRA-55和MERRA-2再分析資料逐6小時(shí)的經(jīng)向風(fēng)和比濕數(shù)據(jù), 估算北極地區(qū)的水汽輸送。為保證使用數(shù)據(jù)的空間分辨率一致性,本文選取的ERA5空間分辨率為0.5°, 并使用線(xiàn)性插值法將JRA-55和MERRA-2的空間分辨率重采樣至0.5°的網(wǎng)格上。水汽輸送在水平方向可分為經(jīng)向和緯向方向, 經(jīng)向方向水汽輸送直接影響北極的局地氣候[29], 所以本文主要研究經(jīng)向方向的水汽輸送。水汽輸送通常以水汽通量的大小進(jìn)行表示, 并規(guī)定水汽通量的正值為向北水汽輸送,負(fù)值為向南。連續(xù)氣壓層pi與pi+1間的水汽通量fqj可以用下式進(jìn)行計(jì)算[19]:

    式中,g是重力加速度, 其值取0.98(m·s-2),q是比濕(kg·kg-1),v是經(jīng)向風(fēng)(m·s-1),p是大氣壓強(qiáng)(hPa),下標(biāo)i表示連續(xù)氣壓層中的某個(gè)氣壓層, 下標(biāo)j指示連續(xù)氣壓層pi和pi+1之間的高度范圍。

    對(duì)地表至500 hPa高度大氣層內(nèi)的整層水汽通量Fq進(jìn)行計(jì)算[19]:

    式中,n是氣壓層數(shù)。需要說(shuō)明的是, 由于北極地區(qū)近98%的大氣水汽均集中在500 hPa以下[30],所以本文沒(méi)有估算高度500 hPa以上的水汽通量。

    1.2 IGRA2數(shù)據(jù)

    無(wú)線(xiàn)電探空觀測(cè)是開(kāi)展北極地區(qū)水汽輸送特征研究的重要手段。為了評(píng)估再分析資料估算水汽通量數(shù)據(jù)的精度, 本文引入了1980—2018年美國(guó)國(guó)家氣候數(shù)據(jù)中心的IGRA2數(shù)據(jù)集進(jìn)行比較驗(yàn)證。IGRA2包括了全球2700多個(gè)分布式站點(diǎn)(無(wú)線(xiàn)電探空儀和探空氣球)的露點(diǎn)溫度、壓強(qiáng)和水汽壓等大氣參數(shù)的逐日數(shù)據(jù)集, 每個(gè)站點(diǎn)數(shù)據(jù)中的壓力、溫度和位勢(shì)高度均經(jīng)過(guò)了嚴(yán)格的質(zhì)量控制, 但是水汽壓和風(fēng)速存在較大偏差。為了保證水汽通量估算結(jié)果的質(zhì)量, 本文以5.4倍標(biāo)準(zhǔn)差(Standard Deviation, STD)作為閾值剔除了水汽壓和風(fēng)速的異常值, 并僅選取站點(diǎn)的有效數(shù)據(jù)次數(shù)高于整個(gè)觀測(cè)周期總數(shù)85%的站點(diǎn)[20]。按照此標(biāo)準(zhǔn), 本文從北極地區(qū)75個(gè)探空測(cè)站中篩選獲得了64個(gè)站點(diǎn)進(jìn)行后續(xù)分析, 篩選出的探空測(cè)站空間分布結(jié)果如圖1所示。

    圖1 本文選用的北極地區(qū)IGRA2站點(diǎn)的空間分布圖。數(shù)字表示探空測(cè)站的世界氣象組織(World Meteorological Organization, WMO)編號(hào)Fig.1.Spatial distribution of IGRA2 sites in the Arctic region selected in this article.The number indicates the World Meteorological Organization (WMO) number of the radiosonde site

    考慮到IGRA2沒(méi)有直接提供比濕數(shù)據(jù), 可以根據(jù)其提供的大氣壓強(qiáng)(p)和水汽壓參數(shù)(e)推算比濕(q), 計(jì)算公式如下[31]:

    式中,w是水汽混合比(kg·kg-1),e是水汽壓(hPa),ε是常量, 取值為0.6219907。

    1.3 再分析資料估算的水汽通量的精度分析

    為了評(píng)估ERA5、JRA-55和MERRA-2估算北極地區(qū)整層水汽通量的精度, 本文將1980—2018年再分析資料與時(shí)空同步的利用IGRA2估算的整層水汽通量結(jié)果進(jìn)行了比較, 兩者之間的相關(guān)系數(shù)(r)和均方根誤差(Root Mean Square Error,RMSE)的空間分布結(jié)果如圖2所示。ERA5和JRA-55與IGRA2水汽通量的相關(guān)性在空間分布上較為一致, 在多數(shù)站點(diǎn)呈現(xiàn)出較強(qiáng)的正相關(guān)性(有49個(gè)站點(diǎn)的r>0.5, 所有站點(diǎn)的平均相關(guān)系數(shù)分別為0.65和0.7), 但ERA5和JRA-55在阿拉斯加、格陵蘭島北側(cè)和加拿大北極群島北側(cè)站點(diǎn)的r較低(圖2a和2b)。MERRA-2與IGRA2估算的水汽通量在大多數(shù)站點(diǎn)上均具有較強(qiáng)的正相關(guān)性(有56個(gè)站點(diǎn)的r>0.5, 所有站點(diǎn)的平均相關(guān)系數(shù)達(dá)到0.78), 僅在西伯利亞北部和北歐個(gè)別測(cè)站的相關(guān)性較低(圖2c), 這說(shuō)明MERRA-2比ERA5和JRA-55更能捕捉到北極大范圍內(nèi)的水汽輸送變化特征。ERA5與JRA-55估算的水汽通量的RMSE空間分布特征相似, 但在阿拉斯加地區(qū)的數(shù)值分別偏高21%和27%, 在挪威海附近分別偏高34%和30%(圖2d和2e), 可能原因是這些地區(qū)的相對(duì)高海拔地形使西風(fēng)偏轉(zhuǎn)為南風(fēng), 造成了低層大氣的水汽通量數(shù)值較大[19], 而再分析資料的同化方案受地形影響未能準(zhǔn)確捕獲這種變化而導(dǎo)致誤差偏大。MERRA-2水汽通量的RMSE在北極相較ERA5和JRA-55偏低5.5%~8%(圖2d和2f), 說(shuō)明MERRA-2表征北極水汽通量的精度更高。此外, 在西伯利亞地區(qū)ERA5和JRA-55的r值(平均值分別為0.8和0.86)高于MERRA-2(平均值為 0.63), 而RMSE值(平均值分別為129 kg·m-1·s-1和122 kg·m-1·s-1)較MERRA-2(平均為137 kg·m-1·s-1)偏低, 這說(shuō)明ERA5和JRA-55的適用性更高??傮w而言, 與ERA5和JRA-55相比, MERRA-2估計(jì)的北極水汽通量與IGRA2具有更高的一致性和精度, 更適用于北極水汽輸送的研究。

    圖2 1980—2018年不同再分析資料估算的北極整層水汽通量精度。a)—c)分別為ERA5、JRA-55、MERRA-2與IGRA2的r值情況; d)—f)分別為ERA5、JRA-55、MERRA-2與IGRA2的RMSE值情況Fig.2.Accuracy of integrated water vapor flux with different reanalysis data over the Arctic from 1980 to 2018.a)—c) indicated the r value between ERA5, JRA-55, MERRA-2 and IGRA2, respectively; d)—f) indicated the RMSE value between ERA5, JRA-55, MERRA-2 and IGRA2, respectively

    2 北極地區(qū)水汽輸送的時(shí)空變化特征

    2.1 北極地區(qū)水汽輸送的時(shí)空變化特征

    圖3給出了基于ERA5、JRA-55和MERAA-2估算的1980—2018年北極極冠邊界處(即沿70°N)多年平均的整層水汽通量的水平分布結(jié)果。整體上, 三種再分析資料的水汽通量分布一致性較高,在東北極均為正值, 在西北極正值和負(fù)值相間分布且數(shù)值較東北極偏高, 這與吳靜[32]的結(jié)論相符。水汽通量的這一水平分布特征說(shuō)明水汽輸送在東北極以向北為主, 而在西北極同時(shí)存在向北和向南的水汽輸送且較東北極更活躍, 這可能是因?yàn)闁|北極多數(shù)區(qū)域?yàn)殛懙? 而西北極區(qū)域內(nèi)海洋和陸地的占比相近, 不同下墊面的差異導(dǎo)致了水汽輸送方向和強(qiáng)度不同?;谘?0°N水汽通量的極值分布, 可以揭示水汽輸送通道所在的位置。水汽通量極大值反映的向北水汽輸送通道,主要位于170°W楚科奇海、50°W巴芬灣、0°北大西洋和90°E中西伯利亞地區(qū), 這與Woods等人[33]的結(jié)論一致, 但比Dufour等人[20]的結(jié)論增加了90°E中西伯利亞地區(qū), 原因可能是Dufour等人[20]僅分析了1—4月和7—10月兩個(gè)時(shí)間段的水汽通量, 而本文與Woods等人[33]則分析了所有月份。水汽通量極小值反映的向南水汽輸送通道, 則主要位于100°W加拿大北極群島和25°W丹麥海峽地區(qū), 這些通道附近的水汽輸送除了與下墊面性質(zhì)有關(guān), 也較大程度地受局地氣象要素的影響。在楚科奇海、巴芬灣和北大西洋等區(qū)域, 水汽通量的大小均超過(guò)18 kg·m-1·s-1, 且范圍超過(guò)50°,這主要是受向北移動(dòng)的氣旋影響[34]; 在加拿大北極群島上空, 向南水汽輸送較為活躍, 這主要是當(dāng)?shù)仄珡?qiáng)的北風(fēng)和較弱的氣旋活動(dòng)共同作用的結(jié)果[35-36]; 丹麥海峽的向南水汽輸送則是冰島低壓產(chǎn)生的近地面北風(fēng)造成的[20]; 在中西伯利亞, 水汽輸送數(shù)值和范圍相對(duì)較小, 主要原因可能是亞歐大陸的風(fēng)速和濕度均值相較海洋偏小。

    圖3 1980—2018年ERA5、JRA-55和MERRA-2沿70°N多年平均的整層水汽通量的水平分布圖Fig.3.Horizontal distribution of multiyear average integrated water vapor flux at 70°N of ERA5, JRA-55, and MERRA-2,respectively, from 1980 to 2018

    為詳細(xì)研究沿70°N水汽輸送的垂直結(jié)構(gòu),圖4給出了該緯度圈多年平均的水汽通量的垂直分布結(jié)果。ERA5和JRA-55所估算的1980—2018年沿70°N多年平均的水汽通量在垂直方向的分布特征較為相似(圖4a和4b), 與Serreze等人[30]的研究結(jié)論相吻合, 這兩種資料的水汽通量正值主要分布于850~950 hPa高度, 并在900 hPa附近出現(xiàn)水汽通量極大值, 這說(shuō)明向北水汽輸送主要集中于850~950 hPa高度, 并在900 hPa較為強(qiáng)烈,這可能是因?yàn)閷?duì)流層內(nèi)風(fēng)速通常隨著高度升高而增大, 濕度隨著高度升高而減小, 兩者在900 hPa高度達(dá)到平衡狀態(tài)。MERRA-2的有效水汽通量結(jié)果在垂直方向上也具有與ERA5和JRA-55相似的空間分布特征。此外, 900 hPa高度附近的水汽通量極大值主要位于170°W楚科奇海、50°W巴芬灣、0°北大西洋和90°E中西伯利亞地區(qū)(圖4a和4b), 與圖3中整層水汽通量極大值的位置相近,這與 Naakka等人[19]的分析結(jié)論一致, 說(shuō)明900 hPa高度水汽通量能夠大致代表整層水汽通量的變化。

    圖4 1980—2018年ERA5、JRA-55和MERRA-2沿70°N多年平均的水汽通量的垂直分布圖。a) ERA5; b) JRA-55;c) MERRA-2。c)圖中的空白區(qū)域是因?yàn)镸ERRA-2的v和q變量在此處存在無(wú)效數(shù)據(jù)Fig.4.Vertical distribution of multiyear average water vapor flux at 70°N of ERA5, JRA-55, and MERRA-2, respectively,from 1980 to 2018.a) ERA5; b) JRA-55; c) MERRA-2.The blank area in the c) is caused by invalid v and q in MERRA-2

    圖5繪制了ERA5、JRA-55和MERRA-2在1980—2018年沿70°N處多年月平均的整層水汽通量的變化結(jié)果。三種再分析資料得到的沿70°N多年月平均水汽通量均為正值, 說(shuō)明向北水汽輸送通量平均高于向南, 北極地區(qū)整體體現(xiàn)為水汽匯入地。就年際變化特性而言, 三種資料基本一致, 均呈單峰型分布, 1—7月水汽通量逐月增長(zhǎng),在7月達(dá)到最高(約10 kg·m-1·s-1), 7—12月逐月降低, 在12月達(dá)到最低(約為4 kg·m-1·s-1), 夏季(6—8月)極大值約為冬季(12月至次年2月)極小值的2倍, 這可能是夏季地表氣溫升高增加了大氣保持水分的能力, 空氣濕度較冬季更高所導(dǎo)致的[30]。

    圖5 1980—2018年ERA5、JRA-55和MERRA-2沿70°N多年月平均的整層水汽通量的變化圖Fig.5.Multiyear monthly average of integrated water vapor flux at 70°N of ERA5, JRA-55 and MERRA-2, respectively, from 1980 to 2018

    圖6繪制了1980—2018年ERA5、JRA-55和MERRA-2三種再分析資料在北極季節(jié)平均的整層水汽通量的空間分布結(jié)果。三種再分析資料均顯示向北水汽輸送的季節(jié)變化比向南輸送更加明顯, 這可能與中緯度與副極地的溫帶氣旋和濕度在一年中的劇烈變化有關(guān)[20]。北極地區(qū)的向北水汽輸送強(qiáng)度在夏(6—8月)、秋(9—11月)、春(3—5月)、冬(12月至次年2月)季依次遞減。夏季向北水汽輸送最為活躍, 北大西洋、楚科奇海和巴芬灣的水汽通量值均達(dá)到50 kg·m-1·s-1以上(圖6c、6g和6k), 這可能是北冰洋及其鄰近海域海冰融化 , 開(kāi)闊水域增多, 海-氣相互作用增強(qiáng)導(dǎo)致的[37]。秋季北極大部分的水汽通量值低于20 kg·m-1·s-1,楚科奇海和巴芬灣區(qū)域約為40 kg·m-1·s-1, 僅北大西洋區(qū)域達(dá)到60 kg·m-1·s-1以上(圖6d、6h和6l), 這可能是當(dāng)?shù)貪穸仁芎1诨绊懚眲≡龆嘣斐傻腫38]。冬季北大西洋海域的水汽通量數(shù)值相比春季偏高20 kg·m-1·s-1, 這可能是因?yàn)樵摰厥嵌練庑闹饕顒?dòng)區(qū)域[39]。除巴芬灣以外,ERA5、JRA-55和MERRA-2的水汽通量在北極的季節(jié)分布特征基本一致。在巴芬灣, MERRA-2較ERA5和JRA-55的水汽通量數(shù)值略小(圖6a—6l), 且MERRA-2在該區(qū)域的精度高于ERA5和JRA-55(圖2), 這說(shuō)明ERA5和JRA-55可能高估了巴芬灣附近的水汽通量。

    圖6 1980—2018年ERA5、JRA-55和MERRA-2在北極多年季節(jié)平均的整層水汽通量的空間分布圖。a)—d) 分別為ERA5在冬、春、夏、秋季的結(jié)果; e)—h) 分別為JRA-55在冬、春、夏、秋季的結(jié)果; i)—l) 分別為MERRA-2在冬、春、夏、秋季的結(jié)果Fig.6.Spatial distribution of multiyear seasonal average integrated water vapor flux over the Arctic from 1980 to 2018.a)—d)indicate winter, spring, summer and autumn results of ERA5 respectively; e)—h) indicate winter, spring, summer and autumn results of JRA55 respectively; i)—l) indicate winter, spring, summer and autumn results of MERRA-2 respectively

    圖7繪制了1980—2018年ERA5、JRA-55和MERRA-2三種資料在北極按季節(jié)平均的整層水汽通量線(xiàn)性趨勢(shì)的空間分布結(jié)果。三種再分析資料均顯示, 在不同季節(jié)水汽通量的變化趨勢(shì)特征差異較大。在夏季, 向北水汽輸送通道附近(0°北大西洋、90°E中西伯利亞、50°W巴芬灣地區(qū))水汽通量的線(xiàn)性增長(zhǎng)速度達(dá)到0.4 kg·m-1·s-1·a-1以上(圖7c、7g和7k)。前人的研究[40-42]表明, 北大西洋和巴芬灣附近水汽通量的增速較快分別是由氣溫升高導(dǎo)致的濕度上升和氣旋活動(dòng)增多導(dǎo)致的, 而東西伯利亞海和中西伯利亞地區(qū)的水汽通量增長(zhǎng)分別與北極偶極子異常和喀拉海位勢(shì)高度負(fù)異常有關(guān)。秋季北大西洋附近海域的水汽通量增速達(dá)到0.8 kg·m-1·s-1·a-1以上(圖7d、7h和7l),這可能與海冰融化導(dǎo)致當(dāng)?shù)貪穸鹊募眲≡黾佑嘘P(guān)[38]。冬季水汽通量增長(zhǎng)最為強(qiáng)烈的區(qū)域主要在北大西洋區(qū)域(圖7a、7e和7i), 春季則主要在西伯利亞區(qū)域(圖7b、7f和7j), 但與夏季相比, 這兩個(gè)區(qū)域的水汽通量增長(zhǎng)趨緩, 均低于0.4 kg·m-1·s-1·a-1。

    圖7 1980—2018年ERA5、JRA-55和MERRA-2在北極多年季節(jié)平均的整層水汽通量的線(xiàn)性趨勢(shì)分布圖。a)—d) 分別為ERA5在冬、春、夏、秋季的結(jié)果; e)—h) 分別為JRA-55在冬、春、夏、秋季的結(jié)果; i)—l) 分別為MERRA-2在冬、春、夏、秋季的結(jié)果Fig.7.Linear trend of multiyear seasonal average integrated water vapor flux over the Arctic from 1980 to 2018.a)—d) indicate winter, spring, summer and autumn results of ERA5 respectively; e)—h) indicate winter, spring, summer and autumn results of JRA55 respectively; i)—l) indicate winter, spring, summer and autumn results of MERRA-2 respectively

    2.2 沿70°N水汽入侵的時(shí)間序列和水平分布

    水汽入侵是指大氣中強(qiáng)烈的水汽輸送事件[33],能夠造成向下長(zhǎng)波輻射和云水路徑正異常[23], 從而引發(fā)北極海冰減退。依據(jù)Woods等人[33]對(duì)水汽入侵的定義, 本文對(duì)沿70°N緯度帶東西方向上9°的范圍內(nèi)水汽通量值超過(guò)總數(shù)90%, 并至少持續(xù)1.5 d的水汽入侵現(xiàn)象進(jìn)行了統(tǒng)計(jì)分析。圖8繪制了分別基于ERA5、JRA-55和MERRA-2數(shù)據(jù), 在1980—2018年沿70°N水汽入侵的年平均次數(shù)。1980—2018年ERA5、JRA-55和MERRA-2的水汽入侵次數(shù)平均為每年27~29次, 2010—2011年為水汽入侵的高值年, ERA5為43次,JRA-55為40次, MERRA-2為40次。從圖8還可以看出, 2008年之后的年平均水汽入侵次數(shù)出現(xiàn)較大波動(dòng), 這與采用滑動(dòng)t檢驗(yàn)方法對(duì)三種再分析資料的水汽入侵年平均次數(shù)的統(tǒng)計(jì)是一致的(檢驗(yàn)結(jié)果未予顯示)。此外, 2008年前后的水汽輸送突變主要是由北極地區(qū)低對(duì)流層濕度的快速異常變化引起的。

    圖8 1980—2018年分別基于ERA5、JRA-55和MERRA-2數(shù)據(jù)的沿70°N水汽入侵次數(shù)的年平均次數(shù)Fig.8.Annual average of water vapor intrusion at 70°N from 1980 to 2018 based on ERA5, JRA-55 and MERRA-2

    圖9繪制了1980—2018年分別基于ERA5、JRA-55和MERRA-2再分析資料中沿70°N處水汽入侵的水平分布結(jié)果。三種再分析資料得到的水汽入侵水平分布相似, 主要集中于0°北大西洋、180°楚科奇海、50°W巴芬灣和20°E—90°E西伯利亞區(qū)域, 與上文中向北水汽輸送的通道位置(圖3)相符, 同時(shí)也與Woods和Caballero[43]的結(jié)論吻合,說(shuō)明了水汽入侵主要發(fā)生在向北水汽輸送通道附近。水汽入侵年平均次數(shù)的極大值出現(xiàn)在北大西洋(圖9), 這可能是當(dāng)?shù)貧庑男眽合到y(tǒng)引起的瞬態(tài)渦旋輸送偏強(qiáng)導(dǎo)致的[29]。此外, 在巴芬灣出現(xiàn)水汽入侵的次數(shù)較低, 且多年平均的整層水汽通量數(shù)值最高(圖3), 說(shuō)明當(dāng)?shù)氐南虮彼斔洼^為平穩(wěn), 這可能是由于巴芬灣附近的濕度受到西風(fēng)急流和大尺度的大氣環(huán)流調(diào)節(jié)而呈穩(wěn)定狀態(tài)[42]。

    圖9 1980—2018年基于ERA5、JRA-55和MERRA-2數(shù)據(jù)的沿70°N水汽入侵的水平分布圖Fig.9.Horizontal distribution of water vapor intrusion at 70°N from 1980 to 2018 based on ERA5, JRA-55 and MERRA-2

    3 結(jié)論

    本文首先以IGRA2無(wú)線(xiàn)電探空數(shù)據(jù)為參照,評(píng)估了1980—2018年ERA5、JRA-55和MERRA-2三種再分析資料估計(jì)的水汽通量結(jié)果在北極地區(qū)的適用性, 并基于三種再分析資料對(duì)沿70°N以及整個(gè)北極地區(qū)的水汽輸送時(shí)空變化特征進(jìn)行了分析, 得到如下主要結(jié)論。

    1.ERA5和JRA-55與IGRA2估計(jì)的北極水汽通量的平均相關(guān)性分別為 0.65和 0.7,MERRA-2與IGRA2的平均相關(guān)性達(dá)到0.78。以IGRA2的估計(jì)結(jié)果為參考, ERA5和JRA-55估計(jì)的水汽通量的RMSE比MERRA-2平均偏高5.5%~8%。整體而言, MERRA-2估計(jì)的北極水汽通量與IGRA2具有更高的一致性, 比ERA5和JRA-55更適用于開(kāi)展北極地區(qū)的水汽輸送研究。

    2.沿70°N水汽通量的空間分布顯示, 在水平方向上, 水汽輸送在東北極基本以向北輸送為主, 在西北極則同時(shí)存在向南和向北輸送且相對(duì)東北極更活躍; 在垂直方向上, 向北水汽輸送在900 hPa高度較為強(qiáng)烈。此外, 沿70°N水汽通量的月平均分布結(jié)果還表明, 向北水汽平均輸送強(qiáng)度高于向南水汽輸送, 北極地區(qū)整體上體現(xiàn)為水汽匯入地。

    3.北極地區(qū)的向北水汽輸送呈現(xiàn)不同的季節(jié)性和空間分布特征, 具體表現(xiàn)為夏秋季的向北水汽輸送較冬春季更加活躍, 水汽通量在北大西洋、楚科奇海和巴芬灣附近均超過(guò)40 kg·m-1·s-1,并分別在向北水汽輸送通道附近和北大西洋海域以超過(guò)0.4 kg·m-1·s-1·a-1和0.8 kg·m-1·s-1·a-1的速度增長(zhǎng); 冬季水汽通量比春季在北大西洋區(qū)域偏高10 kg·m-1·s-1, 同時(shí), 冬春季節(jié)水汽通量在北大西洋和西伯利亞的增長(zhǎng)速度較為平緩, 均低于0.4 kg·m-1·s-1·a-1。

    4.受北極地區(qū)近地表比濕快速增加的影響,2008年以后沿70°N水汽入侵的年均次數(shù)存在明顯的波動(dòng)。同時(shí), 北極地區(qū)水汽入侵的發(fā)生位置與中緯度水汽進(jìn)入北極的輸送通道位置一致, 主要分布在0°的北大西洋、180°的楚科奇海、50°W的巴芬灣和20°E—90°E的西伯利亞區(qū)域。

    猜你喜歡
    北大西洋水汽北極
    青藏高原上空平流層水汽的時(shí)空演變特征
    近60年華北春季干旱特征及其與北大西洋海表溫度的關(guān)系
    北極有個(gè)“放屁湖”
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    北大西洋海浪特征分析
    北極兔乖乖,唱起歌來(lái)
    北極
    心聲歌刊(2017年4期)2017-09-20 11:43:48
    去北極嘍!
    幼兒園(2016年22期)2017-03-01 13:57:52
    1979~2011年間平流層溫度及平流層水汽的演變趨勢(shì)
    與北大西洋接壤的北極海冰和年際氣候變化
    人人澡人人妻人| 亚洲 国产 在线| 国产精品免费一区二区三区在线 | 午夜免费鲁丝| 身体一侧抽搐| 999久久久精品免费观看国产| 色老头精品视频在线观看| 狠狠狠狠99中文字幕| 黄频高清免费视频| 91成人精品电影| 自线自在国产av| 国产精品二区激情视频| 亚洲精品美女久久av网站| 成人免费观看视频高清| 国产欧美亚洲国产| 欧美激情极品国产一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产男女超爽视频在线观看| 91麻豆精品激情在线观看国产 | 一级片'在线观看视频| 亚洲精品国产一区二区精华液| 最近最新中文字幕大全电影3 | 人妻 亚洲 视频| 久久久久久久国产电影| 久久狼人影院| 亚洲精品美女久久av网站| 欧美激情高清一区二区三区| 大香蕉久久网| 国产成人精品久久二区二区91| 黄片大片在线免费观看| 色综合欧美亚洲国产小说| 国产成人啪精品午夜网站| 99久久99久久久精品蜜桃| 国产视频一区二区在线看| 99国产极品粉嫩在线观看| 80岁老熟妇乱子伦牲交| 99re6热这里在线精品视频| 丁香欧美五月| 精品一品国产午夜福利视频| 精品人妻在线不人妻| 久久天堂一区二区三区四区| 中文字幕人妻丝袜一区二区| 日韩三级视频一区二区三区| 亚洲九九香蕉| av福利片在线| 亚洲午夜理论影院| 在线观看一区二区三区激情| 在线看a的网站| 亚洲成国产人片在线观看| 91在线观看av| 热re99久久精品国产66热6| 嫩草影视91久久| 国产真人三级小视频在线观看| 欧美激情 高清一区二区三区| 免费在线观看完整版高清| 国产av一区二区精品久久| 亚洲视频免费观看视频| 热99久久久久精品小说推荐| 天堂√8在线中文| 久久精品国产亚洲av高清一级| 夜夜夜夜夜久久久久| 国产激情欧美一区二区| 曰老女人黄片| 超碰97精品在线观看| 国产亚洲av高清不卡| 12—13女人毛片做爰片一| 国产高清视频在线播放一区| 欧美精品高潮呻吟av久久| avwww免费| 久久中文字幕人妻熟女| 国产精品偷伦视频观看了| 国产高清国产精品国产三级| 在线观看66精品国产| 久久久精品免费免费高清| 精品福利永久在线观看| 两个人看的免费小视频| 看黄色毛片网站| 精品第一国产精品| 天天添夜夜摸| 90打野战视频偷拍视频| 女性生殖器流出的白浆| 天天添夜夜摸| 中国美女看黄片| 久9热在线精品视频| 亚洲综合色网址| 欧美日韩国产mv在线观看视频| 亚洲精品在线观看二区| 亚洲国产欧美一区二区综合| 99精品在免费线老司机午夜| 三上悠亚av全集在线观看| 国产精品九九99| 精品国产一区二区久久| 成年人午夜在线观看视频| 搡老岳熟女国产| 满18在线观看网站| 久久影院123| 热re99久久精品国产66热6| 成年人免费黄色播放视频| 别揉我奶头~嗯~啊~动态视频| 窝窝影院91人妻| 老熟妇仑乱视频hdxx| 国产三级黄色录像| 亚洲第一欧美日韩一区二区三区| 国产精品98久久久久久宅男小说| 啦啦啦 在线观看视频| 精品熟女少妇八av免费久了| 一本大道久久a久久精品| 国产精品久久久av美女十八| 日韩大码丰满熟妇| 少妇猛男粗大的猛烈进出视频| 男人舔女人的私密视频| 日本a在线网址| 国产激情久久老熟女| 黄色成人免费大全| 国产午夜精品久久久久久| 免费在线观看亚洲国产| 久久精品亚洲熟妇少妇任你| 80岁老熟妇乱子伦牲交| 午夜精品在线福利| 亚洲一区二区三区不卡视频| 国产又爽黄色视频| a级毛片在线看网站| 制服人妻中文乱码| 啦啦啦免费观看视频1| 欧美精品人与动牲交sv欧美| 制服人妻中文乱码| av天堂在线播放| 中文字幕人妻丝袜制服| 成人18禁高潮啪啪吃奶动态图| 国产99白浆流出| 亚洲精品粉嫩美女一区| 精品午夜福利视频在线观看一区| 少妇 在线观看| 人妻一区二区av| 69精品国产乱码久久久| 精品久久久精品久久久| 一进一出好大好爽视频| 久9热在线精品视频| 丝袜美足系列| 91麻豆精品激情在线观看国产 | 精品国产一区二区久久| 高清欧美精品videossex| 欧美在线黄色| 亚洲精华国产精华精| 国产aⅴ精品一区二区三区波| 999精品在线视频| 国产精品自产拍在线观看55亚洲 | 亚洲男人天堂网一区| 69精品国产乱码久久久| 香蕉久久夜色| 欧美老熟妇乱子伦牲交| 黄色片一级片一级黄色片| 亚洲专区字幕在线| 精品人妻熟女毛片av久久网站| 三级毛片av免费| 国产精品免费一区二区三区在线 | 老熟妇乱子伦视频在线观看| 欧美日本中文国产一区发布| 久久精品国产99精品国产亚洲性色 | 黄色女人牲交| 精品国产一区二区久久| 一级a爱视频在线免费观看| 在线免费观看的www视频| 在线观看免费午夜福利视频| 欧美不卡视频在线免费观看 | 亚洲国产欧美一区二区综合| 丝袜美足系列| 18禁国产床啪视频网站| 欧美黄色淫秽网站| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 男女午夜视频在线观看| 国产精品久久视频播放| 韩国av一区二区三区四区| 国产淫语在线视频| 国产极品粉嫩免费观看在线| 黄色成人免费大全| 黑丝袜美女国产一区| 法律面前人人平等表现在哪些方面| 激情视频va一区二区三区| 久久精品人人爽人人爽视色| 久久影院123| 久久精品国产a三级三级三级| 狠狠狠狠99中文字幕| 9色porny在线观看| 美女午夜性视频免费| 男人的好看免费观看在线视频 | ponron亚洲| 免费女性裸体啪啪无遮挡网站| 久久久久精品国产欧美久久久| 老熟妇仑乱视频hdxx| 狠狠狠狠99中文字幕| 亚洲在线自拍视频| 在线观看免费视频网站a站| 中国美女看黄片| 在线观看免费高清a一片| 久久国产精品男人的天堂亚洲| 国产成人影院久久av| 久久人妻av系列| 成人av一区二区三区在线看| 在线看a的网站| 制服人妻中文乱码| 国产男靠女视频免费网站| 亚洲七黄色美女视频| 国产日韩一区二区三区精品不卡| 99国产极品粉嫩在线观看| 亚洲七黄色美女视频| 正在播放国产对白刺激| 十八禁高潮呻吟视频| 一夜夜www| 中文字幕最新亚洲高清| 男女午夜视频在线观看| 国产精品偷伦视频观看了| 一级毛片高清免费大全| 在线永久观看黄色视频| 一进一出抽搐动态| 午夜激情av网站| 人妻一区二区av| 久久久精品国产亚洲av高清涩受| 一二三四在线观看免费中文在| 亚洲熟妇中文字幕五十中出 | 婷婷丁香在线五月| 丝袜在线中文字幕| 久久人妻熟女aⅴ| 成人18禁在线播放| 大型黄色视频在线免费观看| 亚洲精品中文字幕一二三四区| 国内久久婷婷六月综合欲色啪| 热re99久久精品国产66热6| 国产精品自产拍在线观看55亚洲 | 精品一品国产午夜福利视频| 女警被强在线播放| 又紧又爽又黄一区二区| 国产一区二区三区在线臀色熟女 | 成人18禁高潮啪啪吃奶动态图| 久久久久久久国产电影| av超薄肉色丝袜交足视频| 50天的宝宝边吃奶边哭怎么回事| 免费少妇av软件| 欧美日韩成人在线一区二区| 欧美人与性动交α欧美软件| 在线免费观看的www视频| 久久这里只有精品19| 757午夜福利合集在线观看| 精品国产美女av久久久久小说| 制服诱惑二区| 天天躁日日躁夜夜躁夜夜| 国产高清激情床上av| 两性夫妻黄色片| 黑人操中国人逼视频| 一区二区日韩欧美中文字幕| 国产精品秋霞免费鲁丝片| 久久久久国产精品人妻aⅴ院 | 婷婷丁香在线五月| 国产精品影院久久| 亚洲人成电影观看| 高清欧美精品videossex| 18禁国产床啪视频网站| 校园春色视频在线观看| 男女午夜视频在线观看| 久久久国产成人免费| 午夜激情av网站| 久久国产精品大桥未久av| xxx96com| 国产不卡av网站在线观看| 亚洲免费av在线视频| 黄色 视频免费看| 色精品久久人妻99蜜桃| 亚洲午夜精品一区,二区,三区| 女人久久www免费人成看片| www.熟女人妻精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 18禁黄网站禁片午夜丰满| 欧美激情久久久久久爽电影 | 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡动漫免费视频| 国产高清视频在线播放一区| 成年人黄色毛片网站| 美女国产高潮福利片在线看| 免费观看人在逋| 国产99久久九九免费精品| av在线播放免费不卡| 午夜老司机福利片| 一进一出抽搐gif免费好疼 | 国产精品一区二区在线不卡| 久久亚洲精品不卡| 精品欧美一区二区三区在线| 国产在线一区二区三区精| 青草久久国产| 亚洲中文字幕日韩| 老鸭窝网址在线观看| 国产99久久九九免费精品| 在线永久观看黄色视频| 国产精品 欧美亚洲| 色在线成人网| e午夜精品久久久久久久| 国产成人欧美| 99riav亚洲国产免费| 久久久久国产精品人妻aⅴ院 | 老熟妇乱子伦视频在线观看| 高清黄色对白视频在线免费看| 亚洲av片天天在线观看| 日韩熟女老妇一区二区性免费视频| 丰满人妻熟妇乱又伦精品不卡| 91成人精品电影| 国产欧美日韩综合在线一区二区| 老汉色∧v一级毛片| 日韩人妻精品一区2区三区| 亚洲专区字幕在线| 亚洲九九香蕉| 搡老熟女国产l中国老女人| 亚洲精品粉嫩美女一区| 两个人看的免费小视频| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 热99久久久久精品小说推荐| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 美女视频免费永久观看网站| 这个男人来自地球电影免费观看| 97人妻天天添夜夜摸| 亚洲精品中文字幕一二三四区| 精品国产国语对白av| 亚洲av美国av| 搡老乐熟女国产| 久久狼人影院| www.精华液| 免费黄频网站在线观看国产| 最新在线观看一区二区三区| 韩国精品一区二区三区| 麻豆国产av国片精品| 一边摸一边抽搐一进一出视频| 亚洲av成人av| 久久久久视频综合| 色综合婷婷激情| 午夜免费鲁丝| ponron亚洲| 69av精品久久久久久| 欧美中文综合在线视频| 亚洲中文字幕日韩| 19禁男女啪啪无遮挡网站| 亚洲精品乱久久久久久| 久久香蕉精品热| 90打野战视频偷拍视频| 亚洲成国产人片在线观看| 亚洲欧美日韩另类电影网站| 99re6热这里在线精品视频| 免费不卡黄色视频| 国产精品国产高清国产av | 日韩有码中文字幕| 精品国产国语对白av| 日韩精品免费视频一区二区三区| 嫁个100分男人电影在线观看| 黄色毛片三级朝国网站| 日韩成人在线观看一区二区三区| 日韩大码丰满熟妇| 男男h啪啪无遮挡| 日韩成人在线观看一区二区三区| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品中文字幕一二三四区| 亚洲午夜理论影院| 欧美日韩亚洲综合一区二区三区_| 91在线观看av| 日韩欧美免费精品| 欧美乱色亚洲激情| 免费女性裸体啪啪无遮挡网站| 91九色精品人成在线观看| 视频区欧美日本亚洲| 国产激情欧美一区二区| 91成年电影在线观看| 久久 成人 亚洲| 亚洲精品粉嫩美女一区| 少妇猛男粗大的猛烈进出视频| 欧美激情极品国产一区二区三区| 一级毛片高清免费大全| 天天添夜夜摸| 99久久99久久久精品蜜桃| 国产免费现黄频在线看| 热re99久久国产66热| av有码第一页| 亚洲国产欧美网| 精品久久久久久久毛片微露脸| 国产高清国产精品国产三级| 中文亚洲av片在线观看爽 | 黄色女人牲交| 久99久视频精品免费| av电影中文网址| 精品一品国产午夜福利视频| 久久九九热精品免费| 日韩免费高清中文字幕av| 欧美黄色片欧美黄色片| 久久久国产成人精品二区 | 超碰97精品在线观看| 在线看a的网站| 亚洲色图综合在线观看| 国产高清激情床上av| 精品一区二区三区四区五区乱码| 麻豆成人av在线观看| 日韩欧美免费精品| 免费在线观看日本一区| 国产av一区二区精品久久| 成人国语在线视频| 日本黄色视频三级网站网址 | 99久久人妻综合| av片东京热男人的天堂| 日本欧美视频一区| 黄色丝袜av网址大全| 久久天堂一区二区三区四区| 黄色a级毛片大全视频| 亚洲熟妇中文字幕五十中出 | 老熟妇仑乱视频hdxx| 99久久国产精品久久久| 看免费av毛片| 国产欧美日韩精品亚洲av| 窝窝影院91人妻| 久久久国产精品麻豆| 搡老熟女国产l中国老女人| 男女之事视频高清在线观看| 下体分泌物呈黄色| 高清毛片免费观看视频网站 | 纯流量卡能插随身wifi吗| 99热国产这里只有精品6| 午夜精品国产一区二区电影| 久久久精品国产亚洲av高清涩受| 亚洲美女黄片视频| 亚洲国产中文字幕在线视频| 在线免费观看的www视频| 黄色a级毛片大全视频| 国产午夜精品久久久久久| 亚洲av成人一区二区三| 国产伦人伦偷精品视频| 欧美日韩精品网址| 国产色视频综合| 日韩免费av在线播放| 免费不卡黄色视频| 午夜福利乱码中文字幕| 黄色a级毛片大全视频| 国产精品九九99| 久久精品国产综合久久久| 国产精品.久久久| 在线国产一区二区在线| 黑丝袜美女国产一区| 欧美性长视频在线观看| 国产黄色免费在线视频| 免费在线观看影片大全网站| 夜夜夜夜夜久久久久| 国产真人三级小视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 人妻一区二区av| 国产成人免费无遮挡视频| 人人妻人人澡人人爽人人夜夜| 精品久久久精品久久久| 欧美黄色片欧美黄色片| 香蕉国产在线看| 两人在一起打扑克的视频| 成人特级黄色片久久久久久久| 欧美日韩乱码在线| 亚洲欧美一区二区三区久久| av有码第一页| 真人做人爱边吃奶动态| 国产免费现黄频在线看| 精品卡一卡二卡四卡免费| 久久精品人人爽人人爽视色| 色综合婷婷激情| 天天添夜夜摸| 精品久久久精品久久久| 黄频高清免费视频| 又黄又爽又免费观看的视频| 亚洲黑人精品在线| 久热爱精品视频在线9| 啪啪无遮挡十八禁网站| 免费在线观看视频国产中文字幕亚洲| 黑丝袜美女国产一区| 又黄又粗又硬又大视频| 亚洲av电影在线进入| 看黄色毛片网站| 俄罗斯特黄特色一大片| 国产在视频线精品| 99精国产麻豆久久婷婷| 午夜福利影视在线免费观看| 91av网站免费观看| 老司机亚洲免费影院| 狠狠婷婷综合久久久久久88av| 人人妻人人添人人爽欧美一区卜| videos熟女内射| 精品国产一区二区久久| 国产欧美日韩一区二区三区在线| 国产成人一区二区三区免费视频网站| 丝瓜视频免费看黄片| 国产精品美女特级片免费视频播放器 | 中文字幕精品免费在线观看视频| 亚洲九九香蕉| 国产成+人综合+亚洲专区| 极品教师在线免费播放| 精品国产国语对白av| 久久久久久久午夜电影 | 欧美 亚洲 国产 日韩一| 午夜亚洲福利在线播放| 又黄又粗又硬又大视频| 免费女性裸体啪啪无遮挡网站| 香蕉久久夜色| 日韩欧美免费精品| 女人爽到高潮嗷嗷叫在线视频| 90打野战视频偷拍视频| 亚洲午夜理论影院| 国产精品自产拍在线观看55亚洲 | 天天添夜夜摸| 欧美乱码精品一区二区三区| 精品久久久久久久毛片微露脸| 精品国产国语对白av| 日韩视频一区二区在线观看| 看黄色毛片网站| 国产精品自产拍在线观看55亚洲 | 人人妻人人澡人人爽人人夜夜| 色播在线永久视频| 亚洲人成77777在线视频| 老司机深夜福利视频在线观看| 极品教师在线免费播放| 亚洲午夜理论影院| 午夜亚洲福利在线播放| 9191精品国产免费久久| 日韩熟女老妇一区二区性免费视频| 中文字幕av电影在线播放| 国产欧美亚洲国产| 中文字幕高清在线视频| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 电影成人av| 少妇 在线观看| 超色免费av| 欧美日韩黄片免| 天天操日日干夜夜撸| 国产精品一区二区精品视频观看| 欧美精品人与动牲交sv欧美| 成年人午夜在线观看视频| 少妇裸体淫交视频免费看高清 | 99久久人妻综合| 免费不卡黄色视频| 亚洲成人免费av在线播放| 亚洲成a人片在线一区二区| 99国产极品粉嫩在线观看| 免费在线观看亚洲国产| 国产在线一区二区三区精| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 精品久久久久久电影网| 少妇粗大呻吟视频| 在线观看免费午夜福利视频| 欧美+亚洲+日韩+国产| 亚洲欧美激情综合另类| 久久久国产欧美日韩av| 欧美日韩亚洲国产一区二区在线观看 | 91精品三级在线观看| av不卡在线播放| 一区在线观看完整版| 国产不卡一卡二| 国产成人免费无遮挡视频| 电影成人av| 亚洲精品在线观看二区| 成人18禁高潮啪啪吃奶动态图| 久久人妻福利社区极品人妻图片| 黄网站色视频无遮挡免费观看| 亚洲成国产人片在线观看| 十八禁网站免费在线| 欧美激情久久久久久爽电影 | 男女之事视频高清在线观看| 国产亚洲精品久久久久久毛片 | 欧美av亚洲av综合av国产av| 成人18禁高潮啪啪吃奶动态图| 午夜亚洲福利在线播放| 国产99白浆流出| av网站在线播放免费| 国产精品九九99| 午夜激情av网站| 99久久国产精品久久久| av欧美777| 99久久99久久久精品蜜桃| 大香蕉久久成人网| 精品一区二区三卡| 亚洲中文日韩欧美视频| 高清毛片免费观看视频网站 | 亚洲视频免费观看视频| 欧美乱色亚洲激情| 久久香蕉激情| 国产av一区二区精品久久| 午夜日韩欧美国产| 麻豆av在线久日| 久久国产乱子伦精品免费另类| 高清在线国产一区| 午夜福利在线免费观看网站| 老司机在亚洲福利影院| 日本wwww免费看| 人妻一区二区av| 欧美日韩福利视频一区二区| 久久久精品免费免费高清| 精品久久久精品久久久| 亚洲熟女精品中文字幕| 久热这里只有精品99| 好看av亚洲va欧美ⅴa在| 国产熟女午夜一区二区三区| 国产免费av片在线观看野外av| 亚洲欧美一区二区三区久久| 国产av又大| 亚洲国产看品久久| 黑人操中国人逼视频| 国产精品偷伦视频观看了| 日本五十路高清| 免费在线观看日本一区| 1024视频免费在线观看| 日韩成人在线观看一区二区三区| 丁香六月欧美| 久久精品成人免费网站|