Metamorphic rocks of the Murmansk domain (Kola Peninsula) as compared with the oldest rock associations of the northeastern Baltic shield,Canada, and Greenland

Rock associations characterized by heterogeneous sets of petrogeochemical parameters were compared by quantifying the degree of their similarity-dissimilarity and searching for discrimination trends between them. Using procedures specially developed for this purpose, it was demonstrated for the first time that the lithotectonic complexes of the Murmansk domain are fundamentally different from those of typical granulite-gneiss terrains and resemble the granite-greenstone terrains of the Baltic shield, Greenland, and Canada. Based on the whole data set, the Murmansk domain can be considered as a deeply eroded Archean granite-greenstone terrain retaining only the tonalite-trondhjemite-gneiss basement with abundant supracrustal enclaves. A trend of the compositional difference between the older and younger rock associations is similar to that between the tholeiitic and boninitic volcanic series. It was suggested that the petrogeochemical “age” trend reflects the initial stage of the compositional evolution of the metamagmatic rocks of the region from metamorphic rocks similar to tholeiites at the early stages to the Paleoproterozoic boninite-like rocks, which are believed to be linked to the unique PGE-bearing province of the northeastern Baltic shield. This implies that the specific metallogenic features of the region emerged already in the Archean, which supported the suggestion on the long duration of geological processes in the Early Precambrian.     

Li-Sr-Nd isotope signatures of the plume and cratonic lithospheric mantle beneath the margin of the rifted Tanzanian craton (Labait)

Lithium concentrations and isotopic compositions of olivine and ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd of coexisting clinopyroxene from peridotite xenoliths from the Quaternary Labait volcano, Tanzania, document the influence of rift-related metasomatism on the ancient cratonic mantle. Olivines show negative correlations between Fo content and both δ⁷Li and Li concentrations. Olivines in iron-rich peridotites (Fo_85–87) have high Li concentrations (3.2–4.8 ppm) and heavy δ⁷Li (+5.2 to +6.6). In contrast, olivines in ancient, refractory peridotites have lower Li concentrations (∼2 ppm) and relatively light δ⁷Li (+2.6 to +3.5). This reflects mixing between ancient, refractory cratonic lithosphere and asthenosphere-derived rift magmas. A uniquely fertile, deformed, high-temperature garnet lherzolite, interpreted to be from the base of the lithosphere, has a ⁸⁷Sr/⁸⁶Sr of 0.7029 and ¹⁴³Nd/¹⁴⁴Nd of 0.51286, similar to HIMU oceanic basalts. It provides the best estimate of the Sr–Nd isotope composition of the upwelling mantle (i.e., plume, sensu lato) underlying this portion of the East African Rift, and is slightly less radiogenic compared to previous estimates of the plume that were based on rift basalts. Although elevated δ⁷Li are not exclusive to HIMU source regions, the data collectively indicate that the plume beneath Labait has HIMU characteristics in Sr, Nd and Li isotope composition. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.

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.     

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.     

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.     

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.     

Late Precambrian rifting in the geological history of the western slope of the South Urals

The rift-related geodynamic setting of the Late Precambrian geological evolution on the western slope of the South Urals is reconstructed on the basis of localization of lithotectonic complexes of this age, their formation conditions, and the geochemistry of rocks. The Early Riphean stage comprises accumulation of coarse-clastic rocks intercalating with alkaline volcanic rocks of the Navysh Complex, which is a constituent of the Ai Formation, and emplacement of doleritic and picritic intrusions of the Shuida Complex and melanocratic dolerite and gabbrodolerite of the Yusha Complex. The Middle Riphean stage is characterized by wide-spread coarse-clastic terrigenous rocks of the Mashak Formation that intercalate with volcanic rocks of the bimodal basalt-rhyolite association, the Berdyaush pluton of rapakivi granite, the Kusa-Kopan layered intrusive complex, the Lapyshta Complex of dolerites and picrites, and numerous occurrences of gabbrodolerites. The terrigenous rocks of the Vendian stage include conglomerate, gravelstone, and sandstone of the Asha Group, while igneous rocks comprise alkaline volcanics of the Arsha Complex, alkali gabbroids of the Miseli Complex, and melanocratic syenite of the Avashla Complex. The geological evolution of the region is distinguished by local (failed or aborted) rifting. The occurrence of lithotectonic complexes is controlled by dynamic conditions of rifting. A certain inheritance in the evolution may be traced for the Early and Middle Riphean and partly for the Late Riphean and Vendian.