Ice Flow Speed Research Articles

Rapid deceleration in the ice-flow velocity of the Shirase Glacier ice tongue and its influence on the velocity field: Observations from Sentinel-1 C-SAR

Temporal and spatial variations in the ice-flow speed of Shirase Glacier ice tongue in East Antarctica between July 2018 and December 2021 were investigated using Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) imagery. We identified pronounced slowdown events in the eastern part of the outer ice tongue, 30–40 km from the grounding line in 2020 and 55 km in 2021. Comparison of ice thickness and bathymetry in areas where the deceleration events occurred suggests that the events were caused by icebergs grounding or landing on the seafloor. The absence of slowdown propagation towards the grounding line demonstrates the ice tongue offers very limited buttressing. This study contributes to a better understanding of the factors influencing glacier dynamics, particularly in the context of grounding events and their localized impacts.

Triggers of the 2022 Larsen B multi-year landfast sea ice breakout and initial glacier response

Abstract. In late March 2011, landfast sea ice (hereafter, “fast ice”) formed in the northern Larsen B embayment and persisted continuously as multi-year fast ice until January 2022. In the 11 years of fast-ice presence, the northern Larsen B glaciers slowed significantly, thickened in their lower reaches, and developed extensive mélange areas, leading to the formation of ice tongues that extended up to 16 km from the 2011 ice fronts. In situ measurements of ice speed on adjacent ice shelf areas spanning 2011 to 2017 show that the fast ice provided significant resistive stress to ice flow. Fast-ice breakout began in late January 2022 and was closely followed by retreat and breakup of both the fast-ice mélange and the glacier ice tongues. We investigate the probable triggers for the loss of fast ice and document the initial upstream glacier responses. The fast-ice breakup is linked to the arrival of a strong ocean swell event (>1.5 m amplitude; wave period waves >5 s) originating from the northeast. Wave propagation to the ice front was facilitated by a 12-year low in sea ice concentration in the northwestern Weddell Sea, creating a near-ice-free corridor to the open ocean. Remote sensing data in the months following the fast-ice breakout reveals an initial ice flow speed increase (>2-fold), elevation loss (9 to 11 m), and rapid calving of floating and grounded ice for the three main embayment glaciers Crane (11 km), Hektoria (25 km), and Green (18 km).

The role of near-terminus conditions in the ice-flow speed of Upernavik Isstrøm in northwest Greenland

Abstract Upernavik Isstrøm, the largest contributor to sea-level rise in northwest Greenland, has experienced complex and contrasting ice-flow-speed changes across its five outlets over the last two decades. In this study, we present a detailed remote-sensing analysis of the ice dynamics at Upernavik's outlets from 2000 to 2021 to evaluate the details of these changes. Previous research suggested that the presence or absence of floating ice tongues strongly influences Upernavik's ice dynamics. We use several lines of evidence to document the presence of floating ice tongues, and find that, while several outlets experienced ice-tongue formation and/or loss during the study period, these changes do not explain observed fluctuations in ice-flow velocity. Further exploration of ice-dynamic forcings using a flowline model suggests that changes in basal slipperiness near the terminus have a strong impact on upstream ice dynamics and can explain the velocity variations. Our results suggest that speed fluctuations at Upernavik's outlets may be seasonally and interannually controlled by bed conditions near the terminus, and highlight the need for further research on the influence of basal conditions on complex tidewater glacier dynamics.

A new model for supraglacial hydrology evolution and drainage for the Greenland Ice Sheet (SHED v1.0)

