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

    南極半島巴赫冰架冰面融水動(dòng)態(tài)遙感監(jiān)測(cè)

    2022-06-11 01:30:26梁相安張聞松李雅陸瑤楊康
    極地研究 2022年2期
    關(guān)鍵詞:冰架融水冰蓋

    梁相安 張聞松 李雅 陸瑤 楊康,2,3

    (1 南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023;2 江蘇省地理信息技術(shù)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023;3 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(珠海),廣東 珠海 519080)

    0 引言

    南極冰蓋物質(zhì)損失是造成全球海平面上升的重要原因之一,研究南極冰蓋物質(zhì)平衡有助于提高海平面上升預(yù)測(cè)精度[1]。近年來,已有研究通過衛(wèi)星遙感影像在南極冰蓋表面觀測(cè)到了規(guī)模龐大形態(tài)復(fù)雜的冰面融水,這些冰面融水會(huì)通過直接或間接的方式影響冰架穩(wěn)定性,進(jìn)而影響南極冰蓋物質(zhì)平衡[1-4]。因此,研究南極冰蓋冰面融水的動(dòng)態(tài)變化,有助于分析冰架穩(wěn)定性,從而更準(zhǔn)確地估計(jì)南極冰蓋物質(zhì)平衡[1]。

    目前,已有研究開展了南極冰蓋典型消融區(qū)冰面融水遙感監(jiān)測(cè)研究。例如,Dell 等[5]利用Landsat-8 和Sentinel-2 影像監(jiān)測(cè)了Nivilisen 冰架冰面融水的動(dòng)態(tài)變化,發(fā)現(xiàn)在消融季內(nèi),大量融水通過冰面或冰下通道向冰架前端輸送;Bell 等[3]通過Landsat-8 OLI 影像發(fā)現(xiàn)Nansen 冰架的冰面融水通過瀑布(waterfall)和冰坑(doline)快速排出冰架,有利于增強(qiáng)冰架穩(wěn)定性;Kingslake 等[4]通過LandsatASTERWorldView 衛(wèi)星影像和航空相片展示了南極冰蓋多個(gè)典型消融區(qū)的冰面水系,指出南極冰蓋表面融水的分布范圍及影響遠(yuǎn)超過去的預(yù)測(cè);璩榆桐等[6]利用Landsat-8 影像實(shí)現(xiàn)了東南極達(dá)爾克冰川季節(jié)性消融信息提取,發(fā)現(xiàn)消融最早出現(xiàn)在地形下降劇烈的觸地線附近。

    南極半島消融劇烈,每年消融期大量融水在冰面存儲(chǔ)與輸送,是研究冰面融水與冰架穩(wěn)定性相互作用的典型區(qū)域。例如,Banwell 等[7]通過Landsat-7 ETM+和ASTER 影像監(jiān)測(cè)了Larsen B冰架崩解前發(fā)生的冰面湖快速排水(fast drainage)事件,發(fā)現(xiàn)融水的存儲(chǔ)和輸送會(huì)壓迫冰架使其發(fā)生彎曲和斷裂;Tuckett 等[8]發(fā)現(xiàn)南極半島地區(qū)的融水會(huì)進(jìn)入冰架內(nèi)部,造成冰流速增加,進(jìn)而減弱冰架穩(wěn)定性;Banwell 等[9]通過微波遙感數(shù)據(jù)和光學(xué)影像數(shù)據(jù)監(jiān)測(cè)了George VI 冰架表面融水長時(shí)間序列的動(dòng)態(tài)變化,發(fā)現(xiàn)2019—2020年消融季,George VI 冰面消融的持續(xù)時(shí)間和范圍急劇增加。然而,目前關(guān)于南極半島冰面融水的研究主要聚焦于南極半島東岸Larsen A/B/C 冰架,南極半島西岸冰面融水的時(shí)空動(dòng)態(tài)尚不明確。此外,這些研究大多分析冰面融水在特定消融時(shí)刻或一個(gè)消融期內(nèi)的分布狀況,對(duì)冰面融水長時(shí)序時(shí)空動(dòng)態(tài)及其驅(qū)動(dòng)機(jī)制的分析較少。

    本研究旨在分析南極半島巴赫冰架(Bach Ice Shelf)上,冰面融水空間分布特征及其長時(shí)序動(dòng)態(tài)變化規(guī)律,并探究遙感反演的冰面融水體積與MERRA-2 模型建模融水徑流量的關(guān)系。

    1 研究區(qū)域與研究數(shù)據(jù)

    1.1 研究區(qū)域

    巴赫冰架(72°W,72°S)位于南極半島西岸的亞歷山大島(Alexander Island,70°W—73°W,71°S—72°S)的南部,長約72 km。巴赫冰架冰面消融十分劇烈,大量融水在冰架表面存儲(chǔ)和輸送,研究其冰面融水的空間分布及動(dòng)態(tài)變化特征,有助于提升對(duì)南極半島冰面融水形成與演化規(guī)律的理解。根據(jù)MEaSUREs Bed Machine Antarctica(Version 1)數(shù)據(jù)集2019年11月5日的數(shù)據(jù),在研究區(qū)域內(nèi),冰架厚度最大處為821 m,最小處為0 m(裸巖),平均冰架厚度 355 m(標(biāo)準(zhǔn)差為172 m)。除少數(shù)裸露的巖石外,該地區(qū)幾乎全部被永久性積雪和冰川所覆蓋。巴赫冰架消融具有季節(jié)性,每年夏季(11月至次年4月)融水通過冰面水系輸送并存儲(chǔ)至冰面湖中,這些冰面湖主要集中在冰架的東北部,因此本文選擇巴赫冰架的東北部(面積約1800 km2)作為研究區(qū)域(圖1)。

    圖1 研究區(qū)域。a)研究區(qū)域位置;b)研究區(qū)域Landsat-8 假彩色合成影像。其中b)圖為a)中紅框區(qū)域,其成像時(shí)間為2020年1月7日;RGB:band 5(近紅外),band 4(紅光),band 3(綠光)Fig.1.The study area.a) the location of the study area;b) Landsat-8 false color composite image of the study area.The red box in a) depicts the area as shown in b).The image was acquired on 7 January 2020;RGB:band 5(NIR),band 4(Red),band 3(Green)

    1.2 研究數(shù)據(jù)

    本文使用的遙感數(shù)據(jù)是Landsat-8 衛(wèi)星搭載的陸地成像儀(Operational Land Imager,OLI)影像,共有9 個(gè)波段,其中2~4 波段分別為藍(lán)光綠光和紅光波段,5 波段為近紅外波段,空間分辨率均為30 m,影像數(shù)據(jù)在研究區(qū)內(nèi)時(shí)空覆蓋較完整。本文選取了2013—2020年的21 期高質(zhì)量(16 期云覆蓋<20%,5 期云覆蓋<40%)Landsat-8 OLI 影像(表1),所有影像通過美國地質(zhì)調(diào)查局(United States Geological Survey,USGS)官方網(wǎng)站(https://www.usgs.gov/)獲得。

    表1 研究中使用的Landsat-8 影像列表Table 1.List of Landsat-8 imagery used in this study

    MERRA-2(Modern-Era Retrospective analysis for Research and Applications,Version 2)是美國國家航空航天局的全球模擬與同化辦公室(Global Modeling and Assimilation Office,GMAO)開發(fā)的全球大氣再分析數(shù)據(jù)集。MERRA-2 可以模擬影響南極冰蓋表面物質(zhì)平衡的日徑流量(runoff)等變量[10],空間分辨率0.5°×0.625°,目前可獲取1980—2020年的數(shù)據(jù)[11]。其中,日徑流量表示當(dāng)天冰面總消融量中可轉(zhuǎn)化為冰面融水徑流的部分,輸出單位為mm·d-1。研究將日徑流量乘以研究區(qū)面積,得到研究區(qū)內(nèi)的日徑流總量(m3·d-1),累積每個(gè)消融季從消融開始到某一指定日期的日徑流總量,可以得到冰面融水總體積模擬結(jié)果。

    本文采用南極冰蓋高空間分辨率數(shù)字高程模型(Reference Elevation Model Antarctica,REMA)定量化描述冰面地形對(duì)冰面融水分布的影響。該產(chǎn)品融合了 2010—2018年獲取的 WorldView-1/2/3 和GeoEye-1 高空間分辨率衛(wèi)星遙感影像,可覆蓋 98%的南極冰蓋區(qū)域,覆蓋范圍延伸到88°S,空間分辨率為8 m,垂直精度為1 m[12]。

    2 研究方法

    2.1 冰面融水范圍遙感信息提取方法

    研究利用歸一化水體指數(shù)(Normalized Difference Water Index,NDWI)提取冰面融水[13]。NDWI利用水體在近紅外波段反射率極低而在可見光波段反射率較高的光譜特征,能夠有效區(qū)分冰面水體與背景冰雪[14],計(jì)算公式如下:

    其中,Green和NIR分別表示綠光和近紅外波段大氣層頂反射率。本文使用NDWI單閾值分割法進(jìn)行冰面融水的水體范圍提取,以0.01 為間隔遍歷0.1 到0.3 內(nèi)所有值后,選取提取效果最好的0.2作為統(tǒng)一閾值。該NDWI閾值與東南極Nansen冰架冰面湖提取(0.25)[3]環(huán)南極冰蓋多個(gè)典型消融區(qū)冰面湖提取(0.19)[15]等研究使用的閾值相近。最終得到21 期二值化水體提取結(jié)果,其中1代表水體像元,0 代表非水體像元,統(tǒng)計(jì)水體像元數(shù)目后,將水體像元數(shù)目與像元面積相乘得到研究區(qū)內(nèi)的水體總面積。進(jìn)一步將21 期冰面融水提取結(jié)果疊加,得到2013—2020年冰面融水最大空間分布范圍。

    2.2 冰面融水體積遙感反演方法

    研究通過光學(xué)影像經(jīng)驗(yàn)公式法反演融水深度。首先,將原始影像中的DN(Digital Number)值轉(zhuǎn)化為大氣層頂反射率(Top of Atmosphere,TOA),然后通過南極冰面湖深度反演經(jīng)驗(yàn)公式反演水深[16],計(jì)算公式如下:

    其中,z為反演得到的冰面融水水深,Blue與Red分別為藍(lán)光與紅光波段的大氣層頂反射率,水深反演均方根誤差(RMSE)為0.38 m[17]。反演得到的水深柵格中,每個(gè)像元的值代表了當(dāng)前位置的水深,單位為m。統(tǒng)計(jì)研究區(qū)內(nèi)水深柵格的像元值之和,并與單個(gè)像元面積(900 m2)相乘,即得到冰面融水總體積。

    3 結(jié)果

    3.1 冰面融水空間分布

    巴赫冰架冰面融水的形態(tài)呈現(xiàn)為冰面湖(0.3±0.2 km2)與河寬較寬(70±30 m)的冰面河。每年12月,融水分布在觸地線附近;在消融較劇烈的1月和2月,融水?dāng)U張至觸地線以外距離較遠(yuǎn)的地區(qū)(圖2)。大部分(>95%)冰面融水分布在低海拔(<40 m)地區(qū)和坡度較小(<1°)的地區(qū)(圖3)。

    圖3 巴赫冰架2013—2020年冰面融水最大空間分布。a) 水體掩膜;b) 水體像元處海拔;c) 水體像元處與觸地線距離;d) 水體像元處坡度Fig.3.Maximum extent of surface meltwater in Bach Ice Shelf during 2013-2020.a) water mask;b) elevation of water pixels;c) distance from water pixels to grounding line;d) slope of water pixels

    冰面融水出現(xiàn)的最高海拔約為120 m,位于觸地線以內(nèi)約5 km 處;最低海拔為6 m,位于觸地線以外16 km 處。此外,91.7%的冰面融水分布在觸地線以外的冰架上,且冰架上>90%的融水與觸地線距離較近(<7 km)。在觸地線另一側(cè)的內(nèi)陸冰蓋表面,冰面融水主要存儲(chǔ)在數(shù)量有限(約30 個(gè))的冰面湖中,彼此孤立并不連通,最遠(yuǎn)的冰面湖與觸地線相距13.5 km(圖23)。

    圖2 巴赫冰架2013—2020年冰面融水提取結(jié)果Fig.2.Supraglacial meltwater extraction results of Bach Ice Shelf during 2013-2020

    本研究利用人工數(shù)字化結(jié)果驗(yàn)證水體遙感提取精度。圖4a4b 是從2017年1月15日Landsat-8影像中選取的兩個(gè)驗(yàn)證區(qū)域,圖4c4d 是從2020年1月17日Landsat-8 影像中選取的兩個(gè)驗(yàn)證區(qū)域。結(jié)果表明本文方法提取的水體邊界與人工數(shù)字化結(jié)果較為一致(R2=0.98),冰面融水遙感提取結(jié)果可靠(圖4)。

    圖4 驗(yàn)證樣本和提取結(jié)果。a)-d) 為驗(yàn)證區(qū)域局部放大圖;e) 精度檢驗(yàn)結(jié)果Fig.4.Validation sample and extraction results.a)-d) partial enlarged views;e) result of accuracy assessment

    3.2 冰面融水年際變化

    本研究利用MERRA-2 模型模擬了研究區(qū)內(nèi)冰面日徑流量,并與遙感影像監(jiān)測(cè)的融水面積和體積進(jìn)行了對(duì)比(圖5)。結(jié)果表明:(1)2013—2020年,巴赫冰架冰面消融強(qiáng)度總體呈上升趨勢(shì),其中2013—2016年消融強(qiáng)度相差不大,2016—2018年消融強(qiáng)度迅速上升,2018—2019年消融強(qiáng)度較弱,2019—2020年消融強(qiáng)度達(dá)到頂峰;(2)在每個(gè)消融季內(nèi),冰面消融開始于12月下旬,于次年1月中旬至下旬達(dá)到峰值,之后開始減弱,直到次年2月中旬結(jié)束,其中2019—2020年消融持續(xù)時(shí)間較長,融水于12月下旬出現(xiàn),直至次年3月中旬結(jié)束;(3)模型模擬的2018年與2020年冰面徑流量峰值與遙感反演的冰面融水面積及體積對(duì)應(yīng)較好,其他年份可能受限于冰面融水遙感觀測(cè)頻次,并未發(fā)現(xiàn)這種對(duì)應(yīng)關(guān)系。

    圖5 冰面融水動(dòng)態(tài)變化。a)模型估計(jì)的日徑流量;b)遙感觀測(cè)融水面積;c)遙感觀測(cè)融水體積Fig.5.Seasonal evolution of surface meltwater.a) model estimated runoff;b) remotely sensed area of meltwater;c) remotely sensed volume of meltwater

    研究進(jìn)一步對(duì)比了遙感反演的融水體積與模型模擬的冰面融水徑流量(圖 6)。結(jié)果表明,MERRA-2 模型高估了融水徑流總量,模型模擬的融水徑流量總體比遙感反演的融水體積高70.7%。其中,12月的樣本點(diǎn)(紅)大多分布于1∶1 線下方;1 份的樣本點(diǎn)(綠)大多分布于1∶1線上方;代表2月融水體積的點(diǎn)(藍(lán))全部分布于1∶1 線上方,且與1∶1 線距離較遠(yuǎn),說明從12月到次年2月,MERRA-2 模型對(duì)冰面融水徑流量的高估逐月增加。

    圖6 模型模擬冰面融水徑流量與遙感反演融水體積對(duì)比Fig.6.Comparison of model estimated volume and remotely sensed meltwater volume

    3.3 冰面融水年內(nèi)變化

    本文進(jìn)一步分析了單個(gè)消融季內(nèi)冰面融水的動(dòng)態(tài)變化(圖7)。消融初期,融水在觸地線附近生成,之后向低高程地區(qū)擴(kuò)展。隨著消融的增強(qiáng),研究區(qū)內(nèi)的冰面湖通過冰面水系相互連通,融水積累匯集成更大的冰面湖,這種現(xiàn)象在研究區(qū)內(nèi)融水總面積大于20 km2時(shí)出現(xiàn)。圖7 展示了消融增強(qiáng)階段,局部融水生成與演化的細(xì)節(jié),其中不同顏色代表了一個(gè)典型消融季(2017—2018年)內(nèi)不同日期的水體提取結(jié)果,可以看到冰面融水向下游擴(kuò)張的連續(xù)過程,新的冰面湖在上游冰面湖附近生成,上下游冰面湖通過冰面水系相互連通,上游冰面湖中的融水通過冰面河道向下游輸送,最終匯入一個(gè)大冰面湖(6.6 km2)中,并在其中存儲(chǔ),不再向冰架前端輸送,冰面融水并未進(jìn)入大洋(圖7)。

    4 討論

    4.1 冰面融水形成與演化的驅(qū)動(dòng)機(jī)制分析

    巴赫冰架冰面融水的動(dòng)態(tài)變化與東南極冰蓋沙克爾頓冰架(Shackleton Ice Shelf)[13]Nivilisen冰架[5]冰面融水相似。在消融初期,冰面融水在觸地線附近生成,而不是在高程更低的冰架前端生成,這可能與南極地區(qū)的下降風(fēng)(katabatic wind)有關(guān)。觸地線位于高海拔地區(qū)(冰蓋)與低海拔地區(qū)(冰架)交界處,高壓氣流受重力作用加速下滑[3-4,17-18],形成強(qiáng)勁的下降風(fēng),吹散觸地線附近的表面積雪,形成反照率較低的藍(lán)冰,加快觸地線附近地區(qū)的消融,從而產(chǎn)生冰面融水[19]。進(jìn)一步,冰面融水的反照率顯著低于冰雪,會(huì)加速周邊區(qū)域冰面融化,形成正反饋機(jī)制,產(chǎn)生的大量融水在地形作用下向下游輸送,冰面融水持續(xù)擴(kuò)張[20-21]。

    融水進(jìn)入冰面湖后,有兩種情景:(1)如果該冰面湖與冰面河道連通,則融水沿河道向下游輸送;(2)如果該冰面湖不與冰面河道連通,則融水將在冰面湖中存儲(chǔ),填滿冰面湖的湖盆后,融水溢出并在重力作用下向下游輸送[22-26]。由于研究區(qū)內(nèi)融水主要分布在坡度較緩(<1°)的地區(qū),上下游高程差小,可轉(zhuǎn)化為水流動(dòng)能的重力勢(shì)能有限。因此,在連接上下游的冰面水系中,融水流速慢徑流量小,在融水向下游擴(kuò)張時(shí),上游的冰面湖縮小不明顯,大部分上游冰面湖面積保持不變或略有上升(圖7)。由此判斷,通過冰面水系輸往下游的融水總量小于或等于上游冰面湖自身位置冰面消融的增量。

    圖7 2017—2018年冰面融水動(dòng)態(tài)變化。a) 空間上的演化過程;b) 體積和面積的動(dòng)態(tài)變化Fig.7.Intra-seasonal(2017-2018) evolution of surface meltwater.a) spatial evolution;b) variation of volume and area

    4.2 MERRA-2 模擬結(jié)果誤差分析

    MERRA-2 模型對(duì)融水模擬的誤差主要有三個(gè)來源。(1)MERRA-2 模型并未考慮融水的蒸發(fā)過程,存儲(chǔ)在冰面湖中的融水,除了因重凍結(jié)作用造成的損失外,還有一部分會(huì)通過蒸發(fā)作用進(jìn)入空氣中[27],MERRA-2 模型忽略了這部分融水損失[11],該誤差會(huì)隨時(shí)間不斷累積,因此從12月到次年2月,模型對(duì)融水的高估呈增加趨勢(shì)(圖6)。(2)MERRA-2 的空間分辨率為0.5°×0.625°(約56 km ×44 km),研究區(qū)的面積較小,只包含1~2 個(gè)MERRA-2 網(wǎng)格單元,在消融初期(12月),由于空間分辨率過粗,MERRA-2 模型可能并沒有捕捉到冰面的消融,導(dǎo)致模型估計(jì)值小于實(shí)際冰面融水總量(圖6)。(3)模型模擬和遙感反演的融水體積是在同一范圍中進(jìn)行統(tǒng)計(jì)的,顯然該統(tǒng)計(jì)范圍的選擇也會(huì)影響結(jié)果,理想的統(tǒng)計(jì)范圍應(yīng)該包含研究區(qū)內(nèi)所有融水及其上游的匯流區(qū)域[28]。然而,本文使用的統(tǒng)計(jì)范圍不僅包含了融水上游的匯流區(qū)域,還包含了部分處于融水下游的區(qū)域,這些區(qū)域產(chǎn)生的消融不會(huì)為研究區(qū)內(nèi)的冰面湖提供融水,但該區(qū)域的冰面消融量也被計(jì)入了模型模擬的融水總量中,造成模型模擬結(jié)果偏大。在未來的對(duì)比中,應(yīng)該使用空間分辨率更高(1~10 km)的區(qū)域氣候模型[29-30],并且考慮冰面融水的蒸發(fā)過程[31]。

    5 結(jié)論

    本文結(jié)合Landsat-8 影像高分辨率數(shù)字高程產(chǎn)品REMA 和再分析數(shù)據(jù)集MERRA-2,分析了2013—2020年南極半島巴赫冰架冰面融水動(dòng)態(tài)變化,對(duì)比了MERRA-2 模型模擬的融水徑流量與遙感反演的融水體積,主要結(jié)論如下。

    1.>90%的冰面融水分布在觸地線以外的浮冰上,與觸地線距離較近(<7 km),融水通過冰面水系在不同冰面湖之間輸送,并未到達(dá)冰架前端進(jìn)入大洋。

    2.2013—2020年冰面融水面積與體積呈上升趨勢(shì),每年 12月中旬融水出現(xiàn)在觸地線周圍,并開始向冰架前端擴(kuò)張,次年1月中下旬達(dá)到最大值,之后冰面融水開始減少,最終完全消失。

    3.冰面融水體積小于MERRA-2 模型模擬的融水徑流量,MERRA-2 模型過高估計(jì)了冰面實(shí)際融水徑流量。

    由于本研究區(qū)可用Landsat-8 影像較少,尚不能連續(xù)監(jiān)測(cè)單個(gè)消融季內(nèi)冰面融水的動(dòng)態(tài)變化。后續(xù)研究應(yīng)融合多源遙感影像,增加影像的時(shí)間覆蓋頻率,實(shí)現(xiàn)更精細(xì)的消融季內(nèi)冰面融水動(dòng)態(tài)監(jiān)測(cè),從而進(jìn)一步提升對(duì)南極冰蓋典型消融區(qū)冰面融水變化規(guī)律的理解。

    猜你喜歡
    冰架融水冰蓋
    軍事文摘(2024年6期)2024-04-30 03:13:59
    南極冰架融化新原因確定
    我在融水唱貝江
    格陵蘭島的冰蓋悄悄融化
    參花(下)(2022年1期)2022-01-15 00:45:01
    融水糯米柚種植栽培管理與技術(shù)探究
    羅斯冰架探險(xiǎn)
    長距離輸水工程的冰期冰蓋數(shù)值模擬研究
    例舉高中數(shù)學(xué)解題切入點(diǎn)的找尋
    廣西苗鄉(xiāng)“蘆笙鬥馬節(jié)”驚險(xiǎn)刺激
    澳門月刊(2014年11期)2014-11-15 15:53:07
    衛(wèi)星數(shù)據(jù)更準(zhǔn)確預(yù)測(cè)冰蓋融化和海平面升高
    精品乱码久久久久久99久播| 在线播放国产精品三级| 亚洲成人久久爱视频| 日韩成人av中文字幕在线观看 | 久久精品国产99精品国产亚洲性色| 国产成人影院久久av| 观看美女的网站| 亚洲国产欧洲综合997久久,| 日韩精品有码人妻一区| 2021天堂中文幕一二区在线观| 久久鲁丝午夜福利片| av黄色大香蕉| 国产精品久久电影中文字幕| 99热6这里只有精品| 中国美女看黄片| 国产精品女同一区二区软件| 日韩三级伦理在线观看| 五月伊人婷婷丁香| 久久精品人妻少妇| av.在线天堂| 国产精品乱码一区二三区的特点| 日韩欧美国产在线观看| 欧美高清性xxxxhd video| 成人亚洲精品av一区二区| 男人狂女人下面高潮的视频| 国产中年淑女户外野战色| 久久久久九九精品影院| 亚洲最大成人av| 国产一区二区亚洲精品在线观看| 亚洲国产欧洲综合997久久,| 国产成年人精品一区二区| 日韩欧美三级三区| 日日撸夜夜添| 日本 av在线| 少妇被粗大猛烈的视频| 国产高清视频在线播放一区| 欧美日韩一区二区视频在线观看视频在线 | 欧美高清性xxxxhd video| 成熟少妇高潮喷水视频| 干丝袜人妻中文字幕| 久久99热6这里只有精品| 久久精品夜色国产| 亚洲国产色片| 99久国产av精品国产电影| 最近手机中文字幕大全| 一边摸一边抽搐一进一小说| 99久久无色码亚洲精品果冻| av卡一久久| 深爱激情五月婷婷| 中国国产av一级| 久久久精品94久久精品| 天堂√8在线中文| 精品久久久久久久久av| 亚洲熟妇中文字幕五十中出| 日日撸夜夜添| 十八禁网站免费在线| 色吧在线观看| 欧美bdsm另类| 搡老岳熟女国产| 欧美日韩精品成人综合77777| 一区二区三区高清视频在线| 亚洲一区二区三区色噜噜| 九九在线视频观看精品| 欧美最新免费一区二区三区| 亚洲成人av在线免费| 小蜜桃在线观看免费完整版高清| 最近最新中文字幕大全电影3| 99热这里只有精品一区| 老司机午夜福利在线观看视频| 亚洲一区高清亚洲精品| 干丝袜人妻中文字幕| 老司机午夜福利在线观看视频| 午夜福利成人在线免费观看| 国产91av在线免费观看| 成年av动漫网址| 99热这里只有精品一区| 九九爱精品视频在线观看| 校园春色视频在线观看| 一进一出抽搐动态| 亚洲国产色片| 精品无人区乱码1区二区| 99热6这里只有精品| 免费在线观看成人毛片| 熟女人妻精品中文字幕| 国产日本99.免费观看| 非洲黑人性xxxx精品又粗又长| 非洲黑人性xxxx精品又粗又长| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩一区二区精品| 久久午夜亚洲精品久久| 国产麻豆成人av免费视频| 国产精品久久久久久久久免| 成年av动漫网址| 国产av一区在线观看免费| 天堂影院成人在线观看| 十八禁国产超污无遮挡网站| 国产精品一区二区免费欧美| 国产高清三级在线| 日本一本二区三区精品| 亚洲在线自拍视频| 久久久a久久爽久久v久久| 99精品在免费线老司机午夜| av专区在线播放| 亚洲精华国产精华液的使用体验 | 99精品在免费线老司机午夜| 午夜日韩欧美国产| 欧美日韩国产亚洲二区| 国产毛片a区久久久久| 亚洲成人av在线免费| 两个人的视频大全免费| 午夜影院日韩av| 亚洲精品日韩av片在线观看| 午夜a级毛片| 免费电影在线观看免费观看| 亚洲av美国av| 哪里可以看免费的av片| 精品久久久久久久久av| a级毛片免费高清观看在线播放| 舔av片在线| 一级毛片久久久久久久久女| 美女xxoo啪啪120秒动态图| 一级毛片我不卡| 丝袜美腿在线中文| 日韩高清综合在线| 高清毛片免费观看视频网站| 久久国内精品自在自线图片| 亚洲欧美清纯卡通| 精品无人区乱码1区二区| 如何舔出高潮| 欧美高清性xxxxhd video| 99热这里只有是精品50| 一进一出好大好爽视频| 日本精品一区二区三区蜜桃| 久久久久久国产a免费观看| 国产成人91sexporn| 免费不卡的大黄色大毛片视频在线观看 | 国产成人a区在线观看| 亚洲在线观看片| 91狼人影院| av在线亚洲专区| 国产成人freesex在线 | 成人漫画全彩无遮挡| 欧美性感艳星| 男女视频在线观看网站免费| 午夜免费激情av| 在现免费观看毛片| 最近的中文字幕免费完整| 成人二区视频| 老女人水多毛片| 一级毛片aaaaaa免费看小| 欧美绝顶高潮抽搐喷水| 波野结衣二区三区在线| av在线蜜桃| 亚洲欧美日韩无卡精品| 日韩欧美三级三区| 激情 狠狠 欧美| 深爱激情五月婷婷| 久久午夜福利片| 亚洲欧美中文字幕日韩二区| 国产精品国产三级国产av玫瑰| av免费在线看不卡| 99热只有精品国产| 亚洲精品粉嫩美女一区| 观看美女的网站| 69人妻影院| 一个人看视频在线观看www免费| 国产精品av视频在线免费观看| 国产精品国产三级国产av玫瑰| 久久久精品94久久精品| 91久久精品国产一区二区三区| 一级毛片久久久久久久久女| 美女cb高潮喷水在线观看| 国产人妻一区二区三区在| 91精品国产九色| 啦啦啦观看免费观看视频高清| 欧美日韩国产亚洲二区| 精品一区二区三区人妻视频| 婷婷精品国产亚洲av在线| 亚洲美女视频黄频| 欧美一区二区国产精品久久精品| 成人午夜高清在线视频| 一本一本综合久久| 女的被弄到高潮叫床怎么办| 精品久久久噜噜| 中文资源天堂在线| 九九在线视频观看精品| 亚洲美女黄片视频| 99热网站在线观看| 亚洲无线在线观看| 日韩高清综合在线| 亚洲成a人片在线一区二区| 久久精品国产鲁丝片午夜精品| 看片在线看免费视频| 人人妻人人澡人人爽人人夜夜 | 啦啦啦啦在线视频资源| 国产综合懂色| 日日摸夜夜添夜夜添av毛片| 全区人妻精品视频| 深夜a级毛片| 99在线人妻在线中文字幕| 一夜夜www| 国产精品嫩草影院av在线观看| 三级毛片av免费| 日本在线视频免费播放| 搡老熟女国产l中国老女人| 国产精品亚洲一级av第二区| 日韩人妻高清精品专区| 草草在线视频免费看| 国产精品一区二区三区四区久久| 日韩欧美免费精品| 悠悠久久av| 亚洲美女黄片视频| 久久久久久久久大av| 国产国拍精品亚洲av在线观看| av卡一久久| 久久久久性生活片| 久久久久久久久久久丰满| 中国美白少妇内射xxxbb| av.在线天堂| 搡老岳熟女国产| 淫秽高清视频在线观看| 春色校园在线视频观看| 亚洲精品456在线播放app| 国产精品野战在线观看| 亚洲最大成人中文| 国产午夜精品久久久久久一区二区三区 | 午夜福利在线观看免费完整高清在 | 亚洲av五月六月丁香网| 美女cb高潮喷水在线观看| 人人妻人人澡欧美一区二区| 少妇熟女aⅴ在线视频| 赤兔流量卡办理| 中文字幕av在线有码专区| 午夜日韩欧美国产| 久久亚洲国产成人精品v| 久久久久国产精品人妻aⅴ院| 一级黄色大片毛片| 成人特级黄色片久久久久久久| 我的老师免费观看完整版| 国产高清视频在线播放一区| 亚洲自拍偷在线| 日韩制服骚丝袜av| 精品久久国产蜜桃| 观看美女的网站| 插阴视频在线观看视频| 精品一区二区三区视频在线| 一级av片app| 最近最新中文字幕大全电影3| 男女那种视频在线观看| 亚洲成av人片在线播放无| 国产成人freesex在线 | a级毛片a级免费在线| 日日摸夜夜添夜夜爱| 一本一本综合久久| 精品人妻熟女av久视频| 久久草成人影院| 成人美女网站在线观看视频| 亚洲国产日韩欧美精品在线观看| 亚洲专区国产一区二区| 又粗又爽又猛毛片免费看| 91av网一区二区| 亚洲18禁久久av| 日韩成人av中文字幕在线观看 | 三级经典国产精品| 国产爱豆传媒在线观看| 久久精品夜夜夜夜夜久久蜜豆| 有码 亚洲区| 欧美zozozo另类| 欧美xxxx性猛交bbbb| 内地一区二区视频在线| 国产伦在线观看视频一区| 最新中文字幕久久久久| 一本久久中文字幕| 国产精品电影一区二区三区| 嫩草影视91久久| 精品久久国产蜜桃| 久久韩国三级中文字幕| 亚洲第一电影网av| 色哟哟哟哟哟哟| 欧美一区二区国产精品久久精品| 99热全是精品| 天堂网av新在线| 老司机福利观看| 国产av一区在线观看免费| 国产亚洲精品av在线| 日韩中字成人| 99热网站在线观看| 99热这里只有是精品在线观看| 国产精品一区二区三区四区免费观看 | 91麻豆精品激情在线观看国产| 精品无人区乱码1区二区| 三级经典国产精品| 日本一本二区三区精品| 国产男靠女视频免费网站| 午夜老司机福利剧场| 欧美绝顶高潮抽搐喷水| 91久久精品国产一区二区三区| 级片在线观看| 成人三级黄色视频| 久久精品夜色国产| 欧美日韩一区二区视频在线观看视频在线 | 村上凉子中文字幕在线| 国产精品美女特级片免费视频播放器| 亚洲中文字幕日韩| 久久午夜福利片| 国产伦在线观看视频一区| 99在线人妻在线中文字幕| 久久精品国产亚洲av天美| 亚洲真实伦在线观看| 国产片特级美女逼逼视频| 国语自产精品视频在线第100页| 久久久精品欧美日韩精品| 日本一本二区三区精品| av免费在线看不卡| 国产精品永久免费网站| 一a级毛片在线观看| aaaaa片日本免费| 寂寞人妻少妇视频99o| 男女视频在线观看网站免费| 日韩国内少妇激情av| 久久久久久久午夜电影| 亚洲丝袜综合中文字幕| 欧美精品国产亚洲| 国产综合懂色| 男人舔女人下体高潮全视频| 中出人妻视频一区二区| 午夜影院日韩av| 亚洲欧美中文字幕日韩二区| 亚洲中文字幕一区二区三区有码在线看| 亚洲av第一区精品v没综合| 国产激情偷乱视频一区二区| 精品午夜福利在线看| 蜜桃久久精品国产亚洲av| 亚洲激情五月婷婷啪啪| 好男人在线观看高清免费视频| 国产黄片美女视频| 卡戴珊不雅视频在线播放| 中文资源天堂在线| 国产成人91sexporn| 超碰av人人做人人爽久久| 九九久久精品国产亚洲av麻豆| 97碰自拍视频| 亚洲精品亚洲一区二区| 亚洲精品在线观看二区| 久久婷婷人人爽人人干人人爱| 日日摸夜夜添夜夜爱| 女的被弄到高潮叫床怎么办| 国产精品99久久久久久久久| 国产精品综合久久久久久久免费| 在线观看一区二区三区| 一个人免费在线观看电影| 亚洲欧美中文字幕日韩二区| 欧美潮喷喷水| 老熟妇仑乱视频hdxx| 国产亚洲91精品色在线| 午夜免费激情av| 国产精品伦人一区二区| 国产伦在线观看视频一区| 亚洲真实伦在线观看| 激情 狠狠 欧美| 久久精品夜色国产| 国产色婷婷99| 国产精品国产高清国产av| 性欧美人与动物交配| 午夜激情福利司机影院| 国产v大片淫在线免费观看| 国产乱人偷精品视频| 国产成人一区二区在线| 三级国产精品欧美在线观看| 男女啪啪激烈高潮av片| 欧美激情国产日韩精品一区| 两个人的视频大全免费| 91麻豆精品激情在线观看国产| 天美传媒精品一区二区| 乱系列少妇在线播放| 亚洲精品久久国产高清桃花| 麻豆精品久久久久久蜜桃| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久久久久免费视频| 久久精品国产亚洲av涩爱 | 人妻少妇偷人精品九色| 99在线人妻在线中文字幕| 国产一区二区三区在线臀色熟女| 一级黄片播放器| 级片在线观看| 中文字幕av成人在线电影| 免费高清视频大片| 久久久久性生活片| 色5月婷婷丁香| 91午夜精品亚洲一区二区三区| 自拍偷自拍亚洲精品老妇| 看非洲黑人一级黄片| 亚洲欧美日韩卡通动漫| 美女高潮的动态| 亚洲人成网站在线观看播放| 精品熟女少妇av免费看| 在线观看一区二区三区| 99九九线精品视频在线观看视频| 三级毛片av免费| 国内久久婷婷六月综合欲色啪| 少妇熟女aⅴ在线视频| 在线观看av片永久免费下载| 无遮挡黄片免费观看| 国产精品精品国产色婷婷| 亚洲一区高清亚洲精品| 国产精品99久久久久久久久| 亚洲成人中文字幕在线播放| 综合色丁香网| 久久久久精品国产欧美久久久| 99热网站在线观看| 亚洲欧美日韩高清在线视频| 免费高清视频大片| 午夜日韩欧美国产| 少妇裸体淫交视频免费看高清| 国产精品亚洲美女久久久| 亚洲av一区综合| 搡老熟女国产l中国老女人| 欧美日韩乱码在线| 免费无遮挡裸体视频| 美女内射精品一级片tv| 亚洲中文字幕日韩| 国产精品一区www在线观看| 欧美+亚洲+日韩+国产| 久久99热6这里只有精品| 美女大奶头视频| 色av中文字幕| 九九久久精品国产亚洲av麻豆| 天堂影院成人在线观看| 波多野结衣高清作品| 可以在线观看的亚洲视频| 青春草视频在线免费观看| 国产视频一区二区在线看| 99热只有精品国产| 干丝袜人妻中文字幕| 黄色欧美视频在线观看| 麻豆久久精品国产亚洲av| 久久久精品94久久精品| 国产美女午夜福利| a级一级毛片免费在线观看| av在线观看视频网站免费| 日日干狠狠操夜夜爽| 亚洲国产精品久久男人天堂| 国产熟女欧美一区二区| 亚洲内射少妇av| 欧美国产日韩亚洲一区| 亚洲精品日韩av片在线观看| 床上黄色一级片| 成人午夜高清在线视频| 在线看三级毛片| 22中文网久久字幕| 美女cb高潮喷水在线观看| 夜夜看夜夜爽夜夜摸| 亚洲七黄色美女视频| 精品国内亚洲2022精品成人| 黄片wwwwww| 亚洲一区高清亚洲精品| 两个人的视频大全免费| 国产精品嫩草影院av在线观看| 国产高清不卡午夜福利| 天美传媒精品一区二区| 中文字幕av在线有码专区| 深夜精品福利| 观看免费一级毛片| 精品人妻偷拍中文字幕| 午夜激情欧美在线| 1000部很黄的大片| 精品人妻熟女av久视频| 欧美一区二区精品小视频在线| 色综合色国产| 男女下面进入的视频免费午夜| 麻豆精品久久久久久蜜桃| 日韩一本色道免费dvd| 精品人妻熟女av久视频| 全区人妻精品视频| 最近手机中文字幕大全| 麻豆一二三区av精品| 久久久久久国产a免费观看| 精品久久久久久久人妻蜜臀av| 麻豆国产97在线/欧美| 可以在线观看毛片的网站| 高清日韩中文字幕在线| 晚上一个人看的免费电影| 五月伊人婷婷丁香| 麻豆av噜噜一区二区三区| av中文乱码字幕在线| 三级毛片av免费| 日韩精品青青久久久久久| 久99久视频精品免费| 亚洲国产高清在线一区二区三| 99riav亚洲国产免费| 欧美区成人在线视频| 免费观看的影片在线观看| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久黄片| 亚洲天堂国产精品一区在线| 99热网站在线观看| 国产精品永久免费网站| 黄片wwwwww| 国产91av在线免费观看| 亚洲精品日韩av片在线观看| 欧美3d第一页| 成人三级黄色视频| 免费看av在线观看网站| 欧美一区二区精品小视频在线| 欧美成人一区二区免费高清观看| 欧美人与善性xxx| 久久久精品94久久精品| 日韩欧美在线乱码| 男女边吃奶边做爰视频| 十八禁国产超污无遮挡网站| 日韩制服骚丝袜av| 在线播放国产精品三级| 少妇人妻一区二区三区视频| 男人狂女人下面高潮的视频| 搡老岳熟女国产| 非洲黑人性xxxx精品又粗又长| 在线观看免费视频日本深夜| 搡老妇女老女人老熟妇| 搡老熟女国产l中国老女人| 日韩强制内射视频| 国产精品人妻久久久久久| 91久久精品国产一区二区三区| 久久久久久久久中文| 精品久久久久久久久av| 99精品在免费线老司机午夜| 国产一区二区在线av高清观看| 国内久久婷婷六月综合欲色啪| 亚洲精品成人久久久久久| 干丝袜人妻中文字幕| 2021天堂中文幕一二区在线观| 在线看三级毛片| 欧美又色又爽又黄视频| 亚洲国产色片| 日韩av不卡免费在线播放| 午夜免费激情av| 日本黄大片高清| 国产人妻一区二区三区在| 黄色一级大片看看| 国内精品美女久久久久久| 美女内射精品一级片tv| av福利片在线观看| 天天躁日日操中文字幕| 91在线观看av| 国产av不卡久久| 久久久a久久爽久久v久久| 国产精品,欧美在线| 成年av动漫网址| 亚洲欧美日韩无卡精品| 国产亚洲欧美98| 国产综合懂色| av天堂在线播放| 女人被狂操c到高潮| 亚洲精品成人久久久久久| 男人狂女人下面高潮的视频| 国产男靠女视频免费网站| 成人鲁丝片一二三区免费| 国产一区二区三区av在线 | 18禁在线无遮挡免费观看视频 | 国产高清有码在线观看视频| 搞女人的毛片| 我要看日韩黄色一级片| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 黄色欧美视频在线观看| 国产av不卡久久| 看黄色毛片网站| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 日日摸夜夜添夜夜爱| 小说图片视频综合网站| 国产成人影院久久av| 国产精品1区2区在线观看.| 91在线精品国自产拍蜜月| 亚洲国产精品合色在线| 美女xxoo啪啪120秒动态图| 日韩在线高清观看一区二区三区| 精品久久久久久久久av| 亚洲人成网站在线观看播放| 一夜夜www| 黄色配什么色好看| 国产亚洲精品久久久com| 成人亚洲精品av一区二区| 久久婷婷人人爽人人干人人爱| 又黄又爽又刺激的免费视频.| 长腿黑丝高跟| 国产精品不卡视频一区二区| 国产高清有码在线观看视频| 久久久久国产网址| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩东京热| 97热精品久久久久久| 久久久久久九九精品二区国产| 国产精品伦人一区二区| 高清日韩中文字幕在线| 国产伦精品一区二区三区四那| 午夜亚洲福利在线播放| 国产熟女欧美一区二区| 在线免费观看不下载黄p国产| 日韩欧美精品v在线| 99热这里只有精品一区| 日韩制服骚丝袜av| 中文字幕久久专区| 黄色视频,在线免费观看| 一个人看视频在线观看www免费| 国模一区二区三区四区视频| 国内少妇人妻偷人精品xxx网站| 欧美丝袜亚洲另类| 午夜福利18| 久久久久久久久久久丰满| 国产亚洲91精品色在线| 国产亚洲精品av在线| 久久韩国三级中文字幕| 午夜福利在线观看免费完整高清在 | 深夜精品福利| 色5月婷婷丁香|