Glacial Velocity Research Articles

The Remotely and Directly Obtained Results of Glaciological Studies on King George Island: A Review

Climate warming has become indisputable, and it is now crucial to increase our understanding of both the mechanisms and consequences of climate change. The Antarctic region is subjected to substantial changes, the trends of which have been recognized for several decades. In the South Shetland Islands, the most visible effect of climate change is progressive deglaciation. The following review focuses on past glaciological studies conducted on King George Island (KGI). The results of collected cryosphere element observations are discussed herein in a comprehensive manner. Our analysis showed that there is a lack of temporal as well as spatial continuity for studies on the basic mass balance parameters on the entire KGI ice dome and only Bellingshausen Dome has a relatively long history of data collection. The methodologies of past work, which have improved over time, are also discussed. When studying the glacier front fluctuations, the authors most frequently use a 1956 aerial photography as reference ice coverage. This was the case for seven papers, while other sources are seldomly mentioned. In other papers as many as 41 other sources were used, and therefore comparison to photos taken up to 60 years later can give misleading trends, as small glaciers may have both advanced and retreated in that time. In the case of glacial velocities there is also an apparent lack of consistency, as different glaciers were indicated as the fastest on KGI. Only Lange, Anna, Crystal, Eldred, and eastern part of Usher glaciers were determined by more than one author as the fastest. Additionally, there are gaps in the KGI Ground Penetrating Radar (GPR) survey area, which includes three ice domes: the Warszawa Icefield, the Krakow Icefield, and eastern part of King George Island. Ideas for further work on the topic are also suggested, allowing for easier access to data and thus contributing to a better understanding of glacier development mechanisms.

Interannual modulation of seasonal glacial velocity variations in the Eastern Karakoram detected by ALOS-1/2 data – ERRATUM

Interannual modulation of seasonal glacial velocity variations in the Eastern Karakoram detected by ALOS-1/2 data – ERRATUM

Recent Surge Behavior of Walsh Glacier Revealed by Remote Sensing Data.

Many surge-type glaciers are present on the St. Elias Mountains, but a detailed study on the surge behavior of the glaciers is still missing. In this study, we used remote sensing data to reveal detailed glacier surge behavior, focusing on the recent surge at Walsh Glacier, which was reported to have surged once in the 1960s. Glacial velocities were derived using a cross-correlation algorithm, and changes in the medial moraines were interpreted based on Landsat images. The digital elevation model (DEM) difference method was applied to Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEMs to evaluate the surface elevation of the glacier. The results showed that the surge initiated near the conjunction of the eastern and northern branches, and then quickly spread downward. The surge period was almost three years, with an active phase of less than two years. The advancing speed of the surge front was much large than the maximum ice velocity of ≈14 m/d observed during the active phase. Summer speed-ups and a winter speed-up in ice velocity were observed from velocity data, with the speed-ups being more obvious during the active phase. Changes in the glacier velocity and the medial moraines suggested that the eastern branch was more affected by the surge. The DEM differencing results showed that the receiving zone thickened up to about 140 m, and the upstream reservoir zone became thinner. These surge behaviors, as characterized by remote sensing data, gave us more detailed insights into the surge dynamics of Walsh Glacier.

Fluctuation analysis in the dynamic characteristics of continental glacier based on Full-Stokes model