Abstract. The Greenland Ice Sheet (GrIS) is losing mass as the climate warms through both increased meltwater runoff and ice discharge at marine-terminating sectors. At the ice sheet surface, meltwater runoff forms a dynamic supraglacial hydrological system which includes stream and river networks and large supraglacial lakes (SGLs). Streams and rivers can route water into crevasses or into supraglacial lakes with crevasses underneath, both of which can then hydrofracture to the ice sheet base, providing a mechanism for the surface meltwater to access the bed. Understanding where, when, and how much meltwater is transferred to the bed is important because variability in meltwater supply to the bed can increase ice flow speeds, potentially impacting the hypsometry of the ice sheet in grounded sectors, and iceberg discharge to the ocean. Here we present a new, physically based, supraglacial hydrology model for the GrIS that is able to simulate (a) surface meltwater routing and SGL filling; (b) rapid meltwater drainage to the ice sheet bed via the hydrofracture of surface crevasses both in and outside of SGLs; (c) slow SGL drainage via overflow in supraglacial meltwater channels; and, by offline coupling with a second model, (d) the freezing and unfreezing of SGLs from autumn to spring. We call the model the Supraglacial Hydrology Evolution and Drainage (or SHED) model. We apply the model to three study regions in southwest Greenland between 2015 and 2019 (inclusive) and evaluate its performance with respect to observed supraglacial lake extents and proglacial discharge measurements. We show that the model reproduces 80 % of observed lake locations and provides good agreement with observations in terms of the temporal evolution of lake extent. Modelled moulin density values are in keeping with those previously published, and seasonal and inter-annual variability in proglacial discharge agrees well with that which is observed, though the observations lag the model by a few days since they include transit time through the subglacial system, while the model does not. Our simulations suggest that lake drainage behaviours may be more complex than traditional models suggest, with lakes in our model draining through a combination of both overflow and hydrofracture and with some lakes draining only partially and then refreezing. This suggests that, in order to simulate the evolution of Greenland's surface hydrological system with fidelity, a model that includes all of these processes needs to be used. In future work, we will couple our model to a subglacial model and an ice flow model and thus use our estimates of where, when, and how much meltwater gets to the bed to understand the consequences for ice flow.

Observing tidal effects on the dynamics of the Ekström Ice Shelf with focus on quarterdiurnal and terdiurnal periods

AbstractAntarctica's ice shelves stabilize the ice sheet and, therefore, understanding processes affecting the mass budgets of ice shelves is important for estimating grounded ice loss. To study the ice shelf dynamics, we analyzed seismological and GNSS data from the Ekström Ice Shelf in Dronning Maud Land. We extracted probabilistic power spectral densities (PPSD) in the frequency band 3.4–6.8 Hz, typical of icequakes, from seismological data and observed pronounced signals in the PPSD with near 3 and 4 cycles per day (cpd) corresponding to tidal overharmonics, in addition to the main tidal constituents near 1 and 2 cpd. GNSS data reveal the same components in ice flow speed but not in vertical displacements. Generally, tide-induced grounding line migration modulates the flow velocity of an entire ice shelf. We find that this velocity modulation causes the increased icequake activity in the tidal overharmonics with 3 and 4 cpd in an ice shear zone where the flow velocity drops to nearly zero.

Glaciological history and structural evolution of the Shackleton Ice Shelf system, East Antarctica, over the past 60 years

Abstract. The discovery of Antarctica's deepest subglacial trough beneath the Denman Glacier, combined with high rates of basal melt at the grounding line, has caused significant concern over its vulnerability to retreat. Recent attention has therefore been focusing on understanding the controls driving Denman Glacier's dynamic evolution. Here we consider the Shackleton system, comprised of the Shackleton Ice Shelf, Denman Glacier, and the adjacent Scott, Northcliff, Roscoe and Apfel glaciers, about which almost nothing is known. We widen the context of previously observed dynamic changes in the Denman Glacier to the wider region of the Shackleton system, with a multi-decadal time frame and an improved biannual temporal frequency of observations in the last 7 years (2015–2022). We integrate new satellite observations of ice structure and airborne radar data with changes in ice front position and ice flow velocities to investigate changes in the system. Over the 60-year period of observation we find significant rift propagation on the Shackleton Ice Shelf and Scott Glacier and notable structural changes in the floating shear margins between the ice shelf and the outlet glaciers, as well as features indicative of ice with elevated salt concentration and brine infiltration in regions of the system. Over the period 2017–2022 we observe a significant increase in ice flow speed (up to 50 %) on the floating part of Scott Glacier, coincident with small-scale calving and rift propagation close to the ice front. We do not observe any seasonal variation or significant change in ice flow speed across the rest of the Shackleton system. Given the potential vulnerability of the system to accelerating retreat into the overdeepened, potentially sediment-filled bedrock trough, an improved understanding of the glaciological, oceanographic and geological conditions in the Shackleton system are required to improve the certainty of numerical model predictions, and we identify a number of priorities for future research. With access to these remote coastal regions a major challenge, coordinated internationally collaborative efforts are required to quantify how much the Shackleton region is likely to contribute to sea level rise in the coming centuries.

Generation and fate of basal meltwater during winter, western Greenland Ice Sheet

Abstract. Basal sliding in the ablation zone of the Greenland Ice Sheet is closely associated with water from surface melt introduced to the bed in summer, yet melting of basal ice also generates subglacial water year-round. Assessments of basal melt rely on modeling with results strongly dependent upon assumptions with poor observational constraints. Here we use surface and borehole measurements to investigate the generation and fate of basal meltwater in the ablation zone of Isunnguata Sermia basin, western Greenland. The observational data are used to constrain estimates of the heat and water balances, providing insights into subglacial hydrology during the winter months when surface melt is minimal or nonexistent. Despite relatively slow ice flow speeds during winter, the basal meltwater generation from sliding friction remains manyfold greater than that due to geothermal heat flux. A steady acceleration of ice flow over the winter period at our borehole sites can cause the rate of basal water generation to increase by up to 20 %. Borehole measurements show high but steady basal water pressure rather than monotonically increasing pressure. Ice and groundwater sinks for water do not likely have sufficient capacity to accommodate the meltwater generated in winter. Analysis of basal cavity dynamics suggests that cavity opening associated with flow acceleration likely accommodates only a portion of the basal meltwater, implying that a residual is routed to the terminus through a poorly connected drainage system. A forcing from cavity expansion at high pressure may explain observations of winter acceleration in western Greenland.

Recent acceleration of Denman Glacier (1972–2017), East Antarctica, driven by grounding line retreat and changes in ice tongue configuration

Abstract. After Totten, Denman Glacier is the largest contributor to sea level rise in East Antarctica. Denman's catchment contains an ice volume equivalent to 1.5 m of global sea level and sits in the Aurora Subglacial Basin (ASB). Geological evidence of this basin's sensitivity to past warm periods, combined with recent observations showing that Denman's ice speed is accelerating and its grounding line is retreating along a retrograde slope, has raised the prospect that its contributions to sea level rise could accelerate. In this study, we produce the first long-term (∼50 years) record of past glacier behaviour (ice flow speed, ice tongue structure and calving) and combine these observations with numerical modelling to explore the likely drivers of its recent change. We find a spatially widespread acceleration of the Denman system since the 1970s across both its grounded (17±4 % acceleration; 1972–2017) and floating portions (36±5 % acceleration; 1972–2017). Our numerical modelling experiments show that a combination of grounding line retreat, ice tongue thinning and the unpinning of Denman's ice tongue from a pinning point following its last major calving event are required to simulate an acceleration comparable with observations. Given its bed topography and the geological evidence that Denman Glacier has retreated substantially in the past, its recent grounding line retreat and ice flow acceleration suggest that it could be poised to make a significant contribution to sea level in the near future.

Hydrological and Kinematic Precursors of the 2017 Calving Event at the Petermann Glacier in Greenland Observed from Multi-Source Remote Sensing Data

Both a decrease of sea ice and an increase of surface meltwater, which may induce ice-flow speedup and frontal collapse, have a significant impact on the stability of the floating ice shelf in Greenland. However, detailed dynamic precursors and drivers prior to a fast-calving process remain unclear due to sparse remote sensing observations. Here, we present a comprehensive investigation on hydrological and kinematic precursors before the calving event on 26 July 2017 of Petermann Glacier in northern Greenland, by jointly using remote sensing observations at high-temporal resolution and an ice-flow model. Time series of ice-flow velocity fields during July 2017 were retrieved with Sentinel-2 observations with a sub-weekly sampling interval. The ice-flow speed quickly reached 30 m/d on 26 July (the day before the calving), which is roughly 10 times quicker than the mean glacier velocity. Additionally, a significant decrease in the radar backscatter coefficient of Sentinel-1 images suggests a rapid transformation from landfast sea ice into open water, associated with a decrease in sea ice extent. Additionally, the area of melt ponds on the floating ice tongue began to increase in mid-May, quickly reached a peak at the end of June and lasted for nearly one month until the calving occurred. We used the ice sheet system model to model the spatial-temporal damage and stress on the floating ice, thereby finding an abnormal stress distribution in a cracked region. It is inferred that this calving event may relate to a weakening of the sea ice, shearing of the tributary glacier, and meltwater infiltrating crevasses.

