COMPUTATION AND MAPPING OF THE ELEMENTS OF THE WATER-ICE BUDGET OF GLACIERS

Two general approaches for compiling maps of components of the water-ice budget in regions of glacier nourishment are compared and the resulting maps are analyzed. The compilation of raw data is based on field measurements rather than remote sensing information. Translated by Edward Torrey, Alexandria, VA 22308 from: Materialy glyatsiologicheskikh issledovaniy, 1991, No. 72, pp. 94-102.     

Multiparametric Cartographic Visualisation of Glacier Rheology

Abstract

In cooperation between remote sensing experts and cartographers interested in glaciology, new types of maps showing the glacier dynamics have been developed. The maps make use of the original phase gradient approach to glacier rheology modelling based on repeat-pass ERTS SAR interferograms. Careful map design and, in particular, colour assignment allow the visualisation of the glacier dynamics in its locally changing velocity with an estimated accuracy of approximately 2.0 cm per day. Two map derivates – a differential interferogram showing the glacier velocity and another product displaying the glacier strain rate – have been designed. Moreover, maps displaying the glacier marginal changes within the space of four years have been generated. The strain rate maps evidence that spots with high values frequently correspond with crevasse-prone areas which are even detectable under thick layers of snow. In this sense, the latter visualisations can be seen as maps of crevasse danger zones. The Svartisen in Norway and the Hintereis Glacier in Austria served as testbeds for the development of these different types of maps which are at the scales of 1:25 000, 1:50 000 and 1:100 000.     

Glacier changes using satellite data and effect of climate in Tirungkhad basin located in western Himalaya

Glaciers are widely recognized as key indicators of climate change, and melt water obtained from them is an important source of fresh water and for hydropower generation. Regular monitoring of a large number of Himalayan glaciers is important for improving our knowledge of glacier response to climate change. In the present study, Survey of India topographical maps (1966) and Landsat datasets as ETM+ (2000, 2006) and TM (2011) have been used to study glacier fluctuations in Tirungkhad basin. A deglaciation of 26.1% (29.1?km²) in terms of area from 1966 to 2011 was observed. Lower altitude small glaciers (area?<?1?km²) lost more ice (34%), while glaciers with an area <10?km² lost less (20%). The percentage of change in glacier length was 26% (31.9?km) from 1966 to 2011. The south-facing glaciers showed high percentages of loss. From 2000 to 2011, debris cover has increased by 1.34%. The analysis of the trend in meteorological data collected from Kalpa and Purbani stations was carried out by Mann Kendall non-parametric method. During the last two decades, the mean annual temperature (T_max and T_min) has increased significantly, accompanied with a fall in snow water equivalent (SWE) and rainfall. The increasing trend in temperature and decreasing trend in SWE were significant at 95% confidence level. This observation shows that the warming of the climate is probably one of the major reasons for the glacier change in the basin.     

Estimation of Snow Cover Area,Snow Physical Properties and Glacier Classification in Parts of Western Himalayas Using C-Band SAR Data

