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

    基于冰川流速的喀喇昆侖地區(qū)典型冰川厚度反演與冰儲(chǔ)量估算

    2023-10-11 09:41:40王紀(jì)碩張齊民孫一丹張豪磊張子彥閆世勇
    關(guān)鍵詞:冰體冰川反演

    王紀(jì)碩, 張齊民, 孫一丹, 張豪磊, 張子彥, 閆世勇

    (1.中國(guó)礦業(yè)大學(xué) 環(huán)境與測(cè)繪學(xué)院,江蘇 徐州 221116; 2.自然資源部國(guó)土環(huán)境與災(zāi)害監(jiān)測(cè)重點(diǎn)實(shí)驗(yàn)室,江蘇 徐州 221116; 3.山東省國(guó)土測(cè)繪院 基礎(chǔ)測(cè)繪中心,山東 濟(jì)南 250102)

    在當(dāng)前全球變暖的背景下[1],以青藏高原為核心的“亞洲水塔”正在發(fā)生以失衡為特征的劇烈變化,青藏高原山地冰川加速退縮,對(duì)人類社會(huì)文明和生態(tài)安全的影響不可低估[2]。青藏高原喀喇昆侖山是全球山岳冰川最發(fā)達(dá)的高大山脈,也是中緯度地區(qū)現(xiàn)代冰川最發(fā)育的地帶,而冰川作為當(dāng)?shù)刂匾牡Y源庫(kù),其厚度與儲(chǔ)量的變化對(duì)周邊區(qū)域的生態(tài)平衡起著至關(guān)重要的作用[3]。

    獲取冰川厚度的方法主要有實(shí)地測(cè)量和間接估計(jì)2類。采用人工鉆孔法進(jìn)行實(shí)地測(cè)量,會(huì)受到地理環(huán)境限制,無(wú)法進(jìn)行大范圍作業(yè)。根據(jù)厚度-面積比例關(guān)系式可間接估計(jì)冰川的平均厚度,但由于不同冰川各自的屬性不同,使得該方法針對(duì)區(qū)域性冰川有效,對(duì)于單一冰川并不適用。20世紀(jì)80年代初期,我國(guó)首次采用探地雷達(dá)(ground penetrating radar,GPR)對(duì)天山烏魯木齊河源1號(hào)冰川進(jìn)行冰厚測(cè)量并取得一些成果,隨著GPR技術(shù)的不斷發(fā)展,該方法被廣泛應(yīng)用在我國(guó)青藏高原山地冰川的冰厚提取與冰儲(chǔ)量估算研究中,但GPR探測(cè)冰厚的準(zhǔn)確性受到電磁波穿透冰體能力和雷達(dá)系統(tǒng)圖像處理能力的制約,并且冰層的細(xì)微變化會(huì)分散GPR的信號(hào),從而降低圖像清晰度,當(dāng)前GPR技術(shù)多應(yīng)用于小型冰川厚度提取,在青藏高原大規(guī)模冰川厚度測(cè)量中的應(yīng)用很少[4]。衛(wèi)星遙感技術(shù)的發(fā)展使得大規(guī)模監(jiān)測(cè)冰川成為可能,在山地冰川應(yīng)用較多的是美國(guó)航天飛機(jī)雷達(dá)地形測(cè)繪使命(Shuttle Radar Topography Mission,SRTM) 數(shù)字高程模型(digital elevation model,DEM)數(shù)據(jù)和Landsat衛(wèi)星光學(xué)影像。文獻(xiàn)[5]根據(jù)淺冰近似(shallow ice approximation,SIA)理論,采用SRTM DEM數(shù)據(jù)進(jìn)行冰川厚度反演。SIA理論是基于冰川中流線平均基底剪切應(yīng)力經(jīng)驗(yàn)公式和冰川中流線平均坡度進(jìn)行冰川厚度反演,然而該方法只能求得冰川的平均冰厚和總體積,并且由于DEM數(shù)據(jù)獲取時(shí)間具有單一性,難以獲取冰川坡度的年際變化,同時(shí)冰川表面運(yùn)動(dòng)、溫度、降水等因素也都嚴(yán)重影響冰川厚度的變化,因而不能為冰川厚度的持續(xù)性變化提供參考。文獻(xiàn)[6]采用Landsat衛(wèi)星影像獲取青藏高原東南部冰川的表面狀態(tài),根據(jù)質(zhì)量平衡分布數(shù)據(jù)、冰川表面速度和冰流原理估算出冰川冰厚的空間分布;文獻(xiàn)[7]采用Landsat 7、Sentinel-2衛(wèi)星影像對(duì)喜馬拉雅山脈冰川的流速進(jìn)行反演,探索冰川動(dòng)態(tài)變化與質(zhì)量守恒的相互作用機(jī)制。但光學(xué)影像在青藏高原冰川監(jiān)測(cè)中的應(yīng)用極大地受限于惡劣的天氣狀況[8],而雷達(dá)影像具有穿透力強(qiáng)、覆蓋范圍廣、不受云雨天氣影響等優(yōu)點(diǎn),利用不同時(shí)期的雷達(dá)影像,基于冰川流速對(duì)冰川厚度進(jìn)行估算,能夠簡(jiǎn)便高效地通過(guò)數(shù)值模擬的方法實(shí)現(xiàn)冰川厚度的反演。本文基于Sentinel-1A衛(wèi)星干涉寬幅(Interferometric Wide swath,IW)模式下的單視復(fù)數(shù)影像和SRTM DEM數(shù)據(jù),采用像素跟蹤算法和層流理論進(jìn)行冰川厚度反演與冰儲(chǔ)量估算,以認(rèn)識(shí)喀喇昆侖地區(qū)冰川厚度空間分布及冰厚與氣候之間的響應(yīng)關(guān)系;將估算結(jié)果與法國(guó)南部比利牛斯天文臺(tái)數(shù)據(jù)服務(wù)中心的全球冰川厚度數(shù)據(jù)產(chǎn)品進(jìn)行對(duì)比分析,從而驗(yàn)證本文方法的可行性。研究結(jié)果可為研究喀喇昆侖地區(qū)冰川厚度變化提供參考。

    1 研究區(qū)概況和數(shù)據(jù)來(lái)源

    喀喇昆侖山是全球山地冰川最發(fā)達(dá)的高大山脈,從阿富汗最東部向東南延伸約480 km,平均海拔超過(guò)6 000 m,山區(qū)終年低溫,年平均氣溫低于0 ℃,由于冬春兩季受到西風(fēng)大環(huán)流的影響,冰川成冰區(qū)域年降水量超過(guò)1 000 mm,獨(dú)特的地理環(huán)境使得喀喇昆侖山有超過(guò)37%的面積被冰川覆蓋。在當(dāng)前氣候變暖的形勢(shì)下,全球絕大多數(shù)冰川出現(xiàn)萎縮現(xiàn)象,而喀喇昆侖地區(qū)冰川對(duì)氣候變化響應(yīng)不夠顯著,在某些情況下出現(xiàn)末端前進(jìn)、面積增加的現(xiàn)象,被稱為“喀喇昆侖異?!盵9],本文選取喀喇昆侖地區(qū)卻哥隆瑪冰川、錫亞琴冰川和巴爾托洛冰川3條規(guī)模較大且具代表性的冰川作為研究對(duì)象。

    為研究喀喇昆侖地區(qū)冰川厚度分布,采用6景Sentinel-1A衛(wèi)星降軌單視復(fù)數(shù)影像,時(shí)間跨度為2019年1月—2020年1月,影像距離向和方位向空間分辨率分別約為5、20 m,觀測(cè)范圍覆蓋研究區(qū)目標(biāo)冰川,影像數(shù)據(jù)全部來(lái)自歐洲航天局(European Space Agency,ESA) (https://search.asf.alaska.edu),具體影像信息見(jiàn)表1所列。另外,美國(guó)地質(zhì)調(diào)查局90 m空間分辨率的SRTM DEM數(shù)據(jù)主要用于研究區(qū)地形位移量補(bǔ)償、地理編碼及坡度計(jì)算;第6版蘭多夫冰川編目 (Randolph Glacier Inventory 6.0,RGI 6.0)用于提取研究區(qū)冰川矢量邊界,數(shù)據(jù)來(lái)源于全球陸地冰空間測(cè)量計(jì)劃(http://www.glims.org/RGI/);冰川厚度數(shù)據(jù)產(chǎn)品(https://doi.org/10.6096/1007)[10]作為驗(yàn)證數(shù)據(jù)用于結(jié)果的精度評(píng)定。

    表1 研究區(qū)影像信息

    2 研究方法

    為獲取目標(biāo)冰川流速,并在其基礎(chǔ)上獲取冰川厚度信息,本研究采用像素跟蹤算法和層流理論進(jìn)行冰川流速提取與冰川厚度的反演,通過(guò)與全球冰川厚度產(chǎn)品進(jìn)行對(duì)比分析,驗(yàn)證該方法的反演精度,相應(yīng)技術(shù)流程如圖1所示。圖1中,NCC表示歸一化互相關(guān)(normalized cross-correlation)。

    圖1 冰川厚度反演技術(shù)流程

    2.1 像素跟蹤算法

    對(duì)于流速較快的山地冰川,像素跟蹤算法具有良好的適用性[11],本研究利用覆蓋研究區(qū)的典型冰川多期Sentinel-1A衛(wèi)星遙感影像,采用改進(jìn)的像素跟蹤算法進(jìn)行位移量計(jì)算,該方法以NCC算法為基礎(chǔ),通過(guò)主、輔影像的非冰川區(qū)穩(wěn)定像素的亞像素級(jí)配準(zhǔn)結(jié)果精確解算主、輔影像之間的幾何變換模型,進(jìn)而估計(jì)冰川表面運(yùn)動(dòng)速度分布。NCC的計(jì)算準(zhǔn)則公式為:

    (1)

    通過(guò)上述復(fù)雜的匹配過(guò)程,能夠獲取研究區(qū)總體位移量,總體位移量可分為冰川運(yùn)動(dòng)位移量和非冰川位移量,其計(jì)算公式為:

    D=Dglacier+Dnonglacier

    (2)

    其中:D為研究區(qū)總體位移量;Dglacier為冰川運(yùn)動(dòng)位移量;Dnonglacier為非冰川位移量,其中包含合成孔徑雷達(dá)(synthetic aperture radar,SAR)影像軌道不同引起的全局位移量、電離層位移量、噪聲位移量和地形高低起伏產(chǎn)生的位移量。為去除SAR影像軌道不同引起的全局位移量,采用隨機(jī)抽樣一致性算法來(lái)估算變換模型,獲得變換模型參數(shù)去除軌道相關(guān)位移量;由于喀喇昆侖地區(qū)處于中低緯度,而電離層引起的位移量對(duì)于高緯度地區(qū)較為明顯,因此電離層影響可以忽略不計(jì);為去除圖像配準(zhǔn)過(guò)程中產(chǎn)生的噪聲點(diǎn),在圖像后處理過(guò)程中采用自適應(yīng)中值濾波進(jìn)行濾波處理;冰川區(qū)地形陡峭引起的位移量受到空間基線長(zhǎng)度的影響,采用DEM數(shù)據(jù)和基線信息對(duì)產(chǎn)生的地形位移量進(jìn)行補(bǔ)償。在考慮非冰川位移量的基礎(chǔ)上進(jìn)行位移合成和地理編碼,即可得到冰川區(qū)表面流速分布。

    2.2 層流理論

    采用層流理論進(jìn)行冰川厚度的反演[12],該方法的基本原理是綜合考慮冰川坡度、基底滑動(dòng)、自身重力等因素對(duì)冰川厚度的影響,根據(jù)輸入?yún)?shù)反演冰川厚度,如圖2所示。

    圖2 層流理論示意圖

    圖2中:H為冰厚;τ為基底剪應(yīng)力;F為冰體沿斜面方向向下滑動(dòng)的力,即重力沿斜面的分力;Flim為冰體靜止時(shí)的最大摩擦限度;m為冰體質(zhì)量。隨著冰體質(zhì)量增加,當(dāng)斜面方向向下滑動(dòng)的力超過(guò)摩擦限度,就會(huì)產(chǎn)生運(yùn)動(dòng)。

    假設(shè)冰體底面積為S,冰體在斜面固定的力為τS,基底剪應(yīng)力等于冰體重力沿斜面方向的分力,即

    τS=mgsinα,m=ρHS

    (3)

    其中:g為重力加速度;α為冰川表面坡度[6];ρ為冰密度。

    冰-巖界面條件變化產(chǎn)生的底部摩擦力不可忽略,需要增加形狀因子f′,則有:

    τ=f′ρgHsinα

    (4)

    其中,f′與山谷冰川橫斷面結(jié)構(gòu)有關(guān),山谷冰川f′取值范圍[13]為0.5~0.9。

    根據(jù)層流理論,冰川表面運(yùn)動(dòng)速率us計(jì)算公式為:

    (5)

    其中:ub為冰川基底運(yùn)動(dòng)速率;A為蠕變參數(shù),與環(huán)境溫度、冰體雜質(zhì)等因素相關(guān)[14];n為格倫定律指數(shù),在1.5~4.2范圍內(nèi)變動(dòng),與冰川所受的應(yīng)力有關(guān)[15]。

    結(jié)合式(4)、式(5)可以得出:

    (6)

    其中,k為基底滑動(dòng)系數(shù)。

    3 數(shù)據(jù)處理與結(jié)果分析

    3.1 冰川流速提取與結(jié)果分析

    利用2019年1月—2020年1月的6景降軌SAR數(shù)據(jù),采用像素跟蹤算法進(jìn)行冰川流速提取,設(shè)置主影像模板窗口的大小為96,輔影像搜索窗口的大小為32,步長(zhǎng)為2,可以在冰川區(qū)達(dá)到較高的配準(zhǔn)精度。為削弱噪聲的影響,采用自適應(yīng)中值濾波進(jìn)行降噪處理,設(shè)置最小濾波半徑為7,最大濾波半徑為15,得到目標(biāo)冰川平均流速如圖3所示。圖3中,冰川運(yùn)動(dòng)方向如紅色箭頭所示。

    圖3 3條冰川流速分布圖

    由于冰川底部基巖本身的差異及冰川受山谷拖滯作用影響,冰川的形態(tài)呈現(xiàn)多樣性、復(fù)雜性。從圖3可以看出:卻哥隆瑪冰川和巴爾托洛冰川最大速率位于冰川上游(約為38 cm/d),隨著冰川坡度不斷減緩,冰川流速逐漸放緩,但均呈上游積累區(qū)相對(duì)于下游消融區(qū)流速較快、冰川中心相對(duì)于冰川邊緣流速較快的特點(diǎn),這是冰體下滑過(guò)程中與周邊山體發(fā)生摩擦導(dǎo)致速度減緩,另外,表面冰磧物的阻熱特性一定程度上減緩了冰川的消融、阻礙冰川運(yùn)動(dòng);而錫亞琴冰川末端出現(xiàn)流速激增的情況(約為40 cm/d),相對(duì)于上游流速明顯增大,期間可能發(fā)生冰川躍動(dòng),這也是喀喇昆侖地區(qū)躍動(dòng)冰川常見(jiàn)現(xiàn)象[16],此外冰川流速變化符合“上游向下游逐漸減緩、中央向兩側(cè)逐漸減緩”的運(yùn)動(dòng)規(guī)律。

    3.2 冰川厚度反演與結(jié)果分析

    為分析冰川厚度分布,在獲取流速的基礎(chǔ)上,采用層流理論進(jìn)行冰川厚度反演,其中蠕變系數(shù)A=2.4×10-24/(Pa3·s),滑動(dòng)系數(shù)k=0.2,冰密度ρ=917 kg/m3,重力加速度[10]g=9.81 m/s2;綜合考慮喀喇昆侖山脈所處氣候條件及冰川結(jié)構(gòu)等影響因素,將格倫定律指數(shù)n取值為2.6,形狀因子f′取值為0.6。目標(biāo)冰川厚度分布如圖4所示。圖4中,a1—a2、b1—b2、c1—c2表示冰川剖面線始末位置。

    圖4 3條冰川厚度反演圖

    由圖4可知:卻哥隆瑪冰川的冰厚最大值出現(xiàn)在支流與干流交匯口,約為380 m,并且大部分區(qū)域厚度為220~250 m,占冰川面積的60.04%;錫亞琴冰川最大冰厚出現(xiàn)在冰川干流下游,約為380 m,大部分區(qū)域冰厚為180~220 m,占冰川面積的42.35%;而巴爾托洛冰川冰厚最大值集中在冰川干流上游,隨著海拔不斷升高,冰厚在210~300 m范圍內(nèi)不斷增加,且冰川大部分區(qū)域冰厚為170~200 m,占冰川面積的57.09%??傮w而言,3條冰川均呈從上游到下游逐漸減薄、支流向干流不斷增厚的趨勢(shì)。

    為進(jìn)一步分析地形因素對(duì)冰川厚度的影響,本文沿著圖4中的剖面線提取3條冰川的冰厚并進(jìn)行擬合,得到剖面圖,如圖5所示。

    圖5 3條冰川干流高程和冰厚變化曲線

    從圖5可以看出,卻哥隆瑪冰川冰厚存在2個(gè)峰值,分析其原因,有以下2點(diǎn):① 可能是冰川支流推動(dòng)作用的影響,使得上游和支流的物質(zhì)在支流交匯處釋放,造成支流交匯處冰體較厚;② 可能是由于冰川的搬運(yùn)作用,使得冰磧物在冰川下游拐角處堆積,一定程度上阻礙了冰體的運(yùn)動(dòng)和消融,促進(jìn)物質(zhì)積累。

    錫亞琴冰川由于冰下水壓力的增加導(dǎo)致冰川末端基底滑動(dòng)增強(qiáng)與冰形變的加速[16],致使上游積累區(qū)物質(zhì)向下游推進(jìn),出現(xiàn)冰川下游相對(duì)上游較厚的現(xiàn)象。

    巴爾托洛冰川為典型的樹(shù)狀型冰川,主干呈東西走向,冰川整體坡度較緩,厚度自東向西逐漸減薄,與高程變化一致,冰川下游海拔較低、氣溫較高等因素促進(jìn)冰雪消融并抑制冰雪積累[17]。

    3.3 冰儲(chǔ)量估算與統(tǒng)計(jì)分析

    為進(jìn)一步估算冰川儲(chǔ)量,將冰川厚度柵格影像重采樣至50 m分辨率,與冰川邊界矢量疊置分析,得到完整的冰川區(qū),利用RGI 6.0全球冰川編目數(shù)據(jù),通過(guò)統(tǒng)計(jì)冰川區(qū)的柵格像元數(shù)目,將柵格像元大小與每個(gè)像元厚度值的乘積之和作為冰川區(qū)冰儲(chǔ)量。卻哥隆瑪冰川、錫亞琴冰川和巴爾托洛冰川統(tǒng)計(jì)面積分別為283.08、1 053.88、781.06 km2,平均厚度值分別為237.27、229.66、172.39 m,卻哥隆瑪冰川、錫亞琴冰川、巴爾托洛冰川的冰儲(chǔ)量估算結(jié)果分別為67.42、242.18、149.72 km3。

    為分析上述結(jié)果的可靠性,采用文獻(xiàn)[10]中的法國(guó)南部比利牛斯天文臺(tái)數(shù)據(jù)服務(wù)中心的冰川厚度產(chǎn)品作為驗(yàn)證數(shù)據(jù),該產(chǎn)品經(jīng)過(guò)精度檢驗(yàn),已被應(yīng)用于全球冰川流速與厚度特征分析等研究中。本文反演結(jié)果與厚度產(chǎn)品對(duì)比見(jiàn)表2所列。

    表2 3條冰川厚度、冰儲(chǔ)量本文結(jié)果與厚度產(chǎn)品對(duì)比

    從表2可以看出,卻哥隆瑪冰川、錫亞琴冰川、巴爾托洛冰川的平均厚度與厚度產(chǎn)品的誤差分別為14.05%、13.86%、11.31%,而冰儲(chǔ)量的誤差分別為14.05%、14.30%、13.85%,表明本文冰厚反演方法具有較好的可行性,可用于輔助分析青藏高原冰川厚度時(shí)空特征差異。

    4 討 論

    本文采用層流理論進(jìn)行冰川厚度反演,該方法對(duì)于光學(xué)影像和雷達(dá)影像的流速結(jié)果都具有良好的適用性,但模型中不同參數(shù)的組合易導(dǎo)致“異參同效”現(xiàn)象的發(fā)生。本研究根據(jù)冰川表面流速、坡度、形狀因子等模型參數(shù)的不確定性,考慮冰川不同區(qū)域的氣象條件差異,結(jié)合不同山谷冰川所處氣候類型、外界應(yīng)力及橫斷面結(jié)構(gòu)等屬性的差異,設(shè)置合理的參數(shù),盡可能減少“異參同效”現(xiàn)象,使得該方法獲得較為理想的結(jié)果。

    模型反演存在不確定性,其中輸入數(shù)據(jù)、模型參數(shù)是造成反演結(jié)果不確定性的主要來(lái)源。

    1) 根據(jù)式(6)可得,若冰川表面流速增加10倍,則會(huì)對(duì)冰厚高估超過(guò)75%[10],會(huì)導(dǎo)致冰厚模擬值在冰川的活躍階段與靜止階段之間不一致,這進(jìn)一步增加了厚度值的不確定性。

    2) 對(duì)于山地冰川,可變參數(shù)τ和形狀因子f′的取值分別與冰川谷壁橫向阻力和冰川橫截面的形狀有關(guān),但實(shí)測(cè)值與理論值往往有很大差異[18]。在區(qū)域尺度上,由于冰川類型、氣候、地形等因素的顯著差異,且統(tǒng)一特征參數(shù)難以確定以及缺少有關(guān)物質(zhì)平衡的區(qū)域性數(shù)據(jù),該方法在區(qū)域研究中具有較大的局限性。

    3) 雖然本文中冰川流速有具體的時(shí)間信息,但由于DEM數(shù)據(jù)的匱乏,無(wú)法獲取目標(biāo)時(shí)間段的DEM數(shù)據(jù),采用了2003年公開(kāi)發(fā)布的SRTM DEM數(shù)據(jù)進(jìn)行坡度提取,時(shí)間的不一致最終也會(huì)直接導(dǎo)致冰厚反演結(jié)果的誤差[19]。

    目前對(duì)于冰體流動(dòng)的相關(guān)研究中,關(guān)于冰體在給定應(yīng)力作用下的變形規(guī)律尚認(rèn)識(shí)不清,且冰體的變形速率與所處溫度和外界施加的壓力密切相關(guān)。雖然眾多研究者設(shè)計(jì)了盡可能接近冰川的冰厚實(shí)驗(yàn)來(lái)建立冰厚模型,但是實(shí)驗(yàn)條件均為理想狀態(tài),相對(duì)于自然條件存在一定差別。

    影響冰川厚度變化的因素有很多,降水和溫度在氣象因子中占有重要地位[20]??龅貐^(qū)降水較為充沛,但3條冰川所在的區(qū)域年均降水量有較大差異,卻哥隆瑪冰川年均降水量約為1 000 mm,而錫亞琴冰川、巴爾托洛冰川年均降水量分別約為600、700 mm,并且卻哥隆瑪冰川、錫亞琴冰川、巴爾托洛冰川年平均氣溫分別為-13.79、-22.93、-20.61 ℃,降水量和氣溫的差異會(huì)影響冰川消融與積累的轉(zhuǎn)化效率,進(jìn)而影響冰川的補(bǔ)給與冰儲(chǔ)量的變化,增加了冰厚反演結(jié)果的不確定性。

    5 結(jié) 論

    本文利用6景Sentinel-1A衛(wèi)星影像和90 m 空間分辨率SRTM DEM數(shù)據(jù),應(yīng)用層流理論模型對(duì)卻哥隆瑪冰川、錫亞琴冰川、巴爾托洛冰川3條冰川進(jìn)行厚度反演及冰儲(chǔ)量估算,分析冰厚空間分布主要特點(diǎn)及其原因,結(jié)論如下:

    1) 基于冰川表面流速和層流理論能夠反演復(fù)雜地形區(qū)山地冰川的厚度信息,通過(guò)與法國(guó)南部比利牛斯天文臺(tái)數(shù)據(jù)服務(wù)中心的冰川厚度數(shù)據(jù)產(chǎn)品對(duì)比驗(yàn)證了結(jié)果的有效性。本文方法可為青藏高原地區(qū)山地冰川厚度反演與冰儲(chǔ)量估算研究提供參考。

    2) 卻哥隆瑪冰川、錫亞琴冰川、巴爾托洛冰川的平均厚度分別為237.27、229.66、172.39 m,同一地區(qū)不同冰川厚度略有差異,冰川的搬運(yùn)作用和堆積作用與冰厚的空間分布存在著必然的聯(lián)系,但冰川整體上呈從上游到下游逐漸減薄、支流向干流不斷增厚的趨勢(shì)。

    3) 層流理論物理意義明確,相對(duì)于傳統(tǒng)經(jīng)驗(yàn)公式更具有穩(wěn)健性,基于層流理論進(jìn)行冰厚估算時(shí),由于模型過(guò)于簡(jiǎn)化,估算結(jié)果誤差偏大[15],實(shí)測(cè)資料的稀缺一定程度上限制了該方法結(jié)果的驗(yàn)證,如何有效利用豐富的冰川編目資料及其他輔助數(shù)據(jù)改進(jìn)厚度估算模型是以后進(jìn)一步研究的方向。

    猜你喜歡
    冰體冰川反演
    彈體高速侵徹冰體研究
    反演對(duì)稱變換在解決平面幾何問(wèn)題中的應(yīng)用
    高速?gòu)楏w侵徹冰材料過(guò)程數(shù)值模擬研究
    為什么冰川會(huì)到處走?
    冰體質(zhì)量和撞擊角度對(duì)船首結(jié)構(gòu)碰撞性能的影響
    冰川會(huì)發(fā)出聲音嗎?
    基于船-水-冰耦合技術(shù)的撞擊參數(shù)對(duì)船冰碰撞性能的影響
    長(zhǎng)途跋涉到冰川
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一级黄色大片毛片| 自拍偷自拍亚洲精品老妇| 午夜老司机福利剧场| 欧美一区二区精品小视频在线| 午夜视频国产福利| 亚洲人成网站高清观看| 国内精品宾馆在线| 国产成人a∨麻豆精品| 欧美日韩一区二区视频在线观看视频在线 | 欧美激情在线99| 国产不卡一卡二| 波野结衣二区三区在线| 久久久久久久久久久丰满| 国产片特级美女逼逼视频| 一级毛片久久久久久久久女| 国内揄拍国产精品人妻在线| 熟妇人妻久久中文字幕3abv| 国产男人的电影天堂91| 日韩一本色道免费dvd| 免费av观看视频| 欧美另类亚洲清纯唯美| 亚洲五月天丁香| 国产精品免费一区二区三区在线| 亚洲国产高清在线一区二区三| 久久久久久久亚洲中文字幕| 亚洲欧美清纯卡通| 久久精品影院6| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av涩爱 | 久久久久国产网址| 国产久久久一区二区三区| 能在线免费看毛片的网站| 少妇被粗大猛烈的视频| 成人特级黄色片久久久久久久| 亚洲av二区三区四区| 日本在线视频免费播放| 午夜激情福利司机影院| 亚洲aⅴ乱码一区二区在线播放| 男女边吃奶边做爰视频| 美女国产视频在线观看| 久久国产乱子免费精品| 国产亚洲精品久久久久久毛片| av在线亚洲专区| 又爽又黄a免费视频| 99热全是精品| 免费搜索国产男女视频| 九九热线精品视视频播放| 欧美bdsm另类| av在线蜜桃| 国产精品.久久久| 老女人水多毛片| 国产美女午夜福利| 97热精品久久久久久| 国产一区二区在线av高清观看| 成人毛片60女人毛片免费| 大香蕉久久网| 亚洲经典国产精华液单| 亚洲三级黄色毛片| 国产av一区在线观看免费| 国产真实乱freesex| 欧美zozozo另类| 欧美日韩乱码在线| 亚洲无线在线观看| 亚洲成人久久性| 国内精品美女久久久久久| 国产伦在线观看视频一区| 精品久久久久久久久久免费视频| 国产精品久久久久久久久免| 亚洲四区av| 国产黄色小视频在线观看| 国产一区亚洲一区在线观看| 成人三级黄色视频| 免费观看在线日韩| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产精品国产精品| 国产精品美女特级片免费视频播放器| 国产黄片视频在线免费观看| 99热网站在线观看| 国产极品天堂在线| 国产精品一区二区性色av| 国产欧美日韩精品一区二区| 久久久欧美国产精品| 尾随美女入室| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕免费在线视频6| 午夜精品国产一区二区电影 | 天堂中文最新版在线下载 | 久久综合国产亚洲精品| 亚洲欧洲国产日韩| 国产午夜福利久久久久久| 少妇的逼水好多| av在线老鸭窝| 国产亚洲av片在线观看秒播厂 | 国产精品电影一区二区三区| 国产精品人妻久久久久久| 亚洲国产精品合色在线| 麻豆久久精品国产亚洲av| 亚洲精品自拍成人| avwww免费| 亚洲美女视频黄频| 草草在线视频免费看| 久久人人爽人人爽人人片va| 五月玫瑰六月丁香| АⅤ资源中文在线天堂| 国产黄a三级三级三级人| 国产成人福利小说| 丰满乱子伦码专区| 成人特级av手机在线观看| 日韩强制内射视频| 日韩av在线大香蕉| 校园人妻丝袜中文字幕| 亚洲精品久久久久久婷婷小说 | 国产精品久久视频播放| 国产精品人妻久久久影院| 中文亚洲av片在线观看爽| 亚洲欧美成人精品一区二区| 国产精品久久久久久久久免| 日韩成人av中文字幕在线观看| 欧美一区二区国产精品久久精品| 欧美三级亚洲精品| 日本一本二区三区精品| 中文字幕免费在线视频6| 久久6这里有精品| www.av在线官网国产| 简卡轻食公司| 丝袜美腿在线中文| 精华霜和精华液先用哪个| 网址你懂的国产日韩在线| 亚洲av中文字字幕乱码综合| 国产美女午夜福利| 国产精品三级大全| 日产精品乱码卡一卡2卡三| 久久亚洲国产成人精品v| 免费av观看视频| 成人永久免费在线观看视频| 中文字幕av在线有码专区| 99精品在免费线老司机午夜| 看黄色毛片网站| 伦精品一区二区三区| av福利片在线观看| 国产极品精品免费视频能看的| 身体一侧抽搐| 欧美+亚洲+日韩+国产| 国产真实伦视频高清在线观看| 五月伊人婷婷丁香| 91麻豆精品激情在线观看国产| 亚洲精品色激情综合| 久久久久久久午夜电影| 激情 狠狠 欧美| 狂野欧美激情性xxxx在线观看| 欧美性猛交╳xxx乱大交人| 久久午夜福利片| 精品无人区乱码1区二区| 亚洲不卡免费看| 赤兔流量卡办理| 波野结衣二区三区在线| 中文字幕av成人在线电影| 26uuu在线亚洲综合色| 色哟哟·www| 欧美成人一区二区免费高清观看| 免费看日本二区| 国产成人a区在线观看| 免费看光身美女| 国产老妇女一区| 12—13女人毛片做爰片一| 深夜a级毛片| 高清日韩中文字幕在线| 18禁在线播放成人免费| 好男人视频免费观看在线| 一进一出抽搐动态| 国内久久婷婷六月综合欲色啪| 观看免费一级毛片| 观看免费一级毛片| 麻豆一二三区av精品| 国内精品一区二区在线观看| 你懂的网址亚洲精品在线观看 | 久久人妻av系列| 日本一本二区三区精品| 欧美最黄视频在线播放免费| 精品久久久久久成人av| 超碰av人人做人人爽久久| 精品午夜福利在线看| 99久久精品热视频| 久久亚洲精品不卡| 日韩成人伦理影院| 久久精品国产亚洲av涩爱 | 日韩欧美精品v在线| 99热只有精品国产| 国产精品日韩av在线免费观看| 亚洲av中文av极速乱| 亚洲精品自拍成人| 久久精品久久久久久噜噜老黄 | 成人三级黄色视频| 亚洲精品日韩在线中文字幕 | 天天躁日日操中文字幕| 精品人妻熟女av久视频| 成人二区视频| 精品久久久久久久久久久久久| 久久这里有精品视频免费| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久久久毛片| 男人舔女人下体高潮全视频| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲91精品色在线| 免费观看人在逋| 久久久国产成人免费| 老司机影院成人| 简卡轻食公司| 亚洲av不卡在线观看| 中文字幕免费在线视频6| 男女那种视频在线观看| 国产激情偷乱视频一区二区| 国产免费男女视频| 国产精品久久久久久av不卡| 中国国产av一级| 18禁裸乳无遮挡免费网站照片| 国产精华一区二区三区| 欧美极品一区二区三区四区| 老师上课跳d突然被开到最大视频| 在线国产一区二区在线| 国内精品美女久久久久久| 日韩中字成人| 亚洲av一区综合| 精品无人区乱码1区二区| 久久精品国产亚洲av天美| 中文字幕人妻熟人妻熟丝袜美| 国产成人午夜福利电影在线观看| 内射极品少妇av片p| 熟女电影av网| av天堂在线播放| 国内揄拍国产精品人妻在线| 亚洲美女视频黄频| 日韩av在线大香蕉| 一进一出抽搐动态| 午夜激情欧美在线| 欧美极品一区二区三区四区| 国产精品麻豆人妻色哟哟久久 | 精华霜和精华液先用哪个| 热99在线观看视频| 欧美性感艳星| 日本欧美国产在线视频| 亚洲欧美中文字幕日韩二区| 国产精品蜜桃在线观看 | 欧美区成人在线视频| 草草在线视频免费看| 听说在线观看完整版免费高清| 日本爱情动作片www.在线观看| 久久久久网色| 日本撒尿小便嘘嘘汇集6| 又黄又爽又刺激的免费视频.| 欧美变态另类bdsm刘玥| 看免费成人av毛片| 熟妇人妻久久中文字幕3abv| 亚洲国产高清在线一区二区三| 婷婷精品国产亚洲av| 免费av不卡在线播放| 国产亚洲精品久久久久久毛片| 久久热精品热| 婷婷六月久久综合丁香| 天堂√8在线中文| 欧美在线一区亚洲| 国产视频内射| 高清毛片免费看| 日韩成人伦理影院| 一区二区三区高清视频在线| 免费看a级黄色片| 久久精品影院6| 狠狠狠狠99中文字幕| 能在线免费看毛片的网站| 久久精品国产99精品国产亚洲性色| 国内久久婷婷六月综合欲色啪| 成人永久免费在线观看视频| 女同久久另类99精品国产91| 22中文网久久字幕| 男女下面进入的视频免费午夜| 狂野欧美激情性xxxx在线观看| 国产精品久久久久久精品电影小说 | 精品人妻一区二区三区麻豆| 啦啦啦韩国在线观看视频| 欧美日韩国产亚洲二区| 国产亚洲5aaaaa淫片| 99热只有精品国产| а√天堂www在线а√下载| 午夜激情欧美在线| 国产精品蜜桃在线观看 | 又粗又爽又猛毛片免费看| 国产国拍精品亚洲av在线观看| 性插视频无遮挡在线免费观看| 国产亚洲精品久久久com| 直男gayav资源| 欧美+日韩+精品| 1000部很黄的大片| 少妇丰满av| 一区二区三区四区激情视频 | 岛国毛片在线播放| 国产 一区精品| 99热只有精品国产| 春色校园在线视频观看| 美女高潮的动态| 99热精品在线国产| 寂寞人妻少妇视频99o| 日韩三级伦理在线观看| 亚洲欧美日韩卡通动漫| 亚洲国产色片| 最近最新中文字幕大全电影3| 亚洲国产日韩欧美精品在线观看| 在线观看av片永久免费下载| 国产精品日韩av在线免费观看| 一边摸一边抽搐一进一小说| 日韩欧美国产在线观看| 国产伦理片在线播放av一区 | 美女黄网站色视频| 免费观看人在逋| 秋霞在线观看毛片| 亚洲在线观看片| 国产精品久久久久久亚洲av鲁大| 日本免费a在线| 老司机影院成人| av在线老鸭窝| 一本一本综合久久| 国内精品久久久久精免费| 一进一出抽搐动态| 国产又黄又爽又无遮挡在线| 久久精品91蜜桃| 婷婷六月久久综合丁香| 国产精品嫩草影院av在线观看| 青春草国产在线视频 | 又爽又黄a免费视频| 寂寞人妻少妇视频99o| 久久精品久久久久久噜噜老黄 | 国产一区二区在线观看日韩| 中文精品一卡2卡3卡4更新| 卡戴珊不雅视频在线播放| 国产黄色视频一区二区在线观看 | 亚洲国产日韩欧美精品在线观看| 亚洲精品日韩在线中文字幕 | 老女人水多毛片| 99视频精品全部免费 在线| 久久精品国产亚洲av香蕉五月| 成人午夜高清在线视频| 欧美另类亚洲清纯唯美| 岛国毛片在线播放| 久久国产乱子免费精品| 日本av手机在线免费观看| 中文字幕免费在线视频6| 亚洲国产精品合色在线| 3wmmmm亚洲av在线观看| 亚洲av第一区精品v没综合| 精品日产1卡2卡| 两个人的视频大全免费| 亚洲国产精品成人综合色| 国产在线男女| 天天躁夜夜躁狠狠久久av| 97在线视频观看| 一进一出抽搐动态| 天堂网av新在线| 精品一区二区三区人妻视频| 亚洲av成人av| 精品欧美国产一区二区三| 国产精品麻豆人妻色哟哟久久 | 久久6这里有精品| 色尼玛亚洲综合影院| 看十八女毛片水多多多| 九九久久精品国产亚洲av麻豆| 18禁在线无遮挡免费观看视频| 免费一级毛片在线播放高清视频| 成人亚洲欧美一区二区av| 丰满人妻一区二区三区视频av| 久久久久九九精品影院| 免费观看精品视频网站| 亚洲真实伦在线观看| 亚洲av熟女| 黄色欧美视频在线观看| 此物有八面人人有两片| 日本av手机在线免费观看| 亚洲第一电影网av| 成人特级黄色片久久久久久久| 床上黄色一级片| 少妇熟女欧美另类| 成人午夜精彩视频在线观看| 日韩av不卡免费在线播放| 日本免费a在线| 欧美极品一区二区三区四区| 黄色日韩在线| 久久99热这里只有精品18| 精品久久久久久久末码| 九九爱精品视频在线观看| 一个人看的www免费观看视频| 草草在线视频免费看| 天天躁日日操中文字幕| 精品久久久噜噜| 人妻夜夜爽99麻豆av| 亚洲内射少妇av| 狂野欧美白嫩少妇大欣赏| 欧美日韩国产亚洲二区| 亚洲成人精品中文字幕电影| 最近中文字幕高清免费大全6| 成人二区视频| 婷婷精品国产亚洲av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 看免费成人av毛片| eeuss影院久久| 99久久无色码亚洲精品果冻| 日韩欧美国产在线观看| ponron亚洲| 成人漫画全彩无遮挡| 亚洲精品456在线播放app| 亚洲人成网站在线播放欧美日韩| 国产精品一二三区在线看| 亚洲最大成人中文| 午夜a级毛片| 久99久视频精品免费| 波多野结衣巨乳人妻| 成人三级黄色视频| 国产精华一区二区三区| 伦精品一区二区三区| 日本色播在线视频| 欧美xxxx黑人xx丫x性爽| 国产伦精品一区二区三区四那| 亚洲乱码一区二区免费版| 亚洲五月天丁香| 波多野结衣巨乳人妻| 高清日韩中文字幕在线| 非洲黑人性xxxx精品又粗又长| 在线观看免费视频日本深夜| 久久久精品94久久精品| 蜜桃久久精品国产亚洲av| 国产爱豆传媒在线观看| 欧美性猛交黑人性爽| 国产精品久久久久久久久免| 日本成人三级电影网站| a级毛片a级免费在线| 内地一区二区视频在线| 亚洲中文字幕日韩| 女人十人毛片免费观看3o分钟| 欧美成人免费av一区二区三区| 国产亚洲精品久久久久久毛片| 色综合站精品国产| 99久国产av精品| av福利片在线观看| 变态另类丝袜制服| 亚洲精品456在线播放app| 永久网站在线| 日日干狠狠操夜夜爽| 国产成人精品一,二区 | 99久久无色码亚洲精品果冻| 久久久久久久亚洲中文字幕| 日韩欧美 国产精品| 日韩强制内射视频| 欧洲精品卡2卡3卡4卡5卡区| 久久久久国产网址| 亚洲自拍偷在线| 欧美在线一区亚洲| 麻豆国产97在线/欧美| 黄色一级大片看看| 成年女人永久免费观看视频| 搡老妇女老女人老熟妇| 人妻系列 视频| 国产成人一区二区在线| 亚洲国产欧美在线一区| 97超视频在线观看视频| 午夜福利在线观看吧| 少妇的逼水好多| 午夜a级毛片| 国产片特级美女逼逼视频| 国产亚洲精品av在线| 草草在线视频免费看| 麻豆国产av国片精品| 伦精品一区二区三区| 看非洲黑人一级黄片| 乱系列少妇在线播放| 看片在线看免费视频| 精品人妻视频免费看| 男女视频在线观看网站免费| 精品久久久噜噜| 国产精品一及| 亚洲av不卡在线观看| 男的添女的下面高潮视频| 最近手机中文字幕大全| .国产精品久久| 日本爱情动作片www.在线观看| 国产精品人妻久久久久久| 亚洲18禁久久av| 久久精品国产亚洲av天美| 日韩中字成人| 亚洲熟妇中文字幕五十中出| 欧美xxxx性猛交bbbb| 寂寞人妻少妇视频99o| 日韩强制内射视频| 特大巨黑吊av在线直播| 精品人妻偷拍中文字幕| 国产一区二区亚洲精品在线观看| 女的被弄到高潮叫床怎么办| 亚洲国产精品合色在线| 亚洲av成人av| 精品日产1卡2卡| av.在线天堂| 欧美性猛交╳xxx乱大交人| 晚上一个人看的免费电影| 成人毛片60女人毛片免费| 少妇裸体淫交视频免费看高清| 亚洲天堂国产精品一区在线| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩卡通动漫| 成人综合一区亚洲| 亚洲成人精品中文字幕电影| 一本久久中文字幕| 麻豆成人av视频| 久久婷婷人人爽人人干人人爱| 麻豆国产av国片精品| av在线亚洲专区| 国产成年人精品一区二区| 欧美成人免费av一区二区三区| ponron亚洲| 性色avwww在线观看| 免费无遮挡裸体视频| or卡值多少钱| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩在线中文字幕 | 婷婷色av中文字幕| 午夜福利在线观看免费完整高清在 | 好男人视频免费观看在线| 国产精品人妻久久久久久| 亚洲人成网站在线播放欧美日韩| 亚洲成人精品中文字幕电影| 一级毛片久久久久久久久女| 99热这里只有是精品在线观看| 少妇的逼水好多| 波多野结衣巨乳人妻| eeuss影院久久| 51国产日韩欧美| 亚洲三级黄色毛片| 蜜臀久久99精品久久宅男| 日本免费一区二区三区高清不卡| 啦啦啦啦在线视频资源| 亚洲国产日韩欧美精品在线观看| 色噜噜av男人的天堂激情| 一个人观看的视频www高清免费观看| 99热精品在线国产| 亚洲欧美成人综合另类久久久 | 国产精品一区二区三区四区久久| 美女高潮的动态| 久久久久免费精品人妻一区二区| 成人毛片a级毛片在线播放| 一级毛片久久久久久久久女| 日韩国内少妇激情av| 久久久久久伊人网av| 日本黄色片子视频| 国产精品久久久久久亚洲av鲁大| 午夜精品在线福利| 乱码一卡2卡4卡精品| 联通29元200g的流量卡| 亚洲美女视频黄频| 特大巨黑吊av在线直播| 精华霜和精华液先用哪个| 桃色一区二区三区在线观看| 国产一级毛片七仙女欲春2| 最近视频中文字幕2019在线8| 青春草亚洲视频在线观看| 日韩欧美国产在线观看| 久久精品国产鲁丝片午夜精品| 国产老妇女一区| 日韩av在线大香蕉| 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 干丝袜人妻中文字幕| 边亲边吃奶的免费视频| 中国美女看黄片| 国产成人freesex在线| 久久久精品欧美日韩精品| 亚洲三级黄色毛片| 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 波多野结衣高清无吗| av卡一久久| 黄色配什么色好看| 一级毛片久久久久久久久女| 91狼人影院| 一级毛片电影观看 | 99热只有精品国产| 亚洲久久久久久中文字幕| 哪个播放器可以免费观看大片| 免费看日本二区| 亚洲国产精品sss在线观看| 国产日本99.免费观看| 村上凉子中文字幕在线| 韩国av在线不卡| 国产三级中文精品| 秋霞在线观看毛片| 午夜福利高清视频| 日韩亚洲欧美综合| av.在线天堂| 亚洲精品国产av成人精品| av天堂中文字幕网| 亚洲五月天丁香| 少妇高潮的动态图| 国产精品久久久久久av不卡| 亚洲色图av天堂| 嫩草影院新地址| 18+在线观看网站| 一级毛片久久久久久久久女| 国内精品一区二区在线观看| 好男人视频免费观看在线| 熟女电影av网| 久久九九热精品免费| 熟妇人妻久久中文字幕3abv| 偷拍熟女少妇极品色| 国产免费男女视频| 岛国在线免费视频观看| 人妻制服诱惑在线中文字幕| 日日摸夜夜添夜夜爱| 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 伦理电影大哥的女人| 国产成人影院久久av|