Spatiotemporal dynamic characteristics of typical temperate glaciers in China

China’s temperate glaciers have a relatively warm and humid climate and hydrothermal conditions at low latitudes. Temperate glaciers, however, have larger ablation, higher ice temperatures, relatively fast movement speeds, and a significant sliding process at the bottom. As a result, these glaciers are more significantly affected by climate change. On the basis of topographic maps, aerial photography, and Landsat OLI images, and combined with existing research results, this paper systematically analyzed the temporal and spatial dynamic characteristics of typical temperate glaciers. The results are as follows: (1) From the 1950s to the 1970s, compared with other types of glaciers, temperate glaciers showed strong retreat and ablation trends in terms of area, length, speed, and mass balance. (2) The reduction rates of glacier areas of Kangri Garpo, Dagu Snow Mountain, Yulong Snow Mountain (YSM), and Meili Snow Mountain (MSM) in China’s temperate glacier areas all exceeded 38%, which was far above the national average of 18% from the 1950s to the 2010s. (3) The recent length retreat rates of Azha Glacier, Kangri Garpo, and Mingyong Glacier, MSM, Hailuogou Glacier (HG), Gongga Snow Mountain (GSM), and Baishui River Glacier No. 1 (BRGN1), YSM were above 22 m/a, which was faster than the retreat rates of other regions. (4) Consistent with glacier retreat, temperate glaciers also had a faster ice flow speed. The ice flow velocities of the BGN1, HG, Parlung River Glaciers No. 4 and 94, and Nyainqêntanglha were, respectively, 6.33–30.78 m/a, 41–205 m/a, 15.1–86.3 m/a, and 7.5–18.4 m/a, which was much faster than the velocity of other types of glaciers. (5) Mass loss of temperate glaciers was most dramatic during the observation period (1959–2015). The annual mass balance from eight typical temperate glaciers fluctuated between − 2.48 and 0.44 m w.e., and the annual average change rate of mass balance (− 0.037 m w.e./a) was much higher than that in China (− 0.015 m w.e./a, p < 0.0001) and globally (− 0.013 m w.e./a, p < 0.0001).

Manifestations and mechanisms of the Karakoram glacier Anomaly.

Global-scale glacier shrinkage is one of the most prominent signs of ongoing climatic change. However, important differences in glacier response exist at the regional scale, and evidence has accumulated that one particular region stands out: the Karakoram. In the past two decades, the region has shown balanced to slightly positive glacier budgets, an increase in glacier ice-flow speeds, stable to partially advancing glacier termini, and widespread glacier surge activity. This is in stark contrast to the rest of High Mountain Asia, where glacier retreat and slowdown dominate, and glacier surging is largely absent. Termed the Karakoram Anomaly, recent observations show that the anomalous glacier behaviour partially extends to the nearby Western Kun Lun and Pamir. Several complementary explanations have now been presented for explaining the Anomaly’s deeper causes, but the understanding is far from being complete. Whether the Anomaly will continue to exist in the coming decades remains unclear, but its long-term persistence seems unlikely in light of the considerable warming anticipated by current projections of future climate.

Subglacial roughness of the Greenland Ice Sheet: relationship with contemporary ice velocity and geology

Abstract. The subglacial environment of the Greenland Ice Sheet (GrIS) is poorly constrained both in its bulk properties, for example geology, the presence of sediment, and the presence of water, and interfacial conditions, such as roughness and bed rheology. There is, therefore, limited understanding of how spatially heterogeneous subglacial properties relate to ice-sheet motion. Here, via analysis of 2 decades of radio-echo sounding data, we present a new systematic analysis of subglacial roughness beneath the GrIS. We use two independent methods to quantify subglacial roughness: first, the variability in along-track topography – enabling an assessment of roughness anisotropy from pairs of orthogonal transects aligned perpendicular and parallel to ice flow and, second, from bed-echo scattering – enabling assessment of fine-scale bed characteristics. We establish the spatial distribution of subglacial roughness and quantify its relationship with ice flow speed and direction. Overall, the beds of fast-flowing regions are observed to be rougher than the slow-flowing interior. Topographic roughness exhibits an exponential scaling relationship with ice surface velocity parallel, but not perpendicular, to flow direction in fast-flowing regions, and the degree of anisotropy is correlated with ice surface speed. In many slow-flowing regions both roughness methods indicate spatially coherent regions of smooth beds, which, through combination with analyses of underlying geology, we conclude is likely due to the presence of a hard flat bed. Consequently, the study provides scope for a spatially variable hard- or soft-bed boundary constraint for ice-sheet models.

