Search

Textural information, such as crystal size distributions (CSDs) or crystal aspect ratios are powerful tools in igneous petrography for interrogating the thermal history of rocks. They facilitate the investigation of crystal nucleation, growth and mixing as well as the cooling rate of the rock. However, they require large volumes of crystal segmentations and measurements that are often obtained with manual methods. Here a deep learning-based computer vision technique, termed instance segmentation, is proposed to automate the pixel-by-pixel detection of each plagioclase crystal in thin section images. Using predictions from a re-trained model the physical properties of the detected crystals, such as size and aspect ratio, can be rapidly generated to provide textural insights. The present segmentations are validated against published results from manual approaches to prove the method's accuracy. The power and efficiency of this automated approach is showcased by analysing an entire sample suite, segmenting over 48,000 crystals in only a matter of days. Widescale use of this method is expected to drive significant developments in the igneous petrography and related fields

There is considerable interest in combining a range of geophysical, geochemical and geomorphic observations with a view to estimating the amplitude, wavelength and depth of mantle thermal anomalies on a global bases. Here, we wish to explore how forward and inverse modelling of major, trace and rare earth elements can be exploited to determine melt fraction as a function of depth for a mantle peridotitic source. Our focus is on an area that includes the Iberian Peninsula where previous work shows that long-wavelength topography is probably generated and maintained by sub-plate thermal anomalies which are manifest by negative shear-wave velocities. Geological and geomorphic studies suggest that this dynamic support is a Neogene phenomenon. 48 newly acquired Neogene basaltic samples from Spain were analyzed and combined with previously published datasets. Both major element thermobarometry and rare earth element inverse modelling are applied to estimate melt fraction as a function of depth. In this way, asthenospheric potential temperature and lithospheric thickness can be gauged. These estimates are compared with those obtained from calibrated shear-wave tomographic models. Our results show that potential temperatures and lithospheric thicknesses are 1250-1300 °C and 65-70 km, respectively. These values broadly agree with calibrated tomographic models which yield values of 1300-1350 °C and 45-70 km. We conclude that a region encompassing Iberia is dynamically supported by a combination of warm asthenosphere and thinned lithosphere. This conclusion broadly agrees with independently obtained residual depth anomalies which indicate that the Western Mediterranean region probably has moderately positive dynamic support.

Textural information, such as crystal size distributions (CSD’s) or crystal aspect ratios are powerful tools in igneous petrography for interrogating the thermal history of rocks and the timescales of processes affecting them [1-3]. Plagioclase feldspar especially has found extensive use as a reliable tracer for igneous thermal history and processes with both the apparent 2D [4] and 3D [5] morphologies shown to vary predictably with crystallization time. However, most textural studies, especially relating to 3D morphologies, require extensive data collection which can be cumbersome and time consuming when performed manually. The aim of this work is to present a holistic and automated workflow to enable the rapid extraction of igneous timescales from plagioclase textures through an automated approach. These developments are vital to better allow petrologists to make timescale estimates that can be used in conjunction with diffusion chronometry and more fully characterise the temperature-time paths of igneous rocks. We propose the use of a deep learning-based computer vision technique, termed instance segmentation [6-7], to automatically detect the exact pixel-by-pixel location of each plagioclase crystal (crystal masks) in thin section images. By re-training the models using a custom set of segmented geological thin section images, one can re-purpose these models for petrographic use, limitations notwithstanding based on the training data. The model outputs can then be used to measure the physical properties of the detected crystals, such as size and aspect ratio, to automate the production of CSDs and aspect ratio distributions which are routinely used to interrogate the timescales of igneous processes. The validity of our method will be showcased using a range of established volcanic and plutonic sample sets that have been previously well-characterised [4,8] through manual segmentation; these will include subglacial pillow basalts from Skuggafjoll and basaltic intrusions such as the Basement Sill in Antarctica and the Karlshamn dyke from Sweden. For sills, we make use of the correlation between plagioclase shape and crystallisation time [4] for rapid timescale determination straight from thin section photomicrographs to complement the information acquired from CSD’s. The vast amounts of data available from the automated segmentation of thin section scans are ripe for 3D shape studies over extensive sample suites to complement traditional textural approaches to timescales. These timescales will be linked to those obtained from diffusion chronometry such as Mg-in-plagioclase diffusion. References: [1] Cashman KV and Marsh BD (1988) Contrib Mineral Petrol 99, 277–291 [2] Higgins MD (2000) American Mineralogist, 85, 1105-1116 [3] Armienti P (2008) Reviews in Mineralogy and Geochemistry. 69. 623-649 [4] Holness MB (2014) Contrib Mineral Petrol 168, 1076 [5] Mangler MF et al. (2022) Contrib Mineral Petrol 177, 64 [6] He K et al. (2017) IEEE International Conference on Computer Vision (ICCV) pp. 2980-2988 [7] Qiao S et al (2021) Proc. IEEE/CVF Conf. CVPR pp. 10213-1022 [8] Neave DA et al. (2014) Crystal Storage and Transfer in Basaltic Systems: The Skuggafjöll Eruption, Iceland, Journal of Petrology, Volume 55 pp.2311–2346

