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: Geochemical data (CH4, SO42−, I−, Cl−, particulate organic carbon (POC), δ13C-CH4, and δ13C-CO2) are presented from the upper 30 m of marine sediment on a tectonic submarine accretionary wedge offshore southwest Taiwan. The sampling stations covered three ridges (Tai-Nan, Yung-An, and Good Weather), each characterized by bottom simulating reflectors, acoustic turbidity, and different types of faulting and anticlines. Sulfate and iodide concentrations varied little from seawater-like values in the upper 1–3 m of sediment at all stations; a feature that is consistent with irrigation of seawater by gas bubbles rising through the soft surface sediments. Below this depth, sulfate was rapidly consumed within 5–10 m by anaerobic oxidation of methane (AOM) at the sulfate-methane transition. Carbon isotopic data imply a mainly biogenic methane source. A numerical transport-reaction model was used to identify the supply pathways of methane and estimate depth-integrated turnover rates at the three ridges. Methane gas ascending from deep layers, facilitated by thrusts and faults, was by far the dominant term in the methane budget at all sites. Differences in the proximity of the sampling sites to the faults and anticlines mainly accounted for the variability in gas fluxes and depth-integrated AOM rates. By comparison, methane produced in situ by POC degradation within the modeled sediment column was unimportant. This study demonstrates that the geochemical trends in the continental margins offshore SW Taiwan are closely related to the different geological settings.

Description: Current estimates suggest that more than 60% of the global seafloor are covered by millions of abyssal hills and mountains. These features introduce spatial fluid-dynamic granularity whose influence on deep-ocean sediment biogeochemistry is unknown. Here we compare biogeochemical surface-sediment properties from a fluid-dynamically well-characterized abyssal hill and upstream plain: (1) In hill sediments, organic-carbon and -nitrogen contents are only about half as high as on the plain while proteinaceous material displays less degradation; (2) on the hill, more coarse-grained sediments (reducing particle surface area) and very variable calcite contents (influencing particle surface charge) are proposed to reduce the extent, and influence compound-specificity, of sorptive organic-matter preservation. Further studies are needed to estimate the representativeness of the results in a global context. Given millions of abyssal hills and mountains, their integrative influence on formation and composition of deep-sea sediments warrants more attention.

Description: Takahe seep, located on the Opouawe Bank, Hikurangi Margin, is characterized by a well-defined subsurface seismic chimney structure ca. 80,500 m2 in area. Sub-seafloor geophysical data based on acoustic anomaly layers indicated the presence of gas hydrate and free gas layers within the chimney structure. Reaction-transport modeling was applied to porewater data from 11 gravity cores to constrain methane turnover rates and benthic methane fluxes in the upper 10 m. Model results show that methane dynamics were highly variable due to transport and dissolution of ascending gas. The dissolution of gas (up to 3761 mmol m−2 yr−1) dwarfed the rate of methanogenesis within the simulated sediment column (2.6 mmol m−2 yr−1). Dissolved methane is mainly consumed by anaerobic oxidation of methane (AOM) at the base of the sulfate reduction zone and trapped by methane hydrate formation below it, with maximum rates in the central part of the chimney (946 and 2420 mmol m−2 yr−1, respectively). A seep-wide methane budget was constrained by combining the biogeochemical model results with geophysical data and led to estimates of AOM rates, gas hydrate formation and benthic dissolved methane fluxes of 3.68 × 104 mol yr−1, 73.85 × 104 mol yr−1and 1.19 × 104 mol yr−1, respectively. A much larger flux of methane probably escapes in gaseous form through focused bubble vents. The approach of linking geochemical model results with spatial geophysical data put forward here can be applied elsewhere to improve benthic methane turnover rates from limited single spot measurements to larger spatial scales.

