Surface complexation of U(VI) on goethite (α-FeOOH)

Sorption of U(VI) to goethite is a fundamental control on the mobility of uranium in soil and groundwater. Here, we investigated the sorption of U on goethite using EXAFS spectroscopy, batch sorption experiments and DFT calculations of the energetics and structures of possible surface complexes. Based on EXAFS spectra, it has previously been proposed that U(VI), as the uranyl cation , sorbs to Fe oxide hydroxide phases by forming a bidentate edge-sharing (E2) surface complex, >Fe(OH)₂UO₂(H₂O)_n. Here, we argue that this complex alone cannot account for the sorption capacity of goethite (α-FeOOH). Moreover, we show that all of the EXAFS signal attributed to the E2 complex can be accounted for by multiple scattering. We propose that the dominant surface complex in CO₂-free systems is a bidentate corner-sharing (C2) complex, (>FeOH)₂UO₂(H₂O)₃ which can form on the dominant {101} surface. However, in the presence of CO₂, we find an enhancement of UO₂ sorption at low pH and attribute this to a (>FeO)CO₂UO₂ ternary complex. With increasing pH, U(VI) desorbs by the formation of aqueous carbonate and hydroxyl complexes. However, this desorption is preceded by the formation of a second ternary surface complex (>FeOH)₂UO₂CO₃. The three proposed surface complexes, (>FeOH)₂UO₂(H₂O)₃, >FeOCO₂UO₂, and (>FeOH)₂UO₂CO₃ are consistent with EXAFS spectra. Using these complexes, we developed a surface complexation model for U on goethite with a 1-pK model for surface protonation, an extended Stern model for surface electrostatics and inclusion of all known UO₂-OH-CO₃ aqueous complexes in the current thermodynamic database. The model gives an excellent fit to our sorption experiments done in both ambient and reduced CO₂ environments at surface loadings of 0.02-2.0 wt% U.     

Seasonal patterns of carbon chemistry and isotopes in tufa depositing groundwaters of southwestern Japan

Annually laminated carbonates, known as tufas, commonly develop in limestone areas and typically record seasonal patterns of oxygen- and carbon-isotope compositions. δ¹⁸O values are principally controlled by seasonal changes of water temperature, whereas δ¹³C values are the result of complex reactions among the gaseous, liquid, and solid sources of carbon in the system. We examined the processes that cause the seasonal patterns of δ¹³C in groundwater systems at three tufa-depositing sites in southwestern Japan by applying model calculations to geochemical data. Underground inorganic carbon species are exchanged with gaseous CO₂, which is mainly introduced to the underground hydrological system by natural atmospheric ventilation and by diffusion of soil air. These processes control the seasonal pattern of δ¹³C, which is low in summer and high in winter. Among the three sites we investigated, we identified two extreme cases of the degree of carbon exchange between liquid and gaseous phases. For the case with high radiocarbon composition (Δ¹⁴C) and low pCO₂, there was substantial carbon exchange because of a large contribution of atmospheric CO₂ and a small water mass. For the other extreme case, which was characterized by low Δ¹⁴C and high pCO₂, the contribution of atmospheric CO₂ was small and the water mass was relatively large. Our results suggest that at two of the three sites water residence time within the soil profile was longer than 1 year. Our results also suggested a short residence time (less than 1 year) of water in the soil profile at the site with the smallest water mass, which is consistent with large seasonal amplitude of the springwater temperature variations. The Δ¹⁴C value of tufas is closely related to the hydrological conditions in which they are deposited. If the initial Δ¹⁴C value of a tufa-depositing system is stable, ¹⁴C-chronology can be used to date paleo-tufas.     

Temperature effects on non-bridging oxygen and aluminum coordination number in calcium aluminosilicate glasses and melts

