The effect of bulk rock composition on the stability of amphibole in the upper mantle: Implications for solidus positions and mantle metasomatism

Summary The stability of pargasitic amphibole in the upper mantle is a function of water content and bulk rock composition, and under water-undersaturated conditions, the stability of amphibole controls the solidus position. Experiments in the system Tinaquillo peridotite +0.2% H₂O, a refractory peridotite under water-undersaturated conditions, show that amphibole is stable to 1030°C and 26 kb. In contrast, pargasitic amphibole is stable to 1150°C and 30 kb in Hawaiian pyrolite, a more fertile peridotite composition. This indicates that under water-undersaturated conditions, the most fertile part of a crystallizing mantle diapir with an inhomogeneous composition will solidify first while a more refractory component will contain an alkali-rich melt which will have the ability to metasomatize adjacent regions. The relative stabilities of amphibole in refractory and fertile bulk compositions may result in increasing rather than diminishing chemical contrasts in high temperature lherzolite, i.e. a process of metamorphic differentiation. Ti, Fe, Al and Na metasomatism can therefore be considered a normal occurrence associated with the upward migration and solidification of an H₂O-bearing mantle diapir.

---

Der Einfluß der Gesamtgesteins-Zusammensetzung auf die Stabilität von Amphibol im oberen Mantel: Bedeutung für Solidus-Positionen und Mantel-Metasomatose

Zusammenfassung Die Stabilität pargasitischer Amphibole im oberen Mantel ist eine Funktion von Wassergehalt und Gesamtgesteins-Zusammensetzung. Unter wasser-untersättigten Bedingungen, kontrolliert die Stabilität von Amphibol die Solidus-Position. Experimente in dem System Tinaquillo Peridotit +0,2% H₂O, einem refraktären Peridotit unter wasser-untersättigten Bedingungen, zeigen daß Amphibol bis 1030°C und 26 Kb stabil ist. Im Gegensatz dazu ist pargasitische Hornblende in einem Hawaii-Pyrolit, von mehr fertiler Peridotit-Zuammensetzung, bis 1150°C und 30 Kb stabil. Das zeigt, daß bei wasser-untersättigten Bedingungen der am meisten produktive Teil eines kristallisierenden Mantel-Diapirs mit inhomogener Zusammensetzung sich zuerst verfestigen wird, während eine mehr refraktäre Komponente eine alkali-reiche Schmelze enthalten wird, die wiederum die Fähigkeit hat, umliegende Bereiche metasomatisch zu beeinflussen. Die relativen Stabilitäten von Amphibol in refraktären und fertilen Gesamtzusammensetzungen können dazu führen, daß die chemischen Gegensätze in Hochtemperaturlherzoliten eher zunehmen als abnehmen, d. h. ein Prozeß metamorpher Differentiation. Ti, Fe, Al und Na Metasomatose können deshalb als ein verbreiteter Vorgang, der mit der Aufwärtsbewegung und Verfestigung eines H₂O-führenden Mantel-Diapirs assoziiert ist, betrachtet werden.

---

---

With 4 Figures     

Late Palaeozoic hydrocarbon migration through the Clair field, West of Shetland, UK Atlantic margin

Geochemical analysis of bitumen- and hydrocarbon-bearing fluid inclusions from the Devonian-Carboniferous Clair field indicates that the reservoirs contain a mixture of oils from different marine and lacustrine sources. Reconstruction of the Clair field oil-charge history using fluid inclusion petrography show that oil-charging occurred at times of K-feldspar, quartz and calcite cementation. Temperature–composition–time data yielded from the integration of fluid inclusion microthermometry with high-resolution Ar–Ar dating, date hydrocarbon-bearing K-feldspar overgrowths at 247 ± 3.3 Ma. These data show that in order for oil to be trapped within primary fluid inclusions in K-feldspar overgrowths, hydrocarbon migration throughout the UK Atlantic margin must have been taking place during the Late Palaeozoic and as such, current industry oil-play models based solely on oil charging from Jurassic-Cretaceous marine sources are clearly incomplete and need revision. Apatite fission track analysis and vitrinite reflectance data were used to reconstruct thermal burial histories and assess potential oil generation from Middle Devonian lacustrine source rocks. Thermal history data from wells along The Rona Ridge adjacent to the Clair field show that the Palaeozoic section was heated to greater than 100 °C at some time between 270 and 230 Ma, confirming that Devonian source rocks were mature and expelling oil during the Late Palaeozoic at the time that authigenic K-feldspar overgrowths were growing in the Clair field.     

