Extraction and Analysis of Temporal-Spatial Variation Characteristics of Surface Velocity of the Gangnalou Glacier
-
摘要: 冰川表面流速是反映冰川动态变化的重要指标，能够为冰川物质平衡提供重要信息。利用2016年的13景Sentinel-1A影像和合成孔径雷达（synthetic aperture radar，SAR）偏移量追踪法测定岗纳楼冰川表面流速场，并根据地表温度是否大于0 ℃将其分为冰封期（10月-次年3月）和消融期（4月-9月）。其中, 冰封期流速场较为稳定，中轴流速快于两侧，边界流速趋于0 m/a；消融期流速场冰舌中下部与冰舌上部差异明显，冰川表面流速具有分区特征。基于SAR观测数据，在冰川积累增厚、消融减薄及冰川表面流动的综合影响下，冰封期冰川表面流速最大为65.43 m/a；消融期冰川表面流速最大为46.28 m/a。通过分析冰川中线流速、冰面高程以及地表温度之间的关系可知, 冰封期内，岩床坡度是影响冰川流速的主要原因，且岗纳楼冰川为“压缩-舒张”流运动模式；消融期内，地表温度升高、冰面海拔变化是冰川表面流速出现分区现象的主要原因。Abstract:Objectives The temporal-spatial variation characteristic of the glacier surface velocity is an important indicator of dynamic change of glacier and can provide important information for glacier mass balance. Surface temperature, glacier bed slope and glacier altitude are all important factors that affect glacier surface velocity, and their respective effects need to be analyzed.Methods In this paper, the offset-tracking method was used to measure the surface velocity of the Gangnalou Glacier, based on high-resolution Sentinel-1A synthetic aperture radar (SAR) images ranging from January 2016 to December 2016. We divided the surface velocity into freezing period and ablation period based on whether surface temperature was greater than 0 ℃ or not. For further analysis, we extracted the central line glacier velocity, and then jointly used surface temperature and digital elevation model (DEM) to analyze the influences of temperature and glacier bed slope.Results The results showed that surface velocity of freezing period was more stable and central line glacier velocity was faster than that of boundary which tended to 0 m/a. SAR measurements were affected by glacier ablation, accumulation and motion, and showed that the maximum surface velocity was 65.43 m/a during freezing period and 46.28 m/a during ablation period. This surface velocity field showed a partitioning phenomenon in which the velocities of the middle and lower parts of glacier tongue were significantly different from that of the upper part during ablation period.Conclusions The surface velocity of the Gangnalou Glacier was affected by many complex factors, of which glacier ablation and glacier bed slope are main factors. The partitioning phenomenon in the ablation period could be better explained by increasing surface temperature and changing glacier altitude, which directly affected glacier ablation and accumulation. Additionally, by analyzing the effect of the glacier bed slope, we found that the Gangnalou Glacier followed the compressive and extending flow mode in the freezing period.
-
作为最为广泛使用的电离层经验模型之一，国际参考电离层模型（international reference ionosphere, IRI）对电离层研究和应用具有十分重要的意义[1]。IRI模型利用了垂测仪、非相干散射雷达、卫星观测及探空火箭数据，可以输出任意指定时间和经纬度的60~2 000 km的电子密度和总电子含量（total electron content, TEC）等多种电离层参数信息[2]。自1972年以来，随着电离层观测数据的积累、观测技术的革新和建模方法的改进，IRI模型已发布了多个版本，目前最新版本为IRI-2016 [3-4]。
峰值密度NmF2和峰值高度hmF2是表征电子密度廓线形态的重要参量。对于NmF2的计算，IRI-2016模型中包含IRI_CCIR和IRI_URSI两种模型选项。对于hmF2的计算，在旧版本模型使用的IRI_CCIR hmF2模型的基础上，IRI-2016模型中另外新增了IRI_AMTB和IRI_Shubin两种模型选项。hmF2的IRI_CCIR模型是基于1954—1958年全球分布的约150个垂测仪观测资料，利用hmF2与参数M(3 000)F2的强相关关系由M(3 000)F2计算得到。其中，M(3 000)F2为3 000 km高度处最大可用频率和F2层临界频率（foF2）的比值。这种间接计算hmF2的方法有诸多缺陷。IRI-2016模型中新增的两种计算hmF2的模型选项直接对hmF2进行建模，不依赖M(3 000)F2。其中，IRI_AMTB模型是根据全球分布的26个垂测仪（其中包括中国区域的武汉站、安阳站和三亚站）在1998—2006年的hmF2观测值，纯粹利用地基垂测仪观测值进行建模[5]；IRI_Shubin模型的建立则使用了挑战性小卫星有效载荷（Challenging Minisatellite Payload，CHAMP）、重力恢复与气候实验（Gravity Recovery and Climate Experiment，GRACE）和气象、电离层和气候卫星联合观测系统(Constellation Observing System for Meteorology, Ionosphere and Climate，COSMIC)3个掩星任务的观测资料以及分布全球的62个垂测仪于1987—2012年的观测数据[6]。
COSMIC由美国国防部和中国台湾地区联合研制并于2006-04-15发射升空，在掩星观测任务中具有里程碑的意义[7]。COSMIC由6颗低轨卫星组成，轨道高度在800 km左右，轨道倾角为72°，并配备能够接收GPS信号的接收机。掩星事件发生时，GPS信号穿过电离层时受电子密度影响而发生折射，产生附加相位延迟。在信号直线传播和球对称假设等前提下，由传播路径上的总电子含量使用Abel积分变换方法可得到电子密度廓线。COSMIC掩星反演结果已经得到了多种独立观测手段的验证[8-10]。
中国区域（70°E~140°E，15°N~55°N）跨越中低纬地区，电离层变化幅度较大。文献[11]利用中国区域的4个垂测仪台站对IRI-2016模型在中国区域的性能进行了评估，但只能对垂测仪台站周围的电离层特征进行验证分析。COSMIC掩星观测具有全球覆盖、高精度、垂直分辨率高等特点，可用于对IRI-2016模型在整个中国区域的表现进行分析。本文的研究内容包括使用2008—2014年的COSMIC掩星数据对地磁平静条件下IRI-2016模型在中国区域的性能进行验证，并对IRI-2016模型在中国区域使用时不同NmF2和hmF2模型选项的输出结果进行统计分析，以确定IRI-2016模型在中国区域输出值的最优选项。
1 数据与模型
1.1 COSMIC掩星数据
COSMIC掩星二级产品“ionPrf”数据可从COSMIC数据存档与分析中心（COSMIC Data Analysis and Archive Center, CDAAC）的网站http://cdaac-www.cosmic.ucar.edu下载得到。该数据包含了掩星事件发生的时间、位置、NmF2、hmF2和电子密度廓线等信息。在分析地磁平静条件下IRI模型的表现时，需要先剔除3 h ap指数大于15的时间段的COSMIC掩星数据[8]。另一方面，尽管CDAAC在反演电子密度廓线时采取了一定的质量检核措施，但通过其检核的产品中仍然存在问题廓线[8]。文献[12]提出根据顶部电离层电子密度廓线的斜率和电子密度波动的平均偏差（mean deviation, MD）来剔除这些问题廓线。本文选择COSMIC电子密度廓线的标准为在420~490 km的斜率小于-0.02；平均偏差不超过1.5。此外，峰值高度hmF2范围应在200~500 km之间。
1.2 IRI-2016模型
IRI-2016模型输出值可从http://irimodel.org/网站获取。该模型中计算NmF2的选项包括IRI_CCIR模型和IRI_URSI模型。其中，IRI_URSI模型为推荐选项。计算峰值高度hmF2的选项包括原有的IRI_CCIR模型以及新增的IRI_AMTB模型和IRI_Shubin模型。其中，IRI_AMTB模型为推荐选项。模型的输入参数包括经纬度和时间。由于掩星观测具有水平漂移，本文选用COSMIC掩星观测的F2层峰值高度处的经纬度和发生时刻作为输入量，通过控制模型选项得到不同选项下输出的NmF2和hmF2。图 1(a)、1(b)分别为当地时间（local time, LT）2018-02-03T13:35和2008-02-08T23:46的COSMIC掩星电子密度廓线和IRI-2016模型不同选项下输出结果的比较图。可以看出，COSMIC电子密度廓线和IRI-2016模型不同选项以及选项之间的输出结果存在着明显差异。
![]()
图 1 COSMIC反演的电子密度廓线和IRI-2016模型不同选项输出值的比较Figure 1. Comparison Among Electron Density Profiles Derived from COSMIC RO and IRI-2016 Model with Different Model Options1.3 F10.7指数
F10.7指数是指10.7 cm（2 800 MHz）波长的太阳辐射通量，单位为sfu，与太阳黑子数之间有强相关性，能很好地表征太阳活动强度。图 2是F10.7指数的月均值在第24太阳活动周2008—2014年期间的变化。可以看出，2008—2010年为太阳活动低年，尤其是2008、2009年，太阳活动达到了极低的水平；之后太阳活动水平逐渐提高，在2011—2014年太阳活动维持在较高水平，为太阳活动高年。
2 结果分析
2.1 IRI-2016模型在中国区域的性能分析
为了了解IRI-2016模型在中国区域的表现，首先采用模型推荐的选项计算NmF2和hmF2，将结果与COSMIC掩星观测进行比较。图 3为2008年和2014年COSMIC掩星观测与IRI-2016模型输出的峰值密度NmF2与峰值高度hmF2的散点图。其中，红色实线代表COSMIC观测值等于IRI预测值，黑色虚线代表使用最小二乘方法得到的拟合线，R和S分别表示两种方法得到的相关系数和拟合线斜率。由图 3(a)和3(b)可以看出，尽管在太阳活动低年2008年和太阳活动高年2014年，COSMIC掩星观测和IRI模型输出的NmF2相关系数非常接近（分别为0.86和0.87），但在2014年数据对较为均匀地围绕红色实线分布，且拟合线和等值线非常接近，斜率高达0.90。而在2008年，数据对较为分散，且较多的数据对分布在等值线下侧，拟合直线与等值线之间存在较大偏差。拟合直线的斜率（0.70）表明在太阳活动低年2008年，相较于COSMIC掩星观测值，IRI-2016模型明显高估了NmF2。与峰值密度NmF2情形类似，图 3(c)和图 3(d)表明, 在2014年IRI-2016模型峰值高度hmF2输出值与COSMIC掩星观测值更为接近，拟合直线几乎与等值线重合。而在2008年，拟合直线的斜率为0.76，远小于2014年的斜率0.97，说明2008年IRI-2016模型的hmF2输出值相对于COSMIC掩星观测存在明显高估。同时，由相关系数可以看出，无论在2008年还是2014年，NmF2的相关系数均高于hmF2，表明IRI-2016模型和COSMIC掩星观测输出的NmF2一致性优于hmF2。
![]()
图 3 COSMIC掩星观测与IRI-2016模型采用推荐选项输出的NmF2与hmF2散点图Figure 3. Scatter Plots of NmF2 and hmF2 from COSMIC RO and IRI-2016 Model with Recommended Model Options为了进一步定量分析太阳活动低年和高年COSMIC掩星观测值和IRI-2016模型输出值的差异，定义了相对中误差绝对值的均值（mean absolute relative error，MARE）M，计算公式为：
$$ M=\frac{\mathop{\sum }^{}\left( \left| I-C \right|/C \right)}{N} $$ 式中，C表示由掩星数据反演的NmF2或者hmF2；I表示对应的由模型计算的NmF2或者hmF2；N为观测数据个数。
表 1给出了2008—2014年的NmF2和hmF2的MARE数值。可以看出，无论是NmF2还是hmF2，太阳活动低年（2008—2010年）的MARE均高于太阳活动高年（2011—2014年），在太阳活动极小年（2008—2009年），IRI-2016模型输出值与COSMIC掩星观测值的差异最大。
表 1 2008—2014年NmF2和hmF2的MARE/%Table 1. MARE of NmF2 and hmF2 During 2008―2014/%统计项 2008年 2009年 2010年 2011年 2012年 2013年 2014年 NmF2 31.112 2 31.481 8 27.836 7 26.7882 21.888 1 24.443 1 25.204 4 hmF2 11.689 0 11.944 2 9.747 9 8.460 8 7.071 9 8.021 6 7.344 6 因此，总体而言，对中国区域，无论是NmF2还是hmF2，以COSMIC掩星数据为参照，IRI-2016模型在太阳活动高年的表现都优于太阳活动低年。
2.2 太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布
为了进一步分析中国区域COSMIC掩星反演的峰值参数与IRI-2016模型输出值间的差异，对比了差异更为明显的太阳活动低年（2008—2010年）峰值参数的空间分布。考虑到电离层活跃水平，分析的时间段选为最具代表性的春秋季的中午时分（LT 12:00—14:00）。在本文中，季节的划分标准为：3月、4月、9月、10月为春秋季，5月—8月为夏季，11月、12月、1月和2月为冬季。图 4为2008—2010年春秋季中午（LT 12:00—14:00）由COSMIC掩星与IRI-2016模型得到的峰值参数NmF2与hmF2的均值在中国区域的空间分布。图 4(a)~4(d)显示，COSMIC掩星观测与IRI模型输出的峰值参数空间分布形态均反映出低纬地区电离层更为活跃，同时观测值与模型输出值之间存在明显差异。图 4(e)、4(g) 表明，对于NmF2，与COSMIC掩星反演结果对比，IRI模型输出值在15°N~30°N的低纬区域明显高估，相对偏差最高可达100%以上；而在30°N~55°N的中纬区域，IRI模型输出的NmF2普遍表现为低估。图 4(f)、4(h)表明，对于hmF2，情况恰恰相反：在30°N~55°N中纬区域，IRI模型的hmF2输出值普遍高于COSMIC观测值，而在15°N~30°N低纬区域，模型输出值普遍低于COSMIC观测值。由图 4（f）可见，在太阳活动低年，IRI-2016模型与COSMIC在中国区域得到的hmF2的绝对偏差一般不超过50 km。
![]()
图 4 2008—2010年春秋季LT 12：00—14：00 COSMIC掩星与IRI-2016模型的NmF2与hmF2在中国区域的空间分布Figure 4. Spatial Distributions of NmF2 and hmF2 over China Region at LT 12:00—14:00 During Equinox Months of 2008—2010 Derived from COSMIC RO and IRI-2016 Model2.3 IRI-2016模型峰值参数不同选项输出结果的比较分析
以上分析中所使用的IRI-2016模型输出值均基于模型的推荐选项得到。除了推荐选项，IRI-2016对NmF2和hmF2的计算均提供其他模型选项，而使用不同选项的模型输出结果存在差异。为了了解不同选项的模型在中国区域的表现，进一步以COSMIC掩星数据为参照，对IRI-2016模型中两种计算NmF2的选项（IRI_CCIR和IRI_URSI模型选项）和3种计算hmF2的选项（IRI_CCIR、IRI_AMTB和IRI_Shubin模型选项）的输出结果进行了比较。图 5给出了不同季节COSMIC掩星与IRI-2016模型两种不同选项得到的NmF2均值的日变化情况，其中图 5(a)~5(c)和图 5(d)~5(f)分别为太阳活动低年和高年的结果。统计中如果掩星事件发生时间与某一整点时刻相差不超过30 min，则该时刻的掩星数据和IRI模型输出结果被用于计算该整点时刻的平均值。可以看出，无论采用哪个选项，NmF2在约LT 04：00达到最小值，之后快速升高，在约LT 14：00达到全天最大值，之后逐渐下降。春秋季NmF2高于其他季节，电离层最为活跃。基于两种不同选项的IRI-2016模型输出值与COSMIC观测值均有相似的日变化趋势，但分别基于两种选项的模型输出结果之间存在一定差异。图 5(a)~5(c)表明采用两种选项的IRI-2016模型输出结果非常接近，但在中午时分，IRI_CCIR模型和IRI_URSI模型均高估了峰值密度，其中IRI_URSI的高估程度更严重。在春秋季和冬季，IRI_CCIR模型和IRI_URSI模型在LT 17：00—05：00则低估了峰值密度，其中IRI_URSI低估程度更严重。图 5(d)~5(f)与图 5(a)~5(c)的对比表明，随着太阳活动水平的增加，NmF2随之增大，但中午时分采用两种选项的IRI-2016模型输出值与COSMIC观测值之间的差异均有减小。特别是在太阳活动高年的夏季，IRI-2016模型和COSMIC观测值全天都符合得很好。而在春秋季LT 16：00—04：00，两种选项的IRI-2016模型输出值均小于COSMIC观测值，其中IRI_URSI的偏离度相对更大。
![]()
图 5 不同季节COSMIC掩星与IRI-2016模型两种不同选项下NmF2的日变化图Figure 5. Diurnal Variations of NmF2 from COSMIC RO and IRI-2016 Model Using Two Different Options for Different Seasons图 6给出了不同季节COSMIC掩星观测与IRI-2016模型3种不同选项输出的hmF2日变化情况。可以看出，3种选项对应的IRI-2016模型输出与COSMIC观测的hmF2呈现相似的日变化趋势，即在LT 06：00―07：00达到最小值，在LT 12：00达到小高峰，并在约LT 24：00LT达到全天最大值。在太阳活动低年，3种计算hmF2的模型选项中，IRI_Shubin模型输出值和COSMIC观测值最为接近，其他两种模型输出值相对于COSMIC观测值均存在高估。IRI_AMTB模型输出值相对于COSMIC观测值的高估现象更为突出，冬季最显著。与太阳活动低年相比较，在太阳活动高年，3种选项对应的hmF2模型输出值与COSMIC观测值的一致性均相对更好。但IRI_AMTB模型输出值在夜间的偏差仍较大。
![]()
图 6 不同季节COSMIC掩星与IRI-2016模型不同选项下hmF2的日变化图Figure 6. Diurnal Variations of hmF2 from COSMIC RO and IRI-2016 Model Using Different Options for Different Seasons3 讨论
本文以COSMIC掩星数据为参照，对IRI-2016模型在中国区域的电离层峰值参数NmF2和hmF2输出值随太阳活动强度、季节和昼夜的变化进行了对比分析。整体而言，模型输出值与COSMIC掩星观测值的相关性较高，但太阳活动低年的相关性低于太阳活动高年。两组数据拟合直线的斜率表明，IRI-2016模型输出值与COSMIC掩星观测值的一致性在太阳活动低年低于太阳活动高年。
文献[13]使用位于磁赤道附近的Jicamarca站垂测仪数据与IRI模型输出值进行验证分析，结果发现模型输出值与垂测仪观测值在太阳活动高年的相关性高于太阳活动低年，这与本文的研究结果一致。这主要是由于太阳活动低年，特别是2008年、2009年为太阳活动非寻常极小年，太阳活动水平极低且持续时间较长，IRI模型输出值会产生较大误差，甚至出现不合理的输出结果[3]。
太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布比较情况显示，虽然两者的空间分布形态基本一致，但仍存在明显差异。在15°N~30°N范围内，IRI-2016模型往往高估NmF2而低估hmF2。在30°N~55°N范围内，IRI-2016模型则存在低估NmF2而高估hmF2的现象。COSMIC掩星与IRI-2016模型之间的差异一方面是由于IRI-2016模型本身的误差，另一方面则是由于COSMIC掩星反演结果存在的误差。在球对称假设的前提下，COSMIC掩星反演的电子密度实际上是掩星切点附近整个掩星区域的平均值。因此，在电子密度梯度较大的区域，COSMIC掩星反演的电子密度与真实大气状态相比可能存在较大差异[14]。中国区域低纬度地区位于赤道电离层异常的驼峰区域附近，电子密度梯度很大，导致COSMIC反演的电子密度廓线相较于实际值可能存在较大偏差。
文献[15]使用中国区域位于同一经度带上的4个垂测仪台站（漠河站、北京站、武汉站和三亚站）数据对COSMIC掩星峰值参数进行了比较验证，发现在漠河站两者符合程度较高，但在武汉站和北京站中午时分COSMIC NmF2观测值高于垂测仪观测值，在三亚站中午时分COSMIC NmF2观测值低于垂测仪观测值。反演存在的误差可能使得在中纬度的武汉站（31.0°N，114.6°E）COSMIC反演的NmF2高于IRI-2016模型值，在低纬度的三亚站（18.3°N，109.6°E）COSMIC反演的NmF2低于IRI-2016模型值。文献[11]使用这4个垂测仪台站数据对IRI-2016模型计算hmF2的不同选项输出值结果进行比较分析，发现IRI-2016模型推荐的IRI_AMTB模型选项在4个台站均高估hmF2。本文结果也表明IRI-2016模型的hmF2输出值在30°N~55°N高于COSMIC掩星观测值。对于15°N~30°N范围，IRI-2016模型的hmF2输出值相对于掩星资料存在低估，原因可能是由于该范围位于赤道电离层异常驼峰附近，掩星反演的hmF2高于实际值[11]。
4 结语
本文对2008—2014年IRI-2016模型在中国区域的电离层峰值密度NmF2和峰值高度hmF2输出值与COSMIC掩星数据进行了对比，并对IRI-2016模型中计算这两个峰值参数的不同选项进行了比较分析。结论如下：
1）IRI-2016模型输出值与COSMIC掩星观测数据的一致性在太阳活动高年比太阳活动低年更好。模型NmF2输出值与COSMIC掩星观测值的相关性高于hmF2输出值与COSMIC掩星观测值的相关性。
2）在太阳活动低年春秋季的中午时段，与COSMIC掩星数据相比，IRI-2016模型在30°N~55°N对NmF2和hmF2分别存在低估和高估现象；在15°N~30°N区域内则恰恰相反。
3）以COSMIC掩星观测值为参照，IRI_CCIR模型和IRI_URSI模型均对NmF2在中午时分存在高估现象，该现象在太阳活动低年更为显著。而在春秋季的LT 17：00—04：00时间段，采用这两个模型选项得到的IRI-2016模型输出值均低于COSMIC观测值。
4）在太阳活动低年，无论任何季节，IRI_CCIR模型和IRI_AMTB模型的hmF2输出值均高于COSMIC掩星观测值，其中IRI_AMTB模型输出值相对于掩星观测的偏差更大，尤其是在冬季。IRI_Shubin模型的hmF2输出值与COSMIC观测值最为接近，即使在太阳活动高年，两者的偏差也小于其他两种模型。这表明COSMIC掩星资料在IRI_Shubin模型改进hmF2输出值中发挥了重要作用。
5）根据本文的研究结果，在中国区域采用IRI-2016计算NmF2和hmF2时，分别推荐使用IRI_CCIR和IRI_Shubin模型，而非官方推荐的IRI_URSI和IRI_AMTB模型。
-
![]()
图 3 2016年岗纳楼冰川地表温度变化
Figure 3. Surface Temperature Changes of the Gangnalou Glacier in 2016
![]()
图 4 2016年岗纳楼冰川表面流速（冰封期）
Figure 4. Surface Velocity of the Gangnalou Glacier in 2016 (Freezing Period)
![]()
图 5 2016年岗纳楼冰川表面流速（消融期）
Figure 5. Surface Velocity of the Gangnalou Glacier in 2016 (Ablation Period)
![]()
图 6 冰川消融与卫星视线向位移关系
Figure 6. Relationship Between Glacier Ablation and Satellite Line-of-Sight Displacement
![]()
图 9 岗纳楼冰川中线表面流速及冰面高程
Figure 9. Surface Velocity and Elevation Along the Central Line of the Gangnalou Glacier
表 1 Sentinel-1A升轨干涉宽幅单视复数影像数据
Table 1 Sentinel-1A Ascending Single Look Complex Images of Interferometric Wide Mode
编号 主影像日期 从影像日期 时间基线/d 垂直基线/m 极化方式 1 2016-01-11 2016-02-04 24 32.04 VV 2 2016-02-04 2016-02-28 24 85.78 VV 3 2016-02-28 2016-04-16 48 3.37 VV 4 2016-04-16 2016-05-10 24 0.17 VV 5 2016-05-10 2016-06-03 24 68.75 VV 6 2016-06-03 2016-07-21 48 67.62 VV 7 2016-07-21 2016-08-14 24 108.32 VV 8 2016-08-14 2016-09-07 24 26.45 VV 9 2016-09-07 2016-10-01 24 55.12 VV 10 2016-10-01 2016-10-25 24 87.04 VV 11 2016-10-25 2016-11-18 24 13.59 VV 12 2016-11-18 2016-12-12 24 107.46 VV 
下载: 导出CSV
表 2 冰川表面流速与均方根误差/(m·a-1)
Table 2 Glacier Surface Velocity and Its RMSE/(m·a-1)
时间跨度 距离向 方位向 流速 流速RMSE 流速 流速RMSE 2016-01-11—2016-02-04 25.89 1.02 40.14 3.50 2016-02-04—2016-02-28 25.22 0.66 62.95 3.35 2016-02-28—2016-04-16 25.90 0.31 25.79 1.92 2016-04-16—2016-05-10 16.34 1.36 62.35 3.83 2016-05-10—2016-06-03 24.41 6.63 92.55 37.50 2016-06-03—2016-07-21 27.34 8.65 43.15 16.46 2016-07-21—2016-08-14 13.83 1.83 78.40 7.75 2016-08-14—2016-09-07 16.74 2.19 40.78 14.77 2016-09-07—2016-10-01 20.98 6.61 64.38 25.71 2016-10-01—2016-10-25 38.65 6.23 42.87 16.21 2016-10-25—2016-11-18 23.57 3.94 46.62 7.22 2016-11-18—2016-12-12 35.00 1.36 47.44 4.87 平均值 24.49 3.40 53, 95 11.92
下载: 导出CSV
-
[1] Cuffey K M, Paterson W S B. The Physics of Glaciers[M]. NewYork, USA: Academic Press, 2010
[2] Bennett M R, Glasser N F. Glacial Geology: Ice Sheets and Landforms[M]. Chichester, UK: John Wiley & Sons, 2011
[3] Jacob T, Wahr J, Pfeffer W T, et al. Recent Contributions of Glaciers and Ice Caps to Sea Level Rise[J]. Nature, 2012, 482(7386): 514-518 doi: 10.1038/nature10847
[4] 周春霞, 邓方慧, 艾松涛, 等. 利用DInSAR的东南极极记录和达尔克冰川冰流速提取与分析[J]. 武汉大学学报·信息科学版, 2014, 39(8): 940-944 doi: 10.13203/j.whugis20130116 Zhou Chunxia, Deng Fanghui, Ai Songtao, et al. Determination of Ice-Flow Velocity at the Polar Record Glacier and Dalk Glacier Using DInSAR[J]. Geomatics and Information Science of Wuhan University, 2014, 39(8): 940-944 doi: 10.13203/j.whugis20130116
[5] Bechor N B D, Zebker H A. Measuring Two-Dimensional Movements Using a Single InSAR Pair[J]. Geophysical Research Letters, 2006, 33(16): L16311 doi: 10.1029/2006GL026883
[6] 王群, 张蕴灵, 范景辉, 等. 利用高分三号影像监测依嘎冰川表面运动[J]. 武汉大学学报·信息科学版, 2020, 45(3): 460-466 doi: 10.13203/j.whugis20190036 Wang Qun, Zhang Yunling, Fan Jinghui, et al. Monitoring the Motion of the Yiga Glacier Using GF-3 Images[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 460-466 doi: 10.13203/j.whugis20190036
[7] 张晓博, 赵学胜, 葛大庆, 等. 利用新型C波段雷达卫星研究南伊内里切克冰川运动特征[J]. 武汉大学学报·信息科学版, 2019, 44(3): 429-435 doi: 10.13203/j.whugis20160538 Zhang Xiaobo, Zhao Xuesheng, Ge Daqing, et al. Motion Characteristics of the South Inilchek Glacier Derived from New C-Band SAR Satellite[J]. Geomatics and Information Science of Wuhan University, 2019, 44(3): 429-435 doi: 10.13203/j.whugis20160538
[8] 邓方慧, 周春霞, 王泽民, 等. 利用偏移量跟踪测定Amery冰架冰流汇合区的冰流速[J]. 武汉大学学报·信息科学版, 2015, 40(7): 901-906 doi: 10.13203/j.whugis20130653 Deng Fanghui, Zhou Chunxia, Wang Zemin, et al. Ice-Flow Velocity Derivation of the Confluence Zone of the Amery Ice Shelf Using Offset-Tracking Method[J]. Geomatics and Information Science of Wuhan University, 2015, 40(7): 901-906 doi: 10.13203/j.whugis20130653
[9] Quincey D J, Luckman A, Benn D. Quantification of Everest Region Glacier Velocities Between 1992 and 2002, Using Satellite Radar Interferometry and Feature Tracking[J]. Journal of Glaciology, 2009, 55(192): 596-606 doi: 10.3189/002214309789470987
[10] Kumar V, Venkataramana G, Høgda K A. Glacier Surface Velocity Estimation Using SAR Interferometry Technique Applying Ascending and Descending Passes in Himalayas[J]. International Journal of Applied Earth Observation and Geoinformation, 2011, 13(4): 545-551 doi: 10.1016/j.jag.2011.02.004
[11] 王群. DInSAR和偏移量跟踪技术在山地冰川流速监测中的应用[D]. 北京: 中国地质大学(北京), 2018 Wang Qun. Application of DInSAR and Offset Tracking Techniques in Monitoring Mountain Glacier Velocity[D]. Beijing: China University of Geosciences(Beijing), 2018
[12] Li J, Li Z W, Wu L X, et al. Deriving a Time Series of 3D Glacier Motion to Investigate Interactions of a Large Mountain Glacial System with Its Glacial Lake: Use of Synthetic Aperture Radar Pixel Offset-Small Baseline Subset Technique[J]. Journal of Hydrology, 2018, 559: 596-608 doi: 10.1016/j.jhydrol.2018.02.067
[13] 丁一凡, 程晓, 程铖, 等. 基于偏移追踪的青藏高原卓琼冰川流速监测及成因分析[J]. 地球物理学报, 2017, 60(5): 1650-1658 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201705004.htm Ding Yifan, Cheng Xiao, Cheng Cheng, et al. Monitoring Ice Velocity by SAR Offset-Tracking and Analysis of Influence Factors for the Kangshung Glacier in the Tibetan Plateau[J]. Chinese Journal of Geophysics, 2017, 60(5): 1650-1658 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201705004.htm
[14] 李佳, 李志伟, 汪长城, 等. SAR偏移量跟踪技术估计天山南依内里切克冰川运动[J]. 地球物理学报, 2013, 56(4): 1226-1236 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201304019.htm Li Jia, Li Zhiwei, Wang Changcheng, et al. Using SAR Offset-Tracking Approach to Estimate Surface Motion of the South Inylchek Glacier in Tianshan[J]. Chinese Journal of Geophysics, 2013, 56(4): 1226-1236 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201304019.htm
[15] Goldstein R M, Engelhardt H, Kamb B, et al. Satellite Radar Interferometry for Monitoring Ice Sheet Motion: Application to an Antarctic Ice Stream[J]. Science, 1993, 262(5139): 1525-1530 doi: 10.1126/science.262.5139.1525
[16] 武安斌. 疏勒南山岗纳楼5号冰川现代冰碛物的沉积组构特征[J]. 兰州大学学报, 1980, 16(3): 119-126 https://www.cnki.com.cn/Article/CJFDTOTAL-LDZK198003013.htm Wu Anbin. The Characteristics of Sedimentary Fabric of Modern till of 5th Gangnalou Glacier, Shuleinanshan[J]. Journal of Lanzhou University, 1980, 16(3): 119-126 https://www.cnki.com.cn/Article/CJFDTOTAL-LDZK198003013.htm
[17] 吴红波, 杨萌, 杨春利, 等. 冰川亏损对哈拉湖流域湖泊水位波动的影响[J]. 水资源与水工程学报, 2016, 27(4): 13-20 https://www.cnki.com.cn/Article/CJFDTOTAL-XBSZ201604003.htm Wu Hongbo, Yang Meng, Yang Chunli, et al. Effect of Glacier Mass Loss on Water Level Fluctuation of Lake in Hala Lake Basin[J]. Journal of Water Resources and Water Engineering, 2016, 27(4): 13-20 https://www.cnki.com.cn/Article/CJFDTOTAL-XBSZ201604003.htm
[18] 刘时银, 郭万钦, 许君利, 等. 中国第二次冰川编目数据集(V1.0)[EB/OL]. [2020-10-30]. https://data.tpdc.ac.cn/zh-hans/data/f92a4346-a33f-497d-9470-2b357ccb4246 Liu Shiyin, Guo Wanqin, Xu Junli, et al. The Second Glacier Inventory Dataset of China (Version 1.0) [EB/OL]. [2020-10-30]. https://data.tpdc.ac.cn/zh-hans/data/f92a4346-a33f-497d-9470-2b357ccb4246
[19] 李鹏, 李振洪, 施闯, 等. 中国地区30 m分辨率SRTM质量评估[J]. 测绘通报, 2016(9): 24-28 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB201609007.htm Li Peng, Li Zhenhong, Shi Chuang, et al. Quality Evaluation of 1 Arc Second Version SRTM DEM in China[J]. Bulletin of Surveying and Mapping, 2016(9): 24-28 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB201609007.htm
[20] 闫世勇. 山地冰川表面运动雷达遥感监测方法研究[D]. 北京: 中国科学院大学, 2013 Yan Shiyong. Research on Extraction of Alpine Glacier Surface Movement by SAR Remote Sensing[D]. Beijing: University of Chinese Academy of Sciences, 2013
[21] Geudtner D, Torres R, Snoeij P, et al. Sentinel-1 System Capabilities and Applications[C]// IEEE Geoscience and Remote Sensing Symposium, Quebec City, Canada, 2014
[22] 吴文豪. 哨兵雷达卫星TOPS模式干涉处理研究[D]. 武汉: 武汉大学, 2016 Wu Wenhao. TOPS Interferometry with Sentinel-1[D]. Wuhan: Wuhan University, 2016
[23] 黄震, 赵双明, 常青. 国产高分三号卫星干涉测量试验[J]. 测绘地理信息, 2021, 46(4): 7-11 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202104003.htm Huang Zhen, Zhao Shuangming, Chang Qing. Interferometric Test Using Domestic GF-3 Satellite[J]. Journal of Geomatics, 2021, 46(4): 7-11 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202104003.htm
[24] Sansosti E, Berardino P, Manunta M, et al. Geometrical SAR Image Registration[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(10): 2861-2870 doi: 10.1109/TGRS.2006.875787
[25] Li J, Li Z W, Ding X L, et al. Investigating Mountain Glacier Motion with the Method of SAR Intensity-Tracking: Removal of Topographic Effects and Analysis of the Dynamic Patterns[J]. Earth-Science Reviews, 2014, 138: 179-195
[26] Wegmuller U, Werner C, Strozzi T, et al. Ionospheric Electron Concentration Effects on SAR and InSAR[C]// IEEE International Symposium on Geoscience and Remote Sensing, Denver, CO, USA, 2006
[27] Strozzi T, Wegmuller U, Werner C, et al. Potential of a C-Band SAR Mission with 12-Day Repeat Cycle to Derive Ice Surface Velocity with Interfero-metry and Offset Tracking[C]// IEEE International Geoscience and Remote Sensing Symposium, Barcelona, Spain, 2007
[28] 王思胜, 江利明, 孙永玲, 等. 基于ALOS PALSAR数据的山地冰川流速估算方法比较: 以喀喇昆仑地区斯克洋坎力冰川为例[J]. 国土资源遥感, 2016, 28(2): 54-61 https://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201602009.htm Wang Sisheng, Jiang Liming, Sun Yongling, et al. Evaluation of Methods for Deriving Mountain Glacier Velocities with ALOS PALSAR Images: A Case Study of Skyang Glacier in Central Karakoram[J]. Remote Sensing for Land & Resources, 2016, 28(2): 54-61 https://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201602009.htm
[29] Nye J F. The Mechanics of Glacier Flow[J]. Journal of Glaciology, 1952, 2(12): 82-93 doi: 10.1017/S0022143000033967
-
期刊类型引用(2)
1. 陆浩楠，叶世榕，张庆兰，夏朋飞，鄂盛龙. 利用COSMIC-1/2 IonPrf产品分析中国区域近15年电离层时空变化特征. 武汉大学学报(信息科学版). 2024(05): 765-774 . 
百度学术
2. 唐丝语，周春红，黄智，袁洪. 国际参考电离层模型输出参数峰值高度性能. 空间科学学报. 2021(06): 898-904 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 817
- HTML全文浏览量: 218
- PDF下载量: 93
- 被引次数: 4
