The geochemistry of Ca,Sr, Ba and Ra sulfates in some deep brines from the Palo Duro Basin,Texas

The geochemistry of Ca, Sr, Ba and Ra sulfates in some deep brines from the Palo Duro Basin of north Texas, was studied to define geochemical controls on radionuclides such as ⁹⁰Sr and ²²⁶Ra. Published solubility data for gypsum, anhydrite, celestite, barite and RaSO₄ were first reevaluated, in most cases using the ion interaction approach of Pitzer, to determine solubility products of the sulfates as a function of temperature and pressure. Ionic strengths of the brines were from 2.9 to 4.8 m, their temperatures and pressures up to 40°C and 130 bars. Saturation indices of the sulfates were computed with the ion-interaction approach in one brine from the arkosic granite wash fades and four from the carbonate Wolfcamp Formation. All five brines are saturated with respect to gypsum, anhydrite and celestite, and three of the five with respect to barite. All are undersaturated by from 5 to 6 orders of magnitude with respect to pure RaSO₄. ²²⁶Ra concentrations in the brines, which ranged from 10^?11.3 to 10^?12.7 m, are not controlled by RaSO₄ solubility or adsorption, but possibly by the solubility of trace Ra solid solutions in sulfates including celestite and barite.     

Studies of equilibrium, structure, and dynamics in the aqueous Al(iii)-oxalate-fluoride system by potentiometry, C and F NMR spectroscopy

The AlOx_1-3 (Ox = oxalate) species were identified in 0.6 M aqueous NaCl by ¹³C nuclear magnetic resonance (NMR). Rate constants and activation parameters for intramolecular cis/trans isomerization of the Werner-type AlOx₂⁻ complex (k(298 K) = 5 s⁻¹, ΔH^# = 67 ± 5 kJ mol⁻¹, ΔS^# = −6 ± 6 J mol⁻¹ K⁻¹, the rate determining step could be the breaking of the Al-O(C=O) bond) and a very slow intermolecular ligand exchange reaction of AlOx₃³⁻ complex and the free ligand (k₃₀(298 K) = 6.6 · 10⁻⁵ s⁻¹, ΔH^# = 164 ± 17 kJ mol⁻¹, ΔS^# = 225 ± 51 J mol⁻¹ K⁻¹, D/I_d mechanism) were determined by dynamic 1D and 2D ¹³C NMR measurements. Mixed complexes, AlFOx, AlFOx₂^2-, AlF₂Ox⁻, and AlF₂Ox₂^3-, with overall stability (logβ) of 11.53 ± 0.03, 15.67 ± 0.03, 15.74 ± 0.02, and 19.10 ± 0.04 were measured by potentiometry using pH- and fluoride-selective electrodes and confirmed by ¹³C and¹⁹F NMR. The role of these complexes in gibbsite dissolution was modeled. The mixed Al(III)-Ox^2--F⁻ complexes have to be considered as the chemical speciation of Al(III) in natural waters is discussed.     

Thermochemistry of yavapaiite KFe(SO4)2: Formation and decomposition

Yavapaiite, KFe(SO₄)₂, is a rare mineral in nature, but its structure is considered as a reference for many synthetic compounds in the alum supergroup. Several authors mention the formation of yavapaiite by heating potassium jarosite above ca. 400°C. To understand the thermal decomposition of jarosite, thermodynamic data for phases in the K-Fe-S-O-(H) system, including yavapaiite, are needed. A synthetic sample of yavapaiite was characterized in this work by X-ray diffraction (XRD), scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR), and thermal analysis. Based on X-ray diffraction pattern refinement, the unit cell dimensions for this sample were found to be a = 8.152 ± 0.001 Å, b = 5.151 ± 0.001 Å, c = 7.875 ± 0.001 Å, and β = 94.80°. Thermal decomposition indicates that the final breakdown of the yavapaiite structure takes place at 700°C (first major endothermic peak), but the decomposition starts earlier, around 500°C. The enthalpy of formation from the elements of yavapaiite, KFe(SO₄)₂, ΔH°_f = −2042.8 ± 6.2 kJ/mol, was determined by high-temperature oxide melt solution calorimetry. Using literature data for hematite, corundum, and Fe/Al sulfates, the standard entropy and Gibbs free energy of formation of yavapaiite at 25°C (298 K) were calculated as S°(yavapaiite) = 224.7 ± 2.0 J.mol⁻¹.K⁻¹ and ΔG°_f = −1818.8 ± 6.4 kJ/mol. The equilibrium decomposition curve for the reaction jarosite = yavapaiite + Fe₂O₃ + H₂O has been calculated, at pH₂O = 1 atm, the phase boundary lies at 219 ± 2°C.     

Tracer diffusion of Mg, Ca, Sr, and Ba in Na-aluminosilicate melts

