Non-modal melting in an upwelling mantle column: Steady-state models with applications to REE depletion in abyssal peridotites and the dynamics of melt migration in the mantle

Simple models for trace element fractionation during concurrent melting and melt migration in an upwelling steady-state mantle were developed. Based on petrologic considerations, we divided the mantle column into two regions: a single-lithology lower region that consists of partially molten garnet and spinel lherzolites and a double-lithology upper region where high-porosity dunite channels or melt-filled fractures are embedded in a porous lherzolite/harzburgite matrix. Analytical solutions for the case of a constant and uniform relative melting suction rate and a linearly variable relative melt suction rate were obtained. Key parameters and the first order characteristics of melting and melt migration in a 1-D steady-state mantle column were examined through forward calculations and Monte Carlo simulations. Melting in the upwelling single-lithology column is equivalent to non-modal batch melting, whereas melting and melt migration in the double-lithology region can be viewed as a nonlinear combination of batch melting and fractional melting, depending on the amount of melt extracted to the channel. The degree of melting (F), the degree of melting at the depth of melt-channel initiation (F_d) and the relative rate of melt suction (R) are important in controlling the extent of depletion of the incompatible trace element in the matrix. Spatially variable R affects the abundance of an incompatible trace element in the melt and residual solid the most in near fractional melting. There is a strong nonlinear trade off among the three parameters. Given F_d, it is possible to constrain F and R from incompatible trace element abundances in residual peridotite.To explore the dynamics of melt migration in the mantle, we used the two melting models developed in this study and published REE and Y abundances in diopside in abyssal peridotites from the Central Indian Ridge to infer their melting and melt migration history. Overall, the degrees of melting inferred from the trace element data are not sensitive to the value of F_d used in the inversion and ranges from 10% to 15%. The relative rate of melt suction depends slightly on the choice of F_d and ranges from 0.85 to 1.0 for F_d = 0.05 and 0.75 to 0.97 for F_d = 0. Further, the estimated R is inversely correlated with F, a robust feature independent of the choice of F_d. The upward decrease of R in an upwelling mantle column can be understood in terms of melt focusing in the lower part of the double-lithology region. And finally, given F and R, we found that the permeability and porosity of the lherzolite/harzburgite matrix also increase as a function of F in the melting column, with melt fractions ranging from 0.2% to 0.7% for a grain size of 5 mm.     

Crustal accretion and dynamic feedback on mantle melting of a ridge centred plume: The Iceland case

A major consequence of the interaction of a plume with an oceanic ridge is the enhanced melt production and associated crust generation. In the case of Iceland crustal thickness as large as 20 to 40 km has been reported. Crustal seismic velocities are high, and have to be explained by thermal or chemical effects. In the first part of the paper we address the question whether extraction of melt out of the plume beneath a slowly spreading ridge and deposition of extracted basalt volumes at the surface produces a dynamic feedback mechanism on mantle melting. To study this question we solve the convection equations for a ridge centred plume with non-Newtonian rheology including melting, melt extraction associated with deposition of cold crust at the surface of the model, and using a simplified approach for compaction. The assumption of cold crust is justified if the thickness of each deposited basaltic layer is less than roughly 1 km. Depending on the buoyancy flux of the plume, crustal thicknesses between 10 and 40 km are modelled, showing characteristic dipping structures resembling the rift-ward dipping basaltic layers of East- and Western Iceland. Comparing the resulting crustal thickness and magma generation rate with models in which the dynamic effect of crust deposition has been suppressed indicates, that melt generation beneath a slowly spreading ridge is considerably damped by the dynamic feedback mechanism if the plume buoyancy flux exceeds 400 to 600 kg/s. Based on the observed crustal thickness of Iceland our models predict a plume buoyancy flux of 1140 kg/s.In the second part we study the accretion of the Icelandic crust by a thermo-mechanical model in more detail based on the Navier–Stokes-, the heat transport and the mass conservation equations including volumetric sources. Hot (1200 °C) molten crustal material is injected into the newly forming crust with a constant rate at different crustal source regions: a) deep, widespread emplacement of dykes and sills including crustal underplating, b) magma chambers at shallow to mid-crustal level, and c) surface extrusions and intrusions in fissure swarms at shallow depth connected to volcanic centres. We identify the material from the different source regions by a marker approach. Varying the relative dominance of these source regions, characteristic crustal structures evolve, showing shallow dipping upper crustal layers with dip angles between 10 and 15°. The thermal structure of the crust varies between cold crust (shallow-source region dominating) and hot crust (deep-source region dominating). We use observations of maximum depth of seismicity to constrain the depth of the 650 °C isotherm and seismological inferences on the lower crust to constrain temperatures in that region. The best agreement with our models is achieved for crust formation dominated by deep dykes and underplating with a considerable influence of magma chamber accretion.     

