Widespread methane release from thawing Arctic gas hydrates is a major concern, yet the processes, sources, and fluxes involved remain unconstrained. We present geophysical data documenting a cluster of kilometer-wide craters and mounds from the Barents Sea floor associated with large-scale methane expulsion. Combined with ice sheet/gas hydrate modeling, our results indicate that during glaciation, natural gas migrated from underlying hydrocarbon reservoirs and was sequestered extensively as subglacial gas hydrates. Upon ice sheet retreat, methane from this hydrate reservoir concentrated in massive mounds before being abruptly released to form craters. We propose that these processes were likely widespread across past glaciated petroleum provinces and that they also provide an analog for the potential future destabilization of subglacial gas hydrate reservoirs beneath contemporary ice sheets.

[1] The roughness of a glacier bed has high importance for the estimation of the sliding velocity and can also provide valuable insights into the dynamics and history of ice sheets, depending on scale. Measurement of basal properties in present-day ice sheets is restricted to ground-penetrating radar and seismics, with surveys retrieving relatively coarse data sets. Deglaciated areas, like the Barents Sea, can be surveyed by shipborne 2-D and 3-D seismics and multibeam sonar and provide the possibility of studying the basal roughness of former ice sheets and ice streams with high resolution. Here, for the first time, we quantify the subglacial roughness of the former Barents Sea ice sheet by estimating the spectral roughness of the basal topography. We also make deductions about the past flow directions by investigating how the roughness varies along a 2-D line as the orientation of the line changes. Lastly, we investigate how the estimated basal roughness is affected by the resolution of the basal topography data set by comparing the spectral roughness along a cross section using various sampling intervals. We find that the roughness typically varies on a similar scale as for other previously marine-inundated areas in West Antarctica, with subglacial troughs having very low roughness, consistent with fast ice flow and high rates of basal erosion. The resolution of the data set seems to be of minor importance when comparing roughness indices calculated with a fixed profile length. A strong dependence on track orientation is shown for all wavelengths, with profiles having higher roughness across former flow directions than along them.

Variations in the flow of ice streams and outlet glaciers are a primary control on ice sheet stability, yet comprehensive understanding of the key processes operating at the ice-bed interface remains elusive. Basal resistance is critical, especially sticky spots-localized zones of high basal traction-for maintaining force balance in an otherwise well-lubricated/high-slip subglacial environment. Here we consider the influence of subglacial gas-hydrate formation on ice stream dynamics, and its potential to initiate and maintain sticky spots. Geophysical data document the geologic footprint of a major palaeo-ice-stream that drained the Barents Sea-Fennoscandian ice sheet approximately 20,000 years ago. Our results reveal a ?250 km2 sticky spot that coincided with subsurface shallow gas accumulations, seafloor fluid expulsion and a fault complex associated with deep hydrocarbon reservoirs. We propose that gas migrating from these reservoirs formed hydrates under high-pressure, low-temperature subglacial conditions. The gas hydrate desiccated, stiffened and thereby strengthened the subglacial sediments, promoting high traction-a sticky spot-that regulated ice stream flow. Deep hydrocarbon reservoirs are common beneath past and contemporary glaciated areas, implying that gas-hydrate regulation of subglacial dynamics could be a widespread phenomenon.Peer reviewe

ABSTRACT. Ice-stream dynamics are strongly controlled by processes taking place at the ice/bed interface where subglacial water both lubricates the base and saturates any existing, underlying sediment. Large parts of the former Eurasian ice sheet were underlain by thick sequences of soft, marine sediments and many areas are imprinted with geomorphological features indicative of fast flow and wet basal conditions. Here, we study the effect of subglacial water on past Eurasian ice-sheet dynamics by incorporating a thin-film model of basal water flow into the ice-sheet model SICOPOLIS and use it to better represent flow in temperate areas. The adjunction of subglacial hydrology results in a smaller icesheet building up over time and generally faster ice velocities, which consequently reduces the total area fraction of temperate basal ice and ice streaming areas. Minima in the hydraulic pressure potential, governing water flow, are used as indicators for potential locations of past subglacial lakes and a probability distribution of lake existence is presented based on estimated lake depth and longevity.

Abstract.As ice flows over a subglacial lake, the drop in bed resistance leads to an increase in ice velocities and a draw down of isochrones and cold ice. The ice surface flattens as it adjusts to the lack of resisting forces at the base. The rapid transition in velocity induces changes in ice viscosity and releases deformation energy that can raise the temperature locally. Recent studies of Antarctic subglacial lakes indicate that many lakes experience very fast and possibly episodic drainage, during which the lake size is rapidly reduced as water flows out. Questions that arise are what effect this would have on internal layers within the ice and whether such past drainage events could be inferred from isochrone structures downstream.Here, we study the effect of a subglacial lake on ice dynamics as well as the influence that such short timescale drainage would have on the internal layers of the ice. To this end, we use a full Stokes, polythermal ice flow model. An enthalpy-gradient method is used to account for the evolution of temperature and water content within the ice.We find that a rapid transition between slow-moving ice outside the lake, and full sliding over the lake, can release considerable amounts of deformational energy, with the potential to form a temperate layer at depth in the transition zone. In addition, we provide an explanation for a characteristic surface feature commonly seen at the edges of subglacial lakes, a hummocky surface depression in the transition zone between little to full sliding. We also conclude that rapid changes in the horizontal extent of subglacial lakes and slippery patches, compared to the average ice column velocity, can create a traveling wave at depth within the isochrone structure that transfers downstream with the advection of ice, thus indicating the possibility of detecting past drainage events with ice penetrating radar.

Abstract. As ice flows over a subglacial lake, the drop in bed resistance leads to an increase in ice velocities and a subsequent draw-down of isochrones and cold ice from the surface. The ice surface flattens as it adjusts to the lack of resisting forces at the base. The rapid transition in velocity induces changes in temperature and ice viscosity, releasing deformation energy which raises the temperature locally. Recent studies of Antarctic subglacial lakes indicate that many lakes experience very fast and possibly episodic drainage, during which the lake size is rapidly reduced as water flows out. A question is what effect this would have on internal layers within the ice, and whether such past events could be inferred from isochrone structures downstream. Here, we study the effect of a subglacial lake on the dynamics of a model ice stream as well as the influence that such short timescale drainage would have on the internal layers of the ice. To this end, we use a Full–Stokes, polythermal ice flow model. An enthalpy gradient method is used to account for the evolution of temperature and water content within the ice. We find that the rapid transition between slow-moving ice outside the lake, and full sliding over the lake, releases large amounts of deformational energy, which has the potential to form a temperate layer at depth in the transition zone. In addition, we provide an explanation for a characteristic surface feature, commonly seen at the edges of subglacial lakes, a hummocky surface depression in the transition zone between little to full sliding. We also conclude that rapid changes in lake geometry or basal friction create a travelling wave at depth within the isochrone structure that transfers downstream with the advection of ice, thus indicating the possibility of detecting past events with ice penetrating radar.