Small Melt Fractions Research Articles

Stabilization and evolution of the Brazilian subcontinental lithospheric mantle: Insights from garnet xenocrysts and peridotite xenoliths of Três Ranchos kimberlite (APIP, Brazil)

The thermal structure and chemical composition of the subcontinental lithospheric mantle beneath the Alto Paranaíba Igneous Province defined by the major and trace element compositions of garnet and clinopyroxene xenocrysts and mantle peridotite xenoliths entrained in the Cretaceous Três Ranchos kimberlite. Estimated temperatures and pressures indicate a typical cratonic geotherm, with a surface heat flow of ∼ 40 mW/m2 slightly hotter than that at the São Francisco Craton margin (35–36 mW/m2), suggesting lateral heterogeneities in the geothermal gradients. A deep mantle portion (160 km) is recorded by a lherzolite xenolith (5 GPa and 1200 °C). This portion has undergone pervasive interstitial melt percolation of an ultra-alkaline silicate melt and late injection of small silicate melt fractions such as the host carbonated kimberlite. This pervasive migration extended into shallower SCLM sectors at 124–100 km (3.9–3.1 GPa and ∼ 940–810 °C) as evidenced by metasomatized (±phlogopite) peridotite identified by sinusoidal REE garnet compositions. This portion is apparently intercalated at 120–110 km depth (3.7–3.3 GPa and ∼ 850 °C) with strongly depleted peridotites showing garnet and clinopyroxene compositions extremely impoverished in moderately incompatible elements, as a result of extreme ancient, melt extraction. Atop the metasomatized and depleted peridotites, at 100 km depth, a layer of fertile lherzolite is identified by garnet xenocrysts, with 'normal', HREE-enriched, steadily fractionated LREE-depleted patterns. Combining our data with those of the Osvaldo França SCLM points to the occurrence of a deeper sector (>100 km depth beneath Trȇs Ranchos) originally formed by extremely depleted peridotites, successively metasomatized by ultra-alkaline melts; and a shallower sector of the São Francisco craton composed of fertile lherzolites unrelated to ultra-alkaline magmatism. Similar fertile lherzolites are found in the uppermost sectors of the Siberian and Kaapval SCLM, suggesting that they may be a ubiquitous, primary feature, possibly related to the stabilization of the Archean cratons.

Slab window–related magmatism as a probe for pyroxenite heterogeneities in the upper mantle

Abstract New high-precision trace-element analyses of magmatic olivines point to a pyroxenite-dominated source for recent alkali basalts erupted above slab windows formed along the Antarctic Peninsula. Melting occurred at ambient mantle temperature, and basalts have geochemical compositions that are indistinguishable from ocean-island basalts (OIBs). We propose that the pyroxenite component originally resided in the upper mantle beneath the subducted slab; formation of a slab window allowed limited decompression and the generation of melts of garnet-pyroxenite, but little or no melting of mantle peridotite. The pyroxenite component in the mantle formed ca. 550 Ma, an age that does not require long-term recycling of subducting slabs to the core-mantle boundary. Enriched mid-ocean ridge basalt (E-MORB) from the adjacent extinct Phoenix Ridge owes its enriched trace-element compositions to mixing between small melt fractions of pyroxenite and peridotite during a period of decreased spreading rate prior to the death of the ridge ca. 3.3 Ma. It is likely that the variable trace-element enrichment seen in East Pacific Rise E-MORB distal from hotspots results from the same process of interactions between small-melt-fraction (<~5%) melts of pyroxenite and peridotite.

Melting, Compaction and Reactive Flow: Controls on Melt Fraction and Composition Change in Crustal Mush Reservoirs

