Diffusion transport of water in basalts

Studying diffusive transport in porous rocks is of fundamental importance in understanding a variety of geochemical processes including: element transfer, primary mineral dissolution kinetics and precipitation of secondary phases. Here we report new findings on the relationship between diffusive transport and textural characteristics of the pore systems on the example of mid-oceanic ridge basalts having different degree of alteration but very similar bulk pore volume. Diffusion processes in porous basalts were studied in situ using H2O -> D2O exchange experiments. The effective diffusion coefficients of water molecules increase systematically from 5.05*10**-11 to 1.19*10**-10 m**2/s for fresh and moderately altered basalts and from 2.40*10**-11 to 6.72*10**-11 m**2/s for completely altered basalt as temperature increases from 5 to 50 °C. The activation energy of the diffusion process increases from 12.29 ± 0.71 kJ/mol for fresh and moderately altered basalts to 14.3 ± 1.33 kJ/mol for completely altered basalt. The results indicate that neither the bulk porosity nor the degree of alteration can be used as proxies for the efficiency of element transport during MORB-water interaction. The formation of secondary phases that replace primary minerals and fill the pore space in the rock leads to the formation of tiny pores and phases with large specific surface area. These factors might have a dominant control on the transport properties of altered basaltic rocks.

Petrophysical measurement of ODP Leg 122 and 123 samples (Tables 1, 2)

The capillary-pressure characteristics of 22 samples of lithified post-Paleozoic Indian-Ocean carbonates were compared to published data from older carbonate rocks (lower Paleozoic Hunton Group of Texas and Oklahoma). The Indian-Ocean samples are considerably more porous than are the Paleozoic samples, yet all of the Indian-Ocean samples fit readily into a descriptive petrofacies scheme previously established for the Hunton Group. The Indian-Ocean samples may be assigned to four petrophysical facies (petrofacies) based on the shapes of their capillary-pressure curves, their pore-throat-size distributions, their estimated recovery efficiency values (for nonwetting fluids), and the visual characteristics of their pore systems, as observed with a scanning-electron microscope. Petrofacies assignments for the Indian-Ocean samples are as follows. Petrofacies I includes six samples collected from the coarse basal portions of event deposits (primarily turbidites). These samples have large throats, leptokurtic throat-size distributions, low- to moderate recovery efficiency values, concave cumulative-intrusion capillary-pressure curves, and high porosity values. Petrofacies II includes two sedimentologically dissimilar samples that have medium-size throats, platykurtic throat-size distributions, moderate- to-high recovery efficiency values, gently sloping cumulative-intrusion capillary-pressure curves, and high porosity values. Petrofacies III includes two polymictic sandstones and a skeletal packstone that have small throats, polymodal throat-size distributions, moderate recovery efficiency values, gently sloping cumulative-intrusion capillary-pressure curves, and high porosity values. Petrofacies IV includes 11 samples, mostly recrystallized neritic carbonates, that have small throats, leptokurtic throat-size distributions, high recovery efficiency values, convex cumulative-intrusion capillary-pressure curves, and low porosity values. Comparison of petrofacies assignment to core-, thin-section-, and smear-slide data, and to inferred depositional setting, suggests that pore systems in most samples from Holes 765C and 766A result from primary depositional features, whereas pore systems in samples from Hole 761C and one sample from Hole 765C have been strongly influenced by diagenetic processes. For Hole 761C, prediction of petrophysical parameters should be most successful if based on diagenetic facies patterns. By contrast, the distribution of favorable reservoir facies and of permeability barriers in less highly altered rocks collected from Holes 765C and 766A is related to depositional patterns. Recovery efficiency is inversely related to both porosity and median throat size for the present data set. This relationship is similar to that observed for carbonates of the lower Paleozoic Hunton Group and the Ordovician Ellenburger dolomite, but opposite of that observed for some other ancient carbonates. The coarse deposits of the massive basal units of turbidites are petrophysically distinct and form a coherent petrophysical group (Petrofacies I) with substantial reservoir potential. Two samples assigned to Petrofacies I have extremely large throats (median throat size at least 4 ?m, and at least six times that of any other sample) and therefore high permeability values. These two samples come from thin, coarse turbidites that lack or have poorly developed fine divisions and are interpreted to have been deposited on channeled suprafan lobes in a proximal mid-fan setting. The restriction of extremely high permeability values to a single depositional facies suggests that careful facies mapping of deep-sea fans in a deliberate search for such coarse turbidites could dramatically enhance the success of exploration for aquifers or hydrocarbon reservoirs. Such reservoirs should have substantial vertical heterogeneity. They should have high lateral permeability values but low vertical permeability values, and reservoir sections should include numerous thin units having widely differing petrophysical characteristics.

Sediment and pore water geochemistry of sediment cores, Potter Cove, King George Island

Redox-sensitive trace metals (Mn, Fe, U, Mo, Re), nutrients and terminal metabolic products (NO3-, NH4+, PO43-, total alkalinity) were for the first time investigated in pore waters of Antarctic coastal sediments. The results of this study reveal a high spatial variability in redox conditions in surface sediments from Potter Cove, King George Island, western Antarctic Peninsula. Particularly in the shallower areas of the bay the significant correlation between sulphate depletion and total alkalinity, the inorganic product of terminal metabolism, indicates sulphate reduction to be the major pathway of organic matter mineralisation. In contrast, dissimilatory metal oxide reduction seems to be prevailing in the newly ice-free areas and the deeper troughs, where concentrations of dissolved iron of up to 700 µM were found. We suggest that the increased accumulation of fine-grained material with high amounts of reducible metal oxides in combination with the reduced availability of metabolisable organic matter and enhanced physical and biological disturbance by bottom water currents, ice scouring and burrowing organisms favours metal oxide reduction over sulphate reduction in these areas. Based on modelled iron fluxes we calculate the contribution of the Antarctic shelf to the pool of potentially bioavailable iron (Feb) to be 6.9x10**3 to 790x10**3 t/yr. Consequently, these shelf sediments would provide an Feb flux of 0.35-39.5/mg/m**2/yr (median: 3.8 mg/m**2/yr) to the Southern Ocean. This contribution is in the same order of magnitude as the flux provided by icebergs and significantly higher than the input by aeolian dust. For this reason suboxic shelf sediments form a key source of iron for the high nutrient-low chlorophyll (HNLC) areas of the Southern Ocean. This source may become even more important in the future due to rising temperatures at the WAP accompanied by enhanced glacier retreat and the accumulation of melt water derived iron-rich material on the shelf.

Stable chlorine isotopes in pore waters from DSDP Hole 79-546 and ODP Hole 131-808C

Stable Cl isotope ratios, measured in marine pore waters associated with the Barbados and Nankai subduction zones, extend significantly (to ~-8 per mil) the range of d37Cl values reported for natural waters. These relatively large negative values, together with geologic and chemical evidence from Barbados and Nankai and recent laboratory data showing that hydrous silicate minerals (i.e., those with structural OH sites) are enriched up to 7.5 per mil in 37Cl relative to seawater, strongly suggest that the isotopic composition of Cl in pore waters from subduction zones reflects diagenetic and metamorphic dehydration and transformation reactions. These reactions involve clays and/or other hydrous silicate phases at depth in the fluid source regions. Chlorine therefore cannot be considered geochemically conservative in these systems. The uptake of Cl by hydrous phases provides a mechanism by which Cl can be cycled into the mantle through subduction zones. Thus, stable Cl isotopes should help in determining the extent to which Cl and companion excess volatiles like H2O and CO2 cycle between the crust and mantle.

Magnesium isotope ratios of pore-fluids and dolomites

Magnesium concentrations in deep-sea sediment pore-fluids typically decrease down core due to net precipitation of dolomite or clay minerals in the sediments or underlying crust. To better characterize and differentiate these processes, we have measured magnesium isotopes in pore-fluids and sediment samples from Ocean Drilling Program sites (1082, 1086, 1012, 984, 1219, and 925) that span a range of oceanographic settings. At all sites, magnesium concentrations decrease with depth. At sites where diagenetic reactions are dominated by the respiration of organic carbon, pore-fluid d26Mg values increase with depth by as much as 2 per mil. Because carbonates preferentially incorporate 24Mg (low d26Mg), the increase in pore-fluid d26Mg values at these sites is consistent with the removal of magnesium in Mg-carbonate (dolomite). In contrast, at sites where the respiration of organic carbon is not important and/or weatherable minerals are abundant, pore-fluid d26Mg values decrease with depth by up to 2 per mil. The decline in pore-fluid d26Mg at these sites is consistent with a magnesium sink that is isotopically enriched relative to the pore-fluid. The identity of this enriched magnesium sink is likely clay minerals. Using a simple 1D diffusion-advection-reaction model of pore-fluid magnesium, we estimate rates of net magnesium uptake/removal and associated net magnesium isotope fractionation factors for sources and sinks at all sites. Independent estimates of magnesium isotope fractionation during dolomite precipitation from measured d26Mg values of dolomite samples from sites 1082 and 1012 are very similar to modeled net fractionation factors at these sites, suggesting that local exchange of magnesium between sediment and pore-fluid at these sites can be neglected. Our results indicate that the magnesium incorporated in dolomite is 2.0-2.7 per mil depleted in d26Mg relative to the precipitating fluid. Assuming local exchange of magnesium is minor at the rest of the studied sites, our results suggest that magnesium incorporated into clay minerals is enriched in d26Mg by 0 per mil to +1.25 per mil relative to the precipitating fluid. This work demonstrates the utility of magnesium isotopes as a tracer for magnesium sources/sinks in low-temperature aqueous systems.

Sulphate concentrations and sulphur isotopic composition in pore fluids from ODP Leg 168 sites

A numerical model of sulfate reduction and isotopic fractionation has been applied to pore fluid SO4**2- and d34S data from four sites drilled during Ocean Drilling Program (ODP) Leg 168 in the Cascadia Basin at 48°N, where basement temperatures reach up to 62°C. There is a source of sulfate both at the top and the bottom of the sediment column due to the presence of basement fluid flow, which promotes bacterial sulfate reduction below the sulfate minimum zone at elevated temperatures. Pore fluid d34S data show the highest values (135 per mil) yet found in the marine environment. The bacterial sulfur isotopic fractionation factor, a, is severely underestimated if the pore fluids of anoxic marine sediments are assumed to be closed systems and Rayleigh fractionation plots yield erroneous values for a by as much as 15 per mil in diffusive and advective pore fluid regimes. Model results are consistent with a = 1.077+/-0.007 with no temperature effect over the range 1.8 to 62°C and no effect of sulfate reduction rate over the range 2 to 10 pmol/ccm/day. The reason for this large isotopic fractionation is unknown, but one difference with previous studies is the very low sulfate reduction rates recorded, about two orders of magnitude lower than literature values that are in the range of µmol/ccm/day to tens of nmol/ccm/day. In general, the greatest 34S depletions are associated with the lowest sulfate reduction rates and vice versa, and it is possible that such extreme fractionation is a characteristic of open systems with low sulfate reduction rates.