We employed the thin source technique to investigate tracer diffusion of Mg, Ca, Sr, and Ba in glasses and supercooled melts of albite (NaAlSi₃O₈) and jadeite (NaAlSi₂O₆) compositions. The experiments were conducted at 1 bar and at temperatures between 645 and 1025°C. Typical run durations ranged between 30 min and 35 days. The analysis of the diffusion profiles was performed with the electron microprobe. Diffusivities of Ca, Sr, and Ba were found to be independent of either duration t of the experiment or tracer concentration M, initially introduced into the sample. Mg exhibits a diffusivity depending on run time and concentration and tracer diffusivity is derived by extrapolation to M/√t = 0. Temperature dependence of the diffusivity D can be represented by an Arrhenius equation D = D_o exp(−E_a/RT), yielding the following least-squares fit parameters (with D in m²/s and E_a in kJ/mol): D_Mg = 1.8 · 10⁻⁵ exp(−234 ± 20/RT), D_Ca = 3.5 · 10⁻⁶ exp(−159 ± 6/RT), D_Sr = 3.6 · 10⁻⁶ exp(−160 ± 6/RT), and D_Ba = 6.0 · 10⁻⁶ exp(−188 ± 12/RT) for albite; and D_Mg = 8.3 · 10⁻⁶ exp(−207 ± 18/RT), D_Ca = 3.8 · 10⁻⁶ exp(−153 ± 4/RT), D_Sr = 2.3 · 10⁻⁶ exp(−150 ± 4/RT), and D_Ba = 3.7 · 10⁻⁵ exp(−198 ± 4/RT) for jadeite composition. Ca and Sr diffusivities agree within error in both compositions and exhibit the fastest diffusivities, whereas Mg reveals the lowest diffusivity. The relationship between activation energy and radius shows a minimum at Ca and Sr for albite and jadeite compositions extending the relationship already observed elsewhere for alkalies. With increasing substitution of Si by (Na + Al), diffusivities increase, whereas activation energies decrease. Furthermore, a simple model modified from that of Anderson and Stuart (Anderson O. L. and Stuart D. A., “Calculation of activation energy of ionic conductivity in silica glasses by classical methods,” J. Am. Ceram. Soc.37, 573-580, 1954) is discussed for calculating the activation energies.     

Cadmium deposition and mobility in the sediments of an acidic oligotrophic lake

A dated core from the profoundal zone in a pristine oligotrophic acidic lake was analyzed for Cd as well as for Al, Ca, Fe, Mg, Mn, Pb, Ti and total carbon and nitrogen. Overlying water and porewater samples were also obtained on six occasions at the same site, and yielded vertical profiles of pH and dissolved Cd, Ca, Fe, Mg, Mn, sulfide, SO₄⁻², organic and inorganic carbon concentrations. These extensive porewater and sediment geochemical data were used, together with information on infaunal benthos, to decipher the sedimentary record of Cd contamination. Depth variation of sediment Ca concentrations indicate that the lake suffered from progressive acidification starting about 1950. The present-day accumulation rate of Cd (J_acc^Cd = 5.4 ± 0.4 × 10⁻¹¹ mol cm⁻² yr⁻¹) in the sediments is the sum of the flux of Cd deposited with settling particles (J_S^Cd = 3.3 ± 0.2 × 10⁻¹¹ mol cm⁻² yr⁻¹) and the fluxes of dissolved Cd across the sediment-water interface due to molecular diffusion (J_D^Cd = 1.8 ± 0.3 × 10⁻¹¹ mol cm⁻² yr⁻¹), bioturbation (J_B^Cd = 1.1 ± 0.2 × 10⁻¹⁴ mol cm⁻² yr⁻¹) and bioirrigation (J_I^Cd = 0.27 ± 0.05 × 10⁻¹¹ mol cm⁻² yr⁻¹). Biological mixing of the sediments was negligible. The shape of the vertical profile of total Cd concentration with depth in the sediment appears to be determined more by its input history than by post-depositional mobilization and redistribution in the sediment column.     

Radium uptake during barite recrystallization at 23 ± 2 °C as a function of solution composition: An experimental Ba and Ra tracer study

High-purity synthetic barite powder was added to pure water or aqueous solutions of soluble salts (BaCl₂, Na₂SO₄, NaCl and NaHCO₃) at 23 ± 2 °C and atmospheric pressure. After a short pre-equilibration time (4 h) the suspensions were spiked either with ¹³³Ba or ²²⁶Ra and reacted under constant agitation during 120-406 days. The pH values ranged from 4 to 8 and solid to liquid (S/L) ratios varied from 0.01 to 5 g/l. The uptake of the radiotracers by barite was monitored through repeated sampling of the aqueous solutions and radiometric analysis. For both ¹³³Ba and ²²⁶Ra, our data consistently showed a continuous, slow decrease of radioactivity in the aqueous phase.Mass balance calculations indicated that the removal of ¹³³Ba activity from aqueous solution cannot be explained by surface adsorption only, as it largely exceeded the 100% monolayer coverage limit. This result was a strong argument in favor of recrystallization (driven by a dissolution-precipitation mechanism) as the main uptake mechanism. Because complete isotopic equilibration between aqueous solution and barite was approached or even reached in some experiments, we concluded that during the reaction all or substantial fractions of the initial solid had been replaced by newly formed barite.The ¹³³Ba data could be successfully fitted assuming constant recrystallization rates and homogeneous distribution of the tracer into the newly formed barite. An alternative model based on partial equilibrium of ¹³³Ba with the mineral surface (without internal isotopic equilibration of the solid) could not reproduce the measured activity data, unless multistage recrystallization kinetics was assumed. Calculated recrystallization rates in the salt solutions ranged from 2.8 × 10⁻¹¹ to 1.9 × 10⁻¹⁰ mol m⁻² s⁻¹ (2.4-16 μmol m⁻² d⁻¹), with no specific trend related to solution composition. For the suspensions prepared in pure water, significantly higher rates (∼5.7 × 10⁻¹⁰ mol m⁻² s⁻¹ or ∼49 μmol m⁻² d⁻¹) were determined.Radium uptake by barite was determined by monitoring the decrease of ²²⁶Ra activity in the aqueous solution with alpha spectrometry, after filtration of the suspensions and sintering. The evaluation of the Ra uptake experiments, in conjunction with the recrystallization data, consistently indicated formation of non-ideal solid solutions, with moderately high Margules parameters (W_AB = 3720-6200 J/mol, a₀ = 1.5-2.5). These parameters are significantly larger than an estimated value from the literature (W_AB = 1240 J/mol, a₀ = 0.5).In conclusion, our results confirm that radium forms solid solutions with barite at fast kinetic rates and in complete thermodynamic equilibrium with the aqueous solutions. Moreover, this study provides quantitative thermodynamic data that can be used for the calculation of radium concentration limits in environmentally relevant systems, such as radioactive waste repositories and uranium mill tailings.     

Joint determination of K decay constants and Ar∗/K for the Fish Canyon sanidine standard, and improved accuracy for Ar/Ar geochronology

⁴⁰Ar/³⁹Ar and K-Ar geochronology have long suffered from large systematic errors arising from imprecise K and Ar isotopic data for standards and imprecisely determined decay constants for the branched decay of ⁴⁰K by electron capture and β⁻ emission. This study presents a statistical optimization approach allowing constraints from ⁴⁰K activity data, K-Ar isotopic data, and pairs of ²³⁸U-²⁰⁶Pb and ⁴⁰Ar/³⁹Ar data for rigorously selected rocks to be used as inputs for estimating the partial decay constants (λ_ε and λ_β) of ⁴⁰K and the ⁴⁰Ar∗/⁴⁰K ratio (κ_FCs) of the widely used Fish Canyon sanidine (FCs) standard. This yields values of κ_FCs = (1.6418 ± 0.0045) × 10⁻³, λ_ε = (0.5755 ± 0.0016) × 10⁻¹⁰ a⁻¹ and λ_β = (4.9737 ± 0.0093) × 10⁻¹⁰ a⁻¹. These results improve uncertainties in the decay constants by a factor of >4 relative to values derived from activity data alone. Uncertainties in these variables determined by our approach are moderately to highly correlated (cov(κ_FCs, λ_ε) = 7.1889 × 10⁻¹⁹, cov(κ_FCs, λ_β) = −7.1390 × 10⁻¹⁹, cov(λ_ε, λ_β) = −3.4497 × 10⁻²⁶) and one must take account of the covariances in error propagation by either linear or Monte Carlo methods. ⁴⁰Ar/³⁹Ar age errors estimated from these results are significantly reduced relative to previous calibrations. Also, age errors are smaller for a comparable level of isotopic measurement precision than those produced by the ²³⁸U/²⁰⁶Pb system, because the ⁴⁰Ar/³⁹Ar system is now jointly calibrated by both the ⁴⁰K and ²³⁸U decay constants, and because λ_ε(⁴⁰K) < λ(²³⁸U). Based on this new calibration, the age of the widely used Fish Canyon sanidine standard is 28.305 ± 0.036 Ma. The increased accuracy of ⁴⁰Ar/³⁹Ar ages is now adequate to provide meaningful validation of high-precision U/Pb or astronomical tuning ages in cases where closed system behavior of K and Ar can be established.     

Diffusion kinetics of Fe and Mg in aluminous spinel: experimental determination and applications

The diffusion coefficients of Fe²⁺ and Mg in aluminous spinel at ∼20 kb, 950 to 1325°C, and at 30 kb, 1125°C have been determined via diffusion couple experiments and numerical modeling of the induced diffusion profiles. The oxygen fugacity, fO₂, was constrained by graphite encapsulating materials. The retrieved self-diffusion coefficients of Fe²⁺ and Mg at ∼20 kb, 950 to 1325°C, fit well the Arrhenian relation, D = D₀exp(−Q/RT), where Q is the activation energy, with D₀(Fe) = 1.8 (±2.8) × 10⁻⁵, D₀(Mg) = 1.9 (±1.4) × 10⁻⁵ cm²/s, Q(Fe) = 198 ± 19, and Q(Mg) = 202 ± 8 kJ/mol. Comparison with the data at 30 kb suggests an activation volume of ∼5 cm³/mol. From analysis of compositional zoning in natural olivine-spinel assemblages in ultramafic rocks, previous reports concluded that D(Fe-Mg) in spinel with Cr/(Cr + Al) ≤0.5 is ∼10 times that in olivine. The diffusion data in spinel and olivine have been applied to the problems of preservation of Mg isotopic inhomogeneity in spinel within the plagioclase-olivine inclusions in Allende meteorite and cooling rates of terrestrial ultramafic rocks.     

