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.     

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.     

Experimental investigation of the effects of temperature and ionic strength on Mo isotope fractionation during adsorption to manganese oxides

The Mo stable isotope system is being applied to study changes in ocean redox. Such applications implicitly assume that Mo isotope fractionation in aqueous systems is relatively insensitive to frequently changing environmental variables such as temperature (T) and ionic strength (I). A major driver of fractionation is the adsorption of Mo to Mn oxyhydroxide surfaces [Barling J. and Anbar A. D. (2004) Molybdenum isotope fractionation during adsorption by manganese oxides. Earth Planet. Sci. Lett.217(3-4), 315-329]. Here, we report the results of experiments that determine the extent to which Mo isotope fractionation during adsorption of Mo to the Mn oxyhydroxide mineral birnessite is sensitive to T and I. The results are compared to new predictions from quantum chemical computations. We measured fractionation from 1 to 50 °C at I = 0.1 m and found that Δ^97/95Mo_{dissolved-adsorbed} varies from 1.9‰ to 1.6‰ over this temperature range. Experiments were also performed at 25 °C in synthetic seawater (I = 0.7); fractionation at this condition was the same within analytical error as in low ionic strength experiments. These findings confirm that the Mo isotope fractionation during adsorption to Mn oxyhydroxides is relatively insensitive to variations and T and I over environmentally relevant ranges. To relate these findings to potential mechanisms of Mo isotope fractionation, we also report results for density functional theory computations of the fractionation between and various possible structures of molybdic acid as a function of temperature. Because no plausible species fractionates from with a magnitude matching the experiments, we are left with three possibilities to explain the fractionation: (1) solvation effects on the vibrational frequencies of aqueous species considered thus far are significant, such that our calculations in vacuo yield inaccurate fractionations; (2) a trace aqueous species not yet considered fractionates from and then adsorbs to birnessite; or (3) a surface complex not present in solution forms on birnessite in which Mo is not tetrahedrally coordinated. Our findings help validate assumptions underlying paleoceanographic applications of the Mo isotope system and also lead us closer to understanding the mechanism of isotope fractionation during adsorption of Mo to Mn oxyhydroxides.     

The iron-isotope fractionation dictated by the carboxylic functional: An ab-initio investigation

The ground-state geometries, electronic energies and vibrational properties of carboxylic complexes of iron were investigated both in vacuo and under the effect of a reaction field, to determine thermodynamic properties of iron-acetates and the role of the carboxylic functional on the isotopic imprinting of this metal in metalorganic complexation. The electronic energy, zero point corrections and thermal corrections of these substances at variational state were investigated at the DFT/B3LYP level of theory with different basis set expansions and the effect of the reaction field on the variational structures was investigated through the Polarized Continuun Model. Thermochemical cycle calculations, combined with solvation energy calculations and appropriate scaling from absolute to conventional properties allowed to compute the Gibbs free energy of formation from the elements of the investigated aqueous species and to select the best procedure to be applied in the successive vibrational analysis. The best compliance with the few existing thermodynamic data for these substances was obtained by coupling the gas phase calculations at DFT/B3LYP level with the [6-31G(d,p)]-[6-31G+(d,p)] (for cations and neutral molecules - anions; respectively) with solvation calculations adopting atomic radii optimized for the HF/6-31G(d) level of theory (UAHF). A vibrational analysis conducted on ⁵⁴Fe, ⁵⁶Fe, ⁵⁷Fe and ⁵⁸Fe gaseous isotopomers yielded reduced partition function ratios which increased not only with the nominal valence of the central cation, as expected, but, more importantly, with the extent of the complexation operated by the organic functional. Coupling thermodynamic data with separative effects it was shown that this last is controlled, as expected, by the relative bond strength of the complex in both aggregation states. Through the Integral Equation Formalism of the Polarized Continuum Model (IEFPCM) the effect of the ionic strength of the solution and of a T-dependent permittivity on the energy and separative effects of the solvated metalorganic complexes were analyzed in detail. The solvent effect in the standard state (hypothetical one-molal solution referred to infinite dilution; T = 298.15 K, P = 1 bar) is a limited reduction of the separative effects of all the isotopomeric couples. With an increase in T (and the concomitant decrease in the dielectric constant of the solvent) this effect diminishes progressively.     

Thermodynamic model for arsenic speciation in sulfidic waters: A novel use of ab initio computations

