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.     

The hydrochemistry of a semi-arid pan basin case study: Sua Pan,Makgadikgadi, Botswana

This study presents results on the fluid and salt chemistry for the Makgadikgadi, a substantial continental basin in the semi-arid Kalahari. The aims of the study are to improve understanding of the hydrology of such a system and to identify the sources of the solutes and the controls on their cycling within pans. Sampling took place against the backdrop of unusually severe flooding as well as significant anthropogenic extraction of subsurface brines. This paper examines in particular the relationship between the chemistry of soil leachates, fresh stream water, salty lake water, surface salts and subsurface brines at Sua Pan, Botswana with the aim of improving the understanding of the system’s hydrology. Occasionally during the short wet season (December–March) surface water enters the saline environment and precipitates mostly calcite and halite, as well as dolomite and traces of other salts associated with the desiccation of the lake. The hypersaline subsurface brine (up to TDS 190,000 mg/L) is homogenous with minor variations due to pumping by BotAsh mine (Botswana Ash (Pty) Ltd.), which extracts 2400 m³ of brine/h from a depth of 38 m. Notable is the decrease in TDS as the pumping rate increases which may be indicative of subsurface recharge by less saline water. Isotope chemistry for Sr (⁸⁷Sr/⁸⁶Sr average 0.722087) and S (δ³⁴S average 34.35) suggests subsurface brines have been subject to a lithological contribution of undetermined origin. Recharge of the subsurface brine from surface water including the Nata River appears to be negligible.     

National-scale estimation of potentially harmful element ambient background concentrations in topsoil using parent material classified soil:stream–sediment relationships

Regulatory authorities require estimates of ambient background concentrations (ABCs) of potentially harmful elements (PHEs) in topsoil; such data are currently not available in many countries. High resolution soil geochemical data exist for only part of England and Wales, whilst stream sediment data cover the entire landscape. A novel methodology is presented for estimating soil equivalent ABCs for PHEs from high-resolution (HR) stream sediment geochemical data grouped by common parent materials (PM), using arsenic (As) as an example. Geometric mean (GM) values for local PM groups are used to investigate different approaches for transforming sediment to soil equivalent concentrations. Holdout validation is used to assess: (i) the optimum number of samples for calculating local GM values, and (ii) the optimum scale at which to group data when using linear regression analysis to estimate GM soil ABCs from local sediment geochemical values. Holdout validation showed that the smallest differences were generally observed when five observations were used to calculate the GM and that these should be grouped over the smallest possible area in order to encompass soils over PMs with elevated GM As concentrations. Geometric mean ABCs are estimated and mapped for As in mineral soil across all of England and Wales within delineations of PM polygons. Errors for the estimation of soil equivalent GM As ABCs based on sediment data for an independent validation set were of a similar magnitude to those from holdout validation applied to the original data suggesting the approach is robust. The estimates of soil equivalent ABCs suggest that As exceeds the regulatory threshold used in risk assessments for residential land use (20 mg kg⁻¹) across 16% of the landscape of England and Wales. The applicability of the method for cognate landscapes, and potential refinements is discussed.