Seismic Noise Interferometry Reveals Transverse Drainage Configuration Beneath the Surging Bering Glacier

AbstractSubglacial drainage systems are known to critically control ice flows, but their spatial configuration and temporal evolution are poorly constrained due to inaccessibility. Here we report a 12‐year‐long monitoring of the drainage underneath Bering Glacier, Alaska, by correlating ambient noise recorded at two seismic stations on the sides of the glacier. We find that the seismic surface waves traveling across Bering Glacier slowed down by 1–2% during its latest 2008–2011 surge, likely due to the switch of the subglacial drainage from a channelized system to a distributed system. In contrast to current models, the relative amplitude of velocity reductions for Rayleigh and Love waves requires the distributed drainage to be highly anisotropic and aligned perpendicular to the ice flow direction. We infer that the subglacial water flow is mainly through a network of transverse basal crevasses during surges and thus can sustain the high water pressure and ice flow speed.

Glacial lake evolution and glacier–lake interactions in the Poiqu River basin, central Himalaya, 1964–2017

ABSTRACTDespite previous studies, glacier–lake interactions and future lake development in the Poiqu River basin, central Himalaya, are still not well understood. We mapped glacial lakes, glaciers, their frontal positions and ice flow from optical remote sensing data, and calculated glacier surface elevation change from digital terrain models. During 1964–2017, the total glacial-lake area increased by ~110%. Glaciers retreated with an average rate of ~1.4 km2a−1between 1975 and 2015. Based on rapid area expansion (&gt;150%), and information from previous studies, eight lakes were considered to be potentially dangerous glacial lakes. Corresponding lake-terminating glaciers showed an overall retreat of 6.0 ± 1.4 to 26.6 ± 1.1 m a−1and accompanying lake expansion. The regional mean glacier elevation change was −0.39 ± 0.13 m a−1while the glaciers associated with the eight potentially dangerous lakes lowered by −0.71 ± 0.05 m a−1from 1974 to 2017. The mean ice flow speed of these glaciers was ~10 m a−1from 2013 to 2017; about double the mean for the entire study area. Analysis of these data along with climate observations suggests that ice melting and calving processes play the dominant role in driving lake enlargement. Modelling of future lake development shows where new lakes might emerge and existing lakes could expand with projected glacial recession.

Ice flow variations at Polar Record Glacier, East Antarctica

ABSTRACTRelatively little is known about the physical mechanisms that drive the dynamics of the East Antarctic outlet glaciers. Here we conduct a remote-sensing investigation of the Polar Record Glacier (PRG), East Antarctica to analyze its ice flow acceleration, ice front variations and ice surface melting. Ice flow speeds at PRG increased by up to 15% from 2005 to 2015, with substantial interannual fluctuations. The ice velocities also showed seasonal variations, accelerating by up to 9% between September and January. Multiple mechanisms contribute to the observed seasonal variations: the initial acceleration may result from the lost back-stress provided by the sea ice in the austral spring and the later speedup relate to the surface meltwater that leads to weakened ice shelf and shear margins. The sensitivity of the PRG to oceanic forcing is confirmed by comparing the secular ice velocity increases with ocean temperatures. These measurements suggest that the dynamics of East Antarctic ice shelves are sensitive to melt at both the surface and base, at a range of timescales.

Ice front and flow speed variations of marine-terminating outlet glaciers along the coast of Prudhoe Land, northwestern Greenland

ABSTRACTSatellite images were analyzed to measure the frontal positions and ice speeds of 19 marine-terminating outlet glaciers along the coast of Prudhoe Land, northwestern Greenland from 1987 to 2014. All the studied glaciers retreated over the study period at a rate of between 12 and 200 m a−1, with a median (mean) retreat rate of 30 (40) m a−1. The glacier retreat began in the year ~2000, which coincided with an increase in summer mean air temperature from 1.4 to 5.5 °C between 1996 and 2000 in this region. Ice speed near the front of the studied glaciers ranged between 20 and 1740 m a−1 in 2014, and many of them accelerated in the early 2000s. In general, the faster retreat was observed at the glaciers that experienced greater acceleration, as represented by Tracy Glacier, which experienced a retreat of 200 m a−1 and a velocity increase of 930 m a−1 during the study period. A possible interpretation of this observation is that flow acceleration induced dynamic thinning near the termini, resulting in enhanced calving and rapid retreat of the studied glaciers. We hypothesize that atmospheric warming conditions in the late 1990s triggered glacier retreat in northwestern Greenland since 2000.

Warm-based basal sediment entrainment and far-field Pleistocene origin evidenced in central Transantarctic blue ice through stable isotopes and internal structures

ABSTRACTStable isotopes of water (δ18O and δ2H) were measured in the debris-laden ice underlying an Antarctic blue ice moraine, and in adjoining Law Glacier in the central Transantarctic Mountains. Air bubble content and morphology were assessed in shallow ice core samples. Stable isotope measurements plot either on the meteoric waterline or are enriched from it. The data cluster in two groups: the ice underlying the moraine has a δ2H:δ18O slope of 5.35 ± 0.92; ice from adjoining portions of Law Glacier has a slope of 6.69 ± 1.39. This enrichment pattern suggests the moraine's underlying blue ice entrained sediment through refreezing processes acting in an open system. Glaciological conditions favorable to warm-based sediment entrainment occur 30–50 km upstream. Basal melting and refreezing are further evidenced by abundant vapor figures formed from internal melting of the ice crystals. Both the moraine ice and Law Glacier are sufficiently depleted of heavy isotopes that their ice cannot be sourced locally, but instead must be derived from far-field interior regions of the higher polar plateau. Modeled ice flow speeds suggest the ice must be at least 80 ka old, with Law Glacier's ice possibly dating to OIS 5 and moraine ice older still.

An analysis of the ice temperature and velocity along the main flowline of Guliya Ice Cap of Western Kunlun Mountains based on glacier dynamical model

Glacier is an important component of the hydrological cycle at different temporal scales. A major mechanism of glacier formation is ″metamorphosis″, by which snow changes into ice. The process is affected by glacier temperature, velocity and mass balance. In recent decades, glaciers in the Western Kunlun Mountains manifested differently from those in other regions. Understanding the mechanism behind their abnormal response to climate change is pressing in order to predict their future trend. However, few simulation researches of Tibetan Plateau glacier change responding to climate warming have been done due to limitation of the availability of glaciological data in situ measurements. Knowledge of present-day ice temperature and velocity is important in order to determine how fast a glacier may respond to climate change. Guliya Ice Cap is located in the Western Kunlun Mountains which lies in the northwestern part of Tibetan Plateau. It is an extremely continental-type (cold) glacier and is characterized by low temperature and precipitation. Using the few observed data for ice temperature and surface velocity, we apply a two-dimensional higher-order thermomechanical flowline dynamical model to simulate the ice temperature and velocity in steady-state along the main flowline of Guliya Ice Cap. The rare but valuable observational data are used to validate the results of numerical simulation of the model in this study. Owing to the uncertainty of parameters′ values, some diagnostic experiments have been conducted to simulate the present-day spatial distributions of ice velocity and temperature along the main flowline of the Guliya Ice Cap. Based on the sensitivity experiments described above, the appropriate parameters on which the model can best simulate the present-day ice velocity and temperature distributions of Guliya are selected. Generally, modeled and observed ice surface velocities exhibit good agreement under the parameters we adopted. Our results show that significant spatial difference of velocity field of the ice cap exists. The ice flow velocities on the terminus and the ablation zone are lower compared to those in the accumulation area. Similar to other alpine land-terminating glaciers, the mean annual horizontal ice flow speeds of Guliya are relatively low (less than 20 m/a). Several velocity extremum regions inside the ice cap can be observed from the simulation results. Along the streamline profile, the glacier surface velocity is higher than that at the bottom, especially in the glacier accumulation area. Modeled ice temperatures agree well with observed values in a borehole drilled to bottom in 1992. Temperature distribution along the glacier′s streamline profile demonstrates an increasing trend from glacier surface to the bed rock, for example −16.2–−1.65°C at a site of 6200 m above sea level (asl). We extract several vertical profiles of temperature distribution at different altitudes. We find that ice temperature near the glacier terminus remains lower than the melting point, and the temperature gradients are various in different altitudes. This phenomenon indicates that physical property of basal ice frozen to the bed constrains the low amplitude response of the temperature field and the terminus dynamics. We speculate that subsequent simulation researches can possibly consider basal sliding contributes in the future by a sensitivity test of the basal temperature under different scenarios. For a large ice cap with few data in situ observations, it is difficult to construct a suitable glacier model with both reliability and simplicity, nevertheless, we can simulate the main glacier dynamic features with much accuracy. The model suggests that the relatively small change in strain rate and low ice temperature near the glacier terminus may be the reason why the Guliya Ice Cap is the most stable glacier found in China up to now.

