Bottom Simulating Reflector Research Articles

Gas hydrate accumulation in shelf break setting: Example from the Qiongdongnan Basin in the northern slope of the South China Sea

Based on interpretation and inversion of high‐resolution seismic data, a previously undocumented occurrence of pore‐filled gas hydrates that are associated with a submarine canyon system was discovered in the shelf break of the north‐east Qiongdongnan Basin (QDNB) in the northern South China Sea (SCS). Continuous bottom simulating reflections (BSRs) were recognized in the submarine striking NW‐SE ridges along the shelf break. The two‐way travel times (TWTs) of the BSRs are about 200–250 ms below the sea floor and the burial depths are mapped at ~180–225 mbsf after time‐depth conversion. The distribution of the enhanced reflections (ERs) and acoustic blanking zone overlying the BSRs indicate possible pore‐filled gas hydrate enrichment. ERs have been commonly observed immediately below the BSRs, indicating the accumulation of free gas, which exhibits distinct amplitude versus offset (AVO) anomalies. The acoustic impedance inversion results show that the gas hydrate saturation is 20–50%, and the thickness of the gas hydrate‐bearing layers is usually ~10–40 m, with a maximum thickness of up to ~90 m in local regions. Large‐scale listric faults originating from the deep sag connect the hydrocarbon source region and the shallow gas hydrate stability zone (GHSZ) and act as migration pathways for biogenic gas and potentially thermogenic gas. In addition, hydrocarbons may also be transported into the GHSZ along the sidewalls of the canyon. Canyon incision is believed to cause the increase in the burial depth of the BSRs on the flank of the ridge toward the base of the channel and thereby to trap the hydrocarbons migrating along the sidewall in permeable beds due to the decrease in permeability caused by the gas hydrate accumulation, resulting in the occurrence of ERs immediately below the BSRs. It is speculated that in the channel‐levee system, the coarse‐grained sediments should be deposited within the GHSZ, and the reservoir's physical properties are conducive to the wide distribution and enrichment of pore‐filled gas hydrates in the shelf break of the northern QDNB.

Reverse Time Migration of Vertical Cable Seismic Data to Image Hydrate-Bearing Sediments With High Resolution

Marine vertical cable seismic (VCS) is a promising survey technique for submarine complex structure imaging and reservoir monitoring, which uses vertical arrays of hydrophones deployed near the seafloor to record seismic wavefields in a quiet environment. Recently, we developed a new type of distributed VCS system for exploration and development of natural gas hydrates preserved in shallow sediments under the seafloor. Using this system and air-gun sources, we accomplished a 3D VCS yield data acquisition for gas hydrates exploration in the Shenhu area, South China Sea. In view of the characteristics of VCS geometry, we implement reverse time migration (RTM) on a common receiver gather to obtain high-resolution images of marine sediments. Due to the unique acquisition method, it is asymmetrical for the reflection path between the sources and the receivers in the VCS survey. Therefore, we apply accurate velocity analysis to common scatter point (CSP) gathers generated from common receiver gathers instead of the conventional velocity analysis based on common depth point gathers. RTM with this reliable velocity model results in high-resolution images of submarine hydrate-bearing sediments in deep water conditions. The RTM imaging section clearly shows the bottom simulating reflector (BSR) and also the reflection characteristics of the hydrate-bearing sediments filled with consolidated hydrates. Moreover, its resolution is relative to that of acoustic logging curves from the nearby borehole, and this imaging section is well consistent with the synthetic seismogram trace generated by the logging data. All these results reveal that VCS is a great potential technology for exploration and production of marine natural gas hydrates.

The hydrocarbon potential of the offshore Talara Basin, Peru

The offshore Talara Basin is the western extension of the hydrocarbon producing onshore fields since the mid-1800s area of Peru and is also located above the subduction zone of the active continental margin of South America. The offshore portion was evaluated using high quality 3D seismic where mapping horizons are all unconformities within the Eocene as well as the unconformities at the top Paleocene and top Cretaceous. Possible source rocks are the Cretaceous black marine shales of the Campanian Redondo Formation, the limestones of the Albian Muerto Formation, and the marine shales of the Paleogene. The primary target offshore is expected to be deep-water turbidites of Paleocene/Eocene age with a depositional source from the northeast from highlands created by the compressional uplift of the Andes. The main seals offshore are expected to be shales of the upper Eocene Lagunitos Formation and shales in the Chacra Formation, which are also seals in the onshore Litoral field. Thermal maturation modeling shows that two hydrocarbon kitchens exist in the offshore portion of the Talara basin, one in the north and one in the south. The probable Cretaceous source rocks reached the onset of maturity (VR = 0.63%) at a depth of 3,250 to 3,285 m (10,663 – 10,778 ft) between 30 and 39 Ma (Late Eocene to Oligocene). Importantly, the Cretaceous source rocks stay within the oil window once they enter it in the late Eocene. Satellite studies show a large offshore present-day oil seep in the southern part of the basin and 3D seismic shows direct hydrocarbon indicators (DHIs) imaged as flat spots and bottom simulating reflectors (BSR). Basin modeling suggests hydrocarbon migration pathways would have been updip (to the east) into the onshore traps and would therefore have first filled the offshore traps along the migration pathway. We conclude that the Talara Basin offshore offers excellent exploration opportunities in a proven productive area where multiple prospects have been mapped.

Seismic Attribute Analyses and Attenuation Applications for Detecting Gas Hydrate Presence

Identifying gas hydrates in the oceanic subsurface using seismic reflection data supported by the presence of a bottom simulating reflector (BSR) is not an easy task, given the wide range of geophysical methods that have been applied to do so. Though the presence of the BSR is attributed to the attenuation response, as seismic waves transition from hydrate-filled sediment within the gas hydrate stability zone (GHSZ) to free gas-bearing sediment below, few studies have applied a direct attenuation measurement. To improve the detection of gas hydrates and associated features, including the BSR and free gas accumulation beneath the gas hydrates, we apply a recently developed method known as Sparse-Spike Decomposition (SSD) that directly measures attenuation from estimating the quality factor (Q) parameter. In addition to performing attribute analyses using frequency attributes and a spectral decomposition method to improve BSR imaging, using a comprehensive analysis of the three methods, we make several key observations. These include the following: (1) low-frequency shadow zones seem to correlate with large values of attenuation; (2) there is a strong relationship between the amplitude strength of the BSR and the increase of the attenuation response; (3) the resulting interpretation of migration pathways of the free gas using the direct attenuation measurement method; and (4) for the data analyzed, the gas hydrates themselves do not give rise to either impedance or attenuation anomalies that fully differentiate them from nearby non-hydrate zones. From this last observation, we find that, although the SSD method may not directly detect in situ gas hydrates, the same gas hydrates often form an effective seal trapping and deeper free gas accumulation, which can exhibit a large attenuation response, allowing us to infer the likely presence of the overlying hydrates themselves.

Study on gas hydrate targets in the Danube Paleo-Delta with a dual polarization controlled-source electromagnetic system

The Black Sea is known to have extensive direct and indirect indicators of methane hydrates. Since the resistivity of the seafloor increases significantly in the presence of gas hydrates or free gas, marine controlled-source electromagnetics (CSEM) is a suitable method for the investigation of hydrates. We have collected CSEM data in a channel-levee system of the Danube paleo-delta (Bulgarian sector) at water depths of about 1500 m. The working area is within the gas hydrates stability zone and seismics suggests the presence of gas hydrates. The CSEM data were acquired with stationary receivers and a novel mobile, dual polarization transmitter system, and interpreted in terms of rotational invariants by means of 1-D inversions at common midpoints (CMP) to generate pseudo 2D resistivity sections. The inversion results reveal two resistive layers at shallow depths of 60–120 mbsf and greater depths of 270–400 mbsf. A comparison with seismics shows a good correlation of the shallow layer with high amplitude reflections and a velocity anomaly. The deeper layer can be tied to the current bottom simulating reflector (BSR). The comparison to a second CSEM experiment, which was conducted by the BGR (Federal Institute for Geosciences and Natural Resources, Germany) shows a good agreement of the derived sections down to a depth of about 350 m. Based on salinity and porosity models derived from boreholes in the Black Sea, we apply Archie's law to estimate potential gas hydrate saturations of up to 23% for the shallow resistor and up to 7% for the deeper layer. Differences are evident at depths greater than 450 m, where we see a conductive layer not evident in the BGR section. This deeper conductor could be evidence for increasing salinities of pore fluids at greater depths, which were previously found in DSDP drilling cores.

Q-compensated reverse-time migration for natural gas hydrate using fractional viscoacoustic wave equation

Natural gas hydrate (NGH) is a potential clean alternative energy source for fossil fuels. In seismic imaging profiles, NGH is often identified by the bottom-simulating reflection (BSR), which is characterized by strong reflection amplitude and negative polarity. High-resolution and amplitude-preserved seismic imaging are demanded for the detection of NGH. However, the traditional acoustic reverse-time migration ignores the attenuation characteristics of the medium, which leads to reduced amplitude and distorted phase of the seismic wave. In particular, when the hydrate saturation is low or the underlying formation does not contain free gas, it is difficult to observe the identifiable BSR in the traditional acoustic imaging profile, which causes difficulties in the identification of NGH. Here, we introduce the fractional viscoacoustic wave equation to perform the Q-compensated reverse-time migration (Q-RTM) for NGH, which can accurately recover the amplitude loss, correct the phase distortion, and provide high-resolution and high-illumination imaging results. Finally, Q-RTM can effectively enhance the BSR, reduce the uncertainty of hydrate identification, help to confirm the location and spatial distribution of the gas hydrate-bearing sediments, further refine the geological properties, and provide some theoretical basis for the exploitation and drilling of hydrate.

Diverse gas composition controls the Moby-Dick gas hydrate system in the Gulf of Mexico

AbstractIn marine basins, gas hydrate systems are usually identified by a bottom simulating reflection (BSR) that parallels the seafloor and coincides with the base of the gas hydrate stability zone (GHSZ). We present a newly discovered gas hydrate system, Moby-Dick, located in the Ship Basin in the northern Gulf of Mexico. In the seismic data, we observe a channel-levee complex with a consistent phase reversal and a BSR extending over an area of ∼14.2 km2, strongly suggesting the presence of gas hydrate. In contrast to classical observations, the Moby-Dick BSR abnormally shoals 150 m toward the seafloor from west to east, which contradicts the northward-shallowing seafloor. We argue that the likely cause of the shoaling BSR is a gradually changing gas mix across the basin, with gas containing heavier hydrocarbons in the west transitioning to methane gas in the east. Our study indicates that such abnormal BSRs can be controlled by gradual changes in the gas mix influencing the shape of the GHSZ over kilometers on a basin scale.

Improved imaging of gas hydrate reservoirs and their plumbing system using 2D elastic full-waveform inversion

Gas hydrates are ice-like crystalline materials that form under submarine environments of moderate pressure and low temperature. Another key factor to their formation is the abundance in gas supply from depth in addition to local biogenic gas. Detailed imaging and velocity analysis of the plumbing system of gas hydrates can provide confidence that amplitude anomalies in seismic data are related to gas hydrate accumulations. We have conducted 2D elastic full-waveform inversion (FWI) along a 14 km long segment of a 2D multichannel seismic profile to obtain a high-resolution velocity model of a hydrate system on the southern Hikurangi margin. We compare the FWI velocity model to previously published semblance- and tomography-based velocity models from the same data to explore how much more can be gained from the FWI. The FWI yielded a structurally more accurate velocity model that better delineated the low-velocity zone associated with free gas beneath the bottom simulating reflector (BSR) compared to the semblance- and tomography-based velocity models. Our results also find a lateral velocity inversion, that is, a narrow low-velocity zone surrounded by bands of higher velocities at a seaward-verging protothrust fault, which the two other methodologies failed to resolve. The FWI provides an improved lateral resolution making it an important tool when imaging the “plumbing” systems of gas hydrate reservoirs. In the southeastern limb of the anticline, our results find that the closely spaced landward-vergent protothrusts provide gas-charged fluids for hydrate formation above the BSR. Moreover, at the center of the anticline, our results find that a seaward-vergent protothrust fault appears to be acting as a conduit for gas-rich fluids into strata, although there is no accumulation of any significant hydrate above the BSR at the apex of the anticline. Our finding emphasizes the significance of densely spaced faults and fractures for providing gas for hydrate formation in the hydrate stability zone.

Estimation of Heat Flow Using a Bottom Simulating Reflection Based on 3D Seismic Data, West Africa

DOI:10.17014/ijog.8.3.297-311A Bottom Simulating Reflection (BSR) has been identified using 3D seismic data from the Cameroon continental slope margin. The BSR covers an area of c. 350 km2 in water depths ranging between 940 m to 1,750 m across an area characterized by high- and low-gradient slopes, gullies, scours, and fans. The thickness of the Gas Hydrate Stability Zone (GHSZ) is ~100 - 250 m, assuming an average velocity of 1,800 m/s. Pockmarks are intensively developed across the slope and most of them are observed in the BSR area. Geothermal estimation is based on hydrate stability conditions for pure methane - seawater system, hydrostatic pressure model, and a range of P-wave velocity models for the GHSZ, ranging from 1,600 -1,800 to 2,000 - 2,200 m/s. Geothermal gradient is calculated showing the range and distribution of thermal gradients in the BSR area from 0.046 C/m to 0.094 oC/m with assumed GHSZ velocity of 1,800 m/s. Thermal gradient anomalies have been observed in association with gullies, vertically stacked channels and in some individual pockmarks. The highest anomalies of 0.08 C/m – 0.094 C/m are found in the depression areas of pockmark trains, within seafloor gullies. These positive anomalies are most likely controlled by active or recently active fluid advection and expulsion through the Cameroon slope.

Seismic response of gas-hydrate zones: a case study from the Oman Sea

An offshore vintage 2D seismic dataset acquired in the Oman Sea recently has been reprocessed to achieve two objectives: mapping gas-hydrate zones and structural imaging beneath the bottom-simulating reflector (BSR). We present results of the analysis of the seismic data associated with three seismic lines in the project area. The sea-bottom water-saturated sediment column (WSZ) has an average thickness of 100 m and an interval velocity range of 1550–1700 m/s. The underlying gas-hydrate stability zone (GHSZ) has an average thickness of 180 m and an interval velocity range of 2000–2200 m/s. The BSR with its low-frequency, large-amplitude character is very prominent in the seismic images derived from prestack time migration (PSTM). The BSR defines the base of the GHSZ and thus is a crucial indicator of the presence of a gas-hydrate zone. Velocities close to the top of the GHSZ interval are slightly lower at some locations along the line traverses – indicative of a transitional boundary between the WSZ and GHSZ. The overthrust tectonics during the Mid-Miocene to Pliocene caused by the counterclockwise rotation of the Arabian Plate is evident in the seismic images. The resulting thrust faults may facilitate upward migration of fluids and gas, thus forming the free-gas zone beneath the BSRs.

A Rapid Numerical Method to Constrain 2D Focused Fluid Flow Rates Along Convergent Margins Using Dense BSR‐Based Temperature Field Data

AbstractEstimates of the sub‐seabed fluid flow rates are important for understanding hydrological budgets, biogeochemical cycles, and physical properties of the sediments. Fluid flow rates and directions, however, are difficult to measure, particularly beneath the seafloor. We developed a rapid method to estimate regional fluid migration rates using an extensive database of seismic reflection profiles taken offshore SW Taiwan. We observe bottom‐simulating reflector (BSR) that deflects toward the seafloor near thrust faults that indicate localized heat flow variations. At these sites, advecting warm pore fluids transport heat to shallower depths and force the BSR shallower. Our 2D steady‐state numerical method quantifies the fluid flow rates required to cause such thermal anomalies. We found that fluid flow rates near the trench of the accretionary wedge range between 0.1 and 16 m3 yr−1 m−1, with slower and faster rates generally associated with slope basin discontinuities and faults, respectively. To evaluate the fluid pattern evolution from subduction to collision, we studied three transects: one along the Manila subduction zone in the south and two in Taiwan's initial collision zone in the north. We quantified the fluid budget and partitioning of fluid flow between focused discharge through faults and diffusive flow through the wedge. Faults in Taiwan's accretionary wedge capture on average 25% of the total dewatering flux in the younger subduction zone and 38.5% in the tectonically mature collision zone. Our method provides estimates of fluid migration rates along convergent plate boundaries, and contributes to our understanding of focused fluid flow processes in many other regions.

Amplitude-Preserved Wave Equation: An Example to Image the Gas Hydrate System

Natural gas hydrate is an important energy source. Therefore, it is extremely important to provide a clear imaging profile to determine its distribution for energy exploration. In view of the problems existing in conventional migration methods, e.g., the limited imaging angles, we proposed to utilize an amplitude-preserved one-way wave equation migration based on matrix decomposition to deal with primary and multiple waves. With respect to seismic data gathered at the Chilean continental margin, a conventional processing flow to obtain seismic records with a high signal-to-noise ratio is introduced. Then, the imaging results of the conventional and amplitude-preserved one-way wave equation migration methods based on primary waves are compared, to demonstrate the necessity of implementing amplitude-preserving migration. Moreover, a simple two-layer model is imaged by using primary and multiple waves, which proves the superiority of multiple waves in imaging compared with primary waves and lays the foundation for further application. For the real data, the imaging sections of primary and multiple waves are compared. We found that multiple waves are able to provide a wider imaging illumination while primary waves fail to illuminate, especially for the imaging of bottom simulating reflections (BSRs), because multiple waves have a longer travelling path and carry more information. By imaging the actual seismic data, we can make a conclusion that the imaging result generated by multiple waves can be viewed as a supplementary for the imaging result of primary waves, and it has some guiding values for further hydrate and in general shallow gas exploration.

Heat Flow Measurements at the Danube Deep-Sea Fan, Western Black Sea

Seafloor heat flow measurements are utilized to determine the geothermal regime of the Danube deep-sea fan in the western Black Sea and are presented in the larger context of regional gas hydrate occurrences. Heat flow data were collected across paleo-channels in water depths of 550–1460 m. Heat flow across levees ranges from 25 to 30 mW m−2 but is up to 65 mW m−2 on channel floors. Gravity coring reveals sediment layers typical of the western Black Sea, consisting of three late Pleistocene to Holocene units, notably red clay within the lowermost unit cored. Heat flow derived from the bottom-simulating reflector (BSR), assumed to represent the base of the gas hydrate stability zone (GHSZ), deviates from seafloor measurements. These discrepancies are linked either to fast sedimentation or slumping and associated variations in sediment physical properties. Topographic effects account of up to 50% of heat flow deviations from average values. Combined with climate-induced variations in seafloor temperature and sea-level since the last glacial maximum large uncertainties in the prediction of the base of the GHSZ remain. A regional representative heat flow value is ~30 mW m−2 for the study region but deviations from this value may be up to 100%.

Дегазация Земли в Арктике: комплексный анализ факторов мощной эмиссии газа в море Лаптевых

Interpretation was conducted for 28 CDP seismic time sections with total length of 5930 km acquired by JSC “MAGE” in the Central Laptev Area, where a zone of powerful gas emission had been discovered earlier. 519 anomalous objects were revealed in near-bottom deposits with an average distance on seismic lines of 11,4 km, potentially connected with accumulations of gas and its migration paths. As a result of comprehensive analysis, for the first time, connection of gas seeps with deep-seated faults in the study area was justified. Highly likely forecast was made that in the area of the discovered seeps (seafloor depths from 50—60 m to 110 m), permafrost and gas hydrates are absent, and the seeps are caused by direct migration of gas from great depths. On the continental slope of the Laptev Sea, a bottom simulating reflector (BSR) was distinguished in CDP seismic sections, associated with the base of gas hydrates.

Velocity modeling and attribute analysis to understand the gas hydrates and free gas system in the Mannar Basin, India

The Mannar basin, located between India and Sri Lanka, is poorly explored for hydrocarbon compared to other petroliferous basins in the Indian offshore. 2D multichannel seismic data was acquired onboard R/V Samudra Ratnakar for exploring gas hydrate deposits in this basin. The bottom-simulating reflectors (BSRs) are observed in the seismic section, indicating the possible gas hydrate deposit. Several pull-down structures are identified below the BSR, suggesting the presence of free gas-bearing sediments. We also observed some chimney-like features, which extend from a deeper horizon up to the seafloor. An automated version of migration velocity analysis was used to estimate P-wave velocities along selected seismic lines. High interval velocities (>1.7 km/s) are observed above the BSR for ~100 ms TWT across the chimney locations, indicating the presence of gas hydrates for the first time in the Mannar basin. Low velocities are observed below the BSR, indicating a widespread region of free gas-bearing sediments. Seismic attributes (instantaneous amplitude and frequency) analysis are used to demarcate the extent of these regions. Chimneys are identified as vertical zones of reduced amplitude, while the gas-bearing sediments below the BSR show high amplitudes and low-frequency. A deeper horizon ~3.3 s TWT shows signatures associated with a deep-seated gas reservoir. We propose that chimneys and small sub-vertical fractures provide pathways for the migration of deep-seated gases to shallower regions. The hydrate deposits of the Mannar basin are unique as the organic-rich sediments are not supplied through any major rivers. Instead, they are derived from the high biomass due to the unusual physiography of the Mannar basin.

Gas hydrate identification and distribution in the Mexican Ridges, Mexico

We have mapped gas hydrates, free gas, and bottom-simulating-reflector (BSR) distributions in an area of the Mexican Ridges, central Gulf of Mexico, Mexico, revealing the relationship between these three elements and the tectonostratigraphy. The three elements are more visible when the host rock is a high-porosity sandstone because there is a large seismic impedance contrast between the solid gas hydrates above and the free gas below, which manifests itself on the seismic as a BSR. Gas hydrates are identified in the well as higher resistivity sandstone layers with a strong positive amplitude. When the host rock has a higher shale content with lower porosity, the impedance contrast is lower and the BSR is weak or not visible. The study area in Mexican Ridges is an anticline where gas hydrates and free gas are trapped on the crest after migrating through the dipping layers and faults from synclines that are generated in calcareous shale. The main seal is mass transport complex deposits from the Pliocene; when they are not deposited at the crest of the anticline, there is gas escape to the seafloor in the form of a gas chimney. In this way, we established a complete petroleum system for gas hydrates and free gas on the Mexican Ridges.

Identification of gas hydrates and bottom-simulating reflectors in far-offset seismic images

The presence of marine gas hydrates is routinely inferred based on the identification of bottom simulating reflectors (BSRs) in common-depth-point (CDP) seismic images. Additional seismic studies such as amplitude variation with offset (AVO) analysis can be applied for corroboration. Although confirmation is needed by drilling and sampling, seismic analysis has proven to be a cost-effective approach to identify the presence of marine gas hydrates. Single-channel far-offset seismic images are investigated for what appears to be a more reliable and cost-effective indicator for the presence of BSRs than traditional CDP processing or AVO analysis. A nontraditional approach for processing seismic data is taken to be more relevant to imaging the gas/gas hydrate contact. Instead of applying the traditional CDP seismic processing workflows from the oil industry, we more carefully review the significant amount of information existing in the data to explore how the character of the data changes as offset angle increases. Three cases from different environments are selected for detailed analysis. These include (1) stratigraphy running parallel with the ocean bottom, (2) a potential BSR, running parallel to the ocean bottom, and cross-cutting dipping reflections, and (3) a suspected thermal intrusion without a recognizable BSR. This investigation considers recently collected multichannel seismic data from the deep waters of the central Aleutian Basin beneath the Bering Sea, the preprocessing of the data sets, and the methodology for processing and display to generate single-channel seismic images. Descriptions are provided for the single-channel near- and far-offset seismic images for the example cases. Results indicate that BSRs related to marine gas hydrates, and originating due to the presence of free gas, are more easily and uniquely identifiable from single-channel displays of far-offset seismic images than from traditional CDP displays.

How gas-hydrate saturation and morphology control seismic attenuation: A case study from the south Hydrate Ridge

Prior studies have shown an ambiguous relationship between gas-hydrate saturation and seismic attenuation in different regions, but the effect of gas-hydrate morphology on seismic attenuation of hydrate-bearing sediments was often overlooked. We have combined seismic data with rock-physics modeling to elucidate how gas-hydrate saturation and morphology may control seismic attenuation. To extract P-wave attenuation, we process the vertical seismic profile data within a frequency range of 30–150 Hz and the sonic logging data within 10–15 kHz from three wells in the south Hydrate Ridge, offshore Oregon (USA), collected during Ocean Drilling Program Leg 204 in 2000. We calculate the P-wave attenuation using spectral matching and centroid frequency shift methods, and we use Archie’s relationship to derive gas-hydrate saturation from the resistivity data above the bottom-simulating reflection at the same wells. To interpret the observed seismic attenuation in terms of the effects of the gas-hydrate saturation and the morphology, we use the hydrate-bearing effective sediment rock-physics model. By comparing the observed and model-predicted attenuation values, we infer that (1) seismic attenuation appears to not be dominated by any single factor — instead, its variation is likely governed by the gas-hydrate saturation and the morphology; (2) the relationship between seismic attenuation and gas-hydrate saturation varies with different hydrate morphologies; (3) the squirt flow, occurring at different compliances of adjacent pores driven by pressure gradients, may be responsible for the significantly large or small attenuation over a broad frequency range.

In-situ borehole temperature measurements confirm dynamics of the gas hydrate stability zone at the upper Danube deep sea fan, Black Sea

Coring, geophysical logging, and in-situ temperature measurements were performed with the MARUM-MeBo200 seafloor rig to characterize gas hydrate occurrences in sediments of the Danube deep sea fan, off Romania, Black Sea. The new drilling data showed no evidence for significant gas hydrate saturations within the sediments but the presence of free gas at the depth of the bottom-simulating reflector (BSR). In-situ temperature and core-derived geochemical data suggest that the current base of the gas hydrate stability zone (BGHSZ) is ∼20 m shallower than the BSR. Investigation of the seismic data around the drill sites shows several locations where free gas previously trapped at a former BGHSZ migrated upwards forming a new reflection above the BSR. This shows that the gas hydrate system in the Danube deep sea fan is still responding to climate changes initiated at the end of the last glacial maximum.

Distribution of Methane Plumes on Cascadia Margin and Implications for the Landward Limit of Methane Hydrate Stability

Nearly 3,500 methane bubble streams, clustered into more than 1,300 methane emission sites, have been identified along the US Cascadia margin, derived both from archived published data and 2011, 2016–2018 dedicated multibeam surveys using co-registered seafloor and water column data. In this study, new multibeam sonar surveys systematically mapped nearly 40% of the US Cascadia margin, extending from the Strait of Juan de Fuca in the north to the Mendocino fracture zone in the south, and bounded East–West by the coast and the base of the accretionary prism. The frequency-depth histogram of the bubble emission sites has a dominant peak at the 500 m isobar, which extends laterally along much of the Cascadia margin off Oregon and Washington. Comparisons with published seismic data on the distribution of bottom simulating reflectors (BSR), which is the acoustic impedance boundary between methane hydrate (solid phase) and free gas phase below, correlates the bottom simulating reflectors upward termination of the feather edge of methane hydrate stability (FEMHS) zone and the newly identified bubble emission sites off Oregon and Washington. The Cascadia margin off northern California, where the BSR ends seaward of the FEMHS, has fewer sites centered on the 500 m isobaths, although data are more limited there. We propose that the peak in bubble emission sites observed near the 500 m isobath results from migration of free gas from beneath the solid phase of the BSR upslope to the FEMHS termination zone, and suggest that this boundary will be useful to monitor for a change in methane release rate potentially related to a warming ocean.