Molecular modeling of water structure in nano-pores between brucite (001) surfaces

Molecular dynamics (MD) computer simulations of liquid water held in one-dimensional nano-confinement by two parallel, electrostatically neutral but hydrophilic surfaces of brucite, Mg(OH)₂, provide greatly increased, atomistically detailed understanding of surface-related effects on the spatial variation in the structural ordering, hydrogen bond (H-bond) organization, and local density of H₂O molecules at this important model hydroxide surface. NVT-ensemble MD simulations (i.e., at constant number of atoms, volume and temperature) were performed for a series of model systems consisting of 3 to 30 Å-thick water layers (containing 35 to 360 H₂O molecules) confined between two 19 Å-thick brucite substrate layers. The results show that the hydrophilic substrate significantly influences the near-surface water structure, with both H-bond donation to the surface oxygen atoms and H-bond acceptance from the surface hydrogen atoms in the first surface layer of H₂O molecules playing key roles. Profiles of oxygen and hydrogen atomic density and H₂O dipole orientation show significant deviation from the corresponding structural properties of bulk water to distances as large as 15 Å (∼5 molecular water layers) from the surface, with the local structural environment varying significantly with the distance from the surface. The water molecules in the first layer at about 2.45 Å from the surface have a two-dimensional hexagonal arrangement parallel to brucite layers, reflecting the brucite surface structure, have total nearest neighbor coordinations of 5 or 6, and are significantly limited in their position and orientation. The greatest degree of the tetrahedral (ice-like) ordering occurs at about 4 Å from the surface. The translational and orientational ordering of H₂O molecules in layers further from the surface become progressively more similar to those of bulk liquid water. A quantitative statistical analysis of the MD-generated instantaneous molecular configurations in terms of local density, molecular orientation, nearest neighbor coordination, and the structural details of the H-bonding network shows that the local structure of interfacial water at the brucite surface results from a combination of “hard wall” (geometric and confinement) effects, highly directional H-bonding, and thermal motion. This structure does not resemble that of bulk water at ambient conditions or at elevated or reduced temperature, but shares some similarities with that of water under higher pressure.     

Pyrite oxidation in moist air

The rate of pyrite oxidation in moist air was determined by measuring, over time, the pressure difference between a sealed chamber containing pyrite plus oxygen and a control. The experiments carried out at 25°C, 96.7% fixed relative humidity, and oxygen partial pressures of 0.21, 0.61, and 1.00 atm showed that the rate of oxygen consumption is a function of oxygen partial pressure and time. The rates of oxygen consumption (r, mol/m²sec) fit the expression

(A)     

Oxygen-exchange pathways in aluminum polyoxocations

Using molecular dynamics simulations and electronic structure methods, we postulate a mechanism to explain the complicated reactivity trends that are observed for oxygen isotope exchange reactions between sites in aluminum polyoxocations of the ε-Keggin type and bulk solution. Experimentally, the molecules have four nonequivalent oxygens that differ considerably in reactivity both within a molecule, and between molecules in the series: Al₁₃, GaAl₁₂, and GeAl₁₂ [MO₄Al₁₂(OH)₂₄(H₂O)₁₂ⁿ⁺(aq); with M = Al(III) for Al₁₃, n = 7; M = Ga(III) for GaAl₁₂, n = 7; M = Ge(IV) for GeAl₁₂, n = 8]. We find that a partly dissociated, metastable intermediate molecule of expanded volume is necessary for exchange of both sets of μ₂-OH and that the steady-state concentration of this intermediate reflects the bond strengths between the central metal and the μ₄-O. Thus the central metal exerts extraordinary control over reactions at hydroxyl bridges, although these are three bonds away.This mechanism not only explains the reactivity trends for oxygen isotope exchange in μ₂-OH and η-OH₂ sites in the ε-Keggin aluminum molecules, but also explains the observation that the reactivities of minerals tend to reflect the presence of highly coordinated oxygens, such as the μ₄-O in boehmite, α-, and γ-Al₂O₃ and their Fe(III) analogs. The partial dissociation of these highly coordinated oxygens, coupled with simultaneous activation and displacement of neighboring metal centers, may be a fundamental process by which metals atoms undergo ligand exchanges at mineral surfaces.     

Fe isotopic fractionation during mineral dissolution with and without bacteria

Fe released into solution is isotopically lighter (enriched in the lighter isotope) than hornblende starting material when dissolution occurs in the presence of the siderophore desferrioxamine mesylate (DFAM). In contrast, Fe released from goethite dissolving in the presence of DFAM is isotopically unchanged. Furthermore, Δ⁵⁶Fe_{solution-hornblende} for Fe released to solution in the presence of ligands varies with the affinity of the ligand for Fe. The extent of isotopic fractionation of Fe released from hornblende also increases when experiments are agitated continuously. The Fe isotope fractionation observed during hornblende dissolution with organic ligands is attributed predominantly to retention of ⁵⁶Fe in an altered surface layer, while the lack of isotopic fractionation during goethite dissolution in DFAM is consistent with the lack of an altered layer. When a siderophore-producing soil bacterium is added to the system (without added organic ligands), Fe released to solution from both hornblende and goethite differs isotopically from Fe in the bulk mineral: Δ⁵⁶Fe_{solution-starting material} = −0.56 ± 0.19 (hornblende) and −1.44 ± 0.16 (goethite). Increased isotopic fractionation is attributed in this case to the fact that as bacterial respiration depletes the system in oxygen and aqueous Fe is reduced, equilibration between aqueous ferrous and ferric iron creates a pool of isotopically heavy ferric iron that is assimilated by bacterial cells. Adsorption of isotopically heavy ferrous iron (Fe(II) enriched in the heavier isotope) or precipitation of isotopically heavy Fe minerals may also contribute to observed fractionations.To test whether these Fe isotope signatures are recorded in natural systems, we also investigated extractions of samples of soils from which the bacteria were isolated. These extractions show variability in the isotopic signatures of exchangeable Fe and Fe oxyhydroxide fractions from one soil sample to another, but exchangeable Fe is observed to be lighter than Fe in soil Fe oxyhydroxides and hornblende. This observation is consistent with isotopically light Fe-organic complexes in soil pore water derived from the Fe-silicate starting materials in the presence of growing microorganisms, as documented in experiments reported here. The contributions from phenomena including organic ligand-promoted nonstoichiometric dissolution of Fe silicates, uptake of ferric iron by organisms, adsorption of isotopically heavy ferrous iron, and precipitation of iron minerals should create complex isotopic signatures in soils. Better understanding of these processes and the timescales over which they contribute to fractionation is needed.     

Carbon isotope effects in carbonate systems

Global carbon cycle models require a complete understanding of the δ¹³C variability of the Earth’s C reservoirs as well as the C isotope effects in the transfer of the element among them. An assessment of δ¹³C changes during CO₂ loss from degassing magmas requires knowledge of the melt-CO₂ carbon isotope fractionation. In order to examine the potential size of this effect for silicate melts of varying composition, ¹³C reduced partition functions were computed in the temperature range 275 to 4000 K for carbonates of varying bond strengths (Mg, Fe, Mn, Sr, Ba, Pb, Zn, Cd, Li, and Na) and the polymorphs of calcite. For a given cation and a given pressure the ¹³C content increases with the density of the carbonate structure. For a given structure the tendency to concentrate ¹³C increases with pressure. The effect of pressure (‰/10 kbar) on the size of the reduced partition function of aragonite varies with temperature; in the pressure range 1 to 10⁵ bars the change is given by:

(1)     

Manganese-chromium isotope systematics of enstatite meteorites

We present results of a study of the ⁵³Mn-⁵³Cr isotope systematics in the enstatite chondrites and achondrites (aubrites). The goal of this study was to explore the capabilities of this isotope system to obtain chronological information on these important classes of meteorites and to investigate the original distribution in the inner solar system of the short-lived radionuclide ⁵³Mn. Our earlier work (Lugmair and Shukolyukov, 1998; Shukolyukov and Lugmair, 2000a) has shown that the asteroid belt bodies are characterized by essentially the same initial ⁵³Mn abundance. However, we have found the presence of a gradient in the abundance of the radiogenic ⁵³Cr between the earth-moon system, Mars, and the asteroid Vesta. If this gradient is considered as a function of the heliocentric distance a linear radial dependence is indicated. This can be explained either by an early, volatility controlled Mn/Cr fractionation in the nebula or by an original radially heterogeneous distribution of ⁵³Mn. The enstatite chondrites are suggested to form in the inner zones of the solar nebula, much closer to the Sun than the ordinary chondrites. Therefore, their investigation may be an important test on the hypothesis on a radial heterogeneity in the initial ⁵³Mn.We have studied the bulk samples of the EH4-chondrites Indarch and Abee and the EL6-chondrite Khairpur. Although these meteorites have essentially the same Mn/Cr ratio as the ordinary chondrites, the relative abundance of the radiogenic ⁵³Cr is three times smaller than in the ordinary chondrites. Because these meteorites are primitive (undifferentiated) and no Mn/Cr fractionation had occurred within their parent bodies, this difference is a strong argument in favor of an initially heterogeneous distribution of ⁵³Mn in the early inner solar system. This finding is also consistent with formation of the enstatite chondrites in the inner zones of the solar nebula. Using the characteristic ⁵³Cr excess of the enstatite chondrites and the observed gradient, their place of origin falls at about 1.4 AU or somewhat closer to the Sun (i.e. >1.0-1.4 AU).We also present chronological results for the enstatite chondrites and achondrites. The ‘absolute’ ⁵³Mn-⁵³Cr ages of the EH4-chondrites are old: ∼4565 Ma. The EL6-chondrite Khairpur is ∼4.5 Ma younger, which is in good agreement with the ¹²⁹I-¹²⁹Xe data from the literature. The age of the aubrite Peña Blanca Spring appears to be similar to those of the enstatite chondrites while that of the aubrite Bishopville is at least ∼10 Ma younger, which is also in agreement with the ¹²⁹I-¹²⁹Xe data. The results from bulk samples of aubrites indicate that the last Mn/Cr fractionation in their parent body occurred ∼ 4563 Ma ago and imply an evolution of the Mn-Cr isotope system in an environment with an higher than chondritic Mn/Cr ratio for several millions of years.     

