Upernavik Isstrøm, the largest contributor to sea-level rise in northwest Greenland, has experienced complex and contrasting ice-flow-speed changes across its five outlets over the last two decades. In this study, we present a detailed remote-sensing analysis of the ice dynamics at Upernavik's outlets from 2000 to 2021 to evaluate the details of these changes. Previous research suggested that the presence or absence of floating ice tongues strongly influences Upernavik's ice dynamics. We use several lines of evidence to document the presence of floating ice tongues, and find that, while several outlets experienced ice-tongue formation and/or loss during the study period, these changes do not explain observed fluctuations in ice-flow velocity. Further exploration of ice-dynamic forcings using a flowline model suggests that changes in basal slipperiness near the terminus have a strong impact on upstream ice dynamics and can explain the velocity variations. Our results suggest that speed fluctuations at Upernavik's outlets may be seasonally and interannually controlled by bed conditions near the terminus, and highlight the need for further research on the influence of basal conditions on complex tidewater glacier dynamics.

Ice-crystal fabric can induce mechanical anisotropy that significantly affects flow, but ice-flow models generally do not include fabric development or its effect upon flow. Here, we incorporate a new spectral expansion of fabric, and more complete description of its evolution, into the ice-flow model Elmer/Ice. This approach allows us to model the effect of both lattice rotation and migration recrystallization on large-scale ice flow. The fabric evolution is coupled to flow using an unapproximated non-linear orthotropic rheology that better describes deformation when the stress and fabric states are misaligned. These improvements are most relevant for simulating dynamically interesting areas, where recrystallization can be important, tuning data are scarce and rapid flow can lead to misalignment between stress and fabric. We validate the model by comparing simulated fabric to ice-core and phase-sensitive radar measurements on a transect across Dome C, East Antarctica. With appropriately tuned rates for recrystallization, the model is able to reproduce observations of fabric. However, these tuned rates differ from those previously derived from laboratory experiments, suggesting a need to better understand how recrystallization acts differently in the laboratory compared to natural settings.

We present new data from the debris-rich basal ice layers of the NEEM ice core (NW Greenland). Using mineralogical observations, SEM imagery, geochemical data from silicates (meteoric 10Be, εNd, 87Sr/86Sr) and organic material (C/N, δ13C), we characterize the source material, succession of previous glaciations and deglaciations and the paleoecological conditions during ice-free episodes. Meteoric 10Be data and grain features indicate that the ice sheet interacted with paleosols and eroded fresh bedrock, leading to mixing in these debris-rich ice layers. Our analysis also identifies four successive stages in NW Greenland: (1) initial preglacial conditions, (2) glacial advance 1, (3) glacial retreat and interglacial conditions and (4) glacial advance 2 (current ice-sheet development). C/N and δ13C data suggest that deglacial environments favored the development of tundra and taiga ecosystems. These two successive glacial fluctuations observed at NEEM are consistent with those identified from the Camp Century core basal sediments over the last 3 Ma. Further inland, GRIP and GISP2 summit sites have remained glaciated more continuously than the western margin, with less intense ice-substratum interactions than those observed at NEEM.

This study provides the first comprehensive reconstruction of the dynamics of Iceberg Glacier, located on western Axel Heiberg Island, and reveals detailed observations of a complete surge for the first time in the Canadian Arctic. Historical aerial photographs, declassified intelligence satellite photographs, optical satellite imagery and synthetic aperture radar data were used to quantify changes in terminus position, ice velocity and glacier thickness since the 1950s. A surge initiated at the terminus in 1981 and terminated in 2003, suggesting a 22-year active phase. High surface velocities, reaching ~2300 m a−1 in 1991, were accompanied by a maximum terminus advance of >7 km and a large transfer of mass down-glacier, causing significant median trunk-wide surface elevation changes attaining >3 ± 1 m a−1. We suggest that the retreat from a pinning point, flotation of the terminus, the removal of sea-ice from the ice front, and an increase in subglacial meltwater availability from relatively high air temperatures in 1981 likely contributed to surge initiation. The ensuing quiescent period has seen a continual decrease in surface flow rates to an average centreline velocity of 11.5 m a−1 in 2020–21, a gradual steepening of the glacier surface and a > 2.5 km terminus retreat.