Snow physical properties, snow cover and glacier facies are important parameters which are used to quantify snowpack characteristics, glacier mass balance and seasonal snow and glacier melt. This study has been done using C-band synthetic aperture radar (SAR) data of Indian radar imaging satellite, radar imaging satellite-1 (RISAT)-1, to estimate the seasonal snow cover and retrieve snow physical properties (snow wetness and snow density), and glacier radar zones or facies classification in parts of North West Himalaya (NWH), India. Additional SAR data used are of Radarsat-2 (RS-2) satellite, which was used for glacier facies classification of Smudra Tapu glacier in Himachal Pradesh. RISAT-1 based snow cover area (SCA) mapping, snow wetness and snow density retrieval and glacier facies classification have been done for the first time in NWH region. SAR-based inversion models were used for finding out wet and dry snow dielectric constant, dry and wet SCA, snow wetness and snow density. RISAT-1 medium resolution scan-SAR mode (MRS) in HV polarization was used for first time in NWH for deriving time series of SCA maps in Beas and Bhagirathi river basins for years 2013–2014. The SAR-based inversion models were implemented separately for RISAT-1 quad pol. FRS2, for wet snow and dry snow permittivity retrieval. Masks for layover and shadow were considered in estimating final snow parameters. The overall accuracy in terms of R² value comes out to be 0.74 for snow wetness and 0.72 for snow density based on the limited ground truth data for subset area of Manali sub-basin of Beas River up to Manali for winter of 2014. Accuracy for SCA was estimated to be 95 % when compared with optical remote sensing based SCA maps with error of ±10 %. The time series data of RISAT-1 MRS and hybrid data in RH/RV mode based decompositions were also used for glacier radar zones classification for Gangotri and Samudra Tapu glaciers. The various glaciers radar zones or facies such as debris covered glacier ice, clean or bare glacier ice radar zone, percolation/refreeze radar zone and wet snow, ice wall etc., were identified. The accuracy of classified maps was estimated using ground truth data collected during 2013 and 2014 glacier field work to Samudra Tapu and Gangotri glaciers and overall accuracy was found to be in range of 82–90 %. This information of various glacier radar zones can be utilized in marking firn line of glaciers, which can be helpful for glacier mass balance studies.     

AN EXPERIMENT IN THE COMPILATION AND USE OF A COMPLEX ENGINEERING-GLACIOLOGICAL MAP

A method for the compilation of 1:3,000,000-scale maps of natural nival-glacial phenomena is proposed for an area of the High Pamir Range in Soviet Central Asia. Such maps depict the distribution, frequency and intensity/magnitude, and potential impacts of such natural hazards as avalanches, mudslides, “dirty snow” avalanches (combined avalanches-mudslides in late spring), icings (naledy), and glacier pulsations, and are used in the planning of construction and the maintenance of roads and other transportation corridors during winter and spring. Hypothetical examples of the use of the maps are provided. Translated from: Materialy glyatsiologicheskikh issledovanly, USSR Geophysical Committee Publication No. 55. Moscow: Mezhduvedomstvenniy geofizicheskiy komitet, 1986, pp. 213–219.     

Evidence for Mountain Glacier Changes in Semi-arid Environments based on Remote Sensing Data

Countries like Iran, which are geographically situated in a rather arid and warm regions, will suffer more from global warming than countries located in humid and semi-humid regions. In such environments, analyzing the variations of mountain glaciers can reveal several aspects of climate change patterns more efficiently in comparison to the other geo-indicators. The present study reports some evidence of changes for Alamkouh glacier between 1955 and 2010 based on several mediums to high-resolution satellite images. Considering that most part of the Alamkouh glacier is covered by debris and delineating its actual area is not possible, planimetric change analysis was restricted to the clean-ice regions. The object-oriented classification approach was used to estimate the clean ice areas. This technique takes into account the shapes of the features along with their spectral patterns. Results revealed that clean ice regions of Alamkouh glacier shrank since 1955 with an overall area reduction of about 59 %. Although the general observed trend is a clean ice area reduction, some advancement was detected over the period from 2000 to 2010. During 1992–2000, the maximum reduction in the clean ice area was observed (0.08 km².a^?1). However, clean ice area of the case study has partially increased about 0.028 km².a^?1 from 2000 to 2010. Supra-glacial lake change analysis illustrated that at the surface of the glacier, lakes have been enlarged remarkably in the past 55 years (about 4.75 times greater). In addition, clean ice area and the surface area of supra-glacial lakes oscillated in compliance with each other. The findings revealed that the maximum expansion of supra-glacial lake occurred during 1992–2000, which demonstrate the glacier maximum reduction during this period. This shrinkage in the Alamkouh glacier caused an extensive glacial lake outburst flood in Jun 2011. The results of this study agree with documented changes in other mountain glaciers located in arid and semi-arid environments and they also confirm the application of mountain glaciers in climate variations monitoring over such regions.     

Estimating surface ice velocity on Chhota Shigri glacier from satellite data using Particle Image Velocimetry (PIV) technique

Information about the surface ice velocity is one of the important parameters for Mass balance and Glacier dynamics. This study estimates the surface ice velocity of Chhota Shigri glacier using Landsat (TM/ETM+) and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) temporal data-sets from a period of 2009 to 2016 and 2006 to 2007, respectively. A correlation based Particle Image Velocimetry (PIV) technique has been used for the estimation of surface ice velocity. This technique uses multiple window sizes in the same data-set. Four window sizes (low, medium, high, very high) are used for each image pair. Estimated results have been compared with the published data. The outcomes attained from the medium window size closely matches with the published results. The estimated mean surface ice velocities of medium window size are 24 and 28.5 myr^?1 for 2009/2010 and 2006/2007 images pair. Highest velocity is observed in middle part of the glacier while lowest in the accumulation zone of the glacier.     

Analyzing ice dynamics using Sentinel-1 data at the Solheimajoküll Glacier,Iceland

ABSTRACT

The climate in southern Iceland has warmed over the last 70 years, resulting in accelerated glacier dynamics at the Solheimajoküll glacier. In this study, we compare glacier terminus locations from 1973 to 2018, to changes in climate across the study area, and we derive ice-surface velocities (2015–2018) from satellite remote-sensing imagery (Sentinel-1) using the offset-tracking method. There have been two regional temperature trends in the study period: cooling (1973–1979) and warming (1980–2018). Our results indicate a time lag of about 20 years between the onset of glacier retreat (?53 m/year since 2000) and the inception of the warming period. Seasonally, the velocity time series suggest acceleration during the summer melt season since 2016, whereas glacier velocities during accumulation months were constant. The highest velocities were observed at high elevations where the ice-surface slope is the steepest. We tested several scenarios to assess the hydrological time response to glacier accelerations, with the highest correlations being found between one and 30 days after the velocity estimates. Monthly correlation analyses indicated inter-annual and intra-annual variability in the glacier dynamics. Additionally, we investigate the linkage between glacier velocities and meltwater outflow parameters as they provide useful information about internal processes in the glacier. Velocity estimates positively correlate with water level and negatively correlate with water conductivity between April and August. There is also a disruption in the correlation trend between water conductivity and ice velocity in June, potentially due to a seasonal release of geothermal water.     

Monitoring thickness and volume changes of the Dongkemadi Ice Field on the Qinghai-Tibetan Plateau (1969–2000) using Shuttle Radar Topography Mission and map data

Abstract

This paper presents the first measurement of multi-decadal thickness and volume changes (1969–2000) of the Dongkemadi Ice Field (DIF) in the Tanggula Mountains, central Qinghai-Tibetan Plateau, China, using multi-source remote sensing data. These include the Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) acquired in February, 2000, a DEM generated by digitising analogue topographic maps from 1969, and Landsat ETM+ imagery from 2000. Digital glacier outlines and GIS-based processing were used to calculate an elevation difference map to evaluate the relative elevation error of these two DEMs over ice-free areas. This method was also used to identify regions of glacier elevation thinning and thickening corresponding to glacier mass loss and gain. Analysis of 67,520 points on flat grass and rock terrain surrounding the DIF, with a slope less than 25°, showed a mean elevation difference of –0.90 m and a standard deviation of 5.58 m. A thickness change error within ±6 m was estimated. Between 1969 and 2000, 76.51% of the whole DIF area appeared to be thinning while 23.49% showed thickening. The average glacier surface thinning was –12.58 m with a standard deviation of 18.29 m and the estimated volume loss was 1.17 km³. The standard deviation of volume change was 0.0006 km³ over the DIF. A thinning rate up to 0.41±0.194 m a^?1 or 0.038 km³ a^?1 for the volume loss was observed for the whole ice field, which seems to be evidence for the ongoing retreat of glaciers on the Qinghai-Tibetan Plateau. It was found that the spatial thickness change pattern derived from the remote sensing method was consistent with the thickness change results of the Small Dongkemadi Glacier (SDG) from field measurements. The estimated error of the annual thickness change rate was on the order of 5%. The relationship between elevation change and absolute glacier elevation over typical glaciers was also analysed, showing considerable variability. These changes have possibly resulted from increased temperature and decreased precipitation in this region.     

基于Landsat8 OLI遥感影像冰湖边界分类叠加提取

青藏高原特殊的地理及其气候环境适合冰湖变化研究。传统的方法由于受到冻湖、冰川以及山体的阴影的影响,提取冰湖边界时会出现大量的错分和误分。本文基于NDWI(Normalized Difference Water Index,水体指数法)、DEM(数字高程模型)、RGI(Randolph Glacier Inventory)6.0全球冰川编目数据提出针对不同状态下的冰湖边界提取的分类叠加提取方法。结果表明,分类叠加提取方法能够实现未结冰、半结冰以及全结冰三种状态下的冰湖边界提取,并且能减少冰川以及山体阴影的影响,使用该方法提取三种状态下的冰湖的错分和漏分误差均小于2%且Kappa系数均大于0.81。     

SATELLITE MEASUREMENTS OF LAND PARAMETERS AS PART OF AN INTERNATIONAL PROGRAM OF CLIMATE RESEARCH

The role of remote sensing in one aspect of the World Climate Research Program is summarized. More specifically, unresolved theoretical and methodological problems are outlined for the use of remote sensing imagery (in conjunction with the International Land Surface Climatology Project) to derive quantitative parameters of land (and ice) surface characteristics. These parameters are raw data for models of the overall circulation of the atmosphere, which can be used to generate maps of albedo, surface roughness, hydrology, etc. A list of preliminary standards for satellite monitoring of land surface characteristics is included. Translated from: Izvestiya Vsesoyuznogo Geograficheskogo Obshchestva, 1988, No. 2, pp. 161-165.     

Contrasting Signals of Glacier Changes in Zanskar Valley,Jammu & Kashmir,India Using Remote Sensing and GIS

The study of advancement and recession of the glaciers in the Himalayas is essential due to their contrasting response towards climatic change. In the present study, Survey of India (SOI) topographical maps of 1962, IRS: LISS-III image of 2001 and LANDSAT-5: TM (Thematic Mapper) image of 2009 were used to analyze the glacier fluctuations in a part of Zanskar valley. The analysis carried out on 212 glaciers indicated decrease of 57 km² (8 %) of glacier area over many glacier which was partly compensated with area increase by 42 km² (6 %) in other glaciers, resulting an overall glacier area decrease by only 15 km² (2 %) from 1962–2001. Due to glacier fragmentation the number of glaciers increased from 212 in 1962 to 238 by 2001. Although majority of glaciers (88 %) exhibited retreat (up to 13 my^?1), minor advancement (<15 my^?1) also took place in few glaciers during this period. Advancement took place mainly in larger glaciers (2–5 km² and >5 km²) located over wider altitudinal range (700 m–1,000 m) whereas smaller glaciers (<2 km²) with narrow altitudinal range (100 m–500 m) exhibited retreat. The supraglacial debris analysis indicated that percentage of debris cover over glaciers ranges from 1.43 % to 18.15 %. Smaller glaciers (<2 km²) were debris free in comparison to the larger glaciers (>5 km²). During 2001–2009 majority of the glaciers were apparently stable in terms of their area and snout position indicating less impact of climate forcing in parts of Zanskar valley as compared to other parts of the Himalaya.     

基于SAR偏移量跟踪法提取岗纳楼冰川流速

冰川动态变化监测有助于反映全球和区域气候演变,保护自然环境和自然资源。近年来,基于SAR数据研究冰川运动已成为主流技术之一。基于SAR提取冰川流速主要包括合成孔径雷达干涉测量、多孔径雷达干涉和偏移量追踪法。本文采用SAR偏移量追踪法中的强度追踪法,提取青海省哈拉湖东北部岗纳楼冰川沿距离向、方位向的冰川流速。试验结果表明,距离向冰川运动速度提取效果较好,最大流速达15.36 m/a,流速从中轴向两侧递减,在冰舌末端趋于0;方位向提取的冰流速最大达18.27 m/a,但因电离层干扰,方位向流速图中存在一些方位向条纹。此外,由于冰流速在方位向分量小等因素的影响,本文研究提取的方位向流速精度低于距离向。     

NATURAL ENVIRONMENT AND RESOURCES OF THE WORLD—A NEW SOVIET ATLAS

The authors employ remote sensing imagery and existing landscape maps to gain greater insight into how pre-glacial tectonic activity contributes to landscape development in old glaciated regions by affecting the form and movements of lobes of past continental ice sheets. Using the Belorussian SSR as a study area, the authors identify several features of current landscape structure which can serve as indicators of past tectonic activity, the most reliable being various characteristics of the landscape family—defined as a combination of natural territorial complexes of similar origin and age. Translated from: Izvestiya Vsesoyuznogo Geograficheskogo Obshchestva, 1990, No. 6, pp. 560–565.     

Elevation change around Dome A region of Antarctica from EnviSat satellite radar altimetry during 2002–2012

The expected responses of ice sheets to climate warming are growth in the thickness of the inland ice areas and thinning near the margins. In recent decades, researchers have identified glacier acceleration along Antarctic ice sheet coastal margins. However, the study of ice sheet interiors where seasonal accumulation eventually balances ice wastage at the lower elevation is poorly understood. In this paper, the ice sheet elevation change around Dome A region is analyzed from 2002 to 2012 using two million elevation change measurements from EnviSat satellite radar altimeter data covering an area of about 7000 km². A declining trend of 0.572 ± 1.31 mm/year which means that the Dome A region was in balance during the last decade can be captured. In addition, two obvious changes in accumulation which divide elevation change time series into three independent equilibration stages are also extracted. In order to explain this phenomenon, two speculations related to snowfall and firn compaction are proposed in this paper.     

USE OF SPACE IMAGERY IN THE STUDY OF MASS EXCHANGE IN GLACIER SYSTEMS

The author surveys the use of remote sensing imagery in the study of mass exchange in glaciers, i.e., glacier dynamics resulting from phase changes in the water which they contain. A program of research now underway at the Institute of Geography (USSR Academy of Sciences) in conjunction with the “Priroda” Remote Sensing Center (Moscow) focuses on improving methods of data collection and mapping of glacier dynamics from space imagery, and particularly on identifying natural glacioclimatic zones believed to represent specific mass exchange conditions or regimes. The boundary of glacier nourishment separating areas of accumulation and ablation represents a key glacioclimatic indicator of mass exchange on satellite imagery and aerial photography. Translated from: Geodeziya i kartografiya, 1986, No. 5, pp. 26–31.     

Light mode and dark mode: Which one is suitable when using public-facing web maps? An experimental evaluation using eye-tracking

Dark maps, which are also called dark modes or dark themes, have gained popularity, but their usability has not been experimentally evaluated. This article presents a user experiment that uses eye-tracking to assess the usability of dark and light maps. Here, two widely used web maps, Google Maps and Gaode Maps, are selected at the street and city scales. Eight map-use tasks are designed and cover four operations (identify, compare, rank, and associate) with space-alone and attributes-in-space distinctions. Four pairs of map-use tasks (light-during-the-day, dark-during-the-day, dark-at-night, and light-at-night) are examined from three aspects: effectiveness, efficiency, and cognitive load. The results provide preliminary evidence that the light-during-the-day performance is generally the best in most cases, followed by the dark-at-night performance; the dark-during-the-day performance is the worst in all cases, followed by the light-at-night performance, which suggests that aligning the map design with the environment (i.e., lighting environment) is critical for better communication. The light-during-the-day performance is the best for space-alone tasks, and the dark-at-night performance is the best for attributes-in-space tasks. Our investigation also indicates that dark maps are far less favored in practice, which suggests that users' preference for using the dark mode of public-facing web maps needs to be shaped. Since light and dark maps are associated with photopic and scotopic vision, respectively, the results indicate the need for future studies on how to leverage scotopic vision to design better dark maps.     