Characterization of calcium isotopes in natural and synthetic barite

The mineral barite (BaSO₄) accommodates calcium in its crystal lattice, providing an archive of Ca-isotopes in the highly stable sulfate mineral. Holocene marine (pelagic) barite samples from the major ocean basins are isotopically indistinguishable from each other (δ^44/40Ca = −2.01 ± 0.15‰) but are different from hydrothermal and cold seep barite samples (δ^44/40Ca = −4.13 to −2.72‰). Laboratory precipitated (synthetic) barite samples are more depleted in the heavy Ca-isotopes than pelagic marine barite and span a range of Ca-isotope compositions, Δ^44/40Ca = −3.42 to −2.40‰. Temperature, saturation state, , and aCa²⁺/aBa²⁺ each influence the fractionation of Ca-isotopes in synthetic barite; however, the fractionation in marine barite samples is not strongly related to any measured environmental parameter. First-principles lattice dynamical modeling predicts that at equilibrium Ca-substituted barite will have much lower ⁴⁴Ca/⁴⁰Ca than calcite, by −9‰ at 0 °C and −8‰ at 25 °C. Based on this model, none of the measured barite samples appear to be in isotopic equilibrium with their parent solutions, although as predicted they do record lower δ^44/40Ca values than seawater and calcite. Kinetic fractionation processes therefore most likely control the extent of isotopic fractionation exhibited in barite. Potential fractionation mechanisms include factors influencing Ca²⁺ substitution for Ba²⁺ in barite (e.g. ionic strength and trace element concentration of the solution, competing complexation reactions, precipitation or growth rate, temperature, pressure, and saturation state) as well as nucleation and crystal growth rates. These factors should be considered when investigating controls on isotopic fractionation of Ca²⁺ and other elements in inorganic and biogenic minerals.     

The effect of dissolved barium on biogeochemical processes at cold seeps

A numerical model was applied to investigate and quantify the biogeochemical processes fueled by the expulsion of barium and methane-rich fluids in the sediments of a giant cold-seep area in the Derugin Basin (Sea of Okhotsk). Geochemical profiles of dissolved Ba²⁺, Sr²⁺, Ca²⁺, SO₄²⁻, HS⁻, DIC, I⁻ and of calcium carbonate (CaCO₃) were fitted numerically to constrain the transport processes and the kinetics of biogeochemical reactions. The model results indicate that the anaerobic oxidation of methane (AOM) is the major process proceeding at a depth-integrated rate of 4.9 μmol cm⁻² a⁻¹, followed by calcium carbonate and strontian barite precipitation/dissolution processes having a total depth-integrated rate of 2.1 μmol cm⁻² a⁻¹. At the low seepage rate prevailing at our study site (0.14 cm a⁻¹) all of the rising barium is consumed by precipitation of barite in the sedimentary column and no benthic barium flux is produced. Numerical experiments were run to investigate the response of this diagenetic environment to variations of hydrological and biogeochemical conditions. Our results show that relatively low rates of fluid flow (<∼5 cm a⁻¹) promote the dispersed precipitation of up to 26 wt% of barite and calcium carbonate throughout the uppermost few meters of the sedimentary column. Distinct and persistent events (several hundreds of years long) of more vigorous fluid flow (from 20-110 cm a⁻¹), instead, result in the formation of barite-carbonate crusts near the sediment surface. Competition between barium and methane for sulfate controls the mineralogy of these sediment precipitates such that at low dissolved methane/barium ratios (<4-11) barite precipitation dominates, while at higher methane/barium ratios sulfate availability is limited by AOM and calcium carbonate prevails. When seepage rates exceed 110 cm a⁻¹, barite precipitation occurs at the seafloor and is so rapid that barite chimneys form in the water column. In the Derugin Basin, spectacular barite constructions up to 20 m high, which cover an area of roughly 22 km² and contain in excess of 5 million tons of barite, are built through this process. In these conditions, our model calculates a flux of barium to the water column of at least 20 μmol cm⁻² a⁻¹. We estimate that a minimum of 0.44 × 10⁶ mol a⁻¹ are added to the bottom waters of the Derugin Basin by cold seep processes, likely affecting the barium cycle in the Sea of Okhotsk.     

Mn-Cr systematics of carbonates in CM chondrites: Implications for the timing and duration of aqueous alteration

CM chondrites contain carbonates and other secondary minerals such as phyllosilicates, sulfides, sulfates, oxides and hydroxides that are believed to have formed by aqueous alteration reactions on their parent asteroid. We report in situ Mn-Cr isotope measurements in the highly aqueously altered CM2.1 chondrites QUE 93005 and ALH 83100 using secondary ion mass spectrometry (Cameca ims-1270 ion microprobe). The ⁵³Cr excesses are correlated with the ⁵³Mn/⁵⁵Mn ratio and result from the in situ decay of ⁵³Mn, a short-lived radioisotope with a half-life of 3.7 Ma. If we assume that carbonate grains in samples QUE 93005 and ALH 83100 are cogenetic, then the excesses define initial ⁵³Mn/⁵⁵Mn ratios ((⁵³Mn/⁵⁵Mn)₀) of (4.1 ± 1.2) × 10⁻⁶ and (5.1 ± 1.7) × 10⁻⁶, respectively. These values are comparable to those in carbonates from other CM chondrites as reported in the literature. Initial ⁵³Mn/⁵⁵Mn ratios for calculated model isochrones for individual carbonate grains range from (3.8 ± 1.4) × 10⁻⁶ to (4.8 ± 2.1) × 10⁻⁶ for QUE 93005 and from (3.1 ± 1.6) × 10⁻⁶ to (1.3 ± 0.5) × 10⁻⁵ for ALH 83100. A possible interpretation for the ranges in (⁵³Mn/⁵⁵Mn)₀ could be that alteration in individual CM chondrites was episodic and occurred over an extended period of time. However, isochrones based on the entire set of carbonate grains in each of the CM chondrites imply that the degree of aqueous alteration is roughly correlated with the age of carbonate formation in CM chondrites of different subtypes and that alteration on the CM parent asteroid started contemporaneously with or shortly after CAI formation and lasted at least 4 Ma.     

Thermodynamics of iron oxides: Part III. Enthalpies of formation and stability of ferrihydrite (∼Fe(OH)3), schwertmannite (∼FeO(OH)3/4(SO4)1/8), and ε-Fe2O3

Enthalpies of formation of ferrihydrite and schwertmannite were measured by acid solution calorimetry in 5 N HCl at 298 K. The published thermodynamic data for these two phases and ε-Fe₂O₃ were evaluated, and the best thermodynamic data for the studied compounds were selected.Ferrihydrite is metastable in enthalpy with respect to α-Fe₂O₃ and liquid water by 11.5 to 14.7 kJ•mol⁻¹ at 298.15 K. The less positive enthalpy corresponds to 6-line ferrihydrite, and the higher one, indicating lesser stability, to 2-line ferrihydrite. In other words, ferrihydrite samples become more stable with increasing crystallinity. The best thermodynamic data set for ferrihydrite of composition Fe(OH)₃ was selected by using the measured enthalpies and (1) requiring ferrihydrite to be metastable with respect to fine-grained lepidocrocite; (2) requiring ferrihydrite to have entropy higher than the entropy of hypothetical, well-crystalline Fe(OH)₃; and (3) considering published estimates of solubility products of ferrihydrite. The ΔG°_f for 2-line ferrihydrite is best described by a range of −708.5±2.0 to −705.2±2.0 kJ•mol⁻¹, and ΔG°_f for 6-line ferrihydrite by −711.0±2.0 to −708.5±2.0 kJ•mol⁻¹.A published enthalpy measurement by acid calorimetry of ε-Fe₂O₃ was re-evaluated, arriving at ΔH°_f (ε-Fe₂O₃) = −798.0±6.6 kJ•mol⁻¹. The standard entropy (S°) of ε-Fe₂O₃ was considered to be equal to S° (γ-Fe₂O₃) (93.0±0.2 J•K⁻¹•mol⁻¹), giving ΔG°_f (ε-Fe₂O₃) = −717.8±6.6 kJ•mol⁻¹. ε-Fe₂O₃ thus appears to have no stability field, and it is metastable with respect to most phases in the Fe₂O₃-H₂O system which is probably the reason why this phase is rare in nature.Enthalpies of formation of two schwertmannite samples are: ΔH°_f (FeO(OH)_0.686(SO₄)_0.157•0.972H₂O) = −884.0±1.3 kJ•mol⁻¹, ΔH°_f (FeO(OH)_0.664(SO₄)_0.168•1.226H₂O) = −960.7±1.2 kJ•mol⁻¹. When combined with an entropy estimate, these data give Gibbs free energies of formation of −761.3 ± 1.3 and −823.3 ± 1.2 kJ•mol⁻¹ for the two samples, respectively. These ΔG_f° values imply that schwertmannite is thermodynamically favored over ferrihydrite over a wide range of pH (2-8) when the system contains even small concentration of sulfate. The stability relations of the two investigated samples can be replicated by schwertmannite of the “ideal” composition FeO(OH)_3/4(SO₄)_1/8 with ΔG°_f = −518.0±2.0 kJ•mol⁻¹.     

Sulfur isotope systematics of the 1982 El Chichón trachyandesite: an ion microprobe study

