北極 Svalbard群島 Austre Lovénbreen和Pedersenbreen冰川表面運(yùn)動(dòng)特征

李鵬 閆明 艾松濤 徐躍通 劉雷保 孫維君

(1中國(guó)極地研究中心,國(guó)家海洋局極地科學(xué)重點(diǎn)實(shí)驗(yàn)室,上海200136;2山東師范大學(xué)人口資源與環(huán)境學(xué)院,山東濟(jì)南250014;3武漢大學(xué)中國(guó)南極測(cè)繪研究中心,湖北武漢430079;4天津市環(huán)境監(jiān)測(cè)中心,天津300191)

0 引言

冰川是氣候的產(chǎn)物,其存在與變化不僅對(duì)氣候變化具有明顯的反饋、調(diào)節(jié)和指示作用,而且對(duì)生態(tài)、水資源變化具有重要影響[1-4]。冰川運(yùn)動(dòng)是冰川區(qū)別于其它自然冰體的主要標(biāo)志之一[5-6],是冰川學(xué)研究的重要內(nèi)容[7-11]。冰川運(yùn)動(dòng)速度與氣候條件、冰川形態(tài)、物質(zhì)平衡、冰下地形、冰川融水、冰內(nèi)應(yīng)力等密切相關(guān)[10,12],不僅是冰川動(dòng)力學(xué)模型和冰川水文模型必不可少的參數(shù),也是進(jìn)行冰川相關(guān)災(zāi)害模擬和預(yù)測(cè)的重要參數(shù)[13-15]。

北極Svalbard群島約59%的地表被冰川(帽)覆蓋,冰川總數(shù)超過2 100條,普遍發(fā)育山谷或冰斗冰川[16]。該群島地處挪威暖流海熱傳輸通道的最北端,冰川對(duì)其波動(dòng)及相應(yīng)的氣候變化十分敏感,是國(guó)際上冰川監(jiān)測(cè)研究的重點(diǎn)區(qū)域[17-25]。2004年7—8月,中國(guó)北極黃河站科學(xué)考察隊(duì)首次在Svalbard群島新奧爾松地區(qū)開展了多學(xué)科聯(lián)合考察,通過對(duì)冰川進(jìn)行實(shí)地勘測(cè)和調(diào)查比較,選定 Austre Lovénbreen和Pedersenbreen兩條冰川作為長(zhǎng)期監(jiān)測(cè)和研究的對(duì)象,開展冰川物質(zhì)平衡、冰川運(yùn)動(dòng)及其對(duì)氣候變化響應(yīng)的研究[8,26]。

1 研究區(qū)概況

Austre Lovénbreen和Pedersenbreen冰川位于北極Svalbard群島新奧爾松地區(qū),距北極黃河站的直線距離分別為6.2和10 km(圖1)。兩條冰川彼此相鄰,均為多溫型山谷冰川,周圍山峰的最高海拔分別為880和1 017 m[8]。依據(jù)挪威極地所2008年出版Svalbard群島大比例尺地形圖A7量算,Austre Lovénbreen 冰川幾何中心在東經(jīng) 12°9′44.4″,北緯78°52′17.3″,面積為 5.69 km2,海拔分布范圍為 75—600 m,Pedersenbreen 冰川幾何中心在12°17′10.1″E,78°51′25.0″N,面積為 6.59 km2,海拔分布范圍為100—700 m。

圖1 Austre Lovénbreen和Pedersenbreen冰川位置示意圖Fig.1.A sketch map showing the location of galciers Austre Lovénbreen and Pedersenbreen

2 數(shù)據(jù)處理

