Retrograde melt–residue interaction and the formation of near‐anhydrous leucosomes in migmatites

Considering physical segregation of melt from its residue, the chemical potentials of the components (oxides) are the same in both when segregation occurs. Then, as P–T conditions change, gradients in chemical potential are established between the melt‐rich domains and residue permitting diffusional interaction to occur. In particular, on cooling, the chemical potential of H₂O becomes higher in the melt segregation than in the residue, particularly when biotite becomes stable in the residue assemblage. Diffusion of water from the melt to the residue promotes crystallization of anhydrous products from the melt and hydrous products in the residue. This diffusive process, when coupled with melt loss from the rocks subsequent to some degree of crystallization, can result in a significant degree of anhydrous leucosome being preserved in a migmatite with only minor retrogression of the residue. If H₂O can diffuse between the melt segregation and all of the residue, then no apparent selvedge between the two will be observed. Alternatively, if H₂O can diffuse between the melt segregation and only part of the residue, then a distinct selvedge may be produced. Diffusion of H₂O into the residue may be in part responsible for the commonly anhydrous nature of leucosomes, especially in granulite facies migmatites. Diffusion of other relatively mobile species such as Na₂O and K₂O has a lesser effect on overall melt crystallization but can change the proportion of quartz, plagioclase and K‐feldspar in the resultant leucosome. The diffusion of H₂O out of the melt results in the enhanced crystallization of the melt in the segregation and increases the amount of resulting anhydrous leucosome relative to the amount produced if melt crystallized in chemical isolation from the residue. For high residue:melt ratios, the proportion of resulting near‐anhydrous leucosome can approach that of the proportion of melt present at the onset of cooling with only minor loss of melt from a given segregation required. Crystallization of melt segregations via the diffusion of H₂O out of them into the host may also play a major role in driving melt‐rich segregations across key rheological transitions that would allow the expulsion of remaining melt from the system.     

Monazite as a monitor for melt‐rock interaction during cooling and exhumation

Granulite facies cordierite–garnet–biotite gneisses from the southeastern Reynolds Range, central Australia, contain both orthopyroxene‐bearing and orthopyroxene‐free quartzofeldspathic leucosomes. Mineral reaction microstructures at the interface of gneiss and leucosome observed in outcrop and petrographically, reflect melt‐rock interaction during crystallization. Accessory monazite, susceptible to fluid alteration, dissolution and recrystallization at high temperature, is tested for its applicability to constrain the chemical and P–T–time evolution of melt‐rock reactions during crystallization upon cooling. Bulk rock geochemistry and phase equilibria modelling constrain peak pressure and temperature conditions to 6.5–7.5 kbar and ~850°C, and U–Pb geochronology constrains the timing of monazite crystallization to 1.55 Ga, coeval with the Chewings Orogeny. Modelling predicts the presence of up to 15 vol.% melt at peak metamorphic conditions. Upon cooling below 800°C, melt extraction and in situ crystallization of melt decrease the melt volume to less than 7%, at which time it becomes entrapped and melt pockets induce replacement reactions in the adjacent host rock. Replacement reactions of garnet, orthopyroxene and K‐feldspar liberate Y, REE, Eu and U in addition to Mg, Fe, Al, Si and K. We demonstrate that distinguishing between monazite varieties solely on the basis of U–Pb ages cannot solve the chronological order of events in this study, nor does it tie monazite to the evolution of melt or stability of rock‐forming minerals. Rather, we argue that analyses of various internal monazite textures, their composition and overprinting relations allow us to identify the chronology of events following the metamorphic peak. We infer that retrograde reactions involving garnet, orthopyroxene and K‐feldspar can be attributed to melt‐rock interaction subsequent to partial melting, which is reflected in the development of compositionally distinct monazite textural domains. Internal monazite textures and their composition are consistent with dissolution and precipitation reactions induced by a high‐T melt. Monazite rims enriched in Y, HREE, Eu and U indicate an increased availability of these elements, consistent with the breakdown of orthopyroxene, garnet and K‐feldspar observed petrographically. Our study indicates that compositional and textural analysis of monazite in relation to major rock‐forming minerals can be used to infer the post‐peak chemical evolution of partial melts during high‐ to ultrahigh‐temperature metamorphism.     

Osumilite–melt interactions in ultrahigh temperature granulites: phase equilibria modelling and implications for the P–T–t evolution of the Eastern Ghats Province,India

The exposed residual crust in the Eastern Ghats Province records ultrahigh temperature (UHT) metamorphic conditions involving extensive crustal anatexis and melt loss. However, there is disagreement about the tectonic evolution of this late Mesoproterozoic–early Neoproterozoic orogen due to conflicting petrological, structural and geochronological interpretations. One of the petrological disputes in residual high Mg–Al granulites concerns the origin of fine‐grained mineral intergrowths comprising cordierite + K‐feldspar ± quartz ± biotite ± sillimanite ± plagioclase. These intergrowths wrap around porphyroblast phases and are interpreted to have formed by the breakdown of primary osumilite in the presence of melt trapped in the equilibration volume by the melt percolation threshold. The pressure (P)–temperature (T) evolution of four samples from three localities across the central Eastern Ghats Province is constrained using phase equilibria modelling in the chemical system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (NCKFMASHTO). Results of the modelling are integrated with published geochronological results for these samples to show that the central Eastern Ghats Province followed a common P–T–t history. This history is characterized by peak UHT metamorphic conditions of 945–955 °C and 7.8–8.2 kbar followed by a slight increase in pressure and close‐to‐isobaric cooling to the conditions of the elevated solidus at 940–900 °C and 8.5–8.3 kbar. In common with other localities from the Eastern Ghats Province, the early development of cordierite before osumilite and the peak to immediate post‐peak retrograde reaction between osumilite and melt to produce the intergrowth features requires that the prograde evolution was one of contemporaneous increasing pressure with increasing temperature. This counter‐clockwise (CCW) evolution is evaluated for one sample using inverse phase equilibria modelling along a schematic P–T path of 150 °C kbar^?1 starting from the low P–T end of the prograde P–T path as constrained by the phase equilibria modelling. The inverse modelling is executed by step‐wise down temperature reintegration of sufficient melt into the residual bulk chemical composition at the P–T point of the 1 mol.% melt isopleth at each step, representing the melt remaining on grain boundaries after each prograde drainage event, to reach the melt connectivity transition (MCT) of 7 mol.%. The procedure is repeated until a plausible protolith composition is recovered. The result demonstrates that clastic sedimentary rocks that followed a CCW P–T evolution could have produced the observed mineral assemblages and microstructures preserved in the central Eastern Ghats Province. This study also highlights the role of melt during UHT metamorphism, particularly its importance to both chemical and physical processes along the prograde and retrograde segments of the P–T path. These processes include: (i) an increase in diffusive length scales during the late prograde to peak evolution, creating equilibration volumes larger than a standard thin section; (ii) the development of retrograde mineral assemblages, which is facilitated if some melt is retained post‐peak; (iii) the presence of melt as a weakening mechanism and the advection of heat by melt, allowing the crust to thicken; and (iv) the effect of melt loss, which makes the deep crust both denser and stronger, and reduces heat production at depth, limiting crustal thickening and facilitating the transition to close‐to‐isobaric cooling.     

胶北地块高压与低压泥质麻粒岩的相平衡关系与p-T演化轨迹

利用最新的内洽性热力学数据库和THERMOCALC3.21程序对胶北地块高压与低压泥质麻粒岩的相平衡关系进行了定量分析。计算了胶北地块高压泥质麻粒岩、低压泥质麻粒岩和夕线石榴黑云片岩等代表性富铝岩石KFMASH(K2O-FeO-MgO-Al2O3-SiO2-H2O)体系的p-T视剖面图,再现了这些岩石随温压条件变化可能出现的各种矿物组合与矿物成分变化,发现原岩成分不同的变质岩石,尽管变质演化过程有所差异,但在麻粒岩相变质条件下所形成的矿物组合一致。通过计算泥质岩石在高压(p=1.0GPa)和低压(p=0.5GPa)条件下的T-X视剖面图,发现极度富铁、贫镁的岩石,在高压麻粒岩相条件下并不会生成含蓝晶石的特征矿物组合,在低压麻粒岩相条件下也不会生成含堇青石的特征矿物组合。将样品实际观测结果与p-T视剖面图的计算结果对比,确定胶北地块高压泥质麻粒岩变质峰期的温压条件为830~860℃,1.25~1.4GPa,峰期后呈现顺时针样式的p-T演化轨迹,反映陆壳先碰撞增厚、后又快速减薄的地质动力学过程;确定胶北地块低压泥质麻粒岩变质峰期的温压条件为790~820℃,0.62~0.68GPa,峰期后呈现近等压冷却的p-T演化轨迹。     

Mineral equilibria constraints on open‐system melting in metamafic compositions

The recent development of activity–composition relations for mineral and melt phases in high‐grade metamafic rocks allows mineral equilibria tools to be used to further aid our understanding of partial melting and the mineralogical consequences of melt segregation in these rocks. We show that bulk compositional data from natural amphibolites cover a wide compositional range, with particular variability in the content and ratios of Ca, Na and K indicating that low‐grade metasomatic alteration can substantially alter the igneous protolith chemistry and potentially affect the volume and composition of melt generated. Mineral equilibria calculations for five samples that span the compositional variability in our data set indicate that melting occurs primarily via the fluid‐absent breakdown of amphibole+quartz to produce a pressure‐sensitive peritectic assemblage of augite, orthopyroxene and/or garnet. The introduction of orthopyroxene at the onset of the amphibolite‐to‐granulite‐facies transition at lower pressure results in an increased rate of melt production until quartz is typically exhausted, and this is similarly seen for the introduction of garnet at higher pressure. Calculated melt compositions are dependent on the protolith composition, but initial solidus melting and biotite breakdown produce 1–3 mol.% of K‐rich granitic melts. As hornblende melting proceeds, 15–20 vol.% of either more granodioritic‐to‐tonalitic or granodioritic‐to‐trondhjemitic melt is produced. Once quartz is exhausted, intermediate to mafic melt compositions are produced at ultrahigh‐temperature conditions. Quartz‐rich lithologies with high Ca coupled to low Na and K are the most fertile under orogenic conditions, yielding up to 25 mol.% of sub‐alkalic granitic melt by 850°C. Such rocks did not experience significant subsolidus alteration. Altered compositions with low Ca and elevated Na and K are not as fertile, yielding less than 15 mol.% of alkalic granitic melt by 850°C. These melt volumes are enough to be segregated, and can make a contribution to granite magmatism and intracrustal differentiation that should not be overlooked.     

Large-scale melt-depletion in granulite terranes: an example from the Archean Ashuanipi Subprovince of Quebec

This study uses field, petrographic and geochemical methods to estimate how much granitic melt was formed and extracted from a granulite facies terrane, and to determine what the grain‐ and outcrop‐scale melt‐flow paths were during the melt segregation process. The Ashuanipi subprovince, located in the north‐eastern Superior Province of Quebec, is a large (90 000 km²) metasedimentary terrane, in which > 85% of the metasediments are of metagreywacke composition, that was metamorphosed at mid‐crustal conditions (820–900 °C and 6–7 kbar) in a late Archean dextral, transpressive orogen. Decrease in modal biotite and quartz as orthopyroxene and plagioclase contents increase, together with preserved former melt textures indicate that anatexis was by the biotite dehydration reaction: biotite + quartz + plagioclase = melt + orthopyroxene + oxides. Using melt/orthopyroxene ratios for this reaction derived from experimental studies, the modal orthopyroxene contents indicate that the metagreywacke rocks underwent an average of 31 vol% partial melting. The metagreywackes are enriched in MgO, CaO and FeO_t and depleted in SiO₂, K₂O, Rb, Cs, and U, have lower Rb/Sr, higher Rb/Cs and Th/U ratios and positive Eu anomalies compared to their likely protolith. These compositions are modelled by the extraction of between 20 and 40 wt %, granitic melt from typical Archean low‐grade metagreywackes. A simple mass balance indicates that about 640 000 km³ of granitic melt was extracted from the depleted granulites. The distribution of relict melt at thin section‐ and outcrop‐scales indicates that in layers without leucosomes melt extraction occurred by a pervasive grain boundary (porous) flow from the site of melting, across the layers and into bedding planes between adjacent layers. In other rocks pervasive grain boundary flow of melt occurred along the layers for a few, to tens of centimetres followed by channelled flow of melt in a network of short interconnected and structurally controlled conduits, visible as the net‐like array of leucosomes in some outcrops. The leucosomes contain very little residual material (< 5% biotite + orthopyroxene) indicating that the melt fraction was well separated from the residuum left in situ as melt‐depleted granulite. Only 1–3 vol percentage melt remained in the melt‐depleted granulites, hence, the extraction of melt generated by biotite dehydration melting in these granulites, was virtually complete under conditions of natural melting and strain rates in a contractional orogen.     

Phase equilibria modelling and zircon age dating of pelitic granulites in Zhaojiayao,from the Jining Group of the Khondalite Belt,North China Craton

The Jining Group occurs as the eastern segment of the Khondalite Belt, North China Craton and is dominated by a series of granulite facies rocks involving ‘normal’ pelitic granulites recording peak temperatures of ~850 °C and ultrahigh‐temperature (UHT) pelitic granulites recording peak temperatures of 950–1100 °C. The P–T paths and ages of these two types of granulites are controversial. Three pelitic granulite samples in the Jining Group comprising two sillimanite–garnet gneiss samples (J1208 and J1210) and one spinel–garnet gneiss sample (J1303) were collected from Zhaojiayao, where ‘normal’ pelitic granulites occur, for determination of their metamorphic evolution and ages. Samples J1208 and J1210 are interpreted to record cooling paths from the T_max stages with P–T conditions respectively of ~870–890 °C/7–8 kbar and >840 °C/>7.5 kbar constrained from the stability fields of the observed mineral assemblages and the isopleths of plagioclase, garnet and biotite compositions in pseudosections. Sample J1303 is interpreted to record pre‐T_max decompression from the kyanite‐stability fields to the T_max stage of 950–1020 °C/8–9 kbar and a post‐T_max cooling path revealed mainly from the stability field of the observed mineral assemblage, the plagioclase zoning and the biotite composition isopleth in pseudosections. The post‐T_max cooling stage can be divided into suprasolidus and subsolidus stages. The suprasolidus cooling may not result in an equilibrium state at the solidus in a rock. Therefore, different minerals may record different P–T conditions along the cooling path; the inferred maximum temperature is commonly higher than the solidus as well as different solidi being recorded for different samples from the same outcrop but experiencing different degrees of melt loss. Plagioclase compositions, especially its zoning in plagioclase‐rich granulites, are predicted to be useful for recording the higher temperature conditions of a granulite's thermal history. The three samples studied seem to record the temperature range covering those of the ‘normal’ and UHT pelitic granulites in the Jining Group, suggesting that UHT conditions may be reached in ‘normal’ granulites without diagnostic UHT indicators. LA‐ICP‐MS zircon U–Pb data provide a continuous trend of concordant ²⁰⁷Pb/²⁰⁶Pb ages from 1.89 to 1.79 Ga for sample J1210, and from 1.94 to 1.80 Ga for sample J1303. These continuous and long age spectrums are interpreted to represent a slow cooling and exhumation process corresponding to the post‐T_max cooling P–T paths recorded by the pelitic granulites, which may have followed the exhumation of deeply buried rocks in a thickened crust region resulted from a collision event at c. 1.95 Ga as suggested by the pre‐T_max decompression P–T path.     

Low‐P melting of metapelitic rocks and the role of H2O: Insights from phase equilibria modelling

Water‐fluxed melting has long been thought to have a minor influence on the thermal and chemical structure of the crust. We report here on amphibolite facies metasedimentary rocks from the 490–450 Ma Famatinian Orogen, in northwest Argentina, that have undergone water‐fluxed incongruent biotite melting at relatively low temperature, which have produced and lost a significant volume of melt. The protoliths consist of the turbiditic Puncoviscana Formation (Neoproterozoic to Early Cambrian). The field area exhibits a condensed metamorphic field gradient, from greenschist to amphibolite facies suprasolidus conditions, recording a low pressure almost isobaric path, reaching peak conditions estimated at 700°C at 4 kbar. Thermodynamic modelling in the MnNCKFMASHTO system is applied to investigate melting at such low pressure as a function of water content. Calculations using a typical turbidite composition show how small amounts of added free H₂O may increase significantly the melt fraction with little or no change in either the melt or residual phase compositions. They indicate negligible difference in normative An–Ab–Or proportions and ferromagnesian contents between melts derived by dehydration and water‐fluxed melts. The same is true for the content of H₂O dissolved in melts, which remains constant and the melt produced is granitic whether or not aqueous fluids are present. Thus, neither the residue nor the melt composition are indicators of the presence of aqueous fluids during anatexis. Recognizing the impact of small additions of H₂O to an anatectic terrane may therefore be difficult. The most significant change related to water‐fluxing is the relative proportions of minerals and melt fraction, rather than the actual mineral assemblage. The modal proportion of feldspar decreased while those of cordierite and biotite increased in the residual assemblages, as <5 mol.% of free H₂O was added. The impact of this addition is to more than double the proportion of water‐undersaturated melt to 25–30 mol.%. We have also developed a simple way to estimate how much melt a residual rock has lost, if the compositional trends of the protoliths are known. In summary, we find that even though the addition of small amounts of free H₂O impacts significantly on rock fertility, there is little obvious record in the field. The combined application of careful petrological investigation and thermodynamic modelling is the key to identify the influence of aqueous fluids, and exploit systems that became open not only to fluid influx but also to the extraction of melt.     

Rapid high‐T decompression recorded by Archean granulites in the northern Wyoming Province: Insights from petrological modelling

This study places new constraints on the pressure–temperature (P–T) path and duration of high‐temperature (HT) metamorphism recorded by Archean granulite facies metasedimentary rocks from the northern Wyoming Province in the eastern Beartooth Mountains, MT and WY, USA. These rocks exist as m‐ to km‐scale xenoliths within a c. 2.8 Ga calc‐alkaline granitoid batholith. Different interpretations of the timing of HT metamorphism relative to batholith intrusion in previous works have led to ambiguity over the mechanism by which these rocks were heated (i.e. batholith intrusion v. a later, cryptic event). The P–T path recorded by these rocks and the duration of this path may be indicative of the heating mechanism but are not currently well constrained. Here, we combine phase equilibria thermobarometry and diffusion modelling of major element zonation in garnet in order to constrain the P–T path of HT metamorphism and the durations of different parts of this path. It is shown that these rocks record a tight, clockwise P–T path characterized by near‐isobaric heating at ~6.5–7 kbar to ?770–800°C, HT decompression to ~6 kbar, 780–800°C, followed by limited decompression while cooling. Diffusion modelling of major element zonation in garnet suggests that HT decompression was brief (likely <1 Ma), and that cooling rates following this decompression were on the order of 10–100°C/Ma. Substantial changes in apparent thermal gradient along this P–T path indicate that the rocks record a significant but short‐lived thermal anomaly that occurred in the Wyoming mid‐crust in the Late Archean.     

Phase equilibria modelling of kyanite‐bearing anatectic paragneisses from the central Grenville Province

Kyanite‐bearing paragneisses from the Manicouagan Imbricate Zone and its footwall (high‐P belt of the central Grenville Province) preserve evidence of partial melting with development of metamorphic textures involving biotite–garnet ± kyanite ± plagioclase ± K‐feldspar–quartz. Garnet in these rocks displays a variety of zoning patterns with respect to Ca. Pseudosection modelling in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system using measured bulk rock compositions accounts for the textural evolution of two aluminous and two sub‐aluminous samples from the presumed thermal peak to conditions at which retained melt solidified. The prograde features are best explained by pseudosections calculated with compositions to account for melt loss. The intersection of isopleths of grossular content and Fe/(Fe + Mg) relating to large porphyroblasts of garnet provide constraints on the P–T conditions of the metamorphic peak. These P–T estimates are considered to be minima because of the potential for diffusional modification of the composition of garnet at high‐T and during the early stages of cooling. However, they are consistent with textural observations and pseudosection topology, with peak assemblages best preserved in rocks for which the calculated pseudosections predict only small changes in mineral proportions in the P–T interval, in which retrograde reactions are inferred to have occurred between the thermal peak and the solidus. Maximum P–T conditions (14.5–15.5 kbar and 840–890 °C) and steep retrograde P–T paths inferred for rocks from the Manicouagan Imbricate Zone are comparable with those determined for mafic rocks from the same area. In contrast, maximum P–T conditions of 12.5–13 kbar and 815–830 °C and flatter P–T paths are inferred for the rocks of the footwall to the Manicouagan Imbricate Zone. The general consistency between textures, mineral compositions and the topologies of the calculated pseudosections suggests that the pseudosection approach is an appropriate tool for inferring the P–T evolution of high‐P anatectic quartzo‐feldspathic rocks.     

Challenges in constraining the P–T conditions of mafic granulites: An example from the northern Trans‐North China Orogen

Some mafic granulites in the Sanggan area of the northern Trans‐North China Orogen (TNCO) have a relatively simple mineralogy with low energy grain shapes that are compatible with an assumption of equilibrium, but the rock‐forming minerals show variations in composition that create challenges for thermobarometry. The mafic granulites, which occur as apparently disrupted dyke‐like bodies in tonalite–trondhjemite–granodiorite gneisses, are divided into two types based on petrography and chemical composition. Type 1 mafic granulites are fine‐ to medium‐grained with an equilibrated texture and an assemblage of plagioclase+clinopyroxene+garnet+magnetite+ilmenite and sometimes minor hornblende±orthopyroxene. Type 2 mafic granulites are coarse‐grained and hornblende bearing with a peak assemblage of garnet+clinopyroxene+plagioclase+hornblende and variably developed coronae and symplectites of plagioclase+hornblende+orthopyroxene partially replacing porphyroblastic garnet±clinopyroxene. SIMS U–Pb dating of metamorphic zircon from two type 1 mafic granulites yields metamorphic ages of c. 1.84 and 1.83 Ga, consistent with published ages of the type 2 mafic granulites. Based on phase equilibrium modelling, we use the common overlap of P–T fields defined by the mineral assemblage limits, and the mole proportion and composition isopleths of different minerals in each sample to quantify the metamorphic conditions. For type 1 granulites, overlap of the mineral proportion and composition fields for each of three samples yields similar P–T conditions of 710–880°C at 0.57–0.79 GPa, 820–850°C at 0.59–0.63 GPa and 800–860°C at 0.59–0.68 GPa. For the type 2 granulites, overlaying the peak assemblage fields for three samples yields common P–T conditions of 870–890°C at 1.1–1.2 GPa. For the retrograde assemblage, overlap of the mineral proportion and composition fields for each sample yields similar P–T conditions of 820–840°C at 0.85–0.88 GPa, 860–880°C at 0.83–0.86 GPa and 880–930°C at 0.89–0.95 GPa. The P–T conditions appear distinct between the two types of mafic granulite, with the mineralogically simple type 1 mafic granulites recording the lowest pressures. However, there are significant uncertainties associated with these results. For the granulites, there are uncertainties related to the determination of modes and composition of the equilibration volume, particularly estimation of O and H₂O contents, and in the phase equilibrium modelling there are uncertainties that propagate through the calculation of mole proportions and mineral compositions. The compound uncertainties on pressure and temperature for high‐T granulites are large and the results of our study show that it may be unwise to rely on P–T conditions determined from the simple intersection of calculated mineral composition isopleths alone. Since the samples in this study are from a limited area—a few hundred square metres—we infer that they record a single P–T path involving both decompression and cooling. However, there is no evidence of the high‐P granulite facies event at 1.93–1.90 Ga that is recorded elsewhere in the TNCO, which suggests that the precursor basic dykes were emplaced late during the assembly of the North China Craton.     

Contrasting metamorphic histories of lenses of high‐pressure rocks and host migmatites with a flat orogenic fabric (Bohemian Massif,Czech Republic): a result of tectonic mixing within horizontal crustal flow?

Migmatites with sub‐horizontal fabrics at the eastern margin of the Variscan orogenic root in the Bohemian Massif host lenses of eclogite, kyanite‐K‐feldspar granulite and marble within a matrix of migmatitic paragneiss and amphibolite. Petrological study and pseudosection modelling have been used to establish whether the whole area experienced terrane‐wide exhumation of lower orogenic crust, or whether smaller portions of higher‐pressure lower crust were combined with a lower‐pressure matrix. Kyanite‐K‐feldspar granulite shows peak conditions of 16.5 kbar and 850 °C with no clear indications of prograde path, whereas in the eclogite the prograde path indicates burial from 10 kbar and 700 °C to a peak of 18 kbar and 800 °C. Two contrasting prograde paths are identified within the host migmatitic paragneiss. The first path is inferred from the presence of staurolite and kyanite inclusions in garnet that contains preserved prograde zoning that indicates burial with simultaneous heating to 11 kbar and 800 °C. The second path is inferred from garnet overgrowths of a flat foliation defined by sillimanite and biotite. Garnet growth in such an assemblage is possible only if the sample is heated at 7–8 kbar to around 700–840 °C. Decompression is associated with strong structural reworking in the flat fabric that involves growth of sillimanite in paragneiss and kyanite‐K‐feldspar granulite at 7–10 kbar and 750–850 °C. The contrasting prograde metamorphic histories indicate that kilometre‐scale portions of high‐pressure lower orogenic crust were exhumed to middle crustal levels, dismembered and mixed with a middle crustal migmatite matrix, with the simultaneous development of a flat foliation. The contrasting P–T paths with different pressure peaks show that tectonic models explaining high‐pressure boudins in such a fabric cannot be the result of heterogeneous retrogression during ductile rebound of the whole orogenic root. The P–T paths are compatible with a model of heterogeneous vertical extrusion of lower crust into middle crust, followed by sub‐horizontal flow.     

Conductively driven,high‐thermal gradient metamorphism in the Anmatjira Range,Arunta region,central Australia

LA–ICP–MS in situ U–Pb monazite geochronology and P–T pseudosections are combined to evaluate the timing and physical conditions of metamorphism in the SE Anmatjira Range in the Aileron Province, central Australia. All samples show age peaks at c. 1580–1555 Ma, with three of five samples showing additional discrete age peaks between c. 1700 and 1630 Ma. P–T phase diagrams calculated for garnet–sillimanite–cordierite–K‐feldspar–ilmenite–melt bearing metapelitic rocks have overlapping peak mineral assemblage stability fields at ~870–920 °C and ~6.5–7.2 kbar. P–T modelling of a fine‐grained spinel–cordierite–garnet–biotite reaction microstructure suggests retrograde P–T conditions evolved down pressure and temperature to ~3–5.5 kbar and ~610–850 °C. The combined geochronological and P–T results indicate the SE Anmatjira Range underwent high‐temperature, low‐pressure metamorphism at c. 1580–1555 Ma, and followed an apparently clockwise retrograde path. The high apparent thermal gradient necessary to produce the estimated P–T conditions does not appear to reflect decompression of high‐P assemblages, nor is there syn‐metamorphic magmatism or structural evidence for extension. Similar to previous workers, we suggest the high‐thermal gradient P–T conditions could have been achieved by heating, largely driven by high heat production from older granites in the region.     

Petrology,phase equilibria modelling and zircon U–Pb geochronology of Paleoproterozoic mafic granulites from the Fuping Complex,North China Craton

The Fuping Complex is one of the important basement terranes within the central segment of the Trans‐North China Orogen (TNCO) where mafic granulites are exposed as boudins within tonalite–trondhjemite–granodiorite (TTG) gneisses. Garnet in these granulites shows compositional zoning with homogeneous cores formed in the peak metamorphic stage, surrounded by thin rims with an increase in almandine and decrease in grossular contents suggesting retrograde decompression and cooling. Petrological and phase equilibria studies including pseudosection calculation using thermocalc define a clockwise P–T path. The peak mineral assemblages comprise garnet+clinopyroxene+amphibole+quartz+plagioclase+K‐feldspar+ilmenite±orthopyroxene±magnetite, with metamorphic P–T conditions estimated at 8.2–9.2 kbar, 870–882 °C (15FP‐02), 9.6–11.3 kbar, 855–870 °C (15FP‐03) and 9.7–10.5 kbar, 880–900 °C (15FP‐06) respectively. The pseudosections for the subsequent retrograde stages based on relatively higher H₂O contents from P/T–M(H₂O) diagrams define the retrograde P–T conditions of <6.1 kbar, <795 °C (15FP‐02), 5.6–5.8 kbar, <795 °C (15FP‐03), and <9 kbar, <865 °C (15FP‐06) respectively. Data from LA‐ICP‐MS zircon U–Pb dating show that the mafic dyke protoliths of the granulite were emplaced at c. 2327 Ma. The metamorphic zircon shows two groups of ages at 1.96–1.90 Ga (peak at 1.93–1.92 Ga) and 1.89–1.80 Ga (peak at 1.86–1.83 Ga), consistent with the two metamorphic events widely reported from different segments of the TNCO. The 1.93–1.92 Ga ages are considered to date the peak granulite facies metamorphism, whereas the 1.86–1.83 Ga ages are correlated with the retrograde event. Thus, the collisional assembly of the major crustal blocks in the North China Craton (NCC) might have occurred during 1.93–1.90 Ga, marking the final cratonization of the NCC.     

An electron‐optical study of melt‐related microstructures in granulite facies rocks from the Torngat Orogen

Cordierite–quartz and plagioclase–quartz intergrowths in a paragneiss from northern Labrador (the Tasiuyak Gneiss) were studied using SEM, STEM and TEM. The gneiss experienced granulite facies conditions and partial melting during both regional and, subsequently, during contact metamorphism. The microstructures examined all results from the contact metamorphism. Cordierite–quartz intergrowths occur on coarse and fine scales. The former sometimes exist as a ‘geometric’ intergrowth in which the interface between cordierite and quartz appears planar at the resolution of the optical microscope and SEM. The latter exists in several microstructural variants. Plagioclase is present as a minor component of the intergrowth in some examples of both the coarse and fine intergrowth. Grain boundaries in cordierite–quartz intergrowths are occupied by amorphous material or a mixture of amorphous material and chlorite. Cordierite and quartz are terminated by crystal faces in contact with amorphous material. Chlorite is sometimes found on cordierite surfaces and penetrating into cordierite grains along defects. Quartz contains (former) fluid inclusions 10–20 nm in maximum dimension. The presence of planar interfaces between cordierite and the amorphous phase is reminiscent of those between crystals and glass in volcanic rocks, but in the absence of compelling evidence that the amorphous material represents former melt, it is interpreted as a reaction product of cordierite. Plagioclase–quartz intergrowths occur in a number of microstructural variants and are commonly associated with cordierite–quartz intergrowths. The plagioclase–quartz intergrowths display simple, non‐planar interfaces between plagioclase and quartz. Quartz contains (former) fluid inclusions of dimensions similar to those observed in cordierite–quartz intergrowths. The boundary between quartz and enclosing K‐feldspar is cuspate, with quartz cusps penetrating a few tens of nanometres into K‐feldspar, commonly along defects in K‐feldspar and sometimes with very low dihedral angles at their tips. This cuspate microstructure is interpreted as melt pseudomorphs. The plagioclase–quartz intergrowths share some features with myrmekite, but differ in some respects: the composition of the plagioclase (An₃₇Ab₆₂Or₁–An₃₈Ab₆₁Or₁); the association with cordierite–quartz intergrowths; and microstructures that are atypical of myrmekite (e.g. quartz vermicules shared with cordierite–quartz intergrowths). It is inferred that the plagioclase–quartz intergrowths may have formed from, or in the presence of, melt. Inferred melt‐related microstructures preserved on the nanometre scale suggest that melt on grain boundaries was more pervasive than is evident from light optical and SEM observations.     

The reaction history of kyanite in high‐P aluminous granulites

Cathodoluminescence (CL) mapping of kyanite in high pressure, aluminous granulites from the central Grenville Province reveals internal structures that are linked to their metamorphic reaction history. In two samples, individual kyanite crystals are shown to be composite porphyroblasts comprising three distinct generations, defined by their CL intensity and Cr (±V, Ti, Fe and Ga) content, and each separated by resorbed interfaces. In contrast, a sub‐aluminous sample contains two types of kyanite, one as resorbed inclusions in garnet and another in the groundmass or replacing garnet. These textural variants of kyanite are interpreted within the framework of phase equilibria modelling. In P–T pseudosections, a first generation of kyanite, which is only present in the most aluminous samples, is potentially linked to staurolite breakdown, and its resorption is consistent with a subsequent increase in pressure. This kyanite represents the earliest remnant of prograde metamorphism identifiable in these rocks. The second generation, present in the porphyroblasts in the same samples and as inclusions in garnet in the sub‐aluminous sample, is interpreted to be the peritectic product of muscovite dehydration melting. Resorption of this kyanite is consistent with subsequent continuous dehydration melting of biotite, which is also inferred based on microstructural considerations. The final generation of kyanite, present as rims on the prograde kyanite porphyroblasts in aluminous samples and as part of the groundmass or replacing garnet in the sub‐aluminous rock, is interpreted to have grown during melt crystallization upon retrogression. The presence of retrograde kyanite implies that the melt crystallized over a wide range of temperatures, and provides an important constraint on the P–T conditions of the metamorphic peak and on the retrograde P–T path. CL mapping is crucial for identifying retrograde kyanite in aluminous samples, as it preferentially overgrows existing kyanite rather than replacing other prograde phases. The scarcity of kyanite in sub‐aluminous rocks allows retrograde kyanite to grow as discrete crystals that can be identified by optical microscopy. This work attests to the potential of unconventional tools such as CL imaging for deciphering the metamorphic history of rocks.     

Genesis of the Kapuskasing (Ontario) migmatitic mafic granulites by dehydration melting of amphibolite: the importance of quartz to reaction progress

Migmatitic, granulite-grade mafic gneisses make up a significant part of the Kapuskasing Structural Zone (KSZ), Ontario. Although they contain a common mineral assemblage [hornblende (Hbl)+plagioclase (Pl)+diopside (Di)±garnet (Grt)+quartz (Qtz)±titanite (Ttn)], the mafic gneisses show wide variations in modal mineralogy from hornblende-rich to diopside+garnet-rich varieties and all gradations between. Up to 25 vol.% segregated plagioclase+quartz-rich (trondhjemitic) leucosome (Tdh) is intimately associated with the mafic gneiss, occurring in a continuum of patches, veins and transecting dykes at scales ranging from decimetres to micrometres. The texture and composition of the leucosome, combined with P-T estimates for the host rocks above the solidus, suggest it represents crystallized trondhjemitic melt. Quartz is mainly restricted to the segregated leucosomes but more rarely occurs in a variety of interstitial textures in the mafic gneiss, suggesting that it crystallized from a melt phase rather than having been present as a solid phase at peak metamorphic conditions. Modal and textural data indicate a reaction relationship of the form: Hbl+Pl(+Qtz?)=Grt+Di+Ttn+leucosome (Tdh), implying that the granulite-forming process involved dehydration melting of an amphibolite protolith. Pressure-temperature estimates from Grt+Di+Pl+Qtz geothermobarometry are 9 kbar and 685-735 °C; however, based on experimental studies of dehydration melting of amphibolite, we estimate that peak conditions were closer to 11 kbar, 850 °C. Mass balance analysis, using the technique of singular value decomposition, and reaction space analysis were used to quantify the reaction and to determine the controls on reaction progress. The following mass balance provides a model for the natural reaction:1.00 Hbl+0.92 Pl+3.76 Qtz=1.14 Grt+1.54 Di+0.21 Ttn+1.49 Tdh+0.14 ‘pg’+0.39 Fe_?1Mg+0.33 NaSiCa_?1Al_?1where ‘pg’ is a pargasite-like exchange. In all model mass balances tested, quartz is a reactant with a large coefficient. We argue that the abundance of quartz in the amphibolite protolith was the primary control on the differing extents of reaction observed. Mineral compositional variation exerted a secondary control on reaction progress, with Fe-richer layers containing An-richer plagioclase and more actinolitic amphibole reacting earliest (i.e. at lowest temperatures). Comparison of the calculated amount of melt produced in the gneisses with that now observed implies expulsion of 5–30% of the melt. These volumes are similar to those predicted from REE modelling of Archaean tonalities and trondhjemites from a garnet amphibolite source, suggesting that the KSZ mafic gneisses may be representative of partially depleted source rocks for trondhjemite-tonalite generation.     

Phase equilibria constraints on melting of stromatic migmatites from Ronda (S. Spain): insights on the formation of peritectic garnet

Stromatic metatexites occurring structurally below the contact with the Ronda peridotite (Ojén nappe, Betic Cordillera, S Spain) are characterized by the mineral assemblage Qtz+Pl+Kfs+Bt+Sil+Grt+Ap+Gr+Ilm. Garnet occurs in low modal amount (2–5 vol.%). Very rare muscovite is present as armoured inclusions, indicating prograde exhaustion. Microstructural evidence of melting in the migmatites includes pseudomorphs after melt films and nanogranite and glassy inclusions hosted in garnet cores. The latter microstructure demonstrates that garnet crystallized in the presence of melt. Re‐melted nanogranites and preserved glassy inclusions show leucogranitic compositions. Phase equilibria modelling of the stromatic migmatite in the MnO–Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂–O₂–C (MnNCaKFMASHOC) system with graphite‐saturated fluid shows P–T conditions of equilibration of 4.5–5 kbar, 660–700 °C. These results are consistent with the complete experimental re‐melting of nanogranites at 700 °C and indicate that nanogranites represent the anatectic melt generated immediately after entering supersolidus conditions. The P–T estimate for garnet and melt development does not, however, overlap with the low‐temperature tip of the pure melt field in the phase diagram calculated for the composition of preserved glassy inclusions in garnet in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (NCKFMASH) system. A comparison of measured melt compositions formed immediately beyond the solidus with results of phase equilibria modelling points to the systematic underestimation of FeO, MgO and CaO in the calculated melt. These discrepancies are present also when calculated melts are compared with low‐T natural and experimental melts from the literature. Under such conditions, the available melt model does not perform well. Given the presence of melt inclusions in garnet cores and the P–T estimates for their formation, we argue that small amounts (<5 vol.%) of peritectic garnet may grow at low temperatures (≤700 °C), as a result of continuous melting reactions consuming biotite.     

Structural,petrological and chemical analysis of syn‐kinematic migmatites: insights from the Western Gneiss Region,Norway

Evidence of melting is presented from the Western Gneiss Region (WGR) in the core of the Caledonian orogen, Western Norway and the dynamic significance of melting for the evolution of orogens is evaluated. Multiphase inclusions in garnet that comprise plagioclase, potassic feldspar and biotite are interpreted to be formed from melt trapped during garnet growth in the eclogite facies. The multiphase inclusions are associated with rocks that preserve macroscopic evidence of melting, such as segregations in mafic rocks, leucosomes and pegmatites hosted in mafic rocks and in gneisses. Based on field studies, these lithologies are found in three structural positions: (i) as zoned segregations found in high‐P (ultra)mafic bodies; (ii) as leucosomes along amphibolite facies foliation and in a variety of discordant structures in gneiss; and (iii) as undeformed pegmatites cutting the main Caledonian structures. Segregations post‐date the eclogite facies foliation and pre‐date the amphibolite facies deformation, whereas leucosomes are contemporaneous with the amphibolite facies deformation, and undeformed pegmatites are post‐kinematic and were formed at the end of the deformation history. The geochemistry of the segregations, leucosomes and pegmatites in the WGR defines two trends, which correlate with the mafic or felsic nature of the host rocks. The first trend with Ca‐poor compositions represents leucosome and pegmatite hosted in felsic gneiss, whereas the second group with K‐poor compositions corresponds to segregation hosted in (ultra)mafic rocks. These trends suggest partial melting of two separate sources: the felsic gneisses and also the included mafic eclogites. The REE patterns of the samples allow distinction between melt compositions, fractionated liquids and cumulates. Melting began at high pressure and affected most lithologies in the WGR before or during their retrogression in the amphibolite facies. During this stage, the presence of melt may have acted as a weakening mechanism that enabled decoupling of the exhuming crust around the peak pressure conditions triggering exhumation of the upward‐buoyant crust. Partial melting of both felsic and mafic sources at temperatures below 800 °C implies the presence of an H₂O‐rich fluid phase at great depth to facilitate H₂O‐present partial melting.     

井间剩余油富集区确定方法研究

剩余油富集部位研究是油田开发后期的重要工作之一,目前有很多成熟的方法和经验。本文提供的方法是利用常规动静态资料快速、准确的进行判断,不需要大型的数值模拟等工作,使油藏管理人员随时监控和调整生产井的工作制度,充分动用地质储量,提高采收率。