Ice thickness has a great influence on glacial movement and ablation. Over the course of the change in thickness, area and external climate, the dynamic process of how glaciers change and whether a glacier’s changes in melting tend to be stable or irregular is a problem that needs to be studied in depth. In our study, the changes in the dynamic process of the No. 8 Glacier in Hei Valley (H8) under the conditions of different thicknesses in 1969 and 2009 were simulated based on the Full-Stokes code Elmer/Ice (http://www.csc.fi/elmer/). The results were as follows: (1) The thickness reduction in glaciers would lead to a decrease in ice surface tension and basal pressure and friction at the bottom, and the resulting extensional and compressional flow played an important role in the variations in glacial velocity. (2) The force at the bottom of the glacier was key to maintaining the overall stress balance, and the glaciers that often melted and collapsed in bedrock were more easily destroyed by the overall force balance and increased change rate of glacial thaw. (3) Temperature changes at different altitudes affected the ice viscous force. The closer the ice surface temperature was to the melting point, the greater the influence of temperature changes on the ice viscous force and ice surface velocity. Finally, we used the RCP 4.8 and 8.5 climate models to simulate the changes in H8 over the next 40 years. The results showed that with some decreases in ice surface compression and tension, the gravity component changes caused by local topography begin to control the ice flow movement on the surface of glacier, and melting of the glacial surface will appear as an irregular change. The simulation results further confirmed that the fluctuation in glacial dynamic characteristics could be attributed to the change in the gravity component caused by ablation.

Monitoring and Analyzing Mountain Glacier Surface Movement Using SAR Data and a Terrestrial Laser Scanner: A Case Study of the Himalayas North Slope Glacier Area

The offset tracking technique based on synthetic aperture radar (SAR) image intensity information can estimate glacier displacement even when glacier velocities are high and the time interval between images is long, allowing for the broad use of this technique in glacier velocity monitoring. Terrestrial laser scanners, a non-contact measuring system, can measure the velocity of a glacier even if there are no control points arranged on a glacier. In this study, six COSMO-SkyMed images acquired between 31 July and 22 December 2016 were used to obtain the glacial movements of five glaciers on the northern slope of the central Himalayas using the offset tracking approach. During the period of image acquirement, a terrestrial laser scanner was used, and point clouds of two periods in a small area at the terminus of the Pingcuoliesa Glacier were obtained. By selecting three fixed areas of the point clouds that have similar shapes across two periods, the displacements of the centers of gravity of the selected areas were calculated by using contrast analyses of feature points. Although the overall low-density point clouds data indicate that the glacial surfaces have low albedos relative to the wavelength of the terrestrial laser scanner and the effect of its application is therefore influenced in this research, the registration accuracy of 0.0023 m/d in the non-glacial areas of the scanner’s measurements is acceptable, considering the magnitude of 0.072 m/d of the minimum glacial velocity measured by the scanner. The displacements from the point clouds broadly agree with the results of the offset tracking technique in the same area, which provides further evidence of the reliability of the measurements of the SAR data in addition to the analyses of the root mean squared error of the velocity residuals in non-glacial areas. The analysis of the movement of five glaciers in the study area revealed the dynamic behavior of these glacial surfaces across five periods. G089972E28213N Glacier, Pingcuoliesa Glacier and Shimo Glacier show increasing surface movement velocities from the terminus end to the upper part with elevations of 1500 m, 4500 m, and 6400 m, respectively. The maximum velocities on the glacial surface profiles were 31.69 cm/d, 62.40 cm/d, and 42.00 cm/d, respectively. In contrast, the maximum velocity of Shie Glacier, 50.60 cm/d, was observed at the glacier’s terminus. For each period, Glacier G090138E28210N exhibited similar velocity values across the surface profile, with a maximum velocity of 39.70 cm/d. The maximum velocities of G089972E28213N Glacier, Pingcuoliesa Glacier, and Shie Glacier occur in the areas where the topography is steepest. In general, glacial surface velocities are higher in the summer than in the winter in this region. With the assistance of a terrestrial laser scanner with optimized wavelengths or other proper ground-based remote sensing instruments, the offset tracking technique based on high-resolution satellite SAR data should provide more reliable and detailed information for local and even single glacial surface displacement monitoring.

Joint geodetic and seismic analysis of surface crevassing near a seasonal glacier-dammed lake at Gornergletscher, Switzerland

ABSTRACTSeasonal lake Gornersee forms at the confluence of Gornergletscher and Grenzgletscher, Switzerland, and experiences outburst floods annually in midsummer. To study the interplay between lake drainage, glacier movement and crevasse activity, high-frequency seismometers and GPS receivers were deployed in networks near Gornersee during the summer ablation seasons of 2004, 2006 and 2007. We use a Rayleigh wave coherence method to locate 3289, 7939 and 4087 icequakes, respectively, primarily along well-defined surface crevasses. We calculate two-dimensional strains from triads of GPS stations and find mean differential strain rates of ~300 × 10−6 d−1 with diurnal variations up to 800 × 10−6 d−1. Crevasse icequake activity and glacial velocity are highest during early season, then decrease as meltwater channels erode and subglacial water pressure decreases. Glacial response to Gornersee drainage varied year-to-year, with icequake activity promoted at some crevasses and inhibited at others, suggesting syn-drainage icequakes may be indicative of local drainage patterns and small-scale features of the stress field. Diurnal pulses in icequake activity exhibit peak activity at different times of day in different locations, coincident with a southeast-to-northwest trending concentrated shear zone near the Gornergletscher–Grenzgletscher confluence, likely due to differences in the timing of peak strain rate in these regions.

Geodetic measurements reveal short-term changes of glacial mass near Jakobshavn Isbræ (Greenland) from 2007 to 2017

The Global Positioning System (GPS) and Gravity Recovery and Climate Experiment (GRACE) provide important geodetic datasets to study glacial mass change. Applying the multichannel singular spectral analysis to the GPS-measured vertical and horizonal crustal displacement and GRACE-derived vertical displacement near Jakobshavn Isbræ (JI) in western Greenland from 2007 to 2017, we reconstruct the short-term loading displacements due to ice mass changes. Both the vertical and east displacements show strong seasonal variability. They also reveal three episodes of transient displacements: downward and eastward motion from late 2007 to around 2010, sustained upward and westward motion from 2010 to early 2013, and downward and eastward motion till late 2016. We also forward model the seasonal and transient displacements caused by surface mass balance (SMB) and glacier dynamics. Our model agrees well with the geodetic observations and provides quantitative insights into the contribution from SMB and ice dynamics to the ice mass changes. We find that SMB is the dominant contributor to the seasonal and transient displacements at three out of four GPS sites (AASI, ILUL, and QEQE). While, at the fourth GPS site (KAGA) that is closest to the glacier, the contributions to the transient displacements from SMB and glacier dynamics are comparable. The forward modeling also suggests that the dynamic mass change in the JI catchment underwent strong seasonal variations and these variations correlated more with the seasonal retreat and advance of the calving front than with the changes of glacial velocities. Our altimetry results reveal that the frontal portion of JI catchment lost 34 Gt in 2012 and this loss of ice declined to only 11 Gt in 2016 due to widespread thickening along the main flowline.

Interannual modulation of seasonal glacial velocity variations in the Eastern Karakoram detected by ALOS-1/2 data

ABSTRACTUnlike in most other regions, Karakoram glaciers are either stable or advancing, a phenomenon known as the Karakoram anomaly. Despite studies of glacier surges and the derivation of surface velocity maps, the spatiotemporal variability of glacier dynamics still remains poorly understood, particularly in the Eastern Karakoram Range. We use Advanced Land Observing Satellite/the Phased Array type L-band Synthetic Aperture Radar (ALOS/PALSAR)-1/2 data from 2007 to 2011 and 2014 to 2015 to examine detailed surface velocity patterns of the Siachen, Baltoro, Kundos, Singkhu and Gasherbrum Glaciers. The first three glaciers show considerable velocity variability (20–350 m a−1), with clear seasonal patterns. Although all glaciers, except for Baltoro, flow slowest in 2015, the velocity structures are individual and vary in space and time. In Gasherbrum Glacier, peak surge-phase velocities are seasonally modulated, with maxima in summers 2006 and 2007, suggesting surface melt plays an important role in maintaining the active phase. Given the relatively close proximity of these glaciers, we assume that surface melt timing and rates are comparable. We therefore argue that the observed spatiotemporal and interannual velocity patterns are determined by local and internal mechanisms, including englacial and subglacial hydrology, thermal processes and tributary configuration of each individual glacier.

Utility of 222Rn as a passive tracer of subglacial distributed system drainage

Water flow beneath the Greenland Ice Sheet (GrIS) has been shown to include slow-inefficient (distributed) and fast-efficient (channelized) drainage systems, in response to meltwater delivery to the bed via both moulins and surface lake drainage. This partitioning between channelized and distributed drainage systems is difficult to quantify yet it plays an important role in bulk meltwater chemistry and glacial velocity, and thus subglacial erosion. Radon-222, which is continuously produced via the decay of 226Ra, accumulates in meltwater that has interacted with rock and sediment. Hence, elevated concentrations of 222Rn should be indicative of meltwater that has flowed through a distributed drainage system network. In the spring and summer of 2011 and 2012, we made hourly 222Rn measurements in the proglacial river of a large outlet glacier of the GrIS (Leverett Glacier, SW Greenland). Radon-222 activities were highest in the early melt season (10–15 dpm L−1), decreasing by a factor of 2–5 (3–5 dpm L−1) following the onset of widespread surface melt. Using a 222Rn mass balance model, we estimate that, on average, greater than 90% of the river 222Rn was sourced from distributed system meltwater. The distributed system 222Rn flux varied on diurnal, weekly, and seasonal time scales with highest fluxes generally occurring on the falling limb of the hydrograph and during expansion of the channelized drainage system. Using laboratory based estimates of distributed system 222Rn, the distributed system water flux generally ranged between 1–5% of the total proglacial river discharge for both seasons. This study provides a promising new method for hydrograph separation in glacial watersheds and for estimating the timing and magnitude of distributed system fluxes expelled at ice sheet margins.

Fast draining lakes on the Greenland Ice Sheet

The rapid drainage of supraglacial lakes around the ablation zone of the Greenland Ice Sheet forms an important link between water at the surface and the ice sheet base, allowing surface meltwater to reach the bed and hence increase glacial velocity. The conduits formed by lake drainages may remain open during the remainder of the melt season providing a pathway for further meltwater to reach the base. We investigated the drainage behavior of lakes from all regions of the Greenland Ice Sheet for the period 2005–2009. We mapped the evolution of 2600 lakes from 3704 MODIS images detecting a mean of 263 drainage events per year, of which 61% occurred in the south‐west region. Only 1% of lake drainages occurred in the rapidly thinning south‐east region. Our results show marked differences between the hydrology of the different regions of the ice sheet, with few lake drainages occurring in the regions where the highest dynamic mass loss is occurring. In the south‐west and north‐east, lake drainages are common and could impact glacier dynamics; in the south‐east they are rare and are thus unlikely to do so.

A DORIS determination of the absolute velocities of the Sorsdal and Mellor glaciers in Antarctica

The Lambert–Amery System is the largest glacier–ice shelf system in East Antarctica, draining a significant portion of the ice sheet. Variation in ice sheet discharge from Antarctica or Greenland has an impact on the rate of change in global mean sea level; which is a manifestation of climate change. In conjunction with a measure of ice thickness change, ice sheet discharge can be monitored by determining the absolute velocities of these glaciers. In order to demonstrate the capability of the DORIS system to determine glacier velocities, Geoscience Australia undertook a Pilot Project under the auspices of the International DORIS Service. A DORIS beacon was deployed on the Sorsdal (November 2001 – January 2002 and November 2003 – January 2004) and Mellor (December 2002 – January 2003) glaciers. The DORIS data, transmitted from the autonomously operating ground beacon for each satellite pass, were stored in the receiver on-board the satellite and later downlinked to the DORIS control centres for processing. This paper describes the campaigns that were conducted at the Sorsdal and Mellor glaciers, the data processing standards for modelling the Doppler measurements, precise orbit determination of the satellites using the data from the globally distributed DORIS network, tracking station position and reference frame modelling, the point positioning mode employed for determining the position and velocities of the transmitting beacon antennas located on the glaciers and provides the velocity estimates that have been determined from the analysis of these tracking data. For the Sorsdal 2001/2002 campaign, using SPOT-4 data only, the measured effective horizontal ice motion was estimated to be 30 ± 0.4 cm/day (azimuth of N246°E ± 0.1°). The inferred velocities for the Sorsdal 2003/2004 campaign, using SPOT-4 and SPOT-5 data, was 5.7 ± 0.8 cm/day (azimuth of N264°E ± 7.5°) for the first eight days and 11.4 ± 1.4 cm/day (azimuth of N241°E ± 1.5°) for the subsequent 21 days. There was a noted decrease in the inferred velocities between the beginning and the end of the observing period. A sub-division of the latter 21 day observing period into three segments showed a decrease in 2D velocity from 18.3 ± 0.7 to 11.2 ± 0.7 cm/day and then to 7.4 ± 0.9 cm/day for the first, second and third segments, respectively. In comparison, a GPS derived velocity over the time-span of the 2001/2002 Sorsdal campaign gave a mean ice flow rate of 31 cm/day. The GPS velocity was derived from two daily position estimates 65 days apart. The DORIS determination from 26 days of continuous SPOT-4 and SPOT-5 data compared well with the GPS derived velocity. For the 2002/2003 Mellor glacier campaign, using SPOT-4 and SPOT-5 data, the estimated average ice velocity was 104 ± 25 cm/day (azimuth of N33°E ± 0.1°); which compared well with an InSAR derived velocity of between 110 and 137 cm/day. The point positioning technique as implemented in this study was further validated and assessed by replicating the computational process to determine the position and velocity of the permanent International DORIS Service site at Terre Adélie, Antarctica. Through these experiments, it has been successfully demonstrated that the DORIS system is capable of determining the velocities of glaciers with an accuracy of a few cm/day over a period of several weeks; operating in remote regions in an autonomous mode. With an increasing number of DORIS-equipped satellites and multiple daily passes, it has the potential to measure glacial velocities at a high temporal resolution (sub-daily).

The morphology of ice streams A, B, and C, west Antarctica, and their environs

Airborne radar soundings of the ice sheet surface made in 1984–1985, together with elevations measured by oversnow traverses between 1957 and 1964 have been used to produce a surface elevation map of ice streams A, B, and C and much of the region around them. The flight positions and elevations were tied to many ground control stations whose positions in three dimensions were fixed by satellite tracking. Ice thicknesses were also measured along the flight lines but have not yet been mapped. The surfaces of active ice streams A and B exhibit a longitudinal ridgetrough topography of uncertain origin. Prominent surface valleys are associated with most, but not all, of their marginal shear zones. There is a deep subglacial trough beneath the grid northeastern side of ice stream A that connects to the subglacial trough beneath Reedy Glacier. The boundary between ice streams Bl and B2, the two tributaries of ice stream B, persists as a suture zone for 250 km downstream from their merger. Ice stream Bl overlies a relatively deep subglacial trough that continued under the corresponding part of ice stream B until it disappears close to the grounding line. Between ice streams Bl and B2 there is a complex zone containing several regions of undisturbed ice separated by bands of disturbed ice, which suggest to us that “rafts” of ice are being incorporated into the ice streams. The previously reported negative mass balance for the drainage system of ice stream B may arise from the expansion of the ice stream into the interior of the West Antarctic ice sheet. Crevasse bands at some places on ridge AB, together with high glacial velocities reported elsewhere, also suggest this expansion. Close to the Ross Ice Shelf, ice stream B widens into a nearly flat area where the mean surface slope is only 3.5 × 10−4, and the surface elevations are only a few tens of meters above hydrostatic equilibrium. Inactive ice stream C differs from active ice streams A and B in surficial and basal characteristics. No elongated ridges and troughs are observed; instead, the ice stream surface exhibits several terraces, including some maxima and minima in elevation. Radar sounding reveals areas where basal echoes are strong and steady, indicating subglacial water, alternating with areas of weaker echoes with short fading lengths. Calculations indicate that the strong reflections require a water layer at least several centimeters to several meters thick, depending on the electrical conductivity of the water. Due to reversals in the surface slope, the hydrologic potential gradient changes sign several times along ice stream C, which implies that subglacial water should collect in some locations.