Debris-rich Ice Research Articles

Hydraulic fracture characterization of debris-rich ice hole for polar geological drilling based on Peridynamics

Polar geological drilling is crucial to investigating the earth's environment and climate evolution. Hydraulic fracturing occurs around the ice holes when drilling into the complex debris-rich ice before entering the bedrock. Based on the nonlocal theory, the ice-debris structure is considered in the stability of the borehole wall, and a peridynamics model is developed to explore fracture initiation and propagation characteristics around the debris-rich ice hole. The integrated equation describes the mechanical behavior of ice holes, and the continuous and discontinuous spaces are described uniformly around debris-rich ice holes, avoiding the singularity of local continuum mechanics at the interface and crack tip. In this paper, the separation behavior of the interface is simulated between ice and debris. The cracking propagation process, the influence of debris size and shape, and horizontal pressure are analyzed. The results show that at a hydraulic pressure of 4.0 MPa, 8 radial cracks are initiated from the wall of the ice hole containing circular and elliptical debris, with the final crack volume of 0.66% and 1.40%, respectively. The cracks near the hole wall propagate a circumferential direction after turning or branching, leading to peeling in the ice hole containing elliptical debris with inclination angles of 30° and 60°. Hydraulic pressure balances stress distribution around debris-rich holes when horizontal pressure ranges from 1.2 to 2.4 MPa. When drilling into debris-rich ice formations at different debris and horizontal pressures, it is imperative to adjust promptly the density of drilling fluid. This adjustment aims to regulate the hydraulic pressure within the hole, minimize crack distribution, and ensure safe access to bedrock for coring.

A thermomechanical model for frost heave and subglacial frozen fringe

Ice-infiltrated sediment, or frozen fringe, is responsible for phenomena such as frost heave, ice lenses and metres of debris-rich ice under glaciers. Understanding frozen fringes is important as frost heave is responsible for damaging infrastructure at high latitudes and sediment freeze-on at the base of glaciers can modulate subglacial friction, influencing the rate of global sea level rise. Here we describe the thermomechanics of liquid water flow and freezing in ice-saturated sediments, focusing on the conditions relevant for subglacial environments. The force balance that governs the frozen fringe thickness depends on the weight of the overlying material, the thermomolecular force between ice and sediments across liquid premelted films and the water pressure required by Darcy flow. We combine this mechanical model with an enthalpy method that conserves energy across phase change interfaces on a fixed computational grid. The force balance and enthalpy model together determine the evolution of the frozen fringe thickness and our simulations predict frost heave rates and ice lens spacing. Our model accounts for premelting at ice–sediment contacts, partial ice saturation of the pore space, water flow through the fringe, the thermodynamics of the ice–water–sediment interface and vertical force balance. We explicitly account for the formation of ice lenses, regions of pure ice that cleave the fringe at the depth where the interparticle force vanishes. Our model results allow us to predict the thickness of a frozen fringe and the spacing of ice lenses in subaerial and subglacial sediments.

Basal debris of the NEEM ice core, Greenland: a window into sub-ice-sheet geology, basal ice processes and ice-sheet oscillations

