Discovery of Late Cretaceous garnet two-pyroxene granulite in the southern Lhasa terrane, Tibet and its tectonic significances

The Lhasa terrane, located between the Bangonghu-Nujiang suture zone and the Indus-Yalung Tsangpo suture zone in the southern Tibetan Plateau, was considered previously as a Precambrian continental block. Mesozoic and Cenozoic tectonic evolution of the Lhasa terrane is closely related to the subduction of the Tethys ocean and the collision between the Indian and European continents; so it is one of the keys to reveal the formation and evolution of the Tibetan plateau. The garnet two-pyroxene granulite which was found at the Nyingtri rock group of the southeastern Lhasa terrene consists of garnet, clinopyroxene, orthopyroxene, labradorite, Ti-rich amphibolite and biotite, with a chemical composition of mafic rock. The metamorphic conditions were estimated to be at T = 747 similar to 834 degrees C and P = 0.90 similar to 1.35GPa, suggesting a formation depth of 45km. The zircon U-Pb dating for the garnet amphibolite and marble associated with the granulite give a metamorphic age of 85 similar to 90Ma. This granulite-facies metamorphic event together with a contemporaneous magmatism demonstrated that the southern Lhasa terrane has undergone an Andean-type orogeny at Late Mesozoic time.

华北克拉通孔玆岩系中石榴子石花岗岩的成因：熔体与残留相混合机制的研究

The most widespread rock associations in the Western Block of North China Craton are khondalites distributed mainly in Jining, Liangcheng and Datong. A large quantitiy of garnet-bearing granites are contained in the khondalites. A great deal of research has been carried out on them by previous researchers. Studies of these garnet-bearing granites consist essentially of structural characteristics, petrography and geochemistry, and finally geochronological determinations. Summing up these researches, it will not be difficult to see that all of these authors have regarded these large numbers of garnets (up to 20%) contained in granites as crystallized products from magmas, but they have not proved this from petrological perspective. Theoretically, there are possibly three kinds of petrogenesis as to these garnets. The first one is that they have been transferred to the granites from khondalites by melt when anatexis happened to khondalites, and they, in essence, are residual metamorphic garnets; The second one is that when the khondalites were being melted, these garnets were produced from biotite dehydration melting, and the newly formed garnets intruded together with the melt and eventually molded the garnet-bearing granites. Garnets of this possible kind either showed independent crystals, or garnets from khondalites took place secondary growth under favorable temperature and pressure conditions for their crystallization; The last possibility is that these garnets were crystallized from magmas in which suitable pressure, temperature and composition were available. These garnets, generally, should be fine-grained. The aim of this study is, through examining the mineral chemistry of the garnets and the whole rock chemistry, to ascertain under which kind of mechanism, in the world, did these garnets form? Besides, we try to calculate the temperatures under which khondalites began melting and reactions of the garnets and the cooled melts happened by garnet-biotite thermometry. The whole rock chemistry analyses of the garnet-bearing granites tell us that all the samples are strongly peraluminous (A/CNK greater than 1.1) on the A/NK vs. A/CNK plot. On the SiO2-K2O plot, the granites are mainly constrained to be high-K calc-alkaline and calc-alkaline series, consistent with previous researches. On the ACF（（Al2O3-Na2O-K2O）-FeO(T)-CaO） discrimination plot, all the six garnet-bearing granite samples drop into the area of S-type granites. The relationship between CaO/Na2O and SiO2 shows that the overwhelming majority of garnet-bearing granites have a CaO/Na2O value over 0.3, revealing that they probably come from metagreywacke precursors or mediate-felsic orthogeneisses compositionally similar to them. Detailed EPMA analyses conducted on the garnets contained in the garnet-bearing granites show that all the garnets are dominated by almandine and pyrope, which occupy 92-96% (Weight Percentage) of each garnet analyzed, typical of granulite facies. Their chemical composition is entirely different from those crystallized in magmas, but extremely similar to those of typical granulite facies metapelites in khondalites and typical granulites, indicating all the garnets to be metamorphogenic. In addition, REEs distribution patterns of the garnets are totally different from typical biotite granites and peraluminous granites. In other words, both LREE and HREE of our garnets are evidently lower than those from these two kinds of rocks. Moreover, compared to the REE pattern of the garnets from typical amphibolites, LREE content of our garnets is obviously higher and HREE content is a little lower. However, REE patterns of our garnets are completely in harmony with those of garnets from typical granulites. So, the REE patterns of garnets, again, prove that all the garnets we studied are metamorphogenic. Biotites appear in two forms, being as inclusions in the garnet and as selvages immediately adjacent to the garnet, respectively. Two reactions and their corresponding temperatures, with the help of petrography and Garnet-Biotite geothermometers, could be obtained, which are Bt+ Pl+ Qtz→Kfs+ Opx+ Grt+ melt as positive reaction and Kfs+ Grt+ melt→Bt+ Pl+ Qtz as reverse reaction, respectively. Summing up the discussion above, we declare that the garnet-bearing granites distributed in the Western Block of North China Craton are the mixture of melts and restites resulted from biotite dehydration melting. The garnets contained in the restites are the products from biotite dehydration melting and restites from the khondalites, respectively.

北大别黄土岭麻粒岩变质演化：对大别造山带演化进程的指示意义

Granulites from Huangtuling in the North Dabie metamorphic core complex in eastern China preserve rare mineralogical and mineral chemical evidence for multistage metamorphism related to Paleoproterozoic metamorphic processes, Triassic continental subduction-collision and Cretaceous collapse of the Dabieshan Orogen. Six stages of metamorphism are established, based on detailed mineralogical and petrological studies: (I) amphibolite facies (6.3–7.0 kbar, 520–550 °C); (II) high-pressure/high-temperature granulite facies (12–15.5 kbar, 920–980 °C); (III) cooling and decompression (4.8–6.0 kbar, 630-700 °C); (IV) medium-pressure granulite facies (7.7–9.0 kbar, 690–790 °C); (V) low-pressure/high-temperature granulite facies (4.0–4.7 kbar, 860–920 °C); (VI) retrograde greenschist facies overprint (1–2 kbar, 340–370 °C). The P–T history derived in this study and existing geochronological data indicate that the Huangtuling granulite records two cycles of orogenic crustal thickening events. The earlier three stages of metamorphism define a clockwise P–T path, implying crustal thickening and thinning events, possibly related to the assembly and breakup of the Columbia Supercontinent ca. 2000 Ma. Stage IV metamorphism indicates another crustal thickening event, which is attributed to the Triassic subduction/collision between the Yangtze and Sino–Korean Cratons. The dry lower crustal granulite persisted metastable during the Triassic subduction/collision due to lack of hydrous fluid and deformation. Stage V metamorphism records the Cretaceous collapse of the Dabieshan Orogen，possibly due to asthenosphere upwelling or removal of the lithospheric mantle resulting in heating of the granulite and partial melting of the North Dabie metamorphic core complex. Comparison of the Huangtuling granulite in North Dabie and the high-pressure (HP)–ultrahigh-pressure (UHP) metamorphic rocks in South Dabie indicates that the subducted upper (South Dabie) and lower (North Dabie) continental crusts underwent contrasting tectonometamorphic evolution during continental subduction–collision and orogenic collapse. High-pressure granulites are generally characterized by the absence of orthopyroxene. However, the Huangtuling felsic granulite rarely preserves the high-pressure granulite facies assemblage of garnet + orthopyroxene + biotite + plagioclase + K-feldspar + quartz. To investigate the effects of bulk rock composition on the stability of orthopyroxene-bearing, high-pressure granulite facies assemblages in the NCKFMASHTO (Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–Fe2O3) system, we constructed a series of P–T–X pseudosections based on the melt-reintegrated composition of the Huangtuling felsic high-pressure granulite. Our calculations demonstrate that the orthopyroxene-bearing, high-pressure granulite facies assemblages are restricted to low XAl [Al2O3/(Na2O + CaO + K2O + FeO + MgO + Al2O3) < 0.35, mole proportion] or high XMg [MgO/(MgO + FeO) > 0.85] felsic–metapelitic rock types. Our study also reveals that the XAl values in the residual felsic–metapelitic, high-pressure granulites could be significantly reduced by a high proportion of melt loss. We suggest that orthopyroxene-bearing high-pressure granulites occur in residual overthickened crustal basement under continental subduction–collision zones and arc–continent collision belts.

云南岩石圈流变学及多震层孕震环境

The Study on rheology of the lithosphere and the environments of the seismogenic layer is currently the basic project of the international earthquake research. Yunnan is the ideal place for studying this project. Through the multi-disciplinary comprehensive study of petrology, geophysics, seismo-geology, rock mechanics, etc., the depth-strength profiles of the lithosphere have been firstly constructed, and the seismogenic layer and its geophysical and tectonic environments in Yunnan have been systematically expounded in this paper. The related results achieved are of the important significances for further understanding the mechanism of strong earthquake generation, dividing the potential foci and exposing recent geodynamical processes in Yunnan. Through the comprehensive contrast of the metamorphic rocks in early and middle Proterozoic outcropping on the surface, DSS data and experimental data of rock seismic velocity under high temperature and high pressure, the petrological structure of the crust and upper mantle has been studied on Yunnan: the upper, middle and lower crust is composed of the metamorphic rocks of greenschist, amphibolite and granulite facies, respectively or granitoids, diorites and gabbros, respectively, and the upper mantle composed of the peridotites. Through the contrast studies of the heat flow and epicenters of the strong earthquakes, the distribution of the geotemperature and the data of focal depth, the relationship of between seismicity and geothermal structure of the lithosphere in Yunnan has been studied: the strong earthquakes with magnitude M ≥ 6.0 mainly take place at the geothermal gradient zone, and the seismic foci densely distribute between 200～500 ℃ isogeotherms. On the basis of studies of the rock properties and constituents of the crust and upper mantle and geothermal structure of the lithosphere, the structure of the rheological stratification of the lithosphere has been studied, and the corresponding depth-strength profiles have been constructed in Yunnan. The lithosphere in majority region of Yunnan has the structure of the rheological stratification, i.e. the brittle regime in the upper crust or upper part of the upper crust, ductile regime in the middle crust or lower part of the upper crust to middle crust, ductile regime in the lower crust and ductile regime in the subcrustal lithosphere. The rheological stratification has the quite marked lateral variations in the various tectonic units. The distributions of the seismogenic layer have been determined by using the high accurate data of focal depth. Through the contrast of the petrological structure, the structure of seismic velocity, electric structure, geotemperature structure, and rheological structure and the study of the focal mechanism in the seismogenic layer, the geophysical environments of the seismogenic layer in Yunnan have been studied. The seismogenic layer in Yunnan is located at the depths of 3 ～ 20 km; the rocks in the seismogenic layer are composed of the metamorphic rocks of greenschist to amphibolite facies (or granites to diorites); the seismogenic layer and its internal focal regions of strong earthquakes have the structure of medium properties with the relatively high seismic velocity, high density and high resistivity; there exists the intracrustal low seismic velocity and high conductivity layer bellow the seismogenic layer, the geotemperature is generally 100～500 ℃ in the depth range in which the seismogenic layer is located. The horizontal stress field predominates in the seismogenic layer, the seismogenic layer corresponds to the brittle regime of the upper crust or brittle regime of the upper crust to semibrittle regime of the middle crust. The formation of the seismogenic layer, preparedness and occurrence of the strong earthquakes is the result of the comprehensive actions of the source fault, rock constituent, structure of the medium properties, distribution of the geotemperature, rheological structure of the seismogenic layer and its external environments. Through the study of the structure, active nature, slip rate, segmentation of the active faults, and seismogenic faults, the tectonic environments of the seismogenic layer in Yunnan have been studied. The source faults of the seismogenic layer in Yunnan are mainly A-type ones and embody mainly the strike slip faults with high dip angle. the source faults are the right-lateral strike slip ones with NW-NNW trend and left-lateral strike slip ones with NE-NEE trend in Southwestern Yunnan, the right-lateral strike slip ones with NNW trend and left-lateral strike slip ones with NNE trend (partially normal ones) in Northwestern Yunnan, the right-lateral strike slip ones with NWW trend in Central Yunnan and left-lateral strike slip ones with NW-NNW trend in Eastern Yunnan. Taking Lijiang earthquake with Ms = 7.0 for example. The generating environments of the strong earthquake and seismogenic mechanical mechanism have been studied: the source region of the strong earthquake has the media structure with the relatively high seismic velocity and high resistivity, there exists the intracrustal low velocity and high conductivity layer bellow it and the strong earthquakes occur near the transitional zone of the crustal brittle to ductile deformation. These characteristics are the generality of the generating environments of strong earthquakes. However, the specific seismogenic tectonic environments and action of the stress field of the seismic source in the various regions, correspondingly constrains the dislocation and rupture mechanical mechanism of source fault of strong earthquake.

内蒙古东部喀喇沁地区早中生代底侵作用

Mafic granulite xenoliths have been extensively concerned over the recent years because they are critical not only to studies of composition and evolution of the deep parts of continental crust but to understanding of the crust-mantle interaction. Detailed petrology, geochemistry and isotope geochronology of the Early Mesozoic mafic-ultramafic cumulate xenoliths and mafic granulite xenoliths and their host diorites from Harqin area, eastern Inner-Mongolia have been studied here. Systematic Rb-Sr isochron, ~(40)Ar-~(39)Ar and K-Ar datings for mafic-ultramafic cumulate xenoliths give ages ranging from 237Ma to 221Ma. Geochemical research and forming temperature and pressure estimates suggest that cumulates are products of the Early Mesozoic mantle-derived magmatic underplating and they formed in the magmatic ponds at the lowermost of the continental crust and are later enclaved by the dioritic magma. Detailed study on the first-discovered mafic granulite xenoliths reveals that their modal composition, mineral chemistry and metamorphic P-T conditions are all different from those of the Precambrian granulite exposed on the earth surface of the North China craton. High-resolution zircon U-Pb dating suggests that the granulite facies metamorphism may take place in 253 ～ 236Ma. Hypersthene single mineral K-Ar dating gives an age of 229Ma, which is believed to represent a cooling age of the granulite. As the host rock of the cumulate and granulite xenoliths, diorites intruded into Archean metamorphic rocks and Permian granite. They are mainly composed of grandodiorite, tonalite and monzogranite and show metaluminous and calc-alkaline features. Whole rock and single mineral K-Ar dating yields age of 221 ～ 223Ma, suggesting a rapid uplift in the forming process of the diorites. Detailed field investigation and geochemical characteristics indicate that these diorites with different rock types are comagmatic rocks, and they have no genetic correlation with cumulate and granulite xenoliths. Geochemical model simulating demonstrates that these diorites in different lithologies are products of highly partial melting of Archean amphibolite. It is considered that the Early Mesozoic underplating induced the intrusion of diorites, and it reflects an extensional geotectonic setting. Compression wave velocity V_P have been measured on 10 representative rock samples from the Early Mesozoic granulite and mafic-ultramafic cumulate xenoliths population as an aid to interpret in-situ seismic velocity data and investigating velocity variation with depth in a mafic lower crust. The experiments have been carried out at constant confining pressures up to 1000MPa and temperatures ranging from 20 ℃ to around 1300 ℃, using the ultrasonic transmission technique. After corrections for estimated in situ crustal pressures and temperatures, elastic wave velocities range from 6.5 ～ 7.4 km s~(-1). On the basis of these experimental data, the Early-Mesozoic continental compression velocity profile has also been reestablished and compared with those of the present and of the different tectonic environments in the world. The result shows that it is similar to the velocity structure of the extensional tectonic area, providing new constraints on the Early Mesozoic continental structure and tectonic evolution of the North-China craton. Combining with some newly advancements about the regional geology, the thesis further proposes some constraints on the Mesozoic geotectonic evolution history, especially the features of deep geology of the North China craton.

牡丹江碰撞带的组成及演化

Jiamusi Massif is an important tectonic unit in Northeast China. It’s significant for understanding the evolution of Paleo-Asian Ocean and reconstruction of the tectonic framework of Northeast China. Mudanjiang area is located in the southern margin of Jiamusi Massif and is the key to understand the evolution of Jiamusi Massif. However, the detailed geological research for Mudanjiang area has long been deficient in many important problems, such as the tectonic components of the Mudanjiang collision zone (MCZ), the age of collisional complexes and the scenario of tectonic evolution. Based on the lithology, geochemistry and the SHRIMP zircon U-Pb geochronology in Mudanjiang area, our new data and results come to some constraints for the tectonic reconstruction of MCZ as follows: 1） It is identified that the former suggestion, which the so-called “Heilongjiang Group” in Mudanjiang area is the vestige of oceanic crust, is correct. The oceanic relics represent the Neo-Proterozoic-Early Paleozoic oceanic basins based on the SHRIMP zircon U-Pb geochronology. 2） One sheet of gabbroic complex with oceanic island-type geochemical signature was discovered by this work in Mudanjiang area. 3） It is verified that the Proterozoic concordant U-Pb ages of the migmatites developed along the southern margin of Jiamusi massif, which represent the events of magmatic intrusion, as the direct evidence for the existence of the Proterozoic crystalline basements of the Jiamusi Massif. Based on geochronology, we suggest that the migmatization and coeval S-type granite magmatism of the southern margin of Jiamusi Massif took place about 490Ma. 4） The island arc complex has been found in the Heilongjiang Group, and the oceanic relics was found distributing on both sides, as provided important constraint for the tectonic reconstruction of the MCZ. 5） ~440Ma metamorphic event and coeval post-collisional granite magmatism have been firmly identified in the MCZ and its southern neighboring area. Together with previous data obtained by other researchers, our conclusions on the reconstruction of the tectonic architecture and evolution of the MCZ as follows: 1） The orogenic assemblages developed in the Mudanjiang collisional zone are featured by a sequence of ancient active continental margins and ensuing orogenic processing. The Mashan Group is the reworking basement of Jiamusi Massif, whereas the Heilongjiang Group represents arc and oceanic complexes characterized by imbricate deep-seated sliced and slivering sheets due to multi-phases of thrusting and nappe stacking. 2） The northern sub-belt of MCZ is probably the arc-continent collisional boundary related to the closure of main oceanic basin. The collisional age can be constrained by the events of syn-orogenic migmatization of migmatite, coeval S-type granite magmatism and the related granulite-facies metamorphism. Therefore, we suggested the collisional age of northern sub-belt is probably Cambrian-Early Ordovician. The extensive granulite-facies metamorphism of the Mashan Group in Jiamusi Massif, as affirmed by former works, was probably related with the collisional event. 3） The southern sub-belt of the MCZ was possibly related with the closure of back-arc basin. We presumed that the collisional age of southern sub-belt is about Ordovician-Early Silurian according to the ~440Ma extensive metamorphism and the occurrence of coeval post-collisional granite magmatism. 4） The extant structural architecture of the MCZ is related to the multi-phases of intra-continental superimposition, which is characterized by the Mesozoic nappe structure.

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

Orientation-based Continuum Damage Models for Rocks

A general formulation of the Helmholtz free energy used in thermodynamics of damage process of rocks is derived within a multi-scale framework. Such a physically-based thermodynamic state potential has a hybrid, discrete/continuum, nature in the sense tha