Abstract Changes in melt fraction and local bulk composition in high-crystallinity, crustal mush reservoirs are essential to produce the large volumes of low-crystallinity, silicic magma that are emplaced to form plutons, or erupted to surface. Heating (and cooling) is well understood and widely invoked in driving melt fraction change, but does not cause chemical differentiation because there is no separation of melt and crystals. Fractional crystallisation at high melt fraction is widely assumed to explain differentiation, but is inconsistent with the evidence that large-scale, long-term magma storage and evolution occurs in high-crystallinity mush reservoirs. Compaction has been suggested to explain melt fraction change and differentiation at low melt fraction, but compaction (and decompaction) causes simple unmixing (and mixing) of melt and solid crystals: to produce very refractory bulk composition by compaction, melt fraction must be driven down to very low values. Yet microstructural evidence demonstrating widespread compaction in crustal mush reservoirs at low melt fraction is lacking. Here we show that melt fraction change can be expressed in terms of heating/cooling and compaction, plus an additional term that we call ‘reactive flow’. Similarly, composition change can be expressed in terms of compaction and reactive flow. Reactive flow changes the local bulk composition, which causes ‘chemical’ melting (dissolution) and freezing (precipitation), distinct from ‘thermal’ melting/freezing caused by changes in enthalpy. We use numerical modelling to show that the contributions of compaction and reactive flow in a crustal magma reservoir are similar in magnitude. However, reactive flow opposes melt fraction and composition changes caused by compaction when compaction occurs in a temperature gradient that increases upwards at, for example, the base of a sill intrusion, or decompaction occurs in a temperature gradient that decreases upwards at, for example, the top of a sill intrusion. If compaction causes melt fraction decrease and creates a more refractory bulk composition, then reactive flow causes melt fraction increase and a more evolved bulk composition, and vice versa. Reactive flow means that very small melt fraction is not required to produce very refractory composition in a crustal magma reservoir, consistent with the relatively scarce microstructural evidence for widespread compaction. The apparent lack of compaction in crustal magma reservoirs, as compared to other natural and engineered systems in which reaction does not occur, is also explained by the contribution of reactive flow. Reactive flow also means that melt loss in compacting regions of a crustal magma reservoir may instead be accompanied by evidence for mineral dissolution, which facilitates ongoing melt fraction loss by preserving connected melt flow paths through the mush pore space. Reactive flow in decompacting regions can also explain why interstitial mineral phases display textures that mimic those of interstitial melt.

Rare elements enrichment in crustal peraluminous magmas: insights from partial melting experiments

Experiments were conducted to explore the behavior of Li, Rb, Nb, Sn, Cs, Ta, W during crustal melting and test the anatectic origin of rare metal-bearing peraluminous granites such as rare metal granites (RMGs). The experiments were performed under fluid-absent conditions at 800 and 850 °C, 400 MPa and moderately reducing fO2 (ΔFMQ = − 0.5 to − 0.8). Starting materials were cores of several millimetres drilled from two natural rocks, a biotite-rich paragneiss (Pg) and a muscovite-rich orthogneiss (Og) enriched in Li, Be, Sn, Cs, W. Both protoliths produced small melt fractions from 8 to 20% vol. Melt distributions were either homogeneously distributed at grain boundaries in the Pg or preferentially associated with muscovite reaction zones in the Og. In the Pg at 800 °C, melting is mainly fluid present, driven by interstitial water at grain boundaries. At 850 °C, biotite dehydration-melting produces peritectic orthopyroxene, hercynitic spinel, ilmenite and alkali feldspar in addition to melt. In the Og, muscovite dehydration-melting generates melt plus peritectic biotite, hercynitic spinel, ilmenite, Al silicates and alkali feldspar. Experimental glasses are nearly homogeneous, silica rich, peraluminous and leucogranitic and their major element compositions differ only little between the two protoliths. In contrast, the trace element concentrations vary as a consequence of chemical and textural heterogeneities in our starting materials. Compared with source rocks, the Og glasses are enriched in Rb, Nb, Ta, W and depleted in Li, Cs and the Pg are enriched in Li, Rb, Cs, W and depleted in Nb, Ta. Mass-balance calculations indicate that during muscovite dehydration-melting, Li, Cs and Rb partition into the melt; whereas Nb, Ta and W are preferentially incorporated in peritectic phases. Li and Cs also partition toward the melt during biotite dehydration-melting. The partitioning behavior of trace elements during crustal melting is a function of the melting reaction and partition coefficients between melt, residual and peritectic phases. Experimental glasses are similar to peraluminous muscovite granites but fail to reproduce RMG compositions. Alternatives to mica dehydration-melting such as fluid-present and residual source melting emphasize the difficulties with an origin of RMGs by purely anatectic processes. Crystallization differentiation might have to be combined with mica dehydration-melting to explain the distinctive geochemical features of RMGs.

Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure: Implications for stability and seismic signature of Fe-rich melts in the mantle

Iron-rich silicate melts play an important role in the magmatic history of the Earth and the Moon. However, their elastic properties at high pressures, especially the sound velocities, are poorly understood. Here we determined the ultrasonic sound velocities of two silicate melts along the hedenbergite (Hd, CaFeSi2O6) – diopside (Di, CaMgSi2O6) join at high pressure and temperature conditions up to 6 GPa and 2329 K, using the high-pressure ultrasonic technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that Fe can significantly reduce the sound velocity while increasing the density of silicate melts. Comparing the melts of Di, Hd, and a mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50), we find that although densities of the Hd-Di melts can be well-described by a linear mixing law empirically at high pressures, sound velocities do not vary with composition linearly. Assuming that the low-velocity regions in Earth's upper mantle are mainly due to the presence of partial melt, we applied our results to study the gravitational stability and seismic signature of Fe-rich silicate melts in the mantle, combined with melt geometry and compaction models. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of seismic velocity reduction can be explained by the presence of a small amount of partial melt distributed in film/band geometry along grain boundaries, silicate melts formed in this depth range are too light to be gravitationally stable, but may be completely entrained in the convecting mantle due to the small melt fraction and low melt-matrix separation velocity. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-rich melts (with FeO>∼10 wt%) distributed in textural equilibrium with the ambient mantle is a plausible explanation.

