• <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ì)
    與北大西洋接壤的北極海冰和年際氣候變化
    自线自在国产av| 一本久久中文字幕| 成人手机av| a级毛片在线看网站| 欧美日韩乱码在线| 久9热在线精品视频| 两人在一起打扑克的视频| 美女扒开内裤让男人捅视频| 极品人妻少妇av视频| 久久久久国产一级毛片高清牌| 婷婷丁香在线五月| 在线av久久热| 久久精品国产亚洲av高清一级| 黄色视频不卡| 在线免费观看的www视频| 亚洲精品久久成人aⅴ小说| 午夜影院日韩av| 国产高清视频在线播放一区| 亚洲精品中文字幕在线视频| 黑人巨大精品欧美一区二区mp4| 亚洲免费av在线视频| 久久性视频一级片| 18禁裸乳无遮挡免费网站照片 | 久久久国产成人精品二区| 亚洲一码二码三码区别大吗| 亚洲久久久国产精品| 精品人妻在线不人妻| 一区二区三区高清视频在线| 欧美乱色亚洲激情| 久久久久久亚洲精品国产蜜桃av| 亚洲精品粉嫩美女一区| 村上凉子中文字幕在线| 99久久久亚洲精品蜜臀av| 女同久久另类99精品国产91| 精品卡一卡二卡四卡免费| 两个人看的免费小视频| 99久久综合精品五月天人人| 后天国语完整版免费观看| 这个男人来自地球电影免费观看| 电影成人av| 国产一区二区三区在线臀色熟女| 91麻豆精品激情在线观看国产| 97人妻精品一区二区三区麻豆 | 国产97色在线日韩免费| 十分钟在线观看高清视频www| 日本 欧美在线| 最新在线观看一区二区三区| 91精品国产国语对白视频| 午夜激情av网站| 欧美一区二区精品小视频在线| av视频免费观看在线观看| 国内毛片毛片毛片毛片毛片| 国产欧美日韩综合在线一区二区| 精品第一国产精品| 韩国av一区二区三区四区| 在线观看免费视频日本深夜| 久久欧美精品欧美久久欧美| 国产精品精品国产色婷婷| 国产成人影院久久av| 国内精品久久久久久久电影| 国产99久久九九免费精品| av免费在线观看网站| 久久久久亚洲av毛片大全| 婷婷六月久久综合丁香| 一夜夜www| 久久香蕉激情| 成人三级做爰电影| 色婷婷久久久亚洲欧美| 涩涩av久久男人的天堂| 琪琪午夜伦伦电影理论片6080| 亚洲精品久久国产高清桃花| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲熟妇少妇任你| 国产av在哪里看| www.999成人在线观看| 日本免费a在线| 别揉我奶头~嗯~啊~动态视频| 性欧美人与动物交配| 免费在线观看日本一区| 动漫黄色视频在线观看| 午夜福利,免费看| 日本免费a在线| 男人的好看免费观看在线视频 | 人人妻人人爽人人添夜夜欢视频| 中国美女看黄片| 国产国语露脸激情在线看| 国产aⅴ精品一区二区三区波| 一本综合久久免费| 老熟妇仑乱视频hdxx| 丰满的人妻完整版| 欧美成人性av电影在线观看| 成年女人毛片免费观看观看9| 91老司机精品| 亚洲精品在线美女| 欧美黄色片欧美黄色片| 欧美日韩黄片免| 日本 av在线| 极品教师在线免费播放| 国产亚洲欧美98| 国产99久久九九免费精品| 成人精品一区二区免费| 亚洲美女黄片视频| 亚洲九九香蕉| 欧美色欧美亚洲另类二区 | 亚洲熟妇中文字幕五十中出| 制服诱惑二区| 在线播放国产精品三级| 手机成人av网站| 国语自产精品视频在线第100页| 亚洲国产毛片av蜜桃av| 成人免费观看视频高清| 国产成人av教育| 高潮久久久久久久久久久不卡| 色老头精品视频在线观看| 99国产精品一区二区三区| av中文乱码字幕在线| 女性生殖器流出的白浆| 久9热在线精品视频| 欧美大码av| 欧美激情久久久久久爽电影 | 99香蕉大伊视频| 国产精品永久免费网站| 曰老女人黄片| 不卡一级毛片| 亚洲少妇的诱惑av| 欧美午夜高清在线| 国产精品九九99| 嫩草影院精品99| 免费看美女性在线毛片视频| 女性生殖器流出的白浆| 久久精品国产清高在天天线| 搡老岳熟女国产| 99re在线观看精品视频| 两个人免费观看高清视频| 搡老岳熟女国产| 女人高潮潮喷娇喘18禁视频| 国产三级在线视频| 99re在线观看精品视频| 国产一卡二卡三卡精品| 亚洲男人天堂网一区| 精品乱码久久久久久99久播| 国产欧美日韩精品亚洲av| 51午夜福利影视在线观看| 国产麻豆69| 可以在线观看的亚洲视频| 涩涩av久久男人的天堂| 欧美 亚洲 国产 日韩一| 精品久久久精品久久久| 亚洲一区中文字幕在线| 欧美亚洲日本最大视频资源| 久久久久精品国产欧美久久久| 亚洲欧美精品综合久久99| √禁漫天堂资源中文www| 狂野欧美激情性xxxx| 国产精品 欧美亚洲| 日韩大码丰满熟妇| 久99久视频精品免费| 可以在线观看的亚洲视频| 一进一出抽搐gif免费好疼| 搡老熟女国产l中国老女人| 精品福利观看| 亚洲欧美日韩另类电影网站| 久久久久久久久免费视频了| 久久久久精品国产欧美久久久| 久久人妻av系列| 国产日韩一区二区三区精品不卡| 国产精品一区二区精品视频观看| 亚洲午夜精品一区,二区,三区| 亚洲电影在线观看av| 18禁国产床啪视频网站| 99riav亚洲国产免费| 一卡2卡三卡四卡精品乱码亚洲| 精品福利观看| 久久久国产欧美日韩av| 国产精品久久久久久人妻精品电影| 高清在线国产一区| 色在线成人网| 亚洲av熟女| 成人三级做爰电影| 成人特级黄色片久久久久久久| 亚洲少妇的诱惑av| 搡老妇女老女人老熟妇| 啦啦啦观看免费观看视频高清 | 午夜福利,免费看| 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频| 999久久久精品免费观看国产| 亚洲男人天堂网一区| tocl精华| 亚洲av片天天在线观看| 欧美日本亚洲视频在线播放| 91国产中文字幕| 午夜免费鲁丝| 男人舔女人的私密视频| 视频区欧美日本亚洲| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品综合一区在线观看 | 国产精品电影一区二区三区| 亚洲全国av大片| 久久久久久国产a免费观看| 国产麻豆69| 咕卡用的链子| 免费少妇av软件| 亚洲成av人片免费观看| 精品久久久精品久久久| 我的亚洲天堂| 亚洲全国av大片| 国产av精品麻豆| www国产在线视频色| 国产激情欧美一区二区| 波多野结衣巨乳人妻| 一级,二级,三级黄色视频| www.www免费av| 99香蕉大伊视频| 国产亚洲欧美精品永久| 啦啦啦韩国在线观看视频| 免费在线观看视频国产中文字幕亚洲| 亚洲电影在线观看av| 亚洲三区欧美一区| 日韩免费av在线播放| 国产精品一区二区三区四区久久 | 黄片大片在线免费观看| 国产区一区二久久| svipshipincom国产片| 日日夜夜操网爽| 激情在线观看视频在线高清| 亚洲第一av免费看| 禁无遮挡网站| 看黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 免费看美女性在线毛片视频| 窝窝影院91人妻| 国产精品久久久久久亚洲av鲁大| 男女之事视频高清在线观看| 九色国产91popny在线| 国产成人欧美| 男女之事视频高清在线观看| svipshipincom国产片| 不卡av一区二区三区| 欧美老熟妇乱子伦牲交| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品免费视频内射| 99在线人妻在线中文字幕| 大码成人一级视频| 欧美精品亚洲一区二区| 亚洲av电影在线进入| 中文字幕色久视频| 久热爱精品视频在线9| 欧美黑人欧美精品刺激| 免费在线观看亚洲国产| 精品不卡国产一区二区三区| 国产成人精品在线电影| 日日夜夜操网爽| 亚洲欧美日韩无卡精品| 18禁裸乳无遮挡免费网站照片 | 午夜久久久在线观看| 在线十欧美十亚洲十日本专区| 国产精品久久久人人做人人爽| 免费搜索国产男女视频| 久久影院123| 黄网站色视频无遮挡免费观看| 51午夜福利影视在线观看| av天堂久久9| 亚洲熟女毛片儿| 亚洲第一青青草原| 可以在线观看的亚洲视频| 国产伦一二天堂av在线观看| 国产高清激情床上av| 亚洲人成77777在线视频| 国产精品亚洲av一区麻豆| www日本在线高清视频| 少妇熟女aⅴ在线视频| avwww免费| 国产精品 国内视频| 最新美女视频免费是黄的| 欧美亚洲日本最大视频资源| 亚洲七黄色美女视频| av天堂久久9| av天堂在线播放| 日韩欧美免费精品| 黄色丝袜av网址大全| 老熟妇仑乱视频hdxx| 51午夜福利影视在线观看| 久久影院123| 神马国产精品三级电影在线观看 | 精品日产1卡2卡| 国产乱人伦免费视频| 丝袜美腿诱惑在线| 亚洲男人的天堂狠狠| 精品电影一区二区在线| 露出奶头的视频| 校园春色视频在线观看| 一区二区三区高清视频在线| 国产精品 国内视频| 性欧美人与动物交配| 国语自产精品视频在线第100页| 免费高清在线观看日韩| 国产人伦9x9x在线观看| 久久国产乱子伦精品免费另类| 最近最新免费中文字幕在线| 国产成人免费无遮挡视频| 两个人免费观看高清视频| 99精品久久久久人妻精品| 中文字幕精品免费在线观看视频| 88av欧美| 亚洲男人的天堂狠狠| 久久久久久人人人人人| 搡老妇女老女人老熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | www.熟女人妻精品国产| 成人国产一区最新在线观看| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 国产av一区二区精品久久| 午夜免费激情av| 国产精品亚洲美女久久久| 亚洲av美国av| 夜夜躁狠狠躁天天躁| 在线永久观看黄色视频| 桃色一区二区三区在线观看| 精品少妇一区二区三区视频日本电影| 午夜福利在线观看吧| 看黄色毛片网站| 国产精品久久久人人做人人爽| 非洲黑人性xxxx精品又粗又长| av免费在线观看网站| 久久香蕉国产精品| 91精品国产国语对白视频| 国产主播在线观看一区二区| 国内久久婷婷六月综合欲色啪| 国产xxxxx性猛交| 岛国在线观看网站| 两个人视频免费观看高清| 欧美最黄视频在线播放免费| 美女午夜性视频免费| 欧美日韩福利视频一区二区| 精品国产超薄肉色丝袜足j| а√天堂www在线а√下载| 欧美中文综合在线视频| 级片在线观看| 伊人久久大香线蕉亚洲五| 一级片免费观看大全| 中文字幕av电影在线播放| 国产主播在线观看一区二区| 久久久国产成人精品二区| 精品国产国语对白av| 亚洲成人精品中文字幕电影| 一级作爱视频免费观看| 老司机在亚洲福利影院| 99国产综合亚洲精品| 天堂影院成人在线观看| 最好的美女福利视频网| 午夜精品久久久久久毛片777| 岛国在线观看网站| 亚洲精品粉嫩美女一区| 美女扒开内裤让男人捅视频| 久久婷婷人人爽人人干人人爱 | 久久婷婷成人综合色麻豆| 免费搜索国产男女视频| 又紧又爽又黄一区二区| 免费搜索国产男女视频| 精品少妇一区二区三区视频日本电影| 中国美女看黄片| 又紧又爽又黄一区二区| 国产色视频综合| 在线观看www视频免费| 国产成人影院久久av| 女同久久另类99精品国产91| 99riav亚洲国产免费| 亚洲熟女毛片儿| 午夜福利视频1000在线观看 | 国产熟女午夜一区二区三区| 天天躁夜夜躁狠狠躁躁| 午夜老司机福利片| 涩涩av久久男人的天堂| 桃红色精品国产亚洲av| 亚洲中文字幕日韩| 免费人成视频x8x8入口观看| 又黄又粗又硬又大视频| 久久久久久大精品| 9191精品国产免费久久| 亚洲精品国产精品久久久不卡| 久久伊人香网站| 女人被躁到高潮嗷嗷叫费观| 亚洲国产欧美一区二区综合| 亚洲精品av麻豆狂野| 久久精品国产清高在天天线| 美女扒开内裤让男人捅视频| 免费av毛片视频| 欧美一级毛片孕妇| 女人高潮潮喷娇喘18禁视频| 1024视频免费在线观看| 亚洲熟女毛片儿| 亚洲黑人精品在线| 欧美性长视频在线观看| 久久亚洲精品不卡| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久av美女十八| av视频在线观看入口| 99国产精品99久久久久| 国产真人三级小视频在线观看| 久久九九热精品免费| 久久精品亚洲精品国产色婷小说| 国产精品九九99| 久久热在线av| 一级a爱片免费观看的视频| 69精品国产乱码久久久| 国产亚洲欧美98| 日韩精品中文字幕看吧| 免费人成视频x8x8入口观看| av免费在线观看网站| 午夜免费激情av| 老汉色∧v一级毛片| 性色av乱码一区二区三区2| 亚洲国产精品合色在线| 亚洲人成电影观看| 9热在线视频观看99| 1024视频免费在线观看| 老司机在亚洲福利影院| 午夜亚洲福利在线播放| 一级作爱视频免费观看| 丝袜在线中文字幕| 香蕉国产在线看| 精品电影一区二区在线| 给我免费播放毛片高清在线观看| 少妇粗大呻吟视频| 欧美在线一区亚洲| 搞女人的毛片| 国产男靠女视频免费网站| 久久草成人影院| 国产成人欧美| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 波多野结衣av一区二区av| 在线观看免费日韩欧美大片| 自线自在国产av| 成人国语在线视频| 色综合亚洲欧美另类图片| 国产av一区在线观看免费| 99国产精品免费福利视频| 亚洲中文日韩欧美视频| 99国产综合亚洲精品| 成人国产综合亚洲| 男人的好看免费观看在线视频 | 男男h啪啪无遮挡| 亚洲中文av在线| 中文字幕精品免费在线观看视频| 18美女黄网站色大片免费观看| 久热爱精品视频在线9| 久久性视频一级片| www.自偷自拍.com| 色av中文字幕| 国产精品一区二区三区四区久久 | av免费在线观看网站| 国产亚洲av嫩草精品影院| 久久 成人 亚洲| 香蕉丝袜av| 丁香六月欧美| 女性生殖器流出的白浆| 国产99久久九九免费精品| 国产亚洲av高清不卡| 久久香蕉国产精品| 青草久久国产| 性欧美人与动物交配| 久久欧美精品欧美久久欧美| 精品欧美国产一区二区三| 不卡av一区二区三区| 岛国视频午夜一区免费看| 在线观看日韩欧美| 色老头精品视频在线观看| 国产精品一区二区三区四区久久 | 一边摸一边抽搐一进一出视频| 99在线人妻在线中文字幕| 欧美中文日本在线观看视频| 女同久久另类99精品国产91| 天天添夜夜摸| 久久人妻av系列| 村上凉子中文字幕在线| 99riav亚洲国产免费| 在线观看舔阴道视频| 久久婷婷人人爽人人干人人爱 | 日韩国内少妇激情av| 88av欧美| 精品久久蜜臀av无| 在线观看免费视频日本深夜| 色综合站精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕人妻熟女乱码| 日韩大尺度精品在线看网址 | 国产91精品成人一区二区三区| 国内毛片毛片毛片毛片毛片| 91大片在线观看| 波多野结衣高清无吗| 久久人妻福利社区极品人妻图片| 亚洲av电影不卡..在线观看| 中文亚洲av片在线观看爽| 欧美黑人欧美精品刺激| 国产成人欧美在线观看| 亚洲成人免费电影在线观看| 欧美激情高清一区二区三区| 亚洲 欧美一区二区三区| 91精品国产国语对白视频| 精品国产一区二区三区四区第35| 精品熟女少妇八av免费久了| 国产精品野战在线观看| 一区福利在线观看| 亚洲人成伊人成综合网2020| 在线观看66精品国产| 国产一区二区三区在线臀色熟女| 法律面前人人平等表现在哪些方面| 成人av一区二区三区在线看| 亚洲中文av在线| 多毛熟女@视频| 亚洲美女黄片视频| 国产成人欧美在线观看| 精品人妻1区二区| 免费搜索国产男女视频| 高潮久久久久久久久久久不卡| 巨乳人妻的诱惑在线观看| 禁无遮挡网站| 99国产精品99久久久久| 乱人伦中国视频| 国产亚洲欧美在线一区二区| 校园春色视频在线观看| 免费一级毛片在线播放高清视频 | 色播亚洲综合网| 欧美成人午夜精品| 一级,二级,三级黄色视频| 国产精品98久久久久久宅男小说| 男女床上黄色一级片免费看| 成年女人毛片免费观看观看9| 黑人巨大精品欧美一区二区mp4| 久久精品人人爽人人爽视色| svipshipincom国产片| 日本五十路高清| 国产成人影院久久av| 老司机深夜福利视频在线观看| 欧美日本中文国产一区发布| 亚洲国产精品久久男人天堂| 久久久久亚洲av毛片大全| 久久久久久亚洲精品国产蜜桃av| 97人妻精品一区二区三区麻豆 | 999久久久精品免费观看国产| 亚洲欧美一区二区三区黑人| 久久国产亚洲av麻豆专区| 国产91精品成人一区二区三区| www国产在线视频色| 91字幕亚洲| 男人的好看免费观看在线视频 | 香蕉丝袜av| 欧美日本视频| 亚洲无线在线观看| 成人永久免费在线观看视频| 99国产精品一区二区蜜桃av| 久久精品亚洲精品国产色婷小说| 一夜夜www| 国语自产精品视频在线第100页| 久久久久九九精品影院| 侵犯人妻中文字幕一二三四区| 日韩一卡2卡3卡4卡2021年| 久久国产精品男人的天堂亚洲| 国产激情久久老熟女| 黄色女人牲交| 国产熟女xx| av视频在线观看入口| 亚洲欧美日韩另类电影网站| 丝袜人妻中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 人妻久久中文字幕网| 久久精品国产亚洲av香蕉五月| 黑人操中国人逼视频| 91国产中文字幕| 国产一区在线观看成人免费| 窝窝影院91人妻| 中国美女看黄片| 在线视频色国产色| 夜夜躁狠狠躁天天躁| 国产欧美日韩精品亚洲av| 亚洲国产日韩欧美精品在线观看 | 久久欧美精品欧美久久欧美| 国产国语露脸激情在线看| 午夜免费观看网址| 99精品在免费线老司机午夜| 叶爱在线成人免费视频播放| 亚洲自偷自拍图片 自拍| 夜夜爽天天搞| 亚洲一区二区三区色噜噜| 中文字幕人妻熟女乱码| 在线永久观看黄色视频| 免费在线观看日本一区| 纯流量卡能插随身wifi吗| 1024视频免费在线观看| 美女扒开内裤让男人捅视频| 欧美绝顶高潮抽搐喷水| 亚洲精品美女久久av网站| 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 中文字幕久久专区| 国产欧美日韩一区二区三区在线| 久久精品国产99精品国产亚洲性色 | 成人欧美大片| 久久午夜综合久久蜜桃| 国产精品一区二区精品视频观看| 动漫黄色视频在线观看| 国产精品久久久久久人妻精品电影| 人妻丰满熟妇av一区二区三区| 国产男靠女视频免费网站| 美女 人体艺术 gogo| 纯流量卡能插随身wifi吗| 18禁裸乳无遮挡免费网站照片 | 丝袜人妻中文字幕| 国产精品1区2区在线观看.| 69av精品久久久久久|