北 極 Svalbard群 島 Austre Lovénbreen和Pedersenbreen冰川表面運(yùn)動(dòng)監(jiān)測(cè)始于2005年7月,在Austre Lovénbreen冰川由A至E斷面布設(shè)16根花桿,在Pedersenbreen冰川沿冰川主流線布設(shè)5根花桿,2007年7月在Austre Lovénbreen冰川增設(shè)F點(diǎn)花桿。2005—2008年,每年7—8月測(cè)量花桿位置以獲得冰川表面年度運(yùn)動(dòng)速度;2009—2011年的每年4—5月和9月分別測(cè)量花桿位置以獲得冰川表面季節(jié)性運(yùn)動(dòng)速度。2005—2008年的每年7—8月份、2009—2011年的每年9月份獲取的是兩條冰川花桿與冰面(很少情況下是雪面)交匯處高程信息,2009—2011年的每年4—5月份獲取的是兩條冰川花桿與雪面交匯處高程信息。由于冰川消融等因素影響,部分年度會(huì)有花桿倒伏且找不到原鉆孔位置,導(dǎo)致冰川運(yùn)動(dòng)數(shù)據(jù)缺失(表1)。本文統(tǒng)計(jì)2005年7月至2011年9月Austre Lovénbreen和Pedersenbreen冰川表面運(yùn)動(dòng)速度,其中 Austre Lovénbreen冰川各觀測(cè)點(diǎn)共有69個(gè)年度數(shù)據(jù),35個(gè)夏季數(shù)據(jù)和32個(gè)冬季數(shù)據(jù);Pedersenbreen冰川各觀測(cè)點(diǎn)共有18個(gè)年度數(shù)據(jù),5個(gè)夏季數(shù)據(jù)和7個(gè)冬季數(shù)據(jù)。

Austre Lovénbreen和Pedersenbreen冰川到黃河站GPS跟蹤站的距離均在10 km以內(nèi),可視為短距離差分計(jì)算,建立地方獨(dú)立坐標(biāo)系統(tǒng)[27-29]:以黃河站GPS跟蹤站為坐標(biāo)原點(diǎn),其正東方為X軸、正北方為Y軸,Z軸垂直于XY平面。每根花桿位置采用Leica 1230 GPS接收機(jī)靜態(tài)測(cè)量50 min以上,測(cè)量精度達(dá)到水平方向優(yōu)于1 cm、高程方向優(yōu)于2 cm的水平。GPS數(shù)據(jù)獲取及處理[27]和冰川表面運(yùn)動(dòng)速度計(jì)算方法[29]詳見文獻(xiàn)[27,29]。本文規(guī)定一個(gè)年度數(shù)據(jù)時(shí)間跨度開始于夏季或夏季末測(cè)量時(shí),截止到來年相同月份。冰川表面水平或垂直運(yùn)動(dòng)速度平均值為統(tǒng)計(jì)時(shí)段內(nèi)所有觀測(cè)點(diǎn)年度平均運(yùn)動(dòng)速度值的平均值,冰川表面夏季(或冬季)水平或者垂直運(yùn)動(dòng)速度平均值為統(tǒng)計(jì)時(shí)段內(nèi)所有觀測(cè)點(diǎn)相應(yīng)季節(jié)內(nèi)水平或者垂直運(yùn)動(dòng)速度值的平均值。冰川表面應(yīng)變率的計(jì)算如下:假設(shè)冰川表面有兩個(gè)點(diǎn)A2和B2,在t年A2到B2點(diǎn)的水平距離為A2B2(t),在t+1年為 A2B2(t+1),則該時(shí)段內(nèi) A2到 B2的應(yīng)變率ε(A2B2,t)為[30]:

表1 2005—2011年冰川表面運(yùn)動(dòng)數(shù)據(jù)缺失情況統(tǒng)計(jì)表Table 1.The data missing points of glacial surface movement from 2005 to 2011

3 結(jié)果與討論

3.1 水平運(yùn)動(dòng)速度特征

Austre Lovénbreen冰川表面水平運(yùn)動(dòng)速度平均為2.14 m·a-1,最大值5.17 m·a-1位于2007/2008年D3,最小值0.27 m·a-1位于2007/2008年A1。相比之下,Pedersenbreen冰川表面水平運(yùn)動(dòng)速度較大,平均為6.28 m·a-1,最大值8.63 m·a-1位于2010/2011年 P4,最小值3.82 m·a-1位于2009/2010年P(guān)2。與新奧爾松地區(qū)的潮水冰川Kronebreen(其冰川下游表面運(yùn)動(dòng)速度大約在600 m·a-1[22])相比,這兩條山谷冰川表面水平運(yùn)動(dòng)速度很小。如圖2所示,Austre Lovénbreen冰川2007/2008年水平運(yùn)動(dòng)速度略高于其他年度,2010/2011年水平運(yùn)動(dòng)速度最小,其余四個(gè)年度水平運(yùn)動(dòng)速度相近,其曲線相互交織在一起。Pedersenbreen冰川2005/2006年水平運(yùn)動(dòng)速度略高于其他年度,2008/2009年水平運(yùn)動(dòng)速度最小。相鄰兩條冰川表面水平運(yùn)動(dòng)速度相差較大的原因可能在于冰川形態(tài)(包括表面坡度)、基底地形和冰下融水狀況等方面的不同,而不同年度之間冰川表面運(yùn)動(dòng)速度的差異則可能與氣溫、物質(zhì)平衡和冰下融水狀況有關(guān)。

圖2 冰川表面水平運(yùn)動(dòng)速度分布Fig.2.The distribution of horizontal velocities on glacial surface

Austre Lovénbreen冰川表面各觀測(cè)點(diǎn)水平運(yùn)動(dòng)速度的變差系數(shù)平均為0.24,最大0.44位于D4,最小0.03位于C2;Pedersenbreen冰川的變差系數(shù)平均為0.14,最大0.28位于 P2,最小0.09位于P5。Austre Lovénbreen冰川表面水平運(yùn)動(dòng)速度變化大于Pedersenbreen冰川,原因可能在于Austre Lovénbreen冰川面積和冰儲(chǔ)量相對(duì)較小——Pedersenbreen冰川面積、冰儲(chǔ)量分別比 Austre Lovénbreen冰川大15.8%和18.6%[31-32]。

從2009年4月至2011年9月,共積累兩條冰川表面三期夏季水平運(yùn)動(dòng)速度和兩期冬季水平運(yùn)動(dòng)速度(圖3)。Austre Lovénbreen冰川表面夏季水平運(yùn)動(dòng)速度平均為0.24 m·mon-1(米/月),最大值0.38 m·mon-1位于2009年夏D2和2011年夏D1,最小值0.02 m·mon-1位于2009年夏D5;冬季水平運(yùn)動(dòng)速度平均為0.22 m·mon-1,最大值0.36 m·mon-1位于2010年冬 D1,最小值 0.02 m·mon-1位于2010年冬A1。Pedersenbreen冰川表面夏季水平運(yùn)動(dòng)速度平均為0.69 m·mon-1,最大值0.82m·mon-1位于2009年夏P4,最小值0.50 m·mon-1位于2009年夏P5;冬季水平運(yùn)動(dòng)速度平均為0.55 m·mon-1,最大值0.67 m·mon-1位于2010年冬P4,最小值0.45 m·mon-1位于2010年冬P1。兩條冰川表面夏季水平運(yùn)動(dòng)速度均略高于冬季值并且變化也強(qiáng)于冬季。Austre Lovénbreen冰川表面夏季水平運(yùn)動(dòng)速度比冬季快9.1%,其極差是冬季值的1.77倍,Pedersenbreen冰川表面夏季水平運(yùn)動(dòng)速度比冬季快25.5%,其極差是冬季值的2.17倍。

圖3 冰川表面夏季和冬季水平運(yùn)動(dòng)速度分布Fig.3.The distribution of horizontal velocities in summer and in winter on glacial surface

冰川表面水平運(yùn)動(dòng)速度分布與海拔高度有較好的相關(guān)性(圖4)。Austre Lovénbreen冰川主流線表面水平運(yùn)動(dòng)速度[y/(m·a-1)]與海拔高度(x/hm)有關(guān)系式:

其相關(guān)系數(shù)rAL為0.9372,滿足顯著性水平0.10。Pedersenbreen冰川主流線表面水平運(yùn)動(dòng)速度[y/(m·a-1)]與海拔高度(x/hm)有關(guān)系式:

其相關(guān)系數(shù)rPD0.854 0,滿足顯著性水平0.20。

