THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

Oxygen isotopic compositions of Allende Type C CAIs: Evidence for isotopic exchange during nebular melting and asteroidal metamorphism

In situ oxygen isotopic measurements of primary and secondary minerals in Type C CAIs from the Allende CV3 chondrite reveal that the pattern of relative enrichments and depletions of ¹⁶O in the primary minerals within each individual CAI are similar to the patterns observed in Types A and B CAIs from the same meteorite. Spinel is consistently the most ¹⁶O-rich (Δ¹⁷O = −25‰ to −15‰), followed by Al,Ti-dioside (Δ¹⁷O = −20‰ to −5‰) and anorthite (Δ¹⁷O = −15‰ to 0‰). Melilite is the most ¹⁶O-depleted primary mineral (Δ¹⁷O = −5‰ to −3‰). We conclude that the original melting event that formed Type C CAIs occurred in a ¹⁶O-rich (Δ¹⁷O  −20‰) nebular gas and they subsequently experienced oxygen isotopic exchange in a ¹⁶O-poor reservoir. At least three of these (ABC, TS26F1 and 93) experienced remelting at the time and place where chondrules were forming, trapping and partially assimilating ¹⁶O-poor chondrule fragments. The observation that the pyroxene is ¹⁶O-rich relative to the feldspar, even though the feldspar preceded it in the igneous crystallization sequence, disproves the class of CAI isotopic exchange models in which partial melting of a ¹⁶O-rich solid in a ¹⁶O-poor gas is followed by slow crystallization in that gas. For the typical (not associated with chondrule materials) Type C CAIs as well for as the Types A and B CAIs, the exchange that produced internal isotopic heterogeneity within each CAI must have occurred largely in the solid state. The secondary phases grossular, monticellite and forsterite commonly have similar oxygen isotopic compositions to the melilite and anorthite they replace, but in one case (CAI 160) grossular is ¹⁶O-enriched (Δ¹⁷O = −10‰ to −6‰) relative to melilite (Δ¹⁷O = −5‰ to −3‰), meaning that the melilite and anorthite must have exchanged its oxygen subsequent to secondary alteration. This isotopic exchange in melilite and anorthite likely occurred on the CV parent asteroid, possibly during fluid-assisted thermal metamorphism.     

Seasonal patterns of carbon chemistry and isotopes in tufa depositing groundwaters of southwestern Japan

Annually laminated carbonates, known as tufas, commonly develop in limestone areas and typically record seasonal patterns of oxygen- and carbon-isotope compositions. δ¹⁸O values are principally controlled by seasonal changes of water temperature, whereas δ¹³C values are the result of complex reactions among the gaseous, liquid, and solid sources of carbon in the system. We examined the processes that cause the seasonal patterns of δ¹³C in groundwater systems at three tufa-depositing sites in southwestern Japan by applying model calculations to geochemical data. Underground inorganic carbon species are exchanged with gaseous CO₂, which is mainly introduced to the underground hydrological system by natural atmospheric ventilation and by diffusion of soil air. These processes control the seasonal pattern of δ¹³C, which is low in summer and high in winter. Among the three sites we investigated, we identified two extreme cases of the degree of carbon exchange between liquid and gaseous phases. For the case with high radiocarbon composition (Δ¹⁴C) and low pCO₂, there was substantial carbon exchange because of a large contribution of atmospheric CO₂ and a small water mass. For the other extreme case, which was characterized by low Δ¹⁴C and high pCO₂, the contribution of atmospheric CO₂ was small and the water mass was relatively large. Our results suggest that at two of the three sites water residence time within the soil profile was longer than 1 year. Our results also suggested a short residence time (less than 1 year) of water in the soil profile at the site with the smallest water mass, which is consistent with large seasonal amplitude of the springwater temperature variations. The Δ¹⁴C value of tufas is closely related to the hydrological conditions in which they are deposited. If the initial Δ¹⁴C value of a tufa-depositing system is stable, ¹⁴C-chronology can be used to date paleo-tufas.     

The age of the martian meteorite Northwest Africa 1195 and the differentiation history of the shergottites

Samarium-neodymium isotopic analyses of unleached and acid-leached mineral fractions from the recently identified olivine-bearing shergottite Northwest Africa 1195 yield a crystallization age of 347 ± 13 Ma and an value of +40.1 ± 0.9. Maskelynite fractions do not lie on the Sm-Nd isochron and appear to contain a martian surface component with low ¹⁴⁷Sm/¹⁴⁴Nd and ¹⁴³Nd/¹⁴⁴Nd ratios that was added during shock. The Rb-Sr system is disturbed and does not yield an isochron. Terrestrial Sr appears to have affected all of the mineral fractions, although a maximum initial ⁸⁷Sr/⁸⁶Sr ratio of 0.7016 is estimated by passing a 347 Ma reference line through the maskelynite fraction that is least affected by contamination. The high initial value and the low initial ⁸⁷Sr/⁸⁶Sr ratio, combined with the geologically young crystallization age, indicate that Northwest Africa 1195 is derived from a source region characterized by a long-term incompatible-element depletion.The age and initial Sr and Nd isotopic compositions of Northwest Africa 1195 are very similar to those of Queen Alexandra Range 94201, indicating these samples were derived from source regions with similar Sr-Nd isotopic systematics. These similarities suggest that these two meteorites share a close petrogenetic relationship and might have been erupted from a common volcano. The meteorites Yamato 980459, Dar al Gani 476, Sayh al Uhaymir 005/008, and Dhofar 019 also have relatively old ages between 474 and 575 Ma and trace element and/or isotopic systematics that are indicative of derivation from incompatible-element-depleted sources. This suggests that the oldest group of meteorites is more closely related to one another than they are to the younger meteorites that are derived from less incompatible-element-depleted sources. Closed-system fractional crystallization of this suite of meteorites is modeled with the MELTS algorithm using the bulk composition of Yamato 980459 as a parent. These models reproduce many of the major element and mineralogical variations observed in the suite. In addition, the rare earth element systematics of these meteorites are reproduced by fractional crystallization using the proportions of phases and extents of crystallization that are calculated by MELTS. Other shergottites that demonstrate enrichments in incompatible-elements and have evolved Sr and Nd isotopic systematics have some geochemical systematics that are similar to those observed in the depleted group. Most notably, although they exhibit a very limited range of incompatible trace element and isotopic compositions, they have highly variable major element compositions. This is also consistent with evolution from a common mantle source region by variable amounts of fractional crystallization. If this scenario is correct, it suggests that the combined effects of source composition and fractional crystallization are likely to account for the major element, trace element, and isotopic diversity of all shergottites.     

Evidence for crustal degassing of CF4 and SF6 in Mojave Desert groundwaters

Dissolved tetrafluoromethane (CF₄) and sulfur hexafluoride (SF₆) concentrations were measured in groundwater samples from the Eastern Morongo Basin (EMB) and Mojave River Basin (MRB) located in the southern Mojave Desert, California. Both CF₄ and SF₆ are supersaturated with respect to equilibrium with the preindustrial atmosphere at the recharge temperatures and elevations of the Mojave Desert. These observations provide the first in situ evidence for a flux of CF₄ from the lithosphere. A gradual basin-wide enhancement in dissolved CF₄ and SF₆ concentrations with groundwater age is consistent with release of these gases during weathering of the surrounding granitic alluvium. Dissolved CF₄ and SF₆ concentrations in these groundwaters also contain a deeper crustal component associated with a lithospheric flux entering the EMB and MRB through the underlying basement. The crustal flux of CF₄, but not of SF₆, is enhanced in the vicinity of local active fault systems due to release of crustal fluids during episodic fracture events driven by local tectonic activity. When fluxes of CF₄ and SF₆ into Mojave Desert groundwaters are extrapolated to the global scale they are consistent, within large uncertainties, with the fluxes required to sustain the preindustrial atmospheric abundances of CF₄ and SF₆.     

Low-oxygen and chemical kinetic constraints on the geochemical niche of neutrophilic iron(II) oxidizing microorganisms

Neutrophilic iron oxidizing bacteria (FeOB) must actively compete with rapid abiotic processes governing Fe(II) oxidation and as a result have adapted to primarily inhabit low-O₂ environments where they can more successfully compete with abiotic Fe(II) oxidation. The spatial distribution of these microorganisms can be observed through the chemical gradients they affect, as measured using in situ voltammetric analysis for dissolved Fe(II), Fe(III), O₂, and FeS_(aq). Field and laboratory determination of the chemical environments inhabited by the FeOB were coupled with detailed kinetic competition studies for abiotic and biotic oxidation processes using a pure culture of FeOB to quantify the geochemical niche these organisms inhabit. In gradient culture tubes, the maximum oxygen levels, which were associated with growth bands of Sideroxydans lithotrophicus (ES-1, a novel FeOB), were 15-50 μM. Kinetic measurements made on S. lithotrophicus compared biotic/abiotic (killed control) Fe oxidation rates. The biotic rate can be a significant and measurable fraction of the total Fe oxidation rate below O₂ concentrations of approximately 50 μM, but biotic Fe(II) oxidation (via the biotic/abiotic rate comparison) becomes difficult to detect at higher O₂ levels. These results are further supported by observations of conditions supporting FeOB communities in field settings. Variablity in cell densities and cellular activity as well as variations in hydrous ferrous oxide mineral quantities significantly affect the laboratory kinetic rates. The microbial habitat (or geochemical niche) where FeOB are active is thus largely controlled by the competition between abiotic and biotic kinetics, which are dependent on Fe(II) concentration, P_O2, temperature and pH in addition to the surface area of hydrous ferric oxide minerals and the cell density/activity of FeOB. Additional field and lab culture observations suggest a potentially important role for the iron-sulfide aqueous molecular cluster, FeS_(aq), in the overall cycling of iron associated with the environments these microorganisms inhabit.     

Geochemistry, petrology and ages of the lunar meteorites Kalahari 008 and 009: New constraints on early lunar evolution

Kalahari 008 and 009 are two lunar meteorites that were found close to each other in Botswana. Kalahari 008 is a typical lunar anorthositic breccia; Kalahari 009 a monomict breccia with basaltic composition and mineralogy. Based on minor and trace elements Kalahari 009 is classified as VLT (very-low-Ti) mare basalt with extremely low contents of incompatible elements, including the REE. The Lu-Hf data define an age of 4286 ± 95 Ma indicating that Kalahari 009 is one of the oldest known basalt samples from the Moon. It provides evidence for lunar basalt volcanism prior to 4.1 Ga (pre-Nectarian) and may represent the first sample from a cryptomare. The very radiogenic initial ¹⁷⁶Hf/¹⁷⁷Hf (εHf = +12.9 ± 4.6), the low REE, Th and Ti concentrations indicate that Kalahari 009 formed from re-melting of mantle material that had undergone strong incompatible trace element depletion early in lunar history. This unusually depleted composition points toward a hitherto unsampled basalt source region for the lunar interior that may represent a new depleted endmember source for low-Ti mare basalt volcanism. Apparently, the Moon became chemically very heterogeneous at an early stage in its history and different cumulate sources are responsible for the diverse mare basalt types.Evidence that Kalahari 008 and 009 may be paired includes the similar fayalite content of their olivine, the identical initial Hf isotope composition, the exceptionally low exposure ages of both rocks and the fact that they were found close to each other. Since cryptomaria are covered by highland ejecta, it is possible that these rocks are from the boundary area, where basalt deposits are covered by highland ejecta. The concentrations of cosmogenic radionuclides and trapped noble gases are unusually low in both rocks, although Kalahari 008 contains slightly higher concentrations. A likely reason for this difference is that Kalahari 008 is a polymict breccia containing a briefly exposed regolith, while Kalahari 009 is a monomict brecciated rock that may never have been at the surface of the Moon.Altogether, the compositions of Kalahari 008 and 009 permit new insight into early lunar evolution, as both meteorites sample lunar reservoirs hitherto unsampled by spacecraft missions. The very low Th and REE content of Kalahari 009 as well as the depletion in Sm and the lack of a KREEP-like signature in Kalahari 008 point to a possible source far from the influence of the Procellarum-KREEP Terrane, possibly the lunar farside.