Possible implications of modal mineralogy for melting in mantle lherzolites

The melting reaction at the solidus of mantle peridotite is commonly peritectic in nature, with liquid and one or more solid phases produced upon melting. In some situations, one of the phases participating on the reactant side of the reaction is present in low abundance. This article explores the possible effects of the low abundance of a reactant phase on the melting behavior of mantle peridotite.For example, spinel lherzolite begins to melt via the peritectic reaction, clinopyroxene + orthopyroxene + spinel = olivine + liquid in the ∼1- to 2-GPa pressure range. In natural spinel lherzolites, spinel is a modally minor mineral and may be infrequently in contact with both clinopyroxene and orthopyroxene. If these mutual contacts are insufficient to generate an interconnected melt, then significant melting may not occur until a combination of minerals that are modally abundant and in contact begin to melt. This scenario could have implications for the physical process of melting and for the timing of formation of an interconnected melt network and separation of the melt from the residue.To begin to investigate this possibility, the spatial relationships between the constituent minerals in two fertile spinel lherzolites were determined by elemental mapping with the electron microprobe. Olivine, orthopyroxene, and clinopyroxene are of similar size, whereas the spinel was smaller and interstitial. Spinel and clinopyroxene are frequently in contact, but mutual contacts of spinel, clinopyroxene, and orthopyroxene are rare. Because of the changes in modal mineralogy anticipated for these lherzolites with increasing temperature, these mutual contacts will be even less common at the solidus. Therefore, an interconnected, potentially extractable, melt may not occur by the solidus spinel + orthopyroxene + clinopyroxene melting reaction.     

The evolution of continental roots in numerical thermo-chemical mantle convection models including differentiation by partial melting

Incorporating upper mantle differentiation through decompression melting in a numerical mantle convection model, we demonstrate that a compositionally distinct root consisting of depleted peridotite can grow and remain stable during a long period of secular cooling. Our modeling results show that in a hot convecting mantle partial melting will produce a compositional layering in a relatively short time of about 50 Ma. Due to secular cooling mantle differentiation finally stops before 1 Ga. The resulting continental root remains stable on a billion year time scale due to the combined effects of its intrinsically lower density and temperature-dependent rheology. Two different parameterizations of the melting phase-diagram are used in the models. The results indicate that during the Archaean melting occurred on a significant scale in the deep regions of the upper mantle, at pressures in excess of 15 GPa. The compositional depths of continental roots extend to 400 km depending on the potential temperature and the type of phase-diagram parameterization used in the model. The results reveal a strong correlation between lateral variations of temperature and the thickness of the continental root. This shows that cold regions in cratons are stabilized by a thick depleted root.     

Effects of mantle upwelling in a compressional setting: the Atlas Mountains of Morocco

We discuss the implications of a lithospheric model of the Moroccan Atlas Mountains based on topography, heat flow, gravity and geoid anomalies, taking into account the regional geology. The NW African cratonic lithosphere, some 160–180 km thick, thins to c. 80 km beneath the Atlas fold-thrust belts, in contrast with the shortening regime prevailing there since the early Cenozoic. This fact explains several geological and geophysical features as high topography with modest tectonic shortening, the occurrence of alkaline magmatism contemporaneous to compression, the absence of large crustal roots to support elevation, the scarce development of foreland basins, and a marked geoid high. The modelled lithosphere thinning is related to a thermal upwelling constrained between the Iberia–Africa convergent plate boundary and the Saharan craton.     

Anhydrous partial melting of an iron-rich mantle II: primary melt compositions at 15 kbar

