Description: The migration of methane through the gas hydrate stability zone (GHSZ) in the marine subsurface is characterized by highly dynamic reactive transport processes coupled to thermodynamic phase transitions between solid gas hydrates, free methane gas, and dissolved methane in the aqueous phase. The marine subsurface is essentially a water-saturated porous medium where the thermodynamic instability of the hydrate phase can cause free gas pockets to appear and disappear locally, causing the model to degenerate. This poses serious convergence issues for the general-purpose nonlinear solvers (e.g., standard Newton), and often leads to extremely small time-step sizes. The convergence problem is particularly severe when the rate of hydrate phase change is much lower than the rate of gas dissolution. In order to overcome this numerical challenge, we have developed an all-at-once Newton scheme tailored to our gas hydrate model, which can handle rate-based hydrate phase change coupled with equilibrium gas dissolution in a mathematically consistent and robust manner.

Description: The exchange of CH4 by CO2 in gas hydrates is of interest for the production of natural gas from methane hydrate with net zero climate gas balance, and for managing risks that are related to sediment destabilization and mobilization after gas-hydrate dissociation. Several experimental studies on the dynamics and efficiency of the process exist, but the results seem to be partly inconsistent. We used confocal Raman spectroscopy to map an area of several tens to hundreds µm of a CH4 hydrate sample during its exposure to liquid and gaseous CO2. On this scale, we could identify and follow different processes in the sample that occur in parallel. Next to guest-molecule exchange, gas-hydrate dissociation also contributes to the release of CH4. During our examination period, about 50% of the CO2 was bound by exchange for CH4 molecules, while the other half was bound by new formation of CO2 hydrates. We evaluated single gas-hydrate grains with confirmed gas exchange and applied a diffusion equation to quantify the process. Obtained diffusion coefficients are in the range of 10−13–10−18 m2/s. We propose to use this analytical diffusion equation for a simple and robust modeling of CH4 production by guest-molecule exchange and to combine it with an additional term for gas-hydrate dissociation.

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

Description: In 1964, exploration drilling in the German Sector of the North Sea hit a gas pocket at ∼2900 m depth below the seafloor and triggered a blowout, which formed a 550 m-wide and up to 38 m deep seafloor crater now known as Figge Maar. Although seafloor craters formed by fluid flow are very common structures, little is known about their formation dynamics. Here, we present 2D reflection seismic, sediment echosounder, and multibeam echosounder data from three geoscientific surveys of the Figge Maar blowout crater, which are used to reconstruct its formation. Reflection seismic data support a scenario in which overpressured gas ascended first through the lower part of the borehole and then migrated along steeply inclined strata and faults towards the seafloor. The focused discharge of gas at the seafloor removed up to 4.8 Mt of sediments in the following weeks of vigorous venting. Eyewitness accounts document that the initial phase of crater formation was characterized by the eruptive expulsion of fluids and sediments cutting deep into the substrate. This was followed by a prolonged phase of sediment fluidization and redistribution widening the crater. After fluid discharge ceased, the Figge Maar acted as a sediment trap reducing the crater depth to ∼12 m relative to the surrounding seafloor in 2018, which corresponds to an average sedimentation rate of ∼22,000 m 3 /yr between 1995 and 2018. Hydroacoustic and geochemical data indicate that the Figge Maar nowadays emits primarily biogenic methane, predominantly during low tide. The formation of Figge Maar illustrates hazards related to the formation of secondary fluid pathways, which can bypass safety measures at the wellhead and are thus difficult to control.