• <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一維反演
    国产老妇女一区| 好看av亚洲va欧美ⅴa在| 脱女人内裤的视频| 18+在线观看网站| 欧美绝顶高潮抽搐喷水| 在线观看免费视频日本深夜| 悠悠久久av| 国产熟女xx| 少妇被粗大猛烈的视频| 色噜噜av男人的天堂激情| 热99在线观看视频| 精品久久久久久久久av| 欧美精品国产亚洲| 国产精品爽爽va在线观看网站| 国产黄色小视频在线观看| 久久性视频一级片| 国产精品一区二区免费欧美| 亚洲成a人片在线一区二区| 美女 人体艺术 gogo| 免费观看人在逋| 免费人成视频x8x8入口观看| 青草久久国产| 国内精品久久久久久久电影| 亚洲久久久久久中文字幕| 日本一本二区三区精品| 色综合婷婷激情| 黄色一级大片看看| 国产极品精品免费视频能看的| 亚洲精品亚洲一区二区| 日韩人妻高清精品专区| 国产精品一区二区三区四区免费观看 | 夜夜爽天天搞| 十八禁人妻一区二区| 午夜精品在线福利| 国产中年淑女户外野战色| 一卡2卡三卡四卡精品乱码亚洲| 天堂影院成人在线观看| 国产色婷婷99| 好看av亚洲va欧美ⅴa在| 久久国产精品影院| 亚洲精品亚洲一区二区| 亚洲av美国av| 99精品久久久久人妻精品| 国产成人a区在线观看| 国产精品久久久久久久电影| 免费电影在线观看免费观看| 国产真实乱freesex| 精品一区二区三区av网在线观看| 日韩人妻高清精品专区| 亚洲av一区综合| 美女高潮的动态| 亚洲成av人片在线播放无| 波多野结衣高清无吗| a级毛片免费高清观看在线播放| 国产成人aa在线观看| 啦啦啦韩国在线观看视频| 日韩中文字幕欧美一区二区| 国产高清视频在线观看网站| 亚洲欧美清纯卡通| 超碰av人人做人人爽久久| 在线观看免费视频日本深夜| 国产精品1区2区在线观看.| 久久久久国内视频| 亚洲熟妇中文字幕五十中出| 美女黄网站色视频| 日韩欧美免费精品| 亚洲五月天丁香| 99热6这里只有精品| 国产成人啪精品午夜网站| 成人亚洲精品av一区二区| 99久久精品国产亚洲精品| 欧美激情在线99| 中国美女看黄片| а√天堂www在线а√下载| 色综合站精品国产| 久久亚洲真实| 他把我摸到了高潮在线观看| 亚洲精品粉嫩美女一区| 搡老妇女老女人老熟妇| 日日干狠狠操夜夜爽| 日本免费a在线| 在现免费观看毛片| 岛国在线免费视频观看| 国产亚洲精品久久久com| 国产欧美日韩一区二区精品| 一级黄片播放器| 国产成人aa在线观看| 欧美成人一区二区免费高清观看| 成人av在线播放网站| 99国产综合亚洲精品| 性插视频无遮挡在线免费观看| 国产黄a三级三级三级人| 亚洲内射少妇av| 中国美女看黄片| 亚洲精品久久国产高清桃花| 国产主播在线观看一区二区| 亚洲av电影在线进入| 啪啪无遮挡十八禁网站| av女优亚洲男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 免费av观看视频| 亚洲中文字幕一区二区三区有码在线看| 精品人妻视频免费看| 长腿黑丝高跟| 欧美xxxx黑人xx丫x性爽| 桃色一区二区三区在线观看| 亚洲欧美日韩高清在线视频| 99久久九九国产精品国产免费| 亚洲 国产 在线| 久久精品久久久久久噜噜老黄 | 丰满的人妻完整版| 午夜激情福利司机影院| 人妻制服诱惑在线中文字幕| 久久亚洲精品不卡| 久久久国产成人免费| 午夜免费男女啪啪视频观看 | 国内少妇人妻偷人精品xxx网站| 欧美一区二区国产精品久久精品| 一级av片app| 一区福利在线观看| 俺也久久电影网| 成人美女网站在线观看视频| 午夜福利在线在线| 亚洲黑人精品在线| 国产av麻豆久久久久久久| 日本五十路高清| 嫩草影院新地址| 国产精品永久免费网站| 波多野结衣高清无吗| 男女之事视频高清在线观看| 久久国产精品人妻蜜桃| 午夜福利18| 亚洲无线在线观看| 淫秽高清视频在线观看| 久久久久国产精品人妻aⅴ院| 亚洲第一区二区三区不卡| 亚洲精品在线观看二区| 男女床上黄色一级片免费看| 99久久无色码亚洲精品果冻| 亚洲一区高清亚洲精品| 可以在线观看毛片的网站| 精品不卡国产一区二区三区| 亚洲精品成人久久久久久| 午夜福利免费观看在线| 毛片女人毛片| 丰满人妻熟妇乱又伦精品不卡| 国产精品乱码一区二三区的特点| 亚洲 国产 在线| 一二三四社区在线视频社区8| 97热精品久久久久久| 成人性生交大片免费视频hd| 又黄又爽又免费观看的视频| 国产精品,欧美在线| 别揉我奶头~嗯~啊~动态视频| 国产精品不卡视频一区二区 | 又爽又黄无遮挡网站| av天堂中文字幕网| 亚洲精品粉嫩美女一区| 久久伊人香网站| 日本黄大片高清| 我要看日韩黄色一级片| 又黄又爽又刺激的免费视频.| 老熟妇仑乱视频hdxx| 日本免费a在线| 精品久久久久久久久久久久久| 日本黄大片高清| 国产精品精品国产色婷婷| 亚洲男人的天堂狠狠| 欧美日韩瑟瑟在线播放| av国产免费在线观看| 国产精品久久视频播放| 日日摸夜夜添夜夜添小说| 美女免费视频网站| 亚洲国产精品成人综合色| 亚洲一区高清亚洲精品| 黄色日韩在线| 精品一区二区免费观看| 欧美日韩福利视频一区二区| 看十八女毛片水多多多| 中文字幕熟女人妻在线| 九色国产91popny在线| 91午夜精品亚洲一区二区三区 | 久久草成人影院| 又紧又爽又黄一区二区| 超碰av人人做人人爽久久| 精品人妻一区二区三区麻豆 | 色综合婷婷激情| avwww免费| 欧美成狂野欧美在线观看| 亚洲成av人片在线播放无| 久久婷婷人人爽人人干人人爱| 麻豆国产av国片精品| 一本综合久久免费| 国产高清视频在线观看网站| 亚洲国产精品久久男人天堂| 国产av在哪里看| 毛片一级片免费看久久久久 | 久久久久久久久中文| 一a级毛片在线观看| 成年人黄色毛片网站| 成人永久免费在线观看视频| 国产高清视频在线观看网站| 91麻豆av在线| 91字幕亚洲| 免费人成视频x8x8入口观看| 欧美日本视频| 日本精品一区二区三区蜜桃| 精品99又大又爽又粗少妇毛片 | АⅤ资源中文在线天堂| 国产精品永久免费网站| av天堂在线播放| 啪啪无遮挡十八禁网站| 中文资源天堂在线| 少妇的逼好多水| 亚洲欧美精品综合久久99| 国产免费一级a男人的天堂| 日本a在线网址| 淫妇啪啪啪对白视频| 精品久久久久久,| 亚洲最大成人手机在线| 网址你懂的国产日韩在线| 有码 亚洲区| 亚洲人成网站在线播| 日韩免费av在线播放| 亚洲国产精品成人综合色| 亚洲av中文字字幕乱码综合| 国产精品,欧美在线| 最新在线观看一区二区三区| 国产高清视频在线播放一区| 最后的刺客免费高清国语| 亚洲在线自拍视频| 一本综合久久免费| 在线免费观看不下载黄p国产 | 人妻久久中文字幕网| 欧美精品啪啪一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 日韩av在线大香蕉| 国内久久婷婷六月综合欲色啪| av福利片在线观看| 他把我摸到了高潮在线观看| 最后的刺客免费高清国语| 男女那种视频在线观看| 日本黄大片高清| 国产国拍精品亚洲av在线观看| 能在线免费观看的黄片| 亚洲成av人片在线播放无| 亚洲欧美日韩无卡精品| 性色av乱码一区二区三区2| 色视频www国产| 免费看光身美女| 成年女人永久免费观看视频| 又黄又爽又刺激的免费视频.| 成人三级黄色视频| 免费在线观看亚洲国产| 男女做爰动态图高潮gif福利片| 在现免费观看毛片| 国产视频一区二区在线看| 欧美日韩乱码在线| 真人做人爱边吃奶动态| 国产91精品成人一区二区三区| 亚洲成a人片在线一区二区| 国产一级毛片七仙女欲春2| 91在线精品国自产拍蜜月| 成人欧美大片| 久久久成人免费电影| 最新在线观看一区二区三区| 久久欧美精品欧美久久欧美| 男女之事视频高清在线观看| 久久精品国产亚洲av天美| 国产精品人妻久久久久久| 亚洲av不卡在线观看| 在线国产一区二区在线| 午夜福利高清视频| 国产精品三级大全| 免费在线观看成人毛片| 伊人久久精品亚洲午夜| 国产高清激情床上av| 俺也久久电影网| 免费在线观看日本一区| 久久久久亚洲av毛片大全| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 国产午夜福利久久久久久| a级毛片免费高清观看在线播放| 男人狂女人下面高潮的视频| 婷婷精品国产亚洲av| 黄片小视频在线播放| 久9热在线精品视频| 好看av亚洲va欧美ⅴa在| 欧美bdsm另类| 国产精品伦人一区二区| 国产av在哪里看| av女优亚洲男人天堂| 国产精品亚洲一级av第二区| 精品久久久久久久末码| 99久久九九国产精品国产免费| 国产精品精品国产色婷婷| 窝窝影院91人妻| 99热只有精品国产| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩av片在线观看| 最新在线观看一区二区三区| 人人妻人人澡欧美一区二区| 无人区码免费观看不卡| 91午夜精品亚洲一区二区三区 | 国产中年淑女户外野战色| 欧美日韩瑟瑟在线播放| 一区二区三区激情视频| 99国产精品一区二区三区| 日本与韩国留学比较| 99在线人妻在线中文字幕| 国产人妻一区二区三区在| 又爽又黄a免费视频| 白带黄色成豆腐渣| 亚洲无线在线观看| 亚洲av不卡在线观看| 国产极品精品免费视频能看的| 综合色av麻豆| 性欧美人与动物交配| 三级男女做爰猛烈吃奶摸视频| 国产精品一区二区三区四区久久| 18禁黄网站禁片午夜丰满| 亚洲国产精品合色在线| 国产精品,欧美在线| 免费在线观看成人毛片| 婷婷亚洲欧美| 99在线人妻在线中文字幕| 免费看日本二区| 在线观看免费视频日本深夜| 综合色av麻豆| 99久久精品一区二区三区| 亚洲av美国av| 高潮久久久久久久久久久不卡| 亚洲av.av天堂| 直男gayav资源| 久久精品夜夜夜夜夜久久蜜豆| 国产v大片淫在线免费观看| 日本五十路高清| 国产中年淑女户外野战色| 在线观看av片永久免费下载| 三级毛片av免费| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 嫩草影视91久久| 欧美日韩福利视频一区二区| 91午夜精品亚洲一区二区三区 | 别揉我奶头 嗯啊视频| 男人舔女人下体高潮全视频| 欧美另类亚洲清纯唯美| 美女免费视频网站| 极品教师在线视频| 精品久久久久久久久亚洲 | 又黄又爽又免费观看的视频| 啪啪无遮挡十八禁网站| 在线a可以看的网站| 夜夜爽天天搞| 九九在线视频观看精品| 成年人黄色毛片网站| 午夜福利在线观看免费完整高清在 | 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 日韩欧美在线乱码| 18禁在线播放成人免费| 国产精品乱码一区二三区的特点| 免费人成视频x8x8入口观看| 在线a可以看的网站| 日韩欧美精品免费久久 | 午夜激情欧美在线| 国产三级在线视频| 首页视频小说图片口味搜索| 五月伊人婷婷丁香| 国产成人啪精品午夜网站| 搡女人真爽免费视频火全软件 | 最后的刺客免费高清国语| 在线a可以看的网站| 在线观看午夜福利视频| 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 1000部很黄的大片| 精品熟女少妇八av免费久了| av天堂在线播放| bbb黄色大片| 精品福利观看| 极品教师在线免费播放| 欧美3d第一页| 亚洲成av人片免费观看| 亚洲午夜理论影院| 久久久成人免费电影| 免费在线观看亚洲国产| 亚洲人成网站在线播放欧美日韩| 嫩草影院精品99| 真实男女啪啪啪动态图| 日韩中字成人| 国产成人av教育| 久久国产精品影院| 欧美性猛交╳xxx乱大交人| 亚洲av成人精品一区久久| 国产不卡一卡二| 无人区码免费观看不卡| 亚洲 欧美 日韩 在线 免费| 天美传媒精品一区二区| 日日夜夜操网爽| 男女做爰动态图高潮gif福利片| 国产亚洲欧美98| 国产爱豆传媒在线观看| 十八禁人妻一区二区| 免费无遮挡裸体视频| 中国美女看黄片| 精品国产亚洲在线| bbb黄色大片| 国产精品女同一区二区软件 | 午夜两性在线视频| 麻豆国产av国片精品| 精品人妻偷拍中文字幕| 成年免费大片在线观看| 最新中文字幕久久久久| 男女那种视频在线观看| 丰满乱子伦码专区| 伊人久久精品亚洲午夜| 在线观看午夜福利视频| 国产一级毛片七仙女欲春2| 成人性生交大片免费视频hd| 日本与韩国留学比较| 欧美成人免费av一区二区三区| 最新中文字幕久久久久| 国产在线精品亚洲第一网站| a级毛片免费高清观看在线播放| 国产熟女xx| 免费观看人在逋| 午夜视频国产福利| 久久久国产一区二区| 亚洲精品国产色婷婷电影| 三级国产精品片| 成年女人看的毛片在线观看| 日本av手机在线免费观看| 精品亚洲乱码少妇综合久久| 国产黄频视频在线观看| 老司机影院成人| 国产高清不卡午夜福利| 国产成人免费无遮挡视频| 尾随美女入室| 久久久精品94久久精品| 欧美97在线视频| 亚洲综合色惰| 女的被弄到高潮叫床怎么办| 久久99蜜桃精品久久| 亚洲精品第二区| 国产高潮美女av| 日本欧美国产在线视频| 国产黄片美女视频| 国产av国产精品国产| 免费看光身美女| 青春草亚洲视频在线观看| 国产在线男女| 免费不卡的大黄色大毛片视频在线观看| 高清在线视频一区二区三区| 国产成人福利小说| 免费看光身美女| 国产美女午夜福利| 亚洲欧美精品专区久久| 亚洲av在线观看美女高潮| 午夜激情久久久久久久| 久久久久久久久久成人| 18禁在线无遮挡免费观看视频| 国产成人精品久久久久久| 亚洲精品第二区| 少妇被粗大猛烈的视频| 国产成人精品一,二区| 国产视频首页在线观看| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 国产精品偷伦视频观看了| 欧美成人午夜免费资源| 久久久久久久久大av| 99热这里只有精品一区| 99久久精品热视频| 成人综合一区亚洲| 久久人人爽av亚洲精品天堂 | .国产精品久久| 女人十人毛片免费观看3o分钟| 五月开心婷婷网| 高清欧美精品videossex| 最近中文字幕2019免费版| 六月丁香七月| av在线蜜桃| 久久精品久久久久久久性| 亚洲怡红院男人天堂| 男女无遮挡免费网站观看| 国产精品女同一区二区软件| 国产老妇女一区| 51国产日韩欧美| 亚洲精品日韩av片在线观看| 真实男女啪啪啪动态图| av在线天堂中文字幕| 欧美最新免费一区二区三区| 亚洲欧美中文字幕日韩二区| 婷婷色综合大香蕉| 国产成人精品一,二区| 欧美性猛交╳xxx乱大交人| 亚洲经典国产精华液单| 国产精品99久久99久久久不卡 | 成人综合一区亚洲| 综合色丁香网| 欧美日韩一区二区视频在线观看视频在线 | 国产日韩欧美在线精品| 青青草视频在线视频观看| 三级男女做爰猛烈吃奶摸视频| 国产美女午夜福利| 国内揄拍国产精品人妻在线| 日韩一区二区视频免费看| 久久99热这里只频精品6学生| 国产有黄有色有爽视频| 韩国高清视频一区二区三区| 我的女老师完整版在线观看| 国产精品不卡视频一区二区| 日韩欧美精品免费久久| 国产极品天堂在线| 亚洲av在线观看美女高潮| 欧美另类一区| 三级经典国产精品| 国产成人aa在线观看| 男女下面进入的视频免费午夜| 青春草视频在线免费观看| 69人妻影院| av在线蜜桃| 国产精品蜜桃在线观看| 国产成人精品婷婷| 亚洲最大成人中文| 男人爽女人下面视频在线观看| 免费观看a级毛片全部| 99热这里只有是精品50| 色视频在线一区二区三区| 免费观看av网站的网址| 成人鲁丝片一二三区免费| 亚洲精品乱久久久久久| 亚洲在久久综合| 毛片女人毛片| 日本午夜av视频| 男人添女人高潮全过程视频| 啦啦啦中文免费视频观看日本| 熟女电影av网| 干丝袜人妻中文字幕| 黑人高潮一二区| 精品少妇黑人巨大在线播放| 亚洲精品日本国产第一区| 亚洲精品成人久久久久久| 男女无遮挡免费网站观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲av欧美aⅴ国产| 中文字幕免费在线视频6| 男女边吃奶边做爰视频| 三级经典国产精品| 99热6这里只有精品| 亚洲av.av天堂| 亚洲性久久影院| 美女cb高潮喷水在线观看| 国产精品偷伦视频观看了| 又黄又爽又刺激的免费视频.| 午夜免费鲁丝| 91久久精品电影网| 麻豆精品久久久久久蜜桃| 精品人妻视频免费看| 天天躁日日操中文字幕| 极品少妇高潮喷水抽搐| 99久久中文字幕三级久久日本| 乱码一卡2卡4卡精品| 亚洲av.av天堂| 亚洲在线观看片| 一级毛片电影观看| 亚洲性久久影院| 丝袜美腿在线中文| 熟女电影av网| 国产精品一区www在线观看| 神马国产精品三级电影在线观看| 干丝袜人妻中文字幕| 久久久久精品久久久久真实原创| 中文在线观看免费www的网站| 又粗又硬又长又爽又黄的视频| 国产成人午夜福利电影在线观看| 超碰av人人做人人爽久久| 大香蕉久久网| 最后的刺客免费高清国语| 永久免费av网站大全| 99热这里只有是精品在线观看| 99热国产这里只有精品6| 直男gayav资源| 大码成人一级视频| 国产午夜精品一二区理论片| 一级毛片久久久久久久久女| 国产白丝娇喘喷水9色精品| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 女人被狂操c到高潮| 久久久成人免费电影| 美女被艹到高潮喷水动态| 午夜亚洲福利在线播放| 成人亚洲欧美一区二区av| 好男人在线观看高清免费视频| av在线观看视频网站免费| 一边亲一边摸免费视频| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 18禁动态无遮挡网站| 我的女老师完整版在线观看| a级一级毛片免费在线观看| 深爱激情五月婷婷| 亚洲,欧美,日韩| 亚洲高清免费不卡视频| 你懂的网址亚洲精品在线观看| 亚洲欧美一区二区三区国产| 欧美日韩视频精品一区| 欧美97在线视频| 中国国产av一级| h日本视频在线播放| 久久久久网色|