Role of molecular oxygen in the dissolution of siderite and rhodochrosite

The dissolution of siderite (FeCO₃) and rhodochrosite (MnCO₃) under oxic and anoxic conditions is investigated at 298 K. The anoxic dissolution rate of siderite is 10^−8.65 mol m⁻² s⁻¹ for 5.5 < pH < 12 and increases as [H⁺]^0.75 for pH < 5.5. The pH dependence is consistent with parallel proton-promoted and water hydrolysis dissolution pathways. Atomic force microscopy (AFM) reveals a change in pit morphology from rhombohedral pits for pH > 4 to pits elongated at one vertex for pH < 4. Under oxic conditions the dissolution rate decreases to below the detection limit of 10⁻¹⁰ mol m⁻² s⁻¹ for 6.0 < pH < 10.3, and hillock precipitation preferential to steps is observed in concurrent AFM micrographs. X-ray photoelectron spectroscopy (XPS) and thermodynamic analysis identify the precipitate as ferrihydrite. At pH > 10.3, the oxic dissolution rate is as high as 10^−7.5 mol m⁻² s⁻¹, which is greater than under the corresponding anoxic conditions. A fast electron transfer reaction between solution O₂ or [Fe³⁺(OH)₄]⁻ species and surficial >Fe^II hydroxyl groups is hypothesized to explain the dissolution kinetics. AFM micrographs do not show precipitation under these conditions. Anoxic dissolution of rhodochrosite is physically observed as rhombohedral pit expansion for 3.7 < pH < 10.3 and is chemically explained by parallel proton- and water-promoted pathways. The dissolution rate law is 10^−4.93[H⁺] + 10^−8.45 mol m⁻² s⁻¹. For 5.8 < pH < 7.7 under oxic conditions, the AFM micrographs show a tabular precipitate growing by preferential expansion along the a-axis, though the macroscopic dissolution rate is apparently unaffected. For pH > 7.7 under oxic conditions, the dissolution rate decreases from 10^−8.45 to 10^−9.0 mol m⁻² s⁻¹. Flattened hillock precipitates grow across the entire surface without apparent morphological influence by the underlying rhodochrosite surface. XPS spectra and thermodynamic calculations implicate the precipitate as bixbyite for 5.8 < pH < 7.7 and MnOOH (possibly feitnkechtite) for pH >7.7.     

A stable isotope-based approach to tropical dendroclimatology

We describe a strategy for development of chronological control in tropical trees lacking demonstrably annual ring formation, using high resolution δ¹⁸O measurements in tropical wood. The approach applies existing models of the oxygen isotopic composition of alpha-cellulose (Roden et al., 2000), a rapid method for cellulose extraction from raw wood (Brendel et al., 2000), and continuous flow isotope ratio mass spectrometry (Brenna et al., 1998) to develop proxy chronological, rainfall and growth rate estimates from tropical trees lacking visible annual ring structure. Consistent with model predictions, pilot datasets from the temperate US and Costa Rica having independent chronological control suggest that observed cyclic isotopic signatures of several permil (SMOW) represent the annual cycle of local rainfall and relative humidity. Additional data from a plantation tree of known age from ENSO-sensitive northwestern coastal Peru suggests that the 1997-8 ENSO warm phase event was recorded as an 8‰ anomaly in the δ¹⁸O of α-cellulose. The results demonstrate reproducibility of the stable isotopic chronometer over decades, two different climatic zones, and three tropical tree genera, and point to future applications in paleoclimatology.     

Thermochemistry of jarosite-alunite and natrojarosite-natroalunite solid solutions

The thermochemistry of jarosite-alunite and natrojarosite-natroalunite solid solutions was investigated. Members of these series were either coprecipitated or synthesized hydrothermally and were characterized by XRD, FTIR, electron microprobe analysis, ICP-MS, and thermal analysis. Partial alkali substitution and vacancies on the Fe/Al sites were observed in all cases, and the solids studied can be described by the general formula K_1-x-yNa_y(H₃O)_xFe_zAl_w(SO₄)₂(OH)_6-3(3-z-w)(H₂O)_3(3-z-w). A strong preferential incorporation of Fe over Al in the jarosite/alunite structure was observed. Heats of formation from the elements, ΔH°_f, were determined by high-temperature oxide melt solution calorimetry. The solid solutions deviate slightly from thermodynamic ideality by exhibiting positive enthalpies of mixing in the range 0 to +11 kJ/mol. The heats of formation of the end members of both solid solutions were derived. The values ΔH°_f = −3773.6 ± 9.4 kJ/mol, ΔH°_f = −4912.2 ± 24.2 kJ/mol, ΔH°_f = −3734.6 ± 9.7 kJ/mol and ΔH°_f = −4979.7 ± 7.5kJ/mol were found for K_0.85(H₃O)_0.15Fe_2.5(SO₄)₂(OH)_4.5(H₂O)_1.5, K_0.85(H₃O)_0.15Al_2.5(SO₄)₂(OH)_4.5(H₂O)_1.5, Na_0.7(H₃O)_0.3Fe_2.7(SO₄)₂(OH)_5.1(H₂O)_0.9, and Na_0.7(H₃O)_0.3Al_2.7(SO₄)₂(OH)_5.1(H₂O)_0.9 respectively. To our knowledge, this is the first experimentally-based report of ΔH°_f for such nonstoichiometric alunite and natroalunite samples. These thermodynamic data should prove helpful to study, under given conditions, the partitioning of Fe and Al between the solids and aqueous solution.     

Neodymium isotopic variations in Northwest Pacific waters

Four vertical profiles of the concentration and isotopic composition of Nd in seawater were obtained in the western North Pacific. Two profiles from the Kuroshio Current regime showed congruently that although the Nd concentration increases gradually with depth, its isotopic composition varies significantly with depth depending upon the water mass occupying the water column. The high-salinity Kuroshio waters originating from the North Pacific Tropical Water (NPTW) carry the least radiogenic Nd (?_Nd = −7.4 to −8.7) to this region at ∼250 m from the western margin continental shelves, most likely from the East China Sea. The Nd isotopic compositions in the North Pacific Intermediate Water (NPIW) that occurs at 600 to 1000 m in the subtropical region are fairly uniform at ?_Nd = −3.7. The profile data from the ∼38° to 40°N Kuroshio/Oyashio mixed water region off Sanriku of Honshu, Japan, also suggest that the newest NPIW with ?_Nd = −3.2 is formed there by the mixing of various source waters, and the radiogenic component of Nd is derived mainly from the Oyashio waters.In the Pacific Deep Water (PDW) below ∼1000 m, the Nd isotopic composition is neither vertically nor horizontally homogeneous, suggesting that it serves as a useful tracer for sluggish deep water circulation as well. Two profiles from the Izu-Ogasawara Trench showed a minimum ?_Nd value at ∼2000 m, suggesting that there exists a horizontal advective flow in the vicinity of Honshu, Japan. There is some evidence from other chemical properties to support this observation. The waters below 4000 m including those within the trench in the subtropical region have ?_Nd values of around −5, suggesting that the deep waters are fed from the south along the western boundary, ultimately from the Antarctic Bottom Water (AABW) in the South Pacific. This extends up to ∼40°N along the Japanese Islands. In the subarctic region (>∼42°N), the waters have more radiogenic Nd with ?_Nd > −4.0 throughout the water column, presumably due to the supply of Nd by weathering in such igneous provinces as the Kuril-Kamchatska-Aleutian Island chain. The lateral inhomogeneity of the Nd isotopic composition in PDW suggests that there may be different circulation and mixing regimes in the North Pacific Basin.     

Reversibility of radiocaesium sorption on illite

