Search

The bulk crystal orientation in ice influences the flow of glaciers and ice streams. The ice c-axes fabric is most reliably derived from ice cores. Because these are sparse, the spatial and vertical distribution of the fabric in the Greenland and Antarctic ice sheets is largely unknown. In recent years, methods have been developed to determine fabric characteristics from polarimetric radar measurements. The aim of this paper is to present an improved method to infer the horizontal fabric asymmetry by precisely determining the travel-time difference using co-polarised phase-sensitive radar data. We applied this method to six radar measurements from the East Greenland Ice-core Project (EastGRIP) drill site on Greenland's largest ice stream to give a proof of concept by comparing the results with the horizontal asymmetry of the bulk crystal anisotropy derived from the ice core. This comparison shows an excellent agreement, which is a large improvement compared to previously used methods. Our approach is particularly useful for determining the vertical profile of the fabric asymmetry in higher resolution and over larger depths than was achievable with previous methods, especially in regions with strong asymmetry.

Ice shelves play a key role in the stability of the Antarctic Ice Sheet due to their buttressing effect. A loss of buttressing as a result of increased basal melting or ice shelf disintegration will lead to increased ice discharge. Some ice shelves exhibit channels at the base that are not yet fully understood. In this study, we present in situ melt rates of a channel which is up to 330 m high and located in the southern Filchner Ice Shelf. Maximum observed melt rates are 2 m yr−1. Melt rates inside the channel decrease in the direction of ice flow and turn to freezing ∼55 km downstream of the grounding line. While closer to the grounding line melt rates are higher within the channel than outside, this relationship reverses further downstream. Comparing the modeled evolution of this channel under present-day climate conditions over 250 years with its present geometry reveals a mismatch. Melt rates twice as large as the present-day values are required to fit the observed geometry. In contrast, forcing the model with present-day melt rates results in a closure of the channel, which contradicts observations. The ice shelf experiences strong tidal variability in vertical strain rates at the measured site, and discrete pulses of increased melting occurred throughout the measurement period. The type of melt channel in this study diminishes in height with distance from the grounding line and is hence not a destabilizing factor for ice shelves.

Curvilinear channels on the surface of an ice shelf indicate the presence of large channels at the base. Modelling studies have shown that where these surface expressions intersect the grounding line, they coincide with the likely outflow of subglacial water. An understanding of the initiation and the ice–ocean evolution of the basal channels is required to understand the present behaviour and future dynamics of ice sheets and ice shelves. Here, we present focused active seismic and radar surveys of a basal channel, ∼950 m wide and ∼200 m high, and its upstream continuation beneath Support Force Glacier, which feeds into the Filchner Ice Shelf, West Antarctica. Immediately seaward from the grounding line, below the basal channel, the seismic profiles show an ∼6.75 km long, 3.2 km wide and 200 m thick sedimentary sequence with chaotic to weakly stratified reflections we interpret as a grounding line fan deposited by a subglacial drainage channel directly upstream of the basal channel. Further downstream the seabed has a different character; it consists of harder, stratified consolidated sediments, deposited under different glaciological circumstances, or possibly bedrock. In contrast to the standard perception of a rapid change in ice shelf thickness just downstream of the grounding line, we find a flat topography of the ice shelf base with an almost constant ice thickness gradient along-flow, indicating only little basal melting, but an initial widening of the basal channel, which we ascribe to melting along its flanks. Our findings provide a detailed view of a more complex interaction between the ocean and subglacial hydrology to form basal channels in ice shelves.

Projections of the contribution of the Greenland ice sheet to future sea-level rise include uncertainties primarily due to the imposed climate forcing and the initial state of the ice sheet model. Several state-of-the-art ice flow models are currently being employed on various grid resolutions to estimate future mass changes in the framework of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). Here we investigate the sensitivity to grid resolution of centennial sea-level contributions from the Greenland ice sheet and study the mechanism at play.We employ the finiteelement higher-order Ice-sheet and Sea-level System Model (ISSM) and conduct experiments with four different horizontal resolutions, namely 4, 2, 1 and 0.75 km. We run the simulation based on the ISMIP6 core climate forcing from the MIROC5 global circulation model (GCM) under the highemission Representative Concentration Pathway (RCP) 8.5 scenario and consider both atmospheric and oceanic forcing in full and separate scenarios. Under the full scenarios, finer simulations unveil up to approximately 5% more sea-level rise compared to the coarser resolution. The sensitivity depends on the magnitude of outlet glacier retreat, which is implemented as a series of retreat masks following the ISMIP6 protocol. Without imposed retreat under atmosphereonly forcing, the resolution dependency exhibits an opposite behaviour with approximately 5% more sea-level contribution in the coarser resolution. The sea-level contribution indicates a converging behaviour below a 1 km horizontal resolution. A driving mechanism for differences is the ability to resolve the bedrock topography, which highly controls ice discharge to the ocean. Additionally, thinning and acceleration emerge to propagate further inland in high resolution for many glaciers. A major response mechanism is sliding, with an enhanced feedback on the effective normal pressure at higher resolution leading to a larger increase in sliding speeds under scenarios with outlet glacier retreat., SCOPUS: ar.j, info:eu-repo/semantics/published

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Pine Island Glacier is one of the fastest changing glaciers of the Antarctic Ice Sheet and therefore of scientific interest. The glacier holds enough ice to raise the global sea level significantly (~ 0.5 m) when fully melted. The question addressed by numerous modelling studies of the glacier focuses on whether the observed changes are a start of an uncontrolled and accelerating retreat. The movement of the glacier is, in the fast-flowing areas, dominated by basal motion. In modelling studies the parametrisation of the basal motion is therefore crucial. Inversion methods are commonly applied to reproduce the complex surface flow structure of Pine Island Glacier by using information of the observed surface velocity field to constrain, among other things, basal sliding. We introduce two different approaches of combining a physical parameter, the basal roughness, with basal sliding parametrisations. This way basal sliding is again connected closer to its original formulation. We show that the basal roughness is an important and helpful parameter to consider and that many features of the flow structure can be reproduced with these approaches.

The ongoing disintegration of large ice shelf parts in Antarctica raise the need for a better understanding of the physical processes that trigger critical crack growth in ice shelves. Finite elements in combination with configurational forces facilitate the analysis of single surface fractures in ice under various boundary conditions and material parameters. The principles of linear elastic fracture mechanics are applied to show the strong influence of different depth dependent functions for the density and the Young's modulus on the stress intensity factor KI at the crack tip. Ice, for this purpose, is treated as an elastically compressible solid and the consequences of this choice in comparison to the predominant incompressible approaches are discussed. The computed stress intensity factors KI for dry and water filled cracks are compared to critical values KIc from measurements that can be found in literature.

The accelerated ice flow of ice streams that reach far into the interior of the ice sheets is associated with lubrication of the ice sheet base by basal meltwater. However, the amount of basal melting under the large ice streams – such as the Northeast Greenland Ice Stream (NEGIS) – is largely unknown. In situ measurements of basal melt rates are important from various perspectives as they indicate the heat budget, the hydrological regime and the relative importance of sliding in glacier motion. The few previous estimates of basal melt rates in the NEGIS region were 0.1 m a−1 and more, based on radiostratigraphy methods. These findings raised the question of the heat source, since even an increased geothermal heat flux could not deliver the necessary amount of heat. Here, we present basal melt rates at the recent deep drill site EastGRIP, located in the centre of NEGIS. Within 2 subsequent years, we found basal melt rates of 0.19±0.04 m a−1 that are based on analysis of repeated phase-sensitive radar measurements. In order to quantify the contribution of processes that contribute to melting, we carried out an assessment of the energy balance at the interface and found the subglacial water system to play a key role in facilitating such high melt rates.