Search

The continental margin off Prins Karls Forland, western Svalbard, is characterized by widespread natural gas seepage into the water column at and upslope of the gas hydrate stability zone. We deployed an ocean bottom seismometer integrated into the MASOX (Monitoring Arctic Seafloor-Ocean Exchange) automated seabed observatory at the pinch-out of this zone at 389 m water depth to investigate passive seismicity over a continuous 297 day period from 13 October 2010. An automated triggering algorithm was applied to detect over 220,000 short duration events (SDEs) defined as having a duration of less than 1 s. The analysis reveals two different types of SDEs, each with a distinctive characteristic seismic signature. We infer that the first type consists of vocal signals generated by moving mammals, likely finback whales. The second type corresponds to signals with a source within a few hundred meters of the seismometer, either due east or west, that vary on short (∼tens of days) and seasonal time scales. Based on evidence of prevalent seafloor seepage and subseafloor gas accumulations, we hypothesize that the second type of SDEs is related to subseafloor fluid migration and gas seepage. Furthermore, we postulate that the observed temporal variations in microseismicity are driven by transient fluid release and due to the dynamics of thermally forced, seasonal gas hydrate decomposition. Our analysis presents a novel technique for monitoring the duration, intensity, and periodicity of fluid migration and seepage at the seabed and can help elucidate the environmental controls on gas hydrate decomposition and release.

Submarine landslides occur along continental margins at depths that often intersect the gas hydrate stability zone, prompting suggestions that slope stability may be affected by perturbations that arise from changes in hydrate stability. Here we develop a numerical model to identify the conditions under which the destabilization of hydrates results in slope failure. Specifically, we focus on high-saturation hydrate anomalies at fine-grained to coarse-grained stratigraphic boundaries that can transmit bridging stresses that decrease the effective stress at sediment contacts and disrupt normal sediment consolidation. We evaluate slope stability before and after hydrate destabilization. Hydrate anomalies act to significantly increase the overall slope stability due to large increases in effective cohesion. However, when hydrate anomalies destabilize there is a loss of cohesion and increase in effective stress that causes the sediment grains to rapidly consolidate and generate pore pressures that can either trigger immediate slope failure or weaken the surrounding sediment until the pore pressure diffuses away. In cases where failure does not occur, the sediment can remain weakened for months. In cases where failure does occur, we quantify landslide dynamics using a rate and state frictional model and find that landslides can display either slow or dynamic (i.e., catastrophic) motion depending on the rate-dependent properties, size of the stress perturbation, and the size of the slip patch relative to a critical nucleation length scale. Our results illustrate the fundamental mechanisms through which the destabilization of gas hydrates can pose a significant geohazard.

The Vestnesa Ridge comprises a > 100 km long sediment drift located between the western continental slope of Svalbard and the Arctic mid-ocean ridges. It hosts a deep-water (>1000 m) gas hydrate and associated seafloor seepage system. Near-seafloor headspace gas compositions and its methane carbon isotopic signature along the ridge indicate a predominance of thermogenic gas sources feeding the system. Prediction of the base of the gas hydrate stability zone for theoretical pressure and temperature conditions and measured gas compositions, results in an unusual underestimation of the observed bottom simulating reflector (BSR) depth. The BSR is up to 60 m deeper than predicted for pure methane and measured gas compositions with > 99 % methane. Models for measured gas compositions with > 4% higher order hydrocarbons result in a better BSR approximation. However, the BSR remains > 20 m deeper than predicted in a region without active seepage. A BSR deeper than predicted is primarily explained by unexpected spatial variations in the geothermal gradient and by larger amounts of thermogenic gas at the base of the gas hydrate stability zone. Hydrates containing higher order hydrocarbons form at greater depths and higher temperatures and contribute with larger amounts of carbons than pure methane hydrates. In thermogenic provinces, this may imply a significant upward revision (up to 50 % in the case of Vestnesa Ridge) of the amount of carbon in gas hydrates.

Our study presents a basin-scale 3D modeling solution, quantifying and exploring gas hydrate accumulations in the marine environment around the Green Canyon (GC955) area, Gulf of Mexico. It is the first modeling study that considers the full complexity of gas hydrate formation in a natural geological system. Overall, it comprises a comprehensive basin re-construction, accounting for depositional and transient thermal history of the basin, source rock maturation, petroleum components generation, expulsion and migration, salt tectonics and associated multi-stage fault development. The resulting 3D gas hydrate distribution in the Green Canyon area is consistent with independent borehole observations. An important mechanism identified in this study and leading to high gas hydrate saturation (> 80 vol. %) at the base of the gas hydrate stability zone (GHSZ), is the recycling of gas hydrate and free gas enhanced by high Neogene sedimentation rates in the region. Our model predicts the rapid development of secondary intra-salt mini-basins situated on top of the allochthonous salt deposits which leads to significant sediment subsidence and an ensuing dislocation of the lower GHSZ boundary. Consequently, large amounts of gas hydrates located in the deepest parts of the basin dissociate and the released free methane gas migrates upwards to recharge the GHSZ. In total, we have predicted the gas hydrate budget for the Green Canyon area that amounts to ∼3,256 Mt of gas hydrate which is equivalent to ∼340 Mt of carbon (∼7 x 1011 m3 of CH4 at STP conditions), and consists mostly of biogenic hydrates.

The goal of this study is to computationally determine the potential distribution patterns of diffusion-driven methane hydrate accumulations in coarse-grained marine sediments. Diffusion of dissolved methane in marine gas hydrate systems has been proposed as a potential transport mechanism through which large concentrations of hydrate can preferentially accumulate in coarse-grained sediments over geologic time. Using one-dimensional compositional reservoir simulations, we examine hydrate distribution patterns at the scale of individual sand layers (1-20 m thick) that are deposited between microbially active fine-grained material buried through the gas hydrate stability zone (GHSZ). We then extrapolate to two-dimensional and basin-scale three-dimensional simulations, where we model dipping sands and multilayered systems. We find that properties of a sand layer including pore size distribution, layer thickness, dip, and proximity to other layers in multilayered systems all exert control on diffusive methane fluxes toward and within a sand, which in turn impact the distribution of hydrate throughout a sand unit. In all of these simulations, we incorporate data on physical properties and sand layer geometries from the Terrebonne Basin gas hydrate system in the Gulf of Mexico. We demonstrate that diffusion can generate high hydrate saturations (upward of 90%) at the edges of thin sands at shallow depths within the GHSZ, but that it is ineffective at producing high hydrate saturations throughout thick (greater than 10 m) sands buried deep within the GHSZ. Furthermore, we find that hydrate in fine-grained material can preserve high hydrate saturations in nearby thin sands with burial.

