Water, oceanic fracture zones and the lubrication of subducting plate boundaries - insights from seismicity

We investigate the relationship between subduction processes and related seismicity for the Lesser Antilles Arc using the Gutenberg-Richter law. This power lawdescribes the earthquakemagnitude distribution, with the gradient of the cumulative magnitude distribution being commonly known as the b-value. The Lesser Antilles Arc was chosen because of its alongstrike variability in sediment subduction and the transition from subduction to strike-slip movement towards its northern and southern ends. The data are derived from the seismicity catalogues from the Seismic Research Centre of The University of the West Indies and the Observatoires Volcanologiques et Sismologiques of the Institut de Physique du Globe de Paris and consist of subcrustal events primarily from the slab interface. The b-value is found using a Kolmogorov-Smirnov test for a maximum-likelihood straight line-fitting routine. We investigate spatial variations in b-values using a grid-search with circular cells as well as an along-arc projection. Tests with different algorithms and the two independent earthquake cataloges provide confidence in the robustness of our results. We observe a strong spatial variability of the b-value that cannot be explained by the uncertainties. Rather than obtaining a simple north-south b-value distribution suggestive of the dominant control on earthquake triggering being water released from the sedimentary cover on the incoming American Plates, or a b-value distribution that correlates with on the obliquity of subduction, we obtain a series of discrete, high b-value 'bull's-eyes' along strike. These bull's-eyes, which indicate stress release through a higher fraction of small earthquakes, coincide with the locations of known incoming oceanic fracture zones on the American Plates. We interpret the results in terms of water being delivered to the Lesser Antilles subduction zone in the vicinity of fracture zones providing lubrication and thus changing the character of the related seismicity. Our results suggest serpentinization around mid-ocean ridge transform faults, which go on to become fracture zones on the incoming plate, plays a significant role in the delivery of water into the mantle at subduction zones

Experimental evidence for polybaric differentiation of primitive arc basalt beneath St. Vincent, Lesser Antilles

Equilibrium crystallization experiments have been performed on a primitive high-MgO basalt (HMB) from Soufrière, St. Vincent, with three initial H2O contents (0·6, 2·3 and 4·5 wt %), at pressures of 0·4, 0·7, 1·0 and 1·3 GPa and temperatures from 1350 to 950°C. Redox conditions, as determined by µXANES analysis of Fe3+ in experimental glasses, were 1–4 log units above the nickel–nickel oxide (NNO) buffer. The aim of the study was to explore the differentiation conditions that gave rise to the observed geochemical variation in lavas and plutonic (cumulate) xenoliths from St. Vincent. An experiment with 4·5 wt % initial H2O is multiply saturated close to its liquidus (1180°C and 1·3 GPa) with a spinel lherzolite assemblage, which is consistent with a primary origin for HMB in the mantle wedge. Multiple saturation of HMB with 2·3 wt % H2O was not observed, but is inferred to occur at pressures >1·3 GPa. The experimental results show that initial H2O content has significant influence on differentiation paths of primary HMB magma, with different lava varieties generated under discrete, well-constrained P–T–H2O conditions. Low-magnesian basalts (LMB) can be generated from HMB with 2·3–4·5 wt % H2O at pressures of 1·0–1·3 GPa, corresponding to Moho depths beneath St. Vincent. The CaO contents of LMB are sensitive to differentiation pressure: high-CaO LMB are produced at pressures >0·5 GPa. Basaltic andesites (BA) can be generated at 0·7–1·0 GPa from HMB with 0·6–2·3 wt % H2O. High-alumina basalts (HAB) are produced at mid- to upper-crustal conditions (≤0·4 GPa) by differentiation of HMB with high initial H2O (≥4 wt %) through delay of plagioclase crystallization and dominant fractionation of olivine, clinopyroxene and spinel. St. Vincent andesites could be produced from relatively dry (≤0·6 wt % H2O) HMB only at lower-crustal conditions. This is suggestive of a partial melting origin from precursor HMB that had solidified at depth to produce gabbros with ∼30% hornblende (i.e. ∼0·6 wt % structurally bound H2O). The experimentally determined differentiation conditions are consistent with polybaric differentiation within a hot zone that extends from the Moho and uppermost mantle to the mid- or upper crust. Within the hot zone differentiation occurs by a combination of crystallization of HMB with 2–5 wt % H2O and partial melting of ancestral HMB gabbros. Although the experimental melts provide an excellent match to erupted lava compositions, experimental crystal compositions do not match either phenocrysts or cumulate crystals, as preserved in xenoliths. The failure to reproduce natural crystal compositions suggests that these are formed as differentiated magmas ascend and attain their H2O-saturated liquidi at shallower pressures. Thus there is a disconnect between the high-pressure phase compositions and assemblages that generate liquid compositional diversity and the low-pressure composition and assemblages that occur as phenocrysts and in cumulate xenoliths. This finding lends support to the idea of cryptic fractionation in the generation of arc magmas

Ascent of volatile-rich felsic magma in dykes:a numerical model applied to deep-sourced porphyry intrusions

Dyke propagation is a mechanism for more rapid ascent of felsic magmas through the crust than is possible via diapirs or percolative flow. As it ascends, the magma undergoes complex physical and chemical transformations induced by decompression and cooling. These processes dramatically change the magma density and viscosity, which in turn affect magma ascent rate and the depth at which the dyke arrests. We present a mathematical model of dyke propagation for silicic magmas taking into account the presence of multiple volatile species (H2O and CO2), bubble growth, heat advection and loss, crystallization and latent heat release. We consider conditions for dykes associated with porphyry ore deposits, which may represent an end-member in rapid ascent of felsic magmas from depth. In particular, we simulate the propagation of dykes launched from a deep (900 MPa), volatile-saturated magma source, testing the effects of the magma H2O/CO2 content, temperature and mass on its ascent rate and final emplacement depth. The model predicts short ascent times (hours to days), with a large increase in viscosity at shallow depth, leading to stagnation and solidification of the dyke. Higher initial water content, higher temperature and larger mass of the magma in the dyke promote faster propagation and shallower arrest. Volatile loss from ascending magma remains limited until the stagnation depth, providing a potential mechanism for transfer of deep volatiles to hypabyssal blind intrusions associated with porphyry ore deposits. Our findings are applicable to the problem of silicic magma ascent through the crust more generally

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

Successive episodes of reactive liquid flow through a layered intrusion (Unit 9, Rum Eastern Layered Intrusion, Scotland)

We present a detailed microstructural and geochemical study of reactive liquid flow in Unit 9 of the Rum Eastern Layered Intrusion, Scotland. Unit 9 comprises an underlying lens-like body of peridotite overlain by a sequence of troctolite and gabbro (termed allivalite), with some local and minor anorthosite. The troctolite is separated from the overlying gabbro by a distinct, sub-horizontal, undulose horizon (the ‘major wavy horizon’). Higher in the stratigraphy is another, similar, horizon (the ‘minor wavy horizon’) that separates relatively clinopyroxene-poor gabbro from an overlying gabbro. To the north of the peridotite lens, both troctolite and gabbro grade into poikilitic gabbro. Clinopyroxene habit in the allivalite varies from thin rims around olivine in troctolite to equigranular crystals in gabbro and to oikocrysts in poikilitic gabbro. The poikilitic gabbros contain multiple generations of clinopyroxene, with Cr-rich (~1.1 wt% Cr2O3) anhedral cores with moderate REE concentrations (core1) overgrown by an anhedral REE-depleted second generation with moderate Cr (~0.7 wt% Cr2O3) (core2). These composite cores are rimmed by Cr-poor (~0.2 wt% Cr2O3) and REE-poor to -moderate clinopyroxene. We interpret these microstructures as a consequence of two separate episodes of partial melting triggered by the intrusion of hot olivine-phyric picrite to form the discontinuous lenses that comprise the Unit 9 peridotite. Loss of clinopyroxene-saturated partial melt from the lower part of the allivalite immediately following the early stages of sill intrusion resulted in the formation of clinopyroxene-poor gabbro. The spatial extent of clinopyroxene loss is marked by the minor wavy horizon. A second partial melting event stripped out almost all clinopyroxene from the lowest allivalite to form a troctolite, with the major wavy horizon marking the extent of melting during this episode. The poikilitic gabbro formed from clinopyroxene-saturated melt moving upwards and laterally through the remobilized cumulate pile and precipitating clinopyroxene en route. This process, called reactive liquid flow, is potentially important in open magma chambers

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