Austre Lovénbreen冰川主流線表面水平運(yùn)動(dòng)速度表現(xiàn)為隨海拔高度升高先增加后減少再增加的趨勢(shì),應(yīng)與冰厚、冰下地形和冰面形態(tài)有關(guān)系:雷達(dá)測(cè)厚數(shù)據(jù)顯示[31]E-F斷面冰厚大幅減小,但冰下地形陡升;且E-F斷面在冰面形態(tài)上表現(xiàn)為巨大的陡坡。與Austre Lovénbreen冰川不同,Pedersenbreen冰川主流線表面水平運(yùn)動(dòng)速度表現(xiàn)為隨海拔高度升高先增加后減少的趨勢(shì),雷達(dá)測(cè)厚數(shù)據(jù)顯示其底部為發(fā)育較為成熟的U型谷,且冰下地形抬升較為平緩[32]。

圖4 冰川主流線表面水平運(yùn)動(dòng)速度與海拔高度關(guān)系Fig.4.The relationship between horizontal velocity on glacial surface and its elevation along themainsteam line

Austre Lovénbreen冰川表面各觀測(cè)點(diǎn)平均水平運(yùn)動(dòng)速度矢量空間分布(圖5)表現(xiàn)出:冰川主流線水平運(yùn)動(dòng)速度高于兩側(cè),并且冰川兩側(cè)水平運(yùn)動(dòng)速度不對(duì)稱,C斷面及其以下主流線東側(cè)水平運(yùn)動(dòng)速度高于西側(cè),而在D斷面主流線東側(cè)水平運(yùn)動(dòng)速度明顯較低。

圖5 冰川表面各觀測(cè)點(diǎn)平均水平和垂直運(yùn)動(dòng)速度矢量圖Fig.5.A vectorgraph showing the distribution of average horizontal and vertical velocities on glacial surface

將分布在冰川中央的花桿稱為主流線花桿,其平均水平運(yùn)動(dòng)速度稱為主流線速度;將分布在Austre Lovénbreen冰川測(cè)量斷面上各花桿的平均水平運(yùn)動(dòng)速度再求平均值,稱其為斷面速度。兩條冰川表面均以(在冰川縱剖面上,冰川主流線為一曲線,用L表示,并假設(shè)原點(diǎn)在冰川源頭,以Uxy表示冰川表面水平運(yùn)動(dòng)速度)為界形成了兩種具有不同運(yùn)動(dòng)狀態(tài)的運(yùn)動(dòng)區(qū):為運(yùn)動(dòng)拉伸區(qū)和＜0為運(yùn)動(dòng)壓縮區(qū)。Austre Lovénbreen冰川從上游至下游依次表現(xiàn)為壓縮、拉伸和壓縮,而Pedersenbreen冰川從上游至下游依次表現(xiàn)為拉伸和壓縮,兩者有所不同。如圖6所示,Austre Lovénbreen冰川主流線速度與斷面速度分布不盡一致,表現(xiàn)為D斷面速度明顯小于C、E斷面,使Austre Lovénbreen冰川E斷面到D斷面表現(xiàn)為運(yùn)動(dòng)壓縮區(qū),原因在于D斷面含蓋Austre Lovénbreen冰川兩個(gè)支流,并且兩個(gè)支流表面水平運(yùn)動(dòng)速度差異懸殊,導(dǎo)致斷面速度值變小。在Austre Lovénbreen冰川 D斷面,D2、D3平均水平速度高于E斷面且低于C斷面。

圖6 冰川表面主流線速度和斷面速度Fig.6.Meanmainstream line and sectional velocities on glacial surface of Austre Lovénbreen and Pedersenbreen

3.2 垂直運(yùn)動(dòng)速度特征