During the past 20 years, multi-channel radar emerged as a key tool for deciphering an ice sheet's internal architecture. To assign ages to radar reflections and connect them over large areas in the ice sheet, the layer genesis has to be understood on a microphysical scale. Synthetic radar trace modelling based on the dielectric profile of ice cores allows for the assignation of observed physical properties’ variations on the decimetre scale to radar reflectors extending from the coring site to a regional or even whole-ice-sheet scale. In this paper we rely on the available dielectric profiling data of the northern Greenland deep ice cores: NGRIP, NEEM and EGRIP. The three records are well suited for assigning an age model to the stratigraphic radar-mapped layers, and linking up the reflector properties to observations in the cores. Our modelling results show that the internal reflections are mainly due to conductivity changes. Furthermore, we deduce fabric characteristics at the EGRIP drill site from two-way-travel-time differences of along and across-flow polarized radarwave reflections of selected horizons (below 980 m). These indicate in deeper parts of the ice column an across-flow concentrated c-axis fabric.

Tiny samples of ancient atmosphere in air bubbles within ice cores contain argon (Ar), which can be used to reconstruct past temperature changes. At a sufficient depth, the air bubbles are compressed by the overburden pressure under low temperature and transform into air-hydrate crystals. While the oxygen (O2) and nitrogen (N2) molecules have indeed been identified in the air-hydrate crystals with Raman spectroscopy, direct observational knowledge of the distribution of Ar at depth within ice sheet and its enclathration has been lacking. In this study, we applied scanning electron microscopy (SEM) and energy-dispersive X-ray spectroscopy (EDS) to five air-hydrate crystals in the Greenland NEEM ice core, finding them to contain Ar and N. Given that Ar cannot be detected by Raman spectroscopy, the method commonly used for O2 and N2, the SEM-EDS measurement method may become increasingly useful for measuring inert gases in deep ice cores.

Subglacial hydrological systems require innovative technological solutions to access and observe. Wireless sensor platforms can be used to collect and return data, but their performance in deep and fast-moving ice requires quantification. We report experimental results from Cryoegg: a spherical probe that can be deployed into a borehole or moulin and transit through the subglacial hydrological system. The probe measures temperature, pressure and electrical conductivity in situ and returns all data wirelessly via a radio link. We demonstrate Cryoegg's utility in studying englacial channels and moulins, including in situ salt dilution gauging. Cryoegg uses VHF radio to transmit data to a surface receiving array. We demonstrate transmission through up to 1.3 km of cold ice – a significant improvement on the previous design. The wireless transmission uses Wireless M-Bus on 169 MHz; we present a simple radio link budget model for its performance in cold ice and experimentally confirm its validity. Cryoegg has also been tested successfully in temperate ice. The battery capacity should allow measurements to be made every 2 h for more than a year. Future iterations of the radio system will enable Cryoegg to transmit data through up to 2.5 km of ice.

Bulk directional enhancement factors are determined for axisymmetric (girdle and single-maximum) orientation fabrics using a transversely isotropic grain rheology with an orientation-dependent non-linear grain fluidity. Compared to grain fluidities that are simplified as orientation independent, we find that bulk strain-rate enhancements for intermediate-to-strong axisymmetric fabrics can be up to a factor of ten larger, assuming stress homogenization over the polycrystal scale. Our work thus extends previous results based on simple basal slip (Schmid) grain rheologies to the transversely isotropic rheology, which has implications for large-scale anisotropic ice-flow modelling that relies on a transversely isotropic grain rheology. In order to derive bulk enhancement factors for arbitrary evolving fabrics, we expand the c-axis distribution in terms of a spherical harmonic series, which allows the rheology-required structure tensors through order eight to easily be calculated and provides an alternative to current structure-tensor-based modelling.

