Strontium concentrations and isotopic compositions of anhydrites from the TAG active mound

Sr isotope analyses have been conducted on anhydrite samples from the TAG (Trans-Atlantic Geotraverse) active hydrothermal mound (26°08?N, Mid-Atlantic Ridge) that have previously been shown to exhibit two distinct patterns of REE behavior when normalized to TAG end-member hydrothermal fluid. Despite differences in REE patterns, the Sr isotope data indicate that all the anhydrites precipitated from fluids with a similar range of hydrothermal fluid and seawater components, and all but one were seawater-dominated (52%-75%). Speciation calculations using the EQ3/6 software package for geochemical modeling of aqueous systems suggest that the REE complexation behavior in different fluid mixing scenarios can explain the variations in the REE patterns. Anhydrites that exhibit relatively flat REE patterns [(La_bs)/(Yb_bs) = 0.8-2.0; subscript bs indicates normalization to end-member black smoker hydrothermal fluid] and a small or no Eu anomaly [(Eu_bs)/(Eu*_bs) = 0.8-2.0] are inferred to have precipitated from mixes of end-member hydrothermal fluid and cold seawater. REE complexes with hard ligands (e.g., fluoride and chloride) are less stable at low temperatures and trivalent Eu has an ionic radius similar to that of Ca2+ and the other REE, and so they behave coherently. In contrast, anhydrites that exhibit slight LREE-depletion [(La_bs)/(Yb_bs) = 0.4-1.4] and a distinct negative anomaly [(Eu_bs)/(Eu*_bs) = 0.2-0.8] are inferred to have precipitated from mixes of end-member hydrothermal fluid and conductively heated seawater. The LREE depletion results from the presence of very stable LREE chloro-complexes that effectively limit the availability of the LREE for partitioning into anhydrite. Above 250°C, Eu is present only in divalent form as chloride complexes, and discrimination against Eu2+ is likely due to both the mismatch in ionic radii between Eu2+ and Ca2+, and the strong chloro-complexation of divalent Eu which promotes stability in the fluid and inhibits partitioning of Eu2+ into precipitating anhydrite. These variations in REE behavior attest to rapid fluctuations in thermal regime, fluid flow and mixing in the subsurface of the TAG mound that give rise to heterogeneity in the formation conditions of individual anhydrite crystals.

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.

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.