The prominent felsic granulites in the southern part of the Bohemian Massif (Gföhl Unit, Moldanubian Zone), with the Variscan (∼340 Ma) high-pressure and high-temperature assemblage garnet+quartz+hypersolvus feldspar ± kyanite, correspond geochemically to slightly peraluminous, fractionated granitic rocks. Compared to the average upper crust and most granites, the U, Th and Cs concentrations are strongly depleted, probably because of the fluid and/or slight melt loss during the high-grade metamorphism (900–1050°C, 1·5–2·0 GPa). However, the rest of the trace-element contents and variation trends, such as decreasing Sr, Ba, Eu, LREE and Zr with increasing SiO2 and Rb, can be explained by fractional crystallisation of a granitic magma. Low Zr and LREE contents yield ∼750°C zircon and monazite saturation temperatures and suggest relatively low-temperature crystallisation. The granulites contain radiogenic Sr (87Sr/86Sr340 = 0·7106–0·7706) and unradiogenic Nd ( = − 4·2 to − 7·5), indicating derivation from an old crustal source. The whole-rock Rb–Sr isotopic system preserves the memory of an earlier, probably Ordovician, isotopic equilibrium.Contrary to previous studies, the bulk of felsic Moldanubian granulites do not appear to represent separated, syn-metamorphic Variscan HP–HT melts. Instead, they are interpreted as metamorphosed (partly anatectic) equivalents of older, probably high-level granites subducted to continental roots during the Variscan collision. Protolith formation may have occurred within an Early Palaeozoic rift setting, which is documented throughout the Variscan Zone in Europe.

No abstract

We present a thermodynamic model for mineral dissolution in aqueous fluids at elevated temperatures and pressures, based on intrinsic thermal properties and variations of volumetric properties of the aqueous solvent. The standard thermodynamic properties of mineral dissolution into aqueous fluid consist of two contributions: one from the energy of transformation from the solid to the hydrated-species state and the other from the compression of solvent molecules during the formation of a hydration shell. The latter contribution has the dimension of the generalized Krichevskii parameter. This approach describes the energetics of solvation more accurately than does the Born electrostatic theory and can be extended beyond the limits of experimental measurements of the dielectric constant of H 2 O. The new model has been calibrated by experimental solubilities of quartz, corundum, rutile, calcite, apatite, fluorite and portlandite in pure H 2 O at temperatures up to 1100°C and pressures up to 20 kbar. All minerals show a steady increase in solubility along constant geothermal gradients or water isochores. By contrast, isobaric solubilities initially increase with rising temperature but then decline above 200-400°C. This retrograde behavior is caused by variations in the isobaric expansivity of the aqueous solvent, which approaches infinity at its critical point. Oxide minerals predominantly dissolve to neutral species; so, their dissolution energetics involve a relatively small contribution from the solvent volumetric properties and their retrograde solubilities are restricted to a relatively narrow window of temperature and pressure near the critical point of water. By contrast, Ca-bearing minerals dissolve to a variety of charged species; so, the energetics of their dissolution reactions involve a comparatively large contribution from volume changes of the aqueous solvent and their isobaric retrograde solubility spans nearly all metamorphic and magmatic conditions. These features correlate with and can be predicted from the standard partial molar volumes of aqueous species.The thermodynamic model can be used over much wider range of settings for terrestrial fluid-rock interaction than has previously been possible. To illustrate, it is integrated with transport theory to show quantitatively that integrated fluid fluxes characteristic of crustal shear zones are capable of precipitating quartz or calcite veins from low-and medium-grade metamorphic conditions, at a geothermal gradient of 20°C km )1 . For subduction zones, modeled by a geotherm of 7°C km )1 , the required fluid fluxes are one to two orders of magnitude lower and predict enhanced efficiency of mass transfer and metasomatic precipitation in comparison with orogenic settings. The new model thus can be applied to shallow hydrothermal, metamorphic, magmatic and subduction fluids, and for retrieval of dependent thermodynamic properties for mass transfer or geodynamic modeling.