Ever since the first deep ice cores were drilled, it has been a challenge to determine their original, in-situ orientation. In general, the orientation of an ice core is lost as the drill is free to rotate during transport to the surface. For shallow ice cores, it is usually possible to match the adjacent core breaks, which preserves the orientation of the ice column. However, this method fails for deep ice cores, such as the EastGRIP ice core in Northeast Greenland. We provide a method to reconstruct ice core orientation using visual stratigraphy and borehole geometry. As the EastGRIP ice core is drilled through the Northeast Greenland Ice Stream, we use information about the directional structures to perform a full geographical re-orientation. We compared the core orientation with logging data from core break matching and the pattern of the stereographic projections of the crystals’ c-axis orientations. Both comparisons agree very well with the proposed orientation method. The method works well for 441 out of 451 samples from a depth of 1375–2120 m in the EastGRIP ice core. It can also be applied to other ice cores, providing a better foundation for interpreting physical properties and understanding the flow of ice.

The detection and monitoring of meltwater within firn presents a significant monitoring challenge. We explore the potential of small wireless sensors (ETracer+, ET+) to measure temperature, pressure, electrical conductivity and thus the presence or absence of meltwater within firn, through tests in the dry snow zone at the East Greenland Ice Core Project site. The tested sensor platforms are small, robust and low cost, and communicate data via a VHF radio link to surface receivers. The sensors were deployed in low-temperature firn at the centre and shear margins of an ice stream for 4 weeks, and a ‘bucket experiment’ was used to test the detection of water within otherwise dry firn. The tests showed the ET+ could log subsurface temperatures and transmit the recorded data through up to 150 m dry firn. Two VHF receivers were tested: an autonomous phase-sensitive radio-echo sounder (ApRES) and a WinRadio. The ApRES can combine high-resolution imaging of the firn layers (by radio-echo sounding) with in situ measurements from the sensors, to build up a high spatial and temporal resolution picture of the subsurface. These results indicate that wireless sensors have great potential for long-term monitoring of firn processes.

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.

We conducted a pit study in July 2009 at the NEEM (North Greenland Eemian Ice Drilling) deep ice-coring site in northwest Greenland. To examine the seasonal variations of snow chemistry and characteristics of the drill site, we collected snow/firn samples from the wall of a 2 m deep pit at intervals of 0.03 m and analyzed them for electric conductivity, pH, Cl–, NO3–, SO42–, CH3SO3– (MSA), Na+, K+, Mg2+, Ca2+ and stable isotopes of water (δ18O and δD). Pronounced seasonal variations in the stable isotopes of water were observed, which indicated that the snow had accumulated regularly during the past 4 years. Concentrations of Na+, Cl– and Mg2+, which largely originate from sea salt, peaked in winter to early spring, while Ca2+, which mainly originates from mineral dust, peaked in late winter to spring, slightly later than Na+, Cl– and Mg2+. Concentrations of NO3– showed double peaks, one in summer and the other in winter to spring, whereas those of SO42– peaked in winter to spring. The winter-to-spring concentrations of NO3– and SO42– seem to have been strongly influenced by anthropogenic inputs. Concentrations of MSA showed double peaks, one in spring and the other in late summer to autumn. Our study confirms that the NEEM deep ice core can be absolutely dated to a certain depth by counting annual layers, using the seasonal variations of stable isotopes of water and those of ions. We calculated the annual surface mass balance for the years 2006–08. The mean annual balance was 176 mm w.e., and the balances for winter-to-summer and summer-to-winter halves of the year were 98 and 78 mm, respectively. Snow deposition during the winter-to-summer half of the year was greater than that during the summer-to-winter half by 10–20mm for all three years covered by this study.