Austre Lovénbreen冰川表面垂直運(yùn)動(dòng)速度平均為0.76 m·a-1,運(yùn)動(dòng)方向垂直向下,最大值1.82 m·a-1位于2008/2009年D1,最小值0.03 m·a-1位于2005/2006年E1。Pedersenbreen冰川表面垂直運(yùn)動(dòng)速度平均為0.90 m·a-1,運(yùn)動(dòng)方向垂直向下,最大值1.77 m·a-1位于2008/2009年P(guān)3,最小值0.01 m·a-1位于2007/2008年P(guān)5。如圖7所示,兩條冰川均表現(xiàn)出2008/2009年冰川表面垂直運(yùn)動(dòng)速度最大,2010/2011年次之,2007/2008年冰川表面垂直運(yùn)動(dòng)速度最小,2005/2006年、2006/2007年和2009/2010年與六年平均冰川表面垂直運(yùn)動(dòng)速度很接近。冰川表面各觀測(cè)點(diǎn)的年際變化趨勢(shì)基本一致,但各觀測(cè)點(diǎn)極差不盡相同,Austre Lovénbreen冰川A斷面觀測(cè)點(diǎn)極差明顯較小。兩條冰川表面垂直運(yùn)動(dòng)速度空間分布最突出的特征是其沿冰川主流線方向逐漸增加(圖5),并且在Austre Lovénbreen冰川兩側(cè)分布不對(duì)稱,東支冰川表面垂直運(yùn)動(dòng)速度明顯小于西支。

圖7 冰川表面垂直運(yùn)動(dòng)速度分布Fig.7.The distribution of vertical velocities on glacial surface

Austre Lovénbreen冰川表面各觀測(cè)點(diǎn)垂直運(yùn)動(dòng)速度的變差系數(shù)平均為0.81,最大1.69位于F,最小0.06位于A3;E2和E1變差系數(shù)也較大,分別為1.61和1.59。Pedersenbreen冰川的變差系數(shù)平均為0.59,最大0.89位于P5,最小0.35位于P3。

Austre Lovénbreen和Pedersenbreen冰川表面夏季垂直運(yùn)動(dòng)速度平均分別為0.14 m·mon-1和0.13 m·mon-1,運(yùn)動(dòng)方向均為垂直向下;冬季垂直運(yùn)動(dòng)速度平均分別為0.13 m·mon-1和0.09 m·mon-1,運(yùn)動(dòng)方向均為垂直向上。兩條冰川表面夏季垂直運(yùn)動(dòng)速度分別高于冬季值7.69%和44.44%,但夏季垂直運(yùn)動(dòng)速度變差系數(shù)分別低于冬季值23.53%和46.67%。

Austre Lovénbreen冰川主流線表面垂直運(yùn)動(dòng)速度[y/(m·a-1)]與海拔高度(x/hm)有關(guān)系式:

其相關(guān)系數(shù)rAL為0.962 7,滿足顯著性水平0.01。Pedersenbreen冰川主流線表面垂直運(yùn)動(dòng)速度[y/(m·a-1)]與海拔高度(x/hm)有關(guān)系式:

其相關(guān)系數(shù)為rPD為0.924 8,滿足顯著性水平0.05。

冰川表面垂直運(yùn)動(dòng)速度(Uz)有三個(gè)方面的來源[33-34]:一是水平運(yùn)動(dòng)速度在垂直方向上的分量Uxy×tga,a為冰川表面坡度,簡(jiǎn)稱為水平分量(H);二為冰川表面物質(zhì)平衡引起的高程變化,簡(jiǎn)稱為物變高程(ΔB);三是冰川冰體內(nèi)部的補(bǔ)償流(V),其運(yùn)動(dòng)方向?yàn)榇怪毕蛏匣蛳蛳隆Ｋ鼈冎g有關(guān)系式:

對(duì)于Austre Lovénbreen冰川表面垂直運(yùn)動(dòng)速度,物變高程占有主導(dǎo)地位,其貢獻(xiàn)率為64%,其次為補(bǔ)償流,貢獻(xiàn)率為19%,最后是水平分量,貢獻(xiàn)率為16%(表2)。對(duì)于Pedersenbreen表面垂直運(yùn)動(dòng)速度,物變高程貢獻(xiàn)率最大,其值為40%,其次為水平分量,貢獻(xiàn)率為35%,最后為補(bǔ)償流,貢獻(xiàn)率為25%。

3.3 表面應(yīng)變率分析