The prevalence, durability and physical significance of crystal mushes in crustal magmatic systems is a topic of current interest in igneous petrology. Fragments of mushes brought to the surface by basaltic eruptions provide a snapshot of the temporal evolution of crustal magmatic systems. Petrographic and geochemical analysis of such fragments give valuable insights into basaltic magma reservoirs, including information about magma storage conditions and possible eruption triggers. A detailed petrological and geochemical study was carried out on gabbroic mush nodules from the Brandur, Fontur and Saxi tuff cones to understand the processes that occur before large fissure eruptions in the Bárðarbunga system, Iceland.Petrographic studies of the mush nodules, from QEMSCAN images, reveal a bimodal phenocryst population in a glassy vesicular groundmass. Probe analyses confirm the bimodal population consists of a primitive and evolved assemblage. The former is composed of large equant crystals of high-anorthite plagioclase (An~88), high-forsterite olivine (Fo~86) and high Mg# clinopyroxene (Mg#~86) forming an interconnected solid framework. The evolved assemblage consists of low-anorthite plagioclase (An~75), low-forsterite olivine (Fo~77) and low Mg# clinopyroxene (Mg#~79) crystallising in the pore space of the mush framework and on the rims of the primitive macrocrysts. The textures and compositions seen suggest the nodules experienced two stages of crystallisation: primitive macrocrysts crystallised first and were stored in crystal mushes. Then a later event caused a change in PTX conditions and triggered relatively rapid crystallisation in the pore-spaces of the mushes.The quenched glass in the pore spaces of the nodules has the composition of a basaltic liquid that in chemical equilibrium with the evolved assemblage of crystals. Thermobarometry based on equilibrium between this liquid and the phases indicates that the final stage of crystallisation occurred at pressures of ~2 kbar. A putative interstitial liquid composition was reconstructed under the assumption of closed system growth of the evolved assemblage by using the QEMSCAN pixel maps to add the evolved crystals to the interstitial glass composition. This reconstructed liquid is far from chemical equilibrium with the primitive crystals in the mush framework, indicating that the assumption of simple closed system crystallisation from an initial mush liquid in equilibrium with the primitive solids is not correct. Therefore, the phase mapping and compositional relationship provide constraints on open-system processes in mushes.The failure of the closed system models to match the observations is significant in two ways. First, the lack of equilibrium between mush liquid and cumulus plagioclase is consistent with the expected sluggish diffusion of NaSi-CaAl in plagioclase. This disequilibrium poses challenges for numerical models of magmatic systems that use the assumption of crystal-melt equilibrium to link temperature, melt fraction and phase compositions. Second, bubble expansion during pre-eruptive ascent forces mush liquid out of solid framework in the nodules and may provide observational constraints on the physics of multiphase flow in deep magmatic systems.

Effective eruption forecasting and volcanic hazard management depend heavily on our ability to detect when a volcanic system switches from a state of unrest into a state of eruption. The 2021 eruption at Fagradalsfjall in SW Iceland, the first deep-sourced eruption on a mid-ocean ridge system monitored with modern instrumentation, presents an ideal opportunity to compare geophysical and petrological datasets to explore processes of deep magma mobilisation and eruption priming. Here we use diffusion chronometry to show that deep magmatic unrest in the roots of volcanic systems can precede apparent geophysical eruption precursors by a few years. Early phases of magma accumulation and reorganisation in the near-Moho plumbing system, part of the priming for eruption, can occur in the absence of significant increases in shallow seismicity (

Supplementary Methods from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Figure 6 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Figure 7 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

S1P2 sphingosine 1-phosphate receptor signaling can regulate proliferation, survival, morphology, and migration in many cell types in vitro. Here, we report that S1P2−/− mice develop clonal B-cell lymphomas with age, such that approximately half of the animals display this neoplasm by 1.5 to 2 years of age. Histologic, immunophenotypic, and molecular analyses revealed a uniform tumor phenotype with features of germinal center (GC)–derived diffuse large B-cell lymphoma (DLBCL). Tumor formation was preceded by increases in GC B cells and CD69+ T cells, as well as an increased formation of spontaneous GCs, suggesting that S1P2 loss may promote lymphomagenesis in part by disrupting GC B-cells homeostasis. With the sole exception of rare lung tumors, the effect of S1P2 gene disruption is remarkably restricted to DLBCL. In humans, 28 of 106 (26%) DLBCL samples were found to harbor multiple somatic mutations in the 5′ sequences of the S1P2 gene. Mutations displayed features resembling those generated by the IgV-associated somatic hypermutation mechanism, but were not detected at significant levels in normal GC B cells, indicating a tumor-associated aberrant function. Collectively, our data suggest that S1P2 signaling may play a critical role in suppressing DLBCL formation in vivo. The high incidence of DLBCL in S1P2−/− mice, its onset at old age, and the relative lack of other neoplasms identify these mice as a novel, and potentially valuable, model for this highly prevalent and aggressive human malignancy. [Cancer Res 2009;69(22):8686–92]

Supplementary Figure 3 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Figure 5 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Figure 1 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Tables 1-3 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Figure 2 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

Supplementary Figure 4 from Targeted Disruption of the S1P2 Sphingosine 1-Phosphate Receptor Gene Leads to Diffuse Large B-Cell Lymphoma Formation

We use a dense seismic network on the Reykjanes Peninsula, Iceland, to image a group of earthquakes at 10–12 km depth, 2 km north-east of 2021 Fagradalsfjall eruption site. These deep earthquakes have a lower frequency content compared to earthquakes located in the upper, brittle crust and are similar to deep long period (DLP) seismicity observed at other volcanoes in Iceland and around the world. We observed several swarms of DLP earthquakes between the start of the study period (June 2020) and the initiation of the 3-week-long dyke intrusion that preceded the eruption in March 2021. During the eruption, DLP earthquake swarms returned 1 km SW of their original location during periods when the discharge rate or fountaining style of the eruption changed. The DLP seismicity is therefore likely to be linked to the magma plumbing system beneath Fagradalsfjall. However, the DLP seismicity occurred ~ 5 km shallower than where petrological modelling places the near-Moho magma storage region in which the Fagradalsfjall lava was stored. We suggest that the DLP seismicity was triggered by the exsolution of CO2-rich fluids or the movement of magma at a barrier to the transport of melt in the lower crust. Increased flux through the magma plumbing system during the eruption likely adds to the complexity of the melt migration process, thus causing further DLP seismicity, despite a contemporaneous magma channel to the surface.

The thermochemical structure of lithospheric and asthenospheric mantle exert primary controls on surface topography and volcanic activity. Volcanic rock compositions and mantle seismic velocities provide indirect observations of this structure. Here, we compile and analyze a global database of the distribution and composition of Neogene-Quaternary intraplate volcanic rocks. By integrating this database with seismic tomographic models, we show that intraplate volcanism is concentrated in regions characterized by slow upper mantle shear-wave velocities and by thin lithosphere (i.e., Not applicable

Chalcophile element concentrations in melt inclusions and matrix glasses may be used to investigate low pressure degassing processes, as well as sulfide saturation during crustal fractionation, and mantle melting. Erupted products from Kīlauea Volcano, Hawaiʻi, record three stages of sulfide saturation (in the mantle, crust, and within lava lakes), separated by episodes of sulfide resorption (i.e., sulfide undersaturation) during ascent through the thick Hawaiian lithosphere, and during syn-eruptive degassing. The presence of residual sulfides in the mantle source throughout the melting interval accounts for the high S concentrations of Kīlauean primary melts (1387–1600 ppm). Residual sulfides retain chalcophile elements during melting, decoupling the variability of these elements in high MgO melts from that of lithophile elements. Decompression associated with magma ascent through the thick Hawaiian lithosphere drives an increase in the sulfide concentration at sulfide saturation (SCSS2-), resulting in shallow storage reservoirs (∼1–5 km depth) being supplied with sulfide-undersaturated melts. A drop in temperature, coupled with major element changes during the fractionation of olivine, causes the SCSS2- to decrease. Combined with an increase in melt S contents during fractionation, this initiates a second stage of sulfide saturation at relatively high MgO contents (∼12 wt% MgO). Syn-eruptive degassing of S drives the resorption of sulfides in contact with the carrier liquid. The covariance structure of Cu, MgO and Ni contents in melt inclusions and matrix glasses indicates that the dissolution of sulfides effectively liberates sulfide-hosted Cu and Ni back into the melt, rather than the vapour phase. The contrasting behaviour of Cu, Ni, Se and S during sulfide resorption indicates that the chalcophile element signature of the Kīlauean plume is largely controlled by silicate melt-vapour partitioning, rather than sulfide-vapour partitioning. The participation of dense sulfide liquids in shallow degassing processes may result from their direct attachment to buoyant vapour bubbles, or olivine crystals which were remobilized prior to eruption. Sulfide resorption obscures the textural and chemical record of sulfide saturation in matrix glasses, but not in melt inclusions, which are isolated from this late-stage release of chalcophile elements. The partitioning of S between the dissolving sulfide, melt and the vapour phase accounts for approximately 20% of the total S release into the atmosphere.

Dynamic topography is the surface expression of sub-plate mantle convective processes. In recent years, there has been considerable interest in combining a wide range of geophysical, geological and geomorphic observations with a view to determining the amplitude, wavelength and depth of mantle thermal anomalies. Here, we are interested in exploring how quantitative modelling of major, trace and rare earth elements can be used to constrain the depth and degree of asthenospheric melting for a mantle peridotitic source. Our focus is on a region that encompasses the Iberian Peninsula where previous research suggests that long-wavelength topography is supported by a significant sub-plate thermal anomaly which is manifest by reduced shear-wave velocities. Stratigraphic and fluvial studies imply that this dynamic support is a Neogene phenomenon. We analyzed 48 Neogene basaltic rocks that were acquired from Iberia in September 2019 and combined these analyses with previously published datasets. Both major element thermobarometry and rare earth element inverse modelling are used to determine the asthenospheric potential temperature and lithospheric thickness. These values are compared with those estimated from calibrated shear-wave tomographic models. Our geochemical results indicate that potential temperatures and lithospheric thicknesses are 1300-1375 °C and 50-80 km, respectively. These values broadly agree with calibrated tomographic models which yield values of 1300-1360 °C and 45-70 km. We conclude that a region encompassing Iberia is dynamically supported by a combination of warm asthenosphere and thinned lithosphere.

Magmas with matrix glass compositions ranging from basalt to dacite erupted from a series of 24 fissures in the first two weeks of the 2018 Lower East Rift Zone (LERZ) eruption of Kīlauea Volcano. Eruption styles ranged from low spattering and fountaining to strombolian activity. Major element trajectories in matrix glasses and melt inclusions hosted by olivine, pyroxene and plagioclase are consistent with variable amounts of fractional crystallization, with incompatible elements (e.g., Cl, F, H2O) becoming enriched by 4-5 times as melt MgO contents evolve from 6 to 0.5 wt%. The high viscosity and high H2O contents (~2 wt%) of the dacitic melts erupting at Fissure 17 account for the explosive Strombolian behavior exhibited by this fissure, in contrast to the low fountaining and spattering observed at fissures erupting basaltic to basaltic-andesite melts. Saturation pressures calculated from melt inclusions CO2-H2O contents indicate that the magma reservoir(s) supplying these fissures was located at ~2-3 km depth, which is in agreement with the depth of a dacitic magma body intercepted during drilling in 2005 (~2.5 km) and a seismically-imaged low Vp/Vs anomaly (~2 km depth). Nb/Y ratios in erupted products are similar to lavas erupted between 1955-1960, indicating that melts were stored and underwent variable amounts of crystallization in the LERZ for >60 years before being remobilized by a dike intrusion in 2018. We demonstrate that extensive fractional crystallization generates viscous and volatile-rich magma with potential for hazardous explosive eruptions, which may be lurking undetected at many ocean island volcanoes.

The Earth’s mantle holds more carbon than its oceans, atmosphere and continents combined, yet the distribution of carbon within the mantle remains uncertain. Our best constraints on the distribution of carbon within the upper mantle are derived from the carbon-trace element systematics of ultra-depleted glasses and melt inclusions from mid-ocean ridge basalts. However, carbon-trace element systematics are susceptible to modification by crustal processes, including concurrent degassing and mixing, and melt inclusion decrepitation. In this study we explore how the influence of these processes varies systematically with both the mantle source and melting process, thereby modulating both global and local carbon-trace element trends. We supplement the existing melt inclusion data from Iceland with four new datasets, significantly enhancing the spatial and geochemical coverage of melt inclusion datasets from the island. Within the combined Iceland dataset there is significant variation in melt inclusion C/Ba ratio, which is tightly correlated with trace element enrichment. The trends in C/Ba-Ba space displayed by our new data coincide with the same trends in data compiled from global ocean islands and mid-ocean ridges, forming a global array. The overall structure of the global C/Ba-Ba array is not a property of the source, instead it is controlled by CO2 vapour loss pre- and post-melt inclusion entrapment; i.e., the array is a consequence of degassing creating near-constant maximum melt-inclusion carbon contents over many orders of magnitude of Ba concentration. On Iceland, extremely high C/Ba ( > 100) and C/Nb ( > 1000) ratios are found in melt inclusions from the most depleted eruptions. The high C/Ba and C/Nb ratios are unlikely to be either analytical artefacts, or to be the product of extreme fractionation of the most incompatible elements during silicate melting. Whilst high C/Ba and C/Nb ratios could be generated by regassing of melt inclusions by CO2 vapour, or by mantle melting occurring in the presence of residual graphite, we suggest the high values most likely derive from an intrinsically high C/Ba and C/Nb mantle component that makes up a small fraction of the Icelandic mantle.

Recent Icelandic rifting events have illuminated the roles of centralized crustal magma reservoirs and lateral magma transport1–4, important characteristics of mid-ocean ridge magmatism1,5. A consequence of such shallow crustal processing of magmas4,5 is the overprinting of signatures that trace the origin, evolution and transport of melts in the uppermost mantle and lowermost crust6,7. Here we present unique insights into processes occurring in this zone from integrated petrologic and geochemical studies of the 2021 Fagradalsfjall eruption on the Reykjanes Peninsula in Iceland. Geochemical analyses of basalts erupted during the first 50 days of the eruption, combined with associated gas emissions, reveal direct sourcing from a near-Moho magma storage zone. Geochemical proxies, which signify different mantle compositions and melting conditions, changed at a rate unparalleled for individual basaltic eruptions globally. Initially, the erupted lava was dominated by melts sourced from the shallowest mantle but over the following three weeks became increasingly dominated by magmas generated at a greater depth. This exceptionally rapid trend in erupted compositions provides an unprecedented temporal record of magma mixing that filters the mantle signal, consistent with processing in near-Moho melt lenses containing 107–108 m3 of basaltic magma. Exposing previously inaccessible parts of this key magma processing zone to near-real-time investigations provides new insights into the timescales and operational mode of basaltic magma systems.

The Hawaiian Island chain in the middle of the Pacific Ocean is a well-studied example of hotspot volcanism caused by an underlying upwelling mantle plume. However, the thermal and compositional na...

In order to reconcile petrological and geophysical observations in the temporal domain, the uncertainties of diffusion timescales need to be rigorously assessed. Here we present a new diffusion chr...

\ud Crystallization temperatures of primitive olivine crystals have been widely used as both a proxy for, or an intermediate step in calculating, mantle temperatures. The olivine-spinel aluminum-exchange thermometer has been applied to samples from mid-ocean ridges and large igneous provinces, yielding considerable variability in olivine crystallization temperatures. We supplement the existing data with new crystallization temperature estimates for Hawaii, between 1282 ± 21 and 1375 ± 19°C. Magmatic temperatures may be linked to mantle temperatures if the thermal changes during melting can be quantified. The magnitude of this temperature change depends on melt fraction, itself controlled by mantle temperature, mantle composition and lithosphere thickness. Both mantle composition and lithosphere thickness vary spatially and temporally, with systematic differences between mid-ocean ridges, ocean islands and large igneous provinces. For crystallization temperatures to provide robust evidence of mantle temperature variability, the controls of lithosphere thickness and mantle lithology on crystallization temperature must be isolated. We develop a multi-lithology melting model for predicting crystallization temperatures of magmas in both intra-plate volcanic provinces and mid-ocean ridges. We find that the high crystallization temperatures seen at mantle plume localities do require high mantle temperatures. In the absence of further constraints on mantle lithology or melt productivity, we cannot robustly infer variable plume temperatures between ocean-islands and large igneous provinces from crystallization temperatures alone; for example, the extremely high crystallization temperatures obtained for the Tortugal Phanerozoic komatiite could derive from mantle of comparable temperature to modern-day Hawaii. This work demonstrates the limit of petrological thermometers when other geodynamic parameters are poorly known.\ud \ud Plain Language Summary\ud The temperature inside the Earth varies a lot. There are many ways of measuring the mantle's temperature in the present-day, but to understand how our planet has changed through time, we need to know how hot its interior was in the past. One of the ways we can estimate mantle temperatures from ancient and modern rocks is from their crystal chemistry. By measuring the aluminum content in crystals of olivine, we can estimate their crystallization temperature. We do this for crystals from Hawaii. To turn the crystallization temperatures into the mantle temperature we need to know the proportions of minerals the mantle is made of. However, we often don't know all of this information. Using a new model of mantle melting we can calculate how uncertain the mantle temperature is when we only have a crystallization temperature. We find that the mantle under Hawaii is 1582 ± 65°C, much hotter than normal mantle, which has a temperature 1364 ± 23°C. We also apply our method to crystallization temperatures from other locations, including ancient volcanic rocks. We find that crystallization temperatures from large igneous provinces, formed by unusually hot mantle, are consistent with their mantle having a similar temperature to Hawaii

The thermochemical structure of lithospheric and asthenospheric mantle exert primary controls on surface topography and volcanic activity. Volcanic rock compositions and mantle seismic velocities provide indirect observations of this structure. Here, we compile and analyze a global database of the distribution and composition of Neogene-Quaternary intraplate volcanic rocks. By integrating this database with seismic tomographic models, we show that intraplate volcanism is concentrated in regions characterized by slow upper mantle shear-wave velocities and by thin lithosphere (i.e. < 100 km). We observe a negative correlation between shear-wave velocities at depths of 125-175 km and melt fractions inferred from volcanic rock compositions. Furthermore, mantle temperature and lithospheric thickness estimates obtained by geochemical modeling broadly agree with values determined from tomographic models that have been converted into temperature. Intraplate volcanism often occurs in regions where uplifted (but undeformed) marine sedimentary rocks are exposed. Regional elevation of these rocks can be generated by a combination of hotter asthenosphere and lithospheric thinning. Therefore, the distribution and composition of intraplate volcanic rocks through geologic time will help to probe past mantle conditions and surface processes.

Understanding accretion and deformation processes at mid-ocean ridges is crucial as they control the resulting oceanic crustal structure, which covers two-thirds of Earth’s surface. The most common tool for observing such dynamic processes within the Earth is seismic anisotropy. Iceland, which is uplifted by a convective mantle plume and has an active spreading ridge system exposed above sea level, offers a unique opportunity for studying this phenomenon. Here we use a high-resolution dataset of Love and Rayleigh wave speeds to constrain the seismic anisotropy in the Icelandic crust. We show that seismic anisotropy in the lower crust is controlled by crystal preferred orientation, providing a direct observation of lower crustal flow. Furthermore, since shear is needed to align the crystals, our results reveal that crustal flow cannot be a simple translation of mass and requires internal deformation. This finding suggests that crustal flow plays an important role in oceanic crustal accretion and deformation where thick, hot oceanic crust is formed, such as at volcanic rifted margins and where there are mantle plume–ridge interactions. The lower oceanic crust beneath Iceland is flowing and internally deforming, according to constraints on seismic anisotropy in the Icelandic crust from an analysis of seismic surface waves.

The 2018 lower East Rift Zone (LERZ) eruption and the accompanying collapse of the summit caldera marked the most destructive episode of activity at Kı̄lauea Volcano in the last 200 years. The eruption was extremely well‐monitored, with extensive real‐time lava sampling as well as continuous geodetic data capturing the caldera collapse. This multi‐parameter dataset provides an exceptional opportunity to determine the reservoir geometry and magma transport paths supplying Kı̄lauea’s LERZ. The forsterite contents of olivine crystals, together with the degree of major element disequilibrium with carrier melts, indicates that two distinct crystal populations were erupted from Fissure 8 (termed High‐ and Low‐Fo). Melt inclusion entrapment pressures reveal that Low‐Fo olivines (close to equilibrium with their carrier melts) crystallized within the Halema’uma’u reservoir ( ∼2 km depth), while many High‐Fo olivines ( > Fo81.5; far from equilibrium with their carrier melts) crystallized within the South Caldera reservoir ( ∼3–5 km depth). Melt inclusions in High‐Fo olivines experienced extensive post‐entrapment crystallization following their incorporation into cooler, more evolved melts. This favoured the growth of a CO2 ‐rich vapor bubble, containing up to 99% of the total melt inclusion CO2 budget (median=93%). If this CO2 ‐rich bubble is not accounted for, entrapment depths are significantly underestimated. Conversely, reconstructions using equation of state methods rather than direct measurements of vapor bubbles overestimate entrapment depths. Overall, we show that direct measurements of melts and vapor bubbles by SIMS and Raman Spectroscopy, combined with a suitable H2 O −CO2 solubility model, is a powerful tool to identify the magma storage reservoirs supplying volcanic eruptions.

Crystallisation temperatures of primitive olivine crystals have been widely used as both a proxy for, or an intermediate step in calculating, mantle temperatures. The olivine-spinel aluminium-exchange thermometer has been applied to many samples from mid-ocean ridges, ocean islands and large igneous provinces, yielding considerable variability in primitive olivine crystallisation temperatures. We supplement the existing data with new crystallisation temperature estimates for Hawaii, in the range 1282±21 - 1375±19°C.Magmatic temperatures may be linked to mantle temperatures if the thermal changes during melting can be quantified. Melting lowers the temperature of co-existing magma and solid mantle, owing to the latent heat of melting. The magnitude of this cooling depends on melt fraction, itself controlled by mantle temperature, mantle lithology and lithosphere thickness. All of these parameters are likely to vary both spatially and temporally. For robust quantification of mantle temperature variability, the controls of lithosphere thickness and mantle lithology on crystallisation temperatures must be isolated.We develop a multi-lithology melting model that can predict crystallisation temperature. The model allows mantle temperature, lithospheric thickness, and fractions of mantle lherzolite, pyroxenite and harzburgite to be varied. Inverting the model using a Bayesian Monte Carlo routine enables assessment of the extent to which crystallisation temperatures require variations in mantle temperature. We find that the high crystallisation temperatures seen at mantle plume localities do require high mantle temperatures. However, in the absence of further constraints on mantle lithology or melt productivity, we cannot robustly infer variable plume temperatures in either the present-day or throughout the phanerozoic. This work demonstrates the limit of petrological thermometers when other geodynamic parameters are poorly known.

Kīlauea Volcano emits large quantities of sulfur and other chalcophile elements into the troposphere as gas and aerosol particles, with widespread implications for regional air quality. The concentration of these elements in erupting melts is controlled by sulfide saturation, as well as the interplay between the exsolved volatile phase, silicate and sulfide liquids on eruption. Analysis of sulfur and other chalcophile elements (e.g, Ni, Cu, Se, As, Bi, Cd) in melt inclusions and matrix glasses allow deconvolution of these various processes. Olivine-hosted melt inclusions have significantly lower Ni and Cu concentrations than matrix glasses, defining trajectories consistent with sulfide saturation. These observations, when interpreted with the latest generation of sulfide saturation models, demonstrate that sulfides saturate at high MgO contents (10-14 wt%), in contrast to the traditional interpretation that sulfide saturation occurs relatively late at Kīlauea (~ 2 wt% MgO). This apparent discrepancy may be reconciled by considering the behaviour of sulfides during syn-eruptive degassing. The release of ~90% of dissolved sulfur into the vapour phase at low pressures leads to previously sulfide-saturated magmas becoming sulfide-undersaturated, driving the resorption of sulfides in contact with the degassing silicate melt. Sulfide resorption releases Cu, Ni, S and other chalcophile elements into the vapour-melt-(sulfide) system. Comparisons of melt inclusion and matrix glasses reveals that significant quantities of S, Se, Bi and As partition into the exsolved volatile phase. Other elements, such as Ni and Cu, remain largely in the melt. The contrasting behaviour of Se and Cu demonstrates that chalcophile element degassing is largely controlled by fluid-melt, rather than sulfide-melt partitioning. Crucially, sulfide resorption obscures the textural and chemical record of sulfide saturation in matrix glass and whole-rocks, but not in melt inclusions, which are isolated from the late-stage release of chalcophile elements from sulfide breakdown. Sulfide resorption during degassing and eruption provides a significant, but previously unquantified flux of sulfur to the atmosphere (as SO2) at Kīlauea. Careful evaluation of melt inclusion Cu-Ni-S systematics reveals that the total S release during eruptions is ~1450 ppm (1.8 x previous estimates).

Distorted olivines of enigmatic origin are ubiquitous in erupted products from a wide range of volcanic systems (e.g., Hawai'i, Iceland, Andean Southern Volcanic Zone). At Kīlauea volcano, distorted olivines are commonly attributed to ductile creep within dunitic bodies located around the central conduit, or within the deep rift zones (~5–9 km depth). However, a recent suggestion that lattice distortions are produced by an early phase of branching dendritic growth, followed by textural ripening and the merging of misoriented crystal buds, has gained considerable traction.A quantitative examination of the microstructures in distorted olivines by electron backscatter diffraction (EBSD) reveals striking similarities to microstructures observed in deformed mantle peridotites, but significant differences to the crystallographic signatures of dendritic growth. This suggests that lattice distortions record the application of differential stresses at high temperatures within the magmatic plumbing system, rather than rapid crystal growth. Previous petrological work has suggested that differential stresses are produced by ductile creep within Kīlauea’s deep rift zones. Crucially, this has fuelled suggestions that significant quantities of magma must travel along these rift zones in order to acquire distorted olivines, despite the paucity of geophysical evidence for these magma transport paths. In contrast, we show that the spatial distribution of eruptions containing distorted olivines is consistent with their derivation from the main magma storage reservoir. This model not only aligns petrological and geophysical observations at Kīlauea, but also accounts for the occurrence of distorted olivines in a wide variety of basaltic systems worldwide (which lack deep rift zones).Application of piezometers developed for mantle peridotites reveals that distorted olivines have experienced differential stresses of ~3–12 MPa. Assuming that mush piles behave as granular materials, and form force chains, these stresses can be generated within cumulate piles of ~180–720 m. Based on available constraints on the magma supply rate and the geometry of Kīlauea’s summit reservoir, these thicknesses accumulate in a few centuries (consistent with residence times inferred from melt inclusion records).Overall, we demonstrate that microstructural investigations of erupted olivine crystals by EBSD generates rich datasets which provide quantitative insights into crystal storage within mush piles. Under the increasingly prevalent view that crustal magmatic systems are mush-dominated, constraining the geometry and dynamics of crystal storage regions is crucial to further our understanding of magmatic plumbing systems. The presence of distorted olivines in many different volcanic settings highlights the global applicability of the methods developed in this study. Furthermore, assessments of deformation conditions using EBSD need not be restricted to olivine-bearing lavas. Microstructural fabrics types in natural and experimental samples have been established for a wide variety of igneous phases (e.g. diopside, plagioclase, hornblende), so similar approaches may be utilized in more evolved volcanic systems.

©2018. American Geophysical Union. All Rights Reserved. Western North America has an average elevation that is ∼2 km higher than cratonic North America. This difference coincides with a westward decrease in average lithospheric thickness from ∼240 to 260 basaltic samples. Forward and inverse modeling of carefully selected major, trace, and rare earth elements were used to determine melt fraction as a function of depth. Basaltic melt appears to have been generated by adiabatic decompression of dry peridotite with asthenospheric potential temperatures of 1340 ± 20 °C. Potential temperatures as high as 1365 °C were obtained for the Snake River Plain. For the youngest (i.e.

In the North Atlantic Ocean, the geometry of diachronous V-shaped features that straddle the Reykjanes Ridge is often attributed to thermal pulses which advect away from the center of the Iceland plume. Recently, two alternative hypotheses have been proposed: rift propagation and buoyant mantle upwelling. Here, we evaluate these different proposals using basin-wide geophysical and geochemical observations. The centerpiece of our analysis is a pair of seismic reflection profiles oriented parallel to flowlines that span the North Atlantic Ocean. V-shaped ridges and troughs are mapped on both Neogene and Paleogene oceanic crust, enabling a detailed chronology of activity to be established for the last 50 million years. Estimates of the cumulative horizontal displacement across normal faults help to discriminate between brittle and magmatic modes of plate separation, suggesting that crustal architecture is sensitive to the changing planform of the plume. Water-loaded residual depth measurements are used to estimate crustal thickness and to infer mantle potential temperature which varies by 25◦C on timescales of 3–8 Ma. This variation is consistent with the range of temperatures inferred from geochemical modeling of dredged basaltic rocks along the ridge axis itself, from changes in Neogene deep-water circulation, and from the regional record of episodic Cenozoic magmatism. We conclude that radial propagation of transient thermal anomalies within an asthenospheric channel that is 150 50 km thick best accounts for the available geophysical and geochemical observations.

Mantle oxygen fugacity (fO2) governs the physico-chemical evolution of the Earth, however current estimates from commonly used basalt redox proxies are often in disagreement. In this study we compare three different potential basalt fO2 proxies: Fe3+/Fetot, V/Sc and V isotopes, determined on the same submarine lavas from a 700 km section of the Reykjanes Ridge, near Iceland. These samples provide a valuable test of the sensitivities of fO2 proxies to basalt petrogenesis, as they formed at different melting conditions and from a mantle that towards Iceland exhibits increasing long-term enrichment of incompatible elements. New trace element data were determined for 63 basalts with known Fe3+/Fetot. A subset of 19 lavas, covering the geographical spread of the ridge transect, was selected for vanadium isotope analyses. Vanadium is a multi-valence element whose isotopic fractionation is theoretically susceptible to redox conditions. Yet, the δ 51 VAA composition of basaltic glasses along the Reykjanes Ridge covers only a narrow range ( δ 51 VAA = −1.09 to −0.86‰; 1SD = 0.02–0.09) and does not co-vary with fractionation-corrected Fe3+/Fetot (0.134–0.151; 1SD = 0.005) or V/Sc (6.6–8.5; 1SD = 0.1-1.3) ratios. However, on a global scale, basaltic δ 51 VAA may be controlled by the extent of melting. The V/Sc compositions of primitive (MgO > 7.5 wt%) basalts show no systematic change along the entire length of the Reykjanes Ridge. Typical peridotite melting models in which source Fe3+/Fetot is constant at 5% and that account for the increased mantle potential temperature nearer the plume center and the fO2 dependent partitioning of V, can reproduce the V/Sc data. However, while these melting models predict that basalt Fe3+/Fetot ratios should decrease with increasing mantle potential temperature towards Iceland, fractionation-corrected Fe3+/Fetot of Reykjanes Ridge lavas remain nearly constant over the ridge length. This discrepancy is explained by source heterogeneity, where an oxidized mantle pyroxenite component contributes to melting with increasing proximity to Iceland. Comparison of observed and modeled Fe3+/Fetot indicate that source variation in fO2 is present under the Reykjanes Ridge, with higher Fe3+/Fetot closer to Iceland. This source variability in fO2 cannot be resolved by V isotopes and redox-sensitive trace element ratios, which instead appear to record magmatic processes.