Submarine groundwater discharge (SGD) is a large-scale, buoyancy-driven, offshore flow of terrestrial groundwater. If SGD occurs within the permafrost-bearing sediments of the circum-Arctic shelf, such fluid circulation may transport large amounts of dissolved methane to the circum-Arctic shelf, aiding the formation of permafrost-associated gas hydrate. We investigate the feasibility of this new permafrost-associated gas hydrate formation mechanism with a 2-D, multiphase fluid flow model, using the Canadian Beaufort Shelf as an example. The numerical model includes freeze/thaw permafrost processes and predicts the unsteady, 2-D methane solubility field for hydrate inventory calculations. Model results show that widespread, low-saturation hydrate deposits accumulate within and below submarine permafrost, even if offshore-flowing groundwater is undersaturated in methane gas. While intrapermafrost hydrate inventory varies widely depending on permafrost extent, subpermafrost hydrate stability remains largely intact across consecutive glacial cycles, allowing widespread subpermafrost accumulation over time. Methane gas escape to the sediment surface (atmosphere) is predicted along the seaward permafrost boundary during the early to middle years of each glacial epoch; however, if free gas is trapped within the forming permafrost layer instead, venting may be delayed until ocean transgression deepens the permafrost table during interglacial periods, and may be related to the spatial distribution of observed pingo-like features (PLFs) on the Canadian Beaufort Shelf. Shallow, gas-charged sediments are predicted above the gas hydrate stability zone at the midshelf to shelf edge and the upper slope, where a gap in hydrate stability allows free gas to accumulate and numerous PLFs have been observed.

The most climate-sensitive methane hydrate deposits occur on upper continental slopes at depths close to the minimum pressure and maximum temperature for gas hydrate stability. At these water depths, small perturbations in intermediate ocean water temperatures can lead to gas hydrate dissociation. The Arctic Ocean has experienced more dramatic warming than lower latitudes, but observational data have not been used to study the interplay between upper slope gas hydrates and warming ocean waters. Here we use (a) legacy seismic data that constrain upper slope gas hydrate distributions on the U.S. Beaufort Sea margin, (b) Alaskan North Slope borehole data and offshore thermal gradients determined from gas hydrate stability zone thickness to infer regional heat flow, and (c) 1088 direct measurements to characterize multidecadal intermediate ocean warming in the U.S. Beaufort Sea. Combining these data with a three-dimensional thermal model shows that the observed gas hydrate stability zone is too deep by 100 to 250 m. The disparity can be partially attributed to several processes, but the most important is the reequilibration (thinning) of gas hydrates in response to significant (~0.5°C at 2σ certainty) warming of intermediate ocean temperatures over 39 years in a depth range that brackets the upper slope extent of the gas hydrate stability zone. Even in the absence of additional ocean warming, 0.44 to 2.2 Gt of methane could be released from reequilibrating gas hydrates into the sediments underlying an area of ~5–7.5 × 103 km2 on the U.S. Beaufort Sea upper slope during the next century.

Natural gas hydrate in sand sediments can increase both the measured compressional velocity and resistivity. The same geophysical signatures occur, however, in low-porosity sand. We investigate the possible occurrence of natural gas hydrate in a sand interval in Alaminos Canyon Block 21 (AC 21) in the Gulf of Mexico, drilled by the U.S. Gas Hydrate Joint Industry Project. The sand interval has an increase in resistivity (~2.2 Ω m) and a strong peak and trough at the top and bottom of the sand on exploration seismic, which has been interpreted as natural gas hydrate. We reexamine the logging data and construct a new reservoir model that matches the measured resistivity, the high-density sublayers in the sand, and the surface seismic trace. Our model shows that the sand interval in AC 21 is most likely water saturated; and the slight increase in resistivity, higher-measured density, and the seismic amplitudes are caused by a reduction in porosity to ~30% in the sand interval relative to a porosity of ~42% in the surrounding marine muds. More broadly, we show that the average depth where the porosity of marine muds becomes lower than sand sediment is 900 mbsf, though it could be as shallow as 600 mbsf for high-porosity sands. In any case, the similar geophysical signatures for water-saturated sand and low saturations of natural gas hydrate in sand probably occur throughout the gas hydrate stability zone at most sites worldwide.

Gas hydrate and free gas accumulation in heterogeneous marine sediment is simulated using a two-dimensional (2-D) numerical model that accounts for mass transfer over geological timescales. The model extends a previously documented one-dimensional (1-D) model such that lateral variations in permeability (k) become important. Various simulations quantitatively demonstrate how focused fluid flow through high-permeability zones affects local hydrate accumulation and saturation. Simulations that approximate a vertical fracture network isolated in a lower permeability shale (kfracture >> kshale) show that focused fluid flow through the gas hydrate stability zone (GHSZ) produces higher saturations of gas hydrate (25–70%) and free gas (30–60%) within the fracture network compared to surrounding shale. Simulations with a dipping, high-permeability sand layer also result in elevated saturations of gas hydrate (60%) and free gas (40%) within the sand because of focused fluid flow through the GHSZ. Increased fluid flux, a deep methane source, or both together increase the effect of flow focusing upon hydrate and free gas distribution and enhance hydrate and free gas concentrations along the high-permeability zones. Permeability anisotropy, with a vertical to horizontal permeability ratio on the order of 10−2, enhances transport of methane-charged fluid to high-permeability conduits. As a result, gas hydrate concentrations are enhanced within these high-permeability zones. The dip angle of these high-permeability structures affects hydrate distribution because the vertical component of fluid flux dominates focusing effects. Hydrate and free gas saturations can be characterized by a local Peclet number (localized, vertical, focused, and advective flux relative to diffusion) relative to the methane solubility gradient, somewhat analogous to such characterization in 1-D systems. Even in lithologically complex systems, local hydrate and free gas saturations might be characterized by basic parameters (local flux and diffusivity).