Uniaxial compression tests were performed on samples of the Greenland Ice Gore Project (GRIP) deep ice core, both in the field and later in a cold-room laboratory, in order to understand the ice-flow behavior of large ice sheets. Experiments were conducted under conditions of constant strain rate (type A) and constant load (type B). Fifty-four uniaxial-compression test specimens from 1327-2922 m were selected. Each test specimen (25 mm x 25 mm x 90 mm) was prepared with its uniaxial stress axis inclined 45° from the core axis in order to examine the flow behavior of strong single-maximum ice-core samples with basal planes parallel to the horizontal plane of the ice sheet. The ice-flow enhancement factors show a gradual increase with depth down to approximately 2000 m. These results can be interpreted in terms of an increase in the fourth-order Schmid factor. Below 2000 m depth, the flow-enhancement factor increases to about 20-30 with a relatively high variability When the Schmid factor was > 0.46, the enhancement factor obtained was higher than expected from the .-axis concentrations measured. The higher values of flow-enhancement factor were obtained from specimens with a cloudy band structure. It was revealed that cloudy bands affect ice-deformation processes, but the details remain unclear.

Detailed measurements of crystal outlines and fabrics have been performed on 35 000 crystals in fifteen 10 × 20 cm2 vertical thin sections from the North Greenland Icecore Project (NorthGRIP) ice core, evenly distributed in the depth interval 115–880m. The crystals exhibit important changes over this period. As the ice gets older the mean crystal area increases towards a constant value, the shape of the crystals becomes increasingly irregular, and the area distribution of crystals develops from a single log-normal distribution into a bimodal lognormal distribution. The c-axis fabric of the ice shows a smooth development of an increasingly stronger vertical fabric with depth, and the formation of a weak vertical girdle. Already in the younger samples the fabric is rather strongly oriented towards vertical. The fabric and the area of individual crystals are found not to correlate. A simple model, which takes into account the vertical strain of the ice, is applied in an attempt to determine the crystal growth rate at NorthGRIP.

A high-resolution study of c-axis fabrics has been performed on the NorthGRIP deep ice core from North Greenland. A newly developed automatic ice-fabric analyzer has been used to measure c-axis orientations on vertical thin sections at 142 different depths between 100 and 2920 m in the core. Detailed comparison studies show that this new method produces results that are in full accordance with those obtained from conventional manual measurements. Fabric development in ice of the NorthGRIP core can be explained by rotation of c axes, controlled by the prevailing stress systems. In the upper 800 mthe c axes appear to distribute uniformly, but in the depth interval 900–2500m they tend to cluster around a vertical plane. This vertical girdle pattern is strengthened with depth, forming the first clear evidence of this fabric type in a Greenland deep ice core. Such fabric development has previously been observed in the Vostok ice core from East Antarctica, where converging flow is believed to occur. the most likely interpretation of the NorthGRIP girdle fabric is that the c axes are rotating away from a horizontal-tension axis across the main ice divide, which runs north-northwest–south-southeast through the NorthGRIP drilling site. This is supported by information available from surface velocity measurements. Within the girdle fabric, increased concentration of c axes parallel to the core axis is observed with increasing depth, indicating combined effects of vertical compression and horizontal tension on the fabric development. from about 2500 m depth, the girdle-type fabric starts to give way to a strong vertical single-maximum fabric, which persists to 2920m depth, where drilling was terminated 150 mabovebedrock in the year 2000. the single maximum seems to suggest that bed-parallel simple shear is exerting a strong influence on the fabric in the lowest part of the ice sheet down to 2920 m.

From temperature measurements down through the 3001 m deep borehole at the North Greenland Icecore Project (NorthGRIP) drill site, it is now clear that the ice at the base, 3080 m below the surface, is at the pressure-melting point. This is supported by the measurements on the ice core where the annual-layer thicknesses show there is bottom melting at the site and upstream from the borehole. Surface velocity measurements, internal radio-echo layers, borehole and ice-core data are used to constrain a time-dependent flow model simulating flow along the north-northwest-trending ice-ridge flow-line, leading to the NorthGRIP site. Also time-dependent melt rates along the flowline are calculated with a heat-flow model. The results show the geothermal heat flow varies from 50 to 200 mW m–2 along the 100km section of the modeled flowline. The melt rate at the NorthGRIP site is 0.75 cm a–1, but the deep ice in the NorthGRIP core originated 50 km upstream and has experienced melt rates as high as 1.1 cm a–1.