Trace Element Partitioning and Accessory Phase Saturation during H2O-Saturated Melting of Basalt with Implications for Subduction Zone Chemical Fluxes

Experimental phase equilibrium and trace element partitioningdata are reported for H₂O-saturated mid-ocean ridge basalt at2·5 GPa, 750–900°C and oxygen fugacities atthe nickel–nickel oxide buffer. Garnet, omphacite andrutile are present at all temperatures. Amphibole and epidotedisappear as residual phases above 800°C; allanite appearsabove 750°C. The Na–Al-rich silicate glass presentin all run products is likely to have quenched from a supercriticalliquid. Trace element analyses of glasses demonstrate the importantcontrol exerted by residual minerals on liquid chemistry. Inaddition to garnet, which controls heavy rare earth elements(HREE) and Sc, and rutile, which controls Ti, Nb and Ta, allanitebuffers the light REE (LREE; La–Sm) contents of liquidsto relatively low levels and preferentially holds back Th relativeto U. In agreement with previous experimental and metamorphicstudies we propose that residual allanite plays a key role inselectively retaining trace elements in the slab during subduction.Experimental data and analyses of allanite-bearing volcanicrocks are used to derive a model for allanite solubility inliquids as a function of pressure, temperature, anhydrous liquidcomposition and LREE content. The large temperature dependenceof allanite solubility is very similar to that previously determinedfor monazite. Our model, fitted to 48 datapoints, retrievesLREE solubility (in ppm) to within a factor of 1· 40over a pressure range of 0–4 GPa, temperature range of700–1200°C and for liquids with anhydrous SiO₂ contentsof 50–84 wt %. This uncertainty in LREE content is equivalentto a temperature uncertainty of only ± 27°C at 1000K, indicating the potential of allanite as a geothermometer.Silicic liquids from either basaltic or sedimentary protolithswill be saturated in allanite except for Ca-poor protolithsor at very high temperatures. For conventional subduction geothermsthe low solubility of LREE (+ Th) in liquids raises questionsabout the mechanism of LREE + Th transport from slab to wedge.It is suggested either that, locally, temperatures experiencedby the slab are high enough to eliminate allanite in the residueor that substantial volumes of H₂O-rich fluids must pass throughthe mantle wedge prior to melting. The solubility of accessoryphases in fluids derived from subducted rocks can provide importantconstraints on subduction zone thermal structure. KEY WORDS: subduction; experimental petrology; allanite; solubility; supercritical liquid; eclogite     

Phase Relations and Melting of Anhydrous K-bearing Eclogite from 1200 to 1600{degrees}C and 3 to 5 GPa

To investigate eclogite melting under mantle conditions, wehave performed a series of piston-cylinder experiments usinga homogeneous synthetic starting material (GA2) that is representativeof altered mid-ocean ridge basalt. Experiments were conductedat pressures of 3·0, 4·0 and 5·0 GPa andover a temperature range of 1200–1600°C. The subsolidusmineralogy of GA2 consists of garnet and clinopyroxene withminor quartz–coesite, rutile and feldspar. Solidus temperaturesare located at 1230°C at 3·0 GPa and 1300°C at5·0 GPa, giving a steep solidus slope of 30–40°C/GPa.Melting intervals are in excess of 200°C and increase withpressure up to 5·0 GPa. At 3·0 GPa feldspar, rutileand quartz are residual phases up to 40°C above the solidus,whereas at higher pressures feldspar and rutile are rapidlymelted out above the solidus. Garnet and clinopyroxene are theonly residual phases once melt fractions exceed 20% and garnetis the sole liquidus phase over the investigated pressure range.With increasing melt fraction garnet and clinopyroxene becomeprogressively more Mg-rich, whereas coexisting melts vary fromK-rich dacites at low degrees of melting to basaltic andesitesat high melt fractions. Increasing pressure tends to increasethe jadeite and Ca-eskolaite components in clinopyroxene andenhance the modal proportion of garnet at low melt fractions,which effects a marked reduction in the Al₂O₃ and Na₂O contentof the melt with pressure. In contrast, the TiO₂ and K₂O contentsof the low-degree melts increase with increasing pressure; thusNa₂O and K₂O behave in a contrasted manner as a function ofpressure. Altered oceanic basalt is an important component ofcrust returned to the mantle via plate subduction, so GA2 maybe representative of one of many different mafic lithologiespresent in the upper mantle. During upwelling of heterogeneousmantle domains, these mafic rock-types may undergo extensivemelting at great depths, because of their low solidus temperaturescompared with mantle peridotite. Melt batches may be highlyvariable in composition depending on the composition and degreeof melting of the source, the depth of melting, and the degreeof magma mixing. Some of the eclogite-derived melts may alsoreact with and refertilize surrounding peridotite, which itselfmay partially melt with further upwelling. Such complex magma-genesisconditions may partly explain the wide spectrum of primitivemagma compositions found within oceanic basalt suites. KEY WORDS: eclogite; experimental petrology; mafic magmatism; mantle melting; oceanic basalts     

