(Table 1) Crystal sizes of measured samples from the Bush Hill region in the Gulf of Mexico and from ODP Leg 204 at the Hydrate Ridge offshore Oregon

The grain sizes of gas hydrate crystallites are largely unknown in natural samples. Single grains are hardly detectable with electron or optical microscopy. For the first time, we have used high-energy synchrotron diffraction to determine grain sizes of six natural gas hydrates retrieved from the Bush Hill region in the Gulf of Mexico and from ODP Leg 204 at the Hydrate Ridge offshore Oregon from varying depth between 1 and 101 metres below seafloor. High-energy synchrotron radiation provides high photon fluxes as well as high penetration depth and thus allows for investigation of bulk sediment samples. Gas hydrate grain sizes were measured at the Beam Line BW 5 at the HASYLAB/Hamburg. A 'moving area detector method', originally developed for material science applications, was used to obtain both spatial and orientation information about gas hydrate grains within the sample. The gas hydrate crystal sizes appeared to be (log-)normally distributed in the natural samples. All mean grain sizes lay in the range from 300 to 600 µm with a tendency for bigger grains to occur in greater depth. Laboratory-produced methane hydrate, aged for 3 weeks, showed half a log-normal curve with a mean grain size value of c. 40 µm. The grains appeared to be globular shaped.

New insights into the transport processes controlling the sulfate-methane-transition-zone near methane vents

Over the past years, several studies have raised concerns about the possible interactions between methane hydrate decomposition and external change. To carry out such an investigation, it is essential to characterize the baseline dynamics of gas hydrate systems related to natural geological and sedimentary processes. This is usually treated through the analysis of sulfate-reduction coupled to anaerobic oxidation of methane (AOM). Here, we model sulfate reduction coupled with AOM as a two-dimensional (2D) problem including, advective and diffusive transport. This is applied to a case study from a deep-water site off Nigeria’s coast where lateral methane advection through turbidite layers was suspected. We show by analyzing the acquired data in combination with computational modeling that a two-dimensional approach is able to accurately describe the recent past dynamics of such a complex natural system. Our results show that the sulfate-methane-transition-zone (SMTZ) is not a vertical barrier for dissolved sulfate and methane. We also show that such a modeling is able to assess short timescale variations in the order of decades to centuries.

用分子动力学模拟甲烷水合物热激法分解

用分子动力学模拟方法研究甲烷水合物热激法分解，系统地研究注入340 K液态水的结构Ⅰ型甲烷水合物的分解机理.模拟显示水合物表层水分子与高温液态水分子接触获得热能，分子运动激烈，摆脱水分子间的氢键束缚，笼状结构被破坏.甲烷分子获得热能从笼中挣脱，向外体系扩散.热能通过分子碰撞从外层传递给内层水分子，水合物逐层分解.对比注入277K液态水体系模拟结果，得出热激法促进水合物分解.

气体水合物生成特性实验研究

实验研究了添加剂对甲烷气体水合物形成过程的影响。发现微量的表面活性剂降低了甲烷气体水合物在静止反应器中形成的诱导时间，并使水合物快速形成和生长，提高了水合物形成过程中的填充密度。阴离子表面活性剂（十二烷基硫酸钠）对水合物生长的促进作用比非离子表面活性剂（烷基多糖苷）强。液态烃环戊烷降低了水合物形成的诱导时间，但环戊烷不能提高水合物的填充密度。

The characteristics of BSRs and their derived heat flow on the profile 973 in the northeastern South China Sea

Processing of a recently acquired seismic line in the northeastern South China Sea by Project 973 has been conducted to study the character and the distribution of gas hydrate Bottom-Simulating Reflectors (BSRs) in the Hengchun ridge. Analysis of different-type seismic profiles shows that the distribution of BSRs can be revealed to some extents by single-channel profile in this area, but seismic data processing plays an important role to resolve the full distribution of BSRs in this area. BSR' s in the northeastern South China Sea have the typical characteristics of BSRs on worldwide continental margins: they cross sediment bed reflections, they are generally parallel to the seafloor and the associated reflections have strong amplitude and a negative polarity. The characteristics of BSRs in this area are obvious and the BSRs indicate the occurrence of gas hydrate-bearing sediments in the northeastern South China Sea. The depth of the base of the gas-hydrate stability zone was calculated using the phase stability boundary curve of methane hydrate and gas hydrate with mixture gas composition and compared with the observed BSR depth. If a single gradient geothermal curve is used for the calculation, the base of the stability zone for methane hydrate or gas hydrate with a gas mixture composition does not correspond to the depth of the BSRs observed along the whole seismic profile. The geothermal gradient therefore changes significantly along the profile. The geothermal gradient and heat flow were estimated from the BSR data and the calculations show that the geothermal gradient and heat flow decrease from west to east, with the increase of the distance from the trench and the decrease of the distance to the island arc. The calculated 2 heat flow changes from 28 to 64 mW/m(2), which is basically consistent with the measured heat flow in southwestern offshore Taiwan.

