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.     

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.     

Microbial mass-dependent fractionation of chromium isotopes

Mass-dependent fractionation of Cr isotopes occurs during dissimilatory Cr(VI) reduction by Shewanella oneidensis strain MR-1. Cells suspended in a simple buffer solution, with various concentrations of lactate or formate added as electron donor, reduced 5 or 10 μM Cr(VI) to Cr(III) over days to weeks. In all nine batch experiments, ⁵³Cr/⁵²Cr ratios of the unreacted Cr(VI) increased as reduction proceeded. In eight experiments covering a range of added donor concentrations up to 100 μM, isotopic fractionation factors were nearly invariant, ranging from 1.0040 to 1.0045, with a mean value somewhat larger than that previously reported for abiotic Cr(VI) reduction (1.0034). One experiment containing much greater donor concentration (10 mM lactate) reduced Cr(VI) much faster and exhibited a lesser fractionation factor (1.0018). These results indicate that ⁵³Cr/⁵²Cr measurements should be effective as indicators of Cr(VI) reduction, either bacterial or abiotic. However, variability in the fractionation factor is poorly constrained and should be studied for a variety of microbial and abiotic reduction pathways.     

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.     

Uranophane dissolution and growth in CaCl2-SiO2(aq) test solutions

The dissolution and growth of uranophane [Ca(UO₂)₂(SiO₃OH)₂·5H₂O] have been examined in Ca- and Si-rich test solutions at low temperatures (20.5 ± 2.0 °C) and near-neutral pH (∼6.0). Uranium-bearing experimental solutions undersaturated and supersaturated with uranophane were prepared in matrices of ∼10⁻² M CaCl₂ and ∼10⁻³ M SiO₂(aq). The experimental solutions were reacted with synthetic uranophane and analyzed periodically over 10 weeks. Interpretation of the aqueous solution data permitted extraction of a solubility constant for the uranophane dissolution reaction and standard state Gibbs free energy of formation for uranophane ( kJ mol⁻¹).     

Calcite precipitation instability under laminar, open-channel flow

We present a 2D numerical model for the growth of calcite from supersaturated aqueous solutions under laminar, open-channel flow conditions. The model couples solution chemistry, precipitation at solution/calcite interfaces, hydrodynamics, diffusion and degassing. The model output is compared with experimental results obtained using an oversaturated calcite solution produced by mixing CaCl₂ and Na₂CO₃. The precipitation rate is observed to increase when the supersaturated solution flows over an obstruction, leading to a growth instability that causes the formation of terraces. At relatively high flow rates, the most important mechanism for this behaviour seems to be hydrodynamic advection of dissolved species either towards or away from the calcite surface, depending on location relative to the obstruction, which deforms the concentration gradients. At lower flow rates, steepening of diffusion gradients around protrusions becomes important. Enhanced degassing over the obstruction due to shallowing and pressure drop is not important on small scales. Diffusion controlled transport close to the calcite surface can lead to a fingering-type growth instability, which generates porous textures. Our results are consistent with existing diffusive boundary layer theory, but for flow over non-smooth surfaces, simple calcite precipitation models that include empirical correlations between fluid flow rate and calcite precipitation rate are inaccurate.     

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.     

The posthumous tectonics and mechanism of exhumation of granitic plutons: The case of the Transbaikal region and the Tien Shan

The study of granitic plutons of the Baikal Highland and the Tien Shan has made it possible to establish new features of their posthumous (after incorporation into the consolidated Earth’s crust) structural reworking and to understand the implications of the cataclastic flow for the exhumation of the crystalline basement in the studied regions. It is shown that granitic plutons undergo appreciable structural transformation at the stages of tectonic reactivation that is significantly separated in time from the moment of formation of plutons as geological bodies. The 3D cataclastic deformation is the main mode of structural reworking of granitic plutons, while the cataclastic flow is the main form of their mobility. Newly recognized slice structures characterize the volumetric deformation of granites.