Description: Highlights • Elongated fault structures are conduits for focused fluid flow. • Gas migration occurs only along a sub-set of faults across Opouawe bank. • Stress state deduced from 3D fault structures appears partially stratigraphically controlled. Abstract High-resolution 2D and 3D seismic data from Opouawe Bank, an accretionary ridge on the Hikurangi subduction margin off New Zealand, show evidence for exceptional gas migration pathways linked to the stress regime of the ridge. Although the ridge has formed by thrusting and folding in response to a sub-horizontal principal compressive stress (σ1), it is clear that local stress conditions related to uplift and extension around the apex of folding (i.e. sub-vertical σ1) are controlling shallow fluid flow. The most conspicuous structural features are parallel and horizontally-elongated extensional fractures that are perpendicular to the ridge axis. At shallower depth near the seafloor, extensional fractures evolve into more concentric structures which ultimately reach the seafloor where they terminate at gas seeps. In addition to the ridge-perpendicular extensional fractures, we also observe both ridge-perpendicular and ridge-parallel normal faults. This indicates that both longitudinal- and ridge-perpendicular extension have occurred in the past. The deepest stratigraphic unit that we image has undergone significant folding and is affected by both sets of normal faults. Shallower stratigraphic units are less deformed and only host the ridge-parallel normal faults, indicating that longitudinal extension was limited to an older phase of ridge evolution. Present-day gas migration has exploited the fabric from longitudinal extension at depth. As the gas ascends to shallower units it ‘self-generates’ its flow pathways through the more concentric structures near the seafloor. This shows that gas migration can evolve from being dependent on inherited tectonic structures at depth, to becoming self-propagating closer to the seafloor.

Description: High-resolution marine seismic data acquisition and subsequent analyses are highly influenced by sea conditions, directly affecting data quality and interpretation. Traditional swell effect correction methods are effective in improving reflector continuity; however, they are less useful for enhancing travel time consistency at intersection points of crossing lines. To develop a robust swell-removal technique for a set of crossing lines multi-beam echo sounder (MBES) data and Chirp sub-bottom profiler (SBP) data were acquired. After generation of a time structure map of the sea-bottom converted from the final processed multi-beam data, a moving average was used to improve the event continuity of the sea-bottom reflection of the Chirp SBP data. Using the position of the Chirp SBP data, the difference between the travel time of the sea-bottom from the smoothed map and the original travel time of the sea-bottom is calculated as a static correction. The static correction method based on the MBES data was compared and verified using three different cases: (i) simple 2D swell effect correction on a line-by-line basis, (ii) comparing the swell corrections at the crossing positions of 2D lines acquired from different dates, and (iii) comparison of ties of intersection points between 2D lines after new swell correction applied. Although a simple 2D swell correction showed great enhancement of reflector continuity, only the full static correction using the newly proposed method using MBES data produced completely corrected reflection events especially at the crossing points of 2D lines.

Description: During the India National Gas Hydrate Program (NGHP) Expedition 01 in 2006 significant sand and gas hydrate were recovered at Site NGHP-01-15 within the Krishna–Godavari Basin, East Coast off India. At the drill site NGHP-01-15, a 5–8 m thick interval was found that is characterized by higher sand content than anywhere else at the site and within the KG Basin. Gas hydrate concentrations were determined to be 20–40% of the pore volume using wire-line electrical resistivity data as well as core-derived pore-fluid freshening trends. The gas hydrate-bearing interval was linked to a prominent seismic reflection observed in the 3D seismic data. This reflection event, mapped for about 1 km2 south of the drill site, is bound by a fault at its northern limit that may act as migration conduit for free gas to enter the gas hydrate stability zone (GHSZ) and subsequently charge the sand-rich layer. On 3D and additional regional 2D seismic data a prominent channel system was imaged mainly by using the seismic instantaneous amplitude attribute. The channel can be clearly identified by changes in the seismic character of the channel fill (sand-rich) and pronounced levees (less sand content than in the fill, but higher than in surrounding mud-dominated sediments). The entire channel sequence (channel fill and levees) has been subsequently covered and back-filled with a more mud-prone sediment sequence. Where the levees intersect the base of the GHSZ, their reflection strengths are significantly increased to 5- to 6-times the surrounding reflection amplitudes. Using the 3D seismic data these high-amplitude reflection edges where linked to the gas hydrate-bearing layer at Site NGHP-01-15. Further south along the channel the same reflection elements representing the levees do not show similarly large reflection amplitudes. However, the channel system is still characterized by several high-amplitude reflection events (a few hundred meters wide and up to ~ 1 km in extent) interpreted as gas hydrate-bearing sand intervals along the length of the channel.

Description: A 1.5 km long, 1 km wide and 70–80 m high carbonate mound was identified on the mid-slope region of the subduction accretionary sedimentary prism offshore Vancouver Island ∼ 3.5 km west of Ocean Drilling Program (ODP) Site 889 and Integrated Ocean Drilling Program (IODP) Site U1327. Seabed-video images show the presence of seafloor carbonate as well as chemosynthetic communities. A high-resolution single channel seismic survey with close line spacing, recording coherent reflectivity down to about 400 m beneath the seafloor, provided acoustic images of this mound and of the gas hydrate bottom-simulating reflector (BSR) beneath it. The mound is interpreted to have developed as a structural topographic high in the hanging wall of a large reverse fault formed at the base of the current seaward slope. The fault zone provides pathways for fluids including gas to migrate to the seafloor where diagenetic carbonate forms and cements the near-surface sediments. To examine the thermal effect of possible upward fluid flow beneath the mound, heat flow at the mound and in the neighbouring region was calculated from the depth of the BSR below the seafloor. These data were combined with heat flow calculated over a broader region from previous multi-channel seismic data. Heat flow within the flattest portion of the surrounding 4 km by 8 km region averages ∼ 74 mW/m2. Taking this value to represent the regional or background heat flow, a simple 2D analytical method was used to calculate theoretical heat flow variations due to topography. Across the mound, most of the variability is explained by topographic effects, including a local 6 mW/m2 negative anomaly over the central mound and a large 20 mW/m2 positive anomaly over the mound steep side slope. However, just south of the mound, there is a 6–7 mW/m2 positive anomaly in a 2-km-long band that has predominantly flat seafloor. Most of this anomaly is probably unrelated to topographic effects, but rather likely due to warm upward fluid flow along faults or fracture zones.

Description: Large uncertainties about the energy resource potential and role in global climate change of gas hydrates result from uncertainty about how much hydrate is contained in marine sediments. During Leg 204 of the Ocean Drilling Program (ODP) to the accretionary complex of the Cascadia subduction zone, we sampled the gas hydrate stability zone (GHSZ) from the seafloor to its base in contrasting geological settings defined by a 3D seismic survey. By integrating results from different methods, including several new techniques developed for Leg 204, we overcome the problem of spatial under-sampling inherent in robust methods traditionally used for estimating the hydrate content of cores and obtain a high-resolution, quantitative estimate of the total amount and spatial variability of gas hydrate in this structural system. We conclude that high gas hydrate content (30–40% of pore space or 20–26% of total volume) is restricted to the upper tens of meters below the seafloor near the summit of the structure, where vigorous fluid venting occurs. Elsewhere, the average gas hydrate content of the sediments in the gas hydrate stability zone is generally 〈2% of the pore space, although this estimate may increase by a factor of 2 when patchy zones of locally higher gas hydrate content are included in the calculation. These patchy zones are structurally and stratigraphically controlled, contain up to 20% hydrate in the pore space when averaged over zones ∼10 m thick, and may occur in up to ∼20% of the region imaged by 3D seismic data. This heterogeneous gas hydrate distribution is an important constraint on models of gas hydrate formation in marine sediments and the response of the sediments to tectonic and environmental change.

Description: During the Indian National Gas Hydrate Program (NGHP) Expedition 01, a series of well logs were acquired at several sites across the Krishna–Godavari (KG) Basin. Electrical resistivity logs were used for gas hydrate saturation estimates using Archie’s method. The measured in situ pore-water salinity, seafloor temperature and geothermal gradients were used to determine the baseline pore-water resistivity. In the absence of core data, Arp’s law was used to estimate in situ pore-water resistivity. Uncertainties in the Archie’s approach are related to the calibration of Archie coefficient (a), cementation factor (m) and saturation exponent (n) values. We also have estimated gas hydrate saturation from sonic P-wave velocity logs considering the gas hydrate in-frame effective medium rock-physics model. Uncertainties in the effective medium modeling stem from the choice of mineral assemblage used in the model. In both methods we assume that gas hydrate forms in sediment pore space. Combined observations from these analyses show that gas hydrate saturations are relatively low (〈5% of the pore space) at the sites of the KG Basin. However, several intervals of increased saturations were observed e.g. at Site NGHP-01-03 (Sh = 15–20%, in two zones between 168 and 198 mbsf), Site NGHP-01-05 (Sh = 35–38% in two discrete zone between 70 and 90 mbsf), and Site NGHP-01-07 shows the gas hydrate saturation more than 25% in two zones between 75 and 155 mbsf. A total of 10 drill sites and associated log data, regional occurrences of bottom-simulating reflectors from 2D and 3D seismic data, and thermal modeling of the gas hydrate stability zone, were used to estimate the total amount of gas hydrate within the KG Basin. Average gas hydrate saturations for the entire gas hydrate stability zone (seafloor to base of gas hydrate stability), sediment porosities, and statistically derived extreme values for these parameters were defined from the logs. The total area considered based on the BSR seismic data covers ∼720 km2. Using the statistical ranges in all parameters involved in the calculation, the total amount of gas from gas hydrate in the KG Basin study area varies from a minimum of ∼5.7 trillion-cubic feet (TCF) to ∼32.1 TCF.

Description: The Gas Hydrate Research and Development Organization (GHDO) of Korea successfully accomplished both coring (hydraulic piston and pressure coring) and logging (logging-while-drilling, LWD, and wireline logging) to investigate the presence of gas hydrate during the first deep drilling expedition in the Ulleung Basin, East Sea of Korea (referred to as UBGH1) in 2007. The LWD data from two sites (UBGH1-9, UBGH1-10) showed elevated electrical resistivity (〉80 Ω-m) and P-wave velocity (〉2000 m/s) values indicating the presence of gas hydrate. During the coring period, the richest gas hydrate accumulation was discovered at these intervals. Based on log data, the occurrence of gas hydrate is primarily controlled by the presence of fractures. The gas hydrate saturation calculated using Archie’s relation shows greater than 60% (as high as ∼90%) of the pore space, although Archie’s equation typically overestimates gas hydrate saturation in near-vertical fractures. The saturation of gas hydrate is also estimated using the modified Biot-Gassmann theory (BGTL) by Lee and Collett (2006). The saturation values estimated rom BGTL are much lower than those calculated from Archie’s equation. Based on log data, the hydrate-bearing sediment section is approximately 70 m (UBGH1-9) to 130 m (UBGH1-10) in thickness at these two sites. This was further directly confirmed by the recovery of gas hydrate samples and pore water freshening collected from deep drilling core during the expedition. LWD data also strongly support the interpretation of the seismic gas hydrate indicators (e.g., vent or chimney structures and bottom-simulating reflectors), which imply the probability of widespread gas hydrate presence in the Ulleung Basin.

Description: Supplies of conventional natural gas and oil are declining fast worldwide, and therefore new, unconventional forms of energy resources are needed to meet the ever-increasing demand. Amongst the many different unconventional natural resources are gas hydrates, a solid, ice-like crystalline compound of methane and water formed under specific low temperature and high pressure conditions. Gas hydrates are believed to exist in large quantities worldwide in oceanic regions of continental margins, as well as associated with permafrost regions in the Arctic. Some studies to estimate the global abundance of gas hydrate suggest that the total volume of natural gas locked up in form of gas hydrates may exceed all known conventional natural gas reserves, although large uncertainties exist in these assessments. Gas hydrates have been intensively studied in the last two decades also due to connections between climate forcing (natural and/or anthropogenic) and the potential large volumes of methane trapped in gas hydrate accumulations. The presence of gas hydrate within unconsolidated sediments of the upper few hundred meters below seafloor may also pose a geo-hazard to conventional oil and gas production. Additionally, climate variability and associated changes in pressure-temperature regimes and thus shifts in the gas hydrate stability zone may cause the occurrence of submarine slope failures. Several large-scale national gas hydrate programs exist especially in countries such as Japan, Korea, Taiwan, China, India, and New Zealand, where large demands of energy cannot be met by domestic supplies from natural resources. The past five years have seen several dedicated deep drilling expeditions and other scientific studies conducted throughout Asia and Oceania to understand gas hydrates off India, China, and Korea. This thematic set of publications is dedicated to summarize the most recent findings and results of geo-scientific studies of gas hydrates in the marginal seas and continental margin of the Asia, and Oceania region.

Description: High resolution multi-channel seismic reflection data acquired in the offshore domain of the Krishna–Godavari (KG) Basin, provide seismic indicators of gas hydrate including bottom-simulating reflectors (BSR), and high reflectivity zones below the BSR as an indication of free gas. The recent recovery of gas hydrate by drilling across various sites in the KG Basin, together with the seismic indicators, suggests overall widespread favourable gas hydrate stability conditions. The BSRs provide indirect estimates of geothermal gradients and heat flow. A simple conductive model is used to calculate heat flow from the depth of the BSR. BSR-derived heat flow values range from about 30 to 102 mW/m2 across the KG Basin. We also modeled the base of gas hydrate stability zone (BGHSZ) to analyse the linkage between gas hydrate occurrences in the basin and its relation to the fault systems in the basin. The influence of bottom water temperature, pore water salinity and geothermal gradient variation on gas hydrate stability zone (GHSZ) thickness is analysed in relation to both geological formations and tectonics. Our analysis suggests a highly variable GHSZ in the KG Basin controlled by local variations in the parameters of stability condition and faulting. The BSR-derived heat flow increases towards deeper water depths across the basin. Consistent local variations were also observed with low heat flow values over prominent topographic highs and high heat flow values over the flanks of the topographic highs. This variation may be due to focusing and defocusing effects of the topography. However, local BSR-depths patterns may also be due to dynamic effects, including the displacement of isotherms by thrust-faulting. Additionally, recent mass-wasting events and normal faulting resulted in the BGHSZ being too shallow relative to the predicted model (which assumed equilibrium conditions). Seismic and heat flow constraints from the Krishna-Godavari Basin gas hydrate system. Mar Geol - ResearchGate. Available from: http://www.researchgate.net/publication/223314484_Seismic_and_heat_flow_constraints_from_the_Krishna-Godavari_Basin_gas_hydrate_system._Mar_Geol [accessed Nov 5, 2015].

Description: Gas hydrate was discovered in the Krishna–Godavari (KG) Basin during the India National Gas Hydrate Program (NGHP) Expedition 1 at Site NGHP-01-10 within a fractured clay-dominated sedimentary system. Logging-while-drilling (LWD), coring, and wire-line logging confirmed gas hydrate dominantly in fractures at four borehole sites spanning a 500 m transect. Three-dimensional (3D) seismic data were subsequently used to image the fractured system and explain the occurrence of gas hydrate associated with the fractures. A system of two fault-sets was identified, part of a typical passive margin tectonic setting. The LWD-derived fracture network at Hole NGHP-01-10A is to some extent seen in the seismic data and was mapped using seismic coherency attributes. The fractured system around Site NGHP-01-10 extends over a triangular-shaped area of ∼2.5 km2 defined using seismic attributes of the seafloor reflection, as well as “seismic sweetness” at the base of the gas hydrate occurrence zone. The triangular shaped area is also showing a polygonal (nearly hexagonal) fault pattern, distinct from other more rectangular fault patterns observed in the study area. The occurrence of gas hydrate at Site NGHP-01-10 is the result of a specific combination of tectonic fault orientations and the abundance of free gas migration from a deeper gas source. The triangular-shaped area of enriched gas hydrate occurrence is bound by two faults acting as migration conduits. Additionally, the fault-associated sediment deformation provides a possible migration pathway for the free gas from the deeper gas source into the gas hydrate stability zone. It is proposed that there are additional locations in the KG Basin with possible gas hydrate accumulation of similar tectonic conditions, and one such location was identified from the 3D seismic data ˜6 km NW of Site NGHP-01-10.