Adsorption of trace amounts of radiocaesium on NH₄-, K-, and Na-saturated Fithian illite and subsequent desorption by 1 M NH₄ showed that a substantial amount of radiocaesium (44%, 46%, and 91% for NH₄-, K-, and Na-illite, respectively) cannot be desorbed after only 5 min of adsorption. Our results suggest that this instantaneous fixation is caused by the collapse of the frayed edges of the clay mineral and the relatively high concentration of radiocaesium building up in solution in the batch desorption experiments. Consequently, commonly applied high-NH₄ extractions underestimate truly exchangeable amounts of radiocaesium in soils and sediments containing illitic clay minerals. The rate of desorption of trace amounts of radiocaesium from the solids using high NH₄ or Cs concentrations has a half-life of about 2 yr, reflecting radiocaesium desorption from (partially) collapsed interlayers. Extraction of radiocaesium from illite after 5 min of contact time with a Cs-selective adsorbent or a 1 × 10^-6 M CsCl solution shows that 100% of the bound radiocaesium is readily available. The desorption rate in the presence of a Cs-selective adsorbent has a half-life of about 0.2 yr. Desorption of radiocaesium from illite using different ammonium concentrations shows that radiocaesium partitioning follows reversible ion-exchange theory if the NH₄ concentration is below 1 × 10^-4 M, and sufficient time (weeks) is allowed for the reaction to proceed. Thus, radiocaesium sorption reversibility in the natural environment is much higher than generally assumed, and equilibrium solid/liquid partitioning may be assumed for the long-term modelling of radiocaesium mobility in the natural environment. In the particular case of anoxic freshwater sediments with very high NH₄ concentrations in the pore waters (up to several mmol.L^-1), collapse of the frayed edges of illite may occur, influencing radiocaesium partitioning. If collapse occurs before radiocaesium adsorbs to illite, high caesium sorption reversibility as measured by high-NH₄ extraction can be expected because further collapse of the frayed edges during the extraction procedure will be limited. This effect has indeed been observed earlier in the extraction of radiocaesium from anoxic freshwater sediments with high-NH₄ solutions and was as yet unexplained.     

Precise and accurate isotopic measurements using multiple-collector ICPMS

New techniques of isotopic measurements by a new generation of mass spectrometers equipped with an inductively-coupled-plasma source, a magnetic mass filter, and multiple collection (MC-ICPMS) are quickly developing. These techniques are valuable because of (1) the ability of ICP sources to ionize virtually every element in the periodic table, and (2) the large sample throughout. However, because of the complex trajectories of multiple ion beams produced in the plasma source whether from the same or different elements, the acquisition of precise and accurate isotopic data with this type of instrument still requires a good understanding of instrumental fractionation processes, both mass-dependent and mass-independent. Although physical processes responsible for the instrumental mass bias are still to be understood more fully, we here present a theoretical framework that allows for most of the analytical limitations to high precision and accuracy to be overcome. After a presentation of unifying phenomenological theory for mass-dependent fractionation in mass spectrometers, we show how this theory accounts for the techniques of standard bracketing and of isotopic normalization by a ratio of either the same or a different element, such as the use of Tl to correct mass bias on Pb. Accuracy is discussed with reference to the concept of cup efficiencies. Although these can be simply calibrated by analyzing standards, we derive a straightforward, very general method to calculate accurate isotopic ratios from dynamic measurements. In this study, we successfully applied the dynamic method to Nd and Pb as examples. We confirm that the assumption of identical mass bias for neighboring elements (notably Pb and Tl, and Yb and Lu) is both unnecessary and incorrect. We further discuss the dangers of straightforward standard-sample bracketing when chemical purification of the element to be analyzed is imperfect. Pooling runs to improve precision is acceptable provided the pooled measurements are shown to be part of a single population. Second-order corrections seem to be able to improve the precision on ¹⁴³Nd/¹⁴⁴Nd measurements. Finally, after discussing a number of potential pitfalls, such as the consequence of peak shape, correlations introduced by counting statistics, and the effect of memory on double-spike methods, we describe an optimal strategy for high-precision and accurate measurements by MC-ICPMS, which involves the repetitive calibration of cup efficiencies and rigorous assessment of mass bias combined with standard-sample bracketing. We suggest that, when these simple guidelines are followed, MC-ICPMS is capable of producing isotopic data precise and accurate to better than 15 ppm.     

Rates of weathering rind formation on Costa Rican basalt