青藏高原西部赛利普中新世火山岩源区:地球化学及Sr-Nd同位素制约

青藏高原拉萨地块西部赛利普地区新生代火山岩依据主量元素可划分为超钾质、钾质和钙碱性系列,主要的岩石类型为粗面安山岩、粗面岩,一个超钾质岩石的40Ar-39Ar年龄为17.58Ma,指示出火山活动为中新世.超钾质、钾质和钙碱性火山岩都显示出富集LREE及LILE(Th、U)、亏损HFSE(Nb、Ta、Ti)的特征.超钾质火山岩具有较高的K2O(6.31%～8.55%)、MgO(6.75%～8.96%)、Cr(270.7×10-6～460.4×10-6)、Ni(142.3×10-6～233.9×10-6)含量,较高的(87Sr/86Sr)i(0.71883～0.72732)和较低的εNd(-14.78～-15.37),指示可能起源于一个前期亏损并经后期俯冲作用改造的富钾的方辉橄榄岩富集地幔源区.钾质火山岩具有比超钾质火山岩低的K2O、MgO、Cr、Ni含量以及高的Ba、Sr含量,初始87Sr/86Sr为0.71553～0.71628,初始143Nd/144Nd为0.51197～0.51198,在空间上与超钾质火山岩共生,可能是前者母岩浆的演化产物.钙碱性火山岩具有较高的Sr(881.7×10-6～1309.2×10-6)、Sr/Y比值(50～108)和较低的Y(12.05×10-6～18.02×10-6),明显亏损重稀土Yb(0.93×10-6～1.30×10-6),类似于典型的埃达克质岩成分特征但相对高钾,并具有相对低的(87Sr/86Sr);(0.70928～0.71374)以及高的εNd(-7.90～-10.91),指示起源于富钾增厚下地壳物质的部分熔融.区域上拉萨地块超钾质岩、钾质岩与N-S向地堑系在空间上共存、时间上相吻合,由此本文认为拉萨地块中新世钾质.超钾质岩和南北向地堑系的形成可能与中新世早期北向俯冲的印度大陆岩石圈断离有关.     

青藏高原东南木里地区二叠纪苦橄岩的Os-Sr-Nd同位素地球化学研究

本文报道了在青藏高原东南木里地区发现的二叠纪苦橄岩和与其共生玄武岩的主微量元素地球化学特征以及Os-Sr-Nd同位素组成。苦橄岩和与其共生玄武岩受地壳混染作用影响较小。根据苦橄岩的Ti/Y比值和初始的Os同位素组成,将木里苦橄岩分为两类:高Ti/Y型苦橄岩和低Ti/Y型苦橄岩,其中高Ti/Y型苦橄岩具有高的γ_(Os)= 5.3～ 10.7和ε_(Nd)= 5.9～ 6.4,与全球典型洋岛玄武岩的Os和Nd同位素组成接近,代表了地幔柱源区的同位素特征;而低Ti/Y型苦橄岩具有低的γ_(Os)=-4.1～ 1.2和ε_(Nd)= 3.2～ 5.0,可能表明受到了SCLM(大陆岩石圈地幔)源区物质的混染。与其共生的玄武岩具有低的γ_(Os)=-3.5～-1.6和ε_(Nd)=-0.6～ 0.7,表明其来自于不同于低Ti/Y型苦橄岩也有异于高Ti/Y型苦橄岩的地幔源区,但是也可能受到了SCLM物质的混染。基于Nd-Os同位素的地幔柱与SCLM的二端元混合模型显示:低Ti/Y型苦橄岩可能是SCLM物质组分与地幔柱起源的苦橄质原始岩浆混合形成的;与苦橄岩共生的玄武岩可能是由地幔柱来源的玄武质岩浆与SCLM小比例熔融的熔体混合形成的。     

冈底斯新近纪钾质火山作用:消减沉积物折返的地球化学与Sr-Nd-Pb同位素证据

在冈底斯岩浆岩带的羊八井地区,林子宗火山岩系上部出现小规模的粗面质火山岩.岩石学与地球化学研究表明粗面岩与时空密切相关的大体积安山质-英安质-流纹质火山岩属于不同的火山岩系:碱性钾质系列与钙碱性系列.详细的地球化学研究证明林子宗晚期小体积钾质岩具有独立的岩浆源区,而非早期钙碱性系统低压岩浆过程演化的产物.粗面质熔岩SiO2含量为62.91%～64.63%,具有高K2O(7.52%～8.05%)、Al2O3(16.64%-17.35%)、低TiO2(0.59%～0.68%)与MsO(0.15%～0.77%).钾质熔岩富集LILE与LREE,亏损HFSE,具有高Rb/Sr(1.1～2.3)、Th/La(0.59～0.80)、Th/U(6.3～10.9)和低Ce/Pb(4.1～6.2)与Ba/Th(4.3～14.5)比值,其同位素组成变化范围分别是87Sr/86Sr(I)=0.7068～0.7075、143Nd/144Nd(I)=0.51241～0.51252、206Pb/204Pb=18.87～18.95、207Pb/204Pb=15.63～15.70和208Pb/204Pb=39.24～39.68.粗面岩的地球化学与同位素特征表明冈底斯新近纪钾质火山岩来源于角闪岩相的消减沉积物熔融.因此,羊八井新近纪钾质火山熔岩提供了印度-亚洲大陆碰撞早期消减沉积物折返的证据.     

Effect of Cyclic Loading Frequency on Undrained Behaviors of Undisturbed Marine Clay

Based on a series of cyclic triaxial tests, the effect of cyclic frequency on the undrained behaviors of undisturbed marine clay is investigated. For a given dynamic stress ratio, the accumulated pore water pressure and dynamic strain increase with the number of cycles. There exists a threshold value for beth the accumulated pore water pressure and dynamic strain, below which the effect of cyclic frequency is very small, but above which the accumulated pore water pressure and dynamic strain increase intensely with the decrease of cyclic frequency for a given number of cycles. The dynamic strength increases with the increase of cyclic frequency, whereas the effect of cyclic frequency on it gradually diminishes to zero when the number of cycles is large enough, and the dynamic strengths at different frequencies tend to the same limiting minimum dynamic strength. The test results demonstrate that the reasons for the frequency effect on the undrained soil behaviors are beth the creep effect induced by the loading rate and the decrease of sample effective confining pressure caused by the accumulated pore water pressure.     

Wave Decay and Its Characteristics in Surf Zone

In accordance with the similarity between breaking waves and hydraulic jumps, the expressions for estimating wave decay and wave energy dissipation in the surf zone are derived based on the fundamental equations of fluid mechanics. Using the numerical solution of cnoidal wave theory, the various kinematic properties of waves in the surf zone, including the relative wave crest height, wave energy, and radiation stress are discussed. The values calculated with the method proposed in this paper are in good agreement with the experimental data gained by other researchers. The present expressions can be used in the studies of sediment transport on gently sloping beaches, especially on muddy beaches.     

A New Pursuit Mode of Nonlinear Wave Surface

The numerical mode of nonlinear wave transformation based on both the Laplace equation for water field and the Bemoulli equation for water surface is a kind of time-domain boundary problem with initial conditions. And the basis for establishing the numerical mode of nonlinear wave in time domain is to trace the position of wave free surface and to calculale the instantaneous surface height and surface potential function. This paper firstly utilizes the ‘0-1‘ combined BEM to separate the boundary by means of discretization of Green‘ s integral equation based on the Laplace equation, then separates the free surface of wave with FEM and derives the FEM equation of wave surface that satisfies the nonlinear boundary conditions. By jointly solving the above BEM and FEM equations, the wave potential and surface height could be obtained with iteration in time domain. Thus a new kind of nonlinear numerical mode is established for calculating wave transformation. The wave test in the numerical wave tank shows that the numerical simulation with this mode is of high accuracy.