Search

The hydrocarbon seepage chronology during deglaciation across the formerly glaciated Barents Sea was established using uranium-thorium (U–Th) dating of seep carbonates. Seep carbonates were sampled with remotely operated vehicles (ROV) from the seafloor at three active hydrocarbon seeps (water depth 156–383 m), located in the north-west (Storfjordrenna), north-central (Storbanken High), and south-west (Loppa High) Barents Sea. Overall, the U–Th dates range from 13.5 to 1.2 thousand years (ka) before present, indicating episodic seep carbonate formation since the late Pleistocene throughout the Holocene. The new U–Th dates indicate protracted post-glacial gas seepage, congruent with previously published seep carbonate ages from the south-west Barents Sea. Gas hydrate dissociation and associated seep carbonate formation occurred at Storfjordrenna between ≈6 and 1.2 ka, and around 13.5 and 6 ka at Storbanken. Early and late Holocene seep carbonate ages from Loppa High support post-glacial seismic activity as potential seepage trigger mechanism.

The Barents Sea, located close to the Arctic Ocean, is a petroleum province featuring an extensive occurrence of gas hydrates and shallow gas in compacted sediments. Glacial erosion and uplift have contributed to the migration of gas originating from deeper rocks to the shallow sediments of this region, resulting in hydrates with higher-order hydrocarbons in addition to methane. This article documents reported gas hydrate indications and major controls on hydrate stability in the Barents Sea.

Submarine groundwater discharge could impact the transport of critical solutes to the ocean. However, its driver(s), significance over geological time scales, and geographical coverage are poorly understood. We characterize a submarine groundwater seep from the continental slope off northern Norway where substantial amount of meteoric water was detected. We reconstruct the seepage history from textural relationships and U‐Th geochronology of authigenic minerals. We demonstrate how glacial‐interglacial dynamics have promoted submarine groundwater circulation more than 100 km offshore and result in high fluxes of critical solutes to the ocean. Such cryosphere‐hydrosphere coupling is likely common in the circum‐Arctic implying that future decay of glaciers and permafrost in a warming Arctic is expected to attenuate such a coupled process and thus decreases the export of critical solutes.

The Harstad Basin is a structural block on the continental shelf of SW Barents Sea where gas hydrates likely occurred below the grounded ice-sheet during the last glaciation and it hosts active gas seepage at numerous seafloor sites. We present an integrated study of fluid flow systems in the Harstad Basin by combining seismic profile interpretations and gas flare mapping data with the geochemical results obtained on seafloor seeping gas and methane-derived carbonate crusts. More than 190 acoustic gas flares were registered in water column, many of them in association with pockmarks and carbonate crust fields. However, weak or absent seepage observed during remotely operated underwater vehicle transects across many pockmarks and crust fields suggests that seepage activity may have decreased since the last deglaciation. In the western Harstad Basin, seeps of microbial methane occur mainly above Tertiary formations that are pinching out below the glacial sediments. High amplitude seismic anomalies suggest the presence of gas pockets at the base of the glacial sediments and within Tertiary deposits. In contrast, gas seeping in the eastern Harstad Basin originates from a biodegraded thermogenic source tentatively connected to the deeply faulted Mesozoic rocks occurring below glacial sediments. This spatial variability in fluid sources is also recorded in the carbon isotope data of seafloor carbonate crusts, with δ13C values typically between −55 and −42‰ and −40 and −20 ‰VPDB for carbonate crusts associated with microbial and thermogenic fluids, respectively. U-Th chronology combined with the stable isotope data suggests that this discrepancy in fluid sources over a distance of about 20 km has been stable since the last glaciation and highlights the significance of regional underlying geology in mediating fluid supply to the seafloor.

