Magma Flux Research Articles

Spreading ridge migration enabled by plume-ridge de-anchoring

It has long been recognised that spreading ridges are kept in place by competing subduction forces that drive plate motions. Asymmetric strain rates pull spreading ridges in the direction of the strongest slab pull force, which partially explains why spreading ridges can migrate vast distances. However, the interaction between mantle plumes and spreading ridges plays a relatively unknown role on the evolution of plate boundaries. Using a numerical model of mantle convection, we show that plumes with high buoyancy flux (>3000 kg/s) can capture spreading ridges within a 1000 km radius and anchor them in place. Exceptionally high buoyancy fluxes may fragment the overriding plate into smaller plates to accommodate more efficient plate motion. If the plume buoyancy flux wanes below 1000 kg/s the ridge may be de-anchored, leading to rapid ridge migration rates when combined with asymmetric plate boundary forces. Our results show that plume-ridge de-anchoring may have contributed to the rapid migration of the SE Indian Ridge from 43 million years ago (Ma) due to waning buoyancy flux from the Kerguelen plume, supported by magma flux estimates and radiogenic isotope geochemistry of eruption products. The plume-ridge de-anchoring mechanism we have identified has global implications for the evolution of plate boundaries near mantle plumes.

Petrogenesis and Ni–Cu–(PGE) prospectivity of the Mount Ayliff Complex in the Karoo Igneous Province: new insights from the Ingeli and Horseshoe lobes

The Mount Ayliff Complex comprises five cognate mafic–ultramafic bodies emplaced along the southern margin of the Kaapvaal Craton. This study examines the Ingeli and Horseshoe lobes and assesses their magmatic sulfide prospectivity with respect to the Insizwa lobe, which hosts massive sulfides at its Waterfall Gorge occurrence. The Ingeli lobe consists of 465 m of olivine–chromite cumulates overlain by 340 m of (olivine-)gabbronorite, while the Horseshoe lobe consists of 45 m of olivine gabbronorite overlain by 5 m of gabbro. The Ingeli lobe possesses sparsely disseminated sulfides at its mafic–ultramafic transition, whereas disseminated sulfides are present throughout the Horseshoe lobe. A petrogenetic model is proposed, where magma accumulated and fractionated nickeliferous olivine and chromite in an upper-crustal staging chamber hosted by Proterozoic basement rocks. Magma then ascended and deposited olivine–chromite cumulates at the level of the complex. Prolonged magma flux in the staging chamber facilitated the assimilation of basement rocks, triggering sulfide saturation and the crystallization of Ni-poor olivine. Contaminated magma then ascended and interacted with pre-existing cumulates, depositing sulfide melt that may have backflowed as magmatic activity waned. Basaltic magma then flowed over the ultramafic cumulates, depositing disseminated sulfides whilst undergoing closed-system fractionation. Basal depressions and underlying feeder structures are the most prospective locations for magmatic sulfide mineralization.

Co-eruptive, endogenous edifice growth, uplift during 4 years of eruption at Sangay Volcano, Ecuador

We report sustained uplift throughout Volcan Sangay's most recent period of eruption (2019–22), moderated only by transient excursions during some of its largest explosions. Volcan Sangay (Amazonia, Ecuador), has been erupting since 2019, impacting both local communities and distant cities with ash fall and lahars. We analyzed ascending and descending Sentinel-1 radar imagery, constructing a robust network of interferograms spanning this eruptive period to measure relative ground displacements across the volcano. Our time series reveals a consistent uplift pattern (∼68 mm/yr) on the western and northern flanks of the volcano, which we attribute to volume increases in a body of magma located within the volcano's edifice beneath its western flank. This source appears to be vertically extensive, and is best fit by a quadrangular magma pathway, dipping towards the west and increasing in volume by 1.1 × 106 m3 between 2019 and 2022. We additionally identify non-magmatic deformation, including subsidence of fresh deposits and downslope displacement (∼50 mm/year) in the southeastern sector of the volcano. Co-eruptive uplift at Sangay is a rare observation of endogenous growth during an eruption and indicates that stratovolcano edifice stability is sensitive to both magma flux into the edifice and shallow controls on eruption rate.

