Laboratory experiments on the weathering of iron meteorites and carbonaceous chondrites by iron‐oxidizing bacteria

Abstract— Batch culture experiments were performed to investigate the weathering of meteoritic material by iron‐oxidizing bacteria. The aerobic, acidophilic iron oxidizer (A. ferrooxidans) was capable of oxidizing iron from both carbonaceous chondrites (Murchison and Cold Bokkeveld) and iron meteorites (York and Casas Grandes). Preliminary iron isotope results clearly show contrasted iron pathways during oxidation with and without bacteria suggesting that a biological role in meteorite weathering could be distinguished isotopically. Anaerobic iron‐oxidizers growing under pH‐neutral conditions oxidized iron from iron meteorites. These results show that rapid biologically‐mediated alteration of extraterrestrial materials can occur in both aerobic and anaerobic environments. These results also demonstrate that iron can act as a source of energy for microorganisms from both iron and carbonaceous chondrites in aerobic and anaerobic conditions with implications for life on the early Earth and the possible use of microorganisms to extract minerals from asteroidal material.     

Changes in Zn speciation during soil formation from Zn-rich limestones

In order to better understand the long-term speciation and fractionation of Zn in soils, we investigated three soils naturally enriched in Zn (237–864 mg/kg Zn) from the weathering of Zn-rich limestones (40–207 mg/kg Zn) using extended X-ray absorption fine structure (EXAFS) spectroscopy and sequential extractions. The analysis of bulk EXAFS spectra by linear combination fitting (LCF) indicated that Zn in the oolitic limestones was mainly present as Zn-containing calcite (at site Dornach), Zn-containing goethite (Gurnigel) and Zn-containing goethite and sphalerite (Liestal). Correspondingly, extraction of the powdered rocks with 1 M NH₄-acetate at pH 6.0 mobilized only minor fractions of Zn from the Gurnigel and Liestal limestones (<30%), but most Zn from the Dornach rock (81%). In the Dornach soil, part of the Zn released from the dissolving limestone was subsequently incorporated into pedogenic hydroxy-interlayered vermiculite (Zn-HIV, 30%) and Zn-containing kaolinite (30%) and adsorbed or complexed by soil organic and inorganic components (40%). The Gurnigel and Liestal soils contained substantial amounts of Zn-containing goethite (50%) stemming from the parent rock, smaller amounts (20%) of Zn-containing kaolinite (and possibly Zn-HIV), as well as adsorbed or complexed Zn-species (30%). In the soil from Liestal, sphalerite was only found in trace amounts, indicating its dissolution during soil formation. In sequential extractions, large percentages of Zn (55–85%) were extracted in recalcitrant extraction steps, confirming that Zn-HIV, Zn-containing kaolinite and Zn-containing goethite are highly resistant to weathering. These Zn-bearing phases thus represent long-term hosts for Zn in soils over thousands of years. The capability of these phases to immobilize Zn in heavily contaminated soils may however be limited by their uptake capacity (especially HIV and kaolinite) or their abundance in soil.     

Multi-metal contaminant dynamics in temporarily flooded soil under sulfate limitation

In many river basins, floodplain soils have accumulated a variety of metal contaminants, which might be released during periods of flooding. We investigated the dynamics of copper, cadmium, lead, zinc, and nickel in a contaminated freshwater floodplain soil under a realistic sulfate-limited flooding regime in microcosm experiments. We found that most contaminants were initially mobilized by processes driven by the reductive dissolution of Fe(III) and Mn(IV, III) (hydr)oxides. Subsequently, bacterial sulfate respiration resulted in the transformation of the entire available sulfate (2.3 mmol/kg) into chromous reducible sulfur (CRS). Cu K-edge X-ray absorption fine structure (XAFS) spectroscopy revealed that the soil Cu speciation changed from predominantly Cu(II) bound to soil organic matter (SOM) intermittently to 14% metallic Cu(0) and subsequently to 66% copper sulfide (Cu_xS). These Cu_xS precipitates accounted for most of the formed CRS, suggesting that Cu_xS was the dominant sulfide phase formed in the flooded soil. Sequential metal extractions, in agreement with CRS results, suggested that easily mobilizable Cd was completely and Pb partially sequestered in sulfide precipitates, controlling their dissolved concentrations to below detection limits. In contrast, Zn and Ni (as well as Fe) were hardly sequestered into sulfide phases, so that micromolar levels of dissolved Zn and Ni (and millimolar dissolved Fe(II)) persisted in the reduced soil. The finding that Cu, Cd, and Pb were sequestered (but hardly any Zn, Ni, and Fe) is consistent with the thermodynamically predicted sulfide ladder following the increasing solubility products of the respective metal sulfides. The observation that Cd and Pb were sequestered in sulfides despite the presence of remaining SOM-bound Cu(II) suggested that the kinetics of Cu(II) desorption, diffusion, and/or Cu_xS precipitation interfered with the sulfide ladder. We conclude that the dynamics of multiple metal contaminants are intimately coupled under sulfate limitation by the relative thermodynamic stabilities and formation kinetics of the respective metal sulfides.     

Reduced sulfur species in a stratified seawater lake (Rogoznica Lake, Croatia); seasonal variations and argument for organic carriers of reactive sulfur

Over a period of a year, Hg⁰-reactive, total reduced sulfur species (RSS_T), as well as a non-volatile fraction that cannot be gas-stripped at pH ∼2 (RSS_NV), have been measured by voltammetry in a stratified, saline lake. In the hypolimnion, RSS_T is dominated by unusually high (up to 5 mM) dissolved divalent sulfur (S^−II), present as H₂S + HS⁻ and as inorganic polysulfides (H_xS_n^x−2). Less abundant RSS_NV is attributed to dissolved zero-valent sulfur (S⁰) in inorganic polysulfides. Assuming negligible contribution of organic S⁰ species in the hypolimnion, the equilibrium distribution of polysulfide ions is calculated; S₅²⁻ is found to predominate. In the epilimnion, all RSS_T consists of RSS_NV within analytical uncertainty. Through spring and summer, RSS_T and RSS_NV display little vertical or seasonal variation, but they increase dramatically when stratification breaks down in autumn. Based on decay rate, RSS during mixing events is attributed to dissolved S₈ from oxidation of sulfide and decomposition of inorganic polysulfides. This hypothesis quantitatively predicts precipitation of elemental sulfur in a year when colloidal sulfur was observed and predicts no precipitation in a year when it was not observed. Except during mixing events, the entire water column is undersaturated with respect to both rhombic sulfur and biologic sulfur, and the limited variations of RSS exclude hydrophobic and volatile aqueous S₈ as a major species. During such periods, RSS (typically 8 nM) may be associated with organic carbon, perhaps as adsorbed S₈ or as covalently bound polysulfanes or polysulfides. The hypolimnion is viewed as a zero-valent sulfur reactor that creates S⁰-containing, dissolved organic macromolecules during stable stratification periods. Some are sufficiently degradation-resistant and hydrophilic to be dispersed throughout the lake during mixing events, subsequently giving rise to ∼10⁻⁸ M RSS in the oxic water column. Voltammetrically determined RSS in oxic natural waters has often been described as “sulfide” or “metal complexed sulfide”, implying an oxidation state of S^−II; we argue that RSS in oxic Rogoznica Lake waters is mainly S⁰.     

Morphological evolution of dissolving feldspar particles with anisotropic surface kinetics and implications for dissolution rate normalization and grain size dependence: A kinetic modeling study

With previous two-dimensional (2D) simulations based on surface-specific feldspar dissolution succeeding in relating the macroscopic feldspar kinetics to the molecular-scale surface reactions of Si and Al atoms ( [Zhang and Lüttge, 2008] and [Zhang and Lüttge, 2009]), we extended our modeling effort to three-dimensional (3D) feldspar particle dissolution simulations. Bearing on the same theoretical basis, the 3D feldspar particle dissolution simulations have verified the anisotropic surface kinetics observed in the 2D surface-specific simulations. The combined effect of saturation state, pH, and temperature on the surface kinetics anisotropy has been subsequently evaluated, found offering diverse options for morphological evolution of dissolving feldspar nanoparticles with varying grain sizes and starting shapes. Among the three primary faces on the simulated feldspar surface, the (1 0 0) face has the biggest dissolution rate across an extensively wide saturation state range and thus acquires a higher percentage of the surface area upon dissolution. The slowest dissolution occurs to either (0 0 1) or (0 1 0) faces depending on the bond energies of Si-(O)-Si (Φ_Si-O-Si/kT) and Al-(O)-Si (Φ_Al-O-Si/kT). When the ratio of Φ_Si-O-Si/kT to Φ_Al-O-Si/kT changes from 6:3 to 7:5, the dissolution rates of three primary faces change from the trend of (1 0 0) > (0 1 0) > (0 0 1) to the trend of (1 0 0) > (0 0 1) > (0 1 0). The rate difference between faces becomes more distinct and accordingly edge rounding becomes more significant. Feldspar nanoparticles also experience an increasing degree of edge rounding from far-from-equilibrium to close-to-equilibrium. Furthermore, we assessed the connection between the continuous morphological modification and the variation in the bulk dissolution rate during the dissolution of a single feldspar particle. Different normalization treatments equivalent to the commonly used mass, cube assumption, sphere assumption, geometric surface area, and reactive surface area normalizations have been used to normalize the bulk dissolution rate. For each of the treatments, time consistence and grain size dependence of the normalized dissolution rate have been evaluated and the results revealed significant dependences on the magnitude of surface kinetic anisotropy under differing environmental conditions. In general, the normalized dissolution rates are strongly dependent on grain size. Time-consistent normalization treatment varies with the investigated condition. The modeling results suggest that the sphere-, cube-, and BET-normalized dissolution rates are appropriate under the far-from-equilibrium conditions at low pH where these normalizations are time-consistent and are slightly dependent on grain size.     

The structure of monomeric and dimeric uranyl adsorption complexes on gibbsite: A combined DFT and EXAFS study

We investigated the structure of uranyl sorption complexes on gibbsite (pH 5.6-9.7) by two independent methods, density functional theory (DFT) calculations and extended X-ray absorption fine structure (EXAFS) spectroscopy at the U-L_III edge. To model the gibbsite surface with DFT, we tested two Al (hydr)oxide clusters, a dimer and a hexamer. Based on polarization, structure, and relaxation energies during geometry optimization, the hexamer cluster was found to be the more appropriate model. An additional advantage of the hexamer model is that it represents both edges and basal faces of gibbsite. The DFT calculations of (monomeric) uranyl sorption complexes show an energetic preference for the corner-sharing versus the edge-sharing configuration on gibbsite edges. The energy difference is so small, however, that possibly both surface species may coexist. In contrast to the edge sites, sorption to basal sites was energetically not favorable. EXAFS spectroscopy revealed in all investigated samples the same interatomic distances of the uranyl coordination environment (R_U-Oax ≈ 1.80 Å, R_U-Oeq ≈ 2.40 Å), and towards the gibbsite surface (R_U-O ≈ 2.87 Å, R_U-Al ≈ 3.38 Å). In addition, two U-U distances were observed, 3.92 Å at pH 9.7 and 4.30 Å at pH 5.6, both with coordination numbers of ∼1. The short U-U distance is close to that of the aqueous uranyl hydroxo dimer, UO₂(OH)₂, reported as 3.875 Å in the literature, but significantly longer than that of aqueous trimers (3.81-3.82 Å), suggesting sorption of uranyl dimers at alkaline pH. The longer U-U distance (4.30 Å) at acidic pH, however, is not in line with known aqueous uranyl polymer complexes. Based on the EXAFS findings we further refined dimeric surface complexes with DFT. We propose two structural models: in the acidic region, the observed long U-U distance can be explained with a distortion of the uranyl dimer to form both a corner-sharing and an edge-sharing linkage to neighboring Al octahedra, leading to R_U-U = 4.150 Å. In the alkaline region, a corner-sharing uranyl dimer complex is the most favorable. The U-O path at ∼2.87 Å in the EXAFS spectra arises from the oxygen atom linking two Al cations in corner-sharing arrangement. The adsorption structures obtained by DFT calculations are in good agreement with the structural parameters from EXAFS analysis: U-Al (3.394 Å), U-U (3.949 Å), and U-O (2.823 Å) for the alkaline pH model, and U-Al (3.279 Å), U-U (4.150 Å), and U-O (2.743 Å) for the acidic pH model. This work shows that by combining EXAFS and DFT, consistent structural models for uranyl sorption complexes can be obtained, which are relevant to predict the migration behavior of uranium at nuclear facilities.     

Oxygen isotopes composition of sapphires from the French Massif Central: implications for the origin of gem corundum in basaltic fields

Alluvial and colluvial gem sapphires are common in the basaltic fields of the French Massif Central (FMC) but sapphire-bearing xenoliths are very rare, found only in the Menet trachytic cone in Cantal. The O-isotope composition of the sapphires ranges between 4.4 and 13.9‰. Two distinct groups have been defined: the first with a restricted isotopic range between 4.4 and 6.8‰ (n = 22; mean δ¹⁸O = 5.6 ± 0.7‰), falls within the worldwide range defined for blue-green-yellow sapphires related to basaltic gem fields (3.0 < δ¹⁸O < 8.2‰, n = 150), and overlaps the ranges defined for magmatic sapphires in syenite (4.4 < δ¹⁸O < 8.3‰, n = 29). A second group, with an isotopic range between 7.6 and 13.9‰ (n = 9), suggests a metamorphic sapphire source such as biotite schist in gneisses or skarns. The δ¹⁸O values of 4.4–4.5‰ for the blue sapphire-bearing anorthoclasite xenolith from Menet is lower than the δ¹⁸O values obtained for anorthoclase (7.7–7.9‰), but suggest that these sapphires were derived from an igneous reservoir in the subcontinental spinel lherzolitic mantle of the FMC. The presence of inclusions of columbite-group minerals, pyrochlore, Nb-bearing rutile, and thorite in these sapphires provides an additional argument for a magmatic origin. In the FMC lithospheric mantle, felsic melts crystallized to form anorthoclasites, the most evolved peraluminous variant of the alkaline basaltic melt. The O-isotopic compositions of the first group suggests that these sapphires crystallized from felsic magmas under upper mantle conditions. The second group of isotopic values, typified for example by the Le Bras sapphire with a δ¹⁸O of 13.9‰, indicates that metamorphic sapphires from granulites were transported to the surface by basaltic magma.     

Rapid and cyclic aeolian deposition during the Last Glacial in European loess: a high-resolution record from Nussloch, Germany

This paper reports the results of an investigation of the Weichselian Upper Pleniglacial loess sequences of Nussloch (Rhine Valley, Germany) based on stratigraphy, palaeopedology, sedimentology, palynology, malacology and geochemistry (δ¹³C), supported by radiocarbon, TL and OSL dating. Grain-size and magnetic susceptibility records are taken at 5 cm intervals from the Upper Pleniglacial (UPG) loess. The data indicate cyclic variations in loess deposition between ca 34 and 17 ka, when the sedimentation rate is especially high (1.0–1.2 m per ka for more than 10 m). The grain-size index (GSI: ratio of coarse silt versus fine silt and clay) shows variations, which are assumed to be an indirect measurement of wind intensity. The sedimentation rate, interpreted from the profiles, indicates high values in loess (Loess events LE-1 to LE-7) and low or negligible values in tundra gley horizons G1 to G8. OSL ages from the loess and ¹⁴C dates from organic matter in the loess show that loess deposition was rapid but was interrupted by shorter periods of reduced aeolian sedimentation. Comparison between the data from Nussloch and other European sequences demonstrates a progressive coarsening of the loess deposits between ca 30 and 22 ka. This coarsening trend ends with a short but major grain-size decrease and is followed by an increase to a new maximum at 20 ± 2 ka (“W” shape). Correlation between the loess GSI and the Greenland ice-core dust records, suggests a global connection between North Atlantic and Western European global atmospheric circulation and wind regimes. In addition, the typical Upper Pleniglacial loess deposition begins at ca 30–31 ka, close to Heinrich event (HE) 3, and the main period of loess sedimentation at about 25 ± 2 ka is coeval to HE 2. Correlation of magnetic susceptibility and grain-size records shows that the periods, characterised by high GSI, coincide with an increase in the amount of ferromagnetic minerals reworked from the Rhine alluvial plain. They suggest enhancement in the frequency of the storms from N–NW. These results are integrated within a palaeogeographical model of dust transport and deposition in Western Europe for the Weichselian Upper Pleniglacial (or Late Pleniglacial).     

IR calibrations for water determination in olivine,r-GeO<Subscript>2</Subscript>, and SiO<Subscript>2</Subscript> polymorphs

Mineral-specific IR absorption coefficients were calculated for natural and synthetic olivine, SiO₂ polymorphs, and GeO₂ with specific isolated OH point defects using quantitative data from independent techniques such as proton–proton scattering, confocal Raman spectroscopy, and secondary ion mass spectrometry. Moreover, we present a routine to detect OH traces in anisotropic minerals using Raman spectroscopy combined with the “Comparator Technique”. In case of olivine and the SiO₂ system, it turns out that the magnitude of ε for one structure is independent of the type of OH point defect and therewith the peak position (quartz ε = 89,000 ± 15,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}), but it varies as a function of structure (coesite ε = 214,000 ± 14,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}; stishovite ε = 485,000 ± 109,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}). Evaluation of data from this study confirms that not using mineral-specific IR calibrations for the OH quantification in nominally anhydrous minerals leads to inaccurate estimations of OH concentrations, which constitute the basis for modeling the Earth’s deep water cycle.     

Oblateness effect of Saturn on periodic orbits in the Saturn-Titan restricted three-body problem

We explore the effect of oblateness of Saturn (more massive primary) on the periodic orbits and the regions of quasi-periodic motion around both the primaries in the Saturn-Titan system in the framework of planar circular restricted three-body problem. First order interior and exterior mean motion resonances are located. The effect of oblateness is studied on the location, nature and size of periodic and quasi-periodic orbits, using the numerical technique of Poincare surface of sections. Some of the periodic orbits change to quasi-periodic orbits due to the effect of oblateness and vice-versa. The stability of the orbits around Saturn, Titan and both varies with the inclusion of oblateness. The centers of the periodic orbits around Titan move towards Saturn, whereas those around Saturn move towards Titan. For the orbit around Titan at C=2.9992, x=0.959494, the apocenter becomes pericenter. By incorporating oblateness effect, the orbit around Titan at C=2.99345, x=0.924938 is captured by Saturn, remains in various trajectories around Saturn, and as time progresses it spirals away around both the primaries.