Sulfur isotopic compositions were determined by ion microprobe for 36 spots on anhydrite crystals in trachyandesitic pumices erupted from El Chichón Volcano in 1982. Individual anhydrite crystals are homogeneous in δ³⁴S, within the ±1‰ (2σ) uncertainty of the method, but crystal-to-crystal variations are large (+2.5 to +10.9‰). The mean δ³⁴S for anhydrite (+6.4 ± 2.1‰, 1σ) is significantly lower than earlier results for bulk anhydrite separates (+9.0 to +9.2‰). The difference between the mean δ³⁴S values in these two populations may reflect a grain-size effect, with heavier sulfur concentrated in smaller anhydrite crystals, few of which were analyzed by ion microprobe. Variations in δ³⁴S show no correlation with complex textures in anhydrite revealed by cathodoluminescence color. Ion-microprobe analyses of δ³⁴S were also obtained on six ovoid-shaped inclusions of pyrrhotite, chalcopyrite, and/or intermediate sulfide solid solution hosted by silicate or oxide crystals, interpreted to be magmatic (δ³⁴S = −0.1 to +2.7‰; mean +0.7‰), and on four irregularly shaped multiphase sulfide fragments in the matrix, interpreted as xenocrystic, which range widely in δ³⁴S (−3.7 to +5.5‰). We evaluate four different mixing scenarios involving (1) magmatic anhydrite and sedimentary sulfate, (2) magmatic anhydrite and hydrothermal anhydrite, and anhydrite and coexisting sulfide crystals precipitated in different domains of a common magma reservoir that were affected by (3) different degrees of degassing or (4) different degrees of crustal sulfur contamination. The model involving physical contamination of sedimentary sulfate is considered untenable. The other three models are considered to be viable, but none of them can explain all observations. The results of this study and other recent investigations prompt a re-evaluation of the sulfur budget for the 1982 El Chichón eruption. We estimate that 2.2 × 10¹³ g of S was emitted, and that 58 wt.% of the sulfur was present as anhydrite prior to eruption, with the remainder in a vapor phase, with H₂S/SO₂ ≈ 9. The bulk magmatic δ³⁴S value for the 1982 El Chichón trachyandesite is estimated as +4.1 to +5.8‰, typical of the relatively heavy sulfur isotopic compositions that characterize subduction-related magmas.     

Modeling diagenesis of lead in sediments of a Canadian Shield lake

Triplicate porewater lead concentration profiles were determined on six occasions in a Canadian Shield lake. Total Pb concentrations were also measured in a dated core obtained at the same site. This information, as well as an extensive dataset comprising ancillary geochemical measurements on porewaters and sediment and the population densities of benthic animals, is used in a one-dimensional transport-reaction diagenetic model to investigate the transport and mobilization of Pb in these sediments. Application of the model consistently indicates the presence of a zone of Pb production to the porewaters that lies above a zone of Pb consumption. The profiles of various porewater constituents and thermodynamic calculations indicate that Pb is mobilized in the zone of production by the reductive dissolution of iron oxyhydroxides, whereas it is removed in the zone of consumption by precipitation as a solid sulfide. Rate constants are estimated for reductive iron dissolution (k_d^Fe(III) = 2.0 ± 0.5 × 10⁻¹ cm³ mol⁻¹ s⁻¹), Pb adsorption on iron oxyhydroxides (k_ads^Pb = 98 ± 55 cm³ mol⁻¹ s⁻¹), and Pb precipitation (k_ppt^Pb = 8 × 10⁻²⁰ mol cm⁻³ s⁻¹ to 16 ± 13 × 10⁻²² mol cm⁻³ s⁻¹, depending on the solubility product assumed for the precipitation of PbS). According to model calculations, diagenetic processes, such as remobilization, molecular diffusion, bioturbation, and bioirrigation have a negligible influence on the solid phase Pb profile. In agreement with this finding, the present-day fluxes of dissolved Pb by diffusion (J_D^Pb = −6.5 × 10⁻¹¹ mol cm⁻² yr⁻¹), bioturbation (J_B^Pb = −1.1 × 10⁻¹³ mol cm⁻² yr⁻¹), and bioirrigation (J_I^Pb = −1.5 × 10⁻¹¹ mol cm⁻² yr⁻¹) are small compared to the flux of Pb deposited with settling particles (J_S^Pb = 5.3 × 10⁻⁹ mol cm⁻² yr⁻¹).     

Flow of river water into a Karstic limestone aquifer. 1. Tracing the young fraction in groundwater mixtures in the Upper Floridan Aquifer near Valdosta, Georgia

The quality of water in the Upper Floridan aquifer near Valdosta, Georgia is affected locally by discharge of Withlacoochee River water through sinkholes in the river bed. Data on transient tracers and other dissolved substances, including Cl⁻, ³H, tritiogenic helium-3 (³He), chlorofluorocarbons (CFC-11, CFC-12, CFC-113), organic C (DOC), O₂ (DO), H₂S, CH₄, δ¹⁸O, δD, and ¹⁴C were investigated as tracers of Withlacoochee River water in the Upper Floridan aquifer. The concentrations of all tracers were affected by dilution and mixing. Dissolved Cl⁻, δ¹⁸O, δD, CFC-12, and the quantity (³H+³He) are stable in water from the Upper Floridan aquifer, whereas DOC, DO, H₂S, CH₄, ¹⁴C, CFC-11, and CFC-113 are affected by microbial degradation and other geochemical processes occurring within the aquifer. Groundwater mixing fractions were determined by using dissolved Cl⁻ and δ¹⁸O data, recognizing 3 end-member water types in the groundwater mixtures: (1) Withlacoochee River water (δ¹⁸O=−2.5±0.3‰, Cl⁻=12.2±2 mg/l), (2) regional infiltration water (δ¹⁸O=−4.2±0.1‰, Cl⁻=2.3±0.1 mg/l), and (3) regional paleowater resident in the Upper Floridan aquifer (δ¹⁸O=−3.4±0.1‰, Cl⁻=2.6±0.1 mg/l) (uncertainties are ±1σ). Error simulation procedures were used to define uncertainties in mixing fractions. Fractions of river water in groundwater range from 0 to 72% and average 10%. The influence of river-water discharge on the quality of water in the Upper Floridan aquifer was traced from the sinkhole area on the Withlacoochee River 25 km SE in the direction of regional groundwater flow. Infiltration of water is most significant to the N and NW of Valdosta, but becomes negligible to the S and SE in the direction of general thickening of post-Eocene confining beds overlying the Upper Floridan aquifer.     