ICE VELOCITY MEASUREMENT IN EAST ANTARCTICA FROM 1960s TO 1980s BASED ON ARGON AND LANDSAT IMAGERY

Abstract. Ice flow velocity is a vital parameter for estimating the ice mass balance of glaciers in Antarctica. Especially long time serial observation of the surface velocity is of great significance to assessing the relationship between Antarctic ice materials and global climate change. However, the existing research on Antarctic ice velocity based on remote sensing data since 1970s due to the harsh climate in Antarctica. This paper presents an ice flow velocity estimating method includes image pre-processing, geometric model reconstruction, image ortho-rectification and feature matching by using ARGON images token in 1963 and Landsat images collected form 1973 to 1989.Considering the temporal-spatial distributions of ARGON images and Landsat images in Antarctica, two different methods respectively based on ortho-photos pair and Non-Ortho photos are adopted in this paper. More specifically, when there exist two stereo pairs taken in different time in the glacier region, after being ortho-rectified, the stereo pairs can be used to calculate ice flow velocity based on feature matching method. Otherwise, a parallax decomposition method that separates the effect of the terrain relief from the ice flow motion is applied when there only exists one stereo pair with a certain time interval. With this method, glacier surface velocity is available in the glacier region lacked enough stereo pairs. The methods mentioned above for estimating ice flow velocity are applied in Totten, Amery and Fimbul, etc. in eastern Antarctica. Furthermore, a 1960-80s ice flow speed map in the main glaciers of East Antarctica is produced for the first time.

A New Analytical Method for Estimating Antarctic Ice Flow in the 1960s From Historical Optical Satellite Imagery

Ice flow velocity is used to estimate ice mass changes in glaciers and is a significant indicator of the stability of the Antarctica ice sheet in global change studies. The existing regional Antarctica ice flow speed maps are usually derived from radar or optical satellite observations of modern satellites since the 1970s. This paper presents a new analytical photogrammetric method for estimating Antarctica ice flow velocity fields by using film-based stereo ARGON photographs collected in the 1960s. The key of the proposed innovative method is a parallax decomposition that separates the effect of the terrain relief from the ice flow motion. An innovative implementation strategy is developed by using a framework that involves key techniques of hierarchical stereo image matching, ice flow direction determination, parallax decomposition, and ice flow speed estimation. This method is applied in the Rayner glacier in eastern Antarctica by using two sets of ARGON images with a two-month interval in 1963. The produced digital terrain model and speed map achieved a ground position accuracy of 61 m and a speed accuracy of 70 m $\text{a}^{-1}$ . A comparison with recent products from 2000 to 2010 shows no significant topographic changes in the study area. Furthermore, the speed around the grounding line remained at the same level, while the speed in the ice shelf front decreased by 73 m $\text{a}^{-1}$ . The ice shelf front advanced by approximately 7 km over more than 40 years. Overall, the observation results indicate favorable conditions for the stability of the Rayner glacier-ice shelf system.