In this study, a modelling approach was used to investigate the cause of the diversion of erratic boulders from Mont Blanc and southern Valais by the Valais Glacier to the Solothurn lobe during the ...Last Glacial Maximum (LGM). Using the Parallel Ice Sheet Model, we simulated the ice flow field during the LGM, and analyzed the trajectories taken by erratic boulders from areas with characteristic lithologies. The main difficulty in this exercise laid with the large uncertainties affecting the paleo climate forcing required as input for the surface mass-balance model. In order to mimic the prevailing climate conditions during the LGM, we applied different temperature offsets and regional precipitation corrections to present-day climate data, and selected the parametrizations, which yielded the best match between the modelled ice extent and the geomorphologically-based ice-margin reconstruction. After running a range of simulations with varying parameters, our results showed that only one parametrization allowed boulders to be diverted to the Solothurn lobe during the LGM. This precipitation pattern supports the existing theory of preferential southwesterly advection of moisture to the alps during the LGM, but also indicates strongly enhanced precipitation over the Mont Blanc massif and enhanced cooling over the Jura Mountains.
After more than a century of geological research, the Cordilleran ice sheet of North America remains among the least understood in terms of its former extent, volume, and dynamics. Because of the ...mountainous topography on which the ice sheet formed, geological studies have often had only local or regional relevance and shown such a complexity that ice-sheet-wide spatial reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (62.2–56.9 ka) and 2 (23.2–16.9 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.7 ka).
Subglacial hydrology plays a key role in many glaciological processes, including ice dynamics via the modulation of basal sliding. Owing to the lack of an overarching theory, however, a variety of ...model approximations exist to represent the subglacial drainage system. The Subglacial Hydrology Model Intercomparison Project (SHMIP) provides a set of synthetic experiments to compare existing and future models. We present the results from 13 participating models with a focus on effective pressure and discharge. For many applications (e.g. steady states and annual variations, low input scenarios) a simple model, such as an inefficient-system-only model, a flowline or lumped model, or a porous-layer model provides results comparable to those of more complex models. However, when studying short term (e.g. diurnal) variations of the water pressure, the use of a two-dimensional model incorporating physical representations of both efficient and inefficient drainage systems yields results that are significantly different from those of simpler models and should be preferentially applied. The results also emphasise the role of water storage in the response of water pressure to transient recharge. Finally, we find that the localisation of moulins has a limited impact except in regions of sparse moulin density.
Abstract
Modelling paleo-glacier networks in mountain ranges on the millennial timescales requires ice flow approximations. Hybrid models calculating ice flow by combining vertical shearing (shallow ...ice approximation) and longitudinal stretching (shallow shelf approximation) have been applied to model paleo-glacier networks on steep terrain, yet their validity has not yet been assessed quantitatively. Moreover, hybrid models consistently yield higher ice thicknesses than Last Glacial Maximum geomorphological reconstructions in the European Alps. Here, we compare results based on the hybrid Parallel Ice Sheet Model (PISM) and the Stokes model Elmer/Ice on the Rhine Glacier, a catchment of the former European Alpine Icefield. For PISM, we also test two magnitudes of flux limitation in a scheme that reduces shearing velocities. We find that the flux limitation typically used in PISM yields significantly reduced shearing speeds and increases ice thicknesses by up to 500 m, partly explaining previous overestimations. However, reducing the ice flux limitation allows the hybrid model to minimize this mismatch and captures sliding speeds, ice thicknesses, ice extent and basal temperatures in close agreement with those obtained with the Stokes model.
Ice mass loss from the Greenland ice sheet is the largest single contributor to sea level rise in the 21st century. The mass loss rate has accelerated in recent decades mainly due to thinning and ...retreat of its outlet glaciers. The diverse calving mechanisms responsible for tidewater glacier retreat are not fully understood yet. Since a tidewater glacier’s sensitivity to external forcings depends on its calving style, detailed insight into calving processes is necessary to improve projections of ice sheet mass loss by calving. As tidewater glaciers are mostly thinning, their calving styles are expected to change. Here, we study calving behaviour changes under a thinning regime at Bowdoin Glacier, north-western Greenland, by combining field and remote-sensing data from 2015 to 2019. Previous studies showed that major calving events in 2015 and 2017 were driven by hydro-fracturing and melt-undercutting. New observations from uncrewed aerial vehicle (UAV) imagery and a GPS network installed at the calving front in 2019 suggest ungrounding and buoyant calving have recently occurred as they show (1) increasing tidal modulation of vertical motion compared to previous years, (2) absence of a surface crevasse prior to calving, and (3) uplift and horizontal surface compression prior to calving. Furthermore, an inventory of calving events from 2015 to 2019 based on satellite imagery provides additional support for a change towards buoyant calving since it shows an increasing occurrence of calving events outside of the melt season. The observed change in calving style could lead to a possible retreat of the terminus, which has been stable since 2013. We therefore highlight the need for high-resolution monitoring to detect changing calving styles and numerical models that cover the full spectrum of calving mechanisms to improve projections of ice sheet mass loss by calving.
All around the margins of the Greenland Ice Sheet, marine-terminating glaciers have recently thinned and accelerated. The reduced basal friction has yielded increased flow velocity, while the rate of ...longitudinal stretching has been limited by ice viscosity, which itself critically depends on temperature. However, ice temperature has rarely been measured on such fast-flowing and heavily crevassed glaciers. Here, we present a 3-year record of englacial temperatures obtained 2 (in 2014) to 1 km (in 2017) from the calving front of Bowdoin Glacier (Kangerluarsuup Sermia), a tidewater glacier in northwestern Greenland. Two boreholes separated by 165 (2014) to 197 m (2017) show significant temperature differences averaging 2.07°C on their entire depth. Englacial warming of up to 0.39°C a−1, an order of magnitude above the theoretical rate of heat diffusion and viscous dissipation, indicates a deep and local heat source within the tidewater glacier. We interpret the heat source as latent heat from meltwater refreezing in crevasses reaching to, or near to, the bed of the glacier, whose localization may be controlled by preferential meltwater infiltration in topographic dips between ogives.
The European Alps, the cradle of pioneering glacial studies, are one of the regions where geological markers of past glaciations are most abundant and well-studied. Such conditions make the region ...ideal for testing numerical glacier models based on simplified ice flow physics against field-based reconstructions and vice versa.
The subglacial hydrology of tidewater glaciers is a key but poorly understood component of the complex ice-ocean system, which affects sea level rise. As it is extremely difficult to access the ...interior of a glacier, our knowledge relies mostly on the observation of input variables such as air temperature, and output variables such as the ice flow velocities reflecting the englacial water pressure, and the dynamics of plumes reflecting the discharge of meltwater into the ocean. In this study we use a cost-effective Vertical Take-Off and Landing (VTOL) Unmanned Aerial Vehicle (UAV) to monitor the daily movements of Bowdoin Glacier, north-west Greenland, and the dynamics of its main plume. Using Structure-from-Motion photogrammetry and feature-tracking techniques, we obtained 22 high-resolution ortho-images and 19 velocity fields at the calving front for 12 days in July 2016. Our results show a two-day-long speed-up event (up to 170%) – caused by an increase in buoyant subglacial forces – with a strong spatial variability revealing that enhanced acceleration is an indication of shallow bedrock. Further, we used the Particle Image Velocimetry (PIV) method to analyze water flow from successive UAV images taken while flying over the main plume of the glacier. We found that PIV successfully captures the area of radially diverging flow of the plume, and provides information on spatial and time variability as no other remote sensing technique can. Most interestingly, the active part of the plume features pulsating water jets at the time scale of seconds, and is 1 to 5 times smaller than its visual footprint defined by the iceberg-free area. Combined with an ice flow model or a non-steady plume model, our approach has the potential to generate a novel set of input data to gather information about the depth of the bedrock, the discharge of meltwater, or the subglacial melting rate of tidewater glaciers.
•First use of a VTOL UAV for glacier monitoring•First use of Particle Image Velocimetry (PIV) to infer meltwater plume flow•A major short-lived ice speed-up event of a calving glacier is captured.•PIV results evidence of the unsteady nature of the plume dynamics.•Ice and plume flows give insights into subglacial hydrological system.