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

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

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

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

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

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°