Description: Mud volcanoes are seafloor expressions of focused fluid flow that are common in compressional tectonic settings. New high-resolution 3-D seismic data from the Mercator mud volcano (MMV) and an adjacent buried mud volcano (BMV) image the internal structure of the top 800 m of sediment at both mud volcanoes, revealing that both are linked and have been active episodically. The total volumes of extruded mud range between 0.15 and 0.35 km3 and 0.02–0.05 km3 for the MMV and the BMV, respectively. The pore water composition of surface sediment samples suggests that halokinesis has played an important role in the evolution of the mud volcanoes. We propose that erosion of the top of the Vernadsky Ridge that underlies the mud volcanoes activated salt movement, triggering deep migration of fluids, dissolution of salt, and sediment liquefaction and mobilization since the end of the Pliocene. Since beginning of mud volcanism in this area, the mud volcanoes erupted four times while there was only one reactivation of salt tectonics. This implies that there are other mechanisms that trigger mud eruptions. The stratigraphic relationship of mudflows from the MMV and BMV indicates that the BMV was triggered by the MMV eruptions. This may either be caused by loading-induced hydrofracturing within the BMV or due to a common feeder system for both mud volcanoes. This study shows that the mud volcanoes in the El Arraiche mud volcano field are long-lived features that erupt with intervals of several tens of thousands of years.

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: Deep-sea mining for polymetallic nodules is expected to have severe environmental impacts because not only nodules but also benthic fauna and the upper reactive sediment layer are removed through the mining operation and blanketed by resettling material from the suspended sediment plume. This study aims to provide a holistic assessment of the biogeochemical recovery after a disturbance event by applying prognostic simulations based on an updated diagenetic background model and validated against novel data on microbiological processes. It was found that the recovery strongly depends on the impact type; complete removal of the reactive surface sediment reduces benthic release of nutrients over centuries, while geochemical processes after resuspension and mixing of the surface sediment are near the pre-impact state 1 year after the disturbance. Furthermore, the geochemical impact in the DISturbance and reCOLonization (DISCOL) experiment area would be mitigated to some degree by a clay-bound Fe(II)-reaction layer, impeding the downward diffusion of oxygen, thus stabilizing the redox zonation of the sediment during transient post-impact recovery. The interdisciplinary (geochemical, numerical and biological) approach highlights the closely linked nature of benthic ecosystem functions, e.g. through bioturbation, microbial biomass and nutrient fluxes, which is also of great importance for the system recovery. It is, however, important to note that the nodule ecosystem may never recover to the pre-impact state without the essential hard substrate and will instead be dominated by different faunal communities, functions and services.

Description: The thriving interest in harvesting deep-sea mineral resources, such as polymetallic nodules, calls for environmental impact studies, and ultimately, for regulations for environmental protection. Industrial-scale deep-sea mining of polymetallic nodules most likely has severe consequences for the natural environment. However, the effects of mining activities on deep-sea ecosystems, sediment geochemistry and element fluxes are still poorly conceived. Predicting the environmental impact is challenging due to the scarcity of environmental baseline studies as well as the lack of mining trials with industrial mining equipment in the deep sea. Thus, currently we have to rely on small-scale disturbances simulating deep-sea mining activities as a first-order approximation to study the expected impacts on the abyssal environment. Here, we investigate surface sediments in disturbance tracks of seven small-scale benthic impact experiments, which have been performed in four European contract areas for the exploration of polymetallic nodules in the Clarion-Clipperton Zone (CCZ). These small-scale disturbance experiments were performed 1 day to 37 years prior to our sampling program in the German, Polish, Belgian and French contract areas using different disturbance devices. We show that the depth distribution of solid-phase Mn in the upper 20 cm of the sediments in the CCZ provides a reliable tool for the determination of the disturbance depth, which has been proposed in a previous study (Paul et al., 2018). We found that the upper 5–15 cm of the sediments were removed during various small-scale disturbance experiments in the different exploration contract areas. Transient transport-reaction modelling for the Polish and German contract areas reveals that the removal of the surface sediments is associated with the loss of reactive labile organic carbon. As a result, oxygen consumption rates decrease significantly after the removal of the surface sediments, and consequently, oxygen penetrates up to tenfold deeper into the sediments inhibiting denitrification and Mn(IV) reduction. Our model results show that the post-disturbance geochemical re-equilibration is controlled by diffusion until the reactive labile TOC fraction in the surface sediments is partly re-established and the biogeochemical processes commence. While the re-establishment of bioturbation is essential, the geochemical re-equilibration of the sediments is ultimately controlled by the burial rates of organic matter. Hence, under current depositional conditions, the new geochemical equilibrium in the sediments of the CCZ is reached only on a millennia scale even for these small-scale disturbances simulating deep-sea mining activities.

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: Free gas migration through the gas hydrate stability zone (GHSZ) and subsequent gas seepage at the seabed are characteristic features in marine gas hydrate provinces worldwide. The biogenic or thermogenic gas is typically transported along faults from deeper sediment strata to the GHSZ. Several mechanisms have been proposed to explain free gas transport through the GHSZ. While inhibition of hydrate formation by elevated salinities and temperatures have been addressed previously in studies simulating unfocused, area-wide upward advection of gas, which is not adequately supported by field observations, the role of focused gas flow through chimney-like structures has been underappreciated in this context. Our simulations suggest that gas migration through the GHSZ is, fundamentally, a result of methane gas supply in excess of its consumption by hydrate formation. The required high gas flux is driven by local overpressure, built up from gas accumulating below the base of the GHSZ that fractures the overburden when exceeding a critical pressure, thereby creating the chimney-like migration pathway. Initially rapid hydrate formation raises the temperature in the chimney structure, thereby facilitating further gas transport through the GHSZ. As a consequence, high hydrate saturations form preferentially close to the seafloor, where temperatures drop to bottom water values, producing a prominent subsurface salinity peak. Over time, hydrates form at a lower rate throughout the chimney structure, while initial temperature elevation and salinity peak dissipate. Thus, our simulations suggest that the near-surface salinity peak and elevated temperatures are a result of transient high-flux gas migration through the GHSZ.