Weathering rind thicknesses were measured on ∼ 200 basaltic clasts collected from three regionally extensive alluvial fill terraces (Qt 1, Qt 2, and Qt 3) preserved along the Pacific coast of Costa Rica. Mass balance calculations suggest that conversion of unweathered basaltic core minerals (plagioclase and augite) to authigenic minerals in the porous rind (kaolinite, allophane, gibbsite, Fe oxyhydroxides) is iso-volumetric and Ti and Zr are relatively immobile. The hierarchy of cation mobility (Ca ≈ Na > K ≈ Mg > Si > Al > Fe ≈ P) is similar to other tropical weathering profiles and is indicative of differential rates of mineral weathering (anorthite > albite ≈ hypersthene > orthoclase ? apatite). Alteration profiles across the cm-thick rinds document dissolution of plagioclase and augite and the growth of kaolinite, with subsequent dissolution of kaolinite and precipitation of gibbsite as weathering rinds age. The rate of weathering rind advance is evaluated using a diffusion-limited model which predicts a parabolic rate law for weathering rind thickness, r_r, as a function of time, t(r_r =), and an interface-limited model which predicts a linear rate law for weathering rind thickness as a function of time (r_r = k_appt). In these rate laws, κ is a diffusion parameter and k_app is an apparent rate constant. The rate of advance is best fit by the interface model.Terrace exposures are confined to the lower reaches of streams draining the Pacific slope near the coast where the stream gradient is less than ∼3 m/km, and terrace deposition is influenced by eustatic sea level fluctuations. Geomorphological evidence is consistent with terrace deposition coincident with sea level maxima when the stream gradient would be lowest. Assigning the most weathered regionally extensive terrace Qt 1 (mean rind thickness 6.9 ± 0. 6cm) to oxygen isotope stage (OIS) 7 (ca. 240 ka), and assuming that at time = 0 rind thickness = 0, it is inferred that terrace Qt 2 (r_r = 2.9 ± 0.1 cm) is coincident with stage 5e (ca. 125 ka) and that Qt 3 (r_r = 0.9 ± 0.1 cm) is consistent with OIS 3 (ca. 37 ka). These assignments yield a value of k_app of 8.6 × 10⁻¹³ cm s⁻¹ (R² = 0.99). Only this value satisfies both the existing age controls and yields ages coincident with sea level maxima. Using this value, elemental weathering release fluxes across a weathering rind from Qt 2 range from 6.0 × 10⁻⁹ mol Si m⁻² s⁻¹ to 2.5 × 10⁻¹¹ mol K m⁻² s⁻¹. The rate of rind advance for the Costa Rican terraces is 2.8 × 10⁻⁷ m yr⁻¹. Basalt rind formation rates in lower temperature settings described in the literature are also consistent with interface-controlled weathering with an apparent activation energy of about 50 kJ mol⁻¹. Rates of rind formation in Costa Rica are an order of magnitude slower than reported for global averages of soil formation rates.     

Oxygen isotope fractionation in synthetic magnesian calcite

Mg-bearing calcite was precipitated at 25°C in closed system free-drift experiments from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The chemical and isotope composition of the solution and precipitate were investigated during time course experiments of 24-h duration. Monohydrocalcite and calcite precipitated early in the experiments (<8 h), while Mg-calcite was the predominant precipitate (>95%) thereafter. Solid collected at the end of the experiments displayed compositional zoning from pure calcite in crystal cores to up to 23 mol% MgCO₃ in the rims. Smaller excursions in Mg were superimposed on this chemical record, which is characteristic of oscillatory zoning observed in synthetic and natural solid-solution carbonates of differing solubility. Magnesium also altered the predominant morphology of crystals over time from the {104} to {100} and {110} growth forms.The oxygen isotope fractionation factor for the magnesian-calcite-water system (as 10³lnα_Mg-cl-H2O) displayed a strong dependence on the mol% MgCO₃ in the solid phase, but quantification of the relationship was difficult due to the heterogeneous nature of the precipitate. Considering only the Mg-content and δ¹⁸O values for the bulk solid, 10³lnα_Mg-cl-H2O increased at a rate of 0.17 ± 0.02 per mol% MgCO₃; this value is a factor of three higher than the single previous estimate (Tarutani T., Clayton R.N., and Mayeda T. K. (1969) The effect of polymorphims and magnesium substitution on oxygen isotope fractionation between calcium carbonate and water. Geochim. Cosmochim. Acta 33, 987-996). Nevertheless, extrapolation of our relationship to the pure calcite end member yielded a value of 27.9 ± 0.02, which is similar in magnitude to published values for the calcite-water system. Although no kinetic effect was observed on 10³lnα_Mg-cl-H2O for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol · m⁻² · h⁻¹, it was impossible to disentangle the potential effect(s) of precipitation rate and Mg-content on 10³lnα_Mg-cl-H2O due to the heterogeneous nature of the solid.The results of this study suggest that paleotemperatures inferred from the δ¹⁸O values of high magnesian calcite (>10 mol% MgCO₃) may be significantly underestimated. Also, the results underscore the need for additional experiments to accurately characterize the effect of Mg coprecipitation on the isotope systematics of calcite from a chemically homogeneous precipitate or a heterogeneous material that is analyzed at the scale of chemical and isotopic zonation.     

