Carbon isotopes and petrography of kerogens in ∼3.5-Ga hydrothermal silica dikes in the North Pole area, Western Australia

More than 600 specimens of ∼3.5 Ga-old hydrothermal silica dikes from the North Pole area, Pilbara craton, Western Australia, have been studied petrographically. The kerogens in 44 samples have been analyzed isotopically (C and N) and chemically (C, N, and H). The silica dikes are composed mainly of fine-grained silica (modal abundance: >97%) and are classified into two types by minor mineral assemblages: B(black)-type and G(gray)-type. The B-type silica dikes contain kerogen (0.37 to 6.72 mgC/g; average 2.44 mgC/g, n = 21) and disseminated sulfides, dominantly pyrite and Fe-poor sphalerite. In some cases, carbonate and apatite are also present. Their silica-dominated and sulfide-poor mineral assemblages suggest precipitation from low-temperature reducing hydrothermal fluid (likely 100-200°C). On the other hand, the G-type silica dikes are sulfide-free and concentrations of kerogen are relatively low (0.05 to 0.41 mgC/g; average 0.17 mgC/g, n = 13). They typically contain Fe-oxide (mainly hematite) which commonly replaces cubic pyrite and rhombic carbonate. Some G-types occur along secondary quartz veins. These textures indicate that the G-type silica dikes were formed by postdepositional metasomatism (oxidation) of the B-types, and that the B-types probably possess premetasomatic signatures. The δ¹³C values of kerogen in the B-types are −38.1 to −33.1‰ (average −35.9‰, n = 21), which are ∼4‰ lower than those of the G-types (−34.5 to −30.0‰; average −32.2‰, n = 19), and ∼6‰ lower than bedded chert (−31.2 to −29.4‰; average −30.5‰, n = 4). This indicates the preferential loss of ¹²C during the metasomatism (estimated fractionation factor: 0.9985). Considering the metasomatic effect on carbon isotopes with probably minor diagenetic and metamorphic overprints, we conclude that the original δ¹³C values of the kerogen in the silica dikes would have been heterogeneous (∼5‰) and at least some material had initial δ¹³C values of ≤ −38‰. The inferred ¹³C-depletions of organic carbon could have been produced by anaerobic chemoautotrophs such as methanogen, but not by aerobic photoautotrophs. This is consistent with the estimated physical and chemical condition of the hydrothermal fluid, which was probably habitable for anaerobic and thermophilic/hyperthermophilic chemoautotrophs. Alternatively, the organic matter may have been possibly produced by abiological reaction such as Fischer-Tropsch Type (FTT) synthesis under the hydrothermal condition. However, the estimated condition is inconsistent with the presence of the effective catalysts for the FTT reaction (i.e., Fe-Ni alloy, magnetite, and hematite). These lines of evidence suggest the possible existence of biosphere in the ∼3.5 Ga sub-seafloor hydrothermal system.     

Oxygen isotope fractionation in synthetic magnesian calcite

Mg-bearing calcite was precipitated at 25°C in closed system free-drift experiments from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The chemical and isotope composition of the solution and precipitate were investigated during time course experiments of 24-h duration. Monohydrocalcite and calcite precipitated early in the experiments (<8 h), while Mg-calcite was the predominant precipitate (>95%) thereafter. Solid collected at the end of the experiments displayed compositional zoning from pure calcite in crystal cores to up to 23 mol% MgCO₃ in the rims. Smaller excursions in Mg were superimposed on this chemical record, which is characteristic of oscillatory zoning observed in synthetic and natural solid-solution carbonates of differing solubility. Magnesium also altered the predominant morphology of crystals over time from the {104} to {100} and {110} growth forms.The oxygen isotope fractionation factor for the magnesian-calcite-water system (as 10³lnα_Mg-cl-H2O) displayed a strong dependence on the mol% MgCO₃ in the solid phase, but quantification of the relationship was difficult due to the heterogeneous nature of the precipitate. Considering only the Mg-content and δ¹⁸O values for the bulk solid, 10³lnα_Mg-cl-H2O increased at a rate of 0.17 ± 0.02 per mol% MgCO₃; this value is a factor of three higher than the single previous estimate (Tarutani T., Clayton R.N., and Mayeda T. K. (1969) The effect of polymorphims and magnesium substitution on oxygen isotope fractionation between calcium carbonate and water. Geochim. Cosmochim. Acta 33, 987-996). Nevertheless, extrapolation of our relationship to the pure calcite end member yielded a value of 27.9 ± 0.02, which is similar in magnitude to published values for the calcite-water system. Although no kinetic effect was observed on 10³lnα_Mg-cl-H2O for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol · m⁻² · h⁻¹, it was impossible to disentangle the potential effect(s) of precipitation rate and Mg-content on 10³lnα_Mg-cl-H2O due to the heterogeneous nature of the solid.The results of this study suggest that paleotemperatures inferred from the δ¹⁸O values of high magnesian calcite (>10 mol% MgCO₃) may be significantly underestimated. Also, the results underscore the need for additional experiments to accurately characterize the effect of Mg coprecipitation on the isotope systematics of calcite from a chemically homogeneous precipitate or a heterogeneous material that is analyzed at the scale of chemical and isotopic zonation.     

Petrology and shock metamorphism of the olivine-phyric shergottite Yamato 980459: Evidence for a two-stage cooling and a single-stage ejection history

The basaltic Martian meteorite Yamato 980459 consists of large olivine phenocrysts and often prismatic pyroxenes set into a fine-grained groundmass of smaller more Fe-rich olivine, chromite, and an interstitial residual material displaying quenching textures of dendritic olivine, chain-like augite and sulfide droplets in a glassy matrix. Yamato 980459 is, thus, the only Martian meteorite without plagioclase/maskelynite. Olivine is compositionally zoned from a Mg-rich core to a Fe-rich rim with the outer few micrometers being especially rich in iron. With Fo₈₄ the cores are the most magnesian olivines found in Martian meteorites so far. Pyroxenes are also mostly composite crystals of large orthopyroxene cores and thin Ca-rich overgrowths. Separate pigeonite and augites are rare. On basis of the mineral compositions, the cooling rates determined from crystal morphologies, and crystal grain size distributions it is deduced that the parent magma of Yamato 980459 initially cooled under near equilibrium conditions e.g., in a magma chamber allowing chromite and the Mg-rich silicates to form as cumulus phases. Fractional crystallization at higher cooling rates and a low degree of undercooling let to the formation of the Ca-, Al-, and Fe-rich overgrowths on olivine and orthopyroxene while the magma was ascending towards the Martian surface. Finally and before plagioclase and also phosphates could precipitate, the magma was very quickly erupted quenching the remaining melt to glass, dendritic silicates and sulfide droplets. The shape preferred orientation of olivine and pyroxene suggests a quick, thin outflow of lava. According to the shock effects found in the minerals of Yamato 980459, the meteorite experienced an equilibration shock pressure of about 20-25 GPa. Its near surface position allowed the ejection from the planet’s surface already by a single impact event and at relatively low shock pressures.     

European vegetation during Marine Oxygen Isotope Stage-3

European vegetation during representative “warm” and “cold” intervals of stage-3 was inferred from pollen analytical data. The inferred vegetation differs in character and spatial pattern from that of both fully glacial and fully interglacial conditions and exhibits contrasts between warm and cold intervals, consistent with other evidence for stage-3 palaeoenvironmental fluctuations. European vegetation thus appears to have been an integral component of millennial environmental fluctuations during stage-3; vegetation responded to this scale of environmental change and through feedback mechanisms may have had effects upon the environment. The pollen-inferred vegetation was compared with vegetation simulated using the BIOME 3.5 vegetation model for climatic conditions simulated using a regional climate model (RegCM2) nested within a coupled global climate and vegetation model (GENESIS-BIOME). Despite some discrepancies in detail, both approaches capture the principal features of the present vegetation of Europe. The simulated vegetation for stage-3 differs markedly from that inferred from pollen analytical data, implying substantial discrepancy between the simulated climate and that actually prevailing. Sensitivity analyses indicate that the simulated climate is too warm and probably has too short a winter season. These discrepancies may reflect incorrect specification of sea surface temperature or sea-ice conditions and may be exacerbated by vegetation-climate feedback in the coupled global model.     

Isotopic exchange of carbon-bound hydrogen over geologic timescales

The increasing popularity of compound-specific hydrogen isotope (D/H) analyses for investigating sedimentary organic matter raises numerous questions about the exchange of carbon-bound hydrogen over geologic timescales. Important questions include the rates of isotopic exchange, methods for diagnosing exchange in ancient samples, and the isotopic consequences of that exchange. This article provides a review of relevant literature data along with new data from several pilot studies to investigate such issues. Published experimental estimates of exchange rates between organic hydrogen and water indicate that at warm temperatures (50-100°C) exchange likely occurs on timescales of 10⁴ to 10⁸ yr. Incubation experiments using organic compounds and D-enriched water, combined with compound-specific D/H analyses, provide a new and highly sensitive method for measuring exchange at low temperatures. Comparison of δD values for isoprenoid and n-alkyl carbon skeletons in sedimentary organic matter provides no evidence for exchange in young (<1 Ma), cool sediments, but strong evidence for exchange in ancient (>350 Ma) rocks. Specific rates of exchange are probably influenced by the nature and abundance of organic matter, pore-water chemistry, the presence of catalytic mineral surfaces, and perhaps even enzymatic activity.Estimates of equilibrium fractionation factors between organic H and water indicate that typical lipids will be depleted in D relative to water by ∼75 to 140‰ at equilibrium (30°C). Thus large differences in δD between organic molecules and water cannot be unambiguously interpreted as evidence against hydrogen exchange. A better approach may be to use changes in stereochemistry as a proxy for hydrogen exchange. For example, estimated rates of H exchange in pristane are similar to predicted rates for stereochemical inversion in steranes and hopanes. The isotopic consequences of this exchange remain in question. Incubations of cholestene with D₂O indicate that the number of D atoms incorporated during structural rearrangements can be far less than the number of C-H bonds that are broken. Sample calculations indicate that, for steranes in immature sediments, the D/H ratio imparted by biosynthesis may be largely preserved in spite of significant structural changes.     

Pyrite oxidation in moist air

The rate of pyrite oxidation in moist air was determined by measuring, over time, the pressure difference between a sealed chamber containing pyrite plus oxygen and a control. The experiments carried out at 25°C, 96.7% fixed relative humidity, and oxygen partial pressures of 0.21, 0.61, and 1.00 atm showed that the rate of oxygen consumption is a function of oxygen partial pressure and time. The rates of oxygen consumption (r, mol/m²sec) fit the expression

(A)     

Oxygen-exchange pathways in aluminum polyoxocations

Using molecular dynamics simulations and electronic structure methods, we postulate a mechanism to explain the complicated reactivity trends that are observed for oxygen isotope exchange reactions between sites in aluminum polyoxocations of the ε-Keggin type and bulk solution. Experimentally, the molecules have four nonequivalent oxygens that differ considerably in reactivity both within a molecule, and between molecules in the series: Al₁₃, GaAl₁₂, and GeAl₁₂ [MO₄Al₁₂(OH)₂₄(H₂O)₁₂ⁿ⁺(aq); with M = Al(III) for Al₁₃, n = 7; M = Ga(III) for GaAl₁₂, n = 7; M = Ge(IV) for GeAl₁₂, n = 8]. We find that a partly dissociated, metastable intermediate molecule of expanded volume is necessary for exchange of both sets of μ₂-OH and that the steady-state concentration of this intermediate reflects the bond strengths between the central metal and the μ₄-O. Thus the central metal exerts extraordinary control over reactions at hydroxyl bridges, although these are three bonds away.This mechanism not only explains the reactivity trends for oxygen isotope exchange in μ₂-OH and η-OH₂ sites in the ε-Keggin aluminum molecules, but also explains the observation that the reactivities of minerals tend to reflect the presence of highly coordinated oxygens, such as the μ₄-O in boehmite, α-, and γ-Al₂O₃ and their Fe(III) analogs. The partial dissociation of these highly coordinated oxygens, coupled with simultaneous activation and displacement of neighboring metal centers, may be a fundamental process by which metals atoms undergo ligand exchanges at mineral surfaces.     

Pt-Re-Os systematics of group IIAB and IIIAB iron meteorites

The Pt-Re-Os isotopic and elemental systematics of 13 group IIAB and 23 group IIIAB iron meteorites are examined. As has been noted previously for iron meteorite groups and experimental systems, solid metal-liquid metal bulk distribution coefficients (D values) for both IIAB and IIIAB systems show D_Os>D_Re>>D_Pt>1 during the initial stages of core crystallization. Assuming closed-system crystallization, the latter stages of crystallization for each core are generally characterized by D_Pt>D_Re>D_Os. The processes governing the concentrations of these elements are much more complex in the IIIAB core relative to the IIAB core. Several crystallization models utilizing different starting parameters and bulk distribution coefficients are considered for the Re-Os pair. Each model has flaws, but in general, the results suggest that the concentrations of these elements were dominated by equilibrium crystallization and subsequent interactions between solid metal and both equilibrium and evolved melts. Late additions of primitive metal to either core were likely minor or nonexistent.The ¹⁸⁷Re-¹⁸⁷Os systematics of the IIAB and IIIAB groups are consistent with generally closed-system behavior for both elements since the first several tens of Ma of the formation of the solar system, consistent with short-lived chronometers. The Re-Os isochron ages for the complete suites of IIAB and IIIAB irons are 4530 ± 50 Ma and 4517 ± 32 Ma, respectively, and are similar to previously reported Re-Os ages for the lower-Ni endmembers of these two groups. Both isochrons are consistent with, but do not require crystallization of the entire groups within 10-30 Ma of the initiation of crystallization.The first high-precision ¹⁹⁰Pt-¹⁸⁶Os isochrons for IIAB and IIIAB irons are presented. The Pt-Os isochron ages for the IIAB and IIIAB irons, calculated using the current best estimate of the λ for ¹⁹⁰Pt, are 4323 ± 80 Ma and 4325 ± 26 Ma respectively. The Re-Os and Pt-Os ages do not overlap within the uncertainties. The younger apparent ages recorded by the Pt-Os system likely reflect error in the ¹⁹⁰Pt decay constant. The slope from the Pt-Os isochron is combined with the age from the Re-Os isochron for the IIIAB irons to calculate a revised λ of 1.415 × 10⁻¹² a⁻¹ for ¹⁹⁰Pt, although additional study of this decay constant is still needed.     

Formation and stability of schwertmannite in acidic mining lakes

Schwertmannite (ideal formula: Fe₈O₈(OH)₆SO₄) is typically found as a secondary iron mineral in pyrite oxidizing environments. In this study, geochemical constraints upon its formation are established and its role in the geochemical cycling of iron between reducing and oxidizing conditions are discussed. The composition of surface waters was analyzed and sediments characterized by X-ray diffraction, FTIR spectroscopy and determination of the Fe:S ratio in the oxalate extractable fraction from 18 acidic mining lakes. The lakes are exposed to a permanent supply of pyritegenous ferrous iron from adjacent ground water. In 3 of the lakes the suspended matter was fractionated using ultra filtration and analyzed with respect to their mineral composition. In addition, stability experiments with synthetic schwertmannite were performed. The examined lake surface waters were O₂-saturated and have sulfate concentrations (10.3 ± 5.5 mM) and pH values (3.0 ± 0.6) that are characteristic for the stability window of schwertmannite. Geochemical modeling implied that i) the waters were saturated with respect to schwertmannite, which controlled the activity of Fe³⁺ and sulfate, and ii) a redox equilibrium exists between Fe²⁺ and schwertmannite. In the uppermost sediment layers (1 to 5 cm depth), schwertmannite was detectable in 16 lakes—in 5 of them by all three methods. FTIR spectroscopy also proved its occurrence in the colloidal fraction (1-10 kDa) in all of the 3 investigated lake surface waters. The stability of synthetic schwertmannite was examined as a function of pH (2-7) by a 1-yr experiment. The transformation rate into goethite increased with increasing pH. Our study suggests that schwertmannite is the first mineral formed after oxidation and hydrolysis of a slightly acidic (pH 5-6), Fe(II)-SO₄ solution, a process that directly affects the pH of the receiving water. Its occurrence is transient and restricted to environments, such as acidic mining lakes, where the coordination chemistry of Fe³⁺ is controlled by the competition between sulfate and hydroxy ions (i.e. mildly acidic).     

Carbon isotope effects in carbonate systems

Global carbon cycle models require a complete understanding of the δ¹³C variability of the Earth’s C reservoirs as well as the C isotope effects in the transfer of the element among them. An assessment of δ¹³C changes during CO₂ loss from degassing magmas requires knowledge of the melt-CO₂ carbon isotope fractionation. In order to examine the potential size of this effect for silicate melts of varying composition, ¹³C reduced partition functions were computed in the temperature range 275 to 4000 K for carbonates of varying bond strengths (Mg, Fe, Mn, Sr, Ba, Pb, Zn, Cd, Li, and Na) and the polymorphs of calcite. For a given cation and a given pressure the ¹³C content increases with the density of the carbonate structure. For a given structure the tendency to concentrate ¹³C increases with pressure. The effect of pressure (‰/10 kbar) on the size of the reduced partition function of aragonite varies with temperature; in the pressure range 1 to 10⁵ bars the change is given by:

(1)     

Experimental study of europium (III) coprecipitation with calcite

Partitioning of Eu(III) in calcite, CaCO₃, was evaluated with the aim of collecting data on partition coefficients and to enhance understanding of the incorporation mechanisms. This information will aid in the interpretation of geological processes from rare Earth element (REE) data and in the use of Eu(III) as a chemical analogue for the trivalent actinides, particularly Am(III) and Cm(III). Coprecipitation experiments were carried out by the constant addition method at 25°C and P_CO2 = 1 atm. Eu(III) was strongly partitioned from the solution into calcite. For dilute solid solutions (X_Eu < 0.001), Eu partition coefficients were estimated to be 770 ± 290 and found to be independent of calcite precipitation rate in the range of 0.02 to 2.7 nmol mg⁻¹ min⁻¹. This could be explained by the approximately equal values of the Eu partition and adsorption coefficients. Several solid solution models were tested. A vacancy model for Eu₂(CO₃)₃-CaCO₃ is consistent with the experimental results and constraints on geometry for Eu fit in the calcite lattice. For low Eu content, vacancy density is independent of Eu concentration in the solid so logarithm of the ion activity product, log (Eu)²(CO₃²⁻)³, depends linearly on log X_Eu². The fit of the data to such a model is good evidence that Eu(III) is taken up as a true solid solution, not simply by physical trapping. A model using EuOHCO₃-CaCO₃ is also consistent with the uptake stoichiometry, but EuOH²⁺ substitution for Ca²⁺ would be expected to distort the calcite structure more than is compatible with such a high K_D. Several other models, including EuNa(CO₃)₂-CaCO₃, were abandoned because their stoichiometric relationships did not fit the experimental data.     

Reversibility of radiocaesium sorption on illite

Adsorption of trace amounts of radiocaesium on NH₄-, K-, and Na-saturated Fithian illite and subsequent desorption by 1 M NH₄ showed that a substantial amount of radiocaesium (44%, 46%, and 91% for NH₄-, K-, and Na-illite, respectively) cannot be desorbed after only 5 min of adsorption. Our results suggest that this instantaneous fixation is caused by the collapse of the frayed edges of the clay mineral and the relatively high concentration of radiocaesium building up in solution in the batch desorption experiments. Consequently, commonly applied high-NH₄ extractions underestimate truly exchangeable amounts of radiocaesium in soils and sediments containing illitic clay minerals. The rate of desorption of trace amounts of radiocaesium from the solids using high NH₄ or Cs concentrations has a half-life of about 2 yr, reflecting radiocaesium desorption from (partially) collapsed interlayers. Extraction of radiocaesium from illite after 5 min of contact time with a Cs-selective adsorbent or a 1 × 10^-6 M CsCl solution shows that 100% of the bound radiocaesium is readily available. The desorption rate in the presence of a Cs-selective adsorbent has a half-life of about 0.2 yr. Desorption of radiocaesium from illite using different ammonium concentrations shows that radiocaesium partitioning follows reversible ion-exchange theory if the NH₄ concentration is below 1 × 10^-4 M, and sufficient time (weeks) is allowed for the reaction to proceed. Thus, radiocaesium sorption reversibility in the natural environment is much higher than generally assumed, and equilibrium solid/liquid partitioning may be assumed for the long-term modelling of radiocaesium mobility in the natural environment. In the particular case of anoxic freshwater sediments with very high NH₄ concentrations in the pore waters (up to several mmol.L^-1), collapse of the frayed edges of illite may occur, influencing radiocaesium partitioning. If collapse occurs before radiocaesium adsorbs to illite, high caesium sorption reversibility as measured by high-NH₄ extraction can be expected because further collapse of the frayed edges during the extraction procedure will be limited. This effect has indeed been observed earlier in the extraction of radiocaesium from anoxic freshwater sediments with high-NH₄ solutions and was as yet unexplained.     

Manganese-chromium isotope systematics of enstatite meteorites

We present results of a study of the ⁵³Mn-⁵³Cr isotope systematics in the enstatite chondrites and achondrites (aubrites). The goal of this study was to explore the capabilities of this isotope system to obtain chronological information on these important classes of meteorites and to investigate the original distribution in the inner solar system of the short-lived radionuclide ⁵³Mn. Our earlier work (Lugmair and Shukolyukov, 1998; Shukolyukov and Lugmair, 2000a) has shown that the asteroid belt bodies are characterized by essentially the same initial ⁵³Mn abundance. However, we have found the presence of a gradient in the abundance of the radiogenic ⁵³Cr between the earth-moon system, Mars, and the asteroid Vesta. If this gradient is considered as a function of the heliocentric distance a linear radial dependence is indicated. This can be explained either by an early, volatility controlled Mn/Cr fractionation in the nebula or by an original radially heterogeneous distribution of ⁵³Mn. The enstatite chondrites are suggested to form in the inner zones of the solar nebula, much closer to the Sun than the ordinary chondrites. Therefore, their investigation may be an important test on the hypothesis on a radial heterogeneity in the initial ⁵³Mn.We have studied the bulk samples of the EH4-chondrites Indarch and Abee and the EL6-chondrite Khairpur. Although these meteorites have essentially the same Mn/Cr ratio as the ordinary chondrites, the relative abundance of the radiogenic ⁵³Cr is three times smaller than in the ordinary chondrites. Because these meteorites are primitive (undifferentiated) and no Mn/Cr fractionation had occurred within their parent bodies, this difference is a strong argument in favor of an initially heterogeneous distribution of ⁵³Mn in the early inner solar system. This finding is also consistent with formation of the enstatite chondrites in the inner zones of the solar nebula. Using the characteristic ⁵³Cr excess of the enstatite chondrites and the observed gradient, their place of origin falls at about 1.4 AU or somewhat closer to the Sun (i.e. >1.0-1.4 AU).We also present chronological results for the enstatite chondrites and achondrites. The ‘absolute’ ⁵³Mn-⁵³Cr ages of the EH4-chondrites are old: ∼4565 Ma. The EL6-chondrite Khairpur is ∼4.5 Ma younger, which is in good agreement with the ¹²⁹I-¹²⁹Xe data from the literature. The age of the aubrite Peña Blanca Spring appears to be similar to those of the enstatite chondrites while that of the aubrite Bishopville is at least ∼10 Ma younger, which is also in agreement with the ¹²⁹I-¹²⁹Xe data. The results from bulk samples of aubrites indicate that the last Mn/Cr fractionation in their parent body occurred ∼ 4563 Ma ago and imply an evolution of the Mn-Cr isotope system in an environment with an higher than chondritic Mn/Cr ratio for several millions of years.     

Experimental study of the stability of aluminate-borate complexes in hydrothermal solutions

Formation of aqueous aluminate-borate complexes was characterized at 25°C using ²⁷Al NMR spectroscopy, and at 50-200°C via measurements of gibbsite and boehmite solubility in the presence of boric acid. ²⁷Al spectra performed at pH = 9 in Al-B solution with m(B) = 0.02 show the presence of two peaks at 80.5 and 74.5 ppm which correspond to Al(OH)₄⁻ and a single Al-substituted Q¹_Al dimer, Al(OH)₃OB(OH)₂⁻, respectively. In 0.08 m and 0.2 m borate solution, a third peak appears at 68.5 ppm which can be assigned to the Q²_Al trimer Al(OH)₂O₂(B(OH)₂)₂⁻. These chemical shifts are close to those measured for Al(OH)₃OSi(OH)₃⁻ and Al(OH)₂O₂(Si(OH)₃)₂⁻ (74 and 69.5 ppm, respectively; Pokrovski et al., Min. Mag.62a (1998), 1194) which demonstrates the similar structure of Al-B and Al-Si complexes formed in alkaline solutions. Gibbsite and boehmite solubility were measured in weakly basic solutions as a function of boric acid concentration at 50°C and 78 to 200°C, respectively. Equilibrium was reached within several days at m(B) = 0.01-0.1, but more slowly at higher boron concentrations, and at 50°C and m(B) = 0.2, Al concentration increased continuously during at least 3 months as a result of the sluggish formation of Al-polyborates. The equilibrium constant of the reaction Al(OH)₄⁻ + B(OH)₃⁰_(aq) = Al(OH)₃OB(OH)₂⁻ + H₂O decreases very slowly with increasing temperature to 200°C. The log K values are 1.58 ± 0.10, 1.46 ± 0.10, 1.52 ± 0.15, and 1.25 ± 0.15 at 50, 78, 150 and 200°C, respectively, which result in the following values of the standard thermodynamic properties for this reaction: Δ_rG⁰ = −9.22 ± 3.25 kJ/mol, Δ_rH⁰ = −4.6 ± 2.5 kJ/mol, Δ_rS⁰ = 15.5 ± 6.9 J/mol K. The thermodynamic data generated in this study indicate that Al-B complexes can dominate aqueous aluminum speciation in solutions containing ≥0.7 g/L of boron at temperature to at least 400°C.     

Role of molecular oxygen in the dissolution of siderite and rhodochrosite

The dissolution of siderite (FeCO₃) and rhodochrosite (MnCO₃) under oxic and anoxic conditions is investigated at 298 K. The anoxic dissolution rate of siderite is 10^−8.65 mol m⁻² s⁻¹ for 5.5 < pH < 12 and increases as [H⁺]^0.75 for pH < 5.5. The pH dependence is consistent with parallel proton-promoted and water hydrolysis dissolution pathways. Atomic force microscopy (AFM) reveals a change in pit morphology from rhombohedral pits for pH > 4 to pits elongated at one vertex for pH < 4. Under oxic conditions the dissolution rate decreases to below the detection limit of 10⁻¹⁰ mol m⁻² s⁻¹ for 6.0 < pH < 10.3, and hillock precipitation preferential to steps is observed in concurrent AFM micrographs. X-ray photoelectron spectroscopy (XPS) and thermodynamic analysis identify the precipitate as ferrihydrite. At pH > 10.3, the oxic dissolution rate is as high as 10^−7.5 mol m⁻² s⁻¹, which is greater than under the corresponding anoxic conditions. A fast electron transfer reaction between solution O₂ or [Fe³⁺(OH)₄]⁻ species and surficial >Fe^II hydroxyl groups is hypothesized to explain the dissolution kinetics. AFM micrographs do not show precipitation under these conditions. Anoxic dissolution of rhodochrosite is physically observed as rhombohedral pit expansion for 3.7 < pH < 10.3 and is chemically explained by parallel proton- and water-promoted pathways. The dissolution rate law is 10^−4.93[H⁺] + 10^−8.45 mol m⁻² s⁻¹. For 5.8 < pH < 7.7 under oxic conditions, the AFM micrographs show a tabular precipitate growing by preferential expansion along the a-axis, though the macroscopic dissolution rate is apparently unaffected. For pH > 7.7 under oxic conditions, the dissolution rate decreases from 10^−8.45 to 10^−9.0 mol m⁻² s⁻¹. Flattened hillock precipitates grow across the entire surface without apparent morphological influence by the underlying rhodochrosite surface. XPS spectra and thermodynamic calculations implicate the precipitate as bixbyite for 5.8 < pH < 7.7 and MnOOH (possibly feitnkechtite) for pH >7.7.     

Source enrichment processes responsible for isotopic anomalies in oceanic island basalts

A new method has been developed to separate the compositional variations in ocean island basalts into those that result from variations in source composition and from the melting process itself. The approach depends on correlations between isotope ratios, which can only come from source inhomogeneities, and elemental concentrations. Analysis of three data sets shows that the inhomogeneities beneath Theistareykir, in NE Iceland, Kilauea and Pitcairn can be produced by subduction of oceanic islands and volcanic ridges. The thicknesses of the lithosphere on which such islands were constructed and potential temperatures of the plumes that produced them can be estimated from the geochemical observations. Model ages are harder to determine, though simple assumptions give about 400 Ma for the Theistareykir source and 1.2 Ga for Kilauea. The model may also provide a physical explanation for the commonly used isotopic classification of ocean island basalts, with the isotopic composition changing from HIMU through EMII to EMI as the melt fraction increases. These results have been obtained from a small number of data sets obtained from ocean island basalts erupted in small areas during short time intervals. More such observations are needed to discover whether geochemical observations from other islands are consistent with the same model.     

Modeling diagenesis of lead in sediments of a Canadian Shield lake

Triplicate porewater lead concentration profiles were determined on six occasions in a Canadian Shield lake. Total Pb concentrations were also measured in a dated core obtained at the same site. This information, as well as an extensive dataset comprising ancillary geochemical measurements on porewaters and sediment and the population densities of benthic animals, is used in a one-dimensional transport-reaction diagenetic model to investigate the transport and mobilization of Pb in these sediments. Application of the model consistently indicates the presence of a zone of Pb production to the porewaters that lies above a zone of Pb consumption. The profiles of various porewater constituents and thermodynamic calculations indicate that Pb is mobilized in the zone of production by the reductive dissolution of iron oxyhydroxides, whereas it is removed in the zone of consumption by precipitation as a solid sulfide. Rate constants are estimated for reductive iron dissolution (k_d^Fe(III) = 2.0 ± 0.5 × 10⁻¹ cm³ mol⁻¹ s⁻¹), Pb adsorption on iron oxyhydroxides (k_ads^Pb = 98 ± 55 cm³ mol⁻¹ s⁻¹), and Pb precipitation (k_ppt^Pb = 8 × 10⁻²⁰ mol cm⁻³ s⁻¹ to 16 ± 13 × 10⁻²² mol cm⁻³ s⁻¹, depending on the solubility product assumed for the precipitation of PbS). According to model calculations, diagenetic processes, such as remobilization, molecular diffusion, bioturbation, and bioirrigation have a negligible influence on the solid phase Pb profile. In agreement with this finding, the present-day fluxes of dissolved Pb by diffusion (J_D^Pb = −6.5 × 10⁻¹¹ mol cm⁻² yr⁻¹), bioturbation (J_B^Pb = −1.1 × 10⁻¹³ mol cm⁻² yr⁻¹), and bioirrigation (J_I^Pb = −1.5 × 10⁻¹¹ mol cm⁻² yr⁻¹) are small compared to the flux of Pb deposited with settling particles (J_S^Pb = 5.3 × 10⁻⁹ mol cm⁻² yr⁻¹).     

The origin and mechanisms of salinization of the lower Jordan river

The chemical and isotopic (⁸⁷Sr/⁸⁶Sr, δ¹¹B, δ³⁴S_sulfate, δ¹⁸O_water, δ¹⁵N_nitrate) compositions of water from the Lower Jordan River and its major tributaries between the Sea of Galilee and the Dead Sea were determined in order to reveal the origin of the salinity of the Jordan River. We identified three separate hydrological zones along the flow of the river:

(1)

A northern section (20 km downstream of its source) where the base flow composed of diverted saline and wastewaters is modified due to discharge of shallow sulfate-rich groundwater, characterized by low ⁸⁷Sr/⁸⁶Sr (0.7072), δ³⁴S_sulfate (−2‰), high δ¹¹B (∼36‰), δ¹⁵N_nitrate (∼15‰) and high δ¹⁸O_water (−2 to-3‰) values. The shallow groundwater is derived from agricultural drainage water mixed with natural saline groundwater and discharges to both the Jordan and Yarmouk rivers. The contribution of the groundwater component in the Jordan River flow, deduced from mixing relationships of solutes and strontium isotopes, varies from 20 to 50% of the total flow.

(2)

A central zone (20-50 km downstream from its source) where salt variations are minimal and the rise of ⁸⁷Sr/⁸⁶Sr and SO₄/Cl ratios reflects predominance of eastern surface water flows.

(3)

