Transcrystalline melt migration in clinopyroxene

Glass inclusions in clinopyroxene phenocrysts from La Sommata (Vulcano Island, Aeolian Arc) were reheated and submitted to a sustained thermal gradient. Each remelted inclusion undergoes a transient textural and chemical reequilibration and concomitantly begins to migrate along a crystallographic direction, at a small angle with the thermal gradient. The completion of morphological evolution requires a characteristic time that is governed by chemical diffusion. Chemical reequilibration results in the formation of a colored halo that delineates the former location and shape of the inclusion after it has migrated away. Transcrystalline migration proceeds by dissolution of the host clinopyroxene ahead and precipitation astern. Its rate is not limited by Fick’s law, but by the crystal-melt interface kinetics. Clinopyroxene dissolution and growth are slower than for olivine in similar conditions but obey the same analytical law, which can be transposed to equally or more sluggish melting or crystallization events in nature. When a gas bubble is initially present, it responds to elastic forces by quickly shifting toward the cold end of the inclusion, where it soon becomes engulfed as an isolated fluid inclusion in the reprecipitated crystal. This study confirms that transcrystalline melt migration, beside its possible implications for small-scale melt segregation and fluid-inclusion generation in the Earth’s mantle, provides an experimental access to interfacial kinetic laws in near-equilibrium conditions.     

The inception and growth of leucosomes: microstructure at the start of melt segregation in migmatites

The beginning stages of melt segregation and the formation of leucosomes are rarely preserved in migmatites. Most arrays of leucosomes record a more advanced stage where flow dominates over segregation. However, the early stages in the formation of leucosomes and the segregation of melt are preserved in a partially melted meta‐argillite from the metatexite zone (>800 °C) of the contact aureole around the Duluth Complex, Minnesota. The rock contains 2.4 modal% leucosome in a matrix consisting of 40.5% in situ neosome and 57.1% cordierite + plagioclase framework. The domainal microstructure in the matrix is a pre‐anatectic feature resulting from the bulk composition. Terminal chlorite reactions produced a large volume of cordierite which, with plagioclase, formed a framework that enclosed patches of biotite + quartz + plagioclase ± K‐feldspar. Upon melting, these fertile domains became patches of in situ neosome. Plagioclase in the neosome is less sodic than in the leucosome, hence segregation of melt occurred during crystallization, not melting. Segregation was delayed because the cordierite + plagioclase framework was strong enough to resist dilatation and compaction until after crystallization started. The leucosomes are small (i.e. they are microleucosomes) and display a systematic progression in morphology as length and aspect ratio increase from ~1 to 19 mm and from ~2.5 to >30 respectively. Small equant micropores form first, and in places these coalesce into small (~1 mm, aspect ratio ~2.5), isolated, blunt‐ended, elliptical microleucosomes. In the next stage, micropores develop ahead of, and at ~45° to the left and right of the blunt tip of a microleucosome; one of these develops into an elliptical leucosome and an en echelon array of either a left‐ or right‐stepping elliptical microleucosome forms. Each elliptical microleucosome in the en echelon arrays is separated by a bridge of matrix. Next, microleucosomes of greater length (>4 mm) and aspect ratio (>5) form when the bridges of cordierite + plagioclase matrix rupture and the elliptical microleucosomes link together to form a zigzag‐shaped microleucosome. Finally, still longer microleucosomes with greater aspect ratios (~30) are formed by the joining of zigzag arrays. Such a progression is characteristic of the way ductile fractures grow. The segregation of melt was driven by the pressure gradient between the dilatant fracture and an adjacent in situ neosome, which drew melt to the growing fracture, thereby creating a microleucosome. The microleucosomes are filled arrays of ductile fractures. Melt was contiguous only between microleucosomes and adjacent patches of in situ neosome. The length‐scale of segregation was ~5 mm, the size of a typical patch of in situ neosome, and restricted by the surrounding impermeable cordierite + plagioclase framework. The melt in the microleucosome was the most fractionated and the last to crystallize. All microleucosomes contain entrained minerals as a consequence of their mechanism of growth. Rupture of the bridges resulted in the entrainment of pre‐anatectic phases. However, microleucosomes that cross patches of in situ neosome are also contaminated with peritectic phases that were transported with the melt.     

Experimental anatexis and genesis of migmatites: A reply

Recentlyvon Platen applied the results of his laboratory experiments to a migmatite studied byHärme. Some comments are presented upon that re-interpretation.     

Continuous vs. discontinuous melt segregation in migmatites: insights from a cellular automaton model

A cellular automaton model is presented that simulates melt extraction from migmatites. The varying parameters are the relative amount of melt, thresholds for melt connection leading to movement and melt escape, and type of deformation (pure and/or simple shear). It is shown that when the melt percentage is about half that of the melt escape threshold, the competition between production and deformation-assisted extraction results in a discontinuous rate of melt extraction. The process can be compared to stick-slip motion during frictional experiments. At higher melt fraction, it is continuously extracted. Deformation assists significantly in segregating melt, particularly when a pure shear (compaction) component is added to simple shear.     

Chemical characteristics of migmatites: accessory phase distribution and evidence for fast melt segregation rates

 Equilibration between melt and solid is inhibited by rapid melt extraction and by restricted equilibration (armouring, slow dissolution). When segregation occurs by channelised migration along high-porosity pathways, melt migration is more rapid than trace element diffusion rates in silicates and faster than accessory phase dissolution rates. Evidence for channelised flow and deformation-enhanced melt segregation into boudin necks, fractures and micro-shears at low melt fractions is present in the Moine Kirtomy Migmatitie Suite (KMS) in Sutherland, Scotland. Melt migration distances are on a metre to tens of metres scale. Concordant leucosomes in stromatic migmatities in the KMS have low Zr contents, low LREE (light rare-earth element) and H (heavy) REE contents and positive Eu anomalies. REE patterns of this type can be produced by removal of leucosome before complete equilibration with source due to the inhibited dissolution of LREE- and HREE-bearing accessory phases in water-undersaturated melts. Melting in the KMS, however, occurred at or near the wet granite solidus, leaving biotite as a residual phase. Detailed back-scattered electron imaging shows that REE-bearing accessory phases remained as residual phases, and were concentrated in the melanosome and at the melanosome-leucosome boundary. Irregularly shaped patches of diatexite contain a small proportion of excess Zr, consistent with entrainment of melanosome-schlieren enriched in zircon. These data indicate that deformation-enhanced melt extraction led to the rapid migration of small melt fractions from the melting site on a time-scale less than that required to saturate the melt in Zr. Leucosomes were thus prevented from equilibrating with accessory phases before extraction. Received: 12 July 1995 / Accepted: 4 March 1996     

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.     

Impact of gneissic layering and localized incipient melting upon melt flow during experimental deformation of migmatites

In this study, we test experimentally the role of compositional layering as a key parameter for controlling melt flow in a natural migmatite during coaxial deformation. We performed in – situ pure-shear experiments on two natural gneisses. The first gneiss is weakly foliated with minerals homogenously distributed. The second gneiss shows a pronounced compositional layering of alternating quartz – feldspar – rich and biotite – muscovite – rich layers. Experimental conditions were selected to obtain homogeneous melt distribution in the homogeneous gneiss and heterogeneous melt distribution in the layered gneiss. Initial melt distribution is not modified by deformation in experiments on the homogeneous gneiss, implying that melting products did not migrate from their initiation sites. In contrast, melt flowed in shear zones or in inter-boudin positions during experimental deformation of the heterogeneous gneiss. These experiments attest to the strong influence of initial gneissic layering on melting pattern, melt segregation and flow during deformation of partially molten rocks.     

pH 值对有机磷在弱透水层中迁移转化的试验模拟

为研究渗透水pH 值对有机磷通过弱透水黏土层迁移转化的影响规律,采用黏性土层充压渗透装置模拟污染潜水中有机磷渗透通过弱透水层的迁移转化过程。研究发现,有机磷组分通过弱透水层的迁移能力总体较弱,不同pH 值条件下有机磷的矿化率为93. 5% ～ 99. 95%; 在进水pH 值为5. 5 ～ 8. 5 时,进水pH 值越高,土层中磷酸酶的活性也越高,降解转化有机磷的能力也越强; 黏土对有机磷的阻滞能力也增强。进水pH 值为5. 5 时,迁移率为6. 5%; 当pH 达到8. 5 时,矿化率明显增加,迁移率＜ 0. 05%,使其不易渗入深层地下水引起污染。     

Roles for fluid and/or melt advection in forming high-P mafic migmatites, Fiordland, New Zealand

A series of striking migmatitic structures occur in rectilinear networks through western Fiordland, New Zealand, involving, for the most part, narrow anorthositic dykes that cut hornblende‐bearing orthogneiss. Adjacent to the dykes, host rocks show patchy, spatially restricted recrystallization and dehydration on a decimetre‐scale to garnet granulite. Although there is general agreement that the migration of silicate melt has formed at least parts of the structures, there is disagreement on the role of silicate melt in dehydrating the host rock. A variety of causal processes have been inferred, including metasomatism due to the ingress of a carbonic, mantle‐derived fluid; hornblende‐breakdown leading to water release and limited partial melting of host rocks; and dehydration induced by volatile scavenging by a migrating silicate melt. Variability in dyke assemblage, together with the correlation between dehydration structures and host rock silica content, are inconsistent with macroscopic metasomatism, and best match open system behaviour involving volatile scavenging by a migrating trondhjemitic liquid.     

Origin of migmatites by deformation-enhanced melt infiltration of orthogneiss: a new model based on quantitative microstructural analysis

A detailed field study reveals a gradual transition from high‐grade solid‐state banded orthogneiss via stromatic migmatite and schlieren migmatite to irregular, foliation‐parallel bodies of nebulitic migmatite within the eastern part of the Gföhl Unit (Moldanubian domain, Bohemian Massif). The orthogneiss to nebulitic migmatite sequence is characterized by progressive destruction of well‐equilibrated banded microstructure by crystallization of new interstitial phases (Kfs, Pl and Qtz) along feldspar boundaries and by resorption of relict feldspar and biotite. The grain size of all felsic phases decreases continuously, whereas the population density of new phases increases. The new phases preferentially nucleate along high‐energy like–like boundaries causing the development of a regular distribution of individual phases. This evolutionary trend is accompanied by a decrease in grain shape preferred orientation of all felsic phases. To explain these data, a new petrogenetic model is proposed for the origin of felsic migmatites by melt infiltration from an external source into banded orthogneiss during deformation. In this model, infiltrating melt passes pervasively along grain boundaries through the whole‐rock volume and changes completely its macro‐ and microscopic appearance. It is suggested that the individual migmatite types represent different degrees of equilibration between the host rock and migrating melt during exhumation. The melt topology mimicked by feldspar in banded orthogneiss forms elongate pockets oriented at a high angle to the compositional banding, indicating that the melt distribution was controlled by the deformation of the solid framework. The microstructure exhibits features compatible with a combination of dislocation creep and grain boundary sliding deformation mechanisms. The migmatite microstructures developed by granular flow accompanied by melt‐enhanced diffusion and/or melt flow. However, an AMS study and quartz microfabrics suggest that the amount of melt present did not exceed a critical threshold during the deformation to allow free movements of grains.     

Petrogenesis of metatexite and diatexite migmatites determined using zircon U–Pb age,trace element and Hf isotope data,Higo metamorphic terrane,central Kyushu,Japan

Metatexite and diatexite migmatites are widely distributed within the upper amphibolite and granulite facies zones of the Higo low‐P/high‐T metamorphic terrane. Here, we report data from an outcrop in the highest grade part of the granulite facies zone, in which diatexite occurs as a 3 m thick layer between 2 m thick layers of stromatic‐structured metatexite within pelitic gneiss. The migmatites and gneiss contain the same peak mineral assemblage of biotite + plagioclase + quartz + garnet + K‐feldspar with retrograde chlorite ± muscovite and some accessory minerals of ilmenite ± rutile ± titanite + apatite + zircon + monazite ± pyrite ± zinc sulphide ± calcite. Calculated metamorphic P–T conditions are 800–900 °C and 9–12 kbar. Zircon in the diatexite forms elongate euhedral crystals with oscillatory zoning, but no core–rim structure. Zircon from the gneiss and metatexite forms euhedral–subhedral grains comprising inherited cores overgrown by thin rims. The overgrowth rims in the metatexite have lower Th/U ratios than zircon in the diatexite and yield a ²⁰⁶Pb/²³⁸U age of 116.0 ± 1.6 Ma, which is older than the 110.1 ± 0.6 Ma ²⁰⁶Pb/²³⁸U age derived from zircon in the diatexite. Zircon from the diatexite has variable REE contents with convex upward patterns and flat normalized HREE, whereas the overgrowth rims in the metatexite and gneiss have steep HREE‐enriched patterns; however, both types have similar positive Ce and negative Eu anomalies. ¹⁷⁶Hf/¹⁷⁷Hf ratios in the overgrowth rims from the metatexite are more variable and generally lower than values from zircon in the diatexite. Based on U–Pb ages, trace element and Hf isotope data, the zircon rims in the metatexite are interpreted to have crystallized from a locally derived melt, following partial dissolution of inherited protolith zircon during anatexis, whereas the zircon in the diatexite is interpreted to have crystallized from a melt that included an externally derived component. By integrating zircon and petrographic data for the migmatites and pelitic gneiss, the metatexite migmatite is interpreted to have formed by in situ partial melting in which the melt did not migrate from the source, whereas the diatexite migmatite included an externally derived juvenile component. The Cretaceous high‐temperature metamorphism of the Higo metamorphic terrane is interpreted to reflect emplacement of mantle‐derived basalts under a volcanic arc along the eastern margin of the Eurasian continent and advection of heat via hybrid silicic melts from the lower crust. Post‐peak crystallization of anatectic melts in a high‐T region at mid‐crustal depths occurred in the interval c. 116–110 Ma, as indicated by the difference in zircon ages from the metatexite and diatexite migmatites.     

Protracted zircon growth in migmatites and In situ melt of Higher Himalayan Crystallines:U-Pb ages from Bhagirathi valley,NW Himalaya,India

The Higher Himalayan Crystallines(HHC), in western Garhwal, Uttarakhand are located in a regionalscale intracontinental ductile shear zone(15-20 km wide) bounded by the Main Central Thrust at the base, and the South Tibetan Detachment System at the top. The migmatite zone in the centre has the highest grade of metamorphism in the NW Himalayas and show evidence of flowage. Zircons extracted from samples of metasediment, migmatite, biotite granite and in situ partial melt(tourmaline-bearing leucogranite) along the Bhagirathi Valley, preserve U-Pb isotopic evidence of magmatic history, magma source and effects of the Himalayan orogeny in the region. Three distinct periods of zircon growth in the leucogranite record the episodic influx of magma between 46 Ma and 20 Ma indicating a time span of more than 25 Ma between the onset of fluid-fluxed partial melting in the mid-crustal intracontinental shear zone and the emplacement of the magma into the upper crust in a post-collisional extensional setting. Metamorphic zircon growth was initiated about 46 Ma, when the partial melts were generated as the migmatite zone was exhumed.     

High-temperature cordierite migmatites in the north-eastern Grenville Province

Abstract In the northeastern part of the Grenville Province, along the gulf of St Lawrence, cordierite is widespread in the migmatites of Baie Jacques Cartier (BJC) and Baie des Ha! Ha! (BHH). In the BJC area, rafts of mesosome occur in a pervasive network of leucosome consisting of cordierite-bearing pegmatite. In BHH, however, the mesosome and leucosome are well segregated and locally separated by thin biotite –hornblende melanosomes. Leucosomes in the BJC area record the highest temperatures (oxide thermometry = 900°C), whereas leucosomes of BHH and mesosomes of both areas indicate peak temperatures around 800°C (oxide thermometry; biotite–garnet thermometry with fluorine-rich biotite). Peak pressures were constrained at 720 MPa using the Ilm-Sil–Qtz–Grt–Rt (GRAIL) equilibrium. The area is thought to have undergone extensive melting under relatively modest pressures. The highest temperatures recorded in the BJC area are probably related to a pervasive impregnation of this terrane by aluminous granitic melts. Most post-peak P–T estimates for the mesosomes fall on a nearly isobaric, clockwise, P–T path (0.6 MPa/°C) with the exception of the high-temperature leucosomes of BJC, which fall about 100°C away from this path; this is additional evidence for the external origin of these leucosomes. The ultimate source of heat that generated the migmatites is thus though to be an underlying plutonic complex (anorthosite?).     

Consequences of diffuse and channelled porous melt migration on uranium series disequilibria

Magmas erupted at mid-ocean ridges (MORB) result from decompression melting of upwelling mantle. However, the mechanism of melt transport from the source region to the surface is poorly understood. It is debated whether melt is transported through melt-filled conduits or cracks on short time scales (<∼ 10³ yrs), or whether there is a significant component of slow, equilibrium porous flow on much longer time scales (>∼ 10³-10⁴ yrs). Radiogenic excess ²²⁶Ra in MORB indicates that melt is transported from the melting region on time scales less than the half life of ²²⁶Ra (∼1600 yrs), and has been used to argue for fast melt transport from the base of the melting column. However, excess ²²⁶Ra can be generated at the bottom of the melt column, during the onset of melting, and at the top of the melt column by reactive porous flow. Determining the depth at which ²²⁶Ra is generated is critical to interpreting the rate and mechanism of magma migration. A recent compilation of high quality U-series isotope data show that in many young basalts, ²²⁶Ra excess in MORB is negatively correlated with ²³⁰Th excess. The data suggest that ²²⁶Ra excess is generated independently of ²³⁰Th excess, and cannot be explained by “dynamic” or fractional melting, where observed radiogenic excesses are all generated at the base of the melt column. One explanation is that the negative correlation of activity ratios is a result of mixing of slow moving melt that has travelled through reactive, low-porosity pathways and relatively fast moving melt that has been transported in unreactive high-porosity channels. We investigate this possibility by calculating U-series disequilibria in a melting column in which high-porosity, unreactive channels form within a low-porosity matrix that is undergoing melting. The results show that the negative correlation of ²²⁶Ra and ²³⁰Th excesses observed in MORB can be produced if ∼60% of the total melt flux travels through the low-porosity matrix. This melt maintains ²²⁶Ra excesses via chromatographic fractionation of Ra and Th during equilibrium transport. Melt that travels through the unreactive, high-porosity channels is not able to maintain significant ²²⁶Ra excesses because Ra and Th are not fractionated from each other during transport and the transport time for melt in the channels to reach the top of the melt column is longer than the time scale for ²²⁶Ra excesses to decay. Mixing of melt from the high porosity channels with melt from the low-porosity matrix at the top of the melting column can produce a negative correlation of ²²⁶Ra and ²³⁰Th excesses with the slope and magnitude observed in MORB. This transport process can also account for other aspects of the geochemistry of MORB, such as correlations between La/Yb, α_Sm/Nd, and Th/U and ²²⁶Ra and ²³⁰Th excess.     

Experimental cpx/melt partitioning of 24 trace elements

Cpx/melt partition coefficients have been determined by ion probe for 24 trace elements at natural levels in an alkali basalt experimentally equilibrated at 1,380°C and 3 GPa. One goal was to intercompare Ds for both high-field-strength elements and rare earth elements (REE) in a single experiment. Relative to the REE spidergram, Hf and Ti show virtually no anomaly, whereas Zr exhibits a major negative anomaly. Other incompatible elements (Ba, K, Nb) fall in the range of published values, as do elements such as Sr, Y, Sc, Cr and V. Pb shows a value intermediate between La and Ce. Values for Be, Li and Ga are reported for the first time, and show that Be is as incompatible as the light REEs whereas Li and Ga are somewhat more compatible than the heavy REE.     

Experimental measurements of zircon/melt trace-element partition coefficients

Zircon was grown from trace-element doped hydrous peralkaline rhyolite melts with buffered oxygen fugacities in cold-seal experiments at 0.1 and 0.2 GPa and 800 °C and piston-cylinder experiments at 1.5 GPa and 900-1300 °C. Zircon and glass were present in all run products, and small monazite crystals were present in eight of the 12 experiments. Average diameters of zircon crystals ranged from 5 to 20 μm at 800 °C to 30-50 μm at 1300 °C. Zircon crystals have thin rims, and adjacent glass has a narrow (∼1 μm thick) compositional boundary layer. Concentrations obtained through in-situ analysis of cores of run product zircon crystals and melt pools were used to calculate trace-element partition coefficients D^zircon/melt for P, Sc, Ti, V, Y, La, Ce, Pr, Nd, Eu, Gd, Ho, Yb, Lu, Hf, Th, and U. In most cases Lu was the most (D 12-105) and La the least (0.06-0.95) compatible elements. D values from this study fall within the range of previously measured values for Rare Earth Elements (REE). However, D values measured experimentally show less fractionation than those recently measured using natural phenocryst/matrix pairs. For example, D_Lu/D_La measured experimentally in this study range between 27 and 206 compared to a value of 706,522 for a natural zircon/dacite pair [Sano, Y., Terada, K., and Fukuoka, T. 2002 High mass resolution ion microprobe analysis of rare earth elements in silicate glass, apatite and zircon: lack of matrix dependency. Chem. Geol.184, 217-230]. Although D values from this study show good agreement with the lattice strain model, D values from natural phenocryst/matrix pairs combined with measured zircon compositions better reproduce host-rock (magma) compositions of igneous rocks. They also yield more reasonable estimates of magma compositions when combined with compositions of ‘‘out-of-context” zircons. For example, compositions of the Hadean detrital zircons from Jack Hills, Australia yield LREE-enriched magmas when combined with D values from phenocryst/matrix pairs yields, but yield LREE-depleted magmas when experimentally determined D values are used. We infer that experimentally measured D^zircon/melt values represent disequilibrium partitioning resulting from rapid zircon growth during short laboratory timescales. Rapid growth causes development of observed diffusive boundary layers in the melt adjacent to zircon crystals. D values from phenocryst/matrix pairs are therefore recommended for petrogenetic modeling.     

Experimental study of magmatic melt oxidation by CO2

The process of CO₂ flashing through hydrous albite-hedenbergite melt was experimentally examined at a temperature of 1100°C and a pressure of 2 kbar. Carbon dioxide was generated when the melt interacted with calcite, and wollastonite was the predominant synthesized phase. Mafic components were introduced into the hydrous albite melt via the dissolution of natural hedenbergite. Raman spectroscopic data on bubbles of the fluid phase in the quench glass indicate that the CO₂/H₂O proportions of the bubbles vary. IR spectroscopic data on the glass prove that the water concentration after CO₂ flashing decreased from 5.5 to approximately 3 wt %. The comparison of the composition of the recrystallized clinopyroxene in contact with melt (with and without CO₂ blowing) indicates that CO₂ oxidizes Fe in the melt. The redox effect of CO₂ is quantified by the empirical clinopyroxene tool for metering oxygen fugacity (oxometer), which was calibrated based on experimental data. The oxygen fugacity in our experiments with CO₂ flashing (estimated by the clinopyroxene oxometer) was NNO + (3.0?C3.5). Our estimates with the application of the clinopyroxene oxometer indicate that the maximum oxygen fugacity in the magmatic chambers of Vesuvius and Stromboli volcanoes (which are bubbled with CO₂) is also close to NNO + (3.5 ± 0.5).     

Retrograde processes in migmatites and granulites revisited

Many migmatites and granulites preserve evidence of a clockwise P–T evolution involving decompression (decrease in P) while close to the thermal peak. The extent of post‐thermal peak reaction is influenced by several factors, including: (1) the P–T path in relation to invariants in the system and the Clapeyron slopes of the equilibria; (2) the rate of cooling; and (3) the availability of fluid (H₂O‐rich volatile phase or melt) for fluid‐consuming reactions. Reaction may occur between products of a prograde (increasing T) fluid‐generating reaction as the same equilibrium is re‐crossed in the retrograde (decreasing T) sense. In general, reaction reversal or ‘back reaction’ requires the P–T path to approximate isobaric heating and cooling, without significant decompression, and evolved fluid to remain within the equilibration volume. The larger the decompression segment in the P–T evolution, the more chance there is of crossing different reactions along the retrograde segment from those crossed along the prograde segment. For common pelite compositions, we may generalize by considering three pressure regimes separated by the [Spl, Ms, H₂O] invariant in KFMASH (approximately 9 kbar) and the intersection of muscovite breakdown with the H₂O‐rich volatile phase‐saturated solidus (approximately 4 kbar). Reaction reversal cannot occur along P–T paths that traverse around one of these points, but may occur along P–T paths confined to one of the three regimes in between. Additionally, above the solidus, melt segregation and loss potentially change the composition of the equilibration volume; and, the size of the equilibration volume shrinks with decreasing T. Since the proportion of melt to residue in the equilibration volume may change with decreasing size, the composition of the equilibration volume may change throughout the supra‐solidus part of the retrograde segment of the P–T evolution. If melt has been lost from the equilibration volume, reaction reversal may not be possible or may be only partial; indeed, the common preservation of close‐to‐peak mineral assemblages in migmatite and granulite demonstrates that extensive reaction with melt is uncommon, which implies melt isolation or loss prior to crossing potential melt‐consuming reactions. Water dissolved in melt is transported through the crust to be exsolved on crystallization at the solidus appropriate to the intrinsic a(H₂O). This recycled water causes retrogression at subsolidus conditions. Consideration of the evidence for supra‐solidus decompression‐dehydration reactions, and review of microstructures that have proven controversial, such as corona and related microstructures, selvage microstructures and ‘late’ muscovite, leads to the conclusion that there is more than one way for these microstructures to form and reminds us that we should always consider multiple working hypotheses!     

What controls partial melting in migmatites?

Abstract The layers of six stromatic migmatites from Northern, Western, and Central Europe display small but systematic chemical and mineralogical differences. At least five of these migmatites do not show any signs of largescale metamorphic differentiation, metasomatism, or segregation of melts. It is concluded, therefore, that the compositional layering observed in most of the investigated migmatites is due to compositional differences inherited from the parent rocks. Almost isochemical partial melting seems to be the most probable process transforming layered paragneisses, metavolcanics, or schists into migmatites.
The formation of neosomes is believed to be caused by higher amounts of partial melts formed due to higher amounts of water moving into these layers. The neosomes have less biotite and more K-feldspar, if K-feldspar is present at all, than the adjacent mesosomes. These differences are small but systematic and seem to control the access of different amounts of water to the various rock portions. Petrographical observations, chemical data, and theoretical considerations indicate a close relationship between rock composition, rock deformation, transport of water, partial melting, and formation of layered migmatites.     

Experimental investigation of contaminant migration in filled fracture

Experiment equipments involved in the single tube fracture and double tube fracture models are designed to research the characteristics of groundwater flow and solute transport in filled fracture. During the experiment, the state of groundwater flow can be characterized as linear flow, and satisfies Darcy’s law. Therefore, based on the pipe flow of hydraulics and Darcy’s law, the flow rate and water flow velocity can be calculated. Also, dispersion parameters were calculated with the fitting of observed data and analytical solution in the single tube fracture model. Furthermore, effects of some factors on solute transport are involved in the double tube fracture model, and length of branch fracture, particles’ diameter and flow rate in water inlet have been discussed. Results show that the arrival time of concentration peak value in the single tube fracture model is faster than that in the double tube fracture model, and two concentration peak values exist in the double tube fracture model. Arrival time of concentration peak value is faster with the increase of branch fracture length. Furthermore, if the branch fracture is longer, arrival time of the first concentration peak value is faster, while arrival time of the second concentration peak value is slower, relative to short branch fracture.