Strontium incorporation into calcite generated by bacterial ureolysis

Strontium incorporation into calcite generated by bacterial ureolysis was investigated as part of an assessment of a proposed remediation approach for ⁹⁰Sr contamination in groundwater. Urea hydrolysis produces ammonium and carbonate and elevates pH, resulting in the promotion of calcium carbonate precipitation. Urea hydrolysis by the bacterium Bacillus pasteurii in a medium designed to mimic the chemistry of the Snake River Plain Aquifer in Idaho resulted in a pH rise from 7.5 to 9.1. Measured average distribution coefficients (D_EX) for Sr in the calcite produced by ureolysis (0.5) were up to an order of magnitude higher than values reported in the literature for natural and synthetic calcites (0.02-0.4). They were also higher than values for calcite produced abiotically by ammonium carbonate addition (0.3). The precipitation of calcite in these experiments was verified by X-ray diffraction. Time-of-flight secondary ion mass spectrometry (ToF SIMS) depth profiling (up to 350 nm) suggested that the Sr was not merely sorbed on the surface, but was present at depth within the particles. X-ray absorption near edge spectra showed that Sr was present in the calcite samples as a solid solution. The extent of Sr incorporation appeared to be driven primarily by the overall rate of calcite precipitation, where faster precipitation was associated with greater Sr uptake into the solid. The presence of bacterial surfaces as potential nucleation sites in the ammonium carbonate precipitation treatment did not enhance overall precipitation or the Sr distribution coefficient. Because bacterial ureolysis can generate high rates of calcite precipitation, the application of this approach is promising for remediation of ⁹⁰Sr contamination in environments where calcite is stable over the long term.     

Barite deposition resulting from phototrophic sulfide-oxidizing bacterial activity

Barite (BaSO₄) deposits generally arise from mixing of soluble barium-containing fluids with sulfate-rich fluids. While the role of biological processes in modulating barium solubility has been shown, no studies have shown that the biological oxidation of sulfide to sulfate leads to barite deposition. Here we present an example of microbially mediated barite deposition in a continental setting. A spring in the Anadarko Basin of southwestern Oklahoma produces water containing abundant barium and sulfide. As emergent water travels down a stream to a nearby creek, sulfate concentration increases from 0.06 mM to 2.2 mM while Ba²⁺ concentration drops from 0.4 mM to less than 7 μM. Stable isotope analysis, microbial activity studies, and in situ experiments provide evidence that as sulfide-rich water flows down the stream, anaerobic, anoxygenic, phototrophic bacteria play a dominant role in oxidizing sulfide to sulfate. Sulfate then precipitates with Ba²⁺ producing barite as travertine, cements, crusts, and accumulations on microbial mats. Our studies suggest that phototrophic sulfide oxidation and concomitant sulfur cycling could prove to be important processes regulating the cycling of barium in continental sulfur-containing systems.     

Transport-controlled kinetics of dissolution and precipitation in the sediments under alkaline and saline conditions

Over 1.6 million liters of radioactive, high-temperature, Al-rich, alkaline and saline high-level waste (HLW) fluids were accidentally discharged from tank leaks onto the sediments at the Hanford Site, Washington. In order to better understand processes that might occur during the migration of HLW through sediments and to estimate their extents, we studied the effects of Al-rich, alkaline and saline solutions on soil mineral dissolution and precipitation during reactive transport. Metal- and glass-free systems were used to conduct miscible-displacement experiments at 50 °C under CO₂ and O₂ free conditions. Results showed significant release of Si, K, Al, Fe, Ca, Mg, and Ba into the aqueous phase. The transport-controlled release of these elements was time dependent as evidenced by its extent varying with the fluid residence time. Silica initial dissolution rates (6.08 × 10^-11 and 5.38 × 10^-13 mol m^-2 s^-1) increased with base concentration, decreased with Al concentration, and decreased with fluid residence time. Aluminum precipitation rates varied in the range from 0.44 to 1.07 × 10^-6 mol s^-1 and were faster in these column experiments than in previous batch studies. The initial rate constant of Al precipitation reaction was 0.07 h^-1 (half-life of 9.9 h at about 3 PV); it increased up to 0.137 h^-1 (half-life of 5.1 h at about 20 PV). The precipitates identified with SEM and suggested from the modeling results were mainly NO₃-cancrinite. SEM analyses also indicated the formation of sodalite when Al was not present in the leaching solution. In addition, results from modeling suggested the precipitation of brucite, goethite and gibbsite; the latter may precipitate in the presence of high Al concentrations. Aqueous and solid phase transformations caused by base-induced dissolution and subsequent secondary phases precipitation should be important determinants of the fate of contaminants and radionuclides in the vadose zone under alkaline and saline conditions.     

