Magnitude–distance relations for earthquake-induced landslides in Greece

A data set of 47 earthquake-induced landslides occurring in Greece from AD 1650 to 1995 has been compiled, and the landslide distribution has been examined. The spatial distribution indicates landslides occurrence almost everywhere in Greece with the exception of the north Greek mainland, which is likely due to the low occurrence frequency of large earthquakes. The earthquake magnitude, M_s, ranges between 5.3 and 7.9, while the main bulk of events falls within the range 7.5≥M_s≥5.3 with peaks at M_s=6.4 and 6.7. The upper bound of the maximum epicentral distance, R_e, from a reported landslide has been approximated by the line log(R_e)=−2.98+0.75M_s for M_s≥5.3, which is compatible with curves proposed world-wide and constitutes a realistic solution at least for hazard assessment purposes. We discuss the importance of the results for assessing hazards from earthquake-induced landslides and for calculating focal parameters of historical earthquakes.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

2019年5月18日松原M5.1地震构造机制分析

基于区域地震台网的数字化波形资料,使用ISOLA方法对2019年5月18日吉林松原M5.1地震进行矩张量反演,研究地震的震源机制,并且收集了地震序列中M_L2.5以上地震的震源机制解,采用FMSI（focal mechanism stress inversion）方法反演震中区构造应力场。结果显示：松原M5.1地震的矩震级为4.9,矩心深度为6 km,双力偶分量为91.5%,主压应力P轴方位角、倾角分别为76°和3°,主张应力T轴方位角、倾角分别为166°和16°,震源机制解显示典型的构造地震特征;震中区构造应力场理论应力轴σ₁方位角、倾伏角分别为88.0°和0.9°,σ₂方位角、倾伏角分别为178.2°和9.6°,σ₃方位角、倾伏角分别为352.5°和80.4°,这一结果与区域构造应力场一致。推断认为区域构造应力场触发了2019年松原M5.1地震活动,地震震源机制解的北西向节面与震中区附近的第二松花江断裂现今活动性质完全一致,认为第二松花断裂可能是松原M5.1地震的发震断层。     

Magnitude-frequency distribution in the European-Mediterranean earthquake regions

The European-Mediterranean earthquake catalogue from 1901 to 1985, which comprises uniformly determined magnitudes M_S and m_B(h ≥ 60 km) of 13300 events, was used in the study of cumulative magnitude-frequency relationships N_c(M) compiled for 75 earthquake regions and 25 larger provinces. In the whole magnitude range observed, the Gutenberg-Richter formula log N_c(M) = a − bM very rarely fits the cumulative (log N_c, M) distributions. The b-values of log-linear segments of N_c(M) vary regionally from b = 0.7 to b = 1.3; averaging of all values leads to (shallow events, M_S and ).

Most distributions pertain to the Mediterranean area (b = 0.86 from the graph for shallow events) and many of them indicate the existence of characteristic earthquakes in accordance with the theoretical single-fault model. Other observed shapes of N_c(M) can be explained by the superposition of populations of different M_max values or by the presence of swarm-type activity. The observed N_c(M) distributions depend very much on the delineation of earthquake regions i.e. on the number and dimension of seismoactive faults in the investigated region.

A premonitory enhancement of medium earthquake activity (M = 4.5–5.5) can be observed only very rarely.     

Strontium diffusion in sanidine and albite, and general comments on strontium diffusion in alkali feldspars

Strontium chemical diffusion has been measured in albite and sanidine under dry, 1 atm, and QFM buffered conditions. Strontium oxide-aluminosilicate powdered sources were used to introduce the diffusant and Rutherford Backscattering Spectroscopy (RBS) used to measure diffusion profiles. For the 1 atm experiments, the following Arrhenius relations were obtained:

Sanidine (Or₆₁), temperature range 725–1075°C, diffusion normal to (001): D=8.4 exp(−450±13 kJ mol⁻¹/RT) m²s⁻¹. Albite (Or₁), temperature range 675–1025°C, diffusion normal to (001): D=2.9 × exp(−224±11 kJ mol⁻¹/RT) m²s⁻¹.

The alkali feldspars in this and earlier work display a broad range of activation energies for Sr diffusion, which may be a consequence of the thermodynamic non-ideality of the alkali feldspar system and/or the mixed alkali effect.     

Basalt magma genesis and fractionation in collision- and extension-related provinces: A comparison between eastern, central and western Anatolia

Experimental studies of synthetic and natural basalt systems suggest that conditions of magma genesis and fractionation depend fundamentally on mantle temperatures and lithospheric stress fields. In general, compressional settings are more conducive to polybaric fractionation than extensional settings and in this regard, the Anatolian magmatic province offers a natural laboratory for comparing near-coeval basalt eruptions as a function of regional tectonics — compressional (collision-related) régimes dominating in eastern Anatolia and extensional tectonics characterizing a western province related to Aegean Sea opening. Projection of Plio-Quaternary basalt normative compositions from the Western Anatolia Extensional Province (WAEP), the Central Anatolian ‘Ova’ Province (CAOP), and Eastern Anatolia Compressional Province (EACP) are projected onto Ol–Ne–Cpx and Pl–Cpx–Ol planes in the simplified basalt system (Ne–Cpx–Ol–Qz), each showing distinctive liquid lines of descent. WAEP basalts are mostly constrained by low-pressure (< 0.5 GPa) cotectics while CAOP and EACP compositions conform to moderate and/or high-pressure (0.8–3.0 GPa) cotectics. Overall, a quasi-linear shift from moderate and/or high-pressure to low-pressure equilibria matches the westward transgression from compressional east Anatolia to the extensional west Anatolian–Aegean region. Comparison of their respective primary (mantle-equilibrated) magmas–simulated by normalizing their compositions to MgO = 15 wt.% (Mg-15)–with parameterized anhydrous and H₂O-undersaturated experimental melts suggests they segregated from spinel- to garnet-lherzolite mantle facies at pressures between c. 2 and 3 GPa (c. 70–100 km depth) under H₂O-undersaturated conditions. Interpolated potential temperatures (T_p) and lithospheric stretching factors (β) range as follows: (1) eastern Anatolian basalts associated with the Arabian foreland show T_p varying between 1250 and 1400 °C (except for the Karacalidag alkali basalts, south of the Bitlis–Zagros fracture zone, for which T_p ranges up to 1450 °C), for β values of 1.2–1.8. T_p values for central Anatolia (e.g. Sivas) range between 1300 and 1375 °C (except for Karapinar, Egrikuyu and Hasandag, which show < 1150 °C), and β values of 1.3–1.4. For western Anatolian basalts, T_p range mostly between 1250 and 1330 °C, except for a single value for Canakkale of 1400 °C and Kula sample showing T_p < 1200 °C, and β values of 1.3–2.0. Variation of these conditions is as great or greater than that between provinces, although there are clearly significant constraints on the inferred polybaric to low-pressure isobaric fractionation régimes. Covariation of total FeO, TiO₂, La/Yb, Ce/Sm, Zr/Y and Zr/Nb reflects small but significant differences in bulk composition and ambient melt fraction while the covariance of Ce/Sm and Sm/Yb is consistent with the segregation of primitive melts at the spinel- to garnet-lherzolite transition.     

Analysis of the Influence of Climate and Vegetation Change on Runoff in the Middle and Upper Reaches of the Pearl River Basin Based on Budyko Hypothesis*

Runoff, which is a key component in the hydrological cycle, is mainly controlled by climate factors and land-surface elements in non-humid regions. The impacts of climate and vegetation changes on runoff based on Budyko hypothesis in the middle and upper reaches of the Pearl River Basin was analyzed in this article. First, the temporal trend of variables in the study area during 1981-2013 was examined by using the Mann-Kendall trend test with trend-free pre-whitening. Second, the relationship of the parameter n in Fu's equation with factors of climate and vegetation coverage was built to reveal the time-variation process of n. Finally, the effects of climatic factors and vegetation coverage on runoff were assessed by analyzing the sensitivity of runoff to each variable. It is found that average temperature (T), maximum temperature (T_max) and minimum temperature (T_min) in the study area present an increasing trend while runoff (Q), precipitation (P), wind speed (u₂) and relative humid (RH) present decreasing trend. The parameter n in Fu's equation is significantly related to both climatic factors (including precipitation (P), average temperature (T), relative humid (RH), sunshine duration (S), wind speed (u₂)) and vegetation coverage index (NDVI). In terms of sensitivity of Runoff (Q) to the variation of each climatic factors and NDVI in the middle and upper reaches of the Pearl River Basin, precipitation (P) and NDVI have the highest sensitivity, followed by other climatic factors. Additionally, the precipitation (P) reduction is the main driving factor to the decline in runoff, while vegetation coverage is another important factor. In general, climate change affects runoff not only by changing the hydrological inputs (precipitation (P) and potential evaporation (PET) but also by altering the watershed characteristics as represented by the parameter n, while the impacts of vegetation coverage on runoff are exerted mainly through the alteration of the watershed characteristics.     

基于区域农业用水量的干旱重现期计算方法

为研究实际水利条件下农业干旱的发生规律,简化农业干旱事件的评估方法,提出基于区域农业用水量的干旱重现期计算方法。通过构建农业用水量距平百分率干旱指标W_A,在基于降雨量距平百分率干旱指标P_A识别干旱事件的基础上,提取W_A干旱指标下的干旱历时和干旱烈度特征变量,并根据以P_A为干旱指标的干旱烈度频率分布曲线F_S（x）和干旱历时频率分布曲线F_D（x）,运用Copula的简化方法计算基于W_A的干旱事件重现期T,最后结合基于P_A的干旱事件重现期T₀,回归分析出T与T₀间关系的计算公式。选取干旱灾害影响严重的亳州市为实证区域开展应用研究,计算得到1975-2007年各场干旱事件的T₀和T以及T₀与T的经验关系式。结果表明:T比T₀更合理地反映区域农业实际受旱状况,重现期T₀和T间存在高度的相关关系,采用T的回归方程可简化计算考虑区域实际抗旱能力下的干旱事件重现期,在区域防旱减灾实践中具有推广应用价值。     

Thermodynamic modelling of C–O–H fluids

H₂O, CO₂, CH₄, CO, H₂ and O₂ are the most important species in crustal fluids. The composition of these C–O–H fluids can be calculated if the pressure, temperature, carbon activity, and either the oxygen fugacity or the atomic H/O ratio of the fluid is known. The calculation methods are discussed and calculation results are illustrated with isobaric T–X_i, P–T, and isobaric–isothermal ternary C–O–H diagrams. Fluid inclusion compositions, in particular, the X_CO2/(X_CO2+X_CH4) ratio, can be used for C–O–H model calculations. However, care should be taken about possible post-entrapment changes, which may have modified the chemical composition of the fluid inclusion.     

Review of the systematics of CO2–H2O fluid inclusions

Aqueous solutions that contain volatile (gas) components are one of the most important types of fluid in the Earth's crust. The record that such fluids have left in the form of fluid inclusions in minerals provides a wealth of insight into the geochemical and petrologic processes in which the fluids participated. This article reviews the systematics of CO₂–H₂O fluid inclusions as a starting point for interpreting the chemically more complex systems. The phase relations of the binary are described with respect to a qualitative P–T–X model, and isoplethic–isochoric paths through this model are used to explain the equilibrium and non-equilibrium behaviour of fluid inclusions during microthermometric heating and cooling. The P–T–X framework is then used to discuss the various modes of fluid inclusion entrapment, and how the resulting assemblage textures can be used to interpret the P–T conditions, phase states, and evolution paths of the parent solutions. Finally, quantitative methods are reviewed by which bulk molar volume and composition of CO₂–H₂O fluid inclusions can be determined from microthermometric observations of phase transitions.     

Provenance of the Neoproterozoic Maricá Formation (Sul-rio-grandense Shield, Southern Brazil): Petrographic and Sm–Nd isotopic constraints

This paper reports the integrated application of petrographic and Sm–Nd isotopic analyses for studying the provenance of the Neoproterozoic Maricá Formation, southern Brazil. This unit encompasses sedimentary rocks of fluvial and marine affiliations. In the lower fluvial succession, sandstones plot in the “craton interior” and “transitional continental” fields of the QFL diagram. Chemical weathering probably caused the decrease of the ¹⁴⁷Sm/¹⁴⁴Nd ratios to 0.0826 and 0.0960, consequently lowering originally > 2.0 Ga T_DM ages to 1.76 and 1.81 Ga. ¹⁴³Nd/¹⁴⁴Nd ratios are also low (0.511521 to 0.511633), corresponding to negative ε_Nd present-day values (− 21.8 and − 19.6). In the intermediate marine succession, sandstones plot in the “dissected arc” field, reflecting the input of andesitic clasts. Siltstones and shales reveal low ¹⁴³Nd/¹⁴⁴Nd ratios (0.511429 to 0.511710), ε_Nd values of − 18.1 and − 23.6, and T_DM ages of 2.16 and 2.37 Ga. Sandstones of the upper fluvial succession have “dissected arc” and “recycled orogen” provenance. ¹⁴³Nd/¹⁴⁴Nd isotopic ratios are also relatively low, from 0.511487 to 0.511560, corresponding to ε_Nd values of − 22.4 and − 21.0 and T_DM of 2.07 Ga. A uniform granite–gneissic basement block of Paleoproterozoic age, with subordinate volcanic rocks, is suggested as the main sediment source of the Maricá Formation.     

Evolution of fluids associated with metasedimentary sequences from Chaves (North Portugal)

In order to identify and characterise fluids associated with metamorphic rocks from the Chaves region (North Portugal), fluid inclusions were studied in quartz veinlets, concordant with the main foliation, in graphitic-rich and nongraphitic-rich lithologies from areas with distinct metamorphic grade. The study indicates multiple fluid circulation events with a variety of compositions, broadly within the C–H–O–N–salt system. Primary fluid inclusions in quartz contain low salinity aqueous–carbonic, H₂O–CH₄–N₂–NaCl fluids that were trapped near the peak of regional metamorphism, which occurred during or immediately after D2. The calculated P–T conditions for the western area of Chaves (CW) is P=300–350 MPa and T500 °C, and for the eastern area (CE), P=200–250 MPa and T=400–450 °C. A first generation of secondary fluid inclusions is restricted to discrete cracks at the grain boundaries of quartz and consists of low salinity aqueous–carbonic, H₂O–CO₂–CH₄–N₂–NaCl fluids. P–T conditions from the fluid inclusions indicate that they were trapped during a thermal event, probably related with the emplacement of the two-mica granites.

A second generation of secondary inclusions occurs in intergranular fractures and is characterised by two types of aqueous inclusions. One type is a low salinity, H₂O–NaCl fluid and the second consists of a high salinity, H₂O–NaCl–CaCl₂ fluid. These fluid inclusions are not related to the metamorphic process and have been trapped after D3 at relatively low P (hydrostatic)–T conditions (P<100 MPa and T<300 °C).

Both the early H₂O–CH₄–N₂–NaCl fluids in quartz from the graphitic-rich lithologies and the later H₂O–CO₂–CH₄–N₂–NaCl carbonic fluid in quartz from graphitic-rich and nongraphitic-rich lithologies seem to have a common origin and evolution. They have low salinity, probably resulting from connate waters that were diluted by the water released from mineral dehydration during metamorphism. Their main component is water, but the early H₂O–CH₄–N₂–NaCl fluids are enriched in CH₄ due to interaction with the C-rich host rocks.

From the early H₂O–CH₄–N₂–NaCl to the later aqueous–carbonic H₂O–CO₂–CH₄–N₂–NaCl fluids, there is an enrichment in CO₂ that is more significant for the fluids associated with nongraphitic-rich lithologies.

The aqueous–carbonic fluids, enriched in H₂O and CH₄, are primarily associated with graphitic-rich lithologies. However, the aqueous–carbonic CO₂-rich fluids were found in both graphitic and nongraphitic-rich units from both the CW and CE studied areas, which are of medium and low metamorphic grade, respectively.     

Petrogenesis of the Maowu pyroxenite–eclogite body from the UHP metamorphic terrane of Dabieshan: chemical and isotopic constraints

The Maowu eclogite–pyroxenite body is a small (250×50 m) layered intrusion that occurs in the ultra-high-pressure (UHP) metamorphic terrane of Dabieshan, China. Like the adjacent Bixiling complex, the Maowu intrusion was initially emplaced at a crustal level, then subducted along with the country gneisses to mantle depths and underwent UHP metamorphism during the collision of the North and South China Blocks in the Triassic. This paper presents the results of a geochemical and isotopic investigation on the metamorphosed Maowu body. The Maowu intrusion has undergone open system chemical and isotopic behavior three times. Early crustal contamination during magmatic differentiation is manifested by high initial ⁸⁷Sr/⁸⁶Sr ratios (0.707–0.708) and inhomogeneous negative _Nd(T) values of −3 to −10 at 500 Ma (probable protolith age). Post-magmatic and pre-UHP metamorphic metasomatism is indicated by sinusoidal REE patterns of garnet orthopyroxenites, lack of whole-rock (WR) Sm–Nd isochronal relationship, low δ¹⁸O values and an extreme enrichment of Th and REE in a clinopyroxenite. Finally, K and Rb depletion during UHP metamorphism is deduced from the high initial ⁸⁷Sr/⁸⁶Sr ratios unsupported by in situ Rb/Sr ratios. Laser ICP-MS spot analyses on mineral grains show that (1) Grt and Cpx attained chemical equilibrium during UHP metamorphism, (2) Cpx/Grt partition coefficients for REE correlate with Ca, and (3) LREE abundances in whole rocks are not balanced by that of the principal phases (Grt and Cpx), implying that the presence of LREE-rich accessory phases, such as monazite and apatite, is required to account for the REE budget.

Sm–Nd isotope analyses of minerals yielded three internal isochrons with ages of 221±5 Ma and (T)=−5.4 for an eclogite, 231±16 Ma and (T)=−6.2 for a garnet websterite, and 236±19 Ma and (T)=−6.9 for a garnet clinopyroxenite. The Cpx/Grt chemical equilibrium and the consistent mineral isochron ages indicate that the metasomatic processes mentioned above must have occurred prior to the UHP metamorphism. These Sm–Nd ages agree with published zircon and monazite U–Pb ages and constrain the time of UHP metamorphism to 220–236 Ma. The Maowu and Bixiling layered intrusions are similar in their in situ tectonic relationship with their country gneisses, but the two bodies are distinguished by their magma-chamber processes. The Bixiling magmas were contaminated by the lower crust, whereas the Maowu magmas were contaminated by the upper crustal rocks during their emplacement and differentiation. The two complexes represent two distinct suites of magmatic rocks, which have resided in the continental crust for about 300–400 Ma before their ultimate subduction to mantle depths, UHP metamorphism and return to the crustal level.     

A correlation for permanent earthquake-induced deformation of earth embankments

Seismic design of earth dams and embankments is mainly controlled by the permanent deformation that may be induced directly or indirectly because of the design earthquake. Simple design charts, such as those developed by Yegian et al. [Yegian, M.K., Marciano, E.A., and Ghahraman, V.G. 1991. Earthquake-induced permanent deformations: Probabilistic approach. J. Geotech. Eng., 117, 35–50.], Hynes-Griffin and Franklin [Hynes-Griffin, M.E., and Franklin, A.G. 1984. Rationalizing the Seismic Coefficient Method. Miscellaneous Paper GL-84-13, US Army Corps of Engineers Waterways Experiment Station, Vicksburg, Mississippi.], Makdisi and Seed [Makdisi, F.I., and Seed, H.B. 1978. Simplified procedure for estimating dam and embankment earthquake-induced deformations. J. Geotech. Eng., 104, 849–867.], and Sarma [Sarma, S.K. 1975. Seismic stability of earth dams and embankments. Géotechnique, 25, 743–761.], are often used to obtain a preliminary estimate of the permanent, earthquake-induced deformation of earth dams and embankments. Comparison of permanent deformations estimated from these procedures with observations from 122 published case histories on performance of earth dams and embankments during past earthquakes indicate that the estimated permanent earthquake-induced deformations were, in general, smaller than the observed deformations. However, the observed permanent deformation, D_avg, was found to relate to the ratio of yield acceleration, a_y, and the peak horizontal ground acceleration, a_max. The scatter in the observational data, upon which the D_avg–a_y/a_max relationship is based, can be partly explained by the variations in the ratio of the fundamental (elastic) period of the earth structure, T_D, and the predominant period of the earthquake ground motion, T_p, and the magnitude of the earthquake, M_W.     

P–T paths reconstruction of a collisional event: The example of the Thiviers-Payzac Unit in the Variscan French Massif Central

Because of late metamorphic and tectonic overprints, the reconstruction of prograde parts of P–T paths is often difficult. In the SW Variscan French Massif Central, the Thiviers-Payzac Unit (TPU) is the uppermost allochthon emplaced above underlying units. The TPU experienced a Barrovian metamorphism coeval with a top-to-the-NW ductile shearing (D2 event) in Early Carboniferous times (ca. 360–350 Ma). The tectonic setting of the D2 event, compression or synconvergence extension, remains unclear. Using the THERMOCALC software and the model system MnNCKFMASH, the peak P–T conditions are estimated from garnet rims and matrix minerals and the prograde evolution is deduced from garnet core compositions. The combination of these two approaches demonstrates that the TPU experienced pressure and temperature increases before reaching peak conditions at 6.6–9.0 +/− 1.2 kbar and 615–655 +/− 35 °C. This kind of P–T path shows that the regional D2 event corresponds to crustal thickening.     

Formation of Al-rich titanite (CaTiSiO4O–CaAlSiO4OH) reaction rims on ilmenite in metamorphic rocks as a function of fH2O and fO2

Reaction rims of titanite on ilmenite are described in samples from four terranes of amphibolite-facies metapelites and amphibolites namely the Tamil Nadu area, southern India; the Val Strona area of the Ivrea-Verbano Zone, northern Italy, the Bamble Sector, southern Norway, and the northwestern Austroalpine Ötztal Complex. The titanite rims, and hence the stability of titanite (CaTiSiO₄O) and Al–OH titanite, i.e. vuaganatite (hypothetical end-member CaAlSiO₄OH), are discussed in the light of fH₂O- and fO₂-buffered equilibria involving clinopyroxene, amphibole, biotite, ilmenite, magnetite, and quartz in the systems CaO–FeO/Fe₂O₃–TiO₂–SiO₂–H₂O–O₂ (CFTSH) and CaO–FeO/Fe₂O₃–Al₂O₃–SiO₂–H₂O–O₂ (CFASH) present in each of the examples. Textural evidence suggests that titanite reaction rims on ilmenite in rocks from Tamil Nadu, Val Strona, and the Bamble Sector originated most likely due to hydration reactions such as clinopyroxene + ilmenite + quartz + H₂O = amphibole + titanite and oxidation reactions such as amphibole + ilmenite + O₂ = titanite + magnetite + quartz + H₂O during amphibolite-facies metamorphism, or, as in the case of the Ötztal Complex, during a subsequent greenschist-facies overprint. Overstepping of these reactions requires fH₂O and fO₂ to be high for titanite formation, which is also in accordance with equilibria involving Al–OH titanite. This study shows that, in addition to P, T, bulk–rock composition and composition of the coexisting fluid, fO₂ and fH₂O also play an important role in the formation of Al-bearing titanite during amphibolite- and greenschist-facies metamorphism.     

Ultra high temperature-metamorphism and its significance in the Central Indian Tectonic Zone

In the present study from the southern margin of the Central Indian Tectonic Zone, it is demonstrated how the metamorphic P–T path of ultrahigh-temperature granulite terranes can be reconstructed using the metamorphic transition in corundum granulites from early biotite melting to later FMAS solid–solid reaction. The extreme metamorphism in these rocks caused two-stage biotite melting, resulting in initial porphyroblastic garnet₁ and later sapphirine–spinel₁ incongruent solid mineral assemblages. During this process, the leucocratic and melanocratic layers in the corundum granulites evolved from an initial silica-oversaturated to a later silica-undersaturated domain. In the melanocratic layer, this allowed localized concentration of sapphirine-spinel₁ and residual sillimanite₁, producing an extremely restitic assemblage, at the culmination of peak metamorphism, BM₁. BM₁ is constrained at  1000 °C at relatively deep crustal levels (P  9 kbar) from the stability of ferroaugite in a co-metamorphosed Iron Formation granulite. During subsequent metamorphism (BM₂), the reaction path and history in the corundum granulites shifted to the restitic domain allowing reacting sapphirine, spinel₁ and sillimanite to produce coronal garnet₂–corundum assemblage via a FMAS univariant reaction. In the final stages of reaction history, biotite₂–sillimanite₂–spinel₂ assemblage was produced after garnet₂–corundum due to localized melt–crystal interaction. The metamorphic sequence, when interpreted with the help of a newly constructed, qualitative KFMASH petrogenetic grid, reveals successive stages of heating, increasing pressure and cooling around the KFMASH invariant point, [Opx,Crd], which is consistent with a counterclockwise metamorphic P–T path. The near isobaric nature of post-peak cooling (ΔT  250–300 °C) is also evident from multistage pyroxene exsolution and by the appearance of lamellar and coronal garnets in the Iron Formation granulites. This study provides the first tight constraint for ultrahigh-T metamorphism along a counter clockwise P–T trajectory in the Central Indian Tectonic zone, and has important bearing for terrane correlations in this part of East Gondwanaland. In addition, the new KFMASH grid allows evaluation of metamorphic phase relations in ultrahigh-T, corundum-bearing and corundum-absent aluminous granulites.     

Reactions and textures in grossular–wollastonite–scapolite calc–silicate granulites from Maligawila, Sri Lanka: evidence for high-temperature isobaric cooling in the meta-sediments of the Highland Complex

Grossular–wollastonite–scapolite calc–silicate granulites from Maligawila in the Buttala klippe, which form part of the overthrusted rocks of the Highland Complex of Sri Lanka, preserve a number of spectacular coronas and replacement textures that could be effectively used to infer their P–T–fluid history. These textures include coronas of garnet, garnet–quartz, and garnet–quartz–calcite at the grain boundaries of wollastonite, scapolite, and calcite as well as calcite–plagioclase and calcite–quartz symplectites or finer grains after scapolite and wollastonite respectively. Other textures include a double rind of coronal scapolite and coronal garnet between matrix garnet and calcite. The reactions that produced these coronas and replacement textures, except those involving clinopyroxene, are modelled in the CaO–Al₂O₃–SiO₂–CO₂ system using the reduced activities. Calculated examples of T–X_CO2 and P–X_CO2 projections indicate that the peak metamorphic temperature of about 900–875 °C at a pressure of 9 kbar and the peak metamorphic fluid composition is constrained to be low in X_CO2 (0.1<X_CO2<0.30). Interpretation of the textural features on the basis of the partial grids revealed that the calc–silicate granulites underwent high-temperature isobaric cooling, from about 900–875 °C to a temperature below 675 °C, following the peak metamorphism. The late-stage cooling was accompanied by an influx of hydrous fluids. The calc–silicate granulites provide evidence for high-temperature isobaric cooling in the meta-sediments of the Highland Complex, earlier considered by some workers to be confined exclusively to the meta-igneous rocks. The coronal scapolite may have formed under open-system metasomatism.     

Magma–host interactions during differentiation and emplacement of a shallow-level, zoned granitic pluton (Tarçouate pluton, Morocco): implications for magma emplacement

The Tarçouate pluton (Anti-Atlas, Morocco) is an inversely zoned laccolith emplaced 583 Ma ago into low-grade metasediments, with the following succession: leucocratic granites, biotite–granodiorites (±monzodiorites), hornblende–granodiorites (±monzodiorites) and monzodiorites syn-plutonic dykes. These rocks form two distinct, chemically coherent, units:

(1) A main unit consists of layered (572<59 wt.%) and homogeneous (632<67%) hornblende–granodiorites, biotite–granodiorites (672<72%) and aplites (702<76%). All these rocks are metaluminous to peraluminous and display fractionated HREE depleted patterns (La/Yb_N=14–61; Yb_N=0.7–6.8). Initial ⁸⁷Sr/⁸⁶Sr ratios (0.7072 to 0.7080) increase, whereas _Nd(t) values (−1.7 to −2.8) decrease from the hornblende– to the biotite–granodiorites. Monzodiorites occur as mafic microgranular enclaves or syn-plutonic dykes.

(2) A subordinate unit consists of leucocratic, distinctly peraluminous, muscovite-bearing granites (722<75%) occurring at the northern edge of the pluton and as dykes in the surrounding schists towards the top of the pluton. These rocks are free of monzodioritic enclaves. They display less fractionated patterns with higher HREE contents (La/Yb_N=2–19; Yb_N=11–18), a distinct _Nd(t) value (−11.8) and a ⁸⁷Sr/⁸⁶Sr initial ratio (0.7480) within those of the surrounding schists (0.7393–0.7819).

Magma–host interactions are closely related to differentiation and occurred at different levels, but mainly before emplacement. Field relationships and petrogenetic modelling show that the bt–granodiorites formed at levels deeper than the level of emplacement, by fractional crystallisation (0.65

These data preclude any significant material transfer process for the emplacement of the Tarçouate pluton, but rather suggest assembly of successive pulses of variably differentiated crystal-poor magmas. These shallow level granitic plutons can be considered as an end-member of magma emplacement with minimum interactions with the country rocks.     

Partial melting of crustal rocks

Shales and graywackes were first metamorphosed at 650°C and then partially melted at 700 and 750°C at 2, 4, 6, and 8 kilobars in the presence of 0.75m NaCl−0.45m KCl and 0.225m CaCl₂−0.750m NaCl solutions. In experiments with shales,KK+Na ratio in the decreases with increasing pressure at 650 and 700°C; however, at 750°C this ratio is equal to 0.5 at all pressures investigated. This suggests that melts at 700°C and at 2 to 8 kilobars pressure may be affected metasomatically whereas melts at 700°C and in the same pressure range will not. Melt composition produced in the shale-KClNaCl experiments is granite at 2, 4 and 6 kilobars pressure, whereas the melt compositions in the shale-CaCl₂NaCl experiments range from quartz monzonite (2–5 kilobars) to granodiorite (above 5 kilobars). Experiments with graywacke-KClNaCl produced melts of trondhjemite composition at 2, 2.5, 4, 6, and 7.5 kilobars.

These results indicate that partial melting of crustal rocks such as metamorphosed shales and graywackes in the deeper parts of the crust can produce large volumes of granitic magmas ranging in composition from true granite to trondhjemite to quartz monzonite and granodiorite.