Reinterpretation of the post-26 ka Taupō Rhyolitic Magmatic System (New Zealand) as Deep and Vertically Extensive Based on Isotope Thermometry and Measured and Modeled Zircon Destinies

Abstract Taupō volcanic zone, the site of the 26 ka Oruanui supereruption, produced ~70 km3 of new rhyolites since 11 ka, culminating in 50 km3 Taupō eruption 1.8 ka. Major phenocrysts decrease from 4 to 1 vol%, and Oruanui and post-Oruanui ignimbrites all have identical high-δ18Omelt values of 7.39 ± 0.1‰ and lack low-δ18O values despite overlapping calderas. The Δ’17O values are −0.07‰, lower than the mantle and indicate source contamination of high-δ18O, low-Δ’17O metasediments, and limited interaction with high-Δ’17O hydrothermally altered crust. Previously published U-Th-Pb zircon ages demonstrate their diversity spanning 104–105 years for each unit. Zircon crystal size distribution shows a decrease in abundance and the mean size, and some units lack small (<~10 um) zircons suggesting that zircons were both growing and dissolving in the coexisting magma generation areas. Isotope thermometry indicates heating of the system from ~812 ± 35°C to 874 ± 36°C past zircon saturation in 1.8 ka eruption. We advocate that a deep vertically continuous and laterally discontinuous silicic magma system at the base of the Taupō rift, rather than a shallow batholith or an evolving mush, drives volcanism at Taupō. To explain the post-Oruanui magma production, rift-base silicic magma origin and moderate (~2 km3/1000 years) rhyodacitic magma flux from a growing and heating liquid magma body creates a sufficient solution for the most recent magmatism.

Magma flux variations triggering shallow-level emplacement of the Takidani pluton (Japan): Insights into the volcanic-plutonic connection

Defining the physico-chemical evolution of felsic magmatic reservoirs in the middle-upper continental crust is crucial for understanding the connection between granitic plutons and silicic volcanic rocks and for research on the trigger mechanisms of volcanic eruptions. In this study, we investigate the formation and evolution of the Hotaka-Takidani volcano-plutonic complex (Japan) using a comprehensive approach that integrates zircon petrochronology with phase-equilibria and thermal modelling. The Takidani pluton is a Pleistocene tilted zoned intrusion of granodiorite-granite composition associated with at least two large caldera-forming eruptions. We have investigated the three uppermost units of the pluton and the oldest and more voluminous Nyukawa eruption, which deposited more than 400 km3 (DRE, dense rock equivalent) of crystal-rich dacitic ignimbrite. Amphibole thermobarometry, U-Pb zircon dates and geochemical data indicate that the Takidani pluton and the Nyukawa dacite were sequentially sourced from a 10–12 km-thick long-lived magmatic reservoir in the middle crust. Thermal evolution of the reservoir led to the extraction of the Nyukawa magma, followed by brief storage in the upper crust and eruption at ca. 1.706 ± 0.091 Ma. For two of the units of the pluton, we combine high-spatial resolution zircon U-Pb geochronology with zircon chemical data. Applying the Ti-in-zircon thermometer, we observe that magma within the reservoir was stored at a relatively constant temperature of ca. 700 °C for 300 kyr and subsequently, starting from 1.4 Ma, the temperature increased progressively reaching 820–850 °C at 1.1 Ma. This temperature shift is inconsistent with a constant magma flux and suggests an increase of the magma injection rate into the reservoir. Using the distribution of zircon ages and thermal simulations, we determined that the reservoir was built by a fairly low average magma flux of ca. 6 × 10−6 km3*km−2*yr−1, shifting to a magma supply rate which is three to four times higher during the last 300 ka of zircon crystallisation. This increase in magma flux caused the final emplacement of the Takidani magma in the upper crust, where it rapidly crystallized. The integration of geochronological data with phase-equilibria and thermo-chemical modelling demonstrates that multiple melt extraction events occurred during the evolution of the Hotaka-Takidani magmatic system. Only the first one triggered a major eruption, while subsequent events resulted in the emplacement of the Takidani pluton in the shallow crust.