Primary melt and coexisting mineral compositions, at increasing degrees of partial melting at 15 kbar, were determined for an iron-rich martian mantle composition, DW. The composition of primary melts near the solidus was determined with basalt-peridotite sandwich experiments. In order to evaluate the approach of the liquids to equilibrium with a DW mantle assemblage, experiments were also performed to establish the liquidus mineralogy of the primary melts. Primary melt compositions produced from an iron-rich mantle are more picritic than those produced from an iron-poor mantle. By increasing the iron content of a model mantle composition (decreasing the mg#, where mg# = atomic [Mg/(Mg+Fe²⁺)^*100]), picritic and komatiitic magmas result at lower percentages of melting and at temperatures closer to the solidus than in an iron-poor mantle. Terrestrial iron-rich primitive volcanics may be the partial melting products of iron-rich, mg# 80, source regions. The DW partial melting results support the conclusion of previous authors that the parent magmas of the SNC (shergottites, nakhlites, chassignites) meteorites were derived from a source region that had been previously depleted in an aluminous phase.     

Early Cretaceous mantle upwelling and melting of juvenile lower crust in the Middle-Lower Yangtze River Metallogenic Belt: Example from Tongshankou Cu-(MoW) ore deposit

The linkage between intracontinental extension and the Early Cretaceous Cu-Mo-W polymetallic metallogenesis in the Middle-Lower Yangtze River Belt (MLYRB) has long been a subject of controversy due to the lack of convincing petrogenetic evidence to identify the nature of magmatic sources and their geological histories during extensional mantle upwelling. Here we present new zircon UPb ages, isotopic and geochemical data for granodiorites, quartz diorites and mafic microgranular enclaves (MMEs) in the Tongshankou area. Comparing the MMEs with their host porphyries, the different ratios of incompatible elements and the similar formation ages, coupled with quenched margins and the xenocrysts in the MMEs, indicate that the MMEs was most likely formed by mixing between mafic magma and their host felsic magma. The MMEs share similar geochemical and isotopic characteristics with the Cretaceous mafic rocks from MLYRB, indicating that MMEs were mostly derived from an enriched lithospheric mantle source without adakitic characteristics. Mixing of a crustal melt derived by melting of an amphibolite bearing juvenile lower crust with a mantle melt derived from melting of enriched lithospheric mantle can account for the generation of the Tongshankou prophyries. The melting of juvenile mafic lower crust and enriched lithospheric mantle is suggested to be caused by upwelling of asthenospheric mantle and the reactivity of trans-lithospheric faults in the intracontinental extensional environment. Our results therefore highlight that juvenile mafic lower crust beneath the Yangtze plate is one of the likely source for ore-forming magmatic rocks in the Early Cretaceous.     

Anhydrous partial melting of an iron-rich mantle I: subsolidus phase assemblages and partial melting phase relations at 10 to 30 kbar

Anhydrous partial melting experiments, at 10 to 30 kbar from solidus to near liquidus temperature, have been performed on an iron-rich martian mantle composition, DW. The DW subsolidus assemblage from 5 kbar to at least 24 kbar is a spinel lherzolite. At 25 kbar garnet is stable at the solidus along with spinel. The clinopyroxene stable on the DW solidus at and above 10 kbar is a pigeonitic clinopyroxene. Pigeonitic clinopyroxene is the first phase to melt out of the spinel lherzolite assemblage at less than 20°C above the solidus. Spinel melts out of the assemblage about 50°C above the solidus followed by a 150° to 200°C temperature interval where melts are in equilibrium with orthopyroxene and olivine. The temperature interval over which pigeonitic clinopyroxene melts out of an iron-rich spinel lherzolite assemblage is smaller than the temperature interval over which augite melts out of an iron-poor spinel lherzolite assemblage. The dominant solidus assemblage in the source regions of the Tharsis plateau, and for a large percentage of the martian mantle, is a spinel lherzolite.     

Electrical conductivity of hydrous basaltic melts: implications for partial melting in the upper mantle

The Earth’s uppermost asthenosphere is generally associated with low seismic wave velocity and high electrical conductivity. The electrical conductivity anomalies observed from magnetotelluric studies have been attributed to the hydration of mantle minerals, traces of carbonatite melt, or silicate melts. We report the electrical conductivity of both H₂O-bearing (0–6 wt% H₂O) and CO₂-bearing (0.5 wt% CO₂) basaltic melts at 2 GPa and 1,473–1,923 K measured using impedance spectroscopy in a piston-cylinder apparatus. CO₂ hardly affects conductivity at such a concentration level. The effect of water on the conductivity of basaltic melt is markedly larger than inferred from previous measurements on silicate melts of different composition. The conductivity of basaltic melts with more than 6 wt% of water approaches the values for carbonatites. Our data are reproduced within a factor of 1.1 by the equation log σ = 2.172 − (860.82 − 204.46 w ^0.5)/(T − 1146.8), where σ is the electrical conductivity in S/m, T is the temperature in K, and w is the H₂O content in wt%. We show that in a mantle with 125 ppm water and for a bulk water partition coefficient of 0.006 between minerals and melt, 2 vol% of melt will account for the observed electrical conductivity in the seismic low-velocity zone. However, for plausible higher water contents, stronger water partitioning into the melt or melt segregation in tube-like structures, even less than 1 vol% of hydrous melt, may be sufficient to produce the observed conductivity. We also show that ~1 vol% of hydrous melts are likely to be stable in the low-velocity zone, if the uncertainties in mantle water contents, in water partition coefficients, and in the effect of water on the melting point of peridotite are properly considered.     

峨眉山大火成岩省：地幔柱活动的证据及其熔融条件

对苦橄岩中橄榄石斑晶及其中熔体包裹体的电子探针分析表明，峨眉山大火山岩省的原始岩浆具高镁（ MgO > 16％）特征。玄武岩的 REE反演计算揭示，参与峨眉山玄武岩岩浆作用的地幔具有异常高的潜能温度（ 1 550℃）。这些特征以及峨眉山玄武岩的大面积分布和一些熔岩所显示的类似于洋岛玄武岩 (OIB)的微量元素和 Sr- Nd同位素特征均为地幔热柱在能量和物质上参与峨眉山溢流玄武岩的形成提供了确凿证据。峨眉山两个主要岩类（高钛和低钛玄武岩）可能是不同地幔源区物质在不同条件下的熔融产物。低钛玄武岩形成于温度最高、岩石圈最薄的地幔柱轴部。地幔（ ISr≈ 0.705,ε Nd(t)≈＋ 2）熔融始于 140 km，并一直延续到较浅的深度（ 60 km,尖晶石稳定区 ),部分熔融程度为 16％，这类岩石可能代表了峨眉山玄武岩的主体。而高钛玄武岩的母岩浆的形成基本局限在石榴子石稳定区（ > 70 km），其源区特征为 : ISr≈ 0.704,ε Nd(t)≈＋ 5，可能代表了热柱边部或消亡期地幔小程度部分熔融（ 1.5％）的产物。     

