Timing of UHP metamorphism in the Hong’an area,western Dabie Mountains,China: evidence from zircon U–Pb age,trace element and Hf isotope composition

The Hong’an area (western Dabie Mountains) is the westernmost terrane in the Qinling-Dabie-Sulu orogen that preserves UHP eclogites. The ages of the UHP metamorphism have not been well constrained, and thus hinder our understanding of the tectonic evolution of this area. LA-ICPMS U–Pb age, trace element and Hf isotope compositions of zircons of a granitic gneiss and an eclogite from the Xinxian UHP unit in the Hong’an area were analyzed to constrain the age of the UHP metamorphism. Most zircons are unzoned or show sector zoning. They have low trace element concentrations, without significant negative Eu anomalies. These metamorphic zircons can be further subdivided into two groups according to their U–Pb ages, and trace element and Lu–Hf isotope compositions. One group with an average age of 239 ± 2 Ma show relatively high and variable HREE contents (527 ≥ Lu_N ≥ 14) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.00008–0.000931), indicating their growth prior to a great deal of garnet growth in the late stage of continental subduction. The other group yields an average age of 227 ± 2 Ma, and shows consistent low HREE contents and ¹⁷⁶Lu/¹⁷⁷Hf ratios, suggesting their growth with concurrent garnet crystallization and/or recrystallization. These two groups of age are taken as recording the time of prograde HP to UHP and retrograde UHP–HP stages, respectively. A few cores have high Th/U ratios, high trace element contents, and a clear negative Eu anomaly. These features support a magmatic origin of these zircon cores. The upper intercept ages of 771 ± 86 and 752 ± 70 Ma for the granitic gneiss and eclogite, respectively, indicate that their protoliths probably formed as a bimodal suite in rifting zones in the northern margin of the Yangtze Block. Young Hf model ages (T _DM1) of magmatic cores indicate juvenile (mantle-derived) materials were involved in their protolith formation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Influence of temperature,composition, silica activity and oxygen fugacity on the H<Subscript>2</Subscript>O storage capacity of olivine at 8 GPa

Olivine crystals were grown in the presence of a hydrous silicate fluid during multi-anvil experiments at 8 GPa and 1,000–1,600°C. Experiments were conducted both in a simple system (FeO–MgO–SiO₂–H₂O) and in a more complex system containing additional elements (CaO–Na₂O–Al₂O₃–Cr₂O₃–TiO₂–FeO–MgO–SiO₂–H₂O). Silica activity was buffered by the presence of either pyroxene (high a _SiO2) or ferropericlase (low a _SiO2), and was buffered by the presence of Ni + NiO or Fe + FeO, or constrained by the presence of Fe₂O₃. Raman spectroscopy was used to identify pyroxene polymorphs in the run products. Clinoenstatite was present in the 1,000°C experiment, and enstatite in experiments at 1,400–1,520°C. The H₂O content of olivine was measured using secondary ion mass spectroscopy, and infrared spectroscopy was used to investigate the nature of hydrous defects. The H₂O storage capacity of olivine decreases with increasing temperature at 8 GPa. In contrast to previous experimental results at ≤2 GPa, no significant effect of varying oxygen fugacity is evident, but H₂O storage capacity is enhanced under conditions of low silica activity. No significant growth of low wavenumber (<3,400 cm⁻¹) peaks, generally associated with high at low pressure, was observed in the FTIR spectra of olivine from the high experiments. Our experiments show that previous high pressure H₂O storage capacity measurements for olivine synthesized under more oxidizing conditions than the Earth’s mantle are not likely to be compromised by the of the experiments. However, the considerable effect of temperature on H₂O storage capacity in olivine must be taken into account to avoid overestimation of the bulk upper mantle H₂O storage capacity.     

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.     

Development and morphometry of sinkholes in coastal plains of Apulia, southern Italy. Preliminary sinkhole susceptibility assessment

Evolution of coastlines in karst areas may be strongly controlled by dissolution processes which favour the development of surface and subsurface landforms. The generation of caves in these environments is commonly favoured by the mixing between fresh and brackish waters. The sinkholes resulting from the upward propagation of the caves may interfere with the anthropogenic environment and cause damage to human elements (property and activities). To highlight the often underestimated importance of karst phenomena in coastal areas, we have analyzed a coastal stretch of Apulia, in southern Italy. The study area, covering an extension of about 6 km², is situated in the Ionian coast, and presents several interesting karst landforms that are generally connected to caves. Tens of sinkholes were mapped through field surveys, multi-year aerial-photographs (dating back to the 1940s) and archival research. We have performed a morphometric analysis of the sinkholes. The analysis describes the main parameters of the sinkholes (area, length, width, and depth), and the control exerted by the main discontinuity systems in the area. The detrimental effects derived from interaction between human environment and these karst landforms is also under consideration. A sinkhole susceptibility map, which may provide useful information for planners, developers and the insurance industry has eventually been produced through application of a decision tree model.     

Formation of Zn–Ca phyllomanganate nanoparticles in grass roots