Thermal Modeling of the Sanbagawa and Ryoke Belts

The Sanbagawa and Ryoke belts were formed in a convergent plate boundary along the eastern margin of Eurasia. Thermal modeling using the geological records of these belts as constraints allows quantitative estimates of both shear heating along the Wadati-Benioff zone and magma fluxes beneath the volcanic arc. In contrast to real-time observations of crustal movement and heat flow, rocks record changes in pressures and temperatures that occur over periods of several million years and can be used to examine conditions from the surface to the mantle. Thermal modeling combined with such geological records helps to bridge the gap in our knowledge between real-time observations of ongoing geological processes and the development of orogenies in convergent plate margins over geological time.

Influence of magma flux on magma storage depths along the Reykjanes Ridge

Quantitative thermomechanical models of mid-ocean ridge magmatic systems focus on spreading rate variations as the principal control on crustal thermal structure and the distribution of melt storage. The Reykjanes Ridge (RR), a slow-spreading portion of the Mid-Atlantic Ridge, has a near-constant spreading rate, but crustal thickness varies by a factor of two. Therefore, magma flux to the ridge varies independently of spreading rate. Here, we investigate the role of magma flux on storage depths along mid-ocean ridges (MORs) and constrain the roles of mantle melt generation rate and spreading rate on oceanic crustal magmatic systems.Using pyOPAM, an open-source Python script written to carry out olivine–plagioclase– augite–melt (OPAM) thermobarometry, with an uncertainty of 1.2 kbar and 9.7°C, samples between 64°N and 52.5°N are compared. The distribution of median OPAM pressure and temperature estimates are compared with calculated magma flux rates across the area of interest.The section of the RR immediately south of the subaerial Reykjanes Peninsula has a median OPAM pressure of 3.7 kbar, and estimated magma flux of 8.82 × 10−6 km3 km−2 y−1. The southern end of the RR has a median storage pressure of 4.8 kbar, and magma flux of 4.58 × 10−6 km3 km−2 y−1. In comparison, The Charlie-Gibbs Fracture Zone to the south has a median OPAM estimate of 6.3 kbar and calculated magma flux of 2.71 × 10−6 km3 km−2 y−1. Despite important local variation, median final magma storage depths and temperatures before eruption increase southwards along the RR. This finding indicates that magma storage depths at mid-ocean ridges are controlled predominantly by magma flux rather than spreading rate alone. Increasing magma flux drives the shallowing of cooler, final magma storage depths, and decreasing magma flux increases magma storage depths and temperature.

Insights on the state of stress in the mantle beneath Pahala, Hawai‘i

Magma supply rates from the mantle to Hawaiian volcanoes serve as an important control on eruptive behavior at the surface. The Pa ̄hala Sill Complex, a collection of magma-bearing, seismogenic structures at 40 km depth beneath Hawai‘i, presents an opportunity to elucidate interactions between stress and magma transport processes in the mantle. We invert for full moment tensors of sill earthquakes and identify predominantly shear mechanisms with persistent tensile faulting components. Slip occurs in-plane with the sill structures. Pressure axes are radially oriented about a point near Mauna Loa, consistent with a stress field generated by a flexural load. Together, these observations suggest that magma flux through the sill structures generates seismicity by increasing pore pressure and promoting slip. Our results suggest that stress changes in mantle structures may enable fluctuations in magma supply rates to the surface over short timescales.

High-precision U-Pb geochronology links magmatism in the Southwestern Laurentia large igneous province and Midcontinent Rift

Abstract The Southwestern Laurentia large igneous province (SWLLIP) comprises voluminous, widespread ca 1.1 Ga magmatism in the southwestern United States and northern Mexico. The timing and tempo of SWLLIP magmatism and its relationship to other late Mesoproterozoic igneous provinces have been unclear due to difficulties in dating mafic rocks at high precision. New precise U-Pb zircon dates for comagmatic felsic segregations within mafic rocks reveal distinct magmatic episodes at ca. 1098 Ma (represented by massive sills in Death Valley, California, the Grand Canyon, and central Arizona) and ca. 1083 Ma (represented by the Cardenas Basalts in the Grand Canyon and a sill in the Dead Mountains, California). The ca. 1098 Ma magmatic pulse was short-lived, lasting 0.25−0.24+0.67 m.y., and voluminous and widespread, evidenced by the ≥100 m sills in Death Valley, the Grand Canyon, and central Arizona, consistent with decompression melting of an upwelling mantle plume. The ca. 1083 Ma magmatism may have been generated by a secondary plume pulse or post-plume lithosphere extension. The ca. 1098 Ma pulse of magmatism in southwestern Laurentia occurred ~2 m.y. prior to an anomalous renewal of voluminous melt generation in the Midcontinent Rift of central Laurentia that is recorded by the ca. 1096 Ma Duluth Complex layered mafic intrusions. Rates of lateral plume spread predicted by mantle plume lubrication theory support a model where a plume derived from the deep mantle impinged near southwestern Laurentia, then spread to thinned Midcontinent Rift lithosphere over ~2 m.y. to elevate mantle temperatures and generate melt. This geodynamic hypothesis reconciles the close temporal relationships between voluminous magmatism across Laurentia and provides an explanation for that anomalous renewal of high magmatic flux within the protracted magmatic history of the Midcontinent Rift.

The magmatic origin of the Columbia River Gorge, USA.

Along subduction zones, high-relief topography is associated with sustained volcanism parallel to the plate margin. However, the relationship between magmatism and mountain building in arcs is poorly understood. Here, we study patterns of surface deformation and correlated fluvial knickpoints in the Columbia River Gorge to link long-term magmatism to the uplift and ensuing topographic development of the Cascade Range. An upwarped paleochannel exposed in the walls of the Gorge constrains unsteady deep magma flux, the ratio of intrusive to extrusive magmatic contributions to topography, and the impact of magmatism on Columbia River incision since 3.5 million years ago. Geophysical data indicate that deep magma influx beneath the arc axis is ongoing and not aligned with the current locations of volcanic edifices, representing a broad regional influence on arc construction.

Lithospheric Structure Controls Sequentially Active Detachment Faulting at the Longqi Segment on the Southwest Indian Ridge

AbstractOceanic detachment faulting, a major mode of seafloor accretion at slow and ultraslow spreading ridges, is thought to occur during magma‐poor phases and be abandoned when magmatism increases. In this framework, detachment faulting is the result of temporal variations in magma flux, which is inconsistent with recent geophysical observations at the Longqi segment on the Southwest Indian Ridge (49°42′E). In this paper, we focus on this sequentially active detachment faulting system that includes an old, inactive detachment fault and a younger, active detachment fault. We investigate the mechanisms controlling the temporal evolution of this tectonomagmatic system by using 2D mid‐ocean ridge spreading models that simulate faulting and magma intrusion into a visco‐elasto‐plastic continuum. Our models show that temporal variations in magma flux alone are insufficient to match the inferred temporal evolution of the sequentially active detachment system. Rather we find that sequentially active detachment faulting spontaneously occurs at the Longqi segment as a function of lithospheric thickness. This finding is in agreement with an analytical model, which shows that a thicker axial lithosphere results in a smaller fault heave and that a flatter angle in lithosphere thickening away from the accretion axis stabilizes the active fault. A thicker axial lithosphere and its flatter off‐axis angle combined have the potential to modulate sequentially active detachment faulting at the Longqi segment. Our results thus suggest that temporal changes of magmatism are not necessary for the development and abandonment of detachment faults at ultraslow spreading ridges.

Decoding the Link Between Magmatic Cyclicity and Episodic Variation of Tectonics and Crustal Thickness in the Overriding Plate

AbstractCyclical change in subduction angle is the favorable mechanism to elucidate the cyclicity of continental arc magmatism, however, the role of episodic tectonics and variation of the lithosphere in overriding plates is much underestimated. Here we focus on structural, magnetic, and gravitational features of the Late Jurassic to Early Cretaceous granites in the Mesozoic Paleo‐Pacific arc system of the North China block. By unraveling the emplacement process and regional tectonics, we establish a three‐staged extension‐contraction cycle with crustal thickness variation controlling the magmatic flux and behavior. The Late Jurassic extension produced high‐flux crustal‐derived magma (1.87 × 103 km2/Myr), but the thick crust >45 km accumulated large granitic batholiths by multi‐feeders emplacement at the middle‐lower crust and prevented magma ascent and eruption. Subsequently, the Latest Jurassic to Earliest Cretaceous contraction resulted in the magmatic lull and thickened crust of ca. 60 km, fueling crustal material for the ensuing magmatism. In the Early Cretaceous, intense crustal extension thinned the crust to 30 km and largely enhanced the magmatic flux (3.03 × 103 km2/Myr). The magma is prone to penetrate the thin crust with an intensive eruption. A small amount of magma was stored, and the emplacement was controlled by ductile detachments or normal faults. Our model emphasizes episodic the deformation of lithosphere and associated crustal thickness variation in controlling magma production, which may shed new light in understanding the magmatic cyclicity under continuous subduction.

The Ferrar Continental Flood Basalt: A ∼1.6 Ma long duration evidenced by high-precision 40Ar/39Ar ages suggest a potential role in the Pliensbachian-Toarcian extinction event

One of the most impassioned topics in large igneous province (LIP) research is how prolonged the duration of these large-scale magmatic events are, as LIP magmatism has considerable impact on models of associated reconstructions, of climate variability or tectonic events. High-precision geochronology is pivotal to LIP basalt emplacement rate, and thus to unravel the role these enormous magmatic events have throughout Earths geological and environmental history. Four high-precision 40Ar/39Ar plagioclase plateau ages for the Tasmanian dolerites (Ferrar) indicate ∼1.6 ± 0.4 Ma of resolvable, continuous magmatic activity; 184.27 ± 0.24 to 182.69 ± 0.54 Ma (2σ). The 40Ar/39Ar results provide evidence of distinctly older intrusions and a more prolonged duration than the observed 182.4-182.9 Ma age range and duration indicated by the main zircon record. Moreover, the precision of our 40Ar/39Ar results coupled with secondary electron microscopy analyses provide evidence of plagioclase crystal inheritance from slightly older magmatism entrained into younger magmatic pulses by exploiting pre-existing conduits. Numerical diffusion models, calculated for a theoretical age spectrum resulting from two slightly different plagioclase ages, provide an excellent match for measured data. Coupling geochemical data to the new age data indicates a silica and incompatible element evolution of the Ferrar magmatic system through time. The older generation of intrusions (ca. Zr: 92 ppm, SiO2: 53.67 wt.%) are seemingly less enriched in incompatible elements and silica than the youngest generation (ca. Zr: 147 ppm, SiO2: 56.5 wt.%). Here, we suggest that the magma chambers differentiated to more incompatible/silica-rich compositions saturating zircon only at evolved magmatic stages. This implies that plagioclase dates the full duration of magmatic Ferrar LIP activity of ca. 1.6 Myrs whilst zircon ages might be naturally biased and restricted to post-Zr saturation stage. The extended duration of Ferrar magmatism indicates that it is coeval with the Pliensbachian-Toarcian boundary. Therefore, we speculate that Ferrar (±Karoo) magmatism triggered the Pliensbachian-Toarcian extinction event and contributed to the Toarcian oceanic anoxic event, from which the environment did not begin to recover until only after the waning and cessation of Ferrar magmatic activity at ∼182 Ma, with zircon crystals recording the final flux of magma.

Structural controls on magma pathways in Bora-Baricha-Tullu Moye (BBTM) volcanic system, Main Ethiopian Rift

The Bora-Baricha-Tullu Moye (BBTM) volcanic complex is located at a transitional zone in the Main Ethiopian Rift where tectonic and volcanic features show complex interplays. We mapped and characterised volcanic and tectonic features using high-resolution digital elevation models and performed morphometric and vent spatial distribution analyses. Structural analysis reveals NNE–SSW, NE–SW, and NW–SE trending faults in the region. The dominant post-caldera volcanic landforms are lava domes, pumice cones, scoria cones, maars, obsidian coulees and lava flows, which have distinct morphological characteristics. Vent elongation and alignment highlight close association between these landforms and the caldera(s) as well as with tectonic structures, suggesting these structures acted as the main magma pathways during the BBTM recent eruptions. We estimate that during the entire BBTM post-caldera phase a total bulk volume of 10.9 km3 of material was erupted. This would represent a time-averaged magma flux of 0.05 km3 ky-1 in the BBTM.

Alkaline magmas in shallow arc plutonic roots: a field and experimental investigation of hydrous cumulate melting in the southern Adamello batholith

Despite the first-order importance of crystallisation–differentiation for arc magma evolution, several other processes contribute to their compositional diversity. Among them is the remelting of partly crystallised magmas, also known as cumulate melting or ‘petrological cannibalism’. The impact of this process on the plutonic record is poorly constrained. We investigate a nepheline-normative dyke suite close to the Blumone gabbros, a large amphibole-gabbro unit of the Tertiary Southern Alpine Adamello igneous complex. The compositions of the studied dykes are characterised by low SiO2 (43–46 wt. %), MgO (5.0–7.2 wt. %), Ni (18–40 μg/g), and high Al2O3 (20.2–22.0 wt. %) contents. Phenocrystic plagioclase in these dykes exhibits major, trace, and Sr isotope compositions similar to Blumone cumulate plagioclase, suggesting a genetic link between the nepheline-normative dykes and the amphibole-gabbro cumulates. We tested this hypothesis by performing saturation experiments on a nepheline-normative dyke composition in an externally heated pressure vessel at 200 MPa between 975 and 1100 °C at fO2 conditions close to the Ni–NiO buffer. Plagioclase and spinel are near-liquidus phases at and above 1050 °C, contrasting with the typical near-liquidus olivine ± spinel assemblage in hydrous calc-alkaline basalts. The alkaline nature of the dykes results from the abundance of amphibole in the protolith, consistent with melting of amphibole-gabbro cumulates. We modelled the heat budget from the repeated injection of basaltic andesite into a partly crystallised amphibole-gabbro cumulate. The results of this model show that no more than 7% of the cumulate pile reaches temperatures high enough to produce nepheline-normative melts. We propose that such nepheline-normative dykes are a hallmark of hydrous cumulate melting in subvolcanic plumbing systems. Therefore, ne-normative dykes in arc batholiths may indicate periods with high magma fluxes.

Tracking changes in magma transport from very-long-period seismic signals at Piton de la Fournaise volcano

Changes in magma properties and transport geometry can have a direct impact on volcanic activity. However, such variations can be difficult to track during eruptions. We report previously undetected very-long-period (VLP) signals at Piton de la Fournaise that can be used to probe changes in magma transport. Source analysis of VLP events during the August-October 2015 eruption indicates a source depth of about 0.9-1.2 km and points to the resonance of the magma dike feeding the eruption. The evolution of the resonance period reveals a shortening of the dike when the magma flux decreases at the end of the eruption. VLP events are actually quite frequent at Piton de la Fournaise: all eruptions analyzed in this study exhibit VLP signals that are indicative of rapid drops in lava discharge. This work encourages the detection of VLP signals to monitor changes in magma flow during volcanic eruptions, and anticipate the corresponding evolution in effusive activity.

Two Distinct Magma Storage Regions at Ambrym Volcano Detected by Satellite Geodesy

AbstractThe flux of eruptible magma into a magmatic plumbing system influences eruption size and timing. If magma transfer is possible between two hydraulically‐connected magma lenses, system destabilization can tap a larger magma volume than stored in any one melt lens. This study identifies two distinct magma reservoirs beneath Ambrym, a basaltic island volcano in Vanuatu, during the time period February 2019 to January 2022. Using InSAR time series and a data assimilation approach, we estimate pressure changes within two reservoirs (located 5–7 and 4–6 km b.s.l.). Furthermore, a theoretical model demonstrates that the reservoirs may not currently be hydraulically connected, despite evidence of physical mixing of magma derived from each reservoir during the December 2018 eruption. These findings further our understanding of how magmatic plumbing systems at basaltic calderas may change after rift‐zone eruptions.

Analysis of magma flux and eruption intensity during the 2021 explosive activity at La Soufrière, St Vincent, West Indies

Abstract Real-time Seismic Amplitude Measurement signals and eruption cloud height measurements were used to estimate peak intensities of 40 explosive events during the 8–22 April 2021 activity of La Soufrière volcano. We estimated magma supply rates and erupted volumes in each explosion, characterized uncertainty by stochastic modelling and identified four eruptive stages. Stage 1 included an intense period of 9.5 hours with 11 explosive events with peak eruption intensity between 2000 and 4000 m 3 s −1 and magma supply rate reaching 828 m 3 s −1 . Twelve high-intensity explosions ( c. 4000 m 3 s −1 ) occurred in Stage 2 with average magma supply rate of 251 m 3 s −1 . Stage 3 involved both declining intensity and magma supply rate and lengthening repose periods between explosions. Stage 4 involved three much weaker explosions. The total erupted volume of magma is estimated at 38.5 × 10 6 m 3 (90% credible interval: [22.0 .. 61.9] × 10 6 m 3 ) consistent with independent estimates from analysis of tephra deposits and volcano subsidence sourced at c. 6 km depth. The 150-fold increase in magma supply rate, from the preceding effusive phase to Stage 1 of the explosive phase, is attributed to replacement of very high-viscosity degassed magma occupying the shallow conduit system with new, lower-viscosity, volatile-rich magma from the magma chamber.

Depth of Magma Storage Under Iceland Controlled by Magma Fluxes

AbstractThe compositions of volcanic materials are sensitive to physical conditions in the underlying magmatic system. When basaltic melts are saturated in olivine‐plagioclase‐augite prior to eruption, their compositions can be used to estimate the pressure at which they last equilibrated. We developed PyOPAM, an open‐source tool that runs in Python, and use this refreshed liquid‐barometer to investigate the relationship between final depths of magma storage and magma flux. We first tested PyOPAM using 312 experimental glasses compiled from literature and found that the 1σ uncertainty is 1.13 kbar (±3 km). PyOPAM was then applied to a data set of 13,400 analyses from Iceland, where suspected controls on magma flux are well constrained. Of these, 3807 analyses return robust pressure estimates, constraining final pre‐eruptive magma storage depths for 23 of the 30 Icelandic volcanic systems. Our results indicate that magma storage pressures on Iceland are linked to melt‐flux from the mantle. This finding is consistent with previous models linking storage depths and spreading rates on the global mid‐ocean ridge system. In addition, we provide clear evidence that the magma flux, rather than spreading rate alone, is the key control on the distribution of melt at spreading centers. Increased melt flux is associated with shoaling of pre‐eruptive storage depths, indicating that mantle melt fluxes dictate the long‐term stabilization of extensive magmatic storage regions at depths shallower than 10 km. Quantitative relationships between mantle melt flux and storage depths can be used to test computational models of transcrustal magmatic systems.

Surges in volcanic activity on the Moon about two billion years ago

The history of mare volcanism critically informs the thermal evolution of the Moon. However, young volcanic eruptions are poorly constrained by remote observations and limited samples, hindering an understanding of mare eruption flux over time. The Chang’e-5 mission returned the youngest lunar basalts thus far, offering a window into the Moon’s late-stage evolution. Here, we investigate the mineralogy and geochemistry of 42 olivine and pyroxene crystals from the Chang’e-5 basalts. We find that almost all of them are normally zoned, suggesting limited magma recharge or shallow-level assimilation. Most olivine grains record a short timescale of cooling. Thermal modeling used to estimate the thickness and volume of the volcanism sampled by Chang’e-5 reveals enhanced magmatic flux ~2 billion years ago, suggesting that while overall lunar volcanic activity may decrease over time, episodic eruptions at the final stage could exhibit above average eruptive fluxes, thus revising models of lunar thermal evolution.