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

    1993—2016年喀喇昆侖山中部Shigar流域冰川物質(zhì)平衡變化空間特征研究

    2021-04-07 16:15:44王璐王寧練李志杰陳安安夏瑋靜
    冰川凍土 2021年1期
    關(guān)鍵詞:冰川海拔高程

    王璐,王寧練,3,李志杰,陳安安,夏瑋靜

    (1.陜西省地表系統(tǒng)與環(huán)境承載力重點實驗室,陜西西安710127;2.西北大學城市與環(huán)境學院地表系統(tǒng)與災(zāi)害研究院,陜西西安710127;3.中國科學院青藏高原地球科學卓越創(chuàng)新中心,北京100101)

    0 引言

    山地冰川是冰凍圈的重要組成部分[1],是氣候變化的靈敏指示器之一[2-3],同時也是重要的淡水資源,被稱為“高山固體水庫”[4]。最新研究結(jié)果表明高亞洲大部分流域冰川融水對河流補給量將在本世紀中葉達到峰值,到本世紀末,冰川融水對河流徑流的補給率將有明顯的下降趨勢[5-8]。因此,冰川變化的長期監(jiān)測和研究在水資源保護和利用、區(qū)域可持續(xù)發(fā)展、災(zāi)害預(yù)警等方面具有重大意義。

    受氣候變暖的影響,全球大多數(shù)山地冰川在過去100年間處于退縮狀態(tài),在近幾十年內(nèi)呈加速退縮趨勢[9-11]。高亞洲地區(qū)冰川亦普遍萎縮,并呈現(xiàn)出自東南向西北萎縮速率逐漸變小的空間格局[12]。Hewitt[13]發(fā)現(xiàn)從20世紀90年代后期開始喀喇昆侖山中部冰川出現(xiàn)了停滯和前進現(xiàn)象,提出了“喀喇昆侖異?!保鹆藝鴥?nèi)外眾多學者的關(guān)注和討論[12-17]。Farinotti等[18]指出喀喇昆侖的兩個顯著異常特征是微弱的正物質(zhì)平衡和大量的躍動冰川存在。Gardelled等[19]基于SRTM DEM和SOPT5立體像對利用大地測量學方法估算得出在喀喇昆侖山中部1999—2011年間冰川呈正物質(zhì)平衡,為(+0.10±0.16)m w.e.·a-1;Ka¨a¨b等[20]基于SRTM DEM和ICESat數(shù)據(jù)估算得出喀喇昆侖山2003—2008年冰川高程變化:冬季為(+0.41±0.04)m·a-1,秋季為(-0.07±0.04)m·a-1。有研究[21-22]指出近二十年來喀喇昆侖山冰川末端前進、退縮和穩(wěn)定等現(xiàn)象同時存在,大量冰川正在發(fā)生躍動或曾出現(xiàn)過躍動。因此,喀喇昆侖山冰川變化存在明顯的空間差異性,其冰川變化的異?,F(xiàn)象主要表現(xiàn)為冰川的前進和躍動。但目前來說,對喀喇昆侖山不同類型冰川的研究,尤其是對不同運動狀態(tài)冰川的表面高程變化和物質(zhì)平衡的研究較少,這種異?,F(xiàn)象值得進一步關(guān)注。

    冰川變化是由氣候、地形、冰川自身規(guī)模等多種因素共同作用導致的。一般情況下,在空間區(qū)域較小的范圍內(nèi),氣候變化的特征通常趨向于一致化[23]。因此對于空間范圍較小的區(qū)域而言,氣候變化是導致冰川變化總體特征的主要因素[24],而地形特征和冰川自身的規(guī)模特征等則是導致冰川變化差異性的主要因素。本研究選取喀喇昆侖山中部冰川變化差異較大的Shigar流域為研究區(qū)域,以探究同一區(qū)域內(nèi)同時存在的不同類型冰川的面積變化、表面高程變化、物質(zhì)平衡變化及其可能的原因。

    1 研究區(qū)概況

    Shigar流域位于印度河上游,喀喇昆侖山中部,介于74°53′~76°45′E和35°19′~36°07′N之間,流域面積約6.92×103km2。流域內(nèi)地形高聳,海拔分布在2 150~8 611 m之間,45%左右的地區(qū)海拔高于5 000 m[25]。Shigar流域受高聳地形的影響,氣溫低且溫差大[26],流域降水主要受西風氣流的影響,陡峭的地形導致氣流的頻繁擾動,使得山地降水相對豐富。流域內(nèi)山地冰川廣泛發(fā)育而且規(guī)模較大,其中巴爾托洛冰川和Biafo冰川的面積分別高達809.11 km2和559.81 km2。

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)來源

    本研究采用的遙感影像數(shù)據(jù)(表1)來源于的Landsat系列衛(wèi)星搭載的TM/ETM+/OLI傳感器,影像均下載于美國地質(zhì)調(diào)查局網(wǎng)站(USGS,http://www.usgs.gov)。本文盡可能選擇消融季末期、云雪覆蓋較少的影像,最終選擇了1990—2016年的影像,其中以1993年、2000年、2016年3個時期的影像共計12景提取冰川邊界,選擇共計14景影像判定躍動冰川的運動特征。所有影像都由USGS進行了輻射校正和地面控制點校正以及地形校正。TM影像數(shù)據(jù)分辨率為30 m,ETM+/OLI影像數(shù)據(jù)采用Gram-Schmidt光譜銳化方法[27-28]進行多光譜與全色波段融合后分辨率為15 m。

    圖1 Shigar流域地理位置及冰川分布Fig.1 Map showing the geographic location of Shigar basin and distribution of glaciers

    本研究中的DEM數(shù)據(jù)采用了SRTM 3 DEM(http://eros.usgs.gov)和TanDEM-X DEM(https://download.geoservice.dlr.de/TDM90/),兩者的空間分辨率都為90 m。由于SRTM DEM是基于EGM96大地水準高的正常高,而TanDEM-X DEM是基于WGS84橢球體的大地高,因此需要通過hwgs84=HEGM96+NEGM96geoid將TanDEM-X DEM從WGS84橢球高轉(zhuǎn)換為EGM96大地水準高。此外,在目視解譯冰川邊界過程中,還參考了Randolph Glacier Inventory 6.0編 目 數(shù) 據(jù)(http://www.glims.org)和Google Earth高分辨率影像。

    2.2 冰川邊界提取

    2.2.1 冰川邊界提取方法

    基于遙感影像提取冰川邊界的方法主要有人工目視解譯方法和監(jiān)督分類、非監(jiān)督分類、比值閾值法、積雪指數(shù)閾值法等自動分類方法[29]。Shigar流域所在的喀喇昆侖山中部地區(qū)冰川末端有大量表磧覆蓋,而當前尚無普遍適用的表磧型冰川邊界自動化提取方法[30],并且計算機自動解譯的結(jié)果通常需要進行大量的人工修正。因此,本文采用目前精度較高的目視解譯方法提取冰川邊界。

    為了減小冰川邊界提取中的誤差,首先利用ENVI軟件對ETM+/OLI影像數(shù)據(jù)進行多光譜與全色波段融合,然后利用ArcGIS 10.4軟件對影像進行假彩色合成,對于裸冰區(qū)而言,直接根據(jù)冰川區(qū)與非冰川區(qū)顯著的色彩差異進行冰川邊界的數(shù)字化,精度可以達到像元級。對于表磧覆蓋區(qū)的識別主要是參考冰川編目數(shù)據(jù)[31],依據(jù)冰川末端的水文特征和兩側(cè)地形條件:一方面是由于冰川的長期融水及融水的沖刷,使得冰川末端形成較明顯的水系出露,這是識別有表磧覆蓋的冰川末端的標志之一;另一方面是冰舌部位表磧覆蓋區(qū)與冰川兩側(cè)地形相比有較大的高程差,可利用DEM數(shù)據(jù)來獲取等高線用于判斷冰川末端表磧覆蓋。相鄰冰川以山脊線為界劃分。在解譯完成后,將目視解譯的冰川邊界與RGI 6.0冰川邊界對比,可看出目視解譯的冰川邊界準確度較高,可有效避免將積雪誤判為冰川[圖2(a)]。同時,利用Google Earth影像對冰川邊界進行檢驗,發(fā)現(xiàn)目視解譯的冰川邊界能將冰川末端的表磧覆蓋區(qū)準確的包括在內(nèi)[圖2(b)]??紤]到遙感影像的分辨率狀況,本文只統(tǒng)計面積大于0.01 km2的冰川[32]。

    表1 遙感影像數(shù)據(jù)Table 1 Remote sensing images

    利用Landsat衛(wèi)星遙感影像(表1)來確定1993年和2016年的冰川末端位置變化情況,用于判斷1993—2016年間冰川是否發(fā)生了前進或退縮,還是保持穩(wěn)定。參考Rankl等[16]提出的定義:當冰川末端至少變化60 m時,該冰川才能被判斷為前進或者退縮。本研究根據(jù)該定義來進行冰川分類:當冰川末端前進超過60 m時為前進冰川;冰川末端退縮超過60 m時為退縮冰川;否則為穩(wěn)定冰川。冰川末端變化利用冰川主流線長度變化來確定,冰川主流線通過自動識別方法和結(jié)合后期人工修正來提取。自動提取方法原理如下:首先利用DEM數(shù)據(jù)識別冰川最低點和局部最高點,然后通過冰川邊界點構(gòu)建泰森多邊形,基于泰森多邊形頂點特性選擇和優(yōu)化路徑,得出冰川主流線,最后選取冰川主流線中最長路徑即為冰川長度。躍動冰川是指周期性發(fā)生快速運動的冰川,當冰川躍動時,冰川會在短時間內(nèi)出現(xiàn)快速前進,其速度往往是平時的數(shù)倍至上百倍,冰川躍動后,可能會出現(xiàn)數(shù)年冰川補償性后退[33]。本研究基于前人的研究成果[21,34-38]來確定躍動冰川,同時結(jié)合冰川末端運動速度、冰川表面高程變化來判斷冰川是否在研究時段內(nèi)發(fā)生過躍動。

    圖2 目視解譯冰川邊界Fig.2 Glacier boundary based on manual visual interpretation:comparison with Randolph 6.0 Glacier inventory(a);the high-resolution Google Earth imagery was used to verify glacier boundary(b)

    2.2.2 冰川邊界精度評價

    冰川邊界人工提取的精度主要受人為判讀誤差、遙感影像質(zhì)量和空間分辨率等因素的影響[39]。在目視解譯工作開始前,工作人員均通過專家培訓系統(tǒng)地學習了冰川學的相關(guān)知識,盡可能減少人為誤判帶來的誤差。有研究對USGS所提供的正射矯正后的Landsat遙感影像進行正射精度檢驗,發(fā)現(xiàn)其檢驗的大部分影像矯正精度在1/2個像元左右,因此USGS提供的Landsat遙感影像具有相當高的校正精度[40],因此本文未對影像進行進一步矯正,僅考慮遙感影像空間分辨率造成的誤差,采用下式計算[41]:

    式中:β為單條冰川面積的不確定性;N為冰川邊界所經(jīng)過的像元數(shù)量;λ為空間分辨率。計算得到1993年、2000年、2016年冰川面積提取的誤差值分別為±160.73 km2、±80.87 km2、±80.69 km2。

    計算冰川面積變化的不確定時,考慮到山地冰川面積變化主要發(fā)生在冰川中下部[42],因此計算誤差時僅涉及冰川末端發(fā)生變化的部分,采用下式計算[43]:

    式中:βAB為兩期冰川面積變化誤差,βA和βB分別為單期冰川面積提取誤差。

    2.3 冰川物質(zhì)平衡計算方法

    2.3.1 DEM配準及校正

    (1)空 間 配 準:SRTM DEM和TanDEM-X DEM的獲取時間分別為2000年和2013年,在這13年內(nèi),非冰川區(qū)的地形是基本保持穩(wěn)定的。因此基于非冰川區(qū)DEM的高程差進行多源DEM間的空間配準。不同DEM之間空間匹配誤差導致的高程差dh與坡度α、坡向φ之間存在明顯的統(tǒng)計關(guān)系,因此可以利用該關(guān)系對偏差進行校正[44]:

    (2)高程殘差校正:多源DEM在經(jīng)過空間匹配后仍會存在一定的高程殘差,高程殘差的存在主要與DEM的空間分辨率有關(guān)。多源DEM之間的高程殘差和地表最大曲率之間存在統(tǒng)計學關(guān)系,這種關(guān)系在冰川區(qū)和非冰川區(qū)表現(xiàn)出良好的一致性[45],因此可利用非冰川區(qū)高程殘差和地表最大曲率之間的關(guān)系對冰川區(qū)的高程殘差進行校正[46]。

    (3)SRTM DEM透射深度校正:由于SRTM C波段數(shù)據(jù)對冰雪有一定的穿透能力,而X波段數(shù)據(jù)對冰雪的穿透可忽略不計。因此本研究利用SRTM X波段數(shù)據(jù)對C波段數(shù)據(jù)的冰雪透射深度進行估算和校 正[22,45,47]。校 正后 的 結(jié) 果 顯示 在Shigar流 域SRTM C波段對雪冰的平均透射深度約為0.85 m。

    (4)剔除異常值:為更真實的反映冰川表面高程變化,根據(jù)冰川表面高程變化圖以及高程變化的柵格值和柵格數(shù)量直方圖得出有95%以上的柵格集中±150 m范圍內(nèi)。因此本文設(shè)定±150 m為閾值范圍來剔除不同DEM間高程差的異常值[19]。

    2.3.2 冰川物質(zhì)平衡估算采用以下公式計算[48]:

    式中:B表示物質(zhì)平衡,ρ為冰川體積向物質(zhì)平衡轉(zhuǎn)換的密度,本文采用Huss推薦的平均密度850 kg·m-3作為轉(zhuǎn)換密度[49]。Sg表示冰川面積,n為冰川區(qū)的像元總數(shù),Δhi表示研究時段內(nèi)像元點的高程差異,Sp表示單個像元的面積。

    2.3.3 誤差評估

    利用平均標準誤差來評估冰川表明高程變化的不確定性[50]:

    式中:σ為不同DEM數(shù)據(jù)間高程變化的相對誤差;SE表示平均標準誤差;N為空間去相關(guān)處理后的像元個數(shù),空間去相關(guān)距離為600 m[51];STDVnoglacier表示非冰川區(qū)高程差的標準差;MED表示非冰川區(qū)的平均高程差。

    研究時段的物質(zhì)平衡不確定性UM用下列公式計算:

    式中:t表示研究時段長度(2000—2013年);Δh代表冰川區(qū)的平均高程差;Δρ為冰密度的不確定性(60 kg·m-3);ρW和ρI分別表示水密度(1 000 kg·m-3)和轉(zhuǎn)換密度(850 kg·m-3)。

    3 結(jié)果與分析

    3.1 冰川面積變化

    2016年Shigar流域內(nèi)分布有冰川516條,總面積為(2 895.30±80.69)km2,約占流域面積的42%。1993—2016年間Shigar流域冰川面積由(2 897.97±160.73)km2減少到(2 895.30±80.69)km2,年均萎縮率為(-0.00±0.02)%·a-1(表2)。1993—2000年和2000—2016年兩個研究時段冰川面積變化率分別為(+0.01±0.02)%·a-1和(-0.01±0.02)%·a-1,表明近20年來Shigar流域冰川面積變化微弱。同時,由于冰川末端前進導致冰川匯合,研究區(qū)內(nèi)冰川數(shù)量由1993年的519條減少為516條。總體看來,1993—2016年間Shigar流域冰川基本處于穩(wěn)定狀態(tài),面積變化非常微弱;部分冰川還存在擴張現(xiàn)象,主要表現(xiàn)為冰川末端的前進和躍動。

    表2 1993—2016年Shigar流域冰川面積變化Table 2 Changes in glacier area of the Shigar basin during 1993—2016

    基于1993—2016年流域內(nèi)冰川末端變化和冰川表面高程變化情況,本研究將冰川分為躍動冰川、前進冰川、穩(wěn)定冰川和退縮冰川進行討論。1993—2016年間,有25條冰川發(fā)生躍動[21,34-38],68條冰川末端發(fā)生前進,50條冰川末端發(fā)生退縮,376條冰川末端保持穩(wěn)定。其中,在68條前進冰川中,有41條在1993—2000年發(fā)生過前進,23條在2000—2016年間發(fā)生過前進,4條冰川在1993—2016年間持續(xù)前進。在此期間,冰川躍動導致流域冰川面積增加1.30 km2,冰川末端前進導致流域冰川面積增加0.86 km2,冰川末端退縮導致流域冰川面積減少3.48 km2,穩(wěn)定冰川面積減少了1.34 km2。

    3.1.1 冰川規(guī)模特征

    躍動冰川中冰川長度最大值為12.58 km(平均值為6.30 km),前進冰川中冰川長度最大值為12.30 km(平均值為3.34 km),退縮冰川中冰川長度最大值為16.80 km(平均值為3.08 km),穩(wěn)定冰川中冰川長度最大值為62.97 km(平均值為1.87 km)。從不同長度區(qū)間內(nèi)各類型冰川數(shù)量占該類冰川總數(shù)量的百分比[圖3(a)]可看出,躍動冰川的冰川長度相比其他3類冰川要集中在較長范圍,60%的躍動冰川長度>5 km;前進冰川比退縮冰川的冰川長度要略長,但這兩類冰川均有超過70%的冰川長度在1~5 km內(nèi);穩(wěn)定冰川長度多集中在<1 km的范圍內(nèi),但流域內(nèi)長度>40 km的冰川共有4條,其中3條為穩(wěn)定冰川。

    圖3 Shigar流域內(nèi)躍動冰川、前進冰川、退縮冰川和穩(wěn)定冰川在不同長度規(guī)模(a)和面積規(guī)模(b)等級內(nèi)的數(shù)量占比Fig.3 Quantity percentage of surge-type,advancing,retreating,and stable glaciers in different glacier length(a)and area(b)in the Shigar basin

    從不同規(guī)模內(nèi)各類型冰川數(shù)量占該類冰川總數(shù)量的百分比[圖3(b)]可看出,躍動冰川多為較大規(guī)模,60%的躍動冰川面積>5 km2,這與冰川長度分布規(guī)律一致;74%的前進冰川分布在規(guī)模0.5~5 km2范圍內(nèi);有50%的退縮冰川面積<1 km2,退縮冰川多為小規(guī)模;穩(wěn)定冰川規(guī)模差異較大,流域內(nèi)規(guī)模>50 km2的6條大冰川中有5條均為穩(wěn)定冰川,此外,穩(wěn)定冰川中有75%的冰川規(guī)模<1 km2。結(jié)合不同規(guī)模中各類冰川數(shù)量占比可看出,在流域內(nèi)氣候條件變化一致的情況下,大規(guī)模冰川對氣候變化的響應(yīng)較小規(guī)模冰川要滯后,更容易保持穩(wěn)定狀態(tài)。

    3.1.2 冰川地形特征

    根據(jù)流域內(nèi)的山谷位置將流域分為四部分,分別是Chogo Lungma區(qū)域(Ch)、Biafo區(qū)域(Bi)、Panmah區(qū)域(Pa)和巴爾托洛區(qū)域(Ba)。從不同類型冰川的空間分布(圖4)可看出,躍動冰川主要分布在Pa區(qū)域和Ba區(qū)域南部山坡的北朝向地區(qū),躍動冰川與前進冰川分布地區(qū)具有一致性;此外,前進冰川還分布在Ch地區(qū);退縮冰川主要分布在Ch區(qū)域和Shigar河干流東部山坡的東朝向地區(qū)。同時還發(fā)現(xiàn)流域內(nèi)的4條特大規(guī)模冰川(巴爾托洛冰川面積為809.11 km2、Biafo冰川面積為559.81 km2、Panmah冰川面積為305.94 km2、Chogo Lungma冰川面積為287.61 km2)在研究時段內(nèi)保持穩(wěn)定狀態(tài),且冰川末端都有大面積表磧覆蓋。

    圖4 Shigar流域內(nèi)躍動冰川、前進冰川、退縮冰川和穩(wěn)定冰川的空間分布Fig.4 Spatial distribution of surge-type,advancing,retreating,and stable glaciers in the Shigar basin

    通過對不同類型冰川的最低海拔統(tǒng)計得出躍動冰川的最低海拔普遍高于其他類型冰川,前進冰川的最低海拔普遍高于退縮冰川,這可能是由于隨著海拔的增高氣溫降水條件改變,使得冰川積累量增加,進而對冰川運動狀態(tài)產(chǎn)生影響。同時,從不同海拔帶內(nèi)各類冰川面積變化率[圖5(a)]可看出,躍動冰川在海拔<5 000 m的范圍內(nèi)冰川面積增加,在海拔>5 000 m的范圍內(nèi)冰川面積變化微弱,面積變化率幾乎接近于0;前進冰川面積增加率隨著海拔升高不斷減少;退縮冰川在各海拔帶內(nèi)的面積均呈減少狀態(tài),并隨著海拔的升高退縮冰川面積減少率逐漸減小;各海拔帶內(nèi)的穩(wěn)定冰川面積變化微弱。同時冰川面積變化主要發(fā)生在低海拔地區(qū)。以上分析說明冰川的躍動和前進在很大程度上彌補了冰川退縮所造成的流域內(nèi)冰川面積減少。將冰川覆蓋范圍內(nèi)的地形平均坡度視為冰川坡度,對冰川坡度進行統(tǒng)計后發(fā)現(xiàn)流域內(nèi)冰川主要分布在坡度20°~30°的范圍內(nèi),但通過統(tǒng)計各坡度范圍內(nèi)各類冰川的數(shù)量占比[圖5(b)]發(fā)現(xiàn)冰川末端變化與冰川坡度并無顯著相關(guān)關(guān)系。

    圖5 躍動冰川、前進冰川、退縮冰川和穩(wěn)定冰川在不同海拔帶內(nèi)冰川面積變化率(a)和不同坡度等級內(nèi)冰川數(shù)量占比(b)Fig.5 Area change rates in different altitudes(a)and quantity percentage in different slopes(b)of surge-type,advancing,retreating,and stable glaciers in the Shigar basin

    3.2 冰川物質(zhì)平衡

    有研究[34]指出TanDEM-X DEM更適用于冰川表面高程變化研究,本文基于SRTM DEM和Tan-DEM-X DEM采用大地學測量方法對Shigar流域的冰川物質(zhì)平衡進行了估算。結(jié)果表明(圖6),在2000—2013年間,Shigar流域內(nèi)冰川表面高程變化速率為(-0.00±0.04)m·a-1,冰川物質(zhì)平衡估算結(jié)果為(-0.00±0.03)m w.e.·a-1,表明2000年以來Shigar流域冰川基本保持穩(wěn)定。其中,2000—2013年躍動冰川物質(zhì)平衡為(+0.17±0.03)m w.e.·a-1,前進冰川物質(zhì)平衡為(-0.01±0.03)m w.e.·a-1,退縮冰川物質(zhì)平衡為(-0.22±0.03)m w.e.·a-1,穩(wěn)定冰川物質(zhì)平衡為(-0.01±0.03)m w.e.·a-1,表明流域內(nèi)冰川物質(zhì)平衡保持穩(wěn)定主要是受躍動冰川物質(zhì)增加影響。

    圖6 2000—2013年Shigar流域冰川高程變化Fig.6 Glacier elevation change in the Shigar basin during 2000—2013

    流域內(nèi)不同海拔帶內(nèi)的冰川表面高程變化(圖7)表明,在海拔4 400~5 600 m的中高海拔地區(qū)冰川表面高程主要呈升高狀態(tài),在海拔<4 400 m的地區(qū)冰川表面高程降低,原因可能是隨著海拔增高氣溫逐漸降低,降水量增加,冰川積累量逐漸增加;而在海拔>5 600 m的高海拔地區(qū)冰川表面高程降低,原因可能有以下兩方面:一方面是在海拔>5 600 m的地區(qū)地形平均坡度>35°,在這樣的陡峭地區(qū)降水不易被保存,同時容易發(fā)生冰崩、雪崩等現(xiàn)象,從而使得冰川物質(zhì)減少,另一方面是在高海拔陡峭地區(qū)不同DEM數(shù)據(jù)之間的配準誤差較大。

    圖7 2000—2013年不同海拔帶內(nèi)冰川面積和表面高程變化Fig.7 The glacier area and elevation change in different altitudes during 2000—2013

    表3 中對不同學者的研究成果進行了對比,發(fā)現(xiàn)在不同的研究中喀喇昆侖山中部地區(qū)冰川均呈微弱正平衡或接近零平衡的狀態(tài)。不同研究結(jié)果的差異可能是研究區(qū)、研究時段、研究方法、轉(zhuǎn)換密度和SRTM C波段透射深度等的選擇不同造成的。但總體來說,研究區(qū)內(nèi)冬季降水量增加和夏季氣溫下降使得冰川積累量增多、消融量減少[17,22],冰川基本保持穩(wěn)定,進一步印證了“喀喇昆侖異常”。

    表3 喀喇昆侖地區(qū)冰川物質(zhì)平衡的不同研究結(jié)果Table 3 Different results of glacier mass balance of the Karakoram

    3.2.1 冰川高程變化特征分析

    對不同類型冰川的表面高程變化隨冰川長度的變化進行歸一化處理(圖中0代表冰川末端,1代表冰川頂部,冰川長度為冰川主流線)后發(fā)現(xiàn)(圖8),四類冰川表現(xiàn)出不同的表面高程變化模式:躍動冰川末端及消融區(qū)冰川表面顯著升高,冰川中上部呈微弱降低,而冰川上部微弱升高;前進冰川末端表面高程升高,但升高程度小于躍動冰川,且前進冰川中部微弱降低,中上部微弱升高,上部呈降低狀態(tài);退縮冰川末端表面高程降低,中部略有升高,上部則降低;穩(wěn)定冰川除上部表面高程降低較明顯外,其他部分冰川表面高程變化十分微弱。

    統(tǒng)計四類冰川在不同海拔帶內(nèi)和不同坡度區(qū)間的表面高程變化(圖9)可看出:躍動冰川在海拔3 600~4 600 m的地區(qū)表面高程明顯升高,在海拔5 000~5 800 m且坡度<45°的地區(qū)表面高程微弱升高,而在海拔>5 800 m的地區(qū)以及海拔4 600~5 800 m且坡度>50°的地區(qū)躍動冰川表面高程降低明顯;前進冰川表面高程升高主要發(fā)生在海拔<3 600 m的冰川末端和海拔4 400~5 600 m的地區(qū),表面高程降低主要發(fā)生在海拔5 600~5 800 m且坡度>45°的冰川上部地區(qū)和海拔3 600~4 200 m的地區(qū);退縮冰川表面高程升高主要發(fā)生在海拔4 600~5 400 m且坡度<40°的地區(qū),表面高程降低主要發(fā)生在海拔<4 400 m的冰川末端和冰川消融區(qū);穩(wěn)定冰川表面高程升高主要發(fā)生在海拔4 400~5 800 m且坡度<40°的地區(qū),表面高程降低主要發(fā)生在海拔<4 400 m且坡度<35°的地區(qū)和海拔>4 600 m且坡度>45°的地區(qū)。同時結(jié)合研究時段內(nèi)冰川表面高程變化隨冰川長度的變化可得出:在相同區(qū)域氣候條件下,躍動冰川物質(zhì)積累量大從而使得冰川表面高程升高明顯,同時推動躍動冰川末端向前運動;前進冰川上部物質(zhì)積累并且積累的物質(zhì)向下運動從而推動冰川末端不斷緩慢前進;退縮冰川物質(zhì)積累量少、冰川末端及消融區(qū)物質(zhì)虧損量大從而影響冰川末端逐漸退縮;穩(wěn)定冰川總體高程變化微弱,但也主要表現(xiàn)在冰川中上部呈物質(zhì)積累,冰川下部呈物質(zhì)虧損。

    圖8 躍動冰川(a)、前進冰川(b)、退縮冰川(c)和穩(wěn)定冰川(d)表面高程變化隨歸一化冰川長度的變化Fig.8 The normalized glaciers,length distribution of elevation changes of surge-type(a),advancing(b),retreating(c),and stable(d)glaciers

    4 結(jié)論

    本研究基于Landsat系列遙感影像、SRTM DEM、和TanDEM-X DEM等資料分析了喀喇昆侖山中部Shigar流域不同類型冰川的面積變化和物質(zhì)平衡,結(jié)果表明:

    (1)1993—2016年間Shigar流域冰川面積減少了(2.67±14.79)km2。其中,躍動冰川面積增加了1.30 km2;前進冰川面積增加了0.86 km2;退縮冰川面積減少了3.48 km2,穩(wěn)定冰川面積減少了1.34 km2。冰川前進和躍動在一定程度上抵消了冰川退縮所造成的冰川面積減少。

    (2)2000—2013年間Shigar流域冰川物質(zhì)平衡估算結(jié)果為(-0.00±0.03)m w.e.·a-1。躍動冰川物質(zhì)平衡為(+0.17±0.03)m w.e.·a-1,前進冰川物質(zhì)平衡為(-0.01±0.03)m w.e.·a-1,退縮冰川物質(zhì)平衡為(-0.22±0.03)m w.e.·a-1,穩(wěn)定冰川物質(zhì)平衡為(-0.01±0.03)m w.e.·a-1。流域內(nèi)冰川保持穩(wěn)定主要是受躍動冰川物質(zhì)增加影響。在同一流域內(nèi),冰川受相同區(qū)域氣候條件的影響,躍動冰川物質(zhì)積累量大、冰川表面高程升高明顯,推動部分躍動冰川末端向前運動;前進冰川主要表現(xiàn)為上部物質(zhì)積累并向下運動,推動冰川末端緩慢向前運動;退縮冰川末端及消融區(qū)物質(zhì)虧損量大,冰川末端逐漸退縮;穩(wěn)定冰川總體保持穩(wěn)定。

    (3)多數(shù)躍動冰川的冰川長度和冰川規(guī)模都集中在較高范圍內(nèi)。超過70%的前進冰川和退縮冰川的冰川長度在1~5 km的范圍內(nèi),73.53%的前進冰川分布在規(guī)模0.5~5 km2范圍內(nèi),50%的退縮冰川面積<1 km2,退縮冰川多為小規(guī)模。流域內(nèi)的4條特大規(guī)模冰川(巴爾托洛冰川、Biafo冰川、Panmah冰川、Chogo Lungma冰川)均為穩(wěn)定冰川,大規(guī)模冰川易保持穩(wěn)定,小規(guī)模冰川易發(fā)生改變。不同類型冰川的空間分布差異較大,但躍動冰川和前進冰川的空間分布具有一致性。海拔高度通過影響氣溫降水條件從而影響冰川運動狀態(tài),冰川面積變化主要發(fā)生在低海拔地區(qū)。流域內(nèi)冰川主要分布在坡度20°~30°的范圍內(nèi),但地形坡度對冰川運動狀態(tài)差異性的影響不顯著。

    猜你喜歡
    冰川海拔高程
    高海拔的甘南
    散文詩(2021年22期)2022-01-12 06:14:04
    8848.86m珠峰新高程
    當代陜西(2020年23期)2021-01-07 09:24:44
    為什么冰川會到處走?
    冰川會發(fā)出聲音嗎?
    南美三大高海拔球場
    海拔4000米的父愛
    長途跋涉到冰川
    GPS控制網(wǎng)的高程異常擬合與應(yīng)用
    SDCORS高程代替等級水準測量的研究
    世界上海拔最高、最雄偉的山脈
    河北遙感(2014年4期)2014-07-10 13:54:58
    又粗又爽又猛毛片免费看| 国产精品三级大全| 久久久久国产精品人妻aⅴ院| 特大巨黑吊av在线直播| 免费大片18禁| 午夜a级毛片| 婷婷精品国产亚洲av| 国产大屁股一区二区在线视频| 99riav亚洲国产免费| 亚洲av.av天堂| 日韩av不卡免费在线播放| 看片在线看免费视频| 成年免费大片在线观看| 色播亚洲综合网| 97超视频在线观看视频| 97碰自拍视频| 伦精品一区二区三区| 最后的刺客免费高清国语| 五月玫瑰六月丁香| 人人妻人人澡欧美一区二区| 成人永久免费在线观看视频| 91精品国产九色| 天天躁日日操中文字幕| 在线免费观看的www视频| 国产真实伦视频高清在线观看| 看片在线看免费视频| 日本撒尿小便嘘嘘汇集6| 免费人成视频x8x8入口观看| 精品乱码久久久久久99久播| 免费观看在线日韩| 日韩在线高清观看一区二区三区| 蜜桃久久精品国产亚洲av| 久久午夜亚洲精品久久| 97超视频在线观看视频| 午夜福利高清视频| 久久久久久久久久成人| 国产大屁股一区二区在线视频| 久久精品夜色国产| 日本爱情动作片www.在线观看 | 国产黄a三级三级三级人| av天堂在线播放| 变态另类成人亚洲欧美熟女| 日韩大尺度精品在线看网址| 精品人妻熟女av久视频| 久久综合国产亚洲精品| 丰满的人妻完整版| 国产精品人妻久久久影院| 18禁在线播放成人免费| 毛片女人毛片| 日本黄色片子视频| 亚洲欧美日韩卡通动漫| 色综合亚洲欧美另类图片| 亚洲av中文av极速乱| 伦理电影大哥的女人| 淫妇啪啪啪对白视频| 成人鲁丝片一二三区免费| 99riav亚洲国产免费| 成人午夜高清在线视频| 噜噜噜噜噜久久久久久91| 18禁黄网站禁片免费观看直播| 蜜臀久久99精品久久宅男| 亚洲专区国产一区二区| 成人欧美大片| 成年av动漫网址| 最近在线观看免费完整版| 亚洲欧美日韩卡通动漫| 日韩制服骚丝袜av| 嫩草影视91久久| 日日摸夜夜添夜夜添av毛片| 亚洲av免费高清在线观看| av免费在线看不卡| 在线观看av片永久免费下载| av福利片在线观看| 小蜜桃在线观看免费完整版高清| 国产在线精品亚洲第一网站| 一级av片app| 精品久久久久久久末码| 狠狠狠狠99中文字幕| 变态另类丝袜制服| 国产三级在线视频| 国产黄片美女视频| 国产激情偷乱视频一区二区| av黄色大香蕉| 欧美极品一区二区三区四区| 久久国内精品自在自线图片| 免费观看人在逋| 久久亚洲精品不卡| 中文亚洲av片在线观看爽| 日韩一区二区视频免费看| 日本一本二区三区精品| 国产单亲对白刺激| 18禁裸乳无遮挡免费网站照片| 亚洲性夜色夜夜综合| 国产精品久久久久久精品电影| 人人妻,人人澡人人爽秒播| 熟女电影av网| 三级男女做爰猛烈吃奶摸视频| 99热精品在线国产| 国产中年淑女户外野战色| 亚洲欧美清纯卡通| 男人的好看免费观看在线视频| 免费人成视频x8x8入口观看| 亚洲无线在线观看| 精品乱码久久久久久99久播| 亚洲精品一区av在线观看| 三级经典国产精品| 亚洲精华国产精华液的使用体验 | 激情 狠狠 欧美| 波多野结衣高清作品| 长腿黑丝高跟| 午夜福利18| 人人妻人人澡欧美一区二区| 成人av一区二区三区在线看| av天堂中文字幕网| 全区人妻精品视频| 亚洲电影在线观看av| 国语自产精品视频在线第100页| 国内精品一区二区在线观看| 国产毛片a区久久久久| 中文亚洲av片在线观看爽| 成人特级黄色片久久久久久久| 国产视频一区二区在线看| 老熟妇乱子伦视频在线观看| 国产综合懂色| АⅤ资源中文在线天堂| 美女免费视频网站| 波多野结衣高清作品| 在线观看一区二区三区| 小说图片视频综合网站| av福利片在线观看| 看十八女毛片水多多多| 直男gayav资源| 韩国av在线不卡| 一进一出抽搐动态| 日韩亚洲欧美综合| 亚洲婷婷狠狠爱综合网| 久久韩国三级中文字幕| 黄色视频,在线免费观看| 国产亚洲91精品色在线| 久久中文看片网| 久久久成人免费电影| 国产蜜桃级精品一区二区三区| 一级a爱片免费观看的视频| 亚洲自拍偷在线| 最好的美女福利视频网| 波多野结衣高清无吗| 成人毛片a级毛片在线播放| 国内少妇人妻偷人精品xxx网站| 欧美日韩综合久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久久久久久久免| 日本黄色片子视频| 熟女人妻精品中文字幕| 成人性生交大片免费视频hd| 少妇熟女欧美另类| 午夜精品国产一区二区电影 | 成年版毛片免费区| 全区人妻精品视频| 欧美激情久久久久久爽电影| 看十八女毛片水多多多| 亚洲中文字幕一区二区三区有码在线看| 国产av在哪里看| 国产男靠女视频免费网站| 无遮挡黄片免费观看| 日韩欧美精品v在线| 久久国产乱子免费精品| eeuss影院久久| 日本色播在线视频| 亚洲av美国av| 国产欧美日韩精品亚洲av| 18禁在线播放成人免费| 亚洲精品影视一区二区三区av| 一区二区三区免费毛片| 国产人妻一区二区三区在| a级毛色黄片| 女人被狂操c到高潮| 免费看日本二区| 老司机福利观看| 99久久精品一区二区三区| 国产激情偷乱视频一区二区| 亚洲无线观看免费| 久久久久免费精品人妻一区二区| 赤兔流量卡办理| 亚洲在线观看片| 国内精品久久久久精免费| 欧美+日韩+精品| 国产女主播在线喷水免费视频网站 | 亚洲自拍偷在线| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品亚洲一区二区| 一区二区三区高清视频在线| 国产精品综合久久久久久久免费| 日韩一本色道免费dvd| 人人妻人人澡欧美一区二区| eeuss影院久久| 日韩欧美三级三区| 久久精品综合一区二区三区| 国产精品女同一区二区软件| 丝袜美腿在线中文| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久av不卡| 高清毛片免费观看视频网站| 99在线视频只有这里精品首页| 简卡轻食公司| 亚洲人与动物交配视频| 日韩亚洲欧美综合| 亚洲人与动物交配视频| 国产大屁股一区二区在线视频| 亚洲av成人精品一区久久| 久久久久久久久中文| 99久久九九国产精品国产免费| 在线观看免费视频日本深夜| 日本黄大片高清| 国产男靠女视频免费网站| 国产蜜桃级精品一区二区三区| 欧美一区二区国产精品久久精品| 国产精品女同一区二区软件| 久久久久国内视频| 亚洲av免费高清在线观看| 中文字幕精品亚洲无线码一区| 三级毛片av免费| 性欧美人与动物交配| 国产大屁股一区二区在线视频| 美女被艹到高潮喷水动态| 99热精品在线国产| 日韩在线高清观看一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 精品福利观看| 亚洲自拍偷在线| 国产一区二区亚洲精品在线观看| 99国产精品一区二区蜜桃av| 精品欧美国产一区二区三| 国产久久久一区二区三区| 国产精品亚洲美女久久久| 成人亚洲欧美一区二区av| 日韩人妻高清精品专区| 97人妻精品一区二区三区麻豆| 高清日韩中文字幕在线| 91久久精品国产一区二区三区| 91久久精品国产一区二区成人| 最后的刺客免费高清国语| 51国产日韩欧美| 午夜福利成人在线免费观看| 亚洲av成人精品一区久久| 熟妇人妻久久中文字幕3abv| 成年女人看的毛片在线观看| 无遮挡黄片免费观看| 99热只有精品国产| 亚洲性夜色夜夜综合| 日韩高清综合在线| 国产视频内射| 免费av毛片视频| 国产精品一区二区三区四区久久| 亚洲性夜色夜夜综合| 天天躁日日操中文字幕| 亚洲国产精品成人久久小说 | 热99在线观看视频| 精品久久国产蜜桃| 性色avwww在线观看| 性色avwww在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲欧美清纯卡通| 少妇人妻精品综合一区二区 | 日本a在线网址| 老司机影院成人| 高清毛片免费看| 午夜视频国产福利| 在线观看免费视频日本深夜| 亚洲真实伦在线观看| 亚洲18禁久久av| 校园人妻丝袜中文字幕| 久久久国产成人免费| 别揉我奶头 嗯啊视频| 中文字幕av在线有码专区| 精品熟女少妇av免费看| 高清日韩中文字幕在线| 欧美成人a在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲av成人精品一区久久| 黄色配什么色好看| 国产高潮美女av| 久久久久久伊人网av| 成人漫画全彩无遮挡| 高清毛片免费观看视频网站| 人妻久久中文字幕网| 国产午夜精品论理片| 人妻夜夜爽99麻豆av| 在线免费观看的www视频| 高清毛片免费观看视频网站| 校园春色视频在线观看| 国产精品一区二区三区四区免费观看 | 国产女主播在线喷水免费视频网站 | 国产探花在线观看一区二区| 亚洲欧美日韩东京热| 99国产极品粉嫩在线观看| 国产私拍福利视频在线观看| 黄色视频,在线免费观看| 午夜精品在线福利| 天堂影院成人在线观看| 九九在线视频观看精品| 国产精品无大码| 免费高清视频大片| 嫩草影院新地址| 别揉我奶头~嗯~啊~动态视频| 九九爱精品视频在线观看| 欧美+日韩+精品| av.在线天堂| 一区福利在线观看| 亚洲国产精品成人综合色| 可以在线观看的亚洲视频| 99久久成人亚洲精品观看| 中国美白少妇内射xxxbb| 一个人看的www免费观看视频| 久久久午夜欧美精品| 久久久精品94久久精品| 性欧美人与动物交配| 日韩成人伦理影院| 国产成人精品久久久久久| 色噜噜av男人的天堂激情| 久久综合国产亚洲精品| 老司机影院成人| 国产淫片久久久久久久久| 一本久久中文字幕| 亚洲精品影视一区二区三区av| 最近中文字幕高清免费大全6| 18禁在线播放成人免费| 噜噜噜噜噜久久久久久91| 国产精品爽爽va在线观看网站| 久久人人爽人人爽人人片va| 五月玫瑰六月丁香| 久久久久久伊人网av| 国产成人91sexporn| 在线观看一区二区三区| eeuss影院久久| www日本黄色视频网| 嫩草影视91久久| 亚洲激情五月婷婷啪啪| 日韩一本色道免费dvd| 在线观看av片永久免费下载| 午夜福利在线观看免费完整高清在 | 久久精品人妻少妇| 国产精品,欧美在线| 俺也久久电影网| 麻豆精品久久久久久蜜桃| 国产 一区 欧美 日韩| 看十八女毛片水多多多| 日韩中字成人| 综合色av麻豆| 亚洲av免费高清在线观看| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 亚洲精品在线观看二区| 色哟哟·www| 午夜精品在线福利| 精品久久久久久久久久久久久| 12—13女人毛片做爰片一| 一区福利在线观看| 国产精品女同一区二区软件| 三级毛片av免费| 在线播放无遮挡| 久久精品人妻少妇| 精品久久国产蜜桃| 成人鲁丝片一二三区免费| 看黄色毛片网站| 亚洲美女黄片视频| 免费大片18禁| 高清毛片免费看| 不卡一级毛片| 国产一区二区三区在线臀色熟女| 人妻少妇偷人精品九色| 美女被艹到高潮喷水动态| 婷婷精品国产亚洲av在线| 亚洲最大成人av| 国产视频内射| 国产伦在线观看视频一区| 亚洲精品一区av在线观看| 舔av片在线| 国产视频内射| 尾随美女入室| 国产三级在线视频| 亚洲av一区综合| 91精品国产九色| 日韩,欧美,国产一区二区三区 | 搡老妇女老女人老熟妇| 欧美丝袜亚洲另类| 亚洲丝袜综合中文字幕| 欧美一区二区精品小视频在线| 99在线人妻在线中文字幕| 成人av在线播放网站| 99精品在免费线老司机午夜| 免费在线观看影片大全网站| 久久久久久国产a免费观看| 国产淫片久久久久久久久| 啦啦啦观看免费观看视频高清| 成年女人看的毛片在线观看| 97超视频在线观看视频| 亚洲四区av| 啦啦啦观看免费观看视频高清| 三级国产精品欧美在线观看| 亚洲18禁久久av| 亚洲欧美中文字幕日韩二区| 天堂av国产一区二区熟女人妻| 久久欧美精品欧美久久欧美| 听说在线观看完整版免费高清| 99热全是精品| 亚洲无线观看免费| 午夜福利18| 午夜亚洲福利在线播放| 中国美女看黄片| 成人性生交大片免费视频hd| 久久精品国产99精品国产亚洲性色| 色尼玛亚洲综合影院| 成人综合一区亚洲| 不卡视频在线观看欧美| 你懂的网址亚洲精品在线观看 | 亚洲美女黄片视频| 十八禁网站免费在线| 嫩草影院精品99| 嫩草影院入口| 美女 人体艺术 gogo| 国产又黄又爽又无遮挡在线| 欧美激情国产日韩精品一区| 久久精品91蜜桃| 精华霜和精华液先用哪个| 亚洲一级一片aⅴ在线观看| 亚洲五月天丁香| 亚洲久久久久久中文字幕| 九九在线视频观看精品| 亚洲国产欧美人成| 啦啦啦啦在线视频资源| 亚洲欧美日韩东京热| 久久精品夜夜夜夜夜久久蜜豆| 午夜老司机福利剧场| 欧美一级a爱片免费观看看| 日日啪夜夜撸| 91久久精品国产一区二区成人| 国产大屁股一区二区在线视频| 国内少妇人妻偷人精品xxx网站| 午夜亚洲福利在线播放| 亚洲美女搞黄在线观看 | 一个人免费在线观看电影| 99久久成人亚洲精品观看| ponron亚洲| 欧美激情在线99| 精品久久久噜噜| 亚洲欧美日韩无卡精品| 一区二区三区免费毛片| 国产成人a区在线观看| 久久欧美精品欧美久久欧美| 少妇猛男粗大的猛烈进出视频 | 亚州av有码| 午夜激情福利司机影院| 欧美不卡视频在线免费观看| 黑人高潮一二区| 色尼玛亚洲综合影院| 老熟妇仑乱视频hdxx| 少妇裸体淫交视频免费看高清| 蜜桃久久精品国产亚洲av| 色哟哟哟哟哟哟| 超碰av人人做人人爽久久| 国产综合懂色| 亚洲精品日韩在线中文字幕 | 中文在线观看免费www的网站| 久久久久久久久大av| 99久久无色码亚洲精品果冻| 国产 一区 欧美 日韩| 日产精品乱码卡一卡2卡三| 性插视频无遮挡在线免费观看| 成人美女网站在线观看视频| 午夜福利18| 午夜免费激情av| 久久久久免费精品人妻一区二区| 国产蜜桃级精品一区二区三区| 国产成人aa在线观看| 亚洲av二区三区四区| 舔av片在线| 可以在线观看毛片的网站| 久久精品国产亚洲av香蕉五月| 春色校园在线视频观看| 成年版毛片免费区| 天堂√8在线中文| 黑人高潮一二区| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| 亚洲国产精品久久男人天堂| 亚洲国产精品sss在线观看| 国产精品伦人一区二区| 人人妻人人澡欧美一区二区| 久久精品国产亚洲网站| 寂寞人妻少妇视频99o| 国产亚洲av嫩草精品影院| 精品福利观看| 亚洲欧美日韩卡通动漫| 小蜜桃在线观看免费完整版高清| 有码 亚洲区| 午夜精品一区二区三区免费看| 国产成人91sexporn| 亚洲国产精品国产精品| www.色视频.com| 亚洲国产精品久久男人天堂| 热99在线观看视频| 女的被弄到高潮叫床怎么办| 免费看a级黄色片| 毛片女人毛片| 小说图片视频综合网站| 久久天躁狠狠躁夜夜2o2o| 日韩国内少妇激情av| 欧美在线一区亚洲| 精品福利观看| 国产成人影院久久av| 成人亚洲欧美一区二区av| 欧美xxxx黑人xx丫x性爽| 国产高清不卡午夜福利| 91在线精品国自产拍蜜月| 尾随美女入室| 精品久久久久久久久av| 国产伦精品一区二区三区四那| 变态另类丝袜制服| 精品国内亚洲2022精品成人| 免费高清视频大片| 成人综合一区亚洲| 国产高清视频在线观看网站| 国内久久婷婷六月综合欲色啪| 日韩欧美一区二区三区在线观看| 特大巨黑吊av在线直播| 国产真实伦视频高清在线观看| 久久人人爽人人片av| 小说图片视频综合网站| 插阴视频在线观看视频| 亚洲av二区三区四区| 亚洲欧美日韩卡通动漫| 欧美+日韩+精品| 亚洲18禁久久av| 女的被弄到高潮叫床怎么办| 欧美bdsm另类| 91精品国产九色| 99国产极品粉嫩在线观看| 日韩亚洲欧美综合| 欧美区成人在线视频| 天美传媒精品一区二区| 国产亚洲精品综合一区在线观看| 亚洲熟妇熟女久久| 国产av麻豆久久久久久久| 国产成人福利小说| 日本一本二区三区精品| 亚洲精品色激情综合| 99久久成人亚洲精品观看| 久久久成人免费电影| 日本在线视频免费播放| avwww免费| 国产v大片淫在线免费观看| 精华霜和精华液先用哪个| 久久久久性生活片| 22中文网久久字幕| 日韩亚洲欧美综合| 不卡一级毛片| 午夜福利视频1000在线观看| 女同久久另类99精品国产91| 免费人成视频x8x8入口观看| 一个人看视频在线观看www免费| 久久综合国产亚洲精品| 欧美日本视频| 国产一区二区三区在线臀色熟女| 少妇的逼水好多| 女生性感内裤真人,穿戴方法视频| 精品国内亚洲2022精品成人| 亚洲国产欧洲综合997久久,| 天堂网av新在线| 国产成年人精品一区二区| 日韩亚洲欧美综合| 国产精品乱码一区二三区的特点| 校园人妻丝袜中文字幕| 成人永久免费在线观看视频| 激情 狠狠 欧美| 在线播放无遮挡| 亚洲性夜色夜夜综合| 日韩成人伦理影院| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| a级毛片a级免费在线| 人人妻人人看人人澡| 综合色av麻豆| 日韩在线高清观看一区二区三区| 久久久久久大精品| 日韩欧美在线乱码| 亚洲最大成人av| 亚洲国产日韩欧美精品在线观看| 免费观看人在逋| 日本熟妇午夜| 欧美在线一区亚洲| 成年版毛片免费区| 欧美xxxx黑人xx丫x性爽| 变态另类成人亚洲欧美熟女| 97热精品久久久久久| 小说图片视频综合网站| 最后的刺客免费高清国语| 国产亚洲精品久久久久久毛片| 一边摸一边抽搐一进一小说| 99视频精品全部免费 在线| 午夜免费激情av| 身体一侧抽搐| 99久久精品国产国产毛片| 国产免费一级a男人的天堂| 麻豆国产av国片精品| 99久久久亚洲精品蜜臀av| 最近2019中文字幕mv第一页| 蜜臀久久99精品久久宅男| 国产69精品久久久久777片| 国产 一区精品| 一卡2卡三卡四卡精品乱码亚洲| 搡女人真爽免费视频火全软件 | 日本三级黄在线观看| www.色视频.com| 老女人水多毛片|