AbstractWe 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 (meteoric10Be, ε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. Meteoric10Be 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.

Cosmogenic nuclide dating of two stacked ice masses: Ong Valley, Antarctica

Abstract. We collected a debris-rich ice core from a buried ice mass in Ong Valley, located in the Transantarctic Mountains in Antarctica. We measured cosmogenic nuclide concentrations in quartz obtained from the ice core to determine the age of the buried ice mass and infer the processes responsible for the emplacement of the debris currently overlaying the ice. Such ice masses are valuable archives of paleoclimate proxies; however, the preservation of ice beyond 800 kyr is rare, and therefore much effort has been recently focused on finding ice that is older than 1 Myr. In Ong Valley, the large, buried ice mass has been previously dated at > 1.1 Ma. Here we provide a forward model that predicts the accumulation of the cosmic-ray-produced nuclides 10Be, 21Ne, and 26Al in quartz in the englacial and supraglacial debris and compare the model predictions to measured nuclide concentrations in order to further constrain the age. Large downcore variation in measured cosmogenic nuclide concentrations suggests that the englacial debris is sourced both from subglacially derived material and recycled paleo-surface debris that has experienced surface exposure prior to entrainment. We find that the upper section of the ice core is 2.95 + 0.18 / −0.22 Myr old. The average ice sublimation rate during this time period is 22.86 + 0.10 / −0.09 m Myr−1, and the surface erosion rate of the debris is 0.206 + 0.013 / −0.017 m Myr−1. Burial dating of the recycled paleo-surface debris suggests that the lower section of the ice core belongs to a separate, older ice mass which we estimate to be 4.3–5.1 Myr old. The ages of these two stacked, separate ice masses can be directly related to glacial advances of the Antarctic ice sheet and potentially coincide with two major global glaciations during the early and late Pliocene epoch when global temperatures and CO2 were higher than present. These ancient ice masses represent new opportunities for gathering ancient climate information.

Melting temperature changes during slip across subglacial cavities drive basal mass exchange

AbstractThe importance of glacier sliding has motivated a rich literature describing the thermomechanical interactions between ice, liquid water and bed materials. Early recognition of the gradient in melting temperature across small bed obstacles led to focused studies of regelation. An appreciation for the limits on ice deformation rates downstream of larger obstacles highlighted a role for cavitation, which has subsequently gained prominence in descriptions of subglacial drainage. Here, we show that the changes in melting temperature that accompany changes in normal stress along a sliding ice interface near cavities and other macroscopic drainage elements cause appreciable supercooling and basal mass exchange. This provides the basis of a novel formation mechanism for widely observed laminated debris-rich basal ice layers.

Extensive glaciation in the Erebus Montes region of Mars

Lobate debris apron (LDA) in the Mars' mid-latitudes substantiate extensive glaciation during the Late Amazonian. Detailed investigation of these landforms is imperative because different areas distributed at different latitudes, distinct geologic settings and varied regional topography may have responded to climate in different ways. In this study, mapping of LDA deposits in Erebus Montes region in the flat, low-lying plains of the northern mid-latitudes has been undertaken and examined at a detail previously not attempted to infer new insights on the history of extensive glaciation. LDA deposits show convex-up, steep terminus profiles consistent with typical down-gradient flow characteristics, and integrated flow patterns akin to the glacial landforms reported along the dichotomy boundary. Evidence for a broad piedmont-like lobe, down-gradient flow within a possible oblique-impact crater, and infilled craters, suggest focused localized flow and glaciation. Lobate flows emanating from small alcoves and superposed on the main LDA are not observed, which likely suggests that there is a lack of multi-stage glaciation facilitated by the alcove microclimatic conditions in the region. Linear-curvilinear ridges on the LDA deposits could be the remnant of the internal flow lineations, and are most likely produced by the sublimation of debris-rich ice. Brain-terrain textures, polygonal cracks and ring-mold craters are ubiquitous on the upper surface of LDA deposits, which provides the morphological evidence for the past accumulation of LDM in the region. Hitherto, radar-based investigations do not provide substantial evidence for the presence of extant water ice beneath LDA deposits. We find that the LDA deposits examined here are more consistent with the cold-based glacial behavior - morphological observations supports existence of the sublimation process in the region. We suggest that the derived best-fit age of ~30 Ma for the LDA deposits indicates age of the debris apron that has been mantled and the mapped LDA deposits in our study should be better represented by a broad age range of ~10–100 Ma. Together, our findings add another well-documented case to support the rapidly accumulating evidences for widespread extensive debris-covered glacial landsystems in the northern mid-latitudes of Mars in the Late Amazonian geological history.

Estimation of Ice Thickness and the Features of Subglacial Media Detected by Ground Penetrating Radar at the Baishui River Glacier No. 1 in Mt. Yulong, China

Using ground-penetrating radar (GPR), we measured and estimated the ice thickness of the Baishui River Glacier No. 1 of Yulong Snow Mountain. According to the position of the reflected media from the GPR image, combined with the radar waveform amplitude and polarity change information, the ice thickness and the changing medium position at the bottom of this temperate glacier were identified. Water paths were found in the measured ice, including ice caves and crevasses. A debris-rich ice layer was found at the bottom of the glacier, which produces strong abrasion and ploughing action at the bedrock surface. This results in the formation of different detrital layers stagnated at the ice-bedrock interface and numerous crevasses on the bedrock surface. Based on the obtained ice thickness and differential GPS data, combined with Landsat images, the kriging interpolation method was used to obtain grid data. The average ice thickness was 52.48 m and between 4740 and 4890 m above sea level, with a maximum depth of 92.83 m. The bedrock topography map of this area was drawn using digital elevation model from the Shuttle Radar Topography Mission. The central part of the glacier was characterized by small ice basins with distributed ice steps and ice ridges at the upper and lower parts.

Antarctic subglacial drilling rig: Part II. Ice and Bedrock Electromechanical Drill (IBED)

AbstractA new, modified version of the cable-suspended Ice and Bedrock Electromechanical Drill (IBED) was designed for drilling in firn, ice, debris-rich ice and rock. The upper part of the drill is almost the same for all drill variants and comprises four sections: cable termination, a slip-ring section, an antitorque system and an electronic pressure chamber. The lower part of the IBED comprises an auger core barrel, reamers, a core barrel for ice/debris-ice drilling and a conventional geological single-tube core barrel or custom-made double-tube core barrel. First, the short and full-scale field versions of the IBED were tested at an outdoor testing stand and a testing facility with a 12.5 m-deep ice well. Then, in the 2018–2019 summer season, the IBED was tested in the field at a site ~12 km south of Zhongshan Station, East Antarctica, and a ~6 cm bedrock core was recovered from a 198 m-deep borehole. A total of 18 d was required to penetrate the ice sheet. The retrieved core samples of blue ice, basal ice and bedrock provided valuable information regarding the Earth's paleo-environment.

Magnetic anisotropy and debris-dependent rheological heterogeneity within stratified basal ice

AbstractBasal ice of glaciers and ice sheets frequently contains a well-developed stratification of distinct, semi-continuous, alternating layers of debris-poor and debris-rich ice. Here, the nature and distribution of shear within stratified basal ice are assessed through the anisotropy of magnetic susceptibility (AMS) of samples collected from Matanuska Glacier, Alaska. Generally, the AMS reveals consistent moderate-to-strong fabrics reflecting simple shear in the direction of ice flow; however, AMS is also dependent upon debris content and morphology. While sample anisotropy is statistically similar throughout the sampled section, debris-rich basal ice composed of semi-continuous mm-scale layers (thestratified facies) possesses well-defined triaxial to oblate fabrics reflecting shear in the direction of ice flow, whereas debris-poor ice containing mm-scale star-shaped silt aggregates (thesuspended facies) possesses nearly isotropic fabrics. Thus, deformation within the stratified basal ice appears concentrated in debris-rich layers, likely the result of decreased crystal size and greater availability of unfrozen water associated with high debris content. These results suggest that variations in debris-content over small spatial scales influence ice rheology and deformation in the basal zone.

Freeze-on limits bed strength beneath sliding glaciers

Discharge from sliding outlet glaciers controls uncertainty in projections for future sea level. Remarkably, over 90% of glacial area is subject to gravitational driving stresses below 150 kPa (median ∼70 kPa). Longstanding explanations that appeal to the shear-thinning rheology of ice tend to overpredict driving stresses and are restricted to areas where ice sheets only deform (roughly 50%). Over the more dynamic portions that slide, driving stresses must be balanced by thermo-mechanical interactions that control basal strength. Here we show that median bed strength is comparable to a threshold effective stress set by ice–liquid surface energy and till pore size. Above this threshold, ice infiltrates sediment to produce basal layers of debris-rich ice, even where net melting takes place. We demonstrate that the narrow range of inferred bed strengths can be explained by the mechanical resistance to sliding where roughness is enhanced by heterogeneous freeze-on.

Assessing the efficiency of carbide drill bits and factors influencing their application to debris-rich subglacial ice

When drilling into subglacial bedrock, drill operators commonly encounter basal ice containing high concentrations of rock debris and melt water. As such conditions can easily damage conventional ice drills, researchers have experimented with carbide, diamond, and polycrystalline diamond compact drill bits, with varying degrees of success. In this study, we analyzed the relationship between drilling speed and power consumption for a carbide drill bit penetrating debris-rich ice. We also assessed drill load, rotation speed, and various performance parameters for the cutting element, as well as the physical and mechanical properties of rock and ice, to construct mathematical models. We show that our modeled results are in close agreement with the experimental data, and that both penetration speed and power consumption are positively correlated with drill speed and load. When used in ice with 30% rock content, the maximum penetration speed of the carbide bit is 3.4 mm/s with a power consumption of ≤0.5 kW, making the bit suitable for use with existing electromechanical drills. Our study also provides a guide for further research into cutting heat and equipment design.

A comprehensive interpretation of the NEEM basal ice build-up using a multi-parametric approach

Abstract. Basal ice is a common expression to describe bottom ice layers of glaciers, ice caps and ice sheets in which the ice is primarily conditioned by processes operating at the bed. It is chemically and/or physically distinct from the ice above and can be characterized by a component of basally derived sediments. The study of basal ice properties provides a rare opportunity to improve our understanding of subglacial environments and processes and to refine ice sheet behaviour modelling. Here, we present and discuss the results of water stable isotopes (δ18O and δD), ice fabrics, debris weight/size distribution and gas content of the basal part of the NEEM (North Greenland Eemian Ice Drilling Project) ice core. Below a depth of 2533.85 m, almost 10 m of basal debris-rich material was retrieved from the borehole, and regular occurrence of frozen sediments with only interstitial ice lenses in the bottom 5 m suggest that the ice–bedrock interface was reached. The sequence is composed of an alternation of three visually contrasting types of ice: clear ice with specks (very small amounts) of particulate inclusions, stratified debris-rich layers and ice containing dispersed debris. The use of water stable isotope signatures (δ18O and δD), together with other parameters, allows discrimination between the different types of ice and to unravel the processes involved in their formation and transformation. The basal debris-rich material presents δ18O values [−39.9 ‰; −34.4 ‰] within the range of the above last 300 m of unaltered meteoric ice [−44.9 ‰; −30.6 ‰] spanning a glacial-interglacial range of values. This rules out the hypothesis of a basal ice layer originating from pre-ice sheet ice overridden by the growing ice sheet, as previously suggested e.g. in the case of GRIP (Greenland Ice Core Project). We show that clear basal ice with specks corresponds to altered meteoric glacial ice where some of the climatic signal could have been preserved. However, the stratified debris-rich layers and the ice containing dispersed debris layers respectively express an “open” or “closed” system melting/refreezing signature, somewhat blurred by mixing processes in the upper part of the sequence. Climatic reconstruction is therefore prohibited from these ice types. We propose a first interpretative framework for the build-up of the NEEM basal ice sequence, based on the origin of the various ice types.

Spatial distribution of surface albedo at the Forni Glacier (Stelvio National Park, Central Italian Alps)

The albedo is a fundamental component of the surface energy balance of glaciers, determining the amount of net solar radiation available for melt. However, its temporal and spatial distribution on the glacier surface is often overlooked, and in melt models constant values of ice and snow albedo are assumed for the whole glacier, leading to inaccurate estimations of glacier melt. In this study, we retrieve the albedo on the surface of the Forni Glacier, Stelvio National Park, Central Italian Alps, using four satellite images from the Landsat 7 ETM+ sensors at 30m resolution, and investigate its spatial distribution during the ablation season between 2011 and 2013. The retrieval method includes correction steps for all processes that influence the relationship between the satellite signal and the albedo: i) radiometric calibration, ii) atmospheric correction using the 6S radiative transfer code, iii) correction for local topographic effects, iv) correction for the anisotropy of the reflected radiation over the hemisphere, and v) narrowband to broadband conversion. In spite of the uncertainties introduced by each of these steps (0.02, 0.01, 0.04, 0.05 and 0.01, respectively), the Landsat-retrieved albedo values have a very good correspondence with field measurements from the permanent Automatic Weather Station (AWS) located over the eastern part of the ablation tongue (±0.01 difference). While meteorological conditions, including snowfalls and rainfalls, exert a short lasting influence on the albedo, under conditions typical of the ablation season the albedo values for the ablation tongue are those of debris-rich ice, confirming the ongoing darkening phenomena previously observed in literature.

Experimental study of the drilling process in debris-rich ice

Debris-rich ice is often encountered when drilling into basal ice and rock glaciers. The standard steel bits used for ice core drilling are not suitable because the cutters are very easily broken by rock particles as their hardness and abrasiveness are higher than that of the ice. The tool steel and tungsten carbide inserts are easily damaged in intermixed ice–rock formations. To obtain high-quality core samples in debris-rich ice, it is necessary to find drill bits that can drill ice–rock mixtures with minimal load and acceptable penetration rate and torque. A special testing stand has been designed and constructed to study both standard and custom-made carbide and polycrystalline diamond compact (PDC) drill bits. The results show that both the carbide and the PDC drill bits can drill with high penetration rates in debris-rich ice containing very hard and abrasive granite particles at low drill loads of 500–1200N. When the rock volume content is 30%, the penetration rates are 4.68m/h, 5.9m/h and 11.12m/h for the standard six-tooth carbide drill bit, a PDC bit with a round compact and a PDC bit with a semi-round compact, respectively, under a drill load of 500N with a rotation speed of 100rpm. Within the range of drill loads of 500 to 1200N and rotation speeds of 50 to 200rpm, the maximum torque is no more than 45Nm, and the power consumption is less than 0.8kW. In addition, the temperature changes of the bit cutters caused by their cutting action were also measured. Results of the preliminary tests show that temperature variations increase from 3.67 to 5.96°C when the drill load increases from 450 to 1200N and from 4.17 to 6.21°C when the rotation speed increases from 50 to 200rpm.

Drilling into debris-rich basal ice at the bottom of the NEEM (Greenland) borehole

AbstractAfter the NEEM (Greenland) deep ice-core drilling was declared terminated with respect to developing stratigraphic climate reconstructions, efforts were turned toward collecting basal ice-sheet debris and, if possible, drilling into the bedrock itself. In 2010, several meters of banded debris-rich ice were obtained under normal ice-drilling operations with the NEEM version of the Hans Tausen (HT) drill, but further penetration was obstructed by a rock in the path of the drill head at 2537.36 m. During short campaigns in 2011 and 2012, attempts were made to penetrate further using various reinforced ice cutters mounted on the HT drill head, tailored to cut through rock. These had some success in penetrating coarse material, but produced severely damaged cutters. Additionally a 51 mm diameter diamond cutting tipped rock drill was adapted to fit the NEEM drill. With this device, several additional meters of core containing subglacial sediments, rocks and rock fragments were collected. With these tools 1.39m of additional material were obtained during the 2011 field season, and 7.1 m during 2012. Subglacial water refreezing into the newly formed borehole hindered further penetration, and the bedrock interface was not reached before final closure of the NEEM Camp.

Quantifying the basal conditions of a mountain glacier using a targeted full-waveform inversion: Bench Glacier, Alaska, USA

AbstractGlacier dynamics are inextricably linked to the basal conditions of glaciers. Seismic reflection methods can image the glacier bed under certain conditions. However, where a seismically thin layer of material is present at the bed, traditional analyses may fail to fully characterize bed properties. We use a targeted full-waveform inversion algorithm to quantify the basal-layer parameters of a mountain glacier: thickness (d), P-wave velocity (α) and density (ρ). We simultaneously invert for the seismic quality factor (Q) of the bulk glacier ice. The inversion seeks to minimize the difference between the data and a one-dimensional reflectivity algorithm using a gradient-based search with starting values initialized from a Monte Carlo scheme. We test the inversion algorithm on four basal layer synthetic data traces with 5% added Gaussian noise. The inversion retrieved thin-layer parameters within 10% of synthetic test parameters with the exception of seismic Q. For the seismic dataset from Bench Glacier, Alaska, USA, inversion results indicate a thin basal layer of debris-rich ice within the study area having mean velocity 4000 ± 700 m s–1, density 1900 ± 200 kg m–3 and thickness 6 ± 1.5 m.

Stable isotope and gas properties of two climatically contrasting (Pleistocene and Holocene) ice wedges from Cape Mamontov Klyk, Laptev Sea, northern Siberia

Abstract. This paper presents and discusses the texture, fabric, water stable isotopes (δ18O, δD) and gas properties (total gas content, O2, N2, Ar, CO2, and CH4 mixing ratios) of two climatically contrasted (Holocene vs. Pleistocene) ice wedges (IW-26 and IW-28) from Cape Mamontov Klyk, Laptev Sea, in northern Siberia. The two ice wedges display contrasting structures: one being of relatively "clean" ice and the other showing clean ice at its centre as well as debris-rich ice on both sides (referred to as "ice-sand wedge"). Our multiparametric approach allows discrimination between three different ice facies with specific signatures, suggesting different climatic and environmental conditions of formation and various intensities and nature of biological activity. More specifically, crystallography, total gas content and gas composition reveal variable levels of meltwater infiltration and contrasting contributions from anaerobic and aerobic conditions to the biological signatures. Stable isotope data are drawn on to discuss changes in paleoenvironmental conditions and in the temporal variation of the different moisture sources for the snow feeding into the ice wedges infillings. Our data set also supports the previous assumption that the ice wedge IW-28 was formed in Pleistocene and the ice wedge IW-26 in Holocene times. This study sheds more light on the conditions of ice wedge growth under changing environmental conditions.

Recent dynamics and climate-glacier relationships of maritime mountain glaciers and their implications for palaeoclimatological interpretation of Holocene glacier chronologies - examples from Norway

Controlled moraines are supraglacial debris concentrations that become hummocky moraine upon de-icing and possess clear linearity due to the inheritance of the former pattern of debris-rich folia in the parent ice. Linearity is most striking wherever glacier ice cores still exist but it increasingly deteriorates with progressive melt-out. As a result, moraine linearity has a low preservation potential in deglaciated terrains but hummocky moraine tracts previously interpreted as evidence of areal stagnation may instead record receding polythermal glacier margins in which debris-rich ice was concentrated in frozen toe zones. Recent applications of modern glaciological analogues to palaeoglaciological reconstructions have implied that: (a) controlled moraine development can be ascribed to a specific process (e.g. englacial thrusting or supercooling); and (b) controlled moraine preservation potential is good enough to imply the occurrence of the specific process in former glacier snouts (e.g. ancient polythermal or supercooled snouts). These assumptions are tested using case studies of controlled moraine construction in which a wide range of debris entrainment and debris-rich ice thickening mechanisms are seen to produce the same geomorphic features. Polythermal conditions are crucial to the concentration of supraglacial debris and controlled moraines in glacier snouts via processes that are most effective at the glacier–permafrost interface. End moraines lie on a process–form continuum constrained by basal thermal regime. The morphological expression of englacial structures in controlled moraine ridges is most striking while the moraines retain ice cores, but the final deposits/landforms tend to consist of discontinuous transverse ridges with intervening hummocks, preserving only a weak impression of the former englacial structure. These are arranged in arcuate zones of hummocky moraine up to 2 km wide containing ice-walled lake plains and lying down flow of streamlined landforms produced by warm-based ice. A variety of debris entrainment mechanisms can produce the same geomorphic signature. Spatial and temporal variability in process–form relationships will lead to the sequential development of different types of end moraines during the recession of a glacier or ice sheet margin.

Formation of ice-shelf moraines by accretion of sea water and marine sediment at the southern margin of the McMurdo Ice Shelf, Antarctica

AbstractA combination of ground-penetrating radar surveys, physical sedimentology and ice composition measurements has been used to characterize ice and sediment accreted to the southern margin of the McMurdo Ice Shelf, Antarctica. The radar data and surface observations show that the ice-shelf margin consists of strongly layered debris-rich ice that contains marine sediment and fossils. A modified Rayleigh-based distillation system has been used to model the isotopic fractionation from sea water to ice in a closed system. The model of ice formation is consistent with formation during almost complete freezing of a sea-water reservoir. By contrast, ice on the upstream side of the grounding line has formed during the early stages of freezing in which a small fraction of the sea-water reservoir has frozen. The model results and the presence of delicate, well-preserved marine fossils are interpreted as evidence of anchor ice formation on the seabed, and rafting of glaciomarine sediment into the bottom of the ice shelf. We argue that repeated accretion of sea water and marine sediment has produced a stacked sequence of ice and glaciomarine debris that forms shore-parallel ice-cored moraines.

Debris entrainment by basal freeze‐on and thrusting during the 1995–1998 surge of Kuannersuit Glacier on Disko Island, west Greenland

AbstractKuannersuit Glacier, a valley glacier on Disko Island in west Greenland, experienced a major surge from 1995 to 1998 where the glacier advanced 10·5 km and produced a ∼65 m thick stacked sequence of debris‐rich basal ice and meteoric glacier ice. The aim of this study is to describe the tectonic evolution of large englacial thrusts and the processes of basal ice formation using a multiproxy approach including structural glaciology, stable isotope composition (δ18O and δD), sedimentology and ground‐penetrating radar. We argue that the major debris layers that can be traced in the terminal zone represent englacial thrusts that were formed early during the surge. Thrust overthrow was at least 200–300 m and this lead to a 30 m thick repetition of basal ice at the ice margin. It is assumed that the englacial thrusting was initiated at the transition between warm ice from the interior and the cold snout. The basal debris‐rich ice was mainly formed after the thrusting phase. Two sub‐facies of stratified basal ice have been identified; a lower massive ice facies (SM) composed of frozen diamict enriched with heavy stable isotopes overlain by laminated ice facies (SL) consisting of millimetre thick lamina of alternating debris‐poor and debris‐rich ice. We interpret the stratified basal ice as a continuum formed mainly by freeze‐on processes and localized regelation. First laminated basal ice is formed and as meltwater is depleted more sediment is entrained and finally the glacier freezes to the base and massive diamict is frozen‐on. The increased ability to entrain sediments may partly be associated with higher basal freezing rates enhanced by loss of frictional heat from cessation of fast flow and conductive cooling through a thin heavily crevassed ice during the final phase of the glacier surge. The dispersed basal ice facies (D) was mainly formed by secondary processes where fine‐grained sediment is mobilized in the vein system of ice. Our results have important implications for understanding the significance of basal ice formation and englacial thrusting beneath fast‐flowing glaciers and it provides new information about the development of landforms during a glacier surge. Copyright © 2010 John Wiley & Sons, Ltd.