Thinning leads to calving-style changes at Bowdoin Glacier, Greenland

This research has been supported by the Alfred and Flora Spälti Fund and the ETH Zurich Foundation (Sun2ice; grant no. ETH-12 16-2); the Swiss National Science Foundation (grant nos. 200021-153179/1 and PP00P2_183719); the SPI Exploratory Grants 2018 awarded to Fabian Walter and Guillaume Jouvet; the Japanese Ministry of Education, Culture, Sports, Science and Technology through the GRENE Arctic Climate Research Project and the Arctic Challenge for Sustainability (ArCS) project; and by NERC (CALISMO: Calving Laws for Ice Sheet Models; grant no. NE/P011365/1).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.Publisher PDFPeer reviewe

Greenland and Canadian Arctic ice temperature profiles database

Here, we present a compilation of 95 ice temperature profiles from 85 boreholes from the Greenland ice sheet and peripheral ice caps, as well as local ice caps in the Canadian Arctic. Profiles from only 31 boreholes (36 %) were previously available in open-access data repositories. The remaining 54 borehole profiles (64 %) are being made digitally available here for the first time. These newly available profiles, which are associated with pre-2010 boreholes, have been submitted by community members or digitized from published graphics and/or data tables. All 95 profiles are now made available in both absolute (meters) and normalized (0 to 1 ice thickness) depth scales and are accompanied by extensive metadata. These metadata include a transparent description of data provenance. The ice temperature profiles span 70 years, with the earliest profile being from 1950 at Camp VI, West Greenland. To highlight the value of this database in evaluating ice flow simulations, we compare the ice temperature profiles from the Greenland ice sheet with an ice flow simulation by the Parallel Ice Sheet Model (PISM). We find a cold bias in modeled near-surface ice temperatures within the ablation area, a warm bias in modeled basal ice temperatures at inland cold-bedded sites, and an apparent underestimation of deformational heating in high-strain settings. These biases provide process level insight on simulated ice temperatures

PISM-SHMIP: PISM scripts for the Subglacial Hydrology Model Intercomparison Project

Initial release used to prepare results for the SHMIP publication: De Fleurian et al., J. Glaciology, http://doi.org/10.1017/jog.2018.78, 2018

Numerical modelling of the Cordilleran ice sheet

This doctoral dissertation presents a study of the glacial history of the North American Cordillera using numerical ice sheet modelling calibrated against field evidence. This area, characterized by the steep topography of several mountain ranges separated by large inter-montane depressions, was once covered by a large-scale ice mass: the former Cordilleran ice sheet. Because of the irregular topography on which the ice sheet formed, geological studies have often had only local or regional relevance, thus leaving the Cordilleran ice sheet least understood among Pleistocene ice sheets in terms of its extent, volume, and dynamics. Here, I present numerical simulations that allow quantitative reconstructions of the former ice sheet evolution based on approximated physics of glacier flow. These simulations show that the geometry of the Last Glacial Maximum Cordilleran ice sheet was largely controlled by sharp contrasts in regional temperature, precipitation, and daily temperature variability associated with the presence of mountain ranges. However, this maximum stage appears short-lived and out of balance with contemporaneous climate. During most of the simulated last glacial cycle, the North American Cordillera is characterized by an intermediate state of glaciation including isolated glaciers and ice caps covering major mountain ranges, the largest of which is located over the Skeena Mountains. The numerically modelled Cordilleran ice sheet appears in constant imbalance with evolving climate conditions, while the complexity of this transient response transcends that encapsulated in two-dimensional, conceptual models of ice sheet growth and decay. This thesis demonstrates the potential of numerical ice sheet modelling to inform on ice sheet history and former climate conditions over a glacial cycle, given that ice sheet models can be calibrated against field constraints.At the time of the doctoral defense, the following paper was unpublished and had a status as follows: Paper 4: Mansucript.</p

Last-glacial-cycle glacier erosion potential in the Alps

The glacial landscape of the Alps has fascinated generations of explorers, artists, mountaineers, and scientists with its diversity, including erosional features of all scales from high-mountain cirques to steep glacial valleys and large overdeepened basins. Using previous glacier modelling results and empirical inferences of bedrock erosion under modern glaciers, we compute a distribution of potential glacier erosion in the Alps over the last glacial cycle from 120 000 years ago to the present. Despite large uncertainties pertaining to the climate history of the Alps and unconstrained glacier erosion processes, the resulting modelled patterns of glacier erosion include persistent features. The cumulative imprint of the last glacial cycle shows a very strong localization of erosion potential with local maxima at the mouths of major Alpine valleys and some other upstream sections where glaciers are modelled to have flowed with the highest velocity. The potential erosion rates vary significantly through the glacial cycle but show paradoxically little relation to the total glacier volume. Phases of glacier advance and maximum extension see a localization of rapid potential erosion rates at low elevation, while glacier erosion at higher elevation is modelled to date from phases of less extensive glaciation. The modelled erosion rates peak during deglaciation phases, when frontal retreat results in steeper glacier surface slopes, implying that climatic conditions that result in rapid glacier erosion might be quite transient and specific. Our results depict the Alpine glacier erosion landscape as a time-transgressive patchwork, with different parts of the range corresponding to different glaciation stages and time periods.</p

Short-lived ice speed-up and plume water flow captured by a VTOL UAV give insights into subglacial hydrological system of Bowdoin Glacier

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

Last glacial maximum precipitation pattern in the alps inferred from glacier Modelling

During the Last Glacial Maximum (LGM), glaciers in the Alps reached a maximum extent, and broad sections of the foreland were covered by ice. In this study, we simulated the alpine ice cap using a glacier flow model to constrain the prevailing precipitation pattern with a geomorphological reconstruction of ice extent. For this purpose we forced the model using different temperature cooling and precipitation reduction factors. The use of the present-day precipitation pattern led to a systematic overestimation of the ice cover on the northern part of the Alps relative to the southern part. To reproduce the LGM ice cap, a more severe decrease in precipitation in the north than in the south was required. This result supports a southwesterly advection of atmospheric moisture to the Alps, sustained by a southward shift of the North Atlantic storm track during the LGM.ISSN:0016-731