Description: Engabreen is an outlet glacier of the Svartisen Ice Cap located in Northern Norway. It is a unique glacier due to the Svartisen Subglacial Laboratory which allows direct access to the glacier bed. In this study, we combine both sub- and supraglacial observations with ice-flow modelling in order to investigate conditions at the bed of Engabreen both spatially and temporally. We use the full-Stokes model Elmer/Ice and satellite-based surface-velocity maps from 2010 and 2014 to infer patterns of basal friction. Direct measurements of basal sliding and deformation of lower layers of the ice are used to adjust the ice viscosity and provide essential input to the setup of our model and influence the interpretation of the results. We find a clear seasonal cycle in the subglacial conditions at the higher elevation region of the study area and discuss this in relation to the subglacial hydrological system. Our results also reveal an area with an overdeepening where basal friction is significantly lower than elsewhere on the glacier all year round. We attribute this to either water pooling at the base, or saturated sediments and increased strain heating at this location which softens the ice further.

Description: Increasing ice flow velocities of marine terminating glaciers are often linked to rising ocean temperatures. Unfortunately, direct comparisons between glacier velocity and ocean temperatures are impeded by the fact that few oceanographic datasets span multiple years or contain temperatures at depth. Here, we use an oceanographic dataset collected in Helheim Fjord over several years (described in Straneo et al., 2011, Nat. Geoscience) in both shallow and deep waters. We compare the water temperatures at different depths with ice flow velocities that have been calculated from feature-tracking of LandSAT 7 and 8 images. Our results cover the period 2009–2013 and show both seasonal and inter-annual variability. We find that the velocity of Helheim glacier is likely influenced by the deep ocean water temperatures, namely the influx of warm Atlantic water, whereas water temperature at shallower depths do not have a significant influence on glacier speed. This is in contrast with findings from, for example, Svalbard. Our study demonstrates the need for multiple–year ocean datasets at different depths, if we are to disentangle the complex interactions between glaciers and ocean.

Description: The North East Greenland ice-stream (NEGIS) is the largest active ice-stream on the Greenland ice-sheet and is a crucial contributor to the ice-sheet mass balance. To investigate the ice-stream dynamics and to gain information about the past climate, a deep ice-core is drilled in the upstream part of the NEGIS, termed the East Greenland ice-core project (EastGRIP). Upstream flow effects introduce non-climatic bias in ice-cores and are particularly strong at EastGRIP due to high ice-flow velocities and the location in an ice-stream on the eastern flank of the Greenland ice-sheet. Understanding and ultimately correcting for such effects requires information on the source area and the local atmospheric conditions at the time of ice deposition. We use a two-dimensional Dansgaard-Johnsen model to simulate ice-flow along three approximated flow-lines between the summit of the ice-sheet and EastGRIP. Model parameters are determined using a Monte Carlo inversion by minimizing the misfit between modeled isochrones and isochrones observed in radio-echo-sounding images. We calculate backward-in-time particle trajectories to determine the source area of ice found in the EastGRIP core today and present estimates of surface elevation and past accumulation-rates at the deposition site. The thinning function and accumulated strain obtained from the modeled velocity field provide useful information on the deformation history in the EastGRIP ice. Our results indicate that increased accumulation in the upstream area is predominantly responsible for the constant annual layer thickness observed in the upper part of the ice column at EastGRIP. Inverted model parameters suggest that the imprint of basal melting and sliding is present in large parts along the flow profiles and that most internal ice deformation happens close to the bedrock. The results of this study can act as a basis for applying upstream corrections to a variety of ice-core measurements, and the model parameters can be useful constraints for more sophisticated modeling approaches in the future.

Description: Crystal anisotropy of ice causes slight birefringence for electromagnetic waves. At the same time, the mechanical anisotropy amounts to several orders of magnitude, thus making fabric properties highly-relevant for internal deformation. To date, bulk anisotropy of glaciers and ice sheets can be determined by geophysical methods, such as polarimetric radar, or direct sampling from ice cores. A shortcoming has been so far that changes of bulk anisotropy could mainly be inferred at single point observations, but less so as continuous profiles. Here, we exploit the effect of birefringence caused by bulk anisotropy in co-polarized airborne radar data to determine the horizontal anisotropy across the North-East Greenland Ice Stream. We base our analysis on the fact that birefringence causes a second-order effect on radar amplitudes, which leads to a beat frequency in the low and medium frequency range (O(100 kHz)), which is proportional to the horizontal anisotropy. Complementing our radar analysis with direct fabric and dielectric property observations we can constrain the range of all three fabric eigenvalues as a function of space across and along the ice stream. Finally, we assess the effect of the inferred fabric distribution on the overall ice rheology in the context of ice stream dynamics. Our overall approach has the advantage that it can be applied to co-polarized radar systems, as commonly used in profiling surveys, and does not require dedicated cross-polarized radar set-up. This provides the opportunity to revisit older data, especially from Greenland and Antarctica, to map fabric anisotropy in ice-dynamically interesting regions.

