Marine‐terminating outlet glaciers of the Greenland Ice Sheet make significant contributions to global sea level rise, yet the conditions that facilitate their fast flow remain poorly constrained owing to a paucity of data. We drilled and instrumented seven boreholes on Store Glacier, Greenland, to monitor subglacial water pressure, temperature, electrical conductivity, and turbidity along with englacial ice temperature and deformation. These observations were supplemented by surface velocity and meteorological measurements to gain insight into the conditions and mechanisms of fast glacier flow. Located 30 km from the calving front, each borehole drained rapidly on attaining ∼600 m depth indicating a direct connection with an active subglacial hydrological system. Persistently high subglacial water pressures indicate low effective pressure (180–280 kPa), with small‐amplitude variations correlated with notable peaks in surface velocity driven by the diurnal melt cycle and longer periods of melt and rainfall. The englacial deformation profile determined from borehole tilt measurements indicates that 63–71% of total ice motion occurred at the bed, with the remaining 29–37% predominantly attributed to enhanced deformation in the lowermost 50–100 m of the ice column. We interpret this lowermost 100 m to be formed of warmer, pre‐Holocene ice overlying a thin (0–8 m) layer of temperate basal ice. Our observations are consistent with a spatially extensive and persistently inefficient subglacial drainage system that we hypothesize comprises drainage both at the ice‐sediment interface and through subglacial sediments. This configuration has similarities to that interpreted beneath dynamically analogous Antarctic ice streams, Alaskan tidewater glaciers, and glaciers in surge.

Outlet glaciers of the Greenland Ice Sheet transport ice from the interior to the ocean and contribute directly to sea level rise because discharge and ablation often exceed the accumulation. To develop a better understanding of these fast‐flowing glaciers, we investigate the basal conditions of Store Glacier, a large outlet glacier flowing into Uummannaq Fjord in west Greenland. We use two crossing seismic profiles acquired near the centerline, 30 km upstream of the calving front, to interpret the physical nature of the ice and bed. We identify one notably englacial and two notably subglacial seismic reflections on both profiles. The englacial reflection represents a change in crystal orientation fabric, interpreted to be the Holocene‐Wisconsin transition. From Amplitude‐Versus‐Angle (AVA) analysis we infer that the deepest ∼80 m of ice of the parallel‐flow profile below this reflection is anisotropic with an enhancement of simple shear of ∼2. The ice is underlain by ∼45 m of unconsolidated sediments, below which there is a strong reflection caused by the transition to consolidated sediments. In the across‐flow profile subglacial properties vary over small scale and the polarity of the ice‐bed reflection switches from positive to negative. We interpret these as patches of different basal slipperiness associated with variable amounts of water. Our results illustrate variability in basal properties, and hence ice‐bed coupling, at a spatial scale of ∼100 m, highlighting the need for direct observations of the bed to improve the basal boundary conditions in ice‐dynamic models.

ABSTRACT. The phase-sensitive radio-echo sounder (pRES) is a powerful new instrument that can measure the depth of internal layers and the glacier bed to millimetre accuracy. We use a stationary 16-antenna pRES array on Store Glacier in West Greenland to measure the three-dimensional orientation of dipping internal reflectors, extending the capabilities of pRES beyond conventional depth sounding. This novel technique portrays the effectiveness of pRES in deriving the orientation of dipping internal layers that may complement profiles obtained through other geophysical surveying methods. Deriving ice vertical strain rates from changes in layer depth as measured by a sequence of pRES observations assumes that the internal reflections come from vertically beneath the antenna. By revealing the orientation of internal reflectors and the potential deviation from nadir of their associated reflections, the use of an antenna array can correct this assumption. While the array configuration was able to resolve the geometry of englacial layers, the same configuration could not be used to accurately image the glacier bed. Here, we use simulations of the performance of different array geometries to identify configurations that can be tailored to study different types of basal geometry for future deployments.