Description: Natural gas hydrates are considered a potential resource for gas production on industrial scales. Gas hydrates contribute to the strength and stiffness of the hydrate-bearing sediments. During gas production, the geomechanical stability of the sediment is compromised. Due to the potential geotechnical risks and process management issues, the mechanical behavior of the gas hydrate-bearing sediments needs to be carefully considered. In this study, we describe a coupling concept that simplifies the mathematical description of the complex interactions occurring during gas production by isolating the effects of sediment deformation and hydrate phase changes. Central to this coupling concept is the assumption that the soil grains form the load-bearing solid skeleton, while the gas hydrate enhances the mechanical properties of this skeleton. We focus on testing this coupling concept in capturing the overall impact of geomechanics on gas production behavior though numerical simulation of a high-pressure isotropic compression experiment combined with methane hydrate formation and dissociation. We consider a linear-elastic stress-strain relationship because it is uniquely defined and easy to calibrate. Since, in reality, the geomechanical response of the hydrate-bearing sediment is typically inelastic and is characterized by a significant shear-volumetric coupling, we control the experiment very carefully in order to keep the sample deformations small and well within the assumptions of poroelasticity. The closely coordinated experimental and numerical procedures enable us to validate the proposed simplified geomechanics-to-flow coupling, and set an important precursor toward enhancing our coupled hydro-geomechanical hydrate reservoir simulator with more suitable elastoplastic constitutive models.

Description: Burial driven recycling is an important process in the natural gas hydrate (GH) systems worldwide, characterized by complex multiphysics interactions like gas migration through an evolving gas hydrate stability zone (GHSZ), competing gas-water-hydrate (i.e. fluid-fluid-solid) phase transitions, locally appearing and disappearing phases, and evolving sediment properties (e.g., permeability, reaction surface area, and capillary entry pressure). Such a recycling process is typically studied in homogeneous or layered sediments. However, there is mounting evidence that structural heterogeneity and anisotropy linked to normal and inclined fault systems or anomalous sediment layers have a strong impact on the GH dynamics. Here, we consider the impacts of such a structurally complex media on the recycling process. To capture the properties of the anomalous layers accurately, we introduce a fully mass conservative, high-order, discontinuous Galerkin (DG) finite element based numerical scheme. Moreover, to handle the rapidly switching thermodynamic phase states robustly, we cast the problem of phase transitions as a set of variational inequalities, and combine our DG discretization scheme with a semismooth Newton solver. Here, we present our new simulator, and demonstrate using synthetic geological scenarios, a) how the presence of an anomalous high-permeability layer, like a fracture or brecciated sediment, can alter the recycling process through flow-localization, and more importantly, b) how an incorrect or incomplete approximation of the properties of such a layer can lead to large errors in the overall prediction of the recycling process.

Description: The presence of gas hydrates (GHs) increases the stiffness and strength of marine sediments. In elasto‐plastic constitutive models, it is common to consider GH saturation (Sh) as key internal variable for defining the contribution of GHs to composite soil mechanical behavior. However, the stress‐strain behavior of GH‐bearing sediments (GHBS) also depends on the microscale distribution of GH and on GH‐sediment fabrics. A thorough analysis of GHBS is difficult, because there is no unique relation between Sh and GH morphology. To improve the understanding of stress‐strain behavior of GHBS in terms of established soil models, this study summarizes results from triaxial compression tests with different Sh, pore fluids, effective confining stresses, and strain histories. Our data indicate that the mechanical behavior of GHBS strongly depends on Sh and GH morphology, and also on the strain‐induced alteration of GH‐sediment fabrics. Hardening‐softening characteristics of GHBS are strain rate‐dependent, which suggests that GH‐sediment fabrics dynamically rearrange during plastic yielding events. We hypothesize that rearrangement of GH‐sediment fabrics, through viscous deformation or transient dissociation and reformation of GHs, results in kinematic hardening, suppressed softening, and secondary strength recovery, which could potentially mitigate or counteract large‐strain failure events. For constitutive modeling approaches, we suggest that strain rate‐dependent micromechanical effects from alterations of the GH‐sediment fabrics can be lumped into a nonconstant residual friction parameter. We propose simple empirical evolution functions for the mechanical properties and calibrate the model parameters against the experimental data. Plain Language Summary Gas hydrates (GHs) are crystalline‐like solids, which are formed from natural gas molecules and water at high pressure and low temperature. GHs, and particularly methane hydrates, are naturally abundant in marine sediments. It is known that the presence of GH increases the mechanical stiffness and strength of sediments, and there is strong effort in analyzing and quantifying these effects in order to understand potential risks of sediment destabilization or slope failure. Based on our experimental results from high‐pressure geotechnical studies, we show that not only the initial amount and distribution of GH are important for the increased strength of GH‐bearing sediments but also the dynamic rearrangement of GH‐sediment fabrics during deformation characterizes the stress‐strain response and enables strength recovery after failure. We propose that different microstructural mechanisms contribute to this rearrangement and strength recovery of GH sediment. However, we consider these complicated processes in a simplified manner in an improved numerical model, which can be applied for geotechnical risk assessment on larger scales.

Description: Abrupt fluid emissions from shallow marine sediments pose a threat to seafloor installations like wind farms and offshore cables. Quantifying such fluid emissions and linking pockmarks, the seafloor manifestations of fluid escape, to flow in the sub-seafloor remains notoriously difficult due to an incomplete understanding of the underlying physical processes. Here, using a compositional multi-phase flow model, we test plausible gas sources for pockmarks in the south-eastern North Sea, which recent observations suggest have formed in response to major storms. We find that the presence of free gas in the subsurface effectively damps storm wave-induced pressure changes due to its high compressibility, so that the mobilization of pre-existing gas pockets is unlikely. Rather, our results point to spontaneous appearance of a free gas phase via storm-induced gas exsolution from pore fluids. This mechanism is primarily driven by the pressure-sensitivity of gas solubility. We show that in highly permeable sediments containing gas-rich pore fluids, wave-induced pressure changes result in the appearance of a persistent gas phase. This suggests that seafloor fluid escape structures are not always proxies for overpressured shallow gas and that periodic seafloor pressure changes can induce persistent free gas phase to spontaneously appear.

Description: Offshore meteoric groundwater (OMG) has long been hypothesised to be a driver of seafloor geomorphic processes in continental margins worldwide. Testing this hypothesis has been challenging because of our limited understanding of the distribution and rate of OMG flow and seepage, and their efficacy as erosive/destabilising agents. Here we carry out numerical simulations of groundwater flow and slope stability using conceptual models and evolving stratigraphy - for passive siliciclastic and carbonate margin cases – to assess whether OMG and its evolution during a late Quaternary glacial cycle can generate the pore pressures required to trigger mechanical instabilities on the seafloor. Conceptual model results show that mechanical instabilities by OMG flow are most likely to occur in the outer shelf to upper slope, at or shortly before the Last Glacial Maximum sea level lowstand. Models with evolving stratigraphy show that OMG flow is a key driver of pore pressure development and instability in the carbonate margin case. In the siliciclastic margin case, OMG flow plays a secondary role in preconditioning the slope to failure. The higher degree of spatial/stratigraphic heterogeneity of carbonate margins, lower shear strengths of their sediments, and limited generation of overpressures by sediment loading may explain the higher susceptibility of carbonate margins, in comparison to siliciclastic margins, to mechanical instability by OMG flow. OMG likely played a more significant role in carbonate margin geomorphology (e.g. Bahamas, Maldives) than currently thought. Key Points Offshore meteoric groundwater (OMG) flow can drive mechanical instabilities in the outer shelf to upper slope Such instabilities occur at, or shortly after, the Last Glacial Maximum sea level lowstand Carbonate margins are more susceptible to mechanical instability by OMG than siliciclastic margins

Description: Subsurface flows, particularly hyporheic exchange fluxes, driven by streambed topography, permeability, channel gradient and dynamic flow conditions provide prominent ecological services such as nitrate removal from streams and aquifers. Stream flow dynamics cause strongly nonlinear and often episodic contributions of nutrient concentrations in river-aquifer systems. Using a fully coupled transient flow and reactive transport model, we investigated the denitrification potential of hyporheic zones during peak-flow events. The effects of streambed permeability, channel gradient and bedform amplitude on the spatio-temporal distribution of nitrate and dissolved organic carbon in streambeds and the associated denitrification potential were explored. Distinct peak-flow events with different intensity, duration and hydrograph shape were selected to represent a wide range of peak-flow scenarios. Our results indicated that the specific hydrodynamic characteristics of individual flow events largely determine the average positive or negative nitrate removal capacity of hyporheic zones, however the magnitude of this capacity is controlled by geomorphological settings (i.e. channel slope, streambed permeability and bedform amplitude). Specifically, events with longer duration and higher intensity were shown to promote higher nitrate removal efficiency with higher magnitude of removal efficiency in the scenarios with higher slope and permeability values. These results are essential for better assessment of the subsurface nitrate removal capacity under the influence of flow dynamics and particularly peak-flow events in order to provide tailored solutions for effective restoration of interconnected river-aquifer systems.

Description: Abrupt fluid emissions from shallow marine sediments pose a threat to seafloor installations like wind farms and offshore cables. Quantifying such fluid emissions and linking pockmarks, the seafloor manifestations of fluid escape, to flow in the sub-seafloor remains notoriously difficult due to an incomplete understanding of the underlying physical processes. Here, using a compositional multi-phase flow model, we test plausible gas sources for pockmarks in the south-eastern North Sea, which recent observations suggest have formed in response to major storms. We find that the mobilization of pre-existing gas pockets is unlikely because free gas, due to its high compressibility, damps the propagation of storm-induced pressure changes deeper into the subsurface. Rather, our results point to spontaneous appearance of a free gas phase via storm-induced gas exsolution from pore fluids. This mechanism is primarily driven by the pressure-sensitivity of gas solubility, and the appearance of free gas is largely confined to sediments in the vicinity of the seafloor. We show that in highly permeable sediments containing gas-rich pore fluids, wave-induced pressure changes result in the appearance of a persistent gas phase. This suggests that seafloor fluid escape structures are not always proxies for overpressured shallow gas and that periodic seafloor pressure changes can induce persistent free gas phase to spontaneously appear. Key Points - Storm-induced pressure changes can lead to spontaneous appearance of free gas phase near the seafloor - This process is driven by pressure-sensitive phase instabilities - This mechanism could help explain elusive gas sources in recently observed pockmarks in the North Sea

Description: Burial driven recycling is an important process in the natural gas hydrate (GH) systems worldwide, characterized by complex multiphysics interactions like gas migration through an evolving gas hydrate stability zone (GHSZ), competing gas-water-hydrate (i.e., fluid-fluid-solid) phase transitions, locally appearing and disappearing phases, and evolving sediment properties (e.g., permeability, reaction surface area, and capillary entry pressure). Such a recycling process is typically studied in homogeneous or layered sediments. However, there is mounting evidence that structural heterogeneity and anisotropy linked to normal and inclined fault systems or anomalous sediment layers have a strong impact on the GH dynamics. Here, we consider the impacts of such a structurally complex media on the recycling process. To capture the properties of the anomalous layers accurately, we introduce a fully mass conservative, high-order, discontinuous Galerkin (DG) finite element based numerical scheme. Moreover, to handle the rapidly switching thermodynamic phase states robustly, we cast the problem of phase transitions as a set of variational inequalities, and combine our DG discretization scheme with a semi-smooth Newton solver. Here, we present our new simulator, and demonstrate using synthetic geological scenarios, (a) how the presence of an anomalous high-permeability layer, like a fracture or brecciated sediment, can alter the recycling process through flow-localization, and more importantly, (b) how an incorrect or incomplete approximation of the properties of such a layer can lead to large errors in the overall prediction of the recycling process. Key Points Structural heterogeneity linked to inclined fault systems or anomalous sediment layers have a strong impact on the gas hydrate dynamics The presence of anomalous high-permeability layers within gas hydrate stability zone alters the recycling process through flow-localization The presented discontinuous Galerkin scheme is able to accurately capture the gas hydrate recycling processes through strongly anisotropic materials

Description: We conducted two‐dimensional numerical simulations to investigate the mechanisms underlying the strong spatiotemporal correlation observed between submarine landslides and gas hydrate dissociation due to glacial sea‐level drops. Our results suggest that potential plastic deformation or slip could occur at localized and small scales in the shallow‐water portion of the gas hydrate stability zone (GHSZ). This shallow‐water portion of the GHSZ typically lies within the area enclosed by three points: the BGHSZ–seafloor intersection, the seafloor at ∼600 m below sea level (mbsl), and the base of the GHSZ (BGHSZ) at ∼1,050 mbsl in low‐latitude regions. The deep BGHSZ (〉1,050 mbsl) could not slip; therefore, the entire BGHSZ was not a complete slip surface. Glacial hydrate dissociation alone is unlikely to cause large‐scale submarine landslides. Observed deep‐water (much greater than 600 mbsl) turbidites containing geochemical evidence of glacial hydrate dissociation potentially formed from erosion or detachment in the GHSZ pinch‐out zone. Plain Language Summary Many submarine landslides spatiotemporally correlate with gas hydrate dissociation. However, direct mechanical evidence supporting whether the overpressure and deformation due to glacial sea‐level drop‐induced hydrate dissociation are adequate for triggering submarine landslides is lacking. Here, we present two‐dimensional thermal‐hydraulic‐chemical and geomechanical models of a gas‐hydrate system in response to glacial sea‐level drops and conduct sensitivity analyses of the model behavior under a wide range of key conditions from a global perspective. Our simulations suggest that glacial hydrate dissociation might induce plastic deformation or slip at localized and small scales only possibly within the shallow‐water portion of the hydrate stability zone. The deep part (〉1,050 m below sea level) of the bottom boundary of the hydrate stability zone could not slip; therefore, the entire bottom boundary of the hydrate stability zone was not a complete slip surface. We demonstrate that glacial hydrate dissociation alone is unlikely to trigger large‐scale submarine landslides. Our work highlights the vicinity of the upper limit of the hydrate stability zone (where the base of the hydrate stability zone intersects the seafloor) as an important area for investigating overpressure and focused fluid flow, localized plastic deformation or slip, and downslope sediment transport related to glacial hydrate dissociation. Key Points Glacial hydrate dissociation might cause potential plastic deformation or slip at localized and small scales in shallow parts of the GHSZ The large deformation surface at the BGHSZ boundary of the potential plastic deformation zone was not a complete slip surface Glacial sea‐level drop‐induced gas hydrate dissociation alone is unlikely to have caused large‐scale submarine landslides