Search

Microbial cells buried in subseafloor sediments comprise a substantial portion of Earth’s biosphere and control global biogeochemical cycles; however, the rate at which they use energy (i.e., power) is virtually unknown. Here, we quantify organic matter degradation and calculate the power utilization of microbial cells throughout Earth’s Quaternary-age subseafloor sediments. Aerobic respiration, sulfate reduction, and methanogenesis mediate 6.9, 64.5, and 28.6% of global subseafloor organic matter degradation, respectively. The total power utilization of the subseafloor sediment biosphere is 37.3 gigawatts, less than 0.1% of the power produced in the marine photic zone. Aerobic heterotrophs use the largest share of global power (54.5%) with a median power utilization of 2.23 × 10 −18 watts per cell, while sulfate reducers and methanogens use 1.08 × 10 −19 and 1.50 × 10 −20 watts per cell, respectively. Most subseafloor cells subsist at energy fluxes lower than have previously been shown to support life, calling into question the power limit to life.

Microbial degradation of organic carbon in marine sediments is a key driver of global element cycles on multiple time scales. However, it is not known to what depth microorganisms alter organic carbon in marine sediments or how microbial rates of organic carbon processing change with depth, and thus time since burial, on a global scale. To better understand the connection between the dynamic carbon cycle and life’s limits in the deep subsurface, we have combined a number of global data sets with a reaction transport model (RTM) describing first, organic carbon degradation in marine sediments deposited throughout the Quaternary Period and second, a bioenergetic model for microbial activity. The RTM is applied globally, recognizing three distinct depositional environments – continental shelf, margin and abyssal zones. The results include the masses of particulate organic carbon, POC, stored in three sediment-depth layers: bioturbated Holocene (1.7 × 10^17 g C), non-bioturbated Holocene (2.6 × 10^18 g C) and Pleistocene (1.4 × 1020 g C) sediments. The global depth-integrated rates of POC degradation have been determined to be 6.8 × 10^13, 1.2 × 10^14 and 1.2 × 10^14 g C yr-1 for the same three layers, respectively. A number of maps depicting the distribution of POC, as well as the fraction that has been degraded have also been generated. Using POC degradation as a proxy for microbial catabolic activity, total heterotrophic processing of POC throughout the Quaternary is estimated to be between 10^-11 – 10^-6 g C cm-3 yr-1, depending on the time since deposition and location. Bioenergetic modeling reveals that laboratory-determined microbial maintenance powers are poor predictors of sediment biomass concentration, but that cell concentrations in marine sediments can be accurately predicted by combining bioenergetic models with the rates of POC degradation determined in this study. Our model can be used to quantitatively describe both the carbon cycle and microbial activi

A site at the gas hydrate stability limit was investigated offshore northwestern Svalbard to study methane transport in sediment. The site was characterized by chemosynthetic communities (sulfur bacteria mats, tubeworms) and gas venting. Sediments were sampled with in‐situ porewater collectors and by gravity coring followed by analyses of porewater constituents, sediment and carbonate geochemistry, and microbial activity, taxonomy, and lipid biomarkers. Sulfide and alkalinity concentrations showed concentration maxima in near‐surface sediments at the bacterial mat and deeper maxima at the gas vent site. Sediments at the periphery of the chemosynthetic field were characterized by two sulfate‐methane transition zones (SMTZ) at ~204 and 45 cm depth, where activity maxima of microbial anaerobic oxidation of methane (AOM) with sulfate were found. Amplicon sequencing and lipid biomarker indicate that AOM at the SMTZs was mediated by ANME‐1 archaea. A 1D numerical transport reaction model suggests that the deeper SMTZ‐1 formed on centennial scale by vertical advection of methane, while the shallower SMTZ‐2 could only be reproduced by non‐vertical methane injections starting on decadal scale. Model results were supported by age distribution of authigenic carbonates, showing youngest carbonates within SMTZ‐2. We propose that non‐vertical methane injection was induced by increasing blockage of vertical transport or formation of sediment fractures. Our study further suggests that the methanotrophic response to the non‐vertical methane injection was commensurate with new methane supply. This finding provides new information about for the response time and efficiency of the benthic methane filter in environments with fluctuating methane transport.

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

Highlights • Total amount of generated biogenic methane is estimated at ~3100 Gt. • Total amount of generated thermogenic methane is estimated at ~1,560 Gt. • The Maykop formation is partially productive in the central basin and not yet fully productive towards the basin peripherals. A new numerical model reconstructing the depositional history (98–0 Ma) of the Western Black Sea sub-basin is presented. The model accounts for changing boundary conditions (i.e. water depth, bottom water temperature, heat flow evolution over time) and estimates the rates and total amounts of the in-situ biogenic methane generation and thermally-driven organic matter maturation in the source rocks. The overall thermogenic and biogenic gas generation predicted by the model is estimated at ~1560 Gt and ~3100 Gt, respectively.

Marine sediments of the Blake Ridge province exhibit clearly defined geophysical indications for the presence of gas hydrates and a free gas phase. Despite being one of the world’s best-studied gas hydrate provinces and having been drilled during Ocean Drilling Program (ODP) Leg 164, discrepancies between previous model predictions and reported chemical profiles as well as hydrate concentrations result in uncertainty regarding methane sources and a possible co-existence between hydrates and free gas near the base of the gas hydrate stability zone (GHSZ). Here, by using a new multi-phase finite element (FE) numerical model, we investigate different scenarios of gas hydrate formation from both single and mixed methane sources (in-situ biogenic formation and a deep methane flux). Moreover, we explore the evolution of the GHSZ base for the past 10 Myr using reconstructed sedimentation rates and non-steady-state P-T solutions. We conclude that (1) the present-day base of the GHSZ predicted by our model is located at the depth of ~450 mbsf, thereby resolving a previously reported inconsistency between the location of the BSR at ODP Site 997 and the theoretical base of the GHSZ in the Blake Ridge region, (2) a single in-situ methane source results in a good fit between the simulated and measured geochemical profiles including the anaerobic oxidation of methane (AOM) zone, and (3) previously suggested 4 vol.%–7 vol.% gas hydrate concentrations would require a deep methane flux of ~170 mM (corresponds to the mass of methane flux of 1.6 × 10−11 kg s−1 m−2) in addition to methane generated in-situ by organic carbon (POC) degradation at the cost of deteriorating the fit between observed and modelled geochemical profiles.

Marine transform faults and associated fracture zones (MTFFZs) cover vast stretches of the ocean floor, where they play a key role in plate tectonics, accommodating the lateral movement of tectonic plates and allowing connections between ridges and trenches. Together with the continental counterparts of MTFFZs, these structures also pose a risk to human societies as they can generate high magnitude earthquakes and trigger tsunamis. Historical examples are the Sumatra-Wharton Basin Earthquake in 2012 (M8.6) and the Atlantic Gloria Fault Earthquake in 1941 (M8.4). Earthquakes at MTFFZs furthermore open and sustain pathways for fluid flow triggering reactions with the host rocks that may permanently change the rheological properties of the oceanic lithosphere. In fact, they may act as conduits mediating vertical fluid flow and leading to elemental exchanges between Earth’s mantle and overlying sediments. Chemicals transported upward in MTFFZs include energy substrates, such as H2 and volatile hydrocarbons, which then sustain chemosynthetic, microbial ecosystems at and below the seafloor. Moreover, up- or downwelling of fluids within the complex system of fractures and seismogenic faults along MTFFZs could modify earthquake cycles and/or serve as “detectors” for changes in the stress state during interseismic phases. Despite their likely global importance, the large areas where transform faults and fracture zones occur are still underexplored, as are the coupling mechanisms between seismic activity, fluid flow, and life. This manuscript provides an interdisciplinary review and synthesis of scientific progress at or related to MTFFZs and specifies approaches and strategies to deepen the understanding of processes that trigger, maintain, and control fluid flow at MTFFZs.

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.

Numerical modeling of natural gas hydrate systems requires an innovative and complex approach. The variability of parameters present in natural geological settings and the lack of wide spread high-quality 3D seismic data are the main factors limiting large-scale numerical simulations. Here, we present the outcome of a joint academic-industry project on testing the feasibility of a newly developed simulation-module included in the commercial software PetroMod TM for modeling the formation of natural gas hydrate deposits at two locations in the Gulf of Mexico. The project aimed at the scientific assessment of required input data quality and validity, choice of the computational methods, and calibration with the field data.

Mud volcanism and fluid seepage are common phenomena on the continental margin in the Gulf of Cadiz, North East Atlantic Ocean. Over the past 2 decades more than 50 mud volcanoes have been discovered and investigated interdisciplinarily. Mud volcano fluids emanating at these sites are sourced at great depths and migration is often mediated by strike slip faults in a seismically active region. The geochemical signals of the mud volcano fluids are affected by widespread various processes such as clay mineral dehydration, but also the recrystallization of ancient carbonate rocks and the alteration of oceanic crust have been suggested (Hensen et al., 2015). We developed a novel fully-coupled, basin-scale, reaction-transport model with an adaptive numerical mesh to simulate the fluid genesis in this region. An advantage of this model is the coupling of a realistic geophysical and geochemical approach, considering a growing sediment column over time together with instant compaction of sediments as well as diffusion and advection of dissolved pore water species and chemical reactions. In this proof of concept study, we looked at various scenarios to identify the processes of fluid genesis for 4 mud volcanoes, representing combinations in different subsurface settings. We can reproduce the fluid signatures (chloride, strontium, 87Sr/86Sr) of all mud volcanoes. Furthermore, we can give additional evidence that alteration of oceanic crust by fluid flow is a likely process affecting the fluid composition.

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 her

Formed under low temperature – high pressure conditions vast amounts of methane hydrates are considered to be locked up in sediments of continental margins including the Arctic shelf regions [1, 2]. Because the Arctic has warmed considerably during the recent decades and because climate models predict accelerated warming if global greenhouse gas emissions continue to rise, it is debated whether shallow Arctic hydrate deposits could be destabilized in the near future [3, 4]. Methane (CH4), a greenhouse gas with a global warming potential about 25 times higher than CO2, could be released from the melting hydrates and enter the water column and atmosphere with uncertain consequences for the environment. Here we present the results of a recent comprehensive study of the future fate of Arctic methane hydrates [5]. Our multi-disciplinary analysis provides a closer look into regional developments of submarine Arctic gas hydrate deposits under future global warming scenarios and reveals where and over which time scales gas hydrates could be destabilized and affect oceanic pH, oxygen, and atmospheric methane. Arctic bottom water temperatures and their future evolution are projected by a climate model. Predicted bottom water warming is spatially inhomogeneous, with strongest impact on shallow regions affected by Atlantic inflow. Within the next 100 years, the warming affects 25% of shallow and mid-depth regions (water depth < 600 m) containing methane hydrates. We have quantified methane release from melting hydrates using transient models resolving the change in stability zone thickness. Due to slow heat diffusion rates, the change in stability zone thickness over the next 100 years is small and methane release limited. Even if these methane emissions were to reach the atmosphere, their climatic impact would be negligible as a climate model run confirms. However, the released methane, if dissolved into the water column, may contribute to ocean acidification and oxygen deplet

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, lithologica

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 continen

Vast amounts of methane hydrates are potentially stored in sediments along the continental margins, owing their stability to low temperature – high pressure conditions. Global warming could destabilize these hydrates and cause a release of methane (CH 4) into the water column and possibly the atmosphere. Since the Arctic has and will be warmed considerably, Arctic bottom water temperatures and their future evolution projected by a climate model were analyzed. The resulting warming is spatially inhomogeneous, with the strongest impact on shallow regions affected by Atlantic inflow. Within the next 100 years, the warming affects 25% of shallow and mid-depth regions containing methane hydrates. Release of methane from melting hydrates in these areas could enhance ocean acidification and oxygen depletion in the water column. The impact of methane release on global warming, however, would not be significant within the considered time span.