Austre Lovénbreen和Pedersenbreen冰川沿主流線的表面應(yīng)變率ε可分為三個(gè)區(qū)域(圖8):(1)ε＞0,分布在冰川上游,為冰川受張力拉伸區(qū),E2F段表面應(yīng)變率平均為0.001 5,D3E2段為0.000 2,P4P5段為0.003 0,冰川上游表面應(yīng)變率具有沿著主流線向下減少的趨勢(shì),表明冰體所受的張力在減少;(2)ε=0,分布在冰川中游,冰川冰既不受擠壓力也不受張力的影響;(3)ε＜0,分布在冰川下游地區(qū),為冰川受擠壓力的壓縮區(qū),Austre Lovénbreen冰川A2D3段表面應(yīng)變率平均為-0.001 5,其中A2B2段為-0.002 2,B2C2段為-0.001 7,C2D3段為-0.000 6,Pedersenbreen冰川P1P4段表面應(yīng)變率平均為-0.001 4,其中P1P2段為-0.002 3,P2P3段為-0.001 8,P3P4段為-0.000 2,冰川下游表面應(yīng)變率表現(xiàn)出沿著主流線向下逐漸增加的趨勢(shì),表明冰體所受的擠壓力在增強(qiáng)。

表2 冰川表面垂直運(yùn)動(dòng)速度分量所占比例表Table 2.The percentage of each component of vertical velocity on glacial surface

圖8 沿冰川主流線表面應(yīng)變率分布Fig.8.The distribution of surface strain rate along glacialmainstream line

Austre Lovénbreen和Pedersenbreen冰川各段表面應(yīng)變率的變差系數(shù)平均分別為0.19和0.15(表3),其變化量可達(dá)到該段表面應(yīng)變率平均值的15%—19%。表面應(yīng)變率的變差系數(shù)表現(xiàn)出沿冰川主流線向下先增加后減少再增加的趨勢(shì);表面應(yīng)變率極差表現(xiàn)為沿冰川主流線向下先減少后增加的趨勢(shì)。

表3 冰川主流線表面應(yīng)變率變化統(tǒng)計(jì)表Table 3.The variation statistics of surface strain rate along galcialmainstream line

4 結(jié)論

本文分析北極Svalbard地區(qū)Austre Lovénbreen和Pedersenbreen冰川表面運(yùn)動(dòng)特征,主要結(jié)論如下:

(1)Austre Lovénbreen和 Pedersenbreen冰川表面水平運(yùn)動(dòng)速度平均分別為2.14 m·a-1和6.28 m·a-1,變差系數(shù)平均分別為0.24和0.14。兩條冰川表面夏季水平運(yùn)動(dòng)速度均略高于冬季并且變化也強(qiáng)于冬季,Pedersenbreen冰川表面水平運(yùn)動(dòng)速度的季節(jié)性變化比Austre Lovénbreen更大。冰川主流線表面水平運(yùn)動(dòng)速度與所處海拔高度具有多項(xiàng)式型相關(guān)性。冰川橫斷面上,Austre Lovénbreen冰川主流線表面水平運(yùn)動(dòng)速度最大,兩側(cè)速度小且不對(duì)稱;冰川縱剖面上,Austre Lovénbreen冰川從源頭至末端表現(xiàn)為運(yùn)動(dòng)的壓縮區(qū)、拉伸區(qū)和壓縮區(qū),而Pedersenbreen冰川從源頭至末端表現(xiàn)為運(yùn)動(dòng)的拉伸區(qū)和壓縮區(qū),兩者有所不同。

(2)Austre Lovénbreen和 Pedersenbreen冰川表面垂直運(yùn)動(dòng)速度平均分別為0.76 m·a-1和0.90 m·a-1,運(yùn)動(dòng)方向均為垂直向下,變差系數(shù)平均分別為0.81和0.59。冰川表面垂直運(yùn)動(dòng)速度年際變化較大且變化趨勢(shì)基本一致。兩條冰川表面夏季垂直運(yùn)動(dòng)速度均大于冬季,且夏季變差系數(shù)比冬季小。冰川主流線表面垂直運(yùn)動(dòng)速度與所處海拔高度具有一元線性相關(guān)性。對(duì) Austre Lovénbreen和Pedersenbreen冰川表面垂直運(yùn)動(dòng)速度構(gòu)成貢獻(xiàn)率最大的因素分別依次為物變高程(64%)、補(bǔ)償流(19%)、水平分量(16%)和物變高程(40%)、水平分量(35%)、補(bǔ)償流(25%)。