Integrated petrography, mineralogy, geochronology and geochemistry of cold seep carbonate crusts and free gas from the Alvheim channel elucidate diagenetic carbonate precipitation and related seepage histories in the central North Sea. Free gas isotope characteristics coupled with carbonate δ13C values as low as − 66‰ VPDB, indicate a predominantly microbial methane source with minor thermogenic contribution. We estimate that ~ 70% of the carbon sequestered into carbonate precipitates was derived from local oxidation of methane. The early stage of crust growth is represented by microcrystalline aragonite and Mg-calcite (10 to 40% mol MgCO3) cementing seafloor sediments consisting of clays, quartz, feldspar, and minor detrital low Mg-calcite and dolomite. Typical association of aragonite cement with coarse-grained detritus may reflect elevated fluid flow and flushing of fine particles prior to cementation close to the seafloor. Middle rare earth element enrichment in early generation microcrystalline cements containing framboidal pyrite indicates diagenetic precipitation within the zone of anaerobic methane oxidation contiguous to iron reduction. The later generation diagenetic phase corresponds to less abundant radial fibrous and botryoidal aragonite which lines cavities developed within the crusts. In contrast to early generation cements, late generation cavity infills have rare earth elements and Y patterns with small negative Ce anomalies similar to seawater, consistent with carbonate precipitation in a more open, seawater dominated system. Aragonite U–Th ages indicate carbonate precipitation between 6.09 and 3.46 kyr BP in the northern part of the channel, whereas in the southern part precipitation occurred between 1.94 and 0.81 kyr BP reflecting regional changes in fluid conduit position.

Gas hydrates stored on continental shelves are susceptible to dissociation triggered by environmental changes. Knowledge of the timescales of gas hydrate dissociation and subsequent methane release are critical in understanding the impact of marine gas hydrates on the ocean–atmosphere system. Here we report a methane efflux chronology from five sites, at depths of 220–400 m, in the southwest Barents and Norwegian seas where grounded ice sheets led to thickening of the gas hydrate stability zone during the last glaciation. The onset of methane release was coincident with deglaciation-induced pressure release and thinning of the hydrate stability zone. Methane efflux continued for 7–10 kyr, tracking hydrate stability changes controlled by relative sea-level rise, bottom water warming and fluid pathway evolution in response to changing stress fields. The protracted nature of seafloor methane emissions probably attenuated the impact of hydrate dissociation on the climate system.

We have estimated the seismic attenuation in gas hydrate and free-gas-bearing sediments from high-resolution P-cable 3D seismic data from the Vestnesa Ridge on the Arctic continental margin of Svalbard. P-cable data have a broad bandwidth (20–300 Hz), which is extremely advantageous in estimating seismic attenuation in a medium. The seismic quality factor (Q), the inverse of seismic attenuation, is estimated from the seismic data set using the centroid frequency shift and spectral ratio (SR) methods. The centroid frequency shift method establishes a relationship between the change in the centroid frequency of an amplitude spectrum and the Q value of a medium. The SR method estimates the Q value of a medium by studying the differential decay of different frequencies. The broad bandwidth and short offset characteristics of the P-cable data set are useful to continuously map the Q for different layers throughout the 3D seismic volume. The centroid frequency shift method is found to be relatively more stable than the SR method. Q values estimated using these two methods are in concordance with each other. The Q data document attenuation anomalies in the layers in the gas hydrate stability zone above the bottom-simulating reflection (BSR) and in the free gas zone below. Changes in the attenuation anomalies correlate with small-scale fault systems in the Vestnesa Ridge suggesting a strong structural control on the distribution of free gas and gas hydrates in the region. We argued that high and spatially limited Q anomalies in the layer above the BSR indicate the presence of gas hydrates in marine sediments in this setting. Hence, our workflow to analyze Q using high-resolution P-cable 3D seismic data with a large bandwidth could be a potential technique to detect and directly map the distribution of gas hydrates in marine sediments.

Author Posting. © Blackwell, 2006. This article is posted here by permission of Blackwell for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 166 (2006): 543–552, doi:10.1111/j.1365-246X.2006.03038.x., The presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments., This work was funded by European Commission under the contract EVK3-CT-2000-00043 (HYDRATECH). JAP and AIB were partly funded by Natural Environment Research Council.

The presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments.