Eclogite thermobarometry: The consistency between conventional thermobarometry and forward phase-equilibrium modelling

Eclogite thermobarometry is crucial for constraining the depths and temperatures to which oceanic and continental crust subduct. However, obtaining the pressure and temperature (P–T) conditions of eclogites is complex as they commonly display high-variance mineral assemblages, and the mineral compositions only vary slightly with P–T. In this contribution, we present a comparison between two independent and commonly used thermobarometric approaches for eclogites: conventional thermobarometry and forward phase-equilibrium modelling. We assess how consistent the thermobarometric calculations are using the garnet–clinopyroxene–phengite barometer and garnet–clinopyroxene thermometer with predictions from forward modelling (i.e. comparing the relative differences between approaches). Our results show that the overall mismatch in methods is typically ±0.2–0.3 GPa and ±29–42°C although differences as large as 80°C and 0.7 GPa are possible for a few narrow ranges of P–T conditions in the forward models. Such mismatch is interpreted as the relative differences among methods, and not as absolute uncertainties or accuracies for either method. For most of the investigated P–T conditions, the relatively minor differences between methods means that the choice in thermobarometric method itself is less important for geological interpretation than careful sample characterization and petrographic interpretation for deriving P–T from eclogites. Although thermobarometry is known to be sensitive to the assumed X_Fe³⁺ of a rock (or mineral), the relative differences between methods are not particularly sensitive to the choice of bulk-rock X_Fe³⁺, except at high temperatures (>650°C, amphibole absent) and for very large differences in assumed X_Fe³⁺ (0–0.5). We find that the most important difference between approaches is the activity–composition (a–x) relations, as opposed to the end-member thermodynamic data or other aspects of experimental calibration. When equivalent a–x relations are used in the conventional barometer, P calculations are nearly identical to phase-equilibrium models (ΔP < 0.1). To further assess the implications of these results for real rocks, we also evaluate common mathematical optimizations of reaction constants used for obtaining the maximum P–T with conventional thermobarometric approaches (e.g. using the highest aGrs² × aPrp in garnet and Si content in phengite, and the lowest aDi in clinopyroxene). These approaches should be used with caution, because they may not represent the compositions of equilibrium mineral assemblages at eclogite facies conditions and therefore systematically bias P–T calculations. Assuming method accuracy, geological meaningful P_max at a typical eclogite facies temperature of ~660°C will be obtained by using the greatest aDi, aCel, and aPrp and lowest aGrs and aMs; garnet and clinopyroxene with the lowest Fe²⁺/Mg ratios may yield geological meaningful T_max at a typical eclogite facies pressure of 2.5 GPa.     

Thermal expansion and high-temperature P21/c–C2/c phase transition in clinopyroxene-type LiFeGe2O6 and comparison to NaFe(Si,Ge)2O6

A synthetic clinopyroxene with composition LiFe³⁺Ge₂O₆, monoclinic s.g. P2₁/c, a = 9.8792(7), b = 8.8095(5), c = 5.3754(3) Å, β = 108.844(6)°, V = 442.75(16) ų, has been studied by in situ low- and high-temperature single-crystal X-ray diffraction. The variation of lattice parameters and the intensity of the b-type reflections (h + k = 2n + 1, only present in the P-symmetry) with increasing temperature showed a displacive phase transition from space group P2₁/c to C2/c at a transition temperature T _tr = 789 K, first order in character, with a sudden volume increase of 1.6% and a decrease of β by 1° at the transition. This spontaneous dilatation is reversible, shows a limited hysteresis of ±10°C, and corresponds to the vanishing of the b-type reflections, thus indicating a symmetry increase to space group C2/c. Below T _tr an expansion is observed for all the cell parameters, while the β angle remained almost constant; at T > T _tr the thermal volume expansion is due to dilatation of the structure in the $(\bar{1}\,0\,1) A synthetic clinopyroxene with composition LiFe³⁺Ge₂O₆, monoclinic s.g. P2₁/c, a = 9.8792(7), b = 8.8095(5), c = 5.3754(3) ?, β = 108.844(6)°, V = 442.75(16) ?³, has been studied by in situ low- and high-temperature single-crystal X-ray diffraction. The variation of lattice parameters and the intensity of the b-type reflections (h + k = 2n + 1, only present in the P-symmetry) with increasing temperature showed a displacive phase transition from space group P2₁/c to C2/c at a transition temperature T _tr = 789 K, first order in character, with a sudden volume increase of 1.6% and a decrease of β by 1° at the transition. This spontaneous dilatation is reversible, shows a limited hysteresis of ±10°C, and corresponds to the vanishing of the b-type reflections, thus indicating a symmetry increase to space group C2/c. Below T _tr an expansion is observed for all the cell parameters, while the β angle remained almost constant; at T > T _tr the thermal volume expansion is due to dilatation of the structure in the ([`1] 0 1)(\bar{1}\,0\,1) plane, mostly along [0 1 0], and pure shear in the (0 1 0) plane due to the decrease of β. From comparison with silicate analogues, the germanate clinopyroxenes are more expansible, while the P2₁/c expands more than the C2/c phase. The evolution of Q ² (calculated as the normalized intensity of b-type reflections) with T in the framework of the Landau theory has been done using a standard expression for a first order phase transition. We observe a jump of Q ₀² = 0.538(2) at T _tr, with T _c of 481(7) K, b/a = −2,290 K, and c/a = 3,192 K, and thus far from being tri-critical point. A closely related composition (LiFe³⁺Si₂O₆) shows an equivalent phase transition at 228 K, which is very close to the tri-critical point and 561 K cooler. This result indicates that a change in the composition of tetrahedral sites can have dramatic effects on the P2₁/c ↔ C2/c displacive phase transition in clinopyroxenes. The major changes observed in the evolution of the crystal structure with T are observed in the M2 polyhedron, with a volume decrease by ca. 13.3%, compared to ca. 1.3% observed in the M1 polyhedron. The tetrahedra behave as rigid units with neither a significant change of volume at T > T _tr (<1‰), nor a change of tilting of the basal plane. No change in coordination is observed at T > T _tr in the M2 polyhedron, which remains sixfold coordinated although a strong deformation of this polyhedron is observed. This deformation is related to a strong change by 51.4° at T _tr of the kinking angle (O3–O3–O3 angle) of the B-chain of tetrahedra, which switches from O-rotated to S-rotated [from 143.3(5)° to 194.7(6)°]. The A-chain is S-rotated at T < T _tr [206.8(5)° at 703 K] and extends by 12° at the transition.     

Dynamics of Saxothuringian subduction channel/wedge constrained by phase‐equilibria modelling and micro‐fabric analysis

Subduction and exhumation dynamics can be investigated through analysis of metamorphic and deformational evolution of associated high‐grade rocks. The Erzgebirge anticline, which forms at the boundary between the Saxothuringian and Teplá‐Barrandian domains of the Bohemian Massif, provides a useful study area for these processes owing to the occurrence of numerous meta‐basites preserving eclogite facies assemblages, and coesite and diamond bearing quartzofeldspathic lithologies indicating subduction to deep mantle depths. The prograde and retrograde evolution of meta‐basite from the Czech portion of the Erzgebirge anticline has been constrained through a combination of thermodynamic modelling and conventional thermobarometry. Garnet growth zoning indicates that the rocks underwent burial and heating to peak conditions of 2.6 GPa and at least 615 °C. Initial exhumation occurred with concurrent cooling and decompression resulting in the growth of amphibole and zoisite poikiloblasts overgrowing and including the eclogite facies assemblage. The development of clinopyroxene–plagioclase–amphibole symplectites after omphacite and Al‐rich rims on matrix amphibole indicate later heating at the base of the lower crust. Omphacite microstructures, in particular grain size analysis and lattice‐preferred orientations, indicate that the prograde evolution was characterized by a constrictional strain geometry transitioning into plane strain and oblate fabrics during exhumation. The initial constrictional strain pattern is interpreted as being controlled by competing slab pull and crustal buoyancy forces leading to necking of the subducting slab. The transition to plane strain and flattening geometries represents transfer of material from the subducting lithosphere into a subduction channel, break‐off of the dense slab and rebound of the buoyant crustal material.     

Wind and freshwater influence over hydrocarbon dispersal on Patos Lagoon, Brazil

The two-dimensional barotropic, hydrodynamic and transport model MOHID is applied to the Patos Lagoon system using a nested modelling approach to reproduce both the lagoon and estuary hydrodynamics. A new Lagrangian oil spill model is presented and used to simulate a hypothetical oil spill in the estuary. Hydrodynamic fields are validated and used to force the oil model. Results show that the hydrodynamics of this system is mainly controlled by the wind and freshwater discharge. The dispersion, concentration and thickness evolution of the oil in the first day after the spill is determined by the equilibrium between these two factors. The freshwater discharge is the major factor controlling the oil dispersion for discharges greater than 5000 m3 while the wind assumes control for lower discharge amounts. The results presented are a first step toward a coastal management tool for the Patos Lagoon.     

Terrestrial heat and volcanism

Summary Although active volcanic territories are not characterised by extreme high heat flow, yet terrestrial heat is one of the main, if not hte first cause of volcanic activities. About 94 per cent of the known active volcanoes are in connection with orogenic zones including ocean ridges, which mean that processes responsible for the evolution of the continents cause the overwhelming majority of volcanism.There stages of the evolution of the continental crust can be distinguished: 1.) Growth of the continental cores by differentiation, 2.) Existence of peripheral growing zones around the first core, 3.) Intergrowth of the existing continents.Volcanism, as well as seismicity are manifestations of the evolution of the crust which is caused by the development of heat sources in the mantle. The heat partly melt the upper mantle and induce slow plastic flow in connection of differentiation. The space, left empty by continental uplift is filled by the slow plastic inflow of mantle from the ocean. This slow redistribution of mantle and curstal rocks disturbs the stress field, the rocks fail along the strained zones and zones of fractures appear sometimes with volcanic activity.     

Relationship between microseisms,geomagnetic activity and ionospheric absorption of radio waves

Summary Relatively good correlation between such remote phenomena as microseisms, geomagnetic activity and ionospheric absorption of radio waves has been observed. These phenomena are probably connected by a considerably increased penetration of energetic particles into the atmosphere under conditions of geomagnetic activity.     

An anthropo-ecological approach to landscape use optimization

A method of landscape evaluation based on the anthropo-ecological approach is worked out. The procedure consists of three stages, viz., general, special and optimization. Within the general stage, the structural arrangement of the landscape is evaluated by means of the concept of geo-ecological sites (GES), which are distributed into four zones (A, B, C, D) according to the degrees of their anthropic disturbance. The zoning performed makes it possible to select those geo-ecological sites where eco-critical relationships can be expected to appear; in order to eliminate such situations, transition zones are calculated for the geo-ecological sites concerned. The outcome of this approach is a proposal of a new structural arrangement of the geo-ecological sites, solving eco-critical situations and serving as a basis for the special and optimization stages of the anthropo-ecological method.