(3)Austre Lovénbreen和 Pedersenbreen冰川由上游至下游,表面應(yīng)變率表現(xiàn)為從大于0逐步減小到0,然后再負(fù)向增加,表明冰川上游為受張力拉伸區(qū),逐步過渡到冰川下游受擠壓力的壓縮區(qū),且冰川下游冰體所受的擠壓力沿著主流線向下逐漸增強(qiáng)。兩條冰川各段表面應(yīng)變率的變差系數(shù)平均分別為0.19和0.15,其變化量可達(dá)該段表面應(yīng)變率平均值的15%—19%。表面應(yīng)變率變差系數(shù)表現(xiàn)出沿冰川主流線向下先增加后減少再增加的趨勢(shì);其極差表現(xiàn)為沿冰川主流線向下先減少后增加的趨勢(shì)。

致謝本文所使用Austre Lovénbreen 和Pedersenbreen冰川運(yùn)動(dòng)及物質(zhì)平衡數(shù)據(jù)是北極黃河站冰川考察隊(duì)全體隊(duì)員艱辛工作的結(jié)果,謹(jǐn)此對(duì)參加該工作的所有人員深表謝意。

1 Dyurgerov M B,Meier M F.Twentieth century climate change:Evidence from small glaciers.Proceedings of the National Academy of Sciences of the United States of America,2000,97(4):1406—1411.

2 Roe G H.What do glaciers tell us about climate variability and climate change.Journal of Glaciology,2011,57(203):567—578.

3 任賈文,葉柏生,丁永建,等.中國(guó)冰凍圈變化對(duì)海平面上升潛在貢獻(xiàn)的初步估計(jì).科學(xué)通報(bào),2011,56(14):1084—1087.

4 姚檀棟,李治國(guó),楊威,等.雅魯藏布江流域冰川分布和物質(zhì)平衡特征及其對(duì)湖泊的影響.科學(xué)通報(bào),2010,55(18):1750—1756.

5 曹敏,李忠勤,李慧林.天山托木爾峰地區(qū)青冰灘72號(hào)冰川表面運(yùn)動(dòng)速度特征研究.冰川凍土,2011,33(1):21—29.

6 黃茂桓,孫作哲,劉宗香.巴托拉冰川末端的運(yùn)動(dòng)學(xué)研究.科學(xué)通報(bào),1981,26(8):494—496.

7 Zwally H J,Li J,Brenner A C,etal.Greenland ice sheetmass balance:distribution of increasedmass losswith climatewarming;2003—2007 versus1992—2002.Journal of Glaciology,2011,57(201):88—102.

8 閆明,任賈文,張占海,等.斯瓦爾巴群島冰川學(xué)研究進(jìn)展與我國(guó)北極冰川監(jiān)測(cè)系統(tǒng)建設(shè).極地研究,2006,18(2):137—147.

9 黃茂桓,孫作哲.我國(guó)大陸型冰川運(yùn)動(dòng)的某些特征.冰川凍土,1982,4(2):35—45.

10 井哲帆,劉力,周在明,等.冰川運(yùn)動(dòng)速度影響因子的強(qiáng)度分析——以祁連山七一冰川為例.冰川凍土,2011,33(6):1222—1228.

11 井哲帆,葉柏生,焦克勤,等.天山奎屯河哈希勒根51號(hào)冰川表面運(yùn)動(dòng)特征分析.冰川凍土,2002,24(5):563—566.

12 Rignot E,Kanagaratnam P.Changes in the velocity structure of the Greenland Ice Sheet.Science,2006,311(5763):986—990.

13 李慧林,李忠勤,秦大河.冰川動(dòng)力學(xué)模式基本原理和參數(shù)觀測(cè)指南.北京:氣象出版社,2009:27—32.