Analytical solutions for Euler‐Bernoulli Beam on Pasternak foundation subjected to arbitrary dynamic loads

In this paper, the dynamic response of an infinite beam resting on a Pasternak foundation and subjected to arbitrary dynamic loads is developed in the form of analytical solution. The beam responses investigated are deflection, velocity, acceleration, bending moment, and shear force. The mechanical resistance of the Pasternak foundation is modeled using two parameters, that is, one accounts for soil resistance due to compressive strains in the soil and the other accounts for the resistance due to shear strains. Because the Winkler model only represents the compressive resistance of soil, comparatively, the Pasternak model is more realistic to consider shear interactions between the soil springs. The governing equation of the beam is simplified into an algebraic equation by employing integration transforms, so that the analytical solution for the dynamic response of the beam can be obtained conveniently in the frequency domain. Both inverse Laplace and inverse Fourier transforms combined with convolution theorem are applied to convert the solution into the time domain. The solutions for several special cases, such as harmonic line loads, moving line loads, and travelling loads are also discussed and numerical examples are conducted to investigate the influence of the shear modulus of foundation on the beam responses. The proposed solutions can be an effective tool for practitioners. Copyright © 2017 John Wiley & Sons, Ltd.     

Porosity of the melting zone and variations in the solid mantle upwelling rate beneath Hawaii: inferences from 238U-230Th-226Ra and 235U-231Pa disequilibria

Measurements of ²³⁸U-²³⁰Th-²²⁶Ra and ²³⁵U-²³¹Pa disequilibria in a suite of tholeiitic-to-basanitic lavas provide estimates of porosity, solid mantle upwelling rate and melt transport times beneath Hawaii. The observation that (²³⁰Th/²³⁸U) > 1 indicates that garnet is required as a residual phase in the magma sources for all of the lavas. Both chromatographic porous flow and dynamic melting of a garnet peridotite source can adequately explain the combined U-Th-Ra and U-Pa data for these Hawaiian basalts. For chromatographic porous flow, the calculated maximum porosity in the melting zone ranges from 0.3–3% for tholeiites and 0.1–1% for alkali basalts and basanites, and solid mantle upwelling rates range from 40 to 100 cm yr⁻¹ for tholeiites and from 1 to 3 cm yr⁻¹ for basanites. For dynamic melting, the escape or threshold porosity is 0.5–2% for tholeiites and 0.1–0.8% for alkali basalts and basanites, and solid mantle upwelling rates range from 10 to 30 cm yr⁻¹ for tholeiites and from 0.1 to 1 cm yr⁻¹ for basanites. Assuming a constant melt productivity, calculated total melt fractions range from 15% for the tholeiitic basalts to 3% for alkali basalts and basanites.     

Formation of the Jurassic Changboshan‐Xieniqishan highly fractionated I‐type granites,northeastern China: implication for the partial melting of juvenile crust induced by asthenospheric mantle upwelling

Zircon U–Pb ages, major and trace elements, and Sr, Nd and Hf isotope compositions of the Changboshan‐Xieniqishan (CX) intrusion from the Great Xing'an Range (GXAR), northeastern China, were studied to investigate its derivation, evolution and geodynamic significance. Laser ablation inductively coupled plasma mass spectrometry (LA‐ICP‐MS) zircon U–Pb dating yields an emplacement age of 161 ± 2 Ma for the CX intrusion. Bulk‐rock analyses show that this intrusion is characterized by high SiO₂, Na₂O and K₂O, but low MgO, CaO and P₂O₅. They are enriched in large‐ion lithophile elements and light rare earth elements, with marked Eu anomalies (mostly from 0.36 to 0.65), and depleted in heavy rare earth elements and high field strength elements. Most samples have relatively low (⁸⁷Sr/⁸⁶Sr)_i values (0.70423–0.70457), with ε_Nd(t) fluctuating between −0.4 and 2.3. The ε_Hf(t) for zircons varies from 5.4 to 8.7. Sr–Nd isotope modelling results, in combination with young Nd and Hf model ages (760–986 and 549–728 Ma, respectively) and the presence of relict zircons, indicate that the CX intrusion may originate from the partial melting of juvenile crust, with minor contamination of recycled crustal components, and then underwent extensive fractional crystallization of K‐feldspar, plagioclase, biotite, sphene, apatite, zircon and allanite. Considering the widespread presence of granitoids with coeval volcanic rocks, we contend that the CX intrusion formed in an extensional environment related to the upwelling of asthenospheric mantle induced by the subduction of the Palaeo‐Pacific plate, rather than a lithospheric delamination model. Copyright © 2013 John Wiley & Sons, Ltd.     

Analytical solution for two-dimensional dynamic consolidation in frequency domain

An analytical solution to the two-dimensional wave propagation in fluid-saturated half-space subjected to a strip load with vertical harmonic oscillation at the surface is presented. The basic equations have been derived on the basis of Biot's linear theory of poro-elasticity and then solved using Fourier complex transform for the horizontal direction. The importance of a number of soil characteristics including compressibility, degree of saturation and soil permeability has been examined. It is shown that the effect of pore fluid is dominant only for fully saturated soils with incompressible solid grains and low permeability. For partially saturated, compressible or very permeable soils, the stresses would be mainly transferred to solid part and there will be considerable reduction in pore pressure amplitude.     

Observations of Li isotopic variations in the Trinity Ophiolite: Evidence for isotopic fractionation by diffusion during mantle melting