笼形水簇和溶解甲烷之间平均力势能的约束型分子动力学计算及其在天然气水合物成核研究中的意义

That the dodecahedral water cluster (DWC) can adsorb dissolved methane molecules, an important phenomenon related to the hydrate nucleation study, has been observed through molecular dynamics simulations, but it has not been explained satisfactorily [Guang-Jun Guo; Yi-Gang Zhang; Hua Liu. J. Phys. Chem. C, 2007, 111, 2595]. In order to explain this phenomenon by using the potential of mean force (PMF) between the DWC and the dissolved methane, we perform several series of constrained molecular dynamics simulations in the methane-water system. The distance between the center of DWC and the methane molecule is constrained from 5 Å to 18 Å by adding 0.2 Å every time. For each fixed distance, we perform 20 independent simulations to improve the statistical precision. We first get the constraint force between the DWC and the dissolved methane in each simulation and then calculate the PMF by integrating these forces. Subsequently, the radial distribution function (RDF) is obtained from the PMF through an equation of statistical mechanics. The results show that the RDF has a sharp peak at about 6.2 Å, successfully explaining why the DWC adsorbs dissolved methane molecules. The preferential binding coefficient is a positive value (=2.05±0.5), indicates that the DWC tends to adsorb dissolved methane rather than water molecules in methane aqueous solutions. The curve of PMF for the DWC encaging a methane almost coincides that for the empty DWC, meaning that it is the DWC rather than the encaged methane who could adsorb dissolved methane molecules. By comparing the curves of PMF for different directions of the DWC relative to the dissolved methane, we find that it is the cage face rather than the cage edge or vertex that plays an essential role when the DWC adsorbing dissolved methane. This research sheds light on the driving force for the methane adsorption, and it is helpful in understanding the nucleation process of methane hydrate.

Carbon and oxygen isotopic compositions of Cibicidoides spp. and fine fraction separates (<10 µm) from ODP Holes 113-689B, 119-738B, and 120-748B

A prominent middle Eocene warming event is identified in Southern Ocean deep-sea cores, indicating that long-term cooling through the middle and late Eocene was not monotonic. At sites on Maud Rise and the Kerguelen Plateau, a distinct negative shift in d18O values (~1.0 per mil) is observed ca. 41.5 Ma. This excursion is interpreted as primarily a temperature signal, with a transient warming of 4°C over 600 k.y. affecting both surface and middle-bathyal deep waters in the Indian-Atlantic region of the Southern Ocean. This isotopic event is designated as the middle Eocene climatic optimum, and is interpreted to represent a significant climatic reversal in the midst of middle to late Eocene deep-sea cooling. The lack of a significant negative carbon isotope excursion, as observed during the Paleocene-Eocene thermal maximum, and the gradual rate of high-latitude warming suggest that this event was not triggered by methane hydrate dissociation. Rather, a transient rise in pCO2 levels is suspected, possibly as a result of metamorphic decarbonation in the Himalayan orogen or increased ridge/arc volcanism during the late middle Eocene.

Stable isotope and geochemical record of sediments from the Bering Sea

here is controversy over the role of marine methane hydrates in atmospheric methane concentrations and climate change during the last glacial period. In this study of two sediment cores from the southeast Bering Sea (700 m and 1467 m water depth), we identify multiple episodes during the last glacial period of intense methane flux reaching the seafloor. Within the uncertainty of the radiocarbon age model, the episodes are contemporaneous in the two cores and have similar timing and duration as Dansgaard-Oeschger events. The episodes are marked by horizons of sediment containing 13C-depleted authigenic carbonate minerals; 13C-depleted archaeal and bacterial lipids, which resemble those found in ANME-1 type anaerobic methane oxidizing microbial consortia; and changes in the abundance and species distribution of benthic foraminifera. The similar timing and isotopic composition of the authigenic carbonates in the two cores is consistent with a region-wide increase in the upward flux of methane bearing fluids. This study is the first observation outside Santa Barbara Basin of pervasive, repeated methane flux in glacial sediments. However, contrary to the "Clathrate Gun Hypothesis" (Kennett et al., 2003), these coring sites are too deep for methane hydrate destabilization to be the cause, implying that a much larger part of the ocean's sedimentary methane may participate in climate or carbon cycle feedback at millennial timescales. We speculate that pulses of methane in these opal-rich sediments could be caused by the sudden release of overpressure in pore fluids that builds up gradually with silica diagenesis. The release could be triggered by seismic shaking on the Aleutian subduction zone caused by hydrostatic pressure increase associated with sea level rise at the start of interstadials.