14 卿文武,陳仁升,劉時(shí)銀.冰川水文模型研究進(jìn)展.水科學(xué)進(jìn)展,2008,19(6):893—902.

15 Mcnabb RW,Hock R,O’Neel S,et al.Using surface velocities to calculate ice thickness and bed topography:a case study at Columbia Glacier,Alaska,USA.Journal of Glaciology,2012,58(212):1151—1164.

16 Hagen JO,Liest?l O,Roland E,et al.Glacier Atlas of Svalbard and Jan Mayen.Oslo:Meddelelser,NR.129,1993:12—35.

17 Mansell D,Luckman A,Murray T.Dynamics of tidewater surge-type glaciers in northwest Svalbard.Journal of Glaciology,2012,58(207):110—118.

18 Nuth C,Schuler T V,Kohler J,etal.Estimating the long-term calving flux of Kronebreen,Svalbard,from geodetic elevation changes andmass-balancemodelling,Journal of Glaciology,2012,58(207):119—133.

19 Nuth C,Moholdt G,Kohler J,et al.Svalbard glacier elevation changes and contribution to sea level rise.Journal of Geophysical Research,2010,115(F1):F01008,doi:10.1029/2008JF001223.

20 Rye C J,Willis IC,Arnold N S,et al.On the need for automated multiobjective optimization and uncertainty estimation of glaciermass balance models.Journal of Geophysical Research,2012,117(F2):F02005,doi:10.1029/2011JF002184.

21 Sund M,Eiken T,Hagen JO,et al.Svalbard surge dynamics derived from geometric changes.Annals of Glaciology,2009,50(52):50—60.

22 K??b A,Lefauconnier B,Melvold K.Flow field of Kronebreen,Svalbard,using repeated Landsat7 and ASTER data.Annals of Glaciology,2005,42(1):7—13.

23 Hagen JO,Kohler J,Melvold K,etal.Glaciers in Svalbard:massbalance,runoff and freshwater flux.Polar Research,2003,22(2):145—159.

24 Rippin D,Willis I,Arnold N,etal.Changes in geometry and subglacial drainage ofMidre Lovénbreen,Svalbard,determined from digital elevation models.Earth Surface Processes and Landforms,2003,28(3):273—298.

25 Svendsen H,Beszczynska-M?ller A,Hagen JO,et al.The physical environment of Kongsfjorden-Krossfjorden,an Arctic fjord system in Svalbard.Polar Research,2002,21(1):133—166.

26 任賈文,閆明.中國(guó)北極黃河站首次科學(xué)考察隊(duì)冰川考察.冰川凍土,2005,27(1):124—127.

27 艾松濤,鄂棟臣,閆明,等.2005年北極冰川首期GPS監(jiān)測(cè).極地研究,2006,18(1):1—8.

28 李鵬,閆明,徐躍通,等.基于GIS冰川末端變化計(jì)算方法研究——以北極Austre Lovénbreen冰川為例,冰川凍土.2012,34(2):367—374.

29 徐明星,閆明,任賈文,等.北極Svalbard地區(qū)Austre Lovénbreen和Pedersenbreen冰川表面物質(zhì)平衡和運(yùn)動(dòng)特征分析.極地研究,2010,22(1):10—22.

30 施雅風(fēng),黃茂桓,任炳輝等.中國(guó)冰川概論.北京:科學(xué)出版社,1988:88—104.

31 Saintenoy A,Friedt J-M,Booth A D,etal.Deriving ice thickness,glacier volume and bedrockmorphology of the Austre Lovénbreen(Svalbard)using Ground-penetrating Radar.Near Surface Geophysics,2013,11(2):253—261.

32 Ai ST,Wang ZM,ED C,et al.Topography,ice thickness and ice volume of the glacier Pedersenbreen in Svalbard,using GPR and GPS.Polar Research,2014,33:18533,doi:10.3402/polar.v33.18533.

33 Cuffey K M,Paterson W SB.The Physics of Glaciers.4th ed.Oxford:Butterworth-Heinemann,2010:285—398.

34 謝自楚,劉潮海.冰川學(xué)導(dǎo)論.上海:上海科學(xué)普及出版社,2010:238—243.