Luminescence centers in anhydrite,barite, celestite and their synthesized analogs

The intrinsic luminescence center (LC) of SO ₄ ^2? appears at 360 nm in all types of sulfates, but is absent in sulfates with large quantities of impurities. Three nonequivalent Gd³⁺ LC have been observed in luminescence spectra of anhydrite. Gd³⁺, Ce³⁺ O ₂ ^? LC were established in celestite. Photoluminescence (PL) bands in the sulfates were assigned as follows: the band at 460 nm is related to Eu²⁺ (a connection with Al³⁺ or donor-acceptor pair is possible also) in sulfates of Ca and Ba; bands at 520 nm and 590 nm are related to VO₄, MoO₄ and TiO₄ substituting SO₄; the band at 660 nm is assigned to Ti³⁺ X-ray luminescence (XL) band at 620 nm is assigned to Ag⁺.     

Synthesis, characterization, and thermochemistry of K-Na-H3O jarosites

The thermochemistry of well-characterized synthetic K-H₃O, Na-H₃O and K-Na-H₃O jarosites was investigated. These phases are solid solutions that obey Vegard’s law. Electron probe microanalyses indicated lower alkali and iron contents than predicted from the theoretical end-member compositions, in agreement with thermal analyses, suggesting the presence of hydronium and “additional” water. The standard enthalpies of formation (ΔH°_f) of K-H₃O, Na-H₃O and K-Na-H₃O jarosites were determined by high-temperature oxide melt solution calorimetry. These enthalpies vary linearly with the K/H₃O, Na/H₃O and K/Na ratio, respectively. The enthalpy of formation of pure hydronium jarosite was also determined experimentally (ΔH°_f = −3741.6 ± 8.3 kJ.mol⁻¹), and it was used to evaluate ΔH°_f for the end-members KFe₃(SO₄)₂(OH)₆ (ΔH°_f = −3829.6 ± 8.3 kJ.mol⁻¹) and NaFe₃(SO₄)₂(OH)₆ (ΔH°_f = −3783.4 ± 8.3 kJ.mol⁻¹). Finally, enthalpies of dehydration (loss of the “additional” water) of some jarosites were determined and found to be near the enthalpy of vaporization of water, suggesting that the “additional” water is weakly bonded in the structure.     

The kinetics of oxygen exchange between the GeO4Al12(OH)24(OH2)12(aq) molecule and aqueous solutions

We report rates of oxygen exchange with bulk solution for an aqueous complex, ^IVGeO₄Al₁₂(OH)₂₄(OH₂)₁₂⁸⁺(aq) (GeAl₁₂), that is similar in structure to both the ^IVAlO₄Al₁₂(OH)₂₄(OH₂)₁₂⁷⁺(aq) (Al₁₃) and ^IVGaO₄Al₁₂(OH)₂₄(OH₂)₁₂⁷⁺(aq) (GaAl₁₂) molecules studied previously. All of these molecules have ε-Keggin-like structures, but in the GeAl₁₂ molecule, occupancy of the central tetrahedral metal site by Ge(IV) results in a molecular charge of +8, rather than +7, as in the Al₁₃ and GaAl₁₂. Rates of exchange between oxygen sites in this molecule and bulk solution were measured over a temperature range of 274.5 to 289.5 K and 2.95 < pH < 4.58 using ¹⁷O-NMR.Apparent rate parameters for exchange of the bound water molecules (η-OH₂) are k_ex²⁹⁸ = 200 (±100) s⁻¹, ΔH^‡ = 46 (±8) kJ · mol⁻¹, and ΔS^‡ = −46 (±24) J · mol⁻¹ K⁻¹ and are similar to those we measured previously for the GaAl₁₂ and Al₁₃ complexes. In contrast to the Al₁₃ and GaAl₁₂ molecules, we observe a small but significant pH dependence on rates of solvolysis that is not yet fully constrained and that indicates a contribution from the partly deprotonated GeAl₁₂ species.The two topologically distinct μ₂-OH sites in the GeAl₁₂ molecule exchange at greatly differing rates. The more labile set of μ₂-OH sites in the GeAl₁₂ molecule exchange at a rate that is faster than can be measured by the ¹⁷O-NMR isotopic-equilibration technique. The second set of μ₂-OH sites have rate parameters of k_ex²⁹⁸ = 6.6 (±0.2) · 10⁻⁴ s⁻¹, ΔH^‡ = 82 (±2) kJ · mol⁻¹, and ΔS^‡ = −29 (±7) J · mol⁻¹ · K⁻¹, corresponding to exchanges ≈40 and ≈1550 times, respectively, more rapid than the less labile μ₂-OH sites in the Al₁₃ and GaAl₁₂ molecules. We find evidence of nearly first-order pH dependence on the rate of exchange of this μ₂-OH site with bulk solution for the GeAl₁₂ molecule, which contrasts with Al₁₃ and GaAl₁₂ molecules.     

