Carbon and deuterium isotopes and C ratios in gas hydrates of ODP Leg 204 sites

Sediments at the southern summit of Hydrate Ridge display two distinct modes of gas hydrate occurrence. The dominant mode is associated with active venting of gas exsolved from the accretionary prism and leads to high concentrations (15%-40% of pore space) of gas hydrate in seafloor or near-surface sediments at and around the topographic summit of southern Hydrate Ridge. These near-surface gas hydrates are mainly composed of previously buried microbial methane but also contain a significant (10%-15%) component of thermogenic hydrocarbons and are overprinted with microbial methane currently being generated in shallow sediments. Focused migration pathways with high gas saturation (>65%) abutting the base of gas hydrate stability create phase equilibrium conditions that permit the flow of a gas phase through the gas hydrate stability zone. Gas seepage at the summit supports rapid growth of gas hydrates and vigorous anaerobic methane oxidation. The other mode of gas hydrate occurs in slope basins and on the saddle north of the southern summit and consists of lower average concentrations (0.5%-5%) at greater depths (30-200 meters below seafloor [mbsf]) resulting from the buildup of in situ-generated dissolved microbial methane that reaches saturation levels with respect to gas hydrate stability at 30-50 mbsf. Net rates of sulfate reduction in the slope basin and ridge saddle sites estimated from curve fitting of concentration gradients are 2-4 mmol/m**3/yr, and integrated net rates are 20-50 mmol/m**2/yr. Modeled microbial methane production rates are initially 1.5 mmol/m**3/yr in sediments just beneath the sulfate reduction zone but rapidly decrease to rates of <0.1 mmol/m**3/yr at depths >100 mbsf. Integrated net rates of methane production in sediments away from the southern summit of Hydrate Ridge are 25-80 mmol/m**2/yr. Anaerobic methane oxidation is minor or absent in cored sediments away from the summit of southern Hydrate Ridge. Ethane-enriched Structure I gas hydrate solids are buried more rapidly than ethane-depleted dissolved gas in the pore water because of advection from compaction. With subsidence beneath the gas hydrate stability zone, the ethane (mainly of low-temperature thermogenic origin) is released back to the dissolved gas-free gas phases and produces a discontinuous decrease in the C1/C2 vs. depth trend. These ethane fractionation effects may be useful to recognize and estimate levels of gas hydrate occurrence in marine sediments.

Oxygen isotopes and hydrocarbon composition of gas hydrate and hydrate water from ODP Leg 165 sites

Ocean Drilling Program (ODP) Leg 164 recovered a number of large solid gas hydrate from Sites 994, 996, and 997 on the Blake Ridge. Sites 994 and 997 samples, either nodular or thick massive pieces, were subjected to laboratory analysis and measurements to determine the structure, molecular and isotopic composition, thermal conductivity, and equilibrium dissociation conditions. X-ray computed tomography (CT) imagery, X-ray diffraction, nuclear magnetic resonance (NMR), and Raman spectroscopy have revealed that the gas hydrates recovered from the Blake Ridge are nearly 100% methane gas hydrate of Structure I, cubic with a lattice constant of a = 11.95 ± 0.05 angström, and a molar ratio of water to gas (hydration number) of 6.2. The d18O of water is 2.67 per mil to 3.51 per mil SMOW, which is 3.5-4.0 heavier than the ambient interstitial waters. The d13C and dD of methane are -66 per mil to -70 per mil and -201 per mil to -206 per mil, respectively, suggesting that the methane was generated through bacterial CO2 reduction. Thermal conductivity values of the Blake Ridge hydrates range from 0.3 to 0.5 W/(m K). Equilibrium dissociation experiments indicate that the three-phase equilibrium for the specimen is 3.27 MPa at 274.7 K. This is almost identical to that of synthetic pure methane hydrate in freshwater.

(Table 3) Run conditions and phase relations and proportions during cooling rate experiments on ODP Interval 163-989B-10R-7,55-59

Equilibrium melting and controlled cooling experiments were undertaken to constrain the crystallization and cooling histories of tholeiitic basalts recovered by the Ocean Drilling Program drilling of Site 989 on the Southeast Greenland continental margin. Isothermal experiments conducted at 1 atm. and at the fayalite-magnetite-quartz buffer using lava sample Section 163-989B-10R-7 yielded the equilibrium appearance sequence with decreasing temperature: olivine at 1184 ± 2ºC; plagioclase at 1177ºC ± 5ºC; augite at 1167 ± 5ºC; and pigeonite at 1113 ± 12ºC. In controlled cooling experiments using the same starting composition and cooling rates between 10ºC/hr and 2000ºC/hr, we find a significant temperature delay in the crystallization of olivine, plagioclase, and augite (relative to the equilibrium appearance temperature); pigeonite does not form under any dynamic crystallization conditions. Olivine exhibits the largest suppression in appearance temperature (e.g., 30º for 10ºC/hr and >190º at 100ºC/hr), while plagioclase shows the smallest (~10ºC at 10ºC/hr; 30ºC at 100ºC/hr, and ~80ºC at 1000ºC/hr). These results are in marked contrast to those obtained on lunar basalts, which generally show a large suppression of plagioclase crystallization and modest suppression of olivine crystallization with an increased cooling rate. The results we report agree well with the petrography of lavas recovered from Site 989. Furthermore, the textural analysis of run products, representing a large range of cooling rates and quench temperatures (1150ºC to 1000ºC), provide a framework for evaluating cooling conditions necessary for glass formation, rates of plagioclase growth, and kinetic factors governing plagioclase growth morphology. Specifically, we use these insights to interpret the textural and mineralogical features of the unusual compound flow recovered at Site 989. We concluded from the analysis that this flow most likely records multiple breakouts from a distal tube at an abrupt break in slope, possibly a fault scarp, resulting in the formation of a lava fan delta. This interpretation implies that normal faulting of the oldest lava sequences (lower and, possibly, middle series) preceded eruption of Site 989 lavas.

Summary of phase compositions of experimental silicate melt, un-normalised analyses

The solubility of Re and Au in haplobasaltic melt has been investigated at 1673-2573 K, 0.1 MPa-2 GPa and IW-1 to +2.5, in both carbon-saturated and carbon-free systems. Results extend the existing, low pressure and temperature, dataset to more accurately predict the results of metal-silicate equilibrium at the base of a terrestrial magma ocean. Solubilities in run-product glasses were measured by laser ablation ICP-MS, which allows for the explicit assessment of contamination by metal inclusions. The Re and Au content of demonstrably contaminant-free glasses increases with temperature, and shows variation with oxygen fugacity (fO2) similar to previous results, although lower valence states for Re (1+, 2+) are suggested by the data. At 2 GPa, and Delta IW of +1.75 to +2, the metal-silicate partition coefficient for Re (DMet/Sil) is defined by the relation LogD[met/sil][Re] = 0.50(±0.022)*10**4/T(K)+3.73(±0.095) For metal-silicate equilibrium to endow Earth's mantle with the observed time-integrated chondritic Re/Os, (and hence 187Os/188Os), DMet/Sil for both elements must converge to a common value. Combined with previously measured DMet/Sil for Os, the estimated temperature at which this convergence occurs is 4500 (±900) K. At this temperature, however, the Re and Os content of the equilibrated silicate is ~100-fold too low to explain mantle abundances. In the same experiments, much lower Dmet/sil values have been determined for Au, and require the metal-silicate equilibration temperature to be <3200 K, as hotter conditions result in an excess of Au in the mantle. Thus, the large disparity in partitioning between Re or Os, and Au at core-forming temperatures argues against their mantle concentrations set solely by metal-silicate equilibrium at the base of a terrestrial magma ocean.

Dissolution of a spatially variable DNAPL (dense non aqueous phase liquid) source

Dissolution of non-aqueous phase liquids (NAPLs) or gases into groundwater is a key process, both for contamination problems originating from organic liquid sources, and for dissolution trapping in geological storage of CO2. Dissolution in natural systems typically will involve both high and low NAPL saturations and a wide range of pore water flow velocities within the same source zone for dissolution to groundwater. To correctly predict dissolution in such complex systems and as the NAPL saturations change over time, models must be capable of predicting dissolution under a range of saturations and flow conditions. To provide data to test and validate such models, an experiment was conducted in a two-dimensional sand tank, where the dissolution of a spatially variable, 5x5 cm**2 DNAPL tetrachloroethene source was carefully measured using x-ray attenuation techniques at a resolution of 0.2x0.2 cm**2. By continuously measuring the NAPL saturations, the temporal evolution of DNAPL mass loss by dissolution to groundwater could be measured at each pixel. Next, a general dissolution and solute transport code was written and several published rate-limited (RL) dissolution models and a local equilibrium (LE) approach were tested against the experimental data. It was found that none of the models could adequately predict the observed dissolution pattern, particularly in the zones of higher NAPL saturation. Combining these models with a model for NAPL pool dissolution produced qualitatively better agreement with experimental data, but the total matching error was not significantly improved. A sensitivity study of commonly used fitting parameters further showed that several combinations of these parameters could produce equally good fits to the experimental observations. The results indicate that common empirical model formulations for RL dissolution may be inadequate in complex, variable saturation NAPL source zones, and that further model developments and testing is desirable.