Recent discoveries demonstrate that the chemistry of arsenic in sulfidic waters is much more complex that previously believed. One implication is that all earlier thermodynamic data on stabilities of As thioanions require revision. Previously used experimental approaches for determining As thioanion stabilities may be inadequate to deal with the full range of complexity. Here we use computational as well as empirical information to construct a provisional model for equilibrium As thioanion distributions in sulfidic waters. Whereas previous authors have argued for either As(III) or As(V) thioanions, the new model predicts that both are important and can occur simultaneously under commonly encountered pH and ΣS^−II conditions. At the order of magnitude level, the model reasonably predicts the solubility of As₂S₃ in sulfidic solutions, provides tentative peak assignments for published Raman spectroscopic data and plausibly accounts for how sulfide modifies the bacterial toxicity of As. The model yields a thermodynamic justification for how sulfide, which is usually regarded as a reducing agent, can counter-intuitively drive oxidation of As(III) to As(V), as has been observed both in the laboratory and in the field. Despite its uncertain accuracy, the model serves as a useful source of new, testable hypotheses about As geochemistry and highlights crucial experimental data needs.     

Exurbia from the bottom-up: Confronting empirical challenges to characterizing a complex system

We describe empirical results from a multi-disciplinary project that support modeling complex processes of land-use and land-cover change in exurban parts of Southeastern Michigan. Based on two different conceptual models, one describing the evolution of urban form as a consequence of residential preferences and the other describing land-cover changes in an exurban township as a consequence of residential preferences, local policies, and a diversity of development types, we describe a variety of empirical data collected to support the mechanisms that we encoded in computational agent-based models. We used multiple methods, including social surveys, remote sensing, and statistical analysis of spatial data, to collect data that could be used to validate the structure of our models, calibrate their specific parameters, and evaluate their output. The data were used to investigate this system in the context of several themes from complexity science, including have (a) macro-level patterns; (b) autonomous decision making entities (i.e., agents); (c) heterogeneity among those entities; (d) social and spatial interactions that operate across multiple scales and (e) nonlinear feedback mechanisms. The results point to the importance of collecting data on agents and their interactions when producing agent-based models, the general validity of our conceptual models, and some changes that we needed to make to these models following data analysis. The calibrated models have been and are being used to evaluate landscape dynamics and the effects of various policy interventions on urban land-cover patterns.     

Late Precambrian rifting in the geological history of the western slope of the South Urals

The rift-related geodynamic setting of the Late Precambrian geological evolution on the western slope of the South Urals is reconstructed on the basis of localization of lithotectonic complexes of this age, their formation conditions, and the geochemistry of rocks. The Early Riphean stage comprises accumulation of coarse-clastic rocks intercalating with alkaline volcanic rocks of the Navysh Complex, which is a constituent of the Ai Formation, and emplacement of doleritic and picritic intrusions of the Shuida Complex and melanocratic dolerite and gabbrodolerite of the Yusha Complex. The Middle Riphean stage is characterized by wide-spread coarse-clastic terrigenous rocks of the Mashak Formation that intercalate with volcanic rocks of the bimodal basalt-rhyolite association, the Berdyaush pluton of rapakivi granite, the Kusa-Kopan layered intrusive complex, the Lapyshta Complex of dolerites and picrites, and numerous occurrences of gabbrodolerites. The terrigenous rocks of the Vendian stage include conglomerate, gravelstone, and sandstone of the Asha Group, while igneous rocks comprise alkaline volcanics of the Arsha Complex, alkali gabbroids of the Miseli Complex, and melanocratic syenite of the Avashla Complex. The geological evolution of the region is distinguished by local (failed or aborted) rifting. The occurrence of lithotectonic complexes is controlled by dynamic conditions of rifting. A certain inheritance in the evolution may be traced for the Early and Middle Riphean and partly for the Late Riphean and Vendian.     

Planktonic foraminifers in the southern Bering Sea: Changes in composition and productivity during the Late Pleistocene-Holocene

Taxonomic composition and distribution of planktonic foraminifers are studied in section of Core GC-11 that penetrated through Upper Quaternary sediments of the Bowers Ridge western slope, the southern Bering Sea. As is shown, structure of foraminiferal assemblage and productivity have varied substantially during the last 32000 calendar years in response to changes in surface water temperatures and water mass circulation in the northern part of the Pacific, the Bering Sea included. The productivity was maximal during deglaciation epoch, being notably lower in the Holocene and minimal at the glaciation time.     

The presence of anomalous trace element levels in present day Jamaican soils and the geochemistry of Late-Miocene or Pliocene phosphorites

High levels of Cd and Zn in Jamaican soils observed in geochemical surveys are related to the presence of phosphorites of possible Late-Miocene or Pliocene age. The trace element and REE geochemistry of the phosphorites, together with SEM studies, indicate a guano origin for the phosphorites. No specific host minerals for Cd could be identified in the fossiliferous phosphorite which is characterized by uniquely high levels of Cd, Zn, Ag, Be, U and Y. However, in the soil Cd is present in lithiophorite and a complex history of pedological development is preserved in the aluminous–goethite present in the soil. The unique guano signature is preserved in the soil despite the fact that guanos themselves have either not been observed or have been destroyed by continuing karst and soil development. The phosphorite geochemical signature can be traced in the data of a 1988 island-wide soil geochemical survey, identifying areas where the Palaeo-environment that supported bird ‘rookeries’ existed in the Late-Miocene or Pliocene.