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.     

Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol + hy-normative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative sequence is exposed in the Dents de Bertol area (center of intrusion). P–T calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 °C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#_cpx = 83–89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment.     

拉萨地体内低速异常成因的探讨--壳内局部熔融的地震层析证据

青藏高原南部天然地震体波的层析反演图像描绘出拉萨地体内在15～90km深度内低速异常的分布特征。在研究速度结构的基础上，根据熔融程度与纵波波速呈倒数衰减的特性，论证了拉萨地体中低速体的出现与壳内的局部熔融有关。壳内局部熔融则来源于板块俯冲、地壳增厚产生的热而非幔源侵入所造成。同时，也发现了在拉萨以东中下地壳低速体有向东发展的趋势。大面积部分熔融层分布的深度范围约在15～35km处。另外，据藏南地区的SKS波分裂时间延迟很小的特点，推测藏南地幔中可能不存在各向异性，也就是没有足够的热源使各向异性矿物的品格重新定位。这从另一方面佐证了部分熔融仅存在于地壳内而不是地幔中。     

Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope     

Continuous Gradations among Primary Carbonatitic, Kimberlitic, Melilititic, Basaltic, Picritic, and Komatiitic Melts in Equilibrium with Garnet Lherzolite at 3-8 GPa

Multianvil melting experiments in the system CaO–MgO–Al₂O₃–SiO₂–CO₂(CMAS–CO₂) at 3–8 GPa, 1340–1800°C, involvingthe garnet lherzolite phase assemblage in equilibrium with CO₂-bearingmelts, yield continuous gradations in melt composition betweencarbonatite, kimberlite, melilitite, komatiite, picrite, andbasalt melts. The phase relations encompass a divariant surfacein P–T space. Comparison of the carbonatitic melts producedat the low-temperature side of this surface with naturally occurringcarbonatites indicates that natural magnesiocarbonatites couldbe generated over a wide range of pressures >2·5 GPa.Melts analogous to kimberlites form at higher temperatures alongthe divariant surface, which suggests that kimberlite genesisrequires more elevated geotherms. However, the amount of waterfound in some kimberlites has the potential to lower temperaturesfor the generation of kimberlitic melts by up to 150°C,provided no hydrous phases are present. Compositions resemblinggroup IB and IA kimberlites are produced at pressures around5–6 GPa and 10 GPa, respectively, whereas the compositionsof some other kimberlites suggest generation at higher pressuresstill. At pressures <4 GPa, an elevated geotherm producesmelilitite-like melt in the CMAS–CO₂ system rather thankimberlite. Even when a relatively CO₂-rich mantle compositioncontaining 0·15 wt % CO₂ is assumed, kimberlites andmelilitites are produced by <1% melting and carbonatitesare generated by even smaller degrees of melting of <0·5%. KEY WORDS: carbonatite; CO₂; kimberlite; melilitite; melt generation     

Origin of Exceptionally Abundant Phonolites on Ua Pou Island (Marquesas, French Polynesia): Partial Melting of Basanites Followed by Crustal Contamination

On the basis of the first systematic mapping of Ua Pou, longknown for its exceptionally abundant phonolites, we estimatethat these rocks cover 65% of the surface of the island whereasmafic lavas cover 27% and intermediate ones 8%. The silica-undersaturatedsuite was erupted in a restricted time span (2·9–2·35Myr), following the emplacement of tholeiites derived from ayoung HIMU-type source at c. 4 Ma. Primitive basanites, derivedfrom a heterogeneous mantle source with a dominant EM II + HIMUsignature, represent likely parental magmas. The series is characterizedby a Daly gap defined by a lack of phonotephrites. We considerthat the most likely model for the origin of evolved lavas ispartial melting at depth of primitive basanites, leaving anamphibole-rich residuum and producing tephriphonolitic magmas.These tephriphonolitic magmas may have evolved by closed-systemfractional crystallization towards Group A phonolites. Threeother groups of phonolites could have been derived from tephriphonoliticmagmas by open-system fractional crystallization processes,characterized respectively by seawater contamination (GroupB), assimilation of nepheline syenite-type materials (GroupC) and extreme fractionation coupled with assimilation of theunderlying oceanic crust (Group D). The prominence of evolvedlavas is a consequence of their origin from partial meltingof mafic precursors followed by crustal contamination. KEY WORDS: Marquesas; French Polynesia; phonolite; partial melting; contamination     

Experimental determination of argon solubility in silicate melts: An assessment of the effects of liquid composition and temperature

The argon solubility of 38 liquids in the system Na₂O-CaO-MgO-Al₂O₃-SiO₂ (NCMAS) has been determined at 1873 K and 1 bar, the argon concentration of presaturated glasses being measured using a static mass spectrometer. For compositions in the subsystem diopside (CaMgSi₂O₆), nepheline (NaAlSiO₄), albite (NaAlSi₃O₈), anorthite (CaAl₂Si₂O₈), argon solubility is generally a linear function of the relative proportion of each end member, solubility being lowest in diopside melt (1.53 10⁻⁵ cm³ STP · g⁻¹ · bar⁻¹) and highest in albite melt (2.88 10⁻⁴ cm³ STP · g⁻¹ · bar⁻¹). For the tectosilicate joins studied (SiO₂-Na₂Al₂O₄, SiO₂-CaAl₂O₄, SiO₂-MgAl₂O₄) solubility decreases with decreasing silica content in all cases, being highest for Na-bearing liquids and lowest for Mg-bearing liquids at constant molar silica content. Where comparison is possible our results are in good agreement with data from the literature. When our data are considered in isolation we find that argon solubility shows an excellent correlation with calculated ionic porosity. The covariation of argon solubility and liquid density is also reasonable, that with molar volume less convincing and that with polymerization state (as defined by the ratio of the number of nonbridging oxygens and tetrahedral network forming cations; NBO/T) nonexistent. However, when our data are combined with those from the literature no well constrained correlation between argon solubility and ionic porosity is apparent. Based upon this observation and consideration of the temperature dependence of noble gas solubility it is concluded that ionic porosity is not a universally applicable parameter which may be used to predict noble gas solubility as a function of composition, temperature and pressure. Two new models for calculating argon solubility are proposed, both employing the notion of partial molar argon solubilities. The first uses oxide components, for which partial molar argon solubility is directly proportional to partial molar ionic porosity calculated at 1873 K, irrespective of the temperature of experimental equilibration. The second model, which offers the best fit to the available data, employs tetrahedral units rather than oxides as the proposed melt components. This latter model successfully accounts for reported argon solubilities in simple Al-free systems, in simple Al-bearing systems and in natural liquids. This is interpreted to infer that argon is incorporated in large sites in the liquid structure (such as the space within rings of n-tetrahedra) although further work is required to understand the quantitative links between melt structure and noble gas solubility.     

Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.     

Micromechanical approach to swelling behavior of capillary-porous media with coupled physics

A detailed multiscale analysis is presented of the swelling phenomenon in unsaturated clay-rich materials in the linear regime through homogenization. Herein, the structural complexity of the material is formulated as a three-scale, triple porosity medium within which microstructural information is transmitted across the various scales, leading ultimately to an enriched stress-deformation relation at the macroscopic scale. As a side note, such derived relationship leads to a tensorial stress partitioning that is reminiscent of a Terzaghi-like effective stress measure. Otherwise, a major result that stands out from previous works is the explicit expression of swelling stress and capillary stress in terms of micromechanical interactions at the very fine scale down to the clay platelet level, along with capillary stress emerging due to interactions between fluid phases at the different scales, including surface tension, pore size, and morphology. More importantly, the swelling stress is correlated with the disjoining forces due to electrochemical effects of charged ions on clay minerals and van der Waals forces at the nanoscale. The resulting analytical expressions also elucidate the role of the various physics in the deformational behavior of clayey material. Finally, the capability of the proposed formulation in capturing salient behaviors of unsaturated expansive clays is illustrated through some numerical examples.     

Constitutive modelling of fine-grained gassy soil: A composite approach

Fine-grained marine sediments containing large undissolved gas bubbles are widely distributed around the world. Presence of the bubbles could degrade the undrained shear strength (s_u ) of the soil, when the gas pressure u_g is relatively high as compared with the effective stress in the saturated soil matrix. Meanwhile, the addition of bubbles may also increase s_u when the difference between u_g and pore water pressure u_w becomes smaller than the water entry value, causing partial water drainage from the saturated matrix into the bubbles (bubble flooding) during globally undrained shearing. A new constitutive model for describing the two competing effects on the stress-strain relationship of fine-grained gassy soil is proposed within the framework of critical state soil mechanics. The gassy soil is considered as a three-phase composite material with compressible cavities, which allows water entry from the saturated matrix. Bubble flooding is modelled by introducing an additional positive volumetric strain increment of the saturated clay matrix, which is dependent on the difference between pore gas and pore water pressure based on experimental observations. A modified hardening law based on that of the modified Cam clay model is employed, which in conjunction with the expression for bubble flooding, can describe both the detrimental and beneficial effects of gas bubbles on soil strength and plastic hardening in shear. Only two extra parameters in addition to those in the modified Cam clay model are used. It is shown that the key features of the stress-strain relationship of three fine-grained gassy soils can be reproduced satisfactorily.