It is now well established that a number of terrestrial and aquatic microorganisms have the capacity to oxidize and precipitate Mn as phyllomanganate. However, this biomineralization has never been shown to occur in plant tissues, nor has the structure of a natural Mn(IV) biooxide been characterized in detail. We show that the graminaceous plant Festuca rubra (red fescue) produces a Zn-rich phyllomanganate with constant Zn:Mn and Ca:Mn atomic ratios (0.46 and 0.38, respectively) when grown on a contaminated sediment. This new phase is so far the Zn-richest manganate known to form in nature (chalcophanite has a Zn:Mn ratio of 0.33) and has no synthetic equivalent. Visual examination of root fragments under a microscope shows black precipitates about ten to several tens of microns in size, and their imaging with backscattered and secondary electrons demonstrates that they are located in the root epidermis. In situ measurements by Mn and Zn K-edge extended X-ray absorption fine structure (EXAFS) spectroscopy and X-ray diffraction (XRD) with a micro-focused beam can be quantitatively described by a single-phase model consisting of Mn(IV) octahedral layers with 22% vacant sites capped with tetrahedral and octahedral Zn in proportions of 3:1. The layer charge deficit is also partly balanced by interlayer Mn and Ca. Diffracting crystallites have a domain radius of 33 Å in the ab plane and contain only 1.2 layers (8.6 Å) on average. Since this biogenic Mn oxide consists mostly of isolated layers, basal 00l reflections are essentially absent despite its lamellar structure. Individual Mn layers are probably held together in the Mn–Zn precipitates by stabilizing organic molecules. Zinc biomineralization by plants likely is a defense mechanism against toxicity induced by excess concentrations of this metal in the rhizosphere.     

Effectiveness of sulfate-reducing passive bioreactors for treating highly contaminated acid mine drainage: I. Effect of hydraulic retention time

Sulfate-reducing passive bioreactors have proved to be an effective technology for the treatment of acid mine drainage (AMD) contaminated waters over relatively short periods of time (1–5 a). However, long-term efficiency can be limited by several factors including problems related to the hydraulic properties of the reactive mixture. In this study, the effect of two hydraulic retention times (HRTs) of 7.3 d and 10 d on the performance of passive bioreactors was evaluated over an 11-month period for the treatment of a highly contaminated AMD. Evolution of the porosity and hydraulic conductivity of the reactive mixture was also evaluated during the 15-month operation of two bioreactors. Results indicated that bioreactors were effective at both HRTs for increasing the pH and alkalinity of contaminated water and for SO₄ and metal removal (60–82% for Fe and up to 99.9% for Cd, Ni and Zn). Although the quality of treated effluent was significantly improved with the 10 d HRT compared to the 7.3 d HRT, results showed that the higher HRT reduced the porosity and the permeability of the reactive mixture which might lead to hydraulic related problems and, eventually, to limited efficiency in long-term operation compared to a shorter HRT. The choice of HRT for a passive bioreactor must therefore consider both the desired quality of treated effluent and the potential for deterioration of hydraulic properties in the reactive mixture.     

Effectiveness of sulfate-reducing passive bioreactors for treating highly contaminated acid mine drainage: II. Metal removal mechanisms and potential mobility

Column bioreactors were used for studying mechanisms of metal removal, assessment of long-term stability of spent reactive mixtures, as well as potential metal mobility after treating highly contaminated acid mine drainage (AMD; pH 2.9–5.7). Several physicochemical, microbiological, and mineralogical analyses were performed on spent reactive mixtures collected from 4 bioreactors, which were tested in duplicate for two hydraulic retention times (7.3d and 10d), with downward flow over an 11-month period. Consistent with the high metal concentrations in the AMD feed, and with low metal concentrations measured in the treated effluent, the physicochemical analyses indicated very high concentrations of metals (Fe, Mn, Cd, Ni, and Zn) in the top and bottom layers of the reactive mixtures from all columns. Moreover, the concentrations of Fe (50.8–57.8 g/kg) and Mn (0.53–0.70 g/kg) were up to twice as high in the bottom layers, whereas the concentrations of Cd (6.77–13.3 g/kg), Ni (1.80–5.19 g/kg) and Zn (2.53–13.2 g/kg) were up to 50-times higher in the top layers. Chemical extractions and elemental analysis gave consistent results, which indicated a low fraction of metals removed as sulfides (up to 15% of total metals recovered in spent reactive mixtures). Moreover, Fe and Mn were found in a more stable chemical form (residual fraction was 42–74% for Mn and 30–77% for Fe) relative to Cd, Ni or Zn, which seemed more weakly bound (oxidisable/reducible fractions) and showed higher potential mobility. Besides identifying (oxy)hydroxide and carbonate minerals, the mineralogical analyses identified metal sulfides containing Fe, Cd, Ni and Zn. Metal removal mechanisms were, therefore, mainly adsorption and other binding mechanisms with organic matter (for Cd, Ni and Zn), and the precipitation as (oxy)hydroxide minerals (for Fe and Mn). After 15 months, however, the column bioreactors did not lose their capacity for removing metals from the AMD. Although the metals were immobile during the bioreactor treatment, their mobility could increase from spent reactive mixtures, if stored inappropriately. Metal recovery by acidic leaching of spent substrates at the end of bioreactor operation could be an alternative.