New evidence for the prograde and retrograde PT-path of high-pressure granulites, Moldanubian Zone, Lower Austria, by Zr-in-rutile thermometry and garnet diffusion modelling

Abstract Compositional zoning in garnet, mineral inclusions and the application of the Zr-in-rutile thermometry on rutile inclusions in garnet in combination with conventional geothermobarometry and thermodynamic modelling allows a reconstruction of the prograde pressure-temperature evolution in felsic and mafic high-pressure granulites from the Moldanubian Zone, Bohemian Massif, Lower Austria. Most garnets in these rocks show homogeneous core compositions with high grossular contents (~30 mol%), while their rim zones have a markedly reduced grossular content. Rutile inclusions in the grossular rich garnet cores have low Zr concentrations (400 to 1300 ppm) indicating a formation temperature of ~810–820 °C which implies that the garnet host grew at these temperature conditions as well. Based on numerous polycrystalline melt inclusions, high Ti-biotite relics and a generally high Ti concentration in garnet cores, the peritectic biotite breakdown reaction is considered to be responsible for a first garnet growth, now observed as high-grossular garnet cores. The corresponding pressure is estimated to be in the range of 1.6 to 2.5 GPa, based on experimentally determined biotite breakdown reactions, thermodynamic modelling and the occurrence of high-Ti biotite in garnet cores. Rutile inclusions in low-Ca garnet rims contain significantly higher Zr concentrations (1700 to 5800 ppm) resulting in ultrahigh temperatures of ~1030 °C. Similar temperature as well as corresponding pressure estimates of 1000 ± 50 °C and 1.60 ± 0.10 GPa were obtained by geothermobarometry and thermodynamic modelling using garnet rim and re-integrated ternary feldspar compositions. These high pressure and ultrahigh temperature conditions are well known from literature for these granulites. The proposed two-phase garnet growth is not only seen in different temperatures obtained from rutile inclusions in garnet core and rim areas, but also in discontinuous trace (Cr, Ga, P, Ti, V, Zr) and heavy rare earth element profiles across garnet porphyroblasts, implying a different reaction mechanism for garnet rim growth. This second phase of garnet growth must have occurred during near isobaric heating to the ultrahigh temperature peak, most likely even at slightly lower pressures compared to the garnet core growth. By applying a binary Fe-Mg diffusion model to strongly zoned garnet grains a maximum timescale of 5–6 million years was estimated for the exhumation and cooling process, assuming a linear cooling path from 1000 °C at 1.6 GPa to 760 °C at 0.8 GPa. This short-lived ultrahigh temperature event corresponds to cooling and exhumation rates of 40–50 °C Ma−1 and 5.3–6.6 mm y−1, respectively.