Gas seepage from marine sediments has implications for understanding feedbacks between the global carbon reservoir, seabed ecology and climate change. Although the relationship between hydrates, gas chimneys and seafloor seepage is well established, the nature of fluid sources and plumbing mechanisms controlling fluid escape into the hydrate zone and up to the seafloor remain one of the least understood components of fluid migration systems. In this study we present the analysis of new three-dimensional high-resolution seismic data acquired to investigate fluid migration systems sustaining active seafloor seepage at Omakere Ridge, on the Hikurangi subduction margin, New Zealand. The analysis reveals at high resolution, complex overprinting fault structures (i.e. protothrusts, normal faults from flexural extension, and shallow (

[1] Identifying the spatial distribution of seabed fluid expulsion features is crucial for understanding the substrate plumbing system of any continental margin. A 1100 km stretch of the U.S. Atlantic margin contains more than 5000 pockmarks at water depths of 120 m (shelf edge) to 700 m (upper slope), mostly updip of the contemporary gas hydrate stability zone (GHSZ). Advanced attribute analyses of high-resolution multichannel seismic reflection data reveal gas-charged sediment and probable fluid chimneys beneath pockmark fields. A series of enhanced reflectors, inferred to represent hydrate-bearing sediments, occur within the GHSZ. Differential sediment loading at the shelf edge and warming-induced gas hydrate dissociation along the upper slope are the proposed mechanisms that led to transient changes in substrate pore fluid overpressure, vertical fluid/gas migration, and pockmark formation.

[1] For the first time, we report the amplitude variation with angle (AVA) pattern of bottom-simulating reflectors (BSRs) beneath fracture-filled gas hydrate deposits when the effective medium is anisotropic. The common depth point (CDP) gathers of two mutually perpendicular multichannel seismic profiles, located in the vicinity of Site NGHP-01-10, are appropriately processed such that they are fit for AVA analysis. AVA analysis of the BSR shows normal-incidence reflection coefficients of −0.04 to −0.11 with positive gradients of 0.04 to 0.31 indicating class IV pattern. The acoustic properties from isotropic rock physics model predict class III AVA pattern which cannot explain the observed class IV AVA pattern in Krishna-Godavari basin due to the anisotropic nature of fracture-filled gas hydrate deposits. We modeled the observed class IV AVA of the BSR by assuming that the gas hydrate bearing sediment can be represented by horizontally transversely isotropic (HTI) medium after accounting for anisotropic wave propagation effects on BSR amplitudes. The effective medium properties are estimated using Backus averaging technique, and the AVA pattern of BSRs is modeled using the properties of overlying HTI and underlying isotropy/HTI media with or without free gas. Anisotropic AVA analysis of the BSR from the inline seismic profile shows 5–30% gas hydrate concentration (equivalent to fracture density) and the azimuth of fracture system (fracture orientation) with respect to the seismic profile is close to 45°. Free gas below the base of gas hydrate stability zone is interpreted in the vicinity of fault system (F1).

The Paleocene/Eocene thermal maximum (PETM) was a time of rapid global warming in both marine and continental realms that has been attributed to a massive methane (CH 4 ) release from marine gas hydrate reservoirs. Previously proposed mechanisms for this methane release rely on a change in deepwater source region(s) to increase water temperatures rapidly enough to trigger the massive thermal dissociation of gas hydrate reservoirs beneath the seafloor. To establish constraints on thermal dissociation, we model heat flow through the sediment column and show the effect of the temperature change on the gas hydrate stability zone through time. In addition, we provide seismic evidence tied to borehole data for methane release along portions of the U.S. continental slope; the release sites are proximal to a buried Mesozoic reef front. Our model results, release site locations, published isotopic records, and ocean circulation models neither confirm nor refute thermal dissociation as the trigger for the PETM methane release. In the absence of definitive evidence to confirm thermal dissociation, we investigate an alternative hypothesis in which continental slope failure resulted in a catastrophic methane release. Seismic and isotopic evidence indicates that Antarctic source deepwater circulation and seafloor erosion caused slope retreat along the western margins of the North Atlantic in the late Paleocene. Continued erosion or seismic activity along the oversteepened continental margin may have allowed methane to escape from gas reservoirs trapped between the frozen hydrate-bearing sediments and the underlying buried Mesozoic reef front, precipitating the Paleocene/Eocene boundary methane release. An important implication of this scenario is that the methane release caused (rather than resulted from) the transient temperature increase of the PETM. Neither thermal dissociation nor mechanical disruption of sediments can be identified unequivocally as the triggering mechanism for methane release with existing data. Further documentation with high-resolution benthic foraminiferal isotopic records and with seismic profiles tied to borehole data is needed to clarify whether erosion, thermal dissociation, or a combination of these two was the triggering mechanism for the PETM methane release.

During Ocean Drilling Program Leg 164, shear sonic velocity and other geophysical logs were acquired in gas hydrate-bearing sediments on the Blake Ridge to characterize the very distinct seismic signature of such formations: anomalous low amplitudes overlying a strong bottom-simulating reflector (BSR). A comparison of the bulk moduli derived from the logs to standard elastic consolidation models shows that the sediments are overconsolidated above the BSR at 440 meters below seafloor (mbsf) because of the presence of hydrates. Below the bottom of the thermodynamic hydrate stability zone at ∼520 mbsf, the high compressibility of the formation and the attenuation of the monopole sonic waveforms are typical of sediments partially saturated with free gas. Between these two depths, gas hydrate and free gas seem to coexist. Within the Gas hydrate stability zone, we estimate the amount of gas hydrates using different models based on theories for wave scattering in multiphase media and for grain cementation. In close agreement with measurements made on discrete in situ samples, the latter describes most accurately the interactions between the matrix, the pore fluids, and the hydrates. This model indicates that 5 to 10% of the pore space is occupied by hydrates deposited uniformly on the surface of the grains. The comparison with Gassmann's model also show that the amount of free gas below the BSR never exceeds 5% of the pore space but is high enough to generate the BSR. The coexistence of free gas and gas hydrates below the BSR may be explained by capillary effects in the smaller pores or by remaining crystalline structures after partial hydrate decomposition.

We offer a first-principle-based effective medium model for elastic-wave velocity in unconsolidated, high porosity, ocean bottom sediments containing gas hydrate. The dry sediment frame elastic constants depend on porosity, elastic moduli of the solid phase, and effective pressure. Elastic moduli of saturated sediment are calculated from those of the dry frame using Gassmann's equation. To model the effect of gas hydrate on sediment elastic moduli we use two separate assumptions: (a) hydrate modifies the pore fluid elastic properties without affecting the frame; (b) hydrate becomes a component of the solid phase, modifying the elasticity of the frame. The goal of the modeling is to predict the amount of hydrate in sediments from sonic or seismic velocity data. We apply the model to sonic and VSP data from ODP Hole 995 and obtain hydrate concentration estimates from assumption (b) consistent with estimates obtained from resistivity, chlorinity and evolved gas data.

[1] The potential impact of future climate change on methane release from oceanic gas hydrates is the subject of much debate. We analyzed World Ocean Database quality controlled data on the Norwegian-Svalbard continental margin from the past 60 years to evaluate the potential effect of ocean temperature variations on continental margin gas hydrate reservoirs. Bottom water temperatures in the Norwegian-Svalbard margin were subject to significant cooling until 1980 (by ∼2°C offshore NW-Svalbard and in the Barents Sea) followed by a general bottom water temperature increase until 2010 (∼0.3°C in deep-water areas offshore NW-Svalbard and mid-Norwegian margin and ∼2°C in the shallow areas of the Barents Sea and Prins Karls Forland). Bottom water warming in the shallow outer shelf areas triggered the Gas Hydrate Stability Zone (GHSZ) retreat toward upper continental slope areas, potentially increasing methane release due to gas hydrate dissociation. GHSZ responses to temperature changes on human time scales occur exclusively in shallow water and only if near-surface gas hydrates exist. The responses are associated with a short time lag of less than 1 year. Temperatures in the bottom water column seem to be partly regulated by the North Atlantic Oscillation (NAO), with positive NAO associated with warm phases. However, cooling events in the surface water offshore NW-Svalbard might be associated with El Nino events of 1976–1977, 1986–1987 and 1997–1998 in the Pacific. Such ocean cooling, if long enough, may delay ocean temperature driven gas hydrate dissociation and potential releases of methane to the ocean.

[1] Evidence for twelve sites with gas bubble emissions causing hydroacoustic anomalies in 18 kHz echosounder records (‘flares’) was obtained at the convergent Makran continental margin. The hydroacoustic anomalies originating from hydrocarbon seeps at water depths between 575 and 2870 m disappeared after rising up to 2000 m in the water column. Dives with the remotely operated vehicle ‘Quest 4000 m’ revealed that several individual bubble vents contributed to one hydroacoustic anomaly. Analyzed gas samples suggest that bubbles were mainly composed of methane of microbial origin. Bubble size distributions and rise velocities were determined and the volume flux was estimated by counting the emitted bubbles and using their average volume. We found that a low volume flux (Flare 1 at 575 mbsl: 90 ml/min) caused a weak hydroacoustic signal in echograms whereas high volume fluxes (Flare 2 at 1027 mbsl: 1590 ml/min; Flare 5 C at 2870 mbsl: 760 ml/min) caused strong anomalies. The total methane bubble flux in the study area was estimated by multiplying the average methane flux causing a strong hydroacoustic anomaly in the echosounder record with the total number of equivalent anomalies. An order-of-magnitude estimate further considers the temporal variability of some of the flares, assuming a constant flux over time, and allows a large range of uncertainty inherent to the method. Our results on the fate of bubbles and the order-of-magnitude estimate suggest that all of the ∼40 ± 32 × 106 mol methane emitted per year within the gas hydrate stability zone remain in the deep ocean.

In fine-grained, faulted sediments, both stratigraphic and fault-induced structural variations can simultaneously determine the gas hydrate distribution. Insights into hydrate distribution can be obtained from P-wave velocity (VP) and attenuation (QP -1 ) character of the gas hydrate stability zone (GHSZ). In this paper, we apply frequency-domain full-waveform inversion (FWI) to surface-towed 2D multichannel seismic data from the Krishna-Godavari (KG) Basin, India, to image the fine-scale (100 X 30 m) VP and QP -1 variations within the GHSZ. We validate the inverted VP model by reconciling it with a sonic log from a nearby (~250 m) well. The VP model shows a patchy distribution of hydrate. Away from the faults-dominated parts of the profile, hydrates demonstrate stratigraphic control which appears to be permeability driven. The QP -1 model suggests that attenuation is relatively suppressed in hydratesbearing sediments. Elevated attenuation in non-hydrate-bearing sediments could be driven by the apparent pore-fluid immiscibility at seismic wavelengths. The VP and the QP -1 models also suggest that fault zones within the GHSZ can be hydrate- or free-gas-rich depending on the relative supply of free gas and water from below the GHSZ.

[1] Global assessment of methane must consider the role of mid-ocean ridges. Fluids from serpentinized mantle and gabbro material are noteworthy on ocean ridges, although they are not very well understood. Only a few sedimented ocean ridges exist worldwide, and they may document past and ongoing serpentinization-driven migration of gas-rich fluids. This study is based on two multichannel reflection seismic profiles acquired across a sedimented segment of the ultraslow spreading Knipovich Ridge offshore NW Svalbard. Seismic data allow suggesting a potential link between inferred areas of serpentinization, transfer of carbon from the deep-seated host rocks through the sediments above by diapirism, and methane capture within the gas hydrate stability zone at the eastern flank of the Knipovich Ridge. The origin of sediment remobilization features can be related to intrusions and the degassing process from mantle serpentinization. These disturbances in sediments overlying the oceanic crust can be observed in seismic data and are interpreted as diapirs. In shallower sediments, at the predicted base of the gas hydrate stability zone, the seismic data show a bright spot with all the characteristics of a gas hydrate related bottom-simulating reflector (BSR), such as enhanced reflection amplitude, phase reversal relative to the seabed reflection, and crosscutting of sedimentary strata. The BSR occurs at about 200 ms two-way time within a sequence of marine sediments. Two-dimensional concentration models of methane hydrate using the differential effective medium theory predict saturations of up to 26% of methane hydrate in the pore space of sediments in the gas hydrate reservoir.

The ongoing warming of bottom water in the Arctic region is anticipated to destabilize some of the gas hydrate present in shallow seafloor sediment, potentially causing the release of methane from dissociating hydrate into the ocean and the atmosphere. Ocean-bottom seismometer (OBS) experiments were conducted along the continental margin of western Svalbard to quantify the amount of methane present as hydrate or gas beneath the seabed. P- and S-wave velocities were modeled for five sites along the continental margin, using ray-trace forward modeling. Two southern sites were located in the vicinity of a 30 km long zone where methane gas bubbles escaping from the seafloor were observed during the cruise. The three remaining sites were located along an E-W orientated line in the north of the margin. At the deepest northern site, Vp anomalies indicate the presence of hydrate in the sediment immediately overlying a zone containing free gas up to 100-m thick. The acoustic impedance contrast between the two zones forms a bottom-simulating reflector (BSR) at approximately 195 m below the seabed. The two other sites within the gas hydrate stability zone (GHSZ) do not show the clear presence of a BSR or of gas hydrate. However, anomalously low Vp, indicating the presence of free gas, was modeled for both sites. The hydrate content was estimated from Vp and Vs, using effective-medium theory. At the deepest northern site, modeling suggests a pore-space hydrate concentration of 7–12%, if hydrate forms as part of a connected framework, and about 22% if it is pore-filling. At the two other northern sites, located between the deepest site and the landward limit of the GHSZ, we suggest that hydrate is present in the sediment as inclusions. Hydrate may be present in small quantities at these two sites (4–5%) of the pore space. The variation in lithology for the three sites indicated by high-resolution seismic profiles may control the distribution, concentration and formation of hydrate and free gas.

[1] The strong coupling between multiphase flow and sediment mechanics determines the spatial distribution and migration dynamics of gas percolating through liquid-filled soft granular media. Here, we investigate, by means of controlled experiments and computed tomography (CT) imaging, the preferential mode of gas migration in three-dimensional samples of water-saturated silica-sand and silica-silt sediments. Our experimental system allowed us to independently control radial and axial confining stresses and pore pressure while performing continuous x-ray CT scanning. The CT image analysis of the three-dimensional gas migration provides the first experimental confirmation that capillary invasion preferentially occurs in coarse-grained sediments whereas grain displacement and conduit openings are dominant in fine-grained sediments. Our findings allow us to rationalize prior field observations and pore-scale modeling results, and provide critical experimental evidence to explain the means by which conduits for the transit of methane gas may be established through the gas hydrate stability zone in oceanic sediments, and cause large episodic releases of carbon into the deep ocean.

The occurrence of gas hydrates at submarine mud volcanoes (MVs) located within the gas hydrate stability zone (GHSZ) is controlled by upward fluid and heat flux associated with MV activity. Determining the spatial distribution of gas hydrates at MVs is crucial to evaluate their sensitivity to known episodic changes in volcanic activity. We determined the hydrocarbon inventory and spatial distribution of hydrates at an individual MV structure. The Hakon Mosby Mud Volcano (HMMV), located at 1,250 m water depth on the Barents Sea slope, was investigated by combined pressure core sampling, heat flow measurements, and pore water chemical analysis. Quantitative pressure core degassing revealed gas-sediment ratios between 3.1 and 25.7, corresponding to hydrate concentrations of up to 21.3% of the pore volume. Hydrocarbon compositions and physicochemical conditions imply that gas hydrates incipiently crystallize as structure I hydrate, with a dissociation temperature of around 13.8 degrees C at this water depth. Based on numerous in situ measurements of the geothermal gradient in the seabed, pore water sulfate profiles and microbathymetric data, we show that the thickness of the GHSZ increases from less than 1 m at the warm center to around 47 m in the outer parts of the HMMV. We estimate the total mass of hydrate-bound methane stored at the HMMV to be about 102.5 kt, of which 2.8 kt are located within the morphological Unit I around the center and thus are likely to be dissociated in the course of a large eruption.

Gas hydrate saturations were estimated using five different methods in silt and silty clay foraminiferous sediments from drill hole SH2 in the South China Sea. Gas hydrate saturations derived from observed pore water chloride values in core samples range from 10 to 45% of the pore space at 190-221 m below seafloor (mbsf). Gas hydrate saturations estimated from resistivity (R-t) using wireline logging results are similar and range from 10 to 40.5% in the pore space. Gas hydrate saturations were also estimated by P wave velocity obtained during wireline logging by using a simplified three-phase equation (STPE) and effective medium theory (EMT) models. Gas hydrate saturations obtained from the STPE velocity model (41.0% maximum) are slightly higher than those calculated with the EMT velocity model (38.5% maximum). Methane analysis from a 69 cm long depressurized core from the hydrate-bearing sediment zone indicates that gas hydrate saturation is about 27.08% of the pore space at 197.5 mbsf. Results from the five methods show similar values and nearly identical trends in gas hydrate saturations above the base of the gas hydrate stability zone at depths of 190 to 221 mbsf. Gas hydrate occurs within units of clayey slit and silt containing abundant calcareous nannofossils and foraminifer, which increase the porosities of the fine-grained sediments and provide space for enhanced gas hydrate formation. In addition, gas chimneys, faults, and fractures identified from three-dimensional (3-D) and high-resolution two-dimensional (2-D) seismic data provide pathways for fluids migrating into the gas hydrate stability zone which transport methane for the formation of gas hydrate. Sedimentation and local canyon migration may contribute to higher gas hydrate saturations near the base of the stability zone.

Gas hydrates are naturally ocurring solids consisting of water molecules forming a lattice of cages, most of which contain a molecule of natural gas, usually methane. The present article discusses three important aspects of gas hydrates: their potential as a fossil fuel resource, their role as a submarine geohazard, and their effects on global climate change. 70 refs., 16 figs., 1 tab.

[1] The distribution of gas hydrate on a continental slope is often characterized as a wedge that pinches out on the seafloor. This part of the hydrate stability zone is particularly relevant for studies of the dynamics of hydrate accumulations, such as processes related to slope stability or hydrate dissociation leading to methane release into the overlying ocean. For regions with very low geothermal gradients, we have produced a series of thermobaric models of the shallow hydrate stability zone that contain an unexpected geometrical distribution of hydrated sediments. In these models, the shallowest part of the stability field thickens and bulges landward. Such a feature is more likely to happen in regions where low geothermal gradients are further lowered by high sedimentation rates. Also, the effect is greater beneath colder oceans. Although a hydrate stability zone bulge would be difficult to image with conventional seismic methods, there are numerous locations around the world where such a system could develop.

[1] Vast quantities of methane are trapped in oceanic hydrate deposits, and there is concern that a rise in the ocean temperature will induce dissociation of these hydrate accumulations, potentially releasing large amounts of methane into the atmosphere. The recent discovery of active methane gas venting along the landward limit of the gas hydrate stability zone (GHSZ) on the shallow continental slope west of Spitsbergen could be an indication of this process, if the source of the methane can be confidently attributed to dissociating hydrates. In the first large-scale simulation study of its kind, we simulate shallow hydrate dissociation in conditions representative of the West Spitsbergen margin to test the hypothesis that the observed gas release originated from hydrates. The simulation results are consistent with this hypothesis, and are in remarkable agreement with the recently published observations. They show that shallow, low-saturation hydrate deposits, when subjected to temperature increases at the seafloor, can release significant quantities of methane, and that the releases will be localized near the landward limit of the top of the GHSZ. These results indicate the possibility that hydrate dissociation and methane release may be both a consequence and a cause of climate change.

[1] A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate-saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate-bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

More than 250 plumes of gas bubbles have been discovered emanating from the seabed of the West Spitsbergen continental margin, in a depth range of 150-400 m, at and above the present upper limit of the gas hydrate stability zone (GHSZ). Some of the plumes extend upward to within 50 m of the sea surface. The gas is predominantly methane. Warming of the northward-flowing West Spitsbergen current by 1°C over the last thirty years is likely to have increased the release of methane from the seabed by reducing the extent of the GHSZ, causing the liberation of methane from decomposing hydrate. If this process becomes widespread along Arctic continental margins, tens of Teragrams of methane per year could be released into the ocean.

[1] High-angle gas hydrate filled fracture planes were identified along a 31 m interval in logging while drilling images in two holes located ∼11 m apart drilled during the Indian National Gas Hydrate Program Expedition 01, offshore India. Using Monte Carlo simulations to account for uncertainty in hole location, hole deviation, strike and dip, we assert with 95% confidence that the fracture planes in the two holes are not the same. The gas hydrate filled fracture planes likely only extend a few meters laterally from each borehole and occur in an isolated interval in the middle of the gas hydrate stability zone. This suggests gas generated microbially within in the gas hydrate stability zone may have supplied the gas hydrate-filled fracture interval. Production of methane from these reservoirs using conventional methods may be quite challenging.

New 3-D seismic investigations carried out across the Sevastopol mud volcano in the Sorokin Trough present 3-D seismic data of a mud volcano in the Black Sea for the first time. The studies allow us to image the complex three-dimensional morphology of a collapse structured mud volcano and to propose an evolution model. The Sevastopol mud volcano is located above a buried diapiric structure with two ridges and controlled by fluid migration along a deep fault system, which developed during the growth of the diapirs in a compressional tectonic system. Overpressured fluids initiated an explosive eruption generating the collapse depression of the Sevastopol mud volcano. Several cones were formed within the depression by subsequent quiet mud extrusions. Although gas hydrates have been recovered at various mud volcanoes in the Sorokin Trough, no gas hydrates were sampled at the Sevastopol mud volcano. A BSR (bottom-simulating reflector) is missing in the seismic data; however, high-amplitude reflections (bright spots) observed above the diapiric ridge near the mud volcano at a relatively constant depth correspond to the approximate depth of the base of the gas hydrate stability zone (BGHSZ). Thus we suggest that gas hydrates are present locally where gas/fluid flow occurs related to mud volcanism, i.e., above the diapir and close to the feeder channel of the mud volcano. Depth variations of the bright spots of up to 200 ms TWT might be caused by temperature variations produced by variable fluid flow.

[1] To elucidate the geological processes associated with hydrate formation and dissociation in the marine environment under a wide range of conditions, we have developed a one-dimensional numerical computer model (simulator). The numerical model can be used to simulate the following aspects of hydrate formation, decomposition, reformation, and distribution: (1) burial history of deep marine sediments and associated phenomena (e.g., sediment compaction and consequent reduction in sediment porosity and permeability, fluid expulsion, time evolution of temperature and pressure, heat flux), (2) in situ generation of biogenic methane from buried organic carbon and methane solubility in formation brine, (3) methane hydrate formation, decomposition, reformation, and distribution in response to changes in gas concentration, pressure, temperature and fluid salinity (the hydrate formation and decomposition are treated as equilibrium processes), (4) influence of sulfate reduction zone under the seafloor on hydrate formation, (5) possible presence of a free gas zone beneath the gas hydrate stability zone, and (6) multiphase (i.e., liquid brine with dissolved gas, free gas, gas hydrate) flow through a deformable porous matrix. The model provides for a reduction/ increase in permeability due to the formation/decomposition of the gas hydrate. Initial applications of the model to study hydrate distribution at the Blake Ridge (site 997) and Hydrate Ridge (site 1249) are described. Model results are compared with chlorinity, sulfate, and hydrate distribution data.

While it is well accepted that gas hydrate dissociation at the base of the Gas Hydrate Stability Zone (GHSZ) can generate high excess pore pressure and leads to sediment deformation, the consequence in terms of pore pressure of the dissolution of the gas hydrate at the top of the Gas Hydrate Occurrence Zone (GHOZ) remains neglected. The purpose of this comment on Xu and Germanovich [2006] article is to demonstrate that gas hydrate dissolution in the GHSZ may generate excess pore pressure and to point out the risk related to hydrate dissolution at the top of the GHOZ.

[1] Methane hydrate bottom-simulating reflectors (BSRs) are widespread on the northern Cascadia margin offshore Vancouver Island. We conducted a three-dimensional tomographic seismic study of the hydrate stability zone in an area around Ocean Drilling Program Site 889 using two deployments of five ocean bottom hydrophones and air gun shots along a series of closely spaced profiles in various orientations. Further constraints on reflector geometry come from coincident single-channel reflection profiles. Travel times of reflected and refracted phases were inverted with a regularized three-dimensional inversion using perturbation ray tracing through smooth isotropic media for the forward step. The seismic data allow us to constrain the velocity structure in a ∼6 km2 area around the drill site. Mean velocities range from 1.50 km s−1 at the seabed to 1.84 km s−1 at the BSR, and velocities at Site 889 match well those measured using a vertical seismic profile. At equivalent depths below the seafloor, velocities vary laterally by typically ∼0.15 km s−1. Close to the seafloor, velocities may be controlled primarily by lithology, but close to the BSR we infer hydrate contents of up to 15% of the pore space from effective medium modeling. The mean hydrate saturation in the well-constrained volume of the velocity model is estimated to be 2.2%. There is no correlation between the seismic velocity above the BSR and the reflection coefficient at the BSR, so the latter is likely controlled primarily by the distribution of free gas beneath the hydrate stability zone.

To constrain the fluxes of methane (CH4) in the water column above the accretionary wedge along the Cascadia continental margin, we measured methane and its stable carbon isotope signature (δ13C-CH4). The studies focused on Hydrate Ridge (HR), where venting occurs in the presence of gas-hydrate-bearing sediments. The vent CH4 has a light δ13C-CH4 biogenic signature (−63 to −66‰ PDB) and forms thin zones of elevated methane concentrations several tens of meters above the ocean floor in the overlying water column. These concentrations, ranging up to 4400 nmol L−1, vary by 3 orders of magnitude over periods of only a few hours. The poleward undercurrent of the California Current system rapidly dilutes the vent methane and distributes it widely within the gas hydrate stability zone (GHSZ). Above 480 m water depth, the methane budget is dominated by isotopically heavier CH4 from the shelf and upper slope, where mixtures of various local biogenic and thermogenic methane sources were detected (−56 to −28‰ PDB). The distribution of dissolved methane in the working area can be represented by mixtures of methane from the two primary source regions with an isotopically heavy background component (−25 to −6‰ PDB). Methane oxidation rates of 0.09 to 4.1% per day are small in comparison to the timescales of advection. This highly variable physical regime precludes a simple characterization and tracing of “downcurrent” plumes. However, methane inventories and current measurements suggest a methane flux of approximately 3 × 104 mol h−1 for the working area (1230 km2), and this is dominated by the shallower sources. We estimate that the combined vent sites on HR produce 0.6 × 104 mol h−1, and this is primarily released in the gas phase rather than dissolved within fluid seeps. There is no evidence that significant amounts of this methane are released to the atmosphere locally.

[1] Wide–angle reflection seismic experiments were performed at the Storegga slide offshore Norway in 2002 with the goal to quantify the amount of gas hydrate and free gas in the sediment. Twenty-two stations with Ocean Bottom Hydrophones (OBH) and Seismometers (OBS) were deployed for a 2D and a 3D experiment. Kirchhoff depth migration is used to transform the seismic wide–angle data into images of the sediment layers and to obtain P wave velocity–depth functions. The gas hydrate and free gas saturations are estimated from the elastic properties of the sediment on the basis of the Frenkel–Gassmann equations. There is 5–15% gas hydrate in the pore space of the sediment in the gas hydrate stability zone (GHSZ). The free gas saturation takes the value of 0.8% for a homogeneous distribution of gas in the pore water and 7% for the model of a patchy gas distribution.

[1] High-resolution three-dimensional (3-D) seismic reflection data acquired on the R/V Thomas G. Thompson in 2000 reveal a pair of bottom simulating reflections (BSRs) across a broad region of southern Hydrate Ridge, offshore Oregon. The primary BSR (BSRp) is a regionally extensive reflection that lies 120–150 m below seafloor and exhibits typical characteristics of a gas hydrate BSR. We also imaged a second weaker BSR (BSRs), 20–40 m below BSRp, with similar characteristics. BSRs is interpreted as a remnant of a BSR that probably formed during the Last Glacial Maximum 18,000 years ago, when the base of the gas hydrate stability zone (GHSZ) was deeper. An increase in bottom water temperatures of 1.75°–2.25° and a corresponding sea level rise of 120 m could have produced the BSR shift. The preservation of BSRs for at least 5000 years, which is the time since subseafloor temperatures stabilized following ocean warming after the Last Glacial Maximum, implies very slow upward advective and diffusive flow of methane (

[1] Gas hydrates inside mud volcanoes have been observed in several locations but are generally found at water depths of 1000 m and deeper. We present the first observation of the base of a gas hydrate stability zone within a shallow mud volcano in the El Arraiche mud volcano field on the Moroccan Atlantic margin. The mud volcano base is located at about 475 m and is over 125 m high. On high-resolution seismics we observed an anomalous but coherent reflection under the slopes of the mud volcano. The event was interpreted as the base of a gas hydrate stability zone because of its inverse polarity and its morphology. Far from the crater, the event is nearly parallel to the seafloor. Closer toward the crater, the event shallows. Inside the mud volcano crater, no event is observed. A stability model using thermogenic gas compositions is applied to local P-T conditions, indicating that thermogenic gas hydrates can be stable at this depth. The high modeled heat flow in the crater of the mud volcano indicates a focused flow of warm fluids. Below the slopes of the mud volcano, the inferred heat flow is also elevated but less high. In areas of thermogenic gas production, gas hydrates can occur at shallow water depths, even in areas with high heat flow. This also suggests that dewatering of the accretionary wedge complex is mainly focused along fault surfaces and through seafloor structures, such as mud volcanoes.

Log and core data document gas saturations as high as 90% in a coarse-grained turbidite sequence beneath the gas hydrate stability zone (GHSZ) at south Hydrate Ridge, in the Cascadia accretionary complex. The geometry of this gas-saturated bed is defined by a strong, negative-polarity reflection in 3D seismic data. Because of the gas buoyancy, gas pressure equals or exceeds the overburden stress immediately beneath the GHSZ at the summit. We conclude that gas is focused into the coarse-grained sequence from a large volume of the accretionary complex and is trapped until high gas pressure forces the gas to migrate through the GHSZ to seafloor vents. This focused flow provides methane to the GHSZ in excess of its proportion in gas hydrate, thus providing a mechanism to explain the observed coexistence of massive gas hydrate, saline pore water and free gas near the summit.

Previous studies have revealed the presence of a widespread Bottom Simulating Reflector (BSR) on the eastern Nankai slope, as well as the occurrence of enigmatic high amplitude reflections that extend well above the BSR. New high-resolution 2-D seismic data were collected on the eastern Nankai slope, during the French-Japanese SFJ cruise in year 2000 and AVA analyses of the enigmatic reflectors are conducted. At the studied location, these analyses suggest that high amplitude anomalies above the BSR delineate the top of gas hydrate rich sediments. Several tens % of the sediment porosity would be filled with gas hydrates between the BSR and a sharp boundary 30 to 60 m above. To account for these observations, we propose that an invasion by free gas of the present day gas hydrate stability zone occurred in the past. Several mechanisms for this intrusion are discussed.

[1] Gas hydrate (GH) stability modeling results explain why some major Holocene submarine landslides along the Norwegian-Barents margin could have been triggered by GH dissociation during the early to middle Holocene, not during the lowest sea levels of the Last Glacial Maximum (LGM). Our model results show that subbottom depths of 170-260 m below the pre-slide continental slope (ca. 350-475 m present water depth) must have passed out of gas hydrate stability zone (GHSZ) by 8.15 ka as the effect of warm bottom water inflow at 11 ka penetrated into the subbottom, overcoming the effects of pressure increase due to sea level rise (SLR). The component of local SLR due to the isostatic response to Fennoscandian deglaciation is shown to be relatively insignificant, particularly for the part of the upper continental slope where the slide probably began. The stability relations show that GH could have formed under the ice sheet before deglaciation, and below deeper shelf areas after sea levels began to rise, but before significant warming near the GHSZ base. To the extent water deeper than 800 m has remained cold (-1° to 0°C) since LGM times, the GHSZ continued to thicken in deep water and GH dissociation could not have triggered Holocene failure in that regime. The present distribution of GH stability is limited to water depths greater than about 400 m in the Storegga slide area, and the thickness of the GHSZ increases with water depth.

[1] Review of the literature reveals that the nature of pore-scale interactions between gas hydrates and porous media remains a matter of controversy. To clarify the situation, nuclear magnetic resonance (NMR) measurements have been made on methane hydrate-bearing sandstones. The samples were synthetically prepared within the gas hydrate stability zone, at or near the seafloor in Monterey Bay, California. The method simulated natural hydrate deposition by gas flows that are not in thermodynamic equilibrium with the surrounding earth. The efficiency of hydrate production was variable, as has been observed elsewhere. When substantial hydrate saturations were achieved, NMR relaxation time measurements indicated that hydrate tended to replace water in the largest pore spaces. The relative permeability to water, as determined by an NMR-based correlation, was significantly reduced. The magnitude of this reduction was also consistent with formation of hydrate in the centers of pores, rather than with hydrate coating the grains. The growth habit suggested by these results is consistent with creation of hydrate nodules and lenses in coarse, unconsolidated sediments. It is also consistent with scenarios in which methane gas is delivered efficiently to the atmosphere as a result of seafloor slope failure, thereby strengthening global warming feedback mechanisms.

Bottom-simulating reflectors (BSR) are observed commonly at a depth of several hundred meters below the seafloor in continental margin sedimentary sections that have undergone recent tectonic consolidation or rapid accumulation. They are believed to correspond to the deepest level at which methane hydrate (clathrate) is stable. We present a model in which BSR hydrate layers are formed through the removal of methane from upward moving pore fluids as they pass into the hydrate stability field. In this model, most of the methane is generated below the level of hydrate stability, but not at depths sufficient for significant thermogenic production; the methane is primarily biogenic in origin. The model requires either a mechanism to remove dissolved methane from the pore fluids or disseminated free gas carried upward with the pore fluid. The model accounts for the evidence that the hydrate is concentrated in a layer at the base of the stability field, for the source of the large amount of methane contained in the hydrate, and for BSRs being common only in special environments. Strong upward fluid expulsion into the hydrate stability field does not occur in normal sediment depositional regimes, so BSRs are uncommon. Upward fluid expulsion does occur as a result of tectonic thickening and loading in subduction zone accretionary wedges and in areas where rapid deposition results in initial undercconsolidation. In these areas hydrate BSRs are common. The most poorly quantified aspect of the model is the efficiency with which methane is removed and hydrate is formed as pore fluids pass into the hydrate stability field. The critical boundary in the phase diagram between the fluid-plus-hydrate and fluid-only fields is not well constrained. However, the amount of methane required to form the hydrate and limited data on methane concentrations in pore fluids from deep-sea boreholes suggest very efficient removal of methane from rising fluid that may contain less than the amount required for free gas production. In most fluid expulsion regimes, the quantity of fluid moved upward to the seafloor is great enough to continually remove the excess chloride and the residue of isotope fractionation resulting from hydrate formation. Thus, as observed in borehole data, there are no large chloride or isotope anomalies remaining in the local pore fluids. The differences in the concentration of methane and probably of CO2 in the pore fluid above and below the base of the stability field may have a significant influence on early sediment diagenetic reactions.

Multichannel seismic reflection data have been analyzed from an area of clear bottom simulating reflectors (BSRs) on the northern Cascadia subduction zone margin off Vancouver Island. The reflector at a depth of about 300 m subbottom is interpreted to represent the base of a layer of methane hydrate or clathrate. The shallow water depth of 1300 m and the 3600-m-long hydrophone array have allowed BSR amplitude-versus-offset and high-resolution velocity analysis, as well as modelling of vertical incidence data. The results of all three types of analysis can be best explained by a 10 to 30-m-thick high-velocity layer located immediately above the BSR about 300 m below the seafloor, having a sharp base and transitional top. In the layer, about one third of the sediment pore spaces must be filled with hydrate “ice”. There is no seismically detectable free gas beneath the BSRs. These results put important constraints on models for the distribution and formation of BSR hydrate.