The Genesis of Intermediate and Silicic Magmas in Deep Crustal Hot Zones

A model for the generation of intermediate and silicic igneous rocks is presented, based on experimental data and numerical modelling. The model is directed at subduction-related magmatism, but has general applicability to magmas generated in other plate tectonic settings, including continental rift zones. In the model mantle-derived hydrous basalts emplaced as a succession of sills into the lower crust generate a deep crustal hot zone. Numerical modelling of the hot zone shows that melts are generated from two distinct sources; partial crystallization of basalt sills to produce residual H2O-rich melts; and partial melting of pre-existing crustal rocks. Incubation times between the injection of the first sill and generation of residual melts from basalt crystallization are controlled by the initial geotherm, the magma input rate and the emplacement depth. After this incubation period, the melt fraction and composition of residual melts are controlled by the temperature of the crust into which the basalt is intruded. Heat and H2O transfer from the crystallizing basalt promote partial melting of the surrounding crust, which can include meta-sedimentary and meta-igneous basement rocks and earlier basalt intrusions. Mixing of residual and crustal partial melts leads to diversity in isotope and trace element chemistry. Hot zone melts are H2O-rich. Consequently, they have low viscosity and density, and can readily detach from their source and ascend rapidly. In the case of adiabatic ascent the magma attains a super-liquidus state, because of the relative slopes of the adiabat and the liquidus. This leads to resorption of any entrained crystals or country rock xenoliths. Crystallization begins only when the ascending magma intersects its H2O-saturated liquidus at shallow depths. Decompression and degassing are the driving forces behind crystallization, which takes place at shallow depth on timescales of decades or less. Degassing and crystallization at shallow depth lead to large increases in viscosity and stalling of the magma to form volcano-feeding magma chambers and shallow plutons. It is proposed that chemical diversity in arc magmas is largely acquired in the lower crust, whereas textural diversity is related to shallow-level crystallizatio

Plagioclase‐Saturated Melt Hygrothermobarometry and Plagioclase‐Melt Equilibria Using Machine Learning

Compositions of plagioclase‐melt pairs are commonly used to constrain temperatures (T), dissolved water contents (H2O) and pressures (P) of pre‐eruptive magma storage and transport. However, previous plagioclase‐based thermometers, hygrometers, and barometers can have significant errors, leading to imprecise reconstructions of conditions during plagioclase growth. Here, we explore whether we can refine existing plagioclase‐based hygrothermobarometers with either plagioclase‐melt or melt‐only chemistry (±T/H2O), calibrated using random forest machine learning on experimental petrology data (n = 1,152). We find that both the plagioclase‐melt and melt‐only models return similar cross‐validation root‐mean‐square errors (RMSEs), as the melt holds most of the P‐T‐H2O information rather than the plagioclase. T/H2O‐dependent melt models have test set RMSEs of 25°C, 0.70 wt.% and 76 MPa for temperature, H2O content and pressure, respectively, while T/H2O‐independent models have RMSEs of 38°C, 0.97 wt.% and 91 MPa. The melt thermometer and hygrometer are applicable to a wide range of plagioclase‐bearing melts at temperatures between 664 and 1355°C, and with H2O concentrations up to 11.2 wt.%, while the melt barometer is suitable for pressures of ≤500 MPa. An updated plagioclase‐melt equilibrium model has also been calibrated, allowing the equilibrium anorthite content to be predicted with an error of 5.8 mol%. The new P‐T‐H2O‐An models were applied to matrix glasses and melt inclusions from the 1980 Mount St Helens (USA) and 2014–2015 Holuhraun (Iceland) eruptions, corroborating previous independent estimates and observations. Models are available at https://github.com/kyra‐cutler/Plag‐saturated‐melt‐P‐T‐H2O‐An, enabling assessment of plagioclase‐melt equilibrium and characterization of last‐equilibrated P‐T‐H2O conditions of plagioclase‐saturated magmas

The sources of granitic melt in Deep Hot Zones

A Deep Hot Zone develops when numerous mafic sills are repeatedly injected at Moho depth or are scattered in the lower crust. The melt generation is numerically modelled for mafic sill emplacement geometries by overaccretion, underaccretion or random emplacement, and for intrusion rates of 2, 5 and 10 mm/yr. After an incubation period, melts are generated by incomplete crystallisation of the mafic magma and by partial melting of the crust. The first melts generated are residual from the mafic magmas that have low solidi due to concentration of H2O in the residual liquids. Once the solidus of the crust is reached, the ratio of crustal partial melt to residual melt increases to a maximum. If wet mafic magma, typical of arc environments, is injected in an amphibolitic crust, the residual melt is dominant over the partial melt, which implies that the generation of I-type granites is dominated by the crystallisation of mafic magma originated from the mantle and not by the partial melting of earlier underplated material. High ratios of crustal partial melt over residual melt are reached when sills are scattered in a metasedimentary crust, allowing the generation of S-type granites. The partial melting of a refractory granulitic crust intruded by dry, high-T mafic magma is limited and subordinate to the production of larger amount of residual melt in the mafic sills. Thus the generation of A-type granites by partial melting of a refractory crust would require a mechanism of selective extraction of the A-type mel

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

An Experimental Study of Trace Element Fluxes from Subducted Oceanic Crust

We have determined experimentally the hydrous phase relations and trace element partitioning behaviour of ocean floor basalt protoliths at pressures and temperatures (3 GPa, 750-1000°C) relevant to melting in subduction zones. To avoid potential complexities associated with trace element doping of starting materials we have used natural, pristine mid-ocean ridge basalt (MORB from Kolbeinsey Ridge) and altered oceanic crust (AOC from Deep Sea Drilling Project leg 46, ∼20°N Atlantic). Approximately 15 wt % water was added to starting materials to simulate fluid fluxing from dehydrating serpentinite underlying the oceanic crust. The vapour-saturated solidus is sensitive to basalt K2O content, decreasing from 825 ± 25°C in MORB (∼0·04 wt % K2O) to ≈750°C in AOC (∼0·25 wt % K2O). Textural evidence indicates that near-solidus fluids are sub-critical in nature. The residual solid assemblage in both MORB and AOC experiments is dominated by garnet and clinopyroxene, with accessory kyanite, epidote, Fe-Ti oxide and rutile (plus quartz-coesite, phengite and apatite below the solidus). Trace element analyses of quenched silica-rich melts show a strong temperature dependence of key trace elements. In contrast to the trace element-doped starting materials of previous studies, we do not observe residual allanite. Instead, abundant residual epidote provides the host for thorium and light rare earth elements (LREE), preventing LREE from being released (ΣLREE 1500 and La/SmPUM (where PUM indicates primitive upper mantle) ∼1, most closely matching the geochemical signal of arc lavas worldwide, were generated from AOC at 800-850°

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

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