Physical properties of 39 sediment cores from the Bengal Fan and northeast Indian Ocean

We obtained sediment physical properties and geochemical data from 47 piston and gravity cores located in the Bay of Bengal, to study the complex history of the Late Pleistocene run-off from the Ganges and Brahmaputra rivers and its imprint on the Bengal Fan. Grain-size parameters were predicted from core logs of density and velocity to infer sediment transport energy and to distinguish different environments along the 3000-km-long transport path from the delta platform to the lower fan. On the shelf, 27 cores indicate rapidly prograding delta foresets today that contain primarily mud, whereas outer shelf sediment has 25% higher silt contents, indicative of stronger and more stable transport regime, which prevent deposition and expose a Late Pleistocene relic surface. Deposition is currently directed towards the shelf canyon 'Swatch of No Ground', where turbidites are released to the only channel-levee system that is active on the fan during the Holocene. Active growth of the channel-levee system occurred throughout sea-level rise and highstand with a distinct growth phase at the end of the Younger Dryas. Coarse-grained material bypasses the upper fan and upper parts of the middle fan, where particle flow is enhanced as a result of flow-restriction in well-defined channels. Sandier material is deposited mainly as sheet-flow deposits on turbidite-dominated plains at the lower fan. The currently most active part of the fan with 10-40 cm thick turbidites is documented for the central channel including inner levees (e.g., site 40). Site 47 from the lower fan far to the east of the active channel-levee system indicates the end of turbidite sedimentation at 300 ka for that location. That time corresponds to the sea-level lowering during late isotopic stage 9 when sediment supply to the fan increased and led to channel avulsion farther upstream, probably indicating a close relation of climate variability and fan activity. Pelagic deep-sea sites 22 and 28 contain a 630-kyear record of climate response to orbital forcing with dominant 21- and 41-kyear cycles for carbonate and magnetic susceptibility, respectively, pointing to teleconnections of low-latitude monsoonal forcing on the precession band to high-latitude obliquity forcing. Upper slope sites 115, 124, and 126 contain a record of the response to high-frequency climate change in the Dansgaard-Oeschger bands during the last glacial cycle with shared frequencies between 0.75 and 2.5 kyear. Correlation of highs in Bengal Fan physical properties to lows in the d18O record of the GISP2 ice-core suggests that times of greater sediment transport energy in the Bay of Bengal are associated with cooler air temperatures over Greenland. Teleconnections were probably established through moisture and other greenhouse-gas forcing that could have been initiated by instabilities in the methane hydrate reservoir in the oceans.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Geochemistry of gas hydrates of DSDP Hole 84-570

During Deep Sea Drilling Project (DSDP) Leg 84 a core 1 m long and 6 cm in diameter of massive gas hydrate was unexpectedly recovered at Site 570 in upper slope sediment of the Middle America Trench offshore of Guatemala. This core contained only 5-7% sediment, the remainder being the solid hydrate composed of gas and water. Samples of the gas hydrate were decomposed under controlled conditions in a closed container maintained at 4°C. Gas pressure increased and asymptotically approached the equilibrium decomposition pressure for an ideal methane hydrate, CH4.5-3/4H2O, of 3930 kPa and approached to this pressure after each time gas was released, until the gas hydrate was completely decomposed. The gas evolved during hydrate decomposition was 99.4% methane, ~0.2% ethane, and ~0.4% CO2. Hydrocarbons from propane to heptane were also present, but in concentrations of less than 100 p.p.m. The carbon-isotopic composition of methane was -41 to -44 per mil, relative to PDB standard. The observed volumetric methane/water ratio was 64 or 67, which indicates that before it was stored and analyzed, the gas hydrate probably had lost methane. The sample material used in the experiments was likely a mixture of methane hydrate and water ice. Formation of this massive gas hydrate probably involved the following processes: (i) upward migration of gas and its accumulation in a zone where conditions favored the growth of gas hydrates, (ii) continued, unusually rapid biological generation of methane, and (iii) release of gas from water solution as pressure decreased due to sea level lowering and tectonic uplift.

Mg/Ca and stable isotope record of benthic foraminifera of the Paleocene-Eocene Thermal Maximum

A rapid increase in greenhouse gas levels is thought to have fueled global warming at the Paleocene-Eocene Thermal Maximum (PETM). Foraminiferal magnesium/calcium ratios indicate that bottom waters warmed by 4° to 5°C, similar to tropical and subtropical surface ocean waters, implying no amplification of warming in high-latitude regions of deep-water formation under ice-free conditions. Intermediate waters warmed before the carbon isotope excursion, in association with downwelling in the North Pacific and reduced Southern Ocean convection, supporting changing circulation as the trigger for methane hydrate release. A switch to deep convection in the North Pacific at the PETM onset could have amplified and sustained warming.

Sedimentology and gas hydrates of IODP Expedition 311 samples

Expedition 311 of the Integrated Ocean Drilling Program (IODP) to northern Cascadia recovered gas-hydrate bearing sediments along a SW-NE transect from the first ridge of the accretionary margin to the eastward limit of gas-hydrate stability. In this study we contrast the gas gas-hydrate distribution from two sites drilled ~ 8 km apart in different tectonic settings. At Site U1325, drilled on a depositional basin with nearly horizontal sedimentary sequences, the gas-hydrate distribution shows a trend of increasing saturation toward the base of gas-hydrate stability, consistent with several model simulations in the literature. Site U1326 was drilled on an uplifted ridge characterized by faulting, which has likely experienced some mass wasting events. Here the gas hydrate does not show a clear depth-distribution trend, the highest gas-hydrate saturation occurs well within the gas-hydrate stability zone at the shallow depth of ~ 49 mbsf. Sediments at both sites are characterized by abundant coarse-grained (sand) layers up to 23 cm in thickness, and are interspaced within fine-grained (clay and silty clay) detrital sediments. The gas-hydrate distribution is punctuated by localized depth intervals of high gas-hydrate saturation, which preferentially occur in the coarse-grained horizons and occupy up to 60% of the pore space at Site U1325 and > 80% at Site U1326. Detailed analyses of contiguous samples of different lithologies show that when enough methane is present, about 90% of the variance in gas-hydrate saturation can be explained by the sand (> 63 µm) content of the sediments. The variability in gas-hydrate occupancy of sandy horizons at Site U1326 reflects an insufficient methane supply to the sediment section between 190 and 245 mbsf.

(Table 1) 87Sr/86Sr ratios and ages of DSDP Site 144A samples

Glassy Turonian foraminifera preserved in clay-rich sediments from the western tropical Atlantic yield the warmest equivalent d18O sea-surface temperatures (SSTs) yet reported for the entire Cretaceous-Cenozoic. We estimate Turonian SSTs that were at least as warm as (conservative mean ~30 °C) to significantly warmer (warm mean ~33 °C) than those in the region today. However, if independent evidence for high middle Cretaceous pCO2 is reliable and resulted in greater isotopic fractionation between seawater and calcite because of lower sea-surface pH, our conservative and warm SST estimates would be even higher (32 and 36°C, respectively). Our new tropical SSTs help reconcile geologic data with the predictions of general circulation models that incorporate high Cretaceous pCO2 and lend support to the hypothesis of a Cretaceous greenhouse. Our data also strengthen the case for a Turonian age for the Cretaceous thermal maximum and highlight a 20-40 m.y. mismatch between peak Cretaceous-Cenozoic global warmth and peak inferred tectonic CO2 production. We infer that this mismatch is either an artifact of a hidden Turonian pulse in global ocean-crust cycling or real evidence of the influence of some other factor on atmospheric CO2 and/or SSTs. A hidden pulse in crust cycling would explain the timing of peak Cretaceous-Cenozoic sea level (also Turonian), but other factors are needed to explain high-frequency (~10-100 k.y.) instability in middle Cretaceous SSTs reported elsewhere.

isotopic composition of gases and interstitial fluides of ODP Site 104-644

Pleistocene- to middle Miocene-age sediment was drilled at Site 341 (67? 20.1'N, 6? 06.6'E) on the inner Voring Plateau during Leg 38 of the Deep Sea Drilling Project (DSDP). In 1985, the Ocean Drilling Program (ODP) returned to the inner Wring Plateau near Site 341 and drilled a new hole at Site 644 (66° 40.7'N, 4° 34.6'E) as part of a transect to study Norwegian Sea paleoenvironments. In Hole 341, gas expansion pockets formed in cores which were recovered from depths below 50 m. This gas was characterized as predominantly methane with delta13C values in the range of -87 to -77 per mil (Morris, 1976, doi:10.2973/dsdp.proc.38.124.1976). At Site 644, sediment gas and pore-water samples were obtained to study the geochemistry of methanogenesis. Of particular interest is the possibility that methane hydrate might be present in these sediments.