Imaging and mechanical property measurements of kerogen via nanoindentation

Most analyses of kerogens rely on samples that have been isolated by dissolving the rock matrix. The properties of the kerogen before and after such isolation may be different and all sample orientation information is lost. We report a method of measuring kerogen mechanical properties in the rock matrix without isolation. An atomic force microscope (AFM) based nanoindenter is used to measure the hardness and reduced modulus of the kerogen within Woodford shale. The same instrument also provides useful images of polished rock sections on a submicrometer scale. Measurements were carried out both parallel and perpendicular to the bedding plane.     

Carbon and hydrogen isotope fractionation by moderately thermophilic methanogens

A series of laboratory studies were conducted to increase understanding of stable carbon (¹³C/¹²C) and hydrogen (D/H) isotope fractionation arising from methanogenesis by moderately thermophilic acetate- and hydrogen-consuming methanogens. Studies of the aceticlastic reaction were conducted with two closely related strains of Methanosaeta thermophila. Results demonstrate a carbon isotope fractionation of only 7‰ (α = 1.007) between the methyl position of acetate and the resulting methane. Methane formed by this process is enriched in ¹³C when compared with other natural sources of methane; the magnitude of this isotope effect raises the possibility that methane produced at elevated temperature by the aceticlastic reaction could be mistaken for thermogenic methane based on carbon isotopic content. Studies of H₂/CO₂ methanogenesis were conducted with Methanothermobacter marburgensis. The fractionation of carbon isotopes between CO₂ and CH₄ was found to range from 22 to 58‰ (1.023 ≤ α ≤ 1.064). Greater fractionation was associated with low levels of molecular hydrogen and steady-state metabolism. The fractionation of hydrogen isotopes between source H₂O and CH₄ was found to range from 127 to 275‰ (1.16 ≤ α ≤ 1.43). Fractionation was dependent on growth phase with greater fractionation associated with later growth stages. The maximum observed fractionation factor was 1.43, independent of the δD-H₂ supplied to the culture. Fractionation was positively correlated with temperature and/or metabolic rate. Results demonstrate significant variability in both hydrogen and carbon isotope fractionation during methanogenesis from H₂/CO₂. The relatively small fractionation associated with deuterium during H₂/CO₂ methanogenesis provides an explanation for the relatively enriched deuterium content of biogenic natural gas originating from a variety of thermal environments. Results from these experiments are used to develop a hypothesis that differential reversibility in the enzymatic steps of the H₂/CO₂ pathway gives rise to variability in the observed carbon isotope fractionation. Results are further used to constrain the overall efficiency of electron consumption by way of the hydrogenase system in M. marburgensis, which is calculated to be less than 55%.     

Long-range vs. short-range ordering in synthetic Cr-substituted goethites

A series of Cr-substituted goethites with (Cr:Fe molar ratio up to 0.12) were prepared. Thermal analysis of the solids indicates the formation of cation-deficient compounds that are more stable towards the transformation to hematite as the Cr content increases. Powder X-ray diffraction (PXRD) and extended X-ray absorption fine structure (EXAFS) techniques were used to assess the structural characteristics of the whole series of the substituted solids. XRD patterns demonstrate that the order around Fe remains typical of a goethite-like structure. Rietveld refinement of X-ray diffraction data indicates that the incorporation of Cr causes a slight decrease in the cell volume with the c-cell parameter following the Vegard’s law. This decrease is accompanied by changes in opposite directions of the various Me-Me distances. EXAFS spectra at the Fe K-edge indicate that the local order around the Fe atom changes slightly upon Cr substitution: Measurements in the Cr K-edge show that the Cr environment remains unchanged in the whole series. All the observed trends in both average Rietveld and local EXAFS distances can be traced back to the differences in the coordination polyhedra around Cr and Fe. The polyhedron around Cr is more symmetric and can be described as Cr(OH_0.5)₆ as opposed to the polyhedron around Fe that contains two distinct sets of ligands, FeO₃(OH)₃. The effects caused by substitution are governed by this difference, rather than by the smaller size of Cr(III) as compared to Fe(III). Simultaneous use of XAS and Rietveld refinement of XRD data permits tracing the trends in the average long range ordering (Me-Me distances) to local changes in distances and angles when Cr³⁺ substitutes Fe³⁺ in goethite. Complex changes in the various interatomic distances and angles may result in deceivingly simple long-range trends. These trends are therefore of limited value as probes for the atomic scale changes. On the other hand, XAS provide direct information on the fundamental, atomic-scale changes.