No continuous record from Greenland of the Eemian interglacial period (130–115 ka BP) currently exists. However, a new ice-core drill site has been suggested at 77.449˚ N, 51.056˚Win north-west Greenland (North Eemian or NEEM). Radio-echo sounding images and flow model investigations indicate that an undisturbed Eemian record may be obtained at NEEM. In this work, a two-dimensional ice flow model with time-dependent accumulation rate and ice thickness is used to estimate the location of the Eemian layer at the new drill site. The model is used to simulate the ice flow along the ice ridge leading to the drill site. Unknown flow parameters are found through a Monte Carlo analysis of the flow model constrained by observed isochrones in the ice. The results indicate that the Eemian layer is approximately 60m thick and that its base is located approximately 100m above bedrock.

The temperature profile in the 1200 m deep Dome Summit South (DSS) borehole near the summit of Law Dome, Antarctica, was measured in 1996, 3 years after the termination of the deep drilling.

The temperature profile contains information on past surface temperature over the last 4ka. This temperature history is determined by the use of a Monte Carlo inverse method in which no constraints are placed on the unknown temperature history and no solution is assumed to be unique. The temperature history is obtained from a selection of equally well-fitting solutions by a statistical treatment.

The results show that solutions covering the last 4ka have a well-developed central value, a most likely temperature history. The temperature record has two well-developed minima at AD 1250 and 1850. From 1850 to the present, temperatures have gradually increased by 0.7 K. The reconstructed temperatures are compared with the stable oxygen isotope (δ18O) from the DSS ice core.

In the mid-1990s, excellent results from the GRIP and GISP2 deep drilling projects in Greenland opened up funding for continued ice-coring efforts in Antarctica (EPICA) and Greenland (NorthGRIP). The Glaciology Group of the Niels Bohr Institute, University of Copenhagen, was assigned the task of providing drilling capability for these projects, as it had done for the GRIP project. The group decided to further simplify existing deep drill designs for better reliability and ease of handling. The drill design decided upon was successfully tested on Hans Tausen Ice Cap, Peary Land, Greenland, in 1995. The 5.0m long Hans Tausen (HT) drill was a prototype for the ~11m long EPICA and NorthGRIP versions of the drill which were mechanically identical to the HT drill except for a much longer core barrel and chips chamber. These drills could deliver up to 4m long ice cores after some design improvements had been introduced. The Berkner Island (Antarctica) drill is also an extended HT drill capable of drilling 2 m long cores. The success of the mechanical design of the HT drill is manifested by over 12 km of good-quality ice cores drilled by the HT drill and its derivatives since 1995.

The North Greenland Icecore Project (NorthGRIP) was initiated in 1995 as a joint international programme involving Denmark, Germany, Japan, Belgium, Sweden, Iceland, the U.S.A., France and Switzerland. the main goal was to obtain undisturbed high-resolution information about the Eemian climatic period (115–130 kyr BP). the records from the Greenland Icecore Project (GRIP) and Greenland Ice Sheet Project 2 (GISP2) in central Greenland are different and disturbed down in the ice covering this period. Internal radio-echo sounding layers show that NorthGRIP, placed 325 km north-northwest of GRIP at the Summit of the Greenland ice sheet, is located on a gently sloping ice ridge with very flat bedrock and internal layers found so high that an undisturbed Eemian record is possible. Internal layers much farther above bedrock than their apparent counter parts at GRIP suggest that conditions are favourable for recovery of an undisturbed Eemian record. So far, a 1351 mdeep ice core (NorthGRIP1) and a 3001 mdeep ice core (NorthGRIP 2) have been recovered. the ice thickness is expected to be 3080 m, and the ice temperature at 3001 m is –5.6°C, so we expect basal melting at the bedrock. Most of the Eemian ice will be melted away, leaving only the last part and the transition between the Eem and the Last Glacial Period. At 3001 m the age of the ice is 110 kyr BP and the annual layers are of the order 1 cm.With modern methods the annual layers can be resolved, resulting in detailed information on the decline of the warm Eemian period into the Last Glacial Period.