The effect of sloped isotherms on melt migration in the shallow mantle: a physical and numerical model based on observations in the Oman ophiolite

Field observations in the Oman ophiolite and petrological data are used to constrain a model of melt segregation at the top of the mantle beneath an oceanic spreading centre. Foliations and lineations in outcrops of mantle-derived peridotites oriented at high angle relative to the crust–mantle boundary have been interpreted as the footprint of a former axial asthenospheric convective upwelling several kilometers in cross-section that reached Moho levels. Basaltic melts migrating through this upwelling reacted with their host harzburgites and suffered fractional crystallization. The mantle–crust transition zone at the top of the upwelling is characterized by an very thick (about 400 m) dunite layer whose detailed structure and composition point to the development by compaction of a former “mantle mush”. The more important structures (in terms of volume of crystallization products) found in the underlying harzburgites are dunitic–troctolitic horizons a few meters thick and of lateral extent reaching 1 km and more. They crystallized at high temperature (>1190 °C) from melts similar to mid-ocean ridge basalts (MORB). These are called “sills” because they are sub-parallel to the crust–mantle boundary, but they can present a moderate dip (15° to 20° at most) relative to this paleo-horizontal surface. These observations have motivated the modelling of melt segregation by compaction within the crystallization domain inside the top convective boundary layer of the mantle upwelling. Two original inputs to the modelling are considered here: (i) the slope of the iso-curves of melt concentration due to the progressive cooling of the mantle in the boundary layer away from the axis of the rising convective flow; (ii) the reduction in permeability caused by the crystallization of the inter-granular melt. Modelling shows that a unique condition is required to generate the troctolite sills and the thick dunite layer nested at the top of the Maqsad diapir: namely a dramatic drop of the interstitial melt concentration at the top of the mantle. Besides, the model developed here allows to scale the time, volume and velocity of the melt segregation.     

Deformation, static recrystallization, and reactive melt transport in shallow subcontinental mantle xenoliths (Tok Cenozoic volcanic field, SE Siberia)

Partial melting and reactive melt transport may change the composition, microstructures, and physical properties of mantle rocks. Here we explore the relations between deformation and reactive melt transport through detailed microstructural analysis and crystallographic orientation measurements in spinel peridotite xenoliths that sample the shallow lithospheric mantle beneath the southeastern rim of the Siberian craton. These xenoliths have coarse-grained, annealed microstructures and show petrographic and chemical evidence for variable degrees of reaction with silicate melts and fluids, notably Fe-enrichment and crystallization of metasomatic clinopyroxene (cpx). Olivine crystal preferred orientations (CPO) range from strong to weak. [010]-fiber patterns, characterized by a point concentration of [010] normal to the foliation and by dispersion of [100] in the foliation plane with a weak maximum parallel to the lineation, predominate relative to the [100]-fiber patterns usually observed in lithospheric mantle xenoliths and peridotite massifs. Variations in olivine CPO patterns or intensity are not correlated with modal and chemical compositions. This, together with the analysis of microstructures, suggests that reactive melt percolation postdated both deformation and static recrystallization. Preferential crystallization of metasomatic cpx along (010) olivine grain boundaries points to an influence of the preexisting deformation fabrics on melt transport, with higher permeability along the foliation. Similarity between orthopyroxene (opx) and cpx CPO suggests that cpx orientations may be inherited from those of opx during melt-rock reaction. As observed in previous studies, reactive melt transport does not weaken olivine CPO and seismic anisotropy in the upper mantle, except in melt accumulation domains. In contrast, recovery and selective grain growth during static recrystallization may lead to development of [010]-fiber olivine CPO and, if foliations are horizontal, result in apparent isotropy for vertically propagating SKS waves, but strong anisotropy for horizontally propagating surface waves.     

The influence of silicate melt composition on distribution of siderophile elements among metal and silicate liquids

Liquid metal-liquid silicate partitioning of Fe, Ni, Co, P, Ge, W and Mo among a carbon-saturated metal and a variety of silicate melts (magnesian-tholeiitic-siliceous-aluminous-aluminosiliceous basalts) depends modestly to strongly upon silicate melt structure and composition. Low valency siderophile elements, Fe, Ni and Co, show a modest influence of silicate melt composition on partitioning. Germanium shows a moderate but consistent preference for the depolymerized magnesian melt. High valency siderophile elements, P, Mo, and W, show more than an order of magnitude decrease in metal-silicate partition coefficients as the silicate melt becomes more depolymerized. Detailed inspection of our and other published W data shows that polymerization state, temperature and pressure are more important controls on W partitioning than oxidation state. For this to be true for a high and variable valence element implies a secondary role in general for oxidation state, even though some role must be present. Equilibrium core segregation through a magma ocean of ‘ultrabasic’ composition can provide a resolution to the ‘excess’ abundances of Ge, P, W and Mo in the mantle, but the mantle composition alone cannot explain the excess abundances of nickel and cobalt in chondritic proportions.     

Silicate and carbonate melt inclusions associated with diamonds in deeply subducted carbonate rocks

Deeply subducted carbonate rocks from the Kokchetav massif (Northern Kazakhstan) recrystallised within the diamond stability field (P = 4.5–6.0 GPa; T ≈ 1000 °C) and preserve evidence for ultra high-pressure carbonate and silicate melts. The carbonate rocks consist of garnet and K-bearing clinopyroxene embedded in a dolomite or magnesian calcite matrix. Polycrystalline magnesian calcite and polyphase carbonate–silicate inclusions occurring in garnet and clinopyroxene show textural features of former melt inclusions. The trace element composition of such carbonate inclusions is enriched in Ba and light rare earth elements and depleted in heavy rare earth elements with respect to the matrix carbonates providing further evidence that the inclusions represent trapped carbonate melt. Polyphase inclusions in garnet and clinopyroxene within a magnesian calcite marble, consisting mainly of a tight intergrowth of biotite + K-feldspar and biotite + zoisite + titanite, are interpreted to represent two different types of K-rich silicate melts. Both melt types show high contents of large ion lithophile elements but contrasting contents of rare earth elements. The Ca-rich inclusions display high REE contents similar to the carbonate inclusions and show a general trace element characteristic compatible with a hydrous granitic origin. Low SiO₂ content in the silicate melts indicates that they represent residual melts after extensive interaction with carbonates. These observations suggest that hydrous granitic melts derived from the adjacent metapelites reacted with dolomite at ultra high-pressure conditions to form garnet, clinopyroxene – a hydrous carbonate melt – and residual silicate melts. Silicate and carbonate melt inclusions contain diamond, providing evidence that such an interaction promotes diamond growth. The finding of carbonate melts in deeply subducted crust might have important consequences for recycling of trace elements and especially C from the slab to the mantle wedge.     

Mineralogy of Antarctic ureilites and a working hypothesis for their origin and evolution

Mineralogical and crystallographic data and electron microprobe analyses of silicate minerals for Antarctic ureilites Y74123, Y74130, Y74659, Y790981, ALHA77257, ALHA78019 and ALH82106 extend the chemical and textural variations among ureilites beyond those known for non-Antarctic ureilites. Low-Ca pyroxenes, possibly orthopyroxene, poikilitically enclosed in olivine and coexisting with Mg-rich pigeonites, have been identified in Y74659 and ALHA77257. Coexisting low-Ca pyroxene, pigeonite, and augite, together with the very Fe-rich olivine in ureilites were observed in Y74130. Almost Fe-free pigeonite and diopside coexist in ALH82106. No systematic differentiation trend has been reconstructed from these data, but the pyroxene pairs yield formation temperatures above 1200°C. A weak anticorrelation between MnO and FeO delineated by these pyroxenes suggests that a reduction reaction is responsible for the chemical variations. A working hypothesis for ureilite formation involves high-temperature episoded and a planetesimal-scale collision, which might have facilitated compaction and recrystallization of mafic silicates. During the episodes, FeNiS eutectic melt and a partial melt rich in Ca, Al and Fe might have been lost from a carbonaceous-chondrite-like source material.     

Nickel partitioning between olivine and silicate melt

Partitioning of Ni between olivine and silicate melt has been determined for compositions in the system Fo-Ab-An (1 atm) for temperatures ranging from 1250°C to 1450°C. Nickel concentrations were determined by electron microprobe; concentration levels in the liquids ranged from 0.1% to 0.5%. Platinum capsules or Pt wire loops were used as containers. Equilibrium was evaluated from kinetic considerations and by variation of run parameters; it was documented in one case by a bracketed reversal. No evidence was found for a dependence of the partition coefficient D (Ni in olivine/Ni in liquid) on Ni concentration. D is strongly dependent on melt composition, varying linearly with (1/MgO) at constant temperature. The intrinsic temperature dependence of D is small; the apparent temperature dependence reported in previous studies is largely related to the variation of melt composition with temperature. Our D values determined in the simple system Fo-An-Ab agree well with those reported by Leeman for natural (Fe-bearing) basalt systems. Overall variation of D in our system (and in natural basalts) can be expressed by the regression: D = (124/MgO) ? 0.9Our data are used to evaluate published Ni-MgO relationships in natural basalt series from Kilauea, Crozet, Cape Verde and Baffin Bay. A combination of olivine accumulation and fractional crystallization processes are sufficient to model these series. Using our data, unique “parental” liquids can be specified for each of these series; the MgO content of these liquids varies from 6% to 13%. Basalts with MgO contents greater than these “parental” liquids must be accumulative. The linear Ni-MgO trends, high absolute Ni concentrations, and large spread of Ni contents for the high-MgO basalts argue convincingly against their being “primary” liquids. Models such as those of O'Hara [6,13] and Clarke [24], based on the assertion of primary high-MgO liquids, must therefore be re-evaluated.Because of the high Si/O ratio and low MgO content of island arc andesites, the Ni partition coefficient D may be quite high. Therefore, the relatively low Ni content of such andesites may not be an argument against their derivation as direct partial melts of the mantle.     

Petrological feature of spinel lherzolite xenolith from Oki-Dogo Island: An implication for variety of the upper mantle peridotite beneath southwestern Japan

Abstract   Spinel lherzolite is a minor component of the deep-seated xenolith suite in the Oki-Dogo alkaline basalts, whereas other types of ultramafic (e.g. pyroxenite and dunite) and mafic (e.g. granulite and gabbro) xenoliths are abundant. All spinel lherzolite xenoliths have spinel with a low Cr number (Cr#; < 0.26). They are anhydrous and are free of modal metasomatism. Their mineral assemblages and microtextures, combined with the high NiO content in olivine, suggest that they are of residual origin. But the Mg numbers of silicate minerals are lower (e.g. down to Fo₈₆) in some spinel lherzolites than in typical upper mantle residual peridotites. The clinopyroxene in the spinel lherzolite shows U-shaped chondrite-normalized rare-earth element (REE) patterns. The abundance of Fe-rich ultramafic and mafic cumulate xenoliths in Oki-Dogo alkali basalts suggests that the later formation of those Fe-rich cumulates from alkaline magma was the cause of Fe- and light REE (LREE)-enrichment in residual peridotite. The similar REE patterns are observed in spinel peridotite xenoliths from Kurose and also in those from the South-west Japan arc, which are non-metasomatized in terms of major-element chemistry (e.g. Fo > 89), and are rarely associated with Fe-rich cumulus mafic and ultramafic xenoliths. This indicates that the LREE-enrichment in mantle rocks has been more prominent and prevalent than Fe and other major-element enrichment during the metasomatism.     

Si- and alkali-rich melt inclusions in minerals of mantle peridotites from eastern China: Implication for lithospheric evolution

CO2 fluid inclusions in mantle minerals are an im-portant source for us to get the information of mantle fluids. Fluid inclusions are mainly composed of CO2, with minor CO, H2O, CH4, N2, H2S, SO2, F, etc., which were demonstrated by lots of Raman spec-trometer analyses in recent years. In contrast, there are very few researches on CO2-bearing melt inclusions since it is more difficult to do so. The available studies have found that the primary CO2-bearing melt inclu-sions are basaltic …     

The ‘small-world’ nature of fracture/conduit networks: Possible implications for disequilibrium transport of magmas beneath mid-ocean ridges

In this contribution we show that natural fracture/conduit networks can be studied by using a new method based on Graph Theory. A number of natural networks at different length scales (from the meter to the millimeter) are analysed and results show that they have typical attributes of ‘small-world’ networks, a special class of networks characterized by high global and local transport efficiency. To our knowledge, this topological feature of natural fracture networks is recognized here for the first time. By starting from results on natural fracture/conduit networks, the possible implications are discussed by focusing on disequilibrium transport of magmas in the upper mantle beneath mid-ocean ridges. Results indicate that the ‘small-world’ topology of natural fracture/conduit networks is an important characteristic to ensure disequilibrium delivery of melts through the upper mantle, thus offering a good explanation of geochemical features of magmas. The remarkable point here is that the modelling of melt migration has been constrained by using real fracture network systems. The results presented in this work may contribute to a better understanding of melt migration in fracture network systems and of the way geochemical features of magmas may be influenced by their transport history.     

A conceptual model for kimberlite emplacement by solitary interfacial mega-waves on the core mantle boundary

If convection in the Earth's liquid outer core is disrupted, degrades to turbulence and begins to behave in a chaotic manner, it will destabilize the Earth's magnetic field and provide the seeds for kimberlite melts via turbulent jets of silicate rich core material which invade the lower mantle. These (proto-) melts may then be captured by extreme amplitude solitary nonlinear waves generated through interaction of the outer core surface with the base of the mantle. A pressure differential behind the wave front then provides a mechanism for the captured melt to ascend to the upper mantle and crust so quickly that emplacement may indirectly promote a type of impact fracture cone within the relatively brittle crust. These waves are very rare but of finite probability. The assumption of turbulence transmission between layers is justified using a simple three-layer liquid model. The core derived melts eventually become frozen in place as localised topographic highs in the Mohorovicic discontinuity (Moho), or as deep rooted intrusive events. The intrusion's final composition is a function of melt contamination by two separate sources: the core contaminated mantle base and subducted Archean crust. The mega-wave hypothesis offers a plausible vehicle for early stage emplacement of kimberlite pipes and explains the age association of diamondiferous kimberlites with magnetic reversals and tectonic plate rearrangements.     

Nd-isotope systematics of ∼2.7 Ga adakites, magnesian andesites, and arc basalts, Superior Province: evidence for shallow crustal recycling at Archean subduction zones

An association of adakite, magnesian andesite (MA), and Nb-enriched basalt (NEB) volcanic flows, which erupted within ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts, has recently been documented in ∼2.7 Ga Wawa greenstone belts. Large, positive initial ?_Nd values (+1.95 to +2.45) of the adakites signify that their basaltic precursors, with a short crustal residence, were derived from a long-term depleted mantle source. It is likely that the adakites represent the melts of subducted late Archean oceanic crust. Initial ?_Nd values in the MA (+0.14 to +1.68), Nb-enriched basalts and andesites (NEBA) (+1.11 to +2.05), and ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts and andesites (+1.44 to +2.44) overlap with, but extend to lower values than, the adakites. Large, tightly clustered ?_Nd values of the adakites, together with Th/Ce and Ce/Yb systematics of the arc basalts that rule out sediment melting, place the enriched source in the sub-arc mantle. Accordingly, isotopic data for the MA, NEBA, and ‘normal’ arc basalts can be explained by melting of an isotopically heterogeneous sub-arc mantle that had been variably enriched by recycling of continental material into the shallow mantle in late Archean subduction zones up to 200 Ma prior to the 2.7 Ga arc. If the late Archean Wawa adakites, MA, and basalts were generated by similar geodynamic processes as their counterparts in Cenozoic arcs, involving subduction of young and/or hot ocean lithosphere, then it is likely that late Archean oceanic crust, and arc crust, were also created and destroyed by modern plate tectonic-like geodynamic processes. This study suggests that crustal recycling through subduction zone processes played an important role for the generation of heterogeneity in the Archean upper mantle. In addition, the results of this study indicate that the Nd-isotope compositions of Archean arc- and plume-derived volcanic rocks are not very distinct, whereas Phanerozoic plumes and intra-oceanic arcs tend to have different Nd-isotopic compositions.     

Melt–wall rock interaction in the mantle shown by silicate melt inclusions in peridotite xenoliths from the central Pannonian Basin (western Hungary)

In this paper we present a detailed textural and geochemical study of two equigranular textured amphibole-bearing spinel lherzolite xenoliths from Szigliget, Bakony–Balaton Highland Volcanic Field (BBHVF, western Hungary) containing abundant primary silicate melt inclusions (SMIs) in clinopyroxene rims and secondary SMIs in orthopyroxene (and rarely spinel) along healed fractures. The SMIs are dominantly composed of silicate glass and CO₂-rich bubbles. Clinopyroxene and orthopyroxene are zoned in both studied xenoliths, especially with respect to Fe, Mg, Na, and Al contents. Cores of clinopyroxenes in both xenoliths show trace element distribution close to primitive mantle. Rims of clinopyroxenes are enriched in Th, U, light rare earth elements (LREEs) and medium REEs (MREEs). Amphiboles in the Szg08 xenolith exhibit elevated Rb, Ba, Nb, Ta, LREE, and MREE contents. The composition of silicate glass in the SMIs covers a wide range from the basaltic trachyandesite and andesite to phonolitic compositions. The glasses are particularly rich in P₂O₅. Both primary and secondary SMIs are strongly enriched in incompatible trace elements (mostly U, Th, La, Zr) and display a slight negative Hf anomaly. The development of zoned pyroxenes, as well as the entrapment of primary SMIs in the clinopyroxene rims, happened after partial melting and subsequent crystallization of clinopyroxenes, most probably due to an interaction between hot volatile-bearing evolved melt and mantle wall-rocks. This silicate melt filled microfractures in orthopyroxenes (and rarely spinels) resulting in secondary SMIs.     

Infrared imagery of ‘breaking’ internal waves

Airborne infrared imagery is shown to provide preliminary evidence of surface thermal expressions associated with internal waves that become unstable and break over the continental shelf. These expressions include a narrow wave front that is warmer than the ambient; a wide, spatially intermittent ‘wake’ that is colder than the ambient; and ∼O (10 m) diameter surface-renewal ‘boils’ that populate the wake. These thermal signatures might be useful in assessing the spatial distribution and structure of breaking internal waves.     

Static compression of hydrous silicate melt and the effect of water on planetary differentiation

High pressure experiments using the sink/float method have bracketed the density of hydrous iron-rich ultrabasic silicate melt from 1.35 to 10.0 GPa at temperatures from 1400 to 1860 °C. The silicate melt composition was a 50–50 mixture of natural komatiite and synthetic fayalite. Water was added in the form of brucite Mg(OH)₂ and was present in the experimental run products at 2 wt.% and 5 wt.% levels as confirmed by microprobe analyses of total oxygen. Buoyancy marker spheres were olivines and garnets of known composition and density. The density of the silicate melt with 5 wt.% water at 2 GPa and 1500 °C is 0.192 g cm^? 3 less than the anhydrous form of this melt at the same P and T. This density difference gives a partial molar volume of water in silicate melt of ～ 7 cm³ mol^? 1, which is similar to previous studies at high pressure. The komatiite–fayalite liquids with 0 and 2 wt.% H₂O, have extrapolated density crossovers with equilibrium liquidus olivine at 8 and 9 GPa respectively, but there is no crossover for the liquid with 5 wt.% H₂O. These results are consistent with the hypothesis that dense hydrous melts could be gravitationally stable atop the 410 km discontinuity in the Earth. The results also support the notion that equilibrium liquidus olivine could float in an FeO-rich hydrous martian magma ocean. Extrapolation of the data suggests that FeO-rich hydrous melt could be negatively buoyant in the Earth's D″-region or atop the core–mantle-boundary (CMB), although experiments at higher pressure are needed to confirm this prediction.     

Magma superflows in the Bering Sea. 2. The mechanism of origin and movement and its evolution in time and space

It is shown that the source for numerous magma superflows (MSFs) with complex structures, whose traces have been identified in the Bering Sea (Melekestsev and Slezin, 2017), was regional mantle plume-like formations (or plumes). The magmatic material propagated from these for some hundreds of kilometers or farther along interfaces at different depths between crustal layers during n × 10⁵ to n × 10⁶ years. The long-continued generation and multi-portion structure of the MSFs is explained by slow pulsations in the rising jet of primary melt in the form of a beaded (“peristalsis-like”) structure of successive bulges (which we shall refer to as “magmons”), which float up as “asthenoliths” as the through flow in the channel becomes lower. The most extensive and longest MSFs were formed at the end of the Eocene through Oligocene, while the youngest and shortest flow was generated in the Pleistocene.     

Subophitic lithologies in KREEP-rich poikilitic impact melt rocks from Cayley Plains,Apollo 16 — remnants of a volcanic highland crust?

KREEP-rich poikilitic impact melt rocks 65777,11, 65015,88, and 62235,66 are the only mafic impact melt rocks from Cayley Plains stations, Apollo 16, from which areas of subophitic texture can be reported.The bulk chemistry of these unique subophitic areas and the surrounding poikilitic matrices, as well as mineral compositions (olivine, plagioclase, pyroxene, Fe-Ni metal) were determined by electron microprobe analysis. All subophitic areas could be undoubtedly identified as impact melt rockclasts. Inclusion 65777,11 II is of uniquely KREEP-rich composition, 62235,66 II can be classified as anorthositic. Therefore our attempt to identify pristine volcanic basement rocks of the Cayley regions among these inclusions of basaltic texture failed.However, the absence of pristine volcanic target rock fragments and the existence of KREEP-rich and anorthositic impact melt clasts in KREEP-rich impact melt rocks from Cayley Plains favors the theory that the Cayley Plains formation is part of the ejecta blanket from a large basin-type impact crater (Imbrium?), which is underlain by anorthositic material (Nectaris ejecta?), and has been reworked by local impacts in post-Imbrian times.     

Experimental studies of thermal grooving in the olivine and albite melt system

Thermal grooving of low angle tilt boundary of San Carlos olivine in the albite melt were experimentally investigated at 1200–1300°C in mixed CO₂ and H₂ gases for 1–20 h. The depth, d, of the thermal groove on (010) of olivine along the (100) sub-boundaries is in the function of time and temperatures as follows; d⁴ = k_o · t · exp(− 190 000/RT), in which R is the gas constant, and k_o is the material constant.The melt shape changes due to the thermal grooving driven by surface tension and deformation of the upper mantle. Compared with the time scales of these two counteracting mechanisms, it is inferred that the melt shape is unstable in the high temperature and low stress conditions, and that the melt shape takes a stable form during progressive deformation in the low temperature and high stress conditions.     

Reaction experiments between tonalitic melt and mantle olivine and their implications for genesis of high-Mg andesites within cratons

High-Mg (Mg#>45) andesites (HMA) within cratons attract great attention from geologists. Their origin remains strongly debated. In order to examine and provide direct evidence for previous assumptions about HMA’s genesis inferred from petrological and geochemical investigations, we performed reaction experiments between tonalitic melt and mantle olivine on a six-anvil apparatus at high-temperature of 1250–1400°C and high-pressure of 2.0–5.0 GPa. Our experiments in this work simulated the interaction between the tonalitic melt derived from partial melting of eclogitized lower crust foundering into the Earth’s mantle and mantle peridotite. The experimental results show that the reacted melts have very similar variations in chemical compositions to the HMAs within the North China Craton. Therefore, this interaction is probably an important process to generate the HMAs within crations.     

Major- and trace-element systematics and isotope geochemistry of Cenozoic mafic volcanic rocks from the Vogelsberg (central Germany): Constraints on the origin of continental alkaline and tholeiitic basalts and their mantle sources

The Cenozoic basaltic province of the Vogelsberg area (central Germany) is mainly composed of intercalated olivine to quartz tholeiites and near-primary nephelinites to basanites. The inferred mantle source for the alkaline and tholeiitic rocks is asthenospheric metasomatized garnet peridotite containing some amphibole as the main hydrous phase. Trace element modelling indicates 2 to 3% partial melting for the alkaline rocks and 5 to 7% partial melting for the olivine tholeiites. Incompatible trace element abundances and ratios as well as Nd and Sr radiogenic isotope compositions lie between plume compositions and enriched mantle compositions and are similar to those measured in Ocean Island Basalts (OIB) and the Central European Volcanic Province elsewhere. The mafic olivine tholeiites have similar Ba/Nb, Ba/La and Nd–Sr isotope ratios to the alkaline rocks indicating derivation of both magma types from chemically comparable mantle sources. However, Zr/Nb ratios are slightly higher in olivine tholeiites than in basanites reflecting some fractionation of Zr relative to Nb during partial melting. Quartz tholeiites have higher Ba/Nb, Zr/Nb, La/Nb, but lower Ce/Pb ratios and lower Nd isotope compositions than the alkaline rocks which can be explained by interaction of the basaltic melt with lower (granulite facies) crustal material or partial melts thereof during stagnation within the lower crust. It appears most likely that upwelling of hot, asthenospheric material results in the generation of primitive alkaline rocks at the base of the lithosphere at depths of 75–90 km. Lithospheric extension together with minor plume activity and probably lower lithosphere erosion induced melting of shallower heterogenous upper mantle generating a spectrum of olivine tholeiitic melts. These olivine tholeiitic rocks evolved via crystal fractionation and probably limited contamination to quartz tholeiites.