To investigate the mechanisms driving recent changes in outlet glaciers in Antarctica, we measured the glacier front position, flow velocity and surface elevation of five outlet glaciers flowing into Lützow-Holm Bay in East Antarctica. After a steady advance from 2008 to 2015, all the glaciers synchronously retreated by 0.4–6.0 km between 2016 and 2018. The initiation of the retreat coincided with the breakup of land-fast sea ice in Lützow-Holm Bay in 2016, which resulted in the largest sea-ice loss in the region since 1998. Similar flow variations and surface elevation changes were observed near the grounding line of Shirase, Skallen and Telen glaciers. The slowdown in 2011–15 (by 13%) and the speedup in 2016–18 (by 7%) coincided with the respective increase and decrease in surface elevation. Simultaneous retreat and acceleration after the land-fast sea-ice breakup implies that sea ice has a significant influence on glacier dynamics. Thickening/thinning observed near the grounding line was attributed to a reduced/enhanced stretching flow regime during the deceleration/acceleration period. Our results demonstrate that land-fast sea ice affects not only terminus positions, but also the flow speed and ice thickness of the Antarctic glaciers.

On the Kamchatka Peninsula in the Russian Far East, 405 glaciers with an estimated total mass of 49 Gt were reported in the 1970s. These have been retreating at an accelerated rate since the start of the 21st century. Because glacier studies in this region are scarce, ice loss and its influence on sea level rise and regional environments is poorly understood. In this study, we analyzed satellite data to quantify glacier mass change from 2000 to 2016 in six major glacier-covered regions on the peninsula. The mean rate of the glacier mass change over the study period was −0.46 ± 0.01 m w.e. a−1 (total mass change was −4.9 ± 0.1 Gt, −304.2 ± 9.1 Mt a−1), which is slightly lower than other regions in mid-latitude and subarctic zones. The mass loss accelerated from >−0.33 ± 0.02 m w.e. a−1 in the period 2000–2006/2010 to <−1.65 ± 0.12 m w.e. a−1 in 2006/2010–2015/16. The increase in mass loss is attributed to a rise in average decadal summer temperatures observed in the region (+0.68°C from 1987–99 to 2000–13). Moreover, a recent trend in Pacific decadal oscillation suggests future acceleration of mass loss due to a decline in winter precipitation.

As a result of climate warming, glacial meltwater discharge has been increasing in Greenland. During the summers of 2015 and 2016, there were rapid increases in discharge from Qaanaaq Glacier in northwestern Greenland. These discharges resulted in floods that destroyed the road linking the settlement of Qaanaaq to Qaanaaq Airport. Field measurements were performed and a numerical model of glacier runoff was developed to quantify these discharges. The high discharge associated with the 2015 flood, estimated at 9.1 m3 s−1 (hourly mean), resulted from intensive glacier melting due to warm air temperature and strong winds, while the high discharge associated with the 2016 flood resulted from heavy rainfall (90 mm d−1) that led to a peak discharge estimated at 19.9 m3 s−1. The developed model, when used to investigate future glacier runoff under warming conditions, revealed a nonlinear increase in glacial melt with increasing temperature. Additionally, the model forecasted a threefold increase in total summer discharge, owing to a 4 °C rise in temperature. Thus, this study quantified the impact of a changing climate on glacier runoff, which gives insight into future risks of flood hazards along the coast of Greenland.

Under the influence of a rapidly warming climate, abrupt changes have been observed along the coast of Greenland. This commentary is based on a Japanese research project initiated in 2012, in which we examined the recent changes in the coastal environment and their impacts on human society in Qaanaaq, a village in northwestern Greenland. Initially, our research sought to quantify the mass loss of glaciers and its interaction with the ocean in the Qaanaaq region. Over the course of the project in collaboration with local communities, we soon realised that the changes in glaciers and the ocean directly impacted the ~600 residents of Qaanaaq. We observed natural disasters triggered by climate change. Environmental changes are also important for local economy and industry because loss of sea ice may lead to growth in transportation, tourism and mineral resource exploration. In order to share the results of our study with the Qaanaaq community, and to gain understanding of local and traditional knowledge, we organised an annual meeting in the village every summer since 2016. Our experience demonstrates the critical importance of performing a long-term multidisciplinary study, including participation of the local communities to understand the changing environment, and to contribute to a sustainable future in Qaanaaq.

Retreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.