极区海冰密集度AMSR-E数据反演算法的试验与验证

海冰密集度是极区海冰监测的重要参数,目前分辨率最高的微波海冰密集度产品为德国Bremen大学发布的针对AMSR-E 89 GHz频段数据利用ASI算法反演的网格数据。为实现中国极区遥感产品从无到有的战略步骤,本文针对AMSR-E 89GHz频段微波数据的ASI算法,进行了插值算法、系点值和天气滤波器一系列试验。针对北极海区,着重对影响反演结果的主要参数——纯冰和纯水的亮温极化差异阈值,即系点值(P₁和P₀)进行了2009年全年的统计分析。研究表明,2009年北极纯冰和纯水的代表区域P₁和P₀年平均值分别为10.0 K和46.67 K;2 K以上的系点值差异引起的海冰密集度差别较为显著;同样的系点值差异在不同极化差异P取值范围对海冰密集度的影响也不同。通过统计确定的系点值推算并修正了海冰密集度反演公式,对2009年全年北极海冰密集度进行了反演,并与Bremen大学产品进行了比较。继而对白令海和楚科奇海12个晴空下MODIS可见光样本数据进行反演,以验证AMSR-E冰密集度反演结果,并对误差原因进行了分析。本研究反演结果与MODIS样本比对的误差略小于Bremen大学的反演产品,空间平均误差为3.84%,空间平均绝对误差10.83%。     

3-Decadal Changes in Sea Ice Melting and Freezing Pattern in Antarctica Using SSM/I Data

The present study aims to estimate the long term changes in sea ice areal extent (SIAE) around Antarctica using Special Sensor Microwave Imager (SSM/I) satellite data (1988–2011). Daily, weekly and monthly brightness temperature (T_B) maps of Antarctica have been generated. Temporal variation in monthly averaged T_B values of sea ice have been analyzed and an increasing trend is observed in T_B. Sea ice concentration (SIC) values are estimated and the concentration class 0.85–0.95 is observed as the main contributor in sea ice areal extent during the Antarctic peak winter month i.e. September and October. SIC changes rapidly from October onwards and shows fast melting of sea ice, however, from April start of freezing has been observed. Continuous increase in SIAE is observed from February to September in all the analyzed years (1988–2011) data. In this observation period maximum SIAE is observed in September 1989 and minimum in February 1993. The Moderate Resolution Imaging Spectroradiometer (MODIS) data is used for the validation of SIAE values. SSM/I and MODIS estimated SIAE values are compared (R²?=?0.98) and a relation between them is established by performing regression analysis. This relation is further used for the estimation of corrected SIAE values. The monthly average sea ice extent maps of last three decades have been analyzed and a shift in onset of freezing of sea ice has been observed. The period of frozen sea has also been observed as decreased.     

非球形冰晶的毫米波k-Ze关系研究

针对毫米波雷达处理数据的实际需要,应用离散偶极子近似法DDA,获得了非球形冰晶的后向散射及衰减截面并进行了参数化,并主要基于细化的冰云模型,假设冰晶粒子谱为Γ分布,通过模拟取样各1330次（代表1330种粒子分布）,分别计算得到了W波段（94 GHz）与Ka波段（35 GHz）毫米波雷达探测的冰云衰减系数k及雷达反射率因子Z_e,而且利用数值模拟的方法,建立了k-Z_e关系的具体表达式。计算表明,非球形和非瑞利散射对W波段毫米波雷达衰减的影响较大,而且在同样滴谱分布条件下,W波段毫米波雷达的衰减比Ka波段毫米波雷达的大几倍,此外细化的冰云模型对k-Z_e关系具有影响。本研究对中纬度非降水性冰云的毫米波雷达的衰减订正具有参考价值,并对中国的毫米波雷达应用具有借鉴作用。