Description: A bottom-simulating reflector (BSR) occurs west of Svalbard in water depths exceeding 600 m, indicating that gas hydrate occurrence in marine sediments is more widespread in this region than anywhere else on the eastern North Atlantic margin. Regional BSR mapping shows the presence of hydrate and free gas in several areas, with the largest area located north of the Knipovich Ridge, a slow-spreading ridge segment of the Mid Atlantic Ridge system. Here, heat flow is high (up to 330 mW m-2), increasing towards the ridge axis. The coinciding maxima in across-margin BSR width and heat flow suggest that the Knipovich Ridge influenced methane generation in this area. This is supported by recent finds of thermogenic methane at cold seeps north of the ridge termination. To evaluate the source rock potential on the western Svalbard margin, we applied 1D petroleum system modeling at three sites. The modeling shows that temperature and burial conditions near the ridge were sufficient to produce hydrocarbons. The bulk petroleum mass produced since the Eocene is at least 5 kt and could be as high as ~0.2 Mt. Most likely, source rocks are Miocene organic-rich sediments and a potential Eocene source rock that may exist in the area if early rifting created sufficiently deep depocenters. Thermogenic methane production could thus explain the more widespread presence of gas hydrates north of the Knipovich Ridge. The presence of microbial methane on the upper continental slope and shelf indicates that the origin of methane on the Svalbard margin varies spatially.

Description: Highlights • BSR position does not match BGHS as predicted based on regional TP conditions. • Use steady state and transient models to determine extent of hydrate stability. • Investigate the influence of topographic focusing on hydrate stability. • Variable thermal properties of sediment impact hydrate stability. The Danube Fan in the western Black Sea shows many features indicating the presence of gas and gas hydrates, including a bottom simulating reflection (BSR), high-amplitude anomalies beneath the BSR and the presence of gas flares at the seafloor. The BSR depth derived from 3D P-cable seismic data of an older slope canyon of the fan (the S2 canyon) suggests that the BSR is not in equilibrium with the present-day topography. The Danube Fan was abandoned ∼7.5 ka, and the S2 canyon was likely incised ∼20 ka, suggesting that the gas hydrate system has had at least 7.5 ka years to equilibrate to the present-day conditions. Here we examine the extent and position of the hydrate stability zone through constructing both steady and transient state models of a 2D profile across the S2 canyon. This was done using inputs from mapping of the 3D P-cable seismic data and geochemical analysis of core samples. Using these models, we investigate the effects of different factors including variable thermal properties of heterogeneous sediments in the vicinity of the canyon and, topographic focusing on the geothermal gradient on the extent of the hydrate stability zone. Our results indicate that both factors have a significant effect and that the hydrate system may actually be in, or approaching equilibrium.

Description: A bottom-simulating reflector (BSR) occurs west of Svalbard in water depths exceeding 600 m, indicating that gas hydrate occurrence in marine sediments is more widespread in this region than anywhere else on the eastern North Atlantic margin. Regional BSR mapping shows the presence of hydrate and free gas in several areas, with the largest area located north of the Knipovich Ridge, a slow-spreading ridge segment of the Mid Atlantic Ridge system. Here, heat flow is high (up to 330 mW m-2), increasing towards the ridge axis. The coinciding maxima in across-margin BSR width and heat flow suggest that the Knipovich Ridge influenced methane generation in this area. This is supported by recent finds of thermogenic methane at cold seeps north of the ridge termination. To evaluate the source rock potential on the western Svalbard margin, we applied 1D petroleum system modeling at three sites. The modeling shows that temperature and burial conditions near the ridge were sufficient to produce hydrocarbons. The bulk petroleum mass produced since the Eocene is at least 5 kt and could be as high as ~0.2 Mt. Most likely, source rocks are Miocene organic-rich sediments and a potential Eocene source rock that may exist in the area if early rifting created sufficiently deep depocenters. Thermogenic methane production could thus explain the more widespread presence of gas hydrates north of the Knipovich Ridge. The presence of microbial methane on the upper continental slope and shelf indicates that the origin of methane on the Svalbard margin varies spatially.

