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.     

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.     

Reaction mechanism of uranyl in the presence of zero-valent iron nanoparticles

The current study provides an investigation of abiotic reduction of an oversaturated uranyl solution driven by iron nanoparticle oxidation. The reactivity of nano-scale zero-valent iron (ZVI) under mildly oxic conditions (1.2% O₂ and 0.0017% CO₂) was studied in 1000 ppm uranyl solution in the pH range 3-7, at reaction times from 10 min to 4 h. Reductive precipitation of UO₂ was observed as the main process responsible for the removal of uranium from solution with the kinetics of reaction becoming increasingly favourable at higher pH. Despite working with an oversaturated uranium solution, the precipitation of UO₂ occurred in preference to precipitation of UO₃·2H₂O (metaschoepite) at reaction times between 1 and 4 h and for uranyl solutions initially set up at pH ?5. Characterisation of both solid and solution phases was performed using X-ray photoelectron spectroscopy (XPS), focused ion beam (FIB) imaging, X-ray diffraction (XRD) and inductively coupled plasma atomic emission spectroscopy (ICP-AES).     

Anaerobic oxidation of methane (AOM) in marine sediments from the Skagerrak (Denmark): II. Reaction-transport modeling

A steady-state reaction-transport model is applied to sediments retrieved by gravity core from two stations (S10 and S13) in the Skagerrak to determine the main kinetic and thermodynamic controls on anaerobic oxidation of methane (AOM). The model considers an extended biomass-implicit reaction network for organic carbon degradation, which includes extracellular hydrolysis of macromolecular organic matter, fermentation, sulfate reduction, methanogenesis, AOM, acetogenesis and acetotrophy. Catabolic reaction rates are determined using a modified Monod rate expression that explicitly accounts for limitation by the in situ catabolic energy yields. The fraction of total sulfate reduction due to AOM in the sulfate-methane transition zone (SMTZ) at each site is calculated. The model provides an explanation for the methane tailing phenomenon which is observed here and in other marine sediments, whereby methane diffuses up from the SMTZ to the top of the core without being consumed. The tailing is due to bioenergetic limitation of AOM in the sulfate reduction zone, because the methane concentration is too low to engender favorable thermodynamic drive. AOM is also bioenergetically inhibited below the SMTZ at both sites because of high hydrogen concentrations (∼3-6 nM). The model results imply there is no straightforward relationship between pore water concentrations and the minimum catabolic energy needed to support life because of the highly coupled nature of the reaction network. Best model fits are obtained with a minimum energy for AOM of ∼11 kJ mol⁻¹, which is within the range reported in the literature for anaerobic processes.     

Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from southern Africa

Fourier transform infrared spectrometry (FTIR) analyses of olivines from peridotite xenoliths found in southern African kimberlites indicate 0 to 80 ppm H₂O concentrations. OH absorbance profiles across olivine grains show homogeneous H contents from core to edge for most samples. In one sample the olivines are H-free, while another has olivines characterized by lower H contents at the grain edges compared to the cores, indicating H loss during transport of the xenolith to the surface. Flat or near-flat H profiles place severe constraints on the duration of H loss from olivine grains, with implications for kimberlite magma ascent rates. Diffusion equations were used to estimate times of H loss of about 4 h for the sample with heterogeneous olivine H contents. Resulting kimberlite ascent rates are calculated to be 5-37 m s⁻¹ minimum, although these estimates are highly dependent on volatile contents and degassing behavior of the host kimberlite magma. Xenolithic olivines from alkali basalts generally have lower H contents and more pronounced H diffusion profiles than do those from kimberlites. This difference is likely caused by higher magma temperatures and lower ascent rates of alkali basalts compared to kimberlites.     

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.     

Phthalic acid complexation and the dissolution of forsteritic glass studied via in situ FTIR and X-ray scattering

Multiple Internal Reflection Fourier Transform Infra-Red (MIR-FTIR) spectroscopy was developed and used for in situ flow-through experiments designed to study the process of organic acid promoted silicate dissolution. In tandem with the FTIR analysis, ex situ X-ray scattering was used to perform detailed analyses of the changes in the surface structure and chemistry resulting from the dissolution process. Phthalic acid and forsteritic glass that had been Chemically Vapour Deposited (CVD) onto an internal reflection element were used as reactants, and the MIR-FTIR results showed that phthalic acid may promote dissolution by directly binding to exposed Mg metal ion centers on the solid surface. Integrated infrared absorption intensity as a function of time shows that phthalic acid attachment apparently follows a t^1/2 dependence, indicating that attachment is a diffusive process. The diffusion coefficient of phthalic acid was estimated to be approximately 7 × 10⁻⁶ cm² s⁻¹ in the solution near the interface with the glass. Shifts in the infrared absorption structure of the phthalate complexed with the surface compared to the solute species indicate that phthalate forms a seven-membered ring chelate complex. This bidentate complex efficiently depletes Mg from the glass surface, such that after reaction as much as 95% of the Mg may be removed. Surface depletion in Mg causes adsorbate density to fall after an initial attachment stage for the organic ligand. In addition, the infrared analysis shows that silica in the near surface polymerizes after Mg removal, presumably to maintain charge balance. X-ray reflectivity shows that the dissolution rate of forsteritic glass at pH 4 based on Mg removal in such flow-through experiments was equal to 4 × 10⁻¹² mol cm⁻² s⁻¹ (geometric surface area normalized). Reflectivity also shows how the surface mass density decreases during reaction from 2.64 g cm⁻³ to 2.2 g cm⁻³, consistent with preferential loss of Mg from the surface. Auxiliary batch experiments with forsteritic glass films deposited onto soda glass were also completed to add further constraints to the mechanism of reaction. By combining reflectivity with diffuse scatter measurements it is shown that the primary interface changes little in terms of atomic-scale roughness even after removal of several hundred angstroms of material. These measurements unequivocally show how a dicarboxylic acid bonds to and may chelate the dissolution of a magnesium-bearing silicate. At the molecular level the solid surface retreat may best be described by a depinning model where Mg is preferentially removed and residual silica tetrahedra polymerize and act to episodically “pin” the surface.     

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.     

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.     

Integrating diverse methods to understand climate-land interactions in East Africa

The questions of how land use change affects climate, and how climate change affects land use, require examination of societal and environmental systems across space at multiple scales, from the global climate to regional vegetative dynamics to local decision making by farmers and herders. It also requires an analysis of causal linkages and feedback loops between systems. These questions and the conceptual approach of the research design of the Climate-Land Interaction Project (CLIP) are rooted in the classical human-environment research tradition in Geography.This paper discusses a methodological framework to quantify the two-way interactions between land use and regional climate systems, using ongoing work by a team of multi-disciplinary scientists examining climate-land dynamics at multiple scales in East Africa. East Africa is a region that is undergoing rapid land use change, where changes in climate would have serious consequences for people’s livelihoods, and requiring new coping and land use strategies. The research involves exploration of linkages between two important foci of global change research, namely, land use/land cover (LULC) and climate change. These linkages are examined through modeling agricultural systems, land use driving forces and patterns, the physical properties of land cover, and the regional climate. Both qualitative and quantitative methods are being used to illustrate a diverse pluralism in scientific discovery.