We present a 3-year record of seasonal variations in ice speed and frontal ablation of 10 marine-terminating outlet glaciers along the coast of Prudhoe Land in northwestern Greenland. The glaciers showed seasonal speedup initiated between late May and early June, and terminated between late June and early July. Ice speed subsequently decreased from July to September. The timing of the speedup coincided with the onset of the air temperature rise to above freezing, suggesting an influence of meltwater availability on the glacier dynamics. No clear relationship was found between the speedup and the terminus position or the sea-ice/ice-mélange conditions. These results suggest that the meltwater input to the glacier bed triggered the summer speedup. The excess of summer speed (June–August) over the mean for the rest of the year accounted for 0.5–13% of the annual ice motion. Several glaciers showed seasonal frontal variations, i.e. retreat in summer and advance in winter. This was not due to ice-speed variations, but was driven by seasonal variations in frontal ablation. The results demonstrate the dominant effect of glacier surface melting on the seasonal speedup, and the importance of seasonal speed patterns on longer-term ice motion of marine-terminating outlet glaciers in Greenland.

To study subglacial hydrological condition and its influence on the glacier dynamics, we drilled Johnsons Glacier on Livingston Island in the Antarctic Peninsula region. Subglacial water pressure was recorded in boreholes at two locations over 2 years, accompanied by high-frequency ice-speed measurements during two summer melt seasons. Water pressure showed two different regimes, namely high frequency and large amplitude variations during the melt season (January–April) and small fluctuations near the overburden pressure the rest of the year. Speed-up events were observed several times in each summer measurement period. Ice motion during these events substantially contributed to total glacier motion, for example, fast ice flow over 1 week accounted for ~70% of the total displacement over a 25-day long measurement period. We did not find a clear relationship between subglacial water pressure and ice speed. This was probably because subglacial hydraulic conditions were spatially inhomogeneous and thus our borehole data did not always represent a large-scale subglacial condition. Ice temperature measurements in the boreholes confirmed the existence of a cold ice layer near the glacier surface. Our data provide a basis to better understand the dynamic and hydrological conditions of relatively unstudied glaciers in the Antarctic Peninsula region.

We recorded the ice motion and icequakes on the floating part of Langhovde Glacier in East Antarctica to better understand the dynamic behavior of ice shelves and floating tongues. Diurnal and semi-diurnal variations in ice motion and seismicity were simultaneously observed at all four global navigation satellite system and three seismic stations over 2 weeks. The short-term along-flow ice motion is explained by the elastic response of the glacier to ocean tide-induced hydrostatic stress variations, which decayed at a rate of 0.8 km−1 toward the grounding line. We observed a large number of icequakes during mid-rising and high tides that covered a broad frequency range and formed two major groups of events centered at 10 and 120 Hz, respectively. The hourly occurrence rates were ~500 events h−1, with the observed seismicity consistent with fracture due to floating tongue bending. We also observed minor secondary peaks at high ice speeds, which could reflect surface cracking due to stretching or basal friction. Our observation demonstrates that tidal-modulation was the main factor to fracture the floating tongue of Langhovde Glacier.

We use the full-Stokes model Elmer/Ice to investigate the present dynamics of Bowdoin Glacier, a marine-terminating outlet glacier in northwestern Greenland. Short-term speed variations of the glacier were observed, correlating with air temperature and precipitation, and with the semi-diurnal ocean tides. We use a control inverse method to determine the distribution of basal friction. This reveals that most of the glacier area is characterized by near-plug-flow conditions, while some sticky spots are also identified. We then conduct experiments to test the sensitivity of the glacier flow to basal lubrication and tidal forcing at the calving front. Reduction of the basal drag by 10–40% produces speed-ups that agree approximately with the observed range of speed-ups that result from warm weather and precipitation events. In agreement with the observations, tidal forcing and surface speed near the calving front are found to be in anti-phase (high tide corresponds to low speed, and vice versa). However, the amplitude of the semi-diurnal variability is underpredicted by a factor ~ 3, which is likely related to either inaccuracies in the surface and bedrock topographies or mechanical weakening due to crevassing.

Calving plays a key role in the recent rapid retreat of glaciers around the world. However, many processes related to calving are poorly understood since direct observations are scarce and challenging to obtain. When calving occurs at a glacier front, surface-water waves arise over the ocean or a lake in front of glaciers. To study calving processes from these surface waves, we performed field observations at Glaciar Perito Moreno, Patagonia. We synchronized time-lapse photography and surface waves record to confirm that glacier calving produces distinct waves compared with local noise. A total of 1074 calving events were observed over the course of 39 d. During austral summer, calving occurred twice more frequently than in spring. The cumulative distribution of calving-interevent time interval followed exponential model, implying random occurrence of events in time. We further investigated wave properties and found that source-to-sensor distance can be estimated from wave dispersion within ~20% error. We also found that waves produced by different calving types showed similar spectra in the same frequency range between 0.05–0.2 Hz, and that the amplitude of surface waves increased with the size of calving. This study demonstrates the potential of surface-wave monitoring for understanding calving processes.

Satellite 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.

As a joint contribution of Japan and Sweden to the International Polar Year 2007–09, a field expedition between Syowa and Wasa stations in East Antarctica was carried out in the 2007/08 austral summer season. Along the 2800 km long expedition route, the dielectric permittivity of the upper 1 m snow layer was measured at intervals of approximately 50 km using a snow fork, a parallel-wire transmission-line resonator. More than 2000 measurements were performed under carefully calibrated conditions, mostly in the interior of Antarctica. The permittivity ε′ was a function of snow density as in previous studies on dry snow, but the values were significantly smaller than those reported before. In the light of the dielectric mixture theory, the relatively smaller ε′ obtained in this study can be attributed to the snow structures characteristic in the studied region. Our data suggest that the permittivity of snow in the Antarctic interior is significantly affected by weak bonding between snow grains, which is due to depth-hoar formation in the extremely low-temperature conditions.

To study the past and future evolution of Rhonegletscher, Switzerland, a flowline model was developed to include valley shape effects more accurately than conventional flowband models. In the model, the ice flux at a gridpoint was computed by a two-dimensional ice-flow model applied to the valley cross-section. The results suggested the underestimation of the accumulation area, which seems to be a general problem of flowline modelling arising from the model’s one-dimensional nature. The corrected mass balance was coupled with the equilibrium-line altitude (ELA) change, which was reconstructed for the period 1878–2003 from temperature and precipitation records, to run the model for the past 125 years. The model satisfactorily reproduced both changes in the terminus position and the total ice volume derived from digital elevation models of the surface obtained by analyses of old maps and aerial photographs. This showed the model’s potential to simulate glacier evolution when an accurate mass balance could be determined. The future evolution of Rhonegletscher was evaluated with three mass-balance conditions: the mean for the period 1994–2003, and the most negative (2003) and positive (1978) mass-balance values for the past 50 years. The model predicted volume changes of –18%, –58% and +38% after 50 years for the three conditions, respectively.

To investigate short-term flow-pattern variations of a temperate glacier, longitudinal surface strain was measured with a wire strainmeter in the ablation area of Koryto glacier, Kamchatka, Russia. Strain-rate anomalies were observed in late summer 2000 that were triggered by a water overflow from a moulin near the measurement site followed by the drainage of accumulated water. The strain event started with (compressive) strain rates of >–10–3 d–1 lasting for 6 hours, which then became tensile. Similar strain-rate variations were observed again on the next day. During the event, basal sliding speed measured at the margin in the lower reach of the glacier fluctuated by about ±50% of the daily mean. Smaller and larger sliding speeds corresponded to the compressive and tensile surface strains, respectively. These measurements suggest that the storage and sudden drainage of water caused spatially non-uniform water-pressure fluctuations along the glacier, changing the sliding regime over short time periods.

To investigate the impact of proglacial lake formation on the dynamics and evolution of glaciers, we measured the ice motion of the terminal part of Rhonegletscher, Switzerland, where a lake formed in 2005. In 2009, the flow velocity near the terminus was >20 m a−1. One of the survey stakes tripled its velocity between 2006 and 2007. Since the lake water pressure was consistently close to the ice overburden pressure, it is likely that the high subglacial water pressure enhanced the basal ice motion. The estimated flow velocity due to ice shearing was negligibly small; almost 100% of the horizontal velocity near the terminus was caused by basal sliding. The longitudinal strain rate was large, 0.064 a–1, indicating that much of the glacier thinning was due to ice dynamics. The region of ice flotation adjacent to the lake expanded between 2008 and 2009 as a result of glacier thinning. Accordingly, a huge uplift of the surface was observed in 2009. It is clear from the vertical ice motion as well as visual observations that the marginal part of the glacier began to float. The ice-thinning rate in the studied area from 2008 to 2009 was 3.4 ma–1, larger than previous estimates.

The effects of spatial and temporal variations in basal lubrication on the englacial strain rate and surface velocity distribution are investigated with a numerical ice-flow model. General aspects of the solutions are compared to measurements made on Lauteraargletscher, Switzerland, in 2001, that showed diurnal fluctuations in both surface velocity and englacial vertical strain. We find that spatial gradients in basal lubrication can set up variations in the deviatoric stress field that increases with distance to the bed and has a maximum value near the glacier surface. This stress field produces a significant strain rate near the surface. The temporal evolution of a slippery zone is identified as a possible cause of the observed diurnal variations in the vertical strain rate. Although general aspects of the measurements can be explained in this way, the calculated vertical strain rates are too small, suggesting that the modeled effective viscosity values using Glen’s flow law are too large near the surface.

The spatial pattern of glacier surface melt was measured with a resolution of 20–100m within a region extending 1 km up-glacier from the terminus of Rhonegletscher, Switzerland. The melt rate was monitored from 6 July to 6 September 2009 using 44 ablation stakes. We also measured the surface albedo near the stakes to investigate the importance of this parameter for the melt-rate distribution. The melt rate varied from 32.8 to 71.9 mm w.e. d–1 in the study area. Our measurements suggest that the spatial variation of the melt rate can be explained by (1) shading of the ice surface by neighbouring mountains, (2) surface albedo and (3) effects of microclimate (e.g. radiation from sidewalls) on the surface energy balance. The observed melt-rate distribution was compared to the results of a temperature-index melt model, which takes into account shading of direct solar illumination but not the other two effects. The model reproduces some important features of the field data, but its spatial variations are generally less than the measured values. Our study shows the importance of albedo and other local conditions in the accurate estimation of the small-scale melt-rate distribution.

As part of a deep ice-core drilling project, GPS and ice-radar measurements were made in the vicinity of Aurora Peak, Alaska, USA. Surface and bed elevation maps were obtained over an area of 400 m×700 m, which includes the drilling site. Surface flow velocity was measured at 18 locations in the same area. The results showed that the drilling site was situated on a gentle north–south ridge with saddle-shaped surface geometry. The horizontal ice-flow velocity at the drilling site was <0.5ma–1 and the ice generally flowed eastward and westward along the surface slope. The glacier bed slopes to the east and has more complex topography than the ice surface. The drilling site was situated above a bedrock dip, with an ice thickness of 252±10 m. These data constitute important information for estimating the effect of ice flow to the ice core retrieved from the depth of up to 180.17 m.

During the Japanese-Swedish Antarctic traverse expedition of 2007/08, we measured the surface snow density at 46 locations along the 2800 km long route from Syowa station to Wasa station in East Antarctica. The mean snow density for the upper 1 (or 0.5) m layer varied from 333 to 439 kg m-3 over a region spanning an elevation range of 365-3800 ma.s.l. The density variations were associated with the elevation of the sampling sites; the density decreased as the elevation increased, moving from the coastal region inland. However, the density was relatively insensitive to the change in elevation along the ridge on the Antarctic plateau between Dome F and Kohnen stations. Because surface wind is weak in this region, irrespective of elevation, the wind speed was suggested to play a key role in the near-surface densification. The results of multiple regression performed on the density using meteorological variables were significantly improved by the inclusion of wind speed as a predictor. The regression analysis yielded a linear dependence between the density and the wind speed, with a coefficient of 13.5 kg m-3 (m s-1)-1. This relationship is nearly three times stronger than a value previously computed from a dataset available in Antarctica. Our data indicate that the wind speed is more important to estimates of the surface snow density in Antarctica than has been previously assumed.

During the outburst flood of a glacier-dammed lake on Gornergletscher, Switzerland, in July 2004, the drained lake water triggered anomalous glacier motion. At the onset of the outburst, the ice-flow direction in the vicinity of the lake became closer to the central flowline. When the lake discharge magnitude decreased, the flow direction altered such that the ice moved back to the azimuth of the initial motion. At one of the survey points, where the ice flows parallel to the central flowline, the ice accelerated along the pre-event flow direction followed by a 180° backward motion that lasted over 2 days. These observations indicate the impact of the lake outburst on the subglacial and englacial stress conditions; however, the reversal in the flow direction is difficult to explain by drawing on our current understanding of glacier mechanics. The timing and the timescale of the flow-direction changes suggest that the elastic glacier motion and its rebound played a role under the rapidly changing stress conditions, but the Young’s modulus of ice is too large to cause the observed ice motion. Other processes, including basal separation and subglacial sediment deformation, are discussed as possible mechanisms for the reversal of the ice motion.