Description: The Northeast Greenland Ice Stream (NEGIS) is the largest active ice stream on the Greenland Ice Sheet (GrIS) and a crucial contributor to the ice-sheet mass balance. To investigate the ice-stream dynamics and to gain information about the past climate, a deep ice core is drilled in the upstream part of the NEGIS, termed the East Greenland Ice-core Project (EastGRIP). Upstream flow can introduce climatic bias into ice cores through the advection of ice deposited under different conditions further upstream. This is particularly true for EastGRIP due to its location inside an ice stream on the eastern flank of the GrIS. Understanding and ultimately correcting for such effects requires information on the atmospheric conditions at the time and location of snow deposition. We use a two-dimensional Dansgaard–Johnsen model to simulate ice flow along three approximated flow lines between the summit of the ice sheet (GRIP) and EastGRIP. Isochrones are traced in radio-echo-sounding images along these flow lines and dated with the GRIP and EastGRIP ice-core chronologies. The observed depth–age relationship constrains the Monte Carlo method which is used to determine unknown model parameters. We calculate backward-in-time particle trajectories to determine the source location of ice found in the EastGRIP ice core and present estimates of surface elevation and past accumulation rates at the deposition site. Our results indicate that increased snow accumulation with increasing upstream distance is predominantly responsible for the constant annual layer thicknesses observed in the upper part of the ice column at EastGRIP, and the inverted model parameters suggest that basal melting and sliding are important factors determining ice flow in the NEGIS. The results of this study form a basis for applying upstream corrections to a variety of ice-core measurements, and the inverted model parameters are useful constraints for more sophisticated modelling approaches in the future.

Description: Anisotropic crystal fabrics in ice sheets develop as a consequence of deformation and hence record information of past ice flow. Simultaneously, the fabric affects the present-day bulk mechanical properties of glacier ice because the susceptibility of ice crystals to deformation is highly anisotropic. This is particularly relevant in dynamic areas such as fast-flowing glaciers and ice streams, where the formation of strong fabrics might play a critical role in facilitating ice flow. Anisotropy is ignored in most state-of-the-art ice sheet models, and while its importance has long been recognized, accounting for fabric evolution and its impact on the ice viscosity has only recently become feasible. Both the application of such models to ice streams and their verification through in-situ observations are still rare. Ice cores provide direct and detailed information on the crystal fabric, but the logistical cost, technical challenges, particularly in fast-flowing ice and shear margins, difficulty in reconstructing the absolute orientation of the core, and their limitation of being a point measurement, make ice cores impractical for a spatially extensive evaluation of the fabric type. Indirect geophysical methods applied from or above the ice surface create the link between the small scale of laboratory experiments and ice–core observations to the large-scale coverage required for ice flow models and the complete understanding of ice stream dynamics. Here, we present a comprehensive analysis of the distribution of the ice fabric in the upstream part of the North-East Greenland Ice Stream (NEGIS). Our results are based on a combination of methods applied to extensive airborne and ground-based radar surveys, ice- and firn-core observations, and numerical ice-flow modelling. They show that in the onset region of NEGIS and around the EGRIP ice core drilling site, the fabric is horizontally strongly anisotropic, forming a horizontal girdle perpendicular to the ice flow, while the horizontal anisotropy reduces quickly over distances of less than five ice thicknesses outside of the ice stream’s shear margins. Downstream of the drill site, the fabric develops into a more vertically symmetric configuration on a time scale of around 2 ka, the first observation of this kind. Our study shows how ice-core based fabric observations, geophysical surveys and ice-flow modelling complement each other to obtain a more comprehensive picture of the spatially strongly varying fabric.

Description: Anisotropic crystal fabrics in ice sheets develop as a consequence of deformation and hence record information of past ice flow. Simultaneously, the fabric affects the present-day bulk mechanical properties of glacier ice because the susceptibility of ice crystals to deformation is highly anisotropic. This is particularly relevant in dynamic areas such as fast-flowing glaciers and ice streams, where the formation of strong fabrics might play a critical role in facilitating ice flow. This fact is ignored in most state-of-the-art ice sheet models, and while their importance has been recognized years ago, accounting for fabrics evolution and their impact on the ice viscosity has only recently become feasible. Both, the application of such models in ice streams as well as their verification through in-situ observations are, however, still rare. We present an extensive dataset of fabric anisotropy derived from radar data recorded in the onset region of the Northeast Greenland Ice Stream by air-borne and ground-based systems. Our methods yield the horizontal anisotropy and are based on travel time anisotropy and splitting as well as birefringence-induced power modulation of radar signals. They complement each other and show good agreement. We compare these in-situ observations with the results obtained from a fabric-evolution model employed along flow tubes in the ice stream onset to discuss the fabric in light of past flow history and its significance for the current flow mechanics of the ice stream.

Description: The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation.