Pore water concentrations and solid phase concentrations of sediment cores from the Ligurian continental slope

In October 1979, a period of heavy rainfall along the French Riviera was followed by the collapse of the Ligurian continental slope adjacent to the airport of Nice, France. A body of slope sediments, which was shortly beforehand affected by construction work south of the airport, was mobilized and traveled hundreds of kilometers downslope into the Var submarine canyon and, eventually, into the deep Ligurian basin. As a direct consequence, the construction was destroyed, seafloor cables were torn, and a small tsunami hit Antibes shortly after the failure. Hypotheses regarding the trigger mechanism include (i) vertical loading by construction of an embankment south of the airport, (ii) failure of a layer of sensitive clay within the slope sequence, and (iii) excess pore fluid pressures from charged aquifers in the underground. Over the previous decades, both the sensitive clay layers and the permeable sand and gravel layers were sampled to detect freshened waters. In 2007, the landslide scar and adjacent slopes were revisited for high-resolution seafloor mapping and systematic sampling. Results from half a dozen gravity and push cores in the shallow slope area reveal a limited zone of freshening (i.e. groundwater influence). A 100-250 m wide zone of the margin shows pore water salinities of 5-50% SW concentration and depletion in Cl, SO4, but Cr enrichment, while cores east or west of the landslide scar show regular SW profiles. Most interestingly, the three cores inside the landslide scar hint towards a complex hydrological system with at least two sources for groundwater. The aquifer system also showed strong freshening after a period of several months without significant precipitation. This freshening implies that charged coarse-grained layers represent a permanent threat to the slope's stability, not just after periods of major rainfall such as in October 1979.

Calcium isotope ratios of pore fluids and solid phase of organic rich sediments

Pore fluid calcium isotope, calcium concentration and strontium concentration data are used to measure the rates of diagenetic dissolution and precipitation of calcite in deep-sea sediments containing abundant clay and organic material. This type of study of deep-sea sediment diagenesis provides unique information about the ultra-slow chemical reactions that occur in natural marine sediments that affect global geochemical cycles and the preservation of paleo-environmental information in carbonate fossils. For this study, calcium isotope ratios (d44/40Ca) of pore fluid calcium from Ocean Drilling Program (ODP) Sites 984 (North Atlantic) and 1082 (off the coast of West Africa) were measured to augment available pore fluid measurements of calcium and strontium concentration. Both study sites have high sedimentation rates and support quantitative sulfate reduction, methanogenesis and anaerobic methane oxidation. The pattern of change of d44/40Ca of pore fluid calcium versus depth at Sites 984 and 1082 differs markedly from that of previously studied deep-sea Sites like 590B and 807, which are composed of nearly pure carbonate sediment. In the 984 and 1082 pore fluids, d44/40Ca remains elevated near seawater values deep in the sediments, rather than shifting rapidly toward the d44/40Ca of carbonate solids. This observation indicates that the rate of calcite dissolution is far lower than at previously studied carbonate-rich sites. The data are fit using a numerical model, as well as more approximate analytical models, to estimate the rates of carbonate dissolution and precipitation and the relationship of these rates to the abundance of clay and organic material. Our models give mutually consistent results and indicate that calcite dissolution rates at Sites 984 and 1082 are roughly two orders of magnitude lower than at previously studied carbonate-rich sites, and the rate correlates with the abundance of clay. Our calculated rates are conservative for these sites (the actual rates could be significantly slower) because other processes that impact the calcium isotope composition of sedimentary pore fluid have not been included. The results provide direct geochemical evidence for the anecdotal observation that the best-preserved carbonate fossils are often found in clay or organic-rich sedimentary horizons. The results also suggest that the presence of clay minerals has a strong passivating effect on the surfaces of biogenic carbonate minerals, slowing dissolution dramatically even in relation to the already-slow rates typical of carbonate-rich sediments.

Geochemistry of pore fluids from ODP Leg 134 sites

Depending on the temperature and the extent of diagenetic alteration of fluid chemistry, fluid flow at convergent margins may transfer important quantities of heat and mass between the crust and seawater, thereby influencing global mass, isotopic and heat budgets. In the North Aoba Basin, an intra-arc basin located at the New Hebrides Island Arc, alteration of volcanic ash to clay minerals and zeolites forms a CaCl2 brine, perhaps in less than 1 to 3 m.y. The brine results from an exchange of Ca for Na, K, and Mg, and an increase in Cl concentrations to a maximum of 1241 mM. The Cl increase is partly due to the transfer of H2O from the pore fluid into authigenic minerals, but water mass balances, d18O-Cl correlations, and Br/Cl ratios suggest that there is a source of Cl in the sediments. Concentration profiles indicate that Li is transferred from the fluid to solid phase at depths <300 meters below seafloor (mbsf), but at greater depths it is transferred from the solid to fluid phase, at temperatures possibly as low as 25°C. In the accretionary wedge extensive fluid flow appears to be confined to highly faulted regions. Although Cl concentrations less than seawater value are common at convergent margins, the New Hebrides margin contains little low-Cl fluid. Br/Cl ratios suggest the low-Cl fluid is from dilution, and d18O values indicate the water may be derived from mineral dehydration and mixing with meteoric water. The New Hebrides margin exhibits few surface manifestations of venting (e.g., sulfide-oxidizing benthic biological communities, carbonate crusts, mud volcanoes) and thus fluid fluxes may be smaller than at many other margins.

Uranium contents and isotope ratios of uranium and strontium in sediments, leachates and pore fluids from ODP Hole 162-984A

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the 234U/238U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%). The pore water 234U/238U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment 234U/238U activity ratios are close to 1.0. The 234U/238U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess 234U to the pore fluid by a-recoil injection of 234Th. The fraction of 238U decays that result in a-recoil injection of 234U to pore fluid is estimated to be 0.10 to 0.20 based on the 234U/238U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 0.0000004 to 0.000002 1/yr. There is significant down-hole variability in pore water 234U/238U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 1000 times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials. The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ~ 0.1 1/age for ages spanning 1000 to 500,000,000 yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.

Barium in sediments and pore waters of several ODP Leg 170 and Leg 205 sites

Barium concentrations were measured on 17 pore fluid and 13 sediment samples from Sites 1253 and 1254 drilled offshore Costa Rica during Ocean Drilling Program (ODP) Leg 205. An additional 83 pore fluid and 29 sediment samples were analyzed for Ba concentrations from Sites 1039 and 1040 drilled during ODP Leg 170 offshore Costa Rica. Sites 1039/1253 and 1040/1254 are part of a transect across the Middle America Trench offshore Nicoya Peninsula. The entire incoming sediment section is being underthrust beneath the margin, providing an ideal setting to examine Ba cycling in the shallow levels of the subduction zone. Results from these analyses indicate that a significant amount of Ba is liberated from the mineral barite (BaSO4) in the uppermost hemipelagic sediments arcward of the trench. The shallow distillation of Ba may impact the amount of sedimentary Ba reaching the deeper subduction zone.

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.

Chemical and isotope composition of poor fluid from gas-hydrate samples of ODP Site 146-892

Although the presence of extensive gas hydrate on the Cascadia margin, offshore from the western U.S. and Canada, has been inferred from marine seismic records and pore water chemistry, solid gas hydrate has only been found at one location. At Ocean Drilling Program (ODP) Site 892, offshore from central Oregon, gas hydrate was recovered close to the sediment-water interface at 2-19 m below the seafloor (mbsf) at 670 m water depth. The gas hydrate occurs as elongated platy crystals or crystal aggregates, mostly disseminated irregularly, with higher concentrations occurring in discrete zones, thin layers, and/or veinlets parallel or oblique to the bedding. A 2- to 3-cm thick massive gas hydrate layer, parallel to bedding, was recovered at ~17 mbsf. Gas from a sample of this layer was composed of both CH4 and H2S. This sample is the first mixed-gas hydrate of CH4-H2S documented in ODP; it also contains ethane and minor amounts of CO2. Measured temperatures of the recovered core ranged from 2 to -1.8°C and are 6 to 8 degrees lower than in-situ temperatures. These temperature anomalies were caused by the partial dissociation of the CH4-H2S hydrate during recovery without a pressure core sampler. During this dissociation, toxic levels of H2S (delta34S, +27.4?) were released. The delta13C values of the CH4 in the gas hydrate, -64.5 to -67.5? (PDB), together with deltaD values of -197 to -199? (SMOW) indicate a primarily microbial source for the CH4. The delta18O value of the hydrate H2O is +2.9? (SMOW), comparable with the experimental fractionation factor for sea-ice. The unusual composition (CH4-H2S) and depth distribution (2-19 mbsf) of this gas hydrate indicate mixing between a methane-rich fluid with a pore fluid enriched in sulfide; at this site the former is advecting along an inclined fault into the active sulfate reduction zone. The facts that the CH4-H2S hydrate is primarily confined to the present day active sulfate reduction zone (2-19 mbsf), and that from here down to the BSR depth (19-68 mbsf) the gas hydrate inferred to exist is a >=99% CH4 hydrate, suggest that the mixing of CH4 and H2S is a geologically young process. Because the existence of a mixed CH4-H2S hydrate is indicative of moderate to intense advection of a methane-rich fluid into a near surface active sulfate reduction zone, tectonically active (faulted) margins with organic-rich sediments and moderate to high sedimentation rates are the most likely regions of occurrence. The extension of such a mixed hydrate below the sulfate reduction zone should reflect the time-span of methane advection into the sulfate reduction zone.

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.

(Table 1) Chemical compositions of beryllium in sediments and pore waters of ODP Hole 131-808

Studies of Be distributions in subduction zone sediments will help to understand questions regarding the enrichments of cosmogenic Be-10 in arc volcanic rocks. Analyses of Be-10 and Be-9 in sediments of Ocean Drilling Program Site 808, Nankai Trough and Be-9 in porewaters of Site 808 and Sites 671 and 672, Barbados ridge complex, show significant decreases in solid phase Be-10 and large increases of porewater Be-9 at the location of the décollement zone and below or at potential flow conduits. These data imply the potential mobilization of Be during pore fluid expulsion upon sediment burial. Experiments involving reaction between a décollement sediment and a synthetic NaCl-CaCl2 solution at elevated pressure and temperatures were conducted in an attempt to mimic early subduction zone processes. The results demonstrate that Be is mobilized under elevated pressure and temperature with a strong pH dependence. The Be mobilization provides an explanation of Be-10 enrichment in arc volcanic rocks and supports the argument of the importance of the fluid processes in subduction zones at convergent margins.

Lithium isotopic composition of pore fluids and sediments in the Costa Rica subduction zone

Pore fluid and sediment Li concentrations and isotopic ratios provide important insights on the hydrology, sediment contribution to the arc volcanoes and fluid-sediment reactions at the dominantly non-accretionary Costa Rica subduction zone. Ocean Drilling Program Site 1039 in the trench axis provides a reference section of 400 m of the incoming sediments, and Site 1040, situated arcward from the trench, consists of a deformed sedimentary wedge and apron sediments, the décollement, and the partially dewatered underthrust sediment section. At the reference site, pore fluids show important isotopic variations (delta6Li=-21.7 to -37.8 per mil), reflecting the interplay of in situ alteration of volcanic material and ion exchange with clay minerals. In the basal section, a reversal of Li concentration and delta6Li toward seawater values is observed, providing supporting evidence for a lateral seawater flow system in the upper oceanic basement underlying this sediment section. At Site 1040, pore fluid of the lower deformed wedge sediments and within the décollement is enriched in Li and the isotopic compositions are relatively light, suggesting infiltration of a deep-seated fluid. The delta6Li value of -22 per mil of this Li-enriched fluid (261 µM), when compared with the delta6Li value of the subducted sediment section (-11 per mil), suggests that the deep source fluid originates from mineral fluid dehydration and transformation reactions at temperatures of 100 to 150°C, consistent with the temperature range of the up-dip seismogenic zone and of transformation of smectite to illite. The distribution of Li and its isotopes in the underthrust section are similar to those at the reference site, indicating near complete subduction of the incoming sediments and that early dewatering of the underthrust sediments occurs predominantly by lateral flow into the ocean. The hemipelagic clay-rich sediment section of the subducting plate carries most of the Li into this subduction zone, and the pelagic diatomaceous and nannofossil calcareous oozes contain little Li. The Li isotopes of both the clay-rich hemipelagic sediments and of the pelagic oozes are, however, similar, with delta6Li values of -9 to -12 per mil. The observations that (1) the delta6Li values of the underthrust sediments are distinctly lower than that of the mantle, and (2) the lavas of the Costa Rican volcanoes are enriched in Li and 7Li, provide an approximation of the contribution of the subducted sediments to the arc volcanoes. A first order mass balance calculation suggests that approximately half of the Li flux delivered by subducted sediments and altered oceanic crust into the Middle American Trench is recycled to the Costa Rican arc and at most a quarter of sedimentary Li is returned into the ocean through thrust faults, primarily the décollement thrust.

Chemical and isotopic compositions of pore fluids from ODP Sites 169-1034 and 169-1037

We have analysed the concentrations of Li, K, Rb, Cs, and B, and the isotopic ratios of Li and B of a suite of pore fluids recovered from ODP Sites 1037 (Leg 169; Escanaba Trough) and 1034 (Leg 169S; Saanich Inlet). In addition, we have analysed dissolved K, Rb, and Cs concentrations for estuarine mixing of the Ganges-Brahmaputra river system. Together, these data sets have been used to assess the role of sediments in the marine geochemical cycles of the alkali elements and boron. Uptake onto clay minerals during estuarine mixing removes 20-30% of the riverine input of dissolved Cs and Rb to the oceans. Prior to this study, the only other recognised sink of Rb and Cs was uptake during low-temperature alteration of the oceanic crust. Even with this additional sink there is an excess of inputs over outputs in their modern oceanic mass balance. Pore fluid data show that Li and Rb are transferred into marine sediments during early diagenesis. However, modeling of the Li isotope systematics of the pore fluids from Site 1037 shows that seawater Li taken up during marine sedimentation can be readily returned to solution in the presence of less hydrated cations, such as NH4+. This process also appears to result in high concentrations of pore fluid Cs (relative to local seawater) due to expulsion of adsorbed Cs from cation exchange sites. Flux calculations based on pore fluid data for a series of ODP sites indicate that early diagenesis of clay sediments removes around 8% of the modern riverine input of dissolved Li. Although NH4+-rich fluids do result in a flux of Cs to the oceans, on the global scale this input only augments the modern riverine Cs flux by ~3%. Nevertheless, this may have implications for the fate of radioactive Cs in the natural environment and waste repositories.

Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones.

The convective instability of pore-fluid flow in inclined and fluid-saturated three-dimensional fault zones has been theoretically investigated in this paper. Due to the consideration of the inclined three-dimensional fault zone with any values of the inclined angle, it is impossible to use the conventional linear stability analysis method for deriving the critical condition (i.e., the critical Rayleigh number) which can be used to investigate the convective instability of the pore-fluid flow in an inclined three-dimensional fault zone system. To overcome this mathematical difficulty, a combination of the variable separation method and the integration elimination method has been used to derive the characteristic equation, which depends on the Rayleigh number and the inclined angle of the inclined three-dimensional fault zone. Using this characteristic equation, the critical Rayleigh number of the system can be numerically found as a function of the inclined angle of the three-dimensional fault zone. For a vertically oriented three-dimensional fault zone system, the critical Rayleigh number of the system can be explicitly derived from the characteristic equation. Comparison of the resulting critical Rayleigh number of the system with that previously derived in a vertically oriented three-dimensional fault zone has demonstrated that the characteristic equation of the Rayleigh number is correct and useful for investigating the convective instability of pore-fluid flow in the inclined three-dimensional fault zone system. The related numerical results from this investigation have indicated that: (1) the convective pore-fluid flow may take place in the inclined three-dimensional fault zone; (2) if the height of the fault zone is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone stabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; (3) if the thickness of the stratum is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone destabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; and that (4) the shape of the inclined three-dimensional fault zone may affect the convective instability of pore-fluid flow in the system. (C) 2004 Published by Elsevier B.V.

(Table 1) Results of pore water analyses, including boron content and boron isotope ratios

Drilling a transect of holes across the Costa Rica forearc during ODP Leg 170 demonstrated the margin wedge to be of continental, non accretionary origin, which is intersected by permeable thrust faults. Pore waters from four drillholes, two of which penetrated the décollement zone and reached the underthrust lower plate sedimentary sequence of the Cocos Plate, were examined for boron contents and boron isotopic signatures. The combined results show dilution of the uppermost sedimentary cover of the forearc, with boron contents lower than half of the present-day seawater values. Pore fluid "refreshening" suggests that gas hydrate water has been mixed with the sediment interstitial water, without profoundly affecting the d11B values. Fault-related flux of a deeply generated fluid is inferred from high B concentration in the interval beneath the décollement, being released from the underthrust sequence with incipient burial. First-order fluid budget calculations over a cross-section across the Costa Rica forearc indicate that no significant fluid transfer from the lower to the upper plate is inferred from boron fluid profiles, at least within the frontal 40 km studied. Expulsed lower plate pore water, which is estimated to be 0.26-0.44 km3 per km trench, is conducted efficiently along and just beneath the décollement zone, indicating effective shear-enhanced compaction. In the upper plate forearc wedge, dewatering occurs as diffuse transport as well as channelled flow. A volume of approximately 2 km3 per km trench is expulsed due to compaction and, to a lesser extent, lateral shortening. Pore water chemistry is influenced by gas hydrate instability, so that it remains unknown whether deep processes like mineral dehydration or hydrocarbon formation may play a considerable role towards the hinterland.