• <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精品久久宅男| 一区二区三区乱码不卡18| av免费观看日本| 亚洲国产高清在线一区二区三| 内射极品少妇av片p| 久久这里有精品视频免费| 日韩视频在线欧美| 亚洲av综合色区一区| 国产精品嫩草影院av在线观看| 一级毛片久久久久久久久女| 久久精品久久久久久噜噜老黄| 少妇被粗大猛烈的视频| 蜜臀久久99精品久久宅男| 男的添女的下面高潮视频| 女性生殖器流出的白浆| 小蜜桃在线观看免费完整版高清| 熟妇人妻不卡中文字幕| 少妇人妻精品综合一区二区| 国产精品嫩草影院av在线观看| 日韩欧美精品免费久久| 视频中文字幕在线观看| 国产精品一区二区在线不卡| 精品久久久久久久久av| 午夜激情久久久久久久| 久久久成人免费电影| 日韩伦理黄色片| 亚洲天堂av无毛| 久久久亚洲精品成人影院| 人妻夜夜爽99麻豆av| 久久久久久久久久久免费av| 亚洲人成网站在线观看播放| 欧美成人午夜免费资源| 韩国av在线不卡| 大香蕉97超碰在线| 久久6这里有精品| 国产精品精品国产色婷婷| 搡老乐熟女国产| 亚洲精品乱久久久久久| 26uuu在线亚洲综合色| freevideosex欧美| 日韩人妻高清精品专区| 成人午夜精彩视频在线观看| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 久久99精品国语久久久| 欧美另类一区| 在线精品无人区一区二区三 | 联通29元200g的流量卡| 亚洲av在线观看美女高潮| 精品亚洲成国产av| 亚洲av在线观看美女高潮| 中文字幕免费在线视频6| 国产精品伦人一区二区| 麻豆精品久久久久久蜜桃| 国产精品99久久99久久久不卡 | 人人妻人人添人人爽欧美一区卜 | 嫩草影院新地址| 大片电影免费在线观看免费| 99热国产这里只有精品6| 中文欧美无线码| 精品久久久久久久久亚洲| av播播在线观看一区| 日本午夜av视频| 精品国产乱码久久久久久小说| 国产精品福利在线免费观看| 日本黄大片高清| 久久久久性生活片| 免费高清在线观看视频在线观看| 久久女婷五月综合色啪小说| 亚洲国产日韩一区二区| 日韩欧美一区视频在线观看 | 国产精品久久久久成人av| 国产成人免费无遮挡视频| 啦啦啦中文免费视频观看日本| 欧美日韩视频高清一区二区三区二| 久久人人爽av亚洲精品天堂 | 久久婷婷青草| 亚洲av成人精品一区久久| 国产欧美日韩精品一区二区| 欧美高清成人免费视频www| 久久久精品免费免费高清| av视频免费观看在线观看| 80岁老熟妇乱子伦牲交| 中文资源天堂在线| 国产探花极品一区二区| 在线观看免费视频网站a站| 老熟女久久久| 少妇精品久久久久久久| 观看免费一级毛片| 国产无遮挡羞羞视频在线观看| 狂野欧美白嫩少妇大欣赏| 少妇 在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲人与动物交配视频| 狂野欧美白嫩少妇大欣赏| 午夜激情福利司机影院| 麻豆精品久久久久久蜜桃| 亚洲精品国产av成人精品| 国产黄片视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 丰满迷人的少妇在线观看| 大片电影免费在线观看免费| 啦啦啦在线观看免费高清www| 国产亚洲欧美精品永久| 国产伦精品一区二区三区视频9| 国产一区二区三区av在线| 国产黄频视频在线观看| 国产精品成人在线| 在线观看av片永久免费下载| 国产乱人视频| 国产精品成人在线| 精品国产乱码久久久久久小说| 久久精品久久久久久噜噜老黄| 国内揄拍国产精品人妻在线| 26uuu在线亚洲综合色| 草草在线视频免费看| 一级毛片 在线播放| www.色视频.com| 日韩制服骚丝袜av| 麻豆国产97在线/欧美| 久久影院123| 精品午夜福利在线看| 在线观看三级黄色| 国产男人的电影天堂91| 色吧在线观看| 我的女老师完整版在线观看| 久久久成人免费电影| 日韩不卡一区二区三区视频在线| 男女边摸边吃奶| 欧美激情国产日韩精品一区| 久久久久久久久久成人| 国产精品久久久久久精品古装| 熟女人妻精品中文字幕| 国内精品宾馆在线| 成人无遮挡网站| 丰满少妇做爰视频| 久久国内精品自在自线图片| av卡一久久| 18禁裸乳无遮挡免费网站照片| 能在线免费看毛片的网站| 日韩欧美精品免费久久| 卡戴珊不雅视频在线播放| av在线播放精品| 嫩草影院入口| 插阴视频在线观看视频| 国产美女午夜福利| 男人爽女人下面视频在线观看| 最后的刺客免费高清国语| 国产精品久久久久成人av| 中文乱码字字幕精品一区二区三区| 欧美精品亚洲一区二区| 乱系列少妇在线播放| 亚洲欧美一区二区三区国产| 久久久久网色| 国产精品一二三区在线看| 亚洲欧美中文字幕日韩二区| 99热这里只有是精品在线观看| 国产视频内射| 久久久久久久国产电影| 乱码一卡2卡4卡精品| 少妇的逼水好多| 少妇高潮的动态图| 欧美xxⅹ黑人| 欧美日韩在线观看h| 性色avwww在线观看| 久久青草综合色| 久久国内精品自在自线图片| 热re99久久精品国产66热6| 久久精品国产鲁丝片午夜精品| 午夜免费男女啪啪视频观看| 久久鲁丝午夜福利片| 肉色欧美久久久久久久蜜桃| 肉色欧美久久久久久久蜜桃| 亚洲精品久久久久久婷婷小说| 99热国产这里只有精品6| xxx大片免费视频| 97在线视频观看| 嫩草影院入口| 99热6这里只有精品| 国产91av在线免费观看| 久热这里只有精品99| 久久女婷五月综合色啪小说| 身体一侧抽搐| 黄片无遮挡物在线观看| 99热网站在线观看| 色哟哟·www| 国产精品久久久久久av不卡| 韩国av在线不卡| 又黄又爽又刺激的免费视频.| 成人综合一区亚洲| 麻豆成人av视频| 你懂的网址亚洲精品在线观看| 久久久亚洲精品成人影院| 看非洲黑人一级黄片| 午夜激情久久久久久久| 欧美日韩视频高清一区二区三区二| 美女内射精品一级片tv| 国内精品宾馆在线| 最近手机中文字幕大全| 日韩 亚洲 欧美在线| 日日啪夜夜爽| 亚洲精华国产精华液的使用体验| 制服丝袜香蕉在线| 国产午夜精品久久久久久一区二区三区| 国产精品一及| 美女xxoo啪啪120秒动态图| 久久午夜福利片| 秋霞在线观看毛片| 亚洲国产精品国产精品| 91精品一卡2卡3卡4卡| 亚洲最大成人中文| 三级国产精品欧美在线观看| 国产精品福利在线免费观看| 久久久久视频综合| 啦啦啦中文免费视频观看日本| 亚洲欧美清纯卡通| 建设人人有责人人尽责人人享有的 | 成人国产av品久久久| 久久精品夜色国产| 在线免费十八禁| 亚洲国产成人一精品久久久| 国产色婷婷99| 少妇的逼水好多| 亚洲av二区三区四区| 亚洲精品自拍成人| 26uuu在线亚洲综合色| 97在线视频观看| 毛片女人毛片| 国产精品久久久久久久久免| 国产色婷婷99| 亚洲av二区三区四区| 欧美亚洲 丝袜 人妻 在线| 亚洲自偷自拍三级| 青春草亚洲视频在线观看| 啦啦啦中文免费视频观看日本| 亚洲欧美清纯卡通| 国产精品一二三区在线看| 久久久久久久亚洲中文字幕| 高清视频免费观看一区二区| 亚洲av日韩在线播放| 又大又黄又爽视频免费| 国产综合精华液| 国产亚洲一区二区精品| 人妻 亚洲 视频| 国产熟女欧美一区二区| 亚洲国产高清在线一区二区三| 少妇人妻一区二区三区视频| 亚洲精品色激情综合| 王馨瑶露胸无遮挡在线观看| 一级av片app| 亚洲欧美成人精品一区二区| 在线看a的网站| 草草在线视频免费看| 欧美日韩视频高清一区二区三区二| 人妻少妇偷人精品九色| 又爽又黄a免费视频| 国产又色又爽无遮挡免| 成年美女黄网站色视频大全免费 | 简卡轻食公司| 直男gayav资源| 精品少妇黑人巨大在线播放| 午夜福利网站1000一区二区三区| 七月丁香在线播放| 99久久综合免费| 色视频www国产| 亚洲精华国产精华液的使用体验| 我的老师免费观看完整版| 又大又黄又爽视频免费| 久久青草综合色| 成年免费大片在线观看| 日韩一区二区视频免费看| 高清不卡的av网站| 嫩草影院入口| 国产v大片淫在线免费观看| 久久这里有精品视频免费| 国产高潮美女av| 少妇人妻久久综合中文| 高清不卡的av网站| 久久精品久久久久久久性| 一级爰片在线观看| 亚洲经典国产精华液单| 国内揄拍国产精品人妻在线| 久久青草综合色| 热re99久久精品国产66热6| 少妇的逼水好多| 男女啪啪激烈高潮av片| 九色成人免费人妻av| 日本一二三区视频观看| 欧美高清性xxxxhd video| 国产又色又爽无遮挡免| 交换朋友夫妻互换小说| 爱豆传媒免费全集在线观看| 成人无遮挡网站| 亚洲不卡免费看| 三级经典国产精品| 最新中文字幕久久久久| 91精品伊人久久大香线蕉| 18禁动态无遮挡网站| 99国产精品免费福利视频| 2022亚洲国产成人精品| 狂野欧美激情性xxxx在线观看| 亚洲四区av| 老熟女久久久| 久久国产精品大桥未久av | 街头女战士在线观看网站| 久久影院123| 色视频www国产| 久久婷婷青草| 日韩一本色道免费dvd| 亚洲内射少妇av| 水蜜桃什么品种好| 在线播放无遮挡| 黄片无遮挡物在线观看| 在线观看美女被高潮喷水网站| 妹子高潮喷水视频| 久久久久国产网址| 精品人妻视频免费看| 国产黄片视频在线免费观看| 国产精品伦人一区二区| 伊人久久精品亚洲午夜| 哪个播放器可以免费观看大片| 校园人妻丝袜中文字幕| 国产v大片淫在线免费观看| 国产探花极品一区二区| 九九久久精品国产亚洲av麻豆| 欧美成人a在线观看| 午夜精品国产一区二区电影| h日本视频在线播放| 搡女人真爽免费视频火全软件| 超碰97精品在线观看| 国产一区有黄有色的免费视频| 免费av不卡在线播放| 99热国产这里只有精品6| 国产成人免费无遮挡视频| 久久久久久人妻| 国产欧美亚洲国产| 男女下面进入的视频免费午夜| 国产成人a∨麻豆精品| 国产精品伦人一区二区| 久久鲁丝午夜福利片| 日韩不卡一区二区三区视频在线| 王馨瑶露胸无遮挡在线观看| 尤物成人国产欧美一区二区三区| 性高湖久久久久久久久免费观看| 婷婷色综合大香蕉| 天天躁日日操中文字幕| 亚洲精品乱码久久久v下载方式| 日本黄色日本黄色录像| 日本午夜av视频| av国产久精品久网站免费入址| 老师上课跳d突然被开到最大视频| 久久久久久久久久人人人人人人| 精品国产乱码久久久久久小说| 毛片女人毛片| 久久精品国产鲁丝片午夜精品| 性高湖久久久久久久久免费观看| 亚洲伊人久久精品综合| 美女cb高潮喷水在线观看| 极品少妇高潮喷水抽搐| 高清午夜精品一区二区三区| 中文精品一卡2卡3卡4更新| 日韩亚洲欧美综合| 黄色欧美视频在线观看| 亚洲av欧美aⅴ国产| 91久久精品国产一区二区三区| 国产视频内射| 成人高潮视频无遮挡免费网站| 草草在线视频免费看| 亚洲精品国产av成人精品| 国产淫语在线视频| 日韩av不卡免费在线播放| 亚洲国产精品一区三区| 国产中年淑女户外野战色| 亚洲久久久国产精品| 日韩不卡一区二区三区视频在线| 一区二区三区免费毛片| 久久精品国产a三级三级三级| 欧美+日韩+精品| 在线观看一区二区三区| 大片免费播放器 马上看| 久久精品国产亚洲av天美| 80岁老熟妇乱子伦牲交| 99re6热这里在线精品视频| 免费看光身美女| 人人妻人人看人人澡| 丝瓜视频免费看黄片| 久久精品久久久久久久性| 97热精品久久久久久| 建设人人有责人人尽责人人享有的 | 春色校园在线视频观看| 伊人久久国产一区二区| 人人妻人人澡人人爽人人夜夜| 美女国产视频在线观看| 美女cb高潮喷水在线观看| 六月丁香七月| 亚洲人成网站在线播| 一个人看的www免费观看视频| www.色视频.com| 亚洲va在线va天堂va国产| 国产免费又黄又爽又色| 交换朋友夫妻互换小说| 一级毛片久久久久久久久女| 亚洲人成网站在线播| 亚洲经典国产精华液单| 亚洲一级一片aⅴ在线观看| 男男h啪啪无遮挡| av在线app专区| 国产免费福利视频在线观看| 亚州av有码| 麻豆成人午夜福利视频| 久久婷婷青草| 亚洲激情五月婷婷啪啪| 中国国产av一级| 最后的刺客免费高清国语| 国产91av在线免费观看| 久久毛片免费看一区二区三区| 永久网站在线| 午夜激情久久久久久久| 国产中年淑女户外野战色| 日韩国内少妇激情av| 国产精品秋霞免费鲁丝片| 国产精品嫩草影院av在线观看| 亚洲人成网站在线观看播放| 欧美日韩视频精品一区| 国模一区二区三区四区视频| 伊人久久国产一区二区| 中文字幕av成人在线电影| 免费大片18禁| 国产乱人偷精品视频| 交换朋友夫妻互换小说| 我要看日韩黄色一级片| 香蕉精品网在线| 成年美女黄网站色视频大全免费 | 人妻 亚洲 视频| 菩萨蛮人人尽说江南好唐韦庄| 女性被躁到高潮视频| 伦精品一区二区三区| 国产精品国产三级国产av玫瑰| 超碰97精品在线观看| 国产免费一级a男人的天堂| 日本爱情动作片www.在线观看| 国产高清三级在线| 在线观看免费日韩欧美大片 | 欧美xxⅹ黑人| 久久综合国产亚洲精品| 亚洲人成网站高清观看| 久久国产乱子免费精品| 亚洲国产色片| 人人妻人人澡人人爽人人夜夜| 欧美国产精品一级二级三级 | 狠狠精品人妻久久久久久综合| 边亲边吃奶的免费视频| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| 99久久精品热视频| av国产免费在线观看| 久久精品久久精品一区二区三区| 日韩电影二区| 国产精品秋霞免费鲁丝片| 在线免费观看不下载黄p国产| 国产高清不卡午夜福利| 美女xxoo啪啪120秒动态图| 中文字幕人妻熟人妻熟丝袜美| 街头女战士在线观看网站| 男人爽女人下面视频在线观看| 看十八女毛片水多多多| 一级毛片电影观看| 尤物成人国产欧美一区二区三区| 日产精品乱码卡一卡2卡三| 国产精品熟女久久久久浪| 大香蕉久久网| 男女国产视频网站| 尤物成人国产欧美一区二区三区| 久久99蜜桃精品久久| 日本欧美国产在线视频| 人人妻人人澡人人爽人人夜夜| 观看美女的网站| 蜜桃在线观看..| 日韩av免费高清视频| 久热久热在线精品观看| 国产高潮美女av| 欧美日韩国产mv在线观看视频 | 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品乱码久久久v下载方式| 国产亚洲5aaaaa淫片| 毛片女人毛片| 国产精品一区二区在线观看99| 青春草亚洲视频在线观看| 亚洲av成人精品一区久久| 3wmmmm亚洲av在线观看| 国产乱来视频区| 高清午夜精品一区二区三区| 欧美zozozo另类| 国产精品不卡视频一区二区| av又黄又爽大尺度在线免费看| 久久久久久久国产电影| 久久99蜜桃精品久久| 国产成人aa在线观看| 精品熟女少妇av免费看| av卡一久久| 欧美人与善性xxx| 激情五月婷婷亚洲| 3wmmmm亚洲av在线观看| 午夜福利影视在线免费观看| 看免费成人av毛片| 我要看日韩黄色一级片| 成人美女网站在线观看视频| 欧美bdsm另类| 久久久久性生活片| 永久免费av网站大全| 午夜老司机福利剧场| 欧美激情国产日韩精品一区| 亚洲国产欧美人成| 亚洲人成网站高清观看| 一区二区av电影网| 如何舔出高潮| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美中文字幕日韩二区| 精品国产一区二区三区久久久樱花 | 亚洲国产精品专区欧美| 国产一级毛片在线| 国产爱豆传媒在线观看| 中文字幕精品免费在线观看视频 | 免费大片黄手机在线观看| 成人漫画全彩无遮挡| 亚洲欧美中文字幕日韩二区| 波野结衣二区三区在线| 久久国产亚洲av麻豆专区| 久久精品国产亚洲av涩爱| 亚洲国产精品一区三区| 国产精品国产av在线观看| 九色成人免费人妻av| 观看免费一级毛片| 又大又黄又爽视频免费| 极品教师在线视频| 国产精品国产三级专区第一集| 国产 精品1| 午夜福利在线观看免费完整高清在| 午夜老司机福利剧场| 国产白丝娇喘喷水9色精品| 黄色怎么调成土黄色| 国产精品不卡视频一区二区| 多毛熟女@视频| 秋霞伦理黄片| 少妇的逼好多水| 联通29元200g的流量卡| 国产亚洲91精品色在线| 97在线人人人人妻| 99热网站在线观看| 高清在线视频一区二区三区| 色婷婷av一区二区三区视频| 亚洲av在线观看美女高潮| 精品国产三级普通话版| 晚上一个人看的免费电影| 中文欧美无线码| 超碰97精品在线观看| 亚洲熟女精品中文字幕| 人妻系列 视频| 国产亚洲5aaaaa淫片| 人妻系列 视频| 免费人妻精品一区二区三区视频| a级毛色黄片| 国产黄频视频在线观看| 天堂俺去俺来也www色官网| 国产亚洲5aaaaa淫片| a级毛色黄片| 一区二区三区精品91| 性色av一级| 一区在线观看完整版| 女的被弄到高潮叫床怎么办| 亚洲天堂av无毛| 亚洲av不卡在线观看| av免费观看日本| 一级av片app| 又大又黄又爽视频免费| 中文字幕制服av| 国产精品av视频在线免费观看| av黄色大香蕉| 久久99热这里只有精品18| 日日啪夜夜爽| 国产精品国产三级国产av玫瑰| 又黄又爽又刺激的免费视频.| 一本色道久久久久久精品综合| 国产日韩欧美亚洲二区| 国产精品国产三级国产av玫瑰| 一级爰片在线观看| 亚洲av福利一区| 中文字幕久久专区| 黑人猛操日本美女一级片| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 91狼人影院| av国产精品久久久久影院| 五月开心婷婷网| 久久 成人 亚洲| 夫妻性生交免费视频一级片| 国产在视频线精品| 成年美女黄网站色视频大全免费 | 色网站视频免费| 能在线免费看毛片的网站| 国产精品久久久久久久电影| 亚洲精品国产av成人精品| 国产女主播在线喷水免费视频网站| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 麻豆精品久久久久久蜜桃| 成人黄色视频免费在线看| 国产爽快片一区二区三区| 一级毛片 在线播放| 夜夜爽夜夜爽视频| 精品久久国产蜜桃| 小蜜桃在线观看免费完整版高清| 丰满迷人的少妇在线观看| 国产探花极品一区二区| 久久韩国三级中文字幕|