The Trinity peridotite (northern CA) contains numerous lithologic sequences consisting of dunite to harzburgite to spinel lherzolite to plagioclase lherzolite. Previous workers have documented geochemical gradients in these sequences consistent with melt-rock reaction processes occurring around dunites, interpreted to reflect conduits for melt ascent. We have undertaken a study of Li isotope compositions of clinopyroxene and some olivine within these sequences using ion probe techniques to test the hypothesis that the geochemical gradients are related to diffusive fluxing of alkali elements into or away from the melt conduit.Results show large variations in ⁷Li/⁶Li occurring in a consistent pattern across three transects from dunite to plagioclase lherzolite within the Trinity peridotite. Specifically, measurements of average δ⁷Li for single thin sections along the traverse reveal a low in δ⁷Li in the harzburgite adjacent to the dunite returning to higher values farther from the dunite with a typical offset of ∼10 per mil in the low δ⁷Li trough. This pattern is consistent with a process whereby Li isotopes are fractionated during diffusion through a melt either from the dunite conduit to the surrounding peridotite, or from the surrounding peridotite into the dunite conduit. The patterns in ⁷Li/⁶Li occur over a length scale similar to the decrease in REE concentration in these same samples. Explaining both the trace element and Li isotopic gradients requires a combined process of alkali diffusion and melt extraction.We develop a numerical model and examine several scenarios of the combined diffusion-extraction process. Using experimentally constrained values for the change in Li diffusion coefficient with isotope mass, large changes in δ⁷Li as a function of distance can be created in year to decade timescales. The addition of the melt extraction term allows larger Li concentration gradients to be developed and thus produces larger isotopic fractionations than diffusion only models. The extraction aspect of the model can also account for the observed decrease in rare earth element concentrations across the transects.     

Experimental phase and melting relations of metapelite in the upper mantle: implications for the petrogenesis of intraplate magmas

We performed a series of piston-cylinder experiments on a synthetic pelite starting material over a pressure and temperature range of 3.0–5.0 GPa and 1,100–1,600°C, respectively, to examine the melting behaviour and phase relations of sedimentary rocks at upper mantle conditions. The anhydrous pelite solidus is between 1,150 and 1,200°C at 3.0 GPa and close to 1,250°C at 5.0 GPa, whereas the liquidus is likely to be at 1,600°C or higher at all investigated pressures, giving a large melting interval of over 400°C. The subsolidus paragenesis consists of quartz/coesite, feldspar, garnet, kyanite, rutile, ±clinopyroxene ±apatite. Feldspar, rutile and apatite are rapidly melted out above the solidus, whereas garnet and kyanite are stable to high melt fractions (>70%). Clinopyroxene stability increases with increasing pressure, and quartz/coesite is the sole liquidus phase at all pressures. Feldspars are relatively Na-rich [K/(K + Na) = 0.4–0.5] at 3.0 GPa, but are nearly pure K-feldspar at 5.0 GPa. Clinopyroxenes are jadeite and Ca-eskolaite rich, with jadeite contents increasing with pressure. All supersolidus experiments produced alkaline dacitic melts with relatively constant SiO₂ and Al₂O₃ contents. At 3.0 GPa, initial melting is controlled almost exclusively by feldspar and quartz, giving melts with K₂O/Na₂O ~1. At 4.0 and 5.0 GPa, low-fraction melting is controlled by jadeite-rich clinopyroxene and K-rich feldspar, which leads to compatible behaviour of Na and melts with K₂O/Na₂O ≫ 1. Our results indicate that sedimentary protoliths entrained in upwelling heterogeneous mantle domains may undergo melting at greater depths than mafic lithologies to produce ultrapotassic dacitic melts. Such melts are expected to react with and metasomatise the surrounding peridotite, which may subsequently undergo melting at shallower levels to produce compositionally distinct magma types. This scenario may account for many of the distinctive geochemical characteristics of EM-type ocean island magma suites. Moreover, unmelted or partially melted sedimentary rocks in the mantle may contribute to some seismic discontinuities that have been observed beneath intraplate and island-arc volcanic regions.     

Fluid inclusions in upper mantle xenoliths from Northern Patagonia,Argentina: Evidence for an upper mantle diapir

Summary Three generations of fluid inclusions can be recognized in upper mantle xenoliths from alkali basalts of the Somoncura Massif, Northern Patagonia, Argentina. The first (early, primary) one consists of dense CO₂ inclusions which were trapped in the mantle-crust boundary zone (22–36 km minimum trapping depth). Their co-genetic relationship with silicate melt inclusions enables us to constrain their minimum trapping temperature at 1200°C, indicating a high temperature event in a cooler environment. The late (pseudosecondary and secondary) generations of fluid inclusions were classified in accordance with their homogenization temperature to liquid CO₂ (L1) and vapor CO₂ (L2) phase. The minimum trapping depth for the first of the late inclusions (L1) is about 16 km. In spite of the uncertainties related to this value, L1 inclusions indicate that the upper mantle rocks, of which samples were delivered by the basalts, had some residence time in the middle crust where they experienced a metasomatic event. The fact that this event did not destroy the earlier inclusions, places severe constraints on its duration. The second late inclusions (L2) are low-pressure CO₂ inclusions with a minimum trapping depth of only 2 km, presumably a shallow magma chamber of the host basalts. The succession of fluid inclusions strongly points toward a fairly fast uprising upper mantle underneath Northern Patagonia. The petrology and mineral chemistry of the peridotitic xenoliths support this view. Extensive partial melting and loss of these melts is indicated by the preponderance of harzburgites in the upper mantle underneath Northern Patagonia, a fairly unusual feature for a continental upper mantle. That depletion event as well as several metasomatic events — including those which left traces of fluid inclusions — are possibly related to a high-speed diapiric uprise of the upper mantle in this area. The path can be traced from the garnet peridotite stability field into the middle crust, a journey which must have been unusually fast. Differences in rock, mineral, and fluid inclusion properties between geographic locations suggest a diffuse and differential type of diapirism. Future studies will hopefully help to map the full extent and the highs and lows of this diapir and elucidate questions related to its origin and future.

---

Fluid-Einschlüsse in Erdmantel-Xenolithen von Nord-Patagonien: Evidenz für einen Diapir im oberen Erdmantel

Zusammenfassung Erdmantel - Xenolithe in Alkali-Basalten des Somoncure Massivs, Nord-Patagonien, Argentinien, führen drei Generationen von Fluid-Einschlüssen. Die erste (frühe, primäre) Generation besteht aus dichten CO₂-Einschlüssen, welche offenbar in der Mantel-Kruste Grenzzone (22–36 km Minimum-Tiefe) eingeschlossen wurden. CO₂-Einschlüsse sind kogenetisch mit Silikat-Schmelzeinschlüssen. Dies erlaubt die Abschätzung der Einschließ-Temperatur mit minimal 1200°C, was auf ein Hochtemperatur-Ereignis in einer deutlich kühleren Umgebung hinweist. Die späten (pseudosekundäre und sekundäre) CO₂- Fluid-Einschlüsse bilden zwei Generationen von denen die eine in die flüssige (L1), die andere in die Dampfphase (L2) homogenisieren. Die minimale Einschließ-Tiefe für die L1 Generation ist etwa 16 km. Dies bedeutet - auch bei Berücksichtigung der mit diesem Wert verbundenen Ungenauigkeit - daß diese Erdmantel-Gesteine einige Zeit in der mittleren Erdkruste verbrachten und ein metasomatisches Ereignis erlebten, bevor sie von den Basalten zur Erdoberfläche gebracht wurden. Die Tatsache, daß dieses Ereignis die frühen Einschlüsse nicht zerstörte, kann nur bedeuten, daß es von kurzer Dauer war. Die L2-Generation besteht aus Niedrigdruck CO₂-Einschlüssen mit einer Minimum-Einschließtiefe von nur 2 km. Dies könnte in einer seichten Magmakammer des Wirt Basaltes geschehen sein.Die Abfolge von Fluid-Einschlüssen deutet auf einen relativ schnell aufsteigenden oberen Erdmantel unterhalb von Patagonien hin. Die Petrologie und Mineralchemie der peridotitischen Xenolithe unterstützen das. Die Vorherrschaft von Harzburgiten im Erdmantel unterhalb von Nord-Patagonien deutet auf umfangreiche Bildung partieller Schmelzen und deren Abfuhr hin — eine für einen kontinentalen Mantel ungewöhnliche Situation. Sowohl die Verarmungsereignisse, als auch die metasomatischen Veränderungen (einschließlich jene, welche Spuren in Form von Fluid Einschlüssen hinterließen) machen das Vorhandensein eines schnell aufsteigenden Daipirs im oberen Erdmantel dieser Gegend wahrscheinlich. Der Aufstieg kann vom Stabilitätsbereich der Granat-Peridotite bis in die mittlere Kruste verfolgt werden und muß daher relativ schnell erfolgt sein. Unterschiede in Gesteins-, Mineral und Fluid-Eigenschaften zwischen verschiedenen Lokalitäten legen einen diffusen und differenziellen Diapirismus nahe. Zukünftige Studien sollten es ermöglichen, das Gesamtausmaß und die unterschiedlichen Aufstiegshöhen des Diapirs zu kartieren und Hinweise auf seine Entstehung und zukünftige Entwicklung zu erhalten.

---

---

With 5 Figures