Configurational changes with temperature are important for the thermodynamic and transport properties of most aluminosilicate melts, but in general are not well understood. Here, we present high-resolution ²⁷Al and ¹⁷O NMR data on several calcium aluminosilicate glasses prepared with varying quench rates and thus with fictive temperatures that span ranges up to about 200 K. In all compositions the content of five-coordinated aluminum increases with fictive temperature, in agreement with recent high temperature NMR data on melts. In a glass of CaAl₂Si₂O₈ (“anorthite”) composition, the content of non-bridging oxygens also increases with temperature; however this effect was not observed in a sample with a much higher CaO/Al₂O₃ ratio. We present a consistent notation for reactions among structural species in these systems that clarify why in some cases, high-coordinated network cations may appear on the same side of the reaction, while in others they occur on the opposite sides: the key difference is in accounting for all coordination changes for oxygens. Mixing of non-bridging oxygens and of high-coordinated aluminum make significant contributions to the overall configurational entropy and heat capacity of the melts, as does the mixing of various bridging oxygens and of tetrahedral network cations. Other, less well known, types of increase in disorder with temperature may be important as well.     

Kinetics of Cd sorption, desorption and fixation by calcite: A long-term radiotracer study

Time-dependent sorption and desorption of Cd on calcite was studied over 210 days utilizing ¹⁰⁹Cd as a tracer to distinguish between ‘labile’ and ‘non-labile’ forms of sorbed Cd. Stabilizing the calcite suspensions for 12 months under atmospheric P_CO2 and controlled temperature was necessary to reliably follow Cd dynamics following initial sorption. Results revealed time-dependant Cd sorption and marked desorption hysteresis by calcite under environmentally relevant conditions. Data obtained were fitted to a first-order kinetic model and a concentric shell diffusion model. Both models described the progressive transfer of Cd²⁺ to a less reactive form within calcite and subsequent desorption of Cd subject to different initial contact times. The kinetic model provided a better fit to the combined sorption and desorption data (R² = 0.992). It differentiates between two ‘pools’ of sorbed Cd²⁺ on calcite, ‘labile’ and ‘non-labile’, in which labile sorbed Cd is in immediate equilibrium with the free Cd²⁺ ion activity in solution whereas non-labile Cd is kinetically restricted. For the diffusion model (R² = 0.959), the rate constants describing Cd dynamics in calcite produced a half-life for Cd desorption of ∼175 d, for release to a ‘zero-sink’ solution. Results from this study allow comment on the likely mechanisms occurring at the calcite surface following long-term Cd sorption.     

The age of the martian meteorite Northwest Africa 1195 and the differentiation history of the shergottites

Samarium-neodymium isotopic analyses of unleached and acid-leached mineral fractions from the recently identified olivine-bearing shergottite Northwest Africa 1195 yield a crystallization age of 347 ± 13 Ma and an value of +40.1 ± 0.9. Maskelynite fractions do not lie on the Sm-Nd isochron and appear to contain a martian surface component with low ¹⁴⁷Sm/¹⁴⁴Nd and ¹⁴³Nd/¹⁴⁴Nd ratios that was added during shock. The Rb-Sr system is disturbed and does not yield an isochron. Terrestrial Sr appears to have affected all of the mineral fractions, although a maximum initial ⁸⁷Sr/⁸⁶Sr ratio of 0.7016 is estimated by passing a 347 Ma reference line through the maskelynite fraction that is least affected by contamination. The high initial value and the low initial ⁸⁷Sr/⁸⁶Sr ratio, combined with the geologically young crystallization age, indicate that Northwest Africa 1195 is derived from a source region characterized by a long-term incompatible-element depletion.The age and initial Sr and Nd isotopic compositions of Northwest Africa 1195 are very similar to those of Queen Alexandra Range 94201, indicating these samples were derived from source regions with similar Sr-Nd isotopic systematics. These similarities suggest that these two meteorites share a close petrogenetic relationship and might have been erupted from a common volcano. The meteorites Yamato 980459, Dar al Gani 476, Sayh al Uhaymir 005/008, and Dhofar 019 also have relatively old ages between 474 and 575 Ma and trace element and/or isotopic systematics that are indicative of derivation from incompatible-element-depleted sources. This suggests that the oldest group of meteorites is more closely related to one another than they are to the younger meteorites that are derived from less incompatible-element-depleted sources. Closed-system fractional crystallization of this suite of meteorites is modeled with the MELTS algorithm using the bulk composition of Yamato 980459 as a parent. These models reproduce many of the major element and mineralogical variations observed in the suite. In addition, the rare earth element systematics of these meteorites are reproduced by fractional crystallization using the proportions of phases and extents of crystallization that are calculated by MELTS. Other shergottites that demonstrate enrichments in incompatible-elements and have evolved Sr and Nd isotopic systematics have some geochemical systematics that are similar to those observed in the depleted group. Most notably, although they exhibit a very limited range of incompatible trace element and isotopic compositions, they have highly variable major element compositions. This is also consistent with evolution from a common mantle source region by variable amounts of fractional crystallization. If this scenario is correct, it suggests that the combined effects of source composition and fractional crystallization are likely to account for the major element, trace element, and isotopic diversity of all shergottites.     

Cr-spinel/olivine and Cr-spinel/liquid nickel partition coefficients from natural samples

The Nernst partition coefficient of nickel (D^Ni) between Cr-spinel and silicate melt in natural systems has been investigated using mid-ocean ridge basalts (MORB) and other volcanic rocks. The Cr-spinel/olivine D^Ni values in volcanic rocks are between 1.2 and 0.3, indicating that the Cr-spinel/liquid D^Ni values vary from slightly higher to significantly lower than the olivine/liquid D^Ni values in natural systems. The Cr-spinel/liquid D^Ni values from the MORB samples vary between 6 and 11, slightly higher than those from the S-bearing experiments of Satari et al. [Satari P., Brenan J. M., Horn I. and McDonough W. F. (2002) Experimental constraints on the sulfide- and chromite-silicate melt partitioning behavior of rhenium and platinum-group elements. Economic Geology97, 385-398]. The results of the MORB samples and the experiments of Satari et al. (2002) indicate a negative correlation between the Cr-spinel/liquid D^Ni and the X_Cr values in Cr-spinels (Cr cation number on the basis of 3 total cations in the spinel structure). Variations of Cr-spinel/liquid D^Ni values with Cr-spinel compositions can be estimated from an empirical equation based on the results of the MORB samples and the experiments by Satari et al. (2002). The choice of Cr-spinel/liquid D^Ni = 10 for numerical modeling by Righter et al. [Righter K., Leeman W. P. and Hervig R. L. (2006) Partitioning of Ni, Co, and V between spinel-structured oxides and silicate melts: importance of spinel composition. Chemical Geology227, 1-25] is reasonable for basaltic systems. For picritic and komatiitic systems a lower value of ∼5 is more appropriate.     

Size scales over which ordinary chondrites and their parent asteroids are homogeneous in oxidation state and oxygen-isotopic composition

Literature data demonstrate that on a global, asteroid-wide scale (plausibly on the order of 100 km), ordinary chondrites (OC) have heterogeneous oxidation states and O-isotopic compositions (represented, respectively, by the mean olivine Fa and bulk Δ¹⁷O compositions of equilibrated samples). Samples analyzed here include: (a) two H5 chondrite Antarctic finds (ALHA79046 and TIL 82415) that have the same cosmic-ray exposure age (7.6 Ma) and were probably within ∼1 km of each other when they were excavated from the H-chondrite parent body, (b) different individual stones from the Holbrook L/LL6 fall that were probably within ∼1 m of each other when their parent meteoroid penetrated the Earth’s atmosphere, and (c) drill cores from a large slab of the Estacado H6 find located within a few tens of centimeters of each other. Our results indicate that OC are heterogeneous in their bulk oxidation state and O-isotopic composition on 100-km-size scales, but homogeneous on meter-, decimeter- and centimeter-size scales. (On kilometer size scales, oxidation state is heterogeneous, but O isotopes appear to be homogeneous.) The asteroid-wide heterogeneity in oxidation state and O-isotopic composition was inherited from the solar nebula. The homogeneity on small size scales was probably caused in part by fluid-assisted metamorphism and mainly by impact-gardening processes (which are most effective at mixing target materials on scales of ?1 m).     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Upwelling, species, and depth effects on coral skeletal cadmium-to-calcium ratios (Cd/Ca)

