Lateral Variation in Crustal Structure along the Lesser Antilles Arc from Petrology of Crustal Xenoliths and Seismic Receiver Functions

We reconstruct crustal structure along the Lesser Antilles island arc using an inversion approach combining constraints from petrology of magmatic crustal xenoliths and seismic receiver functions. Xenoliths show considerable island-to-island variation in xenolith petrology from plagioclase-free ultramafic lithologies to gabbros and gabbronorites with variable proportions of amphibole, indicative of changing magma differentiation depths. Xenoliths represent predominantly cumulate compositions with equilibration depths in the range 5 to 40 km. We use xenolith mineral modes and compositions to calculate seismic velocities () and density at the estimated equilibration depths. We create a five-layer model of crustal structure for testing against receiver functions (RF) from island seismic stations along the arc. Lowermost layer (5) comprises peridotite with physical characteristics of mantle xenoliths from Grenada. Uppermost layer (1) consists of 5 km of volcaniclastics and sediments, whose physical properties are determined via a grid inversion routine. The three middle layers (2) to (4) comprise igneous arc crust with compositions corresponding to the xenoliths sampled at each island. By inversion we obtain a petrological best-fit for the RF on each island to establish the nature and thicknesses of layers (2) to (4). Along the arc we see variations in the depth and strength of both Moho and mid-crustal discontinuity (MCD) on length-scales of tens of km. Moho depths vary from 25 to 37 km; MCD from 11 and 32 km. The Moho is the dominant discontinuity beneath some islands (St. Kitts, Guadeloupe, Martinique, Grenada), whereas the MCD dominates beneath others (Saba, St. Eustatius). Along-arc variability in MCD depth and strength is consistent with variation in estimated magmatic H2O contents and differentiations depths that, in turn, influence xenolith lithologies. A striking feature is steep, along-arc gradients in similar to those observed at other oceanic arcs. These gradients reflect abrupt changes in rates and processes of magma generation in the underlying crust and mantle. We find no evidence for large, interconnected bodies of partial melt beneath the Lesser Antilles. Instead, the crustal velocity structure is consistent with magma differentiation in vertically-extensive, crystal mush-dominated reservoirs. Along-arc variation in crustal structure may reflect heterogeneous upwelling within the mantle wedge, itself driven by variation in slab-derived H2O fluxes.This work was supported by Natural Environment Research Council grants NE/N001966/1, NE/K004883/1, NE/K014978/1 and NE/K010824/1 and ERC Advanced Grant CRITMAG (Blundy

Melting during late-stage rifting in Afar is hot and deep

Investigations of a variety of continental rifts and margins worldwide have revealed that a considerable volume of melt can intrude into the crust during continental breakup, modifying its composition and thermal structure. However, it is unclear whether the cause of voluminous melt production at volcanic rifts is primarily increased mantle temperature or plate thinning. Also disputed is the extent to which plate stretching or thinning is uniform or varies with depth with the entire continental lithospheric mantle potentially being removed before plate rupture. Here we show that the extensive magmatism during rifting along the southern Red Sea rift in Afar, a unique region of sub-aerial transition from continental to oceanic rifting, is driven by deep melting of hotter-than-normal asthenosphere. Petrogenetic modelling shows that melts are predominantly generated at depths greater than 80 kilometres, implying the existence of a thick upper thermo-mechanical boundary layer in a rift system approaching the point of plate rupture. Numerical modelling of rift development shows that when breakup occurs at the slow extension rates observed in Afar, the survival of a thick plate is an inevitable consequence of conductive cooling of the lithosphere, even when the underlying asthenosphere is hot. Sustained magmatic activity during rifting in Afar thus requires persistently high mantle temperatures, which would allow melting at high pressure beneath the thick plate. If extensive plate thinning does occur during breakup it must do so abruptly at a late stage, immediately before the formation of the new ocean basin

Clinopyroxene diversity and magma plumbing system processes in an accreted Pacific ocean island, Panama

Characterising equilibrium and disequilibrium crystal-melt processes is critical in determining the extent of magma mixing and crystallization conditions in the roots of volcanoes. However, these processes remain poorly investigated in most Pacific intraplate ocean settings that are difficult to access and study. To help address this issue, we investigated crystallization conditions of clinopyroxene phenocrysts in an accreted Palaeogene oceanic island in Panama. Petrographic and geochemical observations, petrological modelling of major and trace elements, and liquid-mineral multicomponent equilibrium tests were carried out using basalts, picrites, and hawaiites of the transitional tholeiitic shield to alkaline post-shield volcanic stages of the island. Five types of clinopyroxene crystals were identified, including (1) microphenocrysts with micron-scale oscillatory zoning, (2) primitive, yet resorbed picrite-hosted phenocrysts, (3) chemically homogeneous, anhedral crystals found in the remaining basalts, (4) Ti–rich euhedral hawaiite-hosted phenocrysts, and (5) evolved sector-zoned phenocrysts. Liquid-clinopyroxene multicomponent equilibrium tests in combination with textural analysis show that ~ 74% of the studied clinopyroxenes are in possible major element equilibrium with one of the available whole rock magma compositions, of which only 21% are equilibrated with their carrier liquid. To deconvolute clinopyroxene-melt pairings and determine plumbing system conditions, we combine rhyolite-MELTS modelling, geothermobarometry, and major- and trace-element equilibrium evaluations, limiting crystallization conditions to crustal levels (< 23 km depth). No migration of magmatic reservoirs to deeper levels is observed during the shield- to post-shield transition. These results suggest the occurrence of an extensive crystal mush system during the late shield to post-shield volcanic stages of this intraplate volcanic system, with both primitive and evolved crystallization domains sampled during eruptions

Quantifying garnet-melt trace element partitioning using lattice-strain theory: New crystal-chemical and thermodynamic constraints

Many geochemical models of major igneous differentiation events on the Earth, the Moon, and Mars invoke the presence of garnet or its high-pressure majoritic equivalent as a residual phase, based on its ability to fractionate critical trace element pairs (Lu/Hf, U/Th, heavy REE/light REE). As a result, quantitative descriptions of mid-ocean ridge and hot spot magmatism, and lunar, martian, and terrestrial magma oceans require knowledge of garnet-melt partition coefficients over a wide range of conditions. In this contribution, we present new crystal-chemical and thermodynamic constraints on the partitioning of rare earth elements (REE), Y and Sc between garnet and anhydrous silicate melt as a function of pressure (P), temperature (T), and composition (X). Our approach is based on the interpretation of experimentally determined values of partition coefficients D using lattice-strain theory. In this and a companion paper (Draper and van Westrenen this issue) we derive new predictive equations for the ideal ionic radius of the dodecahedral garnet X-site,

Interaction between protokimberlite melts and mantle lithosphere: evidence from mantle xenoliths from the Dalnyaya kimberlite pipe, Yakutia (Russia)

The Dalnyaya kimberlite pipe (Yakutia, Russia) contains mantle peridotite xenoliths (mostly lherzolites and harzburgites) that show both sheared porphyroclastic (deformed) and coarse granular textures, together with ilmenite and clinopyroxene megacrysts. Deformed peridotites contain high-temperature Fe-rich clinopyroxenes, sometimes associated with picroilmenites, which are products of interaction of the lithospheric mantle with protokimberlite related melts. The orthopyroxene-derived geotherm for the lithospheric mantle beneath Dalnyaya is stepped similar to that beneath the Udachnaya pipe. Coarse granular xenoliths fall on a geotherm of 35 mWm-2 whereas deformed varieties yield a 45 mWm-2 geotherm in the 2–7.5 GPa pressure interval. The chemistry of the constituent minerals including garnet, olivine and clinopyroxene shows trends of increasing Fe# (= Fe/(Fe+Mg) with decreasing pressure. This may suggest that the interaction with fractionating protokimberlite melts occurred at different levels. Two major mantle lithologies are distinguished by the trace element patterns of their constituent minerals, determined by LA-ICP-MS. Orthopyroxenes, some clinopyroxenes and rare garnets are depleted in Ba, Sr, HFSE and MREE and represent relic lithospheric mantle. Re-fertilized garnet and clinopyroxene are more enriched. The distribution of trace elements between garnet and clinopyroxene shows that the garnets dissolved primary orthopyroxene and clinopyroxene. Later high temperature clinopyroxenes related to the protokimberlite melts partially dissolved these garnets. Olivines show decreases in Ni and increases in Al, Ca and Ti from Mg-rich varieties to the more Fe-rich, deformed and refertilized ones. Minerals showing higher Fe# (0.11–0.15) are found within intergrowths of low-Cr ilmenite-clinopyroxene-garnet related to the crystallization of protokimberlite melts in feeder channels. In P-f(O2) diagrams, garnets and Cr-rich clinopyroxenes indicate reduced conditions at the base of the lithosphere at -5 log units below a FMQ buffer. However, Cr-poor clinopyroxenes, together with ilmenite and some Fe-Ca-rich garnets, demonstrate a more oxidized trend in the lower part of lithosphere at -2 to 0 log units relative to FMQ. Clinopyroxenes from xenoliths in most cases show conditions transitional between those determined for garnets and megacrystalline Cr-poor suite. The relatively low diamond grade of Dalnyaya kimberlites is explained by a high degree of interaction with the oxidized protokimberlite melts, which is greater at the base of the lithosphere

First recorded eruption of Nabro volcano, Eritrea, 2011

We present a synthesis of diverse observations of the first recorded eruption of Nabro volcano, Eritrea, which began on 12 June 2011. While no monitoring of the volcano was in effect at the time, it has been possible to reconstruct the nature and evolution of the eruption through analysis of re- gional seismological and infrasound data and satellite remote sensing data, supplemented by petrological analysis of erupted products and brief field surveys. The event is notable for the comparative rarity of recorded historical eruptions in the region and of caldera systems in general, for the prodi- gious quantity of SO2 emitted into the atmosphere and the significant human impacts that ensued notwithstanding the low population density of the Afar region. It is also relevant in understanding the broader magmatic and tectonic signifi- cance of the volcanic massif of which Nabro forms a part and which strikes obliquely to the principal rifting directions in the Red Sea and northern Afar. The whole-rock compositions of Editorial responsibility: G. Giordano the erupted lavas and tephra range from trachybasaltic to trachybasaltic andesite, and crystal-hosted melt inclusions contain up to 3,000 ppm of sulphur by weight. The eruption was preceded by significant seismicity, detected by regional networks of sensors and accompanied by sustained tremor. Substantial infrasound was recorded at distances of hundreds to thousands of kilometres from the vent, beginning at the onset of the eruption and continuing for weeks. Analysis of ground deformation suggests the eruption was fed by a shal- low, NW–SE-trending dike, which is consistent with field and satellite observations of vent distributions. Despite lack of prior planning and preparedness for volcanic events in the country, rapid coordination of the emergency response miti- gated the human costs of the eruption

The evolution and storage of primitive melts in the Eastern Volcanic Zone of Iceland: the 10 ka Grímsvötn tephra series (i.e. the Saksunarvatn ash)

Major, trace and volatile elements were measured in a suite of primitive macrocrysts and melt inclusions from the thickest layer of the 10 ka Grímsvötn tephra series (i.e. Saksunarvatn ash) at Lake Hvítárvatn in central Iceland. In the absence of primitive tholeiitic eruptions (MgO > 7 wt.%) within the Eastern Volcanic Zone (EVZ) of Iceland, these crystal and inclusion compositions provide an important insight into magmatic processes in this volcanically productive region. Matrix glass compositions show strong similarities with glass compositions from the AD 1783–84 Laki eruption, confirming the affinity of the tephra series with the Grímsvötn volcanic system. Macrocrysts can be divided into a primitive assemblage of zoned macrocryst cores (An_78–An_92, Mg#_cpx = 82–87, Fo_79.5–Fo_87) and an evolved assemblage consisting of unzoned macrocrysts and the rims of zoned macrocrysts (An_60–An_68, Mg#_cpx = 71–78, Fo_70–Fo_76). Although the evolved assemblage is close to being in equilibrium with the matrix glass, trace element disequilibrium between primitive and evolved assemblages indicates that they were derived from different distributions of mantle melt compositions. Juxtaposition of disequilibrium assemblages probably occurred during disaggregation of incompatible trace element-depleted mushes (mean La/Yb_melt = 2.1) into aphyric and incompatible trace element-enriched liquids (La/Yb_melt = 3.6) shortly before the growth of the evolved macrocryst assemblage. Post-entrapment modification of plagioclase-hosted melt inclusions has been minimal and high-Mg# inclusions record differentiation and mixing of compositionally variable mantle melts that are amongst the most primitive liquids known from the EVZ. Coupled high field strength element (HFSE) depletion and incompatible trace element enrichment in a subset of primitive plagioclase-hosted melt inclusions can be accounted for by inclusion formation following plagioclase dissolution driven by interaction with plagioclase-undersaturated melts. Thermobarometric calculations indicate that final crystal-melt equilibration within the evolved assemblage occurred at ~1140°C and 0.0–1.5 kbar. Considering the large volume of the erupted tephra and textural evidence for rapid crystallisation of the evolved assemblage, 0.0–1.5 kbar is considered unlikely to represent a pressure of long-term magma accumulation and storage. Multiple thermometers indicate that the primitive assemblage crystallised at high temperatures of 1240–1300°C. Different barometers, however, return markedly different crystallisation depth estimates. Raw clinopyroxene-melt pressures of 5.5–7.5 kbar conflict with apparent melt inclusion entrapment pressures of 1.4 kbar. After applying a correction derived from published experimental data, clinopyroxene-melt equilibria return mid-crustal pressures of 4±1.5 kbar, which are consistent with pressures estimated from the major element content of primitive melt inclusions. Long-term storage of primitive magmas in the mid-crust implies that low CO_2 concentrations measured in primitive plagioclase-hosted inclusions (262–800 ppm) result from post-entrapment CO_2 loss during transport through the shallow crust. In order to reconstruct basaltic plumbing system geometries from petrological data with greater confidence, mineral-melt equilibrium models require refinement at pressures of magma storage in Iceland. Further basalt phase equilibria experiments are thus needed within the crucial 1–7 kbar range.D.A.N. was supported by a Natural Environment Research Council studentship (NE/1528277/1) at the start of this project. SIMS analyses were supported by Natural Environment Research Council Ion Microprobe Facility award (IMF508/1013).This is the final version of the article. It first appeared from Springer via http://dx.doi.org/10.1007/s00410-015-1170-

Variable water input controls evolution of the Lesser Antilles volcanic arc

Oceanic lithosphere carries volatiles, notably water, into the mantle through subduction at convergent plate boundaries. This subducted water exercises control on the production of magma, earthquakes, formation of continental crust and mineral resources. Identifying different potential fluid sources (sediments, crust and mantle lithosphere) and tracing fluids from their release to the surface has proved challenging1. Atlantic subduction zones are a valuable endmember when studying this deep water cycle because hydration in Atlantic lithosphere, produced by slow spreading, is expected to be highly non-uniform2. Here, as part of a multi-disciplinary project in the Lesser Antilles volcanic arc3, we studied boron trace element and isotopic fingerprints of melt inclusions. These reveal that serpentine—that is, hydrated mantle rather than crust or sediments—is a dominant supplier of subducted water to the central arc. This serpentine is most likely to reside in a set of major fracture zones subducted beneath the central arc over approximately the past ten million years. The current dehydration of these fracture zones coincides with the current locations of the highest rates of earthquakes and prominent low shear velocities, whereas the preceding history of dehydration is consistent with the locations of higher volcanic productivity and thicker arc crust. These combined geochemical and geophysical data indicate that the structure and hydration of the subducted plate are directly connected to the evolution of the arc and its associated seismic and volcanic hazards