A southern section (50-100 km downstream of its source) where the total dissolved solids of the Jordan River increase, particularly during the spring (70-80 km) and summer (80-100 km) to values as high as 11.1 g/L. Variations in the chemical and isotopic compositions of river water along the southern section suggest that the Zarqa River (⁸⁷Sr/⁸⁶Sr∼0.70865; δ¹¹B∼25‰) has a negligible affect on the Jordan River. Instead, the river quality is influenced primarily by groundwater discharge composed of sulfate-rich saline groundwater (Cl^-=31-180 mM; SO₄/Cl∼0.2-0.5; Br/Cl∼2-3×10^-3; ⁸⁷Sr/⁸⁶Sr∼0.70805; δ¹¹B∼30‰; δ¹⁵N_nitrate ∼17‰, δ³⁴S_sulfate=4-10‰), and Ca-chloride Rift valley brines (Cl^-=846-1500 mM; Br/Cl∼6-8×10^-3; ⁸⁷Sr/⁸⁶Sr∼0.7080; δ¹¹B>40‰; δ³⁴S_sulfate=4-10‰). Mixing calculations indicate that the groundwater discharged to the river is composed of varying proportions of brines and sulfate-rich saline groundwater. Solute mass balance calculations point to a ∼10% contribution of saline groundwater (Cl⁻=282 to 564 mM) to the river. A high nitrate level (up to 2.5 mM) in the groundwater suggests that drainage of wastewater derived irrigation water is an important source for the groundwater. This irrigation water appears to leach Pleistocene sediments of the Jordan Valley resulting in elevated sulfate contents and altered strontium and boron isotopic compositions of the groundwater that in turn impacts the water quality of the lower Jordan River.

     

Iron isotope fractionation by Fe(II)-oxidizing photoautotrophic bacteria

Photoautotrophic bacteria that oxidize ferrous iron (Fe[II]) under anaerobic conditions are thought to be ancient in origin, and the ferric (hydr)oxide mineral products of their metabolism are likely to be preserved in ancient rocks. Here, two enrichment cultures of Fe(II)-oxidizing photoautotrophs and a culture of the genus Thiodictyon were studied with respect to their ability to fractionate Fe isotopes. Fe isotope fractionations produced by both the enrichment cultures and the Thiodictyon culture were relatively constant at early stages of the reaction progress, where the ⁵⁶Fe/⁵⁴Fe ratios of poorly crystalline hydrous ferric oxide (HFO) metabolic products were enriched in the heavier isotope relative to aqueous ferrous iron (Fe[II]_aq) by ∼1.5 ± 0.2‰. This fractionation appears to be independent of the rate of photoautotrophic Fe(II)-oxidation, and is comparable to that observed for Fe isotope fractionation by dissimilatory Fe(III)-reducing bacteria. Although there remain a number of uncertainties regarding how the overall measured isotopic fractionation is produced, the most likely mechanisms include (1) an equilibrium effect produced by biological ligands, or (2) a kinetic effect produced by precipitation of HFO overlaid upon equilibrium exchange between Fe(II) and Fe(III) species. The fractionation we observe is similar in direction to that measured for abiotic oxidation of Fe(II)_aq by molecular oxygen. This suggests that the use of Fe isotopes to identify phototrophic Fe(II)-oxidation in the rock record may only be possible during time periods in Earth’s history when independent evidence exists for low ambient oxygen contents.     

Thermochemistry of jarosite-alunite and natrojarosite-natroalunite solid solutions

The thermochemistry of jarosite-alunite and natrojarosite-natroalunite solid solutions was investigated. Members of these series were either coprecipitated or synthesized hydrothermally and were characterized by XRD, FTIR, electron microprobe analysis, ICP-MS, and thermal analysis. Partial alkali substitution and vacancies on the Fe/Al sites were observed in all cases, and the solids studied can be described by the general formula K_1-x-yNa_y(H₃O)_xFe_zAl_w(SO₄)₂(OH)_6-3(3-z-w)(H₂O)_3(3-z-w). A strong preferential incorporation of Fe over Al in the jarosite/alunite structure was observed. Heats of formation from the elements, ΔH°_f, were determined by high-temperature oxide melt solution calorimetry. The solid solutions deviate slightly from thermodynamic ideality by exhibiting positive enthalpies of mixing in the range 0 to +11 kJ/mol. The heats of formation of the end members of both solid solutions were derived. The values ΔH°_f = −3773.6 ± 9.4 kJ/mol, ΔH°_f = −4912.2 ± 24.2 kJ/mol, ΔH°_f = −3734.6 ± 9.7 kJ/mol and ΔH°_f = −4979.7 ± 7.5kJ/mol were found for K_0.85(H₃O)_0.15Fe_2.5(SO₄)₂(OH)_4.5(H₂O)_1.5, K_0.85(H₃O)_0.15Al_2.5(SO₄)₂(OH)_4.5(H₂O)_1.5, Na_0.7(H₃O)_0.3Fe_2.7(SO₄)₂(OH)_5.1(H₂O)_0.9, and Na_0.7(H₃O)_0.3Al_2.7(SO₄)₂(OH)_5.1(H₂O)_0.9 respectively. To our knowledge, this is the first experimentally-based report of ΔH°_f for such nonstoichiometric alunite and natroalunite samples. These thermodynamic data should prove helpful to study, under given conditions, the partitioning of Fe and Al between the solids and aqueous solution.