Description: Gas migration pathways in the Gulf of Mexico are strongly influenced by the extensive formation and time evolution of salt canopies, welds and sheets. This multi-level salt system (known as the Louann Salt formation) deposited mostly within Callovian age (upper Middle Jurassic) and mobilized during late Miocene up to Pliocene-Pleistocene times controls the extension and direction of petroleum components migration over the entire history of the basin which, in return, has a major impact on potential gas transportation into the gas hydrate stability zone (GHSZ). In the context of gas hydrate formation, presence of extensive salt deposits tends to bend gas migration pathways from vertical (typical for the Gulf of Mexico region) towards rather horizontal and dispersed. However, amalgamation of two or more salt structures often results in re-focusing of the flow towards the local topographic subsalt heights. Together with the formation of local sediment discontinuity structures such as faults developing at the rims and tops of rootless salt deposits related to further stages of allochthonous salt mobilization, new high-permeability migration pathways develop and act as direct connection for the thermogenic gas to the GHSZ. Our study presents the 3D modeling solution quantifying and exploring the gas hydrate accumulation potential in the marine environment experiencing salt tectonics such as the Green Canyon, Gulf of Mexico. This modeling study evaluates the potential of bio- and thermogenic gas hydrate formation within Pliocene-Pleistocene reservoir layers based on full basin re-construction which accounts for depositional and transient thermal history of the basin, source rock maturation, petroleum generation, expulsion and migration, salt tectonics and associated faults development. Based on a numerical study calibrated with the existing field data, we present a new distribution pattern of gas hydrates attributed to both microbial and thermogenic origin. We present here an explanation for a formation mechanism of large gas hydrate amounts (〉 70 vol. %) wide-spread at the base of the stability zone as a result of the gas hydrate-free gas recycling process enhanced by very high Neogene sedimentation rates in the region. We suggest that the rapid development of secondary intra-salt mini-basins situated on top of the allochthonous salt deposits and following sediment subsidence caused a consequent dislocation of the GHSZ lower boundary and led to efficient gas hydrate dissociation process followed by a free gas re-charge into the GHSZ.

Description: Introduction We have developed a new multi 1-D numerical model to investigate and understand the processes of gas hydrate formation and dissolution in anoxic marine sediments under a wide range of conditions. By this reaction-transport model we are able to investigate a various aspects of gas hydrate dynamics: sediment compaction which results in expulsion of pore fluids containing various chemical species, reduction in porosity and permeability of the sediment matrix due to hydrate formation, time-resolved evolution of pressure and temperature regimes, multiphase flow of compressible pore fluids, gas hydrate, and a free gas, thermal-blanketing effect due to vigorous sedimentation of cold impermeable layers, gas hydrate dissolution as a response to a slowing down sedimentation, and the effects of salinity variations on the thickness of the Gas Hydrate Stability Zone (GHSZ). Numerical model The reaction-transport model contains various chemical compounds (solid organic carbon, dissolved in pore water methane, dissolved inorganic carbon, dissolved sulfates, gas hydrates, and free methane gas). We consider a reference frame which extends from the seafloor to the bottom of the GHSZ (defined as a combination of pressure, temperature, and salinity conditions) plus 50m of Free Gas Zone lying directly beneath. However, the upper part of sediment column (10 cm) is not considered in the model due to strong bioturbation processes which might potentially have an impact on the gradients of dissolved chemical species. Initially, the system is filled by compressible pore fluids of a given salinity (consistent with a value at the sediment-water interface). As the upper boundary conditions, we have applied constant concentrations of dissolved methane, dissolved inorganic carbon, and sulfate according to the mean values in the ocean. At the beginning of each time-step, a new sediment layer is deposited at the top of sediment column according to a given sedimentation rate, lithological type, and initial porosity at the surface. Transport processes have been split into the advection and diffusion part and solved separately for every chemical compound. Multiphase flow of dissolved chemical species and free gas phase has been solved by finite-volumes method according to the Darcy’s law. Molecular diffusion of dissolved species is controlled by changes in concentration gradients and has been solved by finite-elements method. Reaction module contains kinetically controlled rates of methanogenesis, sulfate reduction, methane oxidation, and POC degradation. POC decay via microbial sulfate reduction takes place until the dissolved sulfate pool in ambient pore waters is depleted. Below the sulfate penetration depth, POC is microbially decomposed into methane and CO2. Upward diffusing dissolved methane is consumed by anaerobic oxidation within the sulfate-methane transition zone. This reaction module has been evaluated previously by Wallmann et al., 2006, Marquardt et al., 2010, and Burwicz et al., 2011. Applications Dynamic un-steady state compaction allows us to investigate gas hydrate formation and dissolution in terms of changing parameters (e.g. sedimentation rate or permeability of deposited sediments). By depositing sediment layers of a different grain size (‘sandwich-like’ scenario), we have observed that lithology of potential hydrate-bearing layers (e.g. coarse-grained sands vs. shales) results in preferential hydrate accumulation in the first ones which stays in agreement with field observations. We have also investigated the effect of slowing down sedimentation rates on gas hydrate dissolution. We have concluded that slow deposition of sediment layers at the top of sediment column and, as a result, a decrease in POC input in time, result in undersaturated in CH4 pore waters causing hydrate destabilization. This scenario clearly shows the importance of constraining a time-resolved sedimentation history in gas hydrate simulations which are coupled with climate models. By depositing thick layers of cold low-permeable sediments on top of the column, we have investigated the temperature variations within sediments, known as ‘thermal blanketing’ effect, which has an impact on previously formed hydrates.

Description: This study provides new estimates for the global offshore methane hydrate inventory formed due to microbial CH4 production under Quaternary and Holocene boundary conditions. A multi-1D model for particular organic carbon (POC) degradation, gas hydrate formation and dissolution is presented. The novel reaction-transport model contains an open three-phase system of two solid compounds (organic carbon, gas hydrates), three dissolved species (methane, sulfates, inorganic carbon) and one gaseous phase (free methane). The model computes time-resolved concentration profiles for all compounds by accounting for chemical reactions as well as diffusive and advective transport processes. The reaction module builds upon a new kinetic model of POC degradation which considers a down-core decrease in reactivity of organic matter. Various chemical reactions such as organic carbon decay, anaerobic oxidation of methane, methanogenesis, and sulfate reduction are resolved using appropriate kinetic rate laws and constants. Gas hydrates and free gas form if the concentration of dissolved methane exceeds the pressure, temperature, and salinity-dependent solubility limits of hydrates and/or free gas, with a rate given by kinetic parameters. Global input grids have been compiled from a variety of oceanographic, geological and geophysical data sets including a new parameterization of sedimentation rates in terms of water depth. We find prominent gas hydrate provinces offshore Central America where sediments are rich in organic carbon and in the Arctic Ocean where low bottom water temperatures stabilize methane hydrates. The world’s total gas hydrate inventory is estimated at 0.82 x 10sup13 m3 - 2.10 x 10sup15 m3 CH4 (at STP conditions) or, equivalently, 4.18–995 Gt of methane carbon. The first value refers to present day conditions estimated using the relatively low Holocene sedimentation rates; the second value corresponds to a scenario of higher Quaternary sedimentation rates along continental margins. Our results clearly show that in-situ POC degradation is at present not an efficient hydrate forming process. Significant hydrate deposits in marine settings are more likely to have formed at times of higher sedimentation during the Quaternary or as a consequence of upward fluid transport at continental margins.

Description: 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.

Description: We have developed a multi-1D numerical model of gas hydrate formation and dissolution processes in anoxic marine sediments and, by this model, we have estimated the new global gas hydrate inventory (BURWICZ E. B. et al., 2011). The reaction-transport model contains various chemical compounds (solid organic carbon, dissolved methane, inorganic carbon, and sulfates, gas hydrates, and free methane gas). The rates of POC degradation, anaerobic methane oxidation, sulfate reduction, and methanogenesis are kinetically controlled. Gas hydrate stability zone (GHSZ) is defined as a combination of pressure, temperature, and (to a smaller degree) salinity conditions. The lower boundary of the GHSZ is defined as the intersection of gas hydrate and methane gas solubilities. The diffusion equations are solved using a fully-implicit finite-differences method, while all transport processes are resolved by a Semi-Lagrangian scheme. Global input data sets (1°x1° resolution) were compiled from various oceanographic, geological and geophysical sources. The entire model was implemented in Matlab.

Description: Highlights • There is direct and indirect evidence for hydrate occurrence in several areas around Europe. • Hydrate is particularly widespread offshore Norway and Svalbard and in the Black Sea. • Hydrate occurrence often coincides with conventional thermogenic hydrocarbon provinces. • The regional abundance of hydrate in Europe is poorly known. Abstract Large national programs in the United States and several Asian countries have defined and characterised their marine methane hydrate occurrences in some detail, but European hydrate occurrence has received less attention. The European Union-funded project “Marine gas hydrate – an indigenous resource of natural gas for Europe” (MIGRATE) aimed to determine the European potential inventory of exploitable gas hydrate, to assess current technologies for their production, and to evaluate the associated risks. We present a synthesis of results from a MIGRATE working group that focused on the definition and assessment of hydrate in Europe. Our review includes the western and eastern margins of Greenland, the Barents Sea and onshore and offshore Svalbard, the Atlantic margin of Europe, extending south to the northwestern margin of Morocco, the Mediterranean Sea, the Sea of Marmara, and the western and southern margins of the Black Sea. We have not attempted to cover the high Arctic, the Russian, Ukrainian and Georgian sectors of the Black Sea, or overseas territories of European nations. Following a formalised process, we defined a range of indicators of hydrate presence based on geophysical, geochemical and geological data. Our study was framed by the constraint of the hydrate stability field in European seas. Direct hydrate indicators included sampling of hydrate; the presence of bottom simulating reflectors in seismic reflection profiles; gas seepage into the ocean; and chlorinity anomalies in sediment cores. Indirect indicators included geophysical survey evidence for seismic velocity and/or resistivity anomalies, seismic reflectivity anomalies or subsurface gas escape structures; various seabed features associated with gas escape, and the presence of an underlying conventional petroleum system. We used these indicators to develop a database of hydrate occurrence across Europe. We identified a series of regions where there is substantial evidence for hydrate occurrence (some areas offshore Greenland, offshore west Svalbard, the Barents Sea, the mid-Norwegian margin, the Gulf of Cadiz, parts of the eastern Mediterranean, the Sea of Marmara and the Black Sea) and regions where the evidence is more tenuous (other areas offshore Greenland and of the eastern Mediterranean, onshore Svalbard, offshore Ireland and offshore northwest Iberia). We provide an overview of the evidence for hydrate occurrence in each of these regions. We conclude that around Europe, areas with strong evidence for the presence of hydrate commonly coincide with conventional thermogenic hydrocarbon provinces.