Skeletal cadmium-to-calcium (Cd/Ca) ratios in hermatypic stony corals have been used to reconstruct changes in upwelling over time, yet there has not been a systematic evaluation of this tracer’s natural variability within and among coral species, between depths and across environmental conditions. Here, coral skeletal Cd/Ca ratios were measured in multiple colonies of Pavona clavus, Pavona gigantea and Porites lobata reared at two depths (1 and 7 m) during both upwelling and nonupwelling intervals in the Gulf of Panama (Pacific). Overall, skeletal Cd/Ca ratios were significantly higher during upwelling than during nonupwelling, in shallow than in deep corals, and in both species of Pavona than in P. lobata. P. lobata skeletal Cd/Ca ratios were uniformly low compared to those in the other species, with no significant differences between upwelling and nonupwelling values. Among colonies of the same species, skeletal Cd/Ca ratios were always higher in all shallow P. gigantea colonies during upwelling compared to nonupwelling, though the magnitude of the increase varied among colonies. For P. lobata, P. clavus and deep P. gigantea, changes in skeletal Cd/Ca ratios were not consistent among all colonies, with some colonies having lower ratios during upwelling than during nonupwelling. No statistically significant relationships were found between skeletal Cd/Ca ratios and maximum linear skeletal extension, δ¹³C or δ¹⁸O, suggesting that at seasonal resolution the Cd/Ca signal was decoupled from growth rate, coral metabolism, and ocean temperature and salinity, respectively. These results led to the following conclusions, (1) coral skeletal Cd/Ca ratios are independent of skeletal extension, coral metabolism and ambient temperature/salinity, (2) shallow P. gigantea is the most reliable species for paleoupwelling reconstruction and (3) the average Cd/Ca record of several colonies, rather than of a single coral, is needed to reliably reconstruct paleoupwelling events.     

Isotopic studies of leaf water. Part 1: A physically based two-dimensional model for pine needles

Understanding of isotopic variations in leaf water is important for reconstruction of paleoclimate and assessment of global biochemical processes. We report here a study of isotopic distributions within a single needle of two pine species, Pinus resinosa Ait and Pinus strobes L., with the objective of understanding how isotopic compositions of leaf water are controlled by environmental and physiological variables. A 2D model was developed to simulate along-leaf isotopic variations and bulk leaf water isotopic compositions. In addition to variables common to all leaf water isotopic models, this 2D model also takes into account the specific geometry and dimensions of pine needles and the isotopic transport in xylem and mesophyll. The model can successfully simulate oxygen isotopic variations along a single needle and averaged over a leaf (bulk leaf water). The simulations suggest that isotopic composition of the bulk leaf water does not always depend only upon the average transpiration rate, which in turn raises questions about using leaf water isotopic values to estimate transpiration rates. An unsuccessful attempt to simulate along-needle hydrogen isotopic variations suggests that certain unknown biological process(es) may not have been incorporated into our 2D model, and if so, it calls for a reevaluation of all other models for hydrogen isotopic simulations of leaf water since they too lack these processes.Existing leaf water isotopic models are reviewed in this work. In particular, we evaluate the most frequently used model, the stomatal boundary layer model (also referred to as the Craig-Gordon model). We point out that discrepancy between the boundary layer model and the measured bulk leaf water seems to depend upon relative humidity. Using our 2D model, we show that this humidity dependency is a result of an interplay between environmental and physiological conditions: if the transpiration rate of plant leaves decreases with increasing relative humidity, our 2D model can reproduce the pattern of isotopic discrepancy between boundary layer model predictions and observations, enabling us to understand better the reason behind this discrepancy.