Dissolution and carbonation of a serpentinite: Inferences from acid attack and high P-T experiments performed in aqueous solutions at variable salinity

Dissolution experiments on a serpentinite were performed at 70 °C, 0.1 MPa, in H₂SO₄ solution, in open and closed systems, in order to evaluate the overall dissolution rate of mineral components over different times (4, 9 and 24 h). In addition, the serpentinite powder was reacted with a NaCl-bearing aqueous solution and supercritical CO₂ for 24 h at higher pressures (9-30 MPa) and temperatures (250-300 °C) either in a stirred reactor or in an externally-heated pressure vessel to assess both the dissolution rate of serpentinite minerals and the progress of the carbonation reaction. Results show that, at 0.1 MPa, MgO extraction from serpentinite ranges from 82% to 98% and dissolution rate varies from 8.5 × 10⁻¹⁰ mole m⁻² s⁻¹ to 4.2 × 10⁻⁹ mole m⁻² s⁻¹. Attempts to obtain carbonates from the Mg-rich solutions by increasing their pH failed since Mg- and NH₄- bearing sulfates promptly precipitated. On the other hand, at higher pressures, significant crystallization (5.0-10.4 wt%) of Ca- and Fe-bearing magnesite was accomplished at 30 MPa and 300 °C using 100 g L⁻¹ NaCl aqueous solutions. The corresponding amount of CO₂ sequestered by crystallization of carbonates is 9.4-15.9 mole%. Dissolution rate (from 6.3 × 10⁻¹¹ mole m⁻² s⁻¹ to 1.3 × 10⁻¹⁰ mole m⁻² s⁻¹) is lower than that obtained at 0.1 MPa and 70 °C but it is related to pH values much higher (3.3-4.4) than that (−0.65) calculated for the H₂SO₄ solution.Through a thorough review of previous experimental investigations on the dissolution kinetics of serpentine minerals the authors propose adopting: (i) the log rate [mole m⁻² s⁻¹] value of −12.08 ± 0.16 (1σ), as representative of the neutral dissolution mechanism at 25 °C and (ii) the following relationship for the acidic dissolution mechanism at 25 °C:

log rate=-0.45(±0.09)×pH-10.01(±0.30).     

Effects of seawater carbonate ion concentration and temperature on shell U, Mg, and Sr in cultured planktonic foraminifera

We investigate the sensitivity of U/Ca, Mg/Ca, and Sr/Ca to changes in seawater [CO₃²⁻] and temperature in calcite produced by the two planktonic foraminifera species, Orbulina universa and Globigerina bulloides, in laboratory culture experiments. Our results demonstrate that at constant temperature, U/Ca in O. universa decreases by 25 ± 7% per 100 μmol [CO₃²⁻] kg⁻¹, as seawater [CO₃²⁻] increases from 110 to 470 μmol kg⁻¹. Results from G. bulloides suggest a similar relationship, but U/Ca is consistently offset by ∼+40% at the same environmental [CO₃²⁻]. In O. universa, U/Ca is insensitive to temperature between 15°C and 25°C. Applying the O. universa relationship to three U/Ca records from a related species, Globigerinoides sacculifer, we estimate that Caribbean and tropical Atlantic [CO₃²⁻] was 110 ± 70 μmol kg⁻¹ and 80 ± 40 μmol kg⁻¹ higher, respectively, during the last glacial period relative to the Holocene. This result is consistent with estimates of the glacial-interglacial change in surface water [CO₃²⁻] based on both modeling and on boron isotope pH estimates. In settings where the addition of U by diagenetic processes is not a factor, down-core records of foraminiferal U/Ca have potential to provide information about changes in the ocean’s carbonate concentration.Below ambient pH (pH < 8.2), Mg/Ca decreased by 7 ± 5% (O. universa) to 16 ± 6% (G. bulloides) per 0.1 unit increase in pH. Above ambient pH, the change in Mg/Ca was not significant for either species. This result suggests that Mg/Ca-based paleotemperature estimates for the Quaternary, during which surface-ocean pH has been at or above modern levels, have not been biased by variations in surface-water pH. Sr/Ca increased linearly by 1.6 ± 0.4% per 0.1 unit increase in pH. Shell Mg/Ca increased exponentially with temperature in O. universa, where Mg/Ca = 0.85 exp (0.096*T), whereas the change in Sr/Ca with temperature was within the reproducibility of replicate measurements.