3-D numerical modeling of methane hydrate deposits

Within the German gas hydrate initiative SUGAR, we have developed a new tool for predicting the formation of sub-seafloor gas hydrate deposits. For this purpose, a new 2D/3D module simulating the biogenic generation of methane from organic material and the formation of gas hydrates has been added to the petroleum systems modeling software package PetroMod®. T ypically, PetroMod® simulates the thermogenic generation of multiple hydrocarbon components including oil and gas, their migration through geological strata, and finally predicts the oil and gas accumulation in suitable reservoir formations. We have extended PetroMod® to simulate gas hydrate accumulations in marine and permafrost environments by the implementation of algorithms describing (1) the physical, thermodynamic, and kinetic properties of gas hydrates; and (2) a kinetic continuum model for the microbially mediated, low temperature degradation of particulate organic carbon in sediments. Additionally, the temporal and spatial resolutions of PetroMod® were increased in order to simulate processes on time scales of hundreds of years and within decimeters of spatial extension. As a first test case for validating and improving the abilities of the new hydrate module, the petroleum systems model of the Alaska North Slope developed by IES (currently Shlumberger) and the USGS has been chosen. In this area, gas hydrates have been drilled in several wells, and a field test for hydrate production is planned for 2011/2012. The results of the simulation runs in PetroMod® predicting the thickness of the gas hydrate stability field, the generation and migration of biogenic and thermogenic methane gas, and its accumulation as gas hydrates will be shown during the conference. The predicted distribution of gas hydrates will be discussed in comparison to recent gas hydrate findings in the Alaska North Slope region

The Global Inventory of Methane Hydrate in Marine Sediments: A Theoretical Approach

The accumulation of methane hydrate in marine sediments is controlled by a number of physical and biogeochemical parameters including the thickness of the gas hydrate stability zone (GHSZ), the solubility of methane in pore fluids, the accumulation of particulate organic carbon at the seafloor, the kinetics of microbial organic matter degradation and methane generation in marine sediments, sediment compaction and the ascent of deep-seated pore fluids and methane gas into the GHSZ. Our present knowledge on these controlling factors is discussed and new estimates of global sediment and methane fluxes are provided applying a transport-reaction model at global scale. The modeling and the data evaluation yield improved and better constrained estimates of the global pore volume within the modern GHSZ ( ≥ 44 × 1015 m3), the Holocene POC accumulation rate at the seabed (~1.4 × 1014 g yr−1), the global rate of microbial methane production in the deep biosphere (4−25 × 1012 g C yr−1) and the inventory of methane hydrates in marine sediments ( ≥ 455 Gt of methane-bound carbon)

A quantitative assessment of methane cycling in Hikurangi Margin sediments (New Zealand) using geophysical imaging and biogeochemical modeling

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

An abyssal hill fractionates organic and inorganic matter in deep-sea surface sediments

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

Engineering of Cyclodextrin Product Specificity and pH Optima of the Thermostable Cyclodextrin Glycosyltransferase from Thermoanaerobacterium thermosulfurigenes EM1

The product specificity and pH optimum of the thermostable cyclodextrin glycosyltransferase (CGTase) from Thermoanaerobacterium thermosulfurigenes EM1 was engineered using a combination of x-ray crystallography and site-directed mutagenesis. Previously, a crystal soaking experiment with the Bacillus circulans strain 251 β-CGTase had revealed a maltononaose inhibitor bound to the enzyme in an extended conformation. An identical experiment with the CGTase from T. thermosulfurigenes EM1 resulted in a 2.6-Å resolution x-ray structure of a complex with a maltohexaose inhibitor, bound in a different conformation. We hypothesize that the new maltohexaose conformation is related to the enhanced α-cyclodextrin production of the CGTase. The detailed structural information subsequently allowed engineering of the cyclodextrin product specificity of the CGTase from T. thermosulfurigenes EM1 by site-directed mutagenesis. Mutation D371R was aimed at hindering the maltohexaose conformation and resulted in enhanced production of larger size cyclodextrins (β- and γ-CD). Mutation D197H was aimed at stabilization of the new maltohexaose conformation and resulted in increased production of α-CD. Glu258 is involved in catalysis in CGTases as well as α-amylases, and is the proton donor in the first step of the cyclization reaction. Amino acids close to Glu258 in the CGTase from T. thermosulfurigenes EM1 were changed. Phe284 was replaced by Lys and Asn327 by Asp. The mutants showed changes in both the high and low pH slopes of the optimum curve for cyclization and hydrolysis when compared with the wild-type enzyme. This suggests that the pH optimum curve of CGTase is determined only by residue Glu258.

New metallic allotropes of planar and tubular carbon

We propose a new family of layered sp2-like carbon crystals, incorporating five-, six-, and seven-membered rings in 2D Bravais lattices. These periodic sheets can be rolled so as to generate nanotubes of different diameter and chirality. We demonstrate that these sheets and tubes are metastable and more favorable than C60, and it is also shown that their mechanical properties are similar to those of graphene. Density of states calculations of all structures revealed an intrinsic metallic behavior, independent of orientation, tube diameter, and chirality

Simulating and Quantifying Multiple Natural Subsea CO2 Seeps at Panarea Island (Aeolian Islands, Italy) as a Proxy for Potential Leakage from Subseabed Carbon Storage Sites

Carbon dioxide (CO2) capture and storage (CCS) has been discussed as a potentially significant mitigation option for the ongoing climate warming. Natural CO2 release sites serve as natural laboratories to study subsea CO2 leakage in order to identify suitable analytical methods and numerical models to develop best-practice procedures for the monitoring of subseabed storage sites. We present a new model of bubble (plume) dynamics, advection-dispersion of dissolved CO2, and carbonate chemistry. The focus is on a medium-sized CO2 release from 294 identified small point sources around Panarea Island (South-East Tyrrhenian Sea, Aeolian Islands, Italy) in water depths of about 40–50 m. This study evaluates how multiple CO2 seep sites generate a temporally variable plume of dissolved CO2. The model also allows the overall flow rate of CO2 to be estimated based on field measurements of pH. Simulations indicate a release of ∼6900 t y–1 of CO2 for the investigated area and highlight an important role of seeps located at >20 m water depth in the carbon budget of the Panarea offshore gas release system. This new transport-reaction model provides a framework for understanding potential future leaks from CO2 storage sites

Bloodlines: mammals, leeches, and conservation in southern Asia

Southern Asia is a biodiversity hotspot both for terrestrial mammals and for leeches. Many small-mammal groups are under-studied in this region, while other mammals are of known conservation concern. In addition to standard methods for surveying mammals, it has recently been demonstrated that residual bloodmeals within leeches can be sequenced to find mammals in a given area. While these invertebrate-parasite-derived DNA (iDNA) methods are promising, most of the leech species utilized for this type of survey remain unevaluated, notwithstanding that their diversity varies substantially. Here we examine approximately 750 individual leech specimens in the genus Haemadipsa across a large range in southern Asia (Bangladesh, Cambodia, and China), specifically reviewing the diversity of mammals they feed on and their own genetic structuring. Leeches were found to feed on a considerable variety of mammals, corroborating prior studies. Additionally, leeches were found to have fed both on bats and on birds, neither of which has previously been recorded with this method. The genetic structuring of the leeches themselves revealed 15 distinct clades of which only two precisely corresponded to previously characterized species, indicating that much work is needed to finalize classifications in this genus. Most importantly, with regards to mammal conservation, leeches in these clades appear to feed on a broad range of mammals

Sensing endogenous seasonality in the case of a coffee supply chain

Rogue seasonality, or endogenously generated cyclicality (in variables), is common in supply chains and known to adversely affect performance. This paper explores a technique for sensing rogue seasonality at a supply chain echelon level. A signature and index based on cluster profiles of variables, which are meant to sense echelon-level generation and intensity of rogue seasonality, respectively, are proposed. Their validity is then established on echelons of a downstream coffee supply chain for five stock keeping units (SKUs) with contrasting rogue seasonality generation behaviour. The appropriateness of spectra as the domain for representing variables, data for which is daily sampled, is highlighted. Time-batching cycles which could corrupt the sensing are observed in variables, and the need to therefore filter them out in advance is also highlighted. The knowledge gained about the echelon location, intensity and time of generation of rogue seasonality could enable timely deployment of specific mitigation actions