Lithospheric mantle refertilization by DMM-derived melts beneath the Cameroon Volcanic Line\u2014a case study of the Befang xenolith suite (Oku Volcanic Group, Cameroon)

The origin and evolution of subcontinental lithospheric mantle (SCLM) are important issues of Earth’s chemical and physical evolution. Here, we report detailed textural and chemical analyses on a mantle xenolith suite from Befang (Oku Volcanic Group, Cameroon Volcanic Line), which represents a major tectono-magmatic structure of the African plate. The samples are sourced from spinel-facies mantle and are dominated by lherzolites. Their texture is cataclastic to porphyroclastic, and foliation defined by grain-size variation and alignment of spinel occurs in part of peridotites. Spinel is interstitial and has amoeboidal shape. Clinopyroxene REE patterns are similar to those of Depleted MORB Mantle (DMM) except LREEs, which vary from depleted to enriched. The A-type olivine fabric occurs in the subset of one harzburgite and 7 lherzolites studied by EBSD. Orthopyroxene shows deformation consistent with olivine. The fabric of LREE-enriched clinopyroxene is equivalent to those of orthopyroxene and olivine, whereas spinel and LREE-depleted clinopyroxene are oriented independently of host rock fabric. The textural, chemical and thermobarometric constraints indicate that the Befang mantle section was refertilised by MORB-like melt at pressures of 1.0–1.4 GPa and temperatures slightly above 1200–1275 °C. The olivine-orthopyroxene framework and LREE-enriched clinopyroxene preserve the protolith fabric. In contrast, the LREE-depleted clinopyroxene, showing discordant deformation relative to the olivine-orthopyroxene protolith framework, and amoeboidal spinel crystallized from the infiltrating melt. The major element and REEs composition of minerals forming the Befang peridotites indicate subsequent reequilibration at temperatures 930–1000 °C. This was followed by the formation of websterite veins in the lithospheric mantle, which can be linked to Cenozoic volcanism in the Cameroon Volcanic Line that also brought the xenoliths to the surface. This study therefore supports the origin of fertile SCLM via refertilization rather than by extraction of small melt fractions, and further emphasizes the involvement of depleted melts in this process.

Application of Machine Learning (ML) for Enhancing the Transient Performance of Thermal Energy Storage (TES) Platforms Using Radial Basis Function (RBF)

Thermal energy storage (TES) can be utilized as supplemental platforms for improving operational reliability and systemic efficiency in variety of industries, such as for reducing water usage in power production (food-energy-water/ FEW nexus), chemical and agro-process industries and for improving sustainability (e.g., desalination), etc. Phase change materials (PCMs) can be used in TES due to their high latent heat storage capacity during phase transformation. Inorganic PCMs typically have the highest latent heat capacity and are attractive for their ability to store the larger quantities of thermal energy in small form factors while conferring respectable power ratings (however, they suffer from compromised reliability issues, that often arise from the need for subcooling). Subcooling (also known as supercooling) is a phenomenon where the temperature needs to be reduced substantially below the melting point to initiate solidification. A technique for obviating subcooling issues is to allow a small portion of the PCM to remain un-melted. This allows the PCM to initiate nucleation from the un-melted portion of PCM (this is termed as the “cold finger” technique). Thus, reliability is enhanced at the expense of substantial reduction in storage capacity. A fundamental challenge for using this technique is the inability to reliably predict and control the amount of melt fraction in the total volume of the PCM (such that a target amount of the PCM remains solidified or un-melted at the end of each melt-cycle during repeated melting and solidification of the total mass of PCM). However, using Machine Learning (ML) techniques, this deficiency can be addressed by reliably predicting and thus controlling the amount of melt fraction in the total volume of the PCM with a higher accuracy than conventional techniques (such as using multi-physics-based models or numerical solvers). Conventional techniques for predicting transient characteristics in real time control schemes typically leverage multi-physics-based models that are often effective only for a narrow range of operating conditions with concomitant disadvantages: they are highly sensitive to small variations in the measurement uncertainties and are therefore susceptible to large levels of error in the real time predictions (and are unreliable for implementation in diverse range of operating conditions). In this pioneering study, nearest neighbor search processes (such as radial basis functions) were utilized along with machine learning (ML) algorithm using a training data set to predict the PCM melt fraction and to demonstrate the feasibility (and efficacy) of this approach. This technique is simple to implement and is device independent as well as robust (i.e., it can be deployed successfully even under conditions where the sensors malfunction, such as thermocouples that are off-calibration). This technique was demonstrated successfully for predicting the melt fraction of a PCM with high accuracy and robustness. With this method, the melt fraction of a PCM can be accurately determined, which allows the maximum thermal capacity of a PCM to be utilized while mitigating reliability issues (such as subcooling) and enhancing the thermodynamic efficiencies of the TES platforms. Melting experiments were performed using a digital camera (for video recording) and a graduated cylinder containing PCM for monitoring the transient values of the melt fraction based on the height of the liquid phase of the PCM in the cylinder. An array of 3 thermocouples was mounted at specific heights within the body of the PCM to monitor the temperature transients at these specific location during the propagation of the melt front within the PCM. In the final stages of the melting process, the predictions from the ML algorithm was found to be more accurate (90~95% accuracy) than that of the conventional techniques based on physics-based solvers (~60% accuracy). The accuracy of the ML algorithm was low at smaller melt fractions (~30%) and improved substantially at higher melt fractions (~95%). Furthermore, the accomplishments of this study display the feasibility of a RBF ML method which can be implemented for the accurate prediction and control of a real world stochastic system which can exhibit nonlinear and chaotic dynamics which change over time.

Chalcophile-siderophile element systematics and regional-scale magmatic percolation in the Ronda peridotite massif (Spain)

Eighteen samples from the south-western part of the Ronda lherzolitic massif have been analysed with bulk-rock and in-situ techniques to unravel the fate of Cu, platinum-group element and chalcogens (S, Se) in a middle Proterozoic subcontinental lithospheric mantle segment variably overprinted by kilometre-scale porous flow percolation of asthenosphere-derived silicate melts during the Alpine orogeny. Chalcogen elements (S, Se Cu) and PGE systematics fit well the published data base for orogenic mantle peridotites as a whole. Positive correlations between S, Se, Cu and fertility result from progressive removal of sulfides from the mantle at increasing degrees of melting. A few harzburgites preserved the S-, Se- and Pd-depleted (Pd/Ir < 1)-chondrite normalized PGE patterns of residues after partial melting. The Oligocene thermo-mechanical erosion event triggered secondary partial melting of the sulfide phase liberating a Cu-Ni-rich sulfide melt now identified through Cu-Ni-rich sulfide inclusions (pentlandite + chalcopyrite ± bornite) and a wide range of intergranular Cu-rich sulfides and alloys. Residual monosulfide solid solution (coexisting with refractory platinum group minerals (laurite-erlichmanite) in spinel tectonites) survived throughout all of the tectonometamorphic domains, thus preserving relative and absolute abundances of compatible platinum group elements (Os, Ir, Ru, Rh). Local scale redistribution of the Cu-Ni-rich sulfide melt accounts for the enrichment/depletion trends of Te, Pd, Pt, Au, Ag, Cu identified in granular peridotites. On cooling, the sulfide melt produced pentlandite by reacting with monosulfide solid solution at T < 870 °C while transfering incompatible elements that generated micrometric Pt-Cu-Te-As-bearing platinum-group mineral inclusions (e.g. malanite CuPt2IrS4, moncheite PtTe2, Pt arsenide, secondary Pt–Cu alloys), in addition to superchondritic Pd/Ir and nearly chondritic Se/Te identified in mantle sulfides of metasomatic origin. By contrast, micrometer-sized Au particles identified in Au-enriched spinel harzburgites are clearly primary precipitates from volatile-rich small melt fractions.

The effect of effective rock viscosity on 2-D magmatic porosity waves

Abstract. In source regions of magmatic systems the temperature is above solidus, and melt ascent is assumed to occur predominantly by two-phase flow, which includes a fluid phase (melt) and a porous deformable matrix. Since McKenzie (1984) introduced equations for two-phase flow, numerous solutions have been studied, one of which predicts the emergence of solitary porosity waves. By now most analytical and numerical solutions for these waves used strongly simplified models for the shear- and bulk viscosity of the matrix, significantly overestimating the viscosity or completely neglecting the porosity dependence of the bulk viscosity. Schmeling et al. (2012) suggested viscosity laws in which the viscosity decreases very rapidly for small melt fractions. They are incorporated into a 2-D finite difference mantle convection code with two-phase flow (FDCON) to study the ascent of solitary porosity waves. The models show that, starting with a Gaussian-shaped wave, they rapidly evolve into a solitary wave with similar shape and a certain amplitude. Despite the strongly weaker rheologies compared to previous viscosity laws, the effects on dispersion curves and wave shape are only moderate as long as the background porosity is fairly small. The models are still in good agreement with semi-analytic solutions which neglect the shear stress term in the melt segregation equation. However, for higher background porosities and wave amplitudes associated with a viscosity decrease of 50 % or more, the phase velocity and the width of the waves are significantly decreased. Our models show that melt ascent by solitary waves is still a viable mechanism even for more realistic matrix viscosities.

Extremely young melt infiltration of the sub-continental lithospheric mantle

It has long been inferred that mantle metasomatism and the incompatible element enrichment of the continents both require movement of melts formed by very low degree melting of the mantle. Yet establishing the presence of these melts and whether this process is on-going and continuous, or spatially and temporally restricted, has proved difficult. Here we report large U-Th-Ra disequilibria in metasomatised, mantle xenoliths erupted in very young lavas from the Newer Volcanics Province in southeastern Australia. The 226Ra-230Th disequilibria appear to require reappraisal of previous estimates for the age of eruption that now seems unlikely to be more than a few kyr at most. We propose that infiltration of carbonatitic melts/fluids, combined with crystallization of pargasite, can account for the first order U-series disequilibria observations. Irrespective of the exact details of the complex processes responsible, the half-lives of the nuclides require that some of the chemical and isotopic disturbance was extremely young (« 8 kyr) and potentially on-going at the time of incorporation into the alkali basalts that transported the xenoliths to the surface. This provides evidence for the presence and possibly continuing migration of small melt fractions (~0.02%) in the upper convecting mantle that may contribute to the seismic low velocity zone. By implication, it appears that the asthenosphere must lie close to its solidus, at least in this region. Pressure-temperature estimates indicate that the small degree melts identified could infiltrate as far as 25 km upwards into the sub-continental lithospheric mantle leading to strong incompatible element enrichment and the recent timing of this event this urges a reappraisal of the meaning of 300–500 Ma Nd model ages in mantle xenoliths from this region. In principle, the resultant metasomatised mantle could provide a component for some ocean island basalts, should the sub-continental lithospheric mantle be returned to the asthenosphere by convective removal at some later time.

Deformation, Annealing, Melt‐Rock Interaction, and Seismic Properties of an Old Domain of the Equatorial Atlantic Lithospheric Mantle

AbstractWe analyzed a set of mantle xenoliths from the Fernando de Noronha archipelago to constrain the roles of deformation, cooling, annealing, and melt percolation in the evolution of the mantle lithosphere of the equatorial Atlantic. The peridotites are dominantly lherzolites with coarse granular or porphyroclastic microstructures. Equilibrium temperatures range between 850 and 1,000 °C. Olivine crystal preferred orientations (CPO) have mainly orthorhombic patterns characterized by clear [100] and [010] point maxima or a fiber‐[100] tendency; these patterns imply deformation by dislocation creep with dominant activation of the [100](010) slip system. Olivine fiber‐[010] patterns are less common; they probably result from recrystallization or melt‐rock interaction. Annealing after deformation partially recovered the microstructures and reduced the olivine CPO strength. Pyroxene CPOs often have weak consistency with olivine CPO, implying post‐deformation refertilization. Chemical compositions are systematically more fertile than those from abyssal peridotites from the Mid‐Atlantic Ridge. These compositions together with the microstructures suggest two stages of melt‐rock interaction. An early refertilization of the base of the lithosphere due to continuous percolation of small melt fractions and later, more important, but local chemical changes (both refertilization and dunitization) as well as recrystallization due to reactive melt percolation likely associated with the Cenozoic volcanism. By comparing the calculated seismic properties with existing seismological data we argue that the moderate seismic anisotropy displayed by the Fernando de Noronha mantle xenoliths likely records past horizontal asthenospheric flow parallel to the spreading direction frozen in the lithospheric mantle as the plate cooled.

Influence of water on granite generation: Modeling and perspective

Water influences the physics and chemistry of magmatic differentiation because it reduces the melting point, decreases melt viscosity, modifies phase equilibria, and controls how latent heat is released. Here, we explore how compositional trends of fractionating magmas can be used to evaluate the amount of water in the crystallizing system. Of interest is the role of water in generating silicic magmas such as granites and rhyolites. Water expands the crystallization window over which silicic melts are stable, such that granites are limited to small residual melt fractions (F < 0.2) under dry conditions but stable over a wider range of residual melt fractions (F < 0.6) at high water contents. Variations of SiO2 versus F are relatively insensitive to oxygen fugacity and can be used as a hygrometer for granitoids, where F is estimated through relative enrichments in incompatible elements, such as Th and K. We also show that generation of granitoids in arc environments follow distinctly different Mg#-SiO2 trends compared with anorogenic or intra-plate granites. We confirm that arc magmas differentiate at higher water contents and possibly higher oxygen fugacities than anorogenic/intra-plate magmas. Finally, Archean TTGs (tonalite, trondhjemite, and granodiorite) show similar Mg#-SiO2 systematics as Phanerozoic arc-related granitoids, suggesting similar petrogenetic physical and chemical conditions.

Melt organisation and strain partitioning in the lower crust

Partial melts can form as a result of crustal thickening due to orogenesis. Even small melt fractions weaken the crust, so that partially molten volumes should accumulate significant amounts of strain. However, relatively little is known of how strain partitions in partial melts, and how effective the melt expulsion processes from the partially molten crust are. Using examples from the Western Gneiss Region (WGR), Norway, we consider a case of co-existing migmatites and shear zones. Field, image analysis, and microanalytical methods allow (semi)quantification of melt volume, rock mineralogy and mineral chemistry, and microstructures. Integration of these analyses implies effective syn-melt strain partitioning and subsequent freezing of both the shear zone and migmatite texture. We propose a mechanism that allows i) syn-melt strain localisation at an outcrop scale through stress-driven melt organisation, resulting in significant relative competence differences in a partially molten rock volume; and ii) formation of fine-grained rocks at outcrop that is entirely or mostly syn-melt, without subsequent mylonitic shearing in the solid-state. Syn-melt shear zones that have not acted as effective melt transport channels and/or that have not accumulated post-melt deformation may be more common than conventionally assumed.

Melting of Peridotites through to Granites: A Simple Thermodynamic Model in the System KNCFMASHTOCr

A new set of thermodynamic models is presented for calculating phase relations in bulk compositions extending from peridotite to granite, from 0·001 to 70 kbar and from 650°C to peridotite liquidus temperatures, in the system K2O-Na2O-CaO-FeO-MgO-Al2O3-SiO2-H2O-TiO2-Fe2O3- Cr2O3 (KNCFMASHTOCr). The models may be used to calculate phase equilibria in partial melting of a large range of mantle and crustal compositions. They provide a good fit to experimental phase relation topologies and melt compositions across the compositional range of the model. Compared with the preliminary model of Jennings, E. S. & Holland, T. J. B. (2015) (A simple thermodynamic model for melting of peridotite in the system NCFMASOCr. Journal of Petrology 56, 869-892) for peridotite-basalt melting relations, the inclusion of K2O and TiO2 allows for better modelling of small melt fractions in peridotite melting, and in reproducing rutile-bearing eclogite melting at high pressures. An improved order-disorder model for spinel is now incorporated. Above 10 kbar pressure, wet partial melting relations may be significantly affected by the dissolution of silicates in aqueous fluid, so the set of models includes an aqueous low-density silicatebearing fluid in addition to a high-density H2O-bearing silicate melt. Oxygen fugacity may be readily calculated for the whole range of bulk compositions investigated, and the effect of water content on melt fO2 is assessed.

He, Ne and Ar ‘snapshot’ of the subcontinental lithospheric mantle from CO2 well gases

The subcontinental lithospheric mantle (SCLM) constitutes a significant portion of the upper mantle sourcing magmatic volatiles to the continents above, yet its geochemical signature and evolution remain poorly constrained. Here we present new interpretation of noble gas datasets from two magmatic CO2 fields in the SW US, namely Bravo Dome and Sheep Mountain, which provide a unique insight into the volatile character of the SCLM sourcing the Cenozoic volcanism in the region. We identify that reduction of 3He/4Hemantle ratio within the Sheep Mountain CO2 field can be attributed to radiogenic production within the SCLM. Using a Reduced Chi-Squared minimisation on the variation of derived 4He/21Necrust ratios within samples from the Sheep Mountain field, combined with a radiogenically raised 21Ne/22Nemantle end member, we resolve 3He/4Hemantle ratios of 2.59 ± 0.15 to 3.00 ± 0.18 Ra. These values correspond with a 21Ne/22Nemantle value of 0.136. Using these 3He/4Hemantle end member values with 21Nemantle resolved from Ne three component analysis, we derive the elemental 3He/22Nemantle of 2.80 ± 0.16 and radiogenic 4He/21Ne*mantle range of 1.11 ± 0.11 to 1.30 ± 0.14. A second Reduced Chi-Squared minimisation performed on the variation of 21Ne/40Arcrust ratios has allowed us to also determine both the 4He/40Armantle range of 0.78 to 1.21 and 21Ne/40Armantle of 7.66 ± 1.62 to 7.70 ± 1.54 within the field. Combining these ratios with the known mantle production ranges for 4He/21Ne and 4He/40Ar allows resolution of the radiogenic He/Ne and He/Ar ratios corresponding to the radiogenically lowered 3He/4Hemantle ratios. Comparing these values with those resolved from the Bravo Dome field allows identification of a clear and coherent depletion of He to Ne and He to Ar in both datasets. This depletion can only be explained by partial degassing of small melt fractions of asthenospheric melts that have been emplaced into the SCLM. This is the first time that it has been possible to resolve and account for both the mantle He/Ne and He/Ar ratios within a SCLM source. The data additionally rule out the involvement of a plume component in the mantle source of the two gas fields and hence any plume influence on the Colorado Plateau Uplift event.

Effects of Partial Melting on Seismic Velocity and Attenuation: A New Insight from Experiments

The effects of partial melting on seismic velocity and attenuation have long been studied by focusing on the direct effects of melt, such as the poroelastic effect. The direct effects are generally very small for a very small melt fraction. Because geochemical studies have shown that the melt fraction during partial melting is very small (∼0.1%), it is difficult to explain upper-mantle low-velocity regions by the direct effects of melt. Recent experimental studies, by using a rock analog, have captured a significant enhancement of polycrystal anelasticity just before partial melting in the absence of melt. This newly recognized effect enables us to interpret seismological and geochemical observations consistently. The new anelasticity model significantly changes the interpretation of upper-mantle seismic structures. This review summarizes the recent progress in the understanding of polycrystal anelasticity, starting from a basic knowledge of linear anelasticity.

Asthenospheric kimberlites: Volatile contents and bulk compositions at 7 GPa

During ascent, kimberlites react with the lithospheric mantle, entrain and assimilate xenolithic material, loose volatiles and suffer from syn- and post-magmatic alteration. Consequently, kimberlite rocks deviate heavily from their primary melt. Experiments at 7 GPa, 1300–1480 °C, 10–30 wt% CO2 and 0.46 wt% H2O on a proposed primitive composition from the Jericho kimberlite show that saturation with a lherzolitic mineral assemblage occurs only at 1300–1350 °C for a carbonatitic melt with <8 wt% SiO2 and >35 wt% CO2. At asthenospheric temperatures of >1400 °C, where the Jericho melt stays kimberlitic, this composition saturates only in low-Ca pyroxene, garnet and partly olivine.We hence forced the primitive Jericho kimberlite into multiple saturation with a lherzolitic assemblage by adding a compound peridotite. Saturation in olivine, low- and high-Ca pyroxene and garnet was obtained at 1400–1650 °C (7 GPa), melts are kimberlitic with 18–29 wt% SiO2 + Al2O3, 22.1–24.6 wt% MgO, 15–27 wt% CO2 and 0.4–7.1 wt% H2O; with a trade-off of H2O vs. CO2 and temperature. Melts in equilibrium with high-Ca pyroxene with typical mantle compositions have ≥2.5 wt% Na2O, much higher than the commonly proposed 0.1–0.2 wt%.The experiments allow for a model of kimberlite origin in the convective upper mantle, which only requires mantle upwelling that causes melting at the depth where elemental carbon (in metal, diamond or carbide) converts to CO2 (at ∼250 km). If primary melts leading to kimberlites contain a few wt% H2O, then adiabatic temperatures of 1400–1500 °C would yield asthenospheric mantle melts that are kimberlitic (>18 wt% SiO2 + Al2O3) but not carbonatitic (<10 wt% SiO2 + Al2O3) in composition, carbonatites only forming 100–200 °C below the adiabat. These kimberlites represent small melt fractions concentrating CO2 and H2O and then acquire part of their chemical signature by assimilation/fractionation during ascent in the subcratonic lithosphere.

Polycrystal anelasticity at near‐solidus temperatures

AbstractElasticity, anelasticity, and viscosity of polycrystalline aggregates were measured at the near‐solidus temperatures ranging from below to above the solidus temperature (Tm). The result shows that the mechanical effects of the partial melting are twofold; changes just below the solidus temperature in the absence of melt and changes at the solidus temperature due to the onset of partial melting. As homologous temperature (T/Tm) increases from about 0.92 to 1, high‐frequency part of the attenuation spectrum significantly grows. Viscosity of the grain boundary diffusion creep is also reduced in this temperature range. These changes are caused by a solid‐state mechanism and have a large amplitude even for the samples which can generate very small amounts of melt at the solidus temperature. At the onset of melting, further increases in the elastic, anelastic, and viscous compliances occur. These changes are caused by the direct effects of the melt phase and are very small for the samples with very small melt fractions. Mechanical properties of a partially molten aggregate are determined by these twofold changes, and when melt fraction is small, the former changes are dominant. We performed a parameterization of the present experimental results and applied the obtained empirical formula to the seismic tomographic data in the upper mantle. The present model explains well the steep reduction of the seismic shear wave velocity in the oceanic lithosphere just below the solidus temperature.

Melting and dynamic metasomatism of mixed harzburgite + glimmerite mantle source: Implications for the genesis of orogenic potassic magmas

Tectonically young, orogenic settings are commonly the sites of post-collisional silica-rich ultrapotassic magmas with extreme K2O-contents of up to 9wt% and K2O/Na2O>2. Many experimental studies investigating the generation of these melts have concentrated on melting of homogenous phlogopite bearing peridotites, whereas geochemical signatures indicate the involvement of at least two types of source rocks: ultra-depleted and K and trace elements-enriched ones. We report the results of melting experiments at 1–2GPa of mixed glimmerite and harzburgite, in which these rock types make up two halves each capsule. Melting begins in the glimmerite, and its metasomatic effects on the harzburgite are apparent at 1100°C even before melt pools are visible. The first melts are Na-rich, seen in zoning of olivines and as growth of clinopyroxene in the harzburgite, but change at higher degrees of melting to produce a typical lamproite-like melt with K2O>10wt%. A major advantage of this study is the preservation of distinct melts in different parts of the capsule, which reflect a process of dynamic metasomatism: within the harzburgite matrix, the infiltrating melt derived from melting of the glimmerite changes consistently with the distance of travel through the harzburgite, enabling quantification of the metasomatic effects as an increase in SiO2 and K2O. This results principally from assimilation of orthopyroxene, which increases the Ol/Opx ratio of the residual harzburgite. The effects of quench olivine growth are recognizable and can be quantified due to a step-change in composition at the glimmerite/harzburgite border: the large total surface area of olivine and small melt fraction mean that the amount of quench olivine is high within the harzburgite, but negligible in the almost completely molten glimmerite. Melts of the glimmerite contain up to 8–10wt% K2O and 53wt% SiO2, which increase to 55–56wt% after interaction with the harzburgite. Mediterranean lamproites resemble melts of glimmerite, whereas melts that have interacted with harzburgite are more similar to less potassic, but more SiO2-rich shoshonites of the Mediterranean region.

Fractional crystallization of Si-undersaturated alkaline magmas leading to unmixing of carbonatites on Brava Island (Cape Verde) and a general model of carbonatite genesis in alkaline magma suites

The carbonatites of Brava Island, Cape Verde hot spot, allow to investigate whether they represent small mantle melt fractions or form through extreme fractionation and/or liquid immiscibility from CO2-bearing silicate magmas. The intrusive carbonatites on Brava Island are part of a strongly silica-undersaturated pyroxenite, ijolite, nephelinite, nepheline syenite, combeite–foiditite, carbonatite series. The major and trace element composition of this suite is reproduced by a model fractionating olivine, clinopyroxene, perovskite, biotite, apatite, titanite, sodalite and FeTi oxides, all present as phenocrysts in the rocks corresponding to their fractionation interval. Fractionation of ~90 wt% crystals reproduces the observed geochemical trend from the least evolved ultramafic dikes (bulk X Mg = 0.64) to syenitic compositions. The modelled fractional crystallization leads to alkali enrichment, driving the melt into the carbonatite–silicate miscibility gap. An initial CO2 content of 4000 ppm is sufficient to saturate in CO2 at the point where the rock record suggests continuing unmixing carbonatites from nephelinites to nepheline syenites after 61 wt% fractionation. Such immiscibility is also manifested in carbonatite and silicate domains on a hand-specimen scale. Furthermore, almost identical primary clinopyroxene, biotite and carbonate compositions from carbonatites and nephelinites to nepheline syenites substantiate their conjugate character and our unmixing model. The modelled carbonatite compositions correspond to the natural ones except for their much higher alkali contents. The alkali-poor character of the carbonatites on Brava and elsewhere is likely a consequence of the release of alkali-rich CO2 + H2O fluids during final crystallization, which cause fenitization in adjacent rocks. We propose a general model for carbonatite generation during alkaline magmatism, where the fractionation of heavily Si-undersaturated, alkaline parent melts results in alkali and CO2 enrichment in the evolving melt, ultimately leading to immiscibility between carbonatites and evolved Si-undersaturated alkaline melts. Early saturation in feldspathoids or feldspars would limit alkali enrichment preventing the formation of carbonatites. The complete and continuous fractionation line from almost primitive melts to syenitic compositions on Brava underlines the possibly important role of intrusives for hot spot volcanism.