Magma Ascent Research Articles

Magmatic to solid-state evolution of a shallow emplaced agpaitic tinguaite (the Suc de Sara dyke, Velay volcanic province, France): implications for peralkaline melt segregation and extraction in ascending magmas

Abstract. In the last decades the mush model has been generalized to the complete trans-crustal magmatic system in which differentiation would be driven by segregation and extraction of trapped melts from crystal-rich mushes. Melt extraction processes involved are porous flow and strain localization, the latter being regarded as the main process acting during transfer through dykes and necks along which high differential stresses are acting on. We combine structural measurements together with petrological analyses and textural observations to constrain the model of emplacement and finally emphasize how shear deformation and strain localization structures promoted the residual melt segregation that occurred in a shallow silica-undersaturated peralkaline intrusion (Suc de Sara, Velay volcanic province, French Massif Central). In this study, we demonstrate that segregation and subsequent extraction of the CO2-rich residual melt occurred during magma ascent and final emplacement of the Suc de Sara tinguaite. Contrasting features of shear deformation between the margins that exhibited different permeabilities highlight that melt segregation started by compaction as a loose packing of emerging microlites and continued with melt filling of an anastomosed C/C′ band network developing in the crystal-rich mush subjected to high shear strain. Subsequent melt extraction throughout the country rock was controlled by the permeability of the hanging wall. Along the western hanging wall of the intrusion, extraction of the residual melt was prevented by the 15 cm thick chilled margin. In contrast, segregated melt circulated through the highly porous and permeable eastern margin, causing the fenitization of the country rock.

On the occurrence of kalsilite in melilite-bearing ultrapotassic lavas from the Roman Province (Vulsini Mts., central Italy)

The Montefiascone Volcanic Complex belongs to the Vulsini Mts. volcanic district, one of the several leucite-bearing ultrapotassic Pleistocene volcanoes of the Roman Province (Central Italy). It is characterized by abundant pyroclastic successions and minor lava flows emplaced during a caldera formation phase occurred ∼0.3–0.2 Ma. The volcanic products span from leucite-basanites and leucite-tephrites to rarer leucitites and melilitites, occasionally associated with very rare kalsilite melilitolite ejecta. Here we report for the first time the presence of primary kalsilite in two lava flows in Forcinella and Feltricci localities, identifying a kamafugitic component among the Montefiascone volcanics.Detailed petrographic, SEM and EPMA studies on ten samples from the two localities revealed, in addition to kalsilite, variable amounts of clinopyroxene, melilite, leucite and nepheline plus rarer olivine coupled with carbonate-bearing segregation pockets in the groundmass. Geochemical and isotopic data confirm the exotic nature of the analysed rocks, showing ultrabasic/basic (SiO2 = 42.15–45.94 wt%) and ultrapotassic (K2O = 6.17–8.48 wt%; K2O/Na2O = 5.3–8.7) compositions, as well as strongly radiogenic 87Sr/86Sr (0.71038–0.711068), poorly radiogenic 143Nd/144Nd (0.512057.512091) and uniform middle radiogenic 206Pb/204Pb (18.76–18.78). Boron isotopes (δ11B ranges from 7.36 to 8.82‰) almost overlaps the typical MORB-OIB range values.These geochemical features recall those of the nearby kamafugite rocks from the Intra-Apennine Province. Major oxide and trace element modelling support fractionation links among the Italian kamafugite types. The overall geochemical and isotopic signatures speak for a common strongly subduction-metasomatized and orthopyroxene-free lithospheric mantle source. As already proposed for some nearby kamafugite centres in central Italy, limestone assimilation during magma ascent is a viable option to explain the mineralogical and whole-rock composition of Montefiascone lavas (i.e., the presence of carbonate plagues in the lava groundmass). This is particularly relevant considering the thick sedimentary sequence pierced by the limited volume of the upwelling magma. On the other hand, many other key features (such as the general lack of crustal xenoliths, the overall incompatible element fractionation and the radiogenic Sr isotopic ratios) clearly require also derivation from subduction-modified mantle sources.

Transition of dome formation to sudden explosive eruptions at Popocatépetl, Mexico: magnetic indicators

Transitions from effusive to explosive activity can increase hazards making it crucial to define early indicators such as changes in the magnetic signals. After more than 80 cycles of crater-dome extrusion and destruction from 1996 on, Popocatépetl volcano (Mexico) experienced changes in its behavior from March 15 to 18 July 2019, when no lava domes were observed. Some of the domes behaved as contained lava flows within the crater floor (pancakes) while others were more irregular-shaped. Activity decreased considerably over this 2019 interval except for the unexpected explosions in March and June, that produced ash plumes reaching up to 14,000 m a.s.l. In order to investigate the causes of the transition from effusive to explosive behavior in March and June, we analyzed the time series from the magnetic monitoring network at Popocatépetl volcano between October 2018 and December 2019. The raw signals were analyzed by weighted differences (WD) based on the elimination of non-local changes from the total intensity values of the geomagnetic field and the discrete-time continuous wavelet transform was used to evaluate the local variations of energy within the time series. The high energy periods (linked to negative magnetic anomalies) are induced by magma ascent associated with movement within the conduit. They indicate that the sudden explosions were due to the ascent of several magma batches that were slowed during ascent and were not able to reach the surface. Changes in the rheology of the lava are linked to the influx of several batches of magma with different compositions as well as to compaction by gas loss when ascending andesitic magma pushed out overlying more viscous degassed magma clearing the conduit, which can explain why these sudden explosions were more energetic. Several geophysical data sets as well as tephra compositions were integrated to support this conclusion. The correlated multiparameters also confirm that geomagnetic volcano monitoring has been essential in understanding the processes that drive the observed changes in eruptive behavior. We present new evidence for the detection of transient events produced by magma ascent and changes in the feeding system of Popocatépetl volcano with wavelet analysis. Detailed vulcanomagnetic processing, especially when it is correlated with other monitoring parameters, provides information on ascending magma and several conduit processes that would otherwise be camouflaged. Ascending batches may precede an eruption but they can also ascend in several pulses indicating how dome growth occurs.

Petrographic and Geochemical Evidence for a Complex Magmatic Plumbing System beneath Bagana Volcano, Papua New Guinea

Abstract Bagana is a persistently active stratovolcano located on Bougainville Island, Papua New Guinea. Characteristic activity consists of prolonged lava effusion over months to years, with occasional shifts to explosive vulcanian or sub Plinian eruptions that threaten surrounding communities. Satellite observations have shown that Bagana is a major SO2 emitter, particularly during eruptive intervals. Despite persistent and potentially hazardous activity, no previous geophysical, petrological, or geochemical studies have constrained the magma storage conditions and reservoir processes at Bagana. To address this knowledge gap, we present new bulk rock major, trace element, and radiogenic isotope data, plus mineral phase major element compositions, for Bagana lavas erupted in 2005 and 2012 and ash erupted in 2016. We use our new data to understand the magmatic processes controlling the typical effusive activity and provide the first estimates of magma storage conditions beneath Bagana. The basaltic andesite bulk rock compositions (56–58 wt% SiO2) of our Bagana lavas reflect accumulation of a plagioclase + clinopyroxene + amphibole + magnetite + orthopyroxene crystal cargo by andesitic-dacitic (57–66 wt% SiO2) carrier melts. Constraints from clinopyroxene and amphibole thermobarometry, amphibole hygrometry, and experimental petrology suggest that the high-An plagioclase + clinopyroxene + amphibole + magnetite assemblage crystallizes from basaltic-basaltic andesite parental magmas with >4 wt% H2O, over a temperature interval of ~1100–900°C, at pressures of ~130–570 MPa, corresponding to ~5–21 km depth. Continued crystallization in the magma storage region at ~5–21 km depth produces andesitic to dacitic residual melts, which segregate and ascend towards the surface. These ascending melts entrain a diverse crystal cargo through interaction with melt-rich and mushy magma bodies. Degassing of carrier melts during ascent results in crystallization of low-An plagioclase and the formation of amphibole breakdown rims. The radiogenic isotope and trace element compositions of Bagana lavas suggest that parental magmas feeding the system derive from an enriched mantle source modified by both slab fluids and subducted sediments. Our findings suggest that the prolonged lava effusion and persistently high gas emissions that characterise Bagana’s activity in recent decades are sustained by a steady state regime of near-continuous ascent and degassing of magmas from the crustal plumbing system. Our characterisation of the Bagana magmatic plumbing system during effusive activity provides a valuable framework for interpreting ongoing monitoring data, and for identifying any differences in magmatic processes during any future shift to explosive activity.

Pleistocene to recent evolution of Mocho-Choshuenco volcano during growth and retreat of the Patagonian Ice Sheet

Mocho-Choshuenco volcano (39.9°S, 72.0°W) produced ∼75 explosive eruptions following retreat of the >1.5-km-thick Patagonian Ice Sheet associated with the local Last Glacial Maximum (LGM, from 35 to 18 ka). Here, we extend this record of volcanic evolution to include pre- and syn-LGM lavas that erupted during the Pleistocene. We establish a long-term chronology of magmatic and volcanic evolution and evaluate the relationship between volcanism and loading/unloading of the Patagonian Ice Sheet via twenty-four 40Ar/39Ar and two 3He age determinations integrated with stratigraphy and whole-rock compositions of lava flows and glass compositions of tephra. Our findings reveal that the edifice is much younger than previously thought and preserves 106 km3 of eruptive products, of which 50% were emplaced immediately following the end of the penultimate glaciation and 20% after the end of the LGM. A period of volcanic inactivity between 37 and 26 ka, when glaciers expanded, was followed by the eruption of incompatible element-rich basaltic andesites. Several of these syn-LGM lavas dated between 26 and 16 ka, which crop out at 1500−1700 m above sea level, show ice contact features that are consistent with emplacement against a 1400- to 1600-m-thick Patagonian Ice Sheet. Small volume dacitic eruptions and two explosive rhyolitic eruptions dominate the volcanic output from 18 to 8 ka, when the Patagonian Ice Sheet began to retreat rapidly. We hypothesize that increased lithostatic loading as the Patagonian Ice Sheet grew prohibited dike propagation, thus stalling the ascent of magma, promoting growth of at least three discrete magma reservoirs, and enhancing minor crustal assimilation to generate incompatible element-rich basaltic andesitic to dacitic magmas that erupted between 26 and 17 ka. From an adjacent reservoir, incompatible element-poor dacites erupted from 17 to 12 ka. These lava flows were followed by the caldera-forming eruption at 11.5 ka of 5.3 km3 of rhyolite from a deeper reservoir atop which a silicic melt lens had formed and expanded. Subsequent eruptions of oxidized dacitic magmas from the Choshuenco cone from 11.5 to 8 ka were followed by andesitic to dacitic eruptions at the more southerly Mocho cone, as well as small flank vent eruptions of basaltic andesite at 2.5 and 0.5 ka. This complex history reflects a multi-reservoir plumbing system beneath Mocho-Choshuenco, which is characterized by depths of magma storage, oxidation states, and trace element compositions that vary over short periods of time (<2 k.y.).

Volcanic Unrest After the 2021 Eruption of La Palma

AbstractLa Palma, Canary Islands, had its largest historical eruption in 2021. From January 2022 to May 2023 there were >2,100 seismic events, primarily at depths ≤20 km, prompting us to update the deformation and modeling study, using interferometric synthetic aperture radar observations and a last generation interpretation tool. We detect the evolution of the remaining magmatic body in the SW portion of the island, with arrival of new magma moving into the oceanic crust out to sea, and a pressurized zone in the central‐eastern area, at regions of structural weakness. The current source characteristics have some similarities to the early stage dynamics prior to the 2021 eruption. Operational and multidisciplinary studies must continue to monitor either their stabilization or growth and destabilization. The ability to identify magma ascent using only deformation data over short time periods allows us to characterize unrest patterns and provide new insights into volcanic processes.

The Chain of Processes Forming Porphyry Copper Deposits—An Invited Paper

Abstract Porphyry-related mineral deposits are giant geochemical anomalies in the Earth’s crust with orders-of-magnitude differences in the content and proportion of the three main ore metals Cu, Au, and Mo. Deposit formation a few kilometers below surface is the product of a chain of geologic processes operating at different scales in space and time. This paper explores each process in this chain with regard to optimizing the chances of forming these rare anomalies. On the lithosphere scale, deposits with distinct metal ratios occur in provinces that formed during brief times of change in plate motions. Similar metal ratios of several deposits in such provinces compared with global rock reservoirs suggest preceding enrichment of Au or Mo in lithospheric regions giving rise to distinct ore provinces. The largest Cu-dominated deposits and provinces are traditionally explained by selective removal of Au during generation or subsequent evolution of mantle magmas, but the possibility of selective Cu pre-enrichment of lithosphere regions by long-term subduction cannot be dismissed, even though its mechanism remains speculative. Evolution of hydrous basaltic melts to fertile magmas forming porphyry Cu deposits requires fractionation toward more H2O-rich magmas in the lower crust, as shown by their adakite-like trace element composition. The prevailing interpretation that this fractionation leads to significant loss of chalcophile ore metals by saturation and removal of magmatic sulfide might be inverted to a metal enrichment step, if the saturating sulfides are physically entrained with the melt fraction of rapidly ascending magmas. Ascent of fertile magma delivers a large mass of H2O-rich ore fluid to the upper crust, along points of weakness in an overall compressive stress regime, within a limited duration as required by mass and heat balance constraints. Two mechanisms of rapid magma ascent are in debate: (1) wholesale emplacement of highly fractionated and volatile-rich granitic melt into a massive transcrustal channelway, from which fluids are exsolved by decompression starting in the lower crust, or (2) partly fractionated magmas filling a large upper crustal magma chamber, from which fluids are expelled by cooling and crystallization. Transfer of ore-forming components to a hydrothermal ore fluid is optimized if the first saturating fluid is dense and Cl rich. This can be achieved by fluid saturation at high pressure, or after a moderately H2O rich intermediate-composition melt further crystallizes in an upper crustal reservoir before reaching fluid saturation. In either case, metals and S (needed for later hydrothermal sulfide precipitation) are transferred to the fluid together, no matter whether ore components are extracted from the silicate melt or liberated to the ore fluid by decomposition of magmatic sulfides. Production and physical focusing of fluids in a crystallizing upper crustal magma chamber are controlled by the rate of heat loss to surrounding rocks. Fluid focusing, requiring large-scale lateral flow, spontaneously occurs in mushy magma because high water content and intermediate melt/crystal ratio support a network of interconnected tubes at the scale of mineral grains. Calculated cooling times of such fluid-producing magma reservoirs agree with the duration of hydrothermal ore formation measured by high-precision zircon geochronology, and both relate to the size of ore deposits. Ore mineral precipitation requires controlled flow of S- and metal-rich fluids through a vein network, as shown by fluid inclusion studies. The degree of hydrothermal metal enrichment is optimized by the balance between fluid advection and the efficiency of cooling of the magmatic fluid plume by heat loss to convecting meteoric water. The depth of fluid production below surface controls the pressure-temperature (P-T) evolution along the upflow path of magmatic fluids. Different evolution paths controlling density, salinity, and phase state of fluids contribute to selective metal precipitation: porphyry Au deposits can form at shallow subvolcanic levels from extremely saline brine or salt melt; high-grade Au-Cu coprecipitation from coexisting and possibly rehomogenizing brine and vapor is most efficient at a depth of a few kilometers; whereas fluids cooling at greater depth tend to precipitate Cu ± Mo but transport Au selectively to shallower epithermal levels. Exhumation and secondary oxidation and enrichment by groundwater finally determine the economics of a deposit, as well as the global potential of undiscovered metal resources available for future mining.

Oligocene magmatism in the eastern Sakarya Zone, NE Turkey: Geodynamic implications

During the Cenozoic, the eastern Sakarya Zone, northeastern Turkey, has witnessed voluminous Eocene and subsidiary Miocene magmatism. By contrast, except in one previous study, Oligocene magmatism has not been reported. Here, we present new data about an Oligocene dyke swarm from the Artvin district in the easternmost part of the Sakarya Zone. Our study improves our understanding of the Oligocene magmatic and geodynamic evolution of the central Tethyan orogenic belt.The Oligocene dykes consist of two distinct rock types: quartz-bearing trachytic and adakitic trachyandesite/andesite, dated at 31.0 ± 0.2 Ma and 31.1 ± 0.2 Ma to 29.2 ± 0.3 Ma by U-Pb zircon geochronology, respectively. Geochemical analysis reveals enrichment in LREE and MREE for quartz-bearing trachytic dykes, whereas adakitic trachyandesite/andesite dykes show depletion in MREE and HREE. Radiogenic isotopic ratios (87Sr/86Sr(i) and 143Nd/144Nd(i)) and trace element data suggest that the source of Oligocene magmatism was linked to low degree (≈10–15 %) partial melting of a metasomatised lithospheric mantle previously enriched by slab-derived melts and fluids. During the magma ascent in the crust, parental magmas of the trachyandesite/andesite dykes have assimilated crustal rocks and undergone fractional crystallization dominated by amphibole. By contrast, quartz-bearing trachytic dykes have recorded plagioclase fractionation.Local extension, in conjunction with continuous foundering of delaminated thickened lower crust beneath the eastern Sakarya Zone, which likely triggered asthenospheric mantle upwelling and resulted in sporadic Oligocene magmatism.Our study highlights the distinct tectonic and magmatic setting of the rare Oligocene magmatism in the eastern Sakarya Zone compared to Oligocene magmatism in the Tauride-Anatolide and South Armenian blocks.

The Andean Southern Volcanic Zone: a review on the legacy of the latest volcanic eruptions

The Andean Southern Volcanic Zone (SVZ) concentrates many of the most active volcanoes of the Andean continental arc, as well as the region’s most recent and impactful volcanic eruptions. In this contribution, we briefly revise the general characteristics of the SVZ volcanism and provide a synthesis of the scientific findings related to the latest volcanic eruptions (430 peer-reviewed publications with over 9,000 citations, with large-magnitude (VEI 4-5) eruptions being the most studied. Our study shows that SVZ research has been primarily focused on environmental and atmospheric impacts (29%), eruption descriptions and physical volcanology (20%), volcanic hazard and risk assessments (15%), and other investigations complementary to volcanology. Whereas the least silicic eruptions (e.g., Llaima 2008-2009 and Villarrica 2015) shed light on magma replenishment and degassing dynamics controlling eruption styles, intermediate eruptions (andesitic-dacitic) offered clues on either rapid or slow eruption initiation, with relevant findings on phreatic-to-magmatic style transitions and eruption triggering mechanisms. On the other hand, silicic (i.e., rhyolite-rhyodacite) eruptions provided unique observations on rapid magma ascent, high-rate magma extrusion, rheology, fragmentation processes, and style transitions. These recent eruptions have also inspired a new generation of tephrochronological, tephrostratigraphical, and physical volcanology studies, aimed at assessing the long-term (kyr-scale) evolution of the volcanic systems and their associated hazards. We debate how the knowledge gained from research and the long-term human coexistence with volcanoes are relevant to reducing volcanic risk in the SVZ. Finally, we discuss how challenges and opportunities emerging from other disciplines can complement our understanding of volcanism in this active region.

Stressful crystal histories recorded around melt inclusions in volcanic quartz

Magma ascent and eruption are driven by a set of internally and externally generated stresses that act upon the magma. We present microstructural maps around melt inclusions in quartz crystals from six large rhyolitic eruptions using synchrotron Laue X-ray microdiffraction to quantify elastic residual strain and stress. We measure plastic strain using average diffraction peak width and lattice misorientation, highlighting dislocations and subgrain boundaries. Quartz crystals across studied magma systems preserve similar and relatively small magnitudes of elastic residual stress (mean 53–135 MPa, median 46–116 MPa) in comparison to the strength of quartz (~ 10 GPa). However, the distribution of strain in the lattice around inclusions varies between samples. We hypothesize that dislocation and twin systems may be established during compaction of crystal-rich magma, which affects the magnitude and distribution of preserved elastic strains. Given the lack of stress-free haloes around faceted inclusions, we conclude that most residual strain and stress was imparted after inclusion faceting. Fragmentation may be one of the final strain events that superimposes stresses of ~ 100 MPa across all studied crystals. Overall, volcanic quartz crystals preserve complex, overprinted deformation textures indicating that quartz crystals have prolonged deformation histories throughout storage, fragmentation, and eruption.

Deciphering magmatic processes from plagioclase, clinopyroxene and amphibole chemistry and textures: A case study of a basaltic lava flow in the Basque-Cantabrian Basin (Northern Spain)

This paper addresses the study of a pillow lava interbedded with Late Albian-Early Cenomanian sediments that crops out in Armintza (Bizkaia, Northern Spain). The lava flow is an alkaline basalt with abundant macrocrysts of clinopyroxene, kaersutite, Ca-rich plagioclase (50-86% An) and ilmenite, which display a variety of textures and complex zoning patterns indicative of open-system magmatic behaviour. Macrocryst cores are likely to be inherited antecrysts that underwent complex processes under deep pre-eruptive conditions (≈ 700-800 MPa). Microcrysts and macrocryst rims formed during magmatic ascent and emplacement at shallower levels (≈ 35 MPa). Hypothetical melts in equilibrium with clinopyroxenes and amphiboles have trace element compositions like metasomatic vein melts containing amphibole, and their patterns overlap with those of the Armintza pillow lava. This suggests a metasomatised lithospheric mantle with amphibole-rich veins as a potential source for the alkaline basaltic melt. It is even conceivable that the Armintza pillow lava and other alkaline volcanic manifestations of the Basque-Cantabrian Basin were part of the same magma plumbing system through which a series of time-limited eruptions of different batches of magma ascended from the lithospheric mantle to the upper crust during the Albian to the Santonian.

Time scales of olivine storage and transport as revealed by diffusion chronometry at Waitomokia Volcanic Complex, Auckland Volcanic Field, New Zealand

Detailed knowledge of the pre-eruptive time scales associated with magma storage and transport is vital to improve volcanic hazard forecasting in active volcanic regions. However, quantification of the timescales of volcanic processes at mafic volcanic centres in continental intraplate settings is challenging, despite them being a source of significant hazards for human populations and infrastructure due to their limited predictability in space and time. We conducted a detailed petrological study to investigate the time scales of olivine storage and transfer throughout the eruption sequence of Waitomokia Volcanic Complex, a tuff ring and scoria cone complex in the Auckland Volcanic Field. Olivine crystal textures and compositions were determined from stratigraphically-constrained samples of the volcanic complex, from the initial phreatomagmatic phase to the final magmatic phase. Olivine crystals are typically <300 μm in length and characterised by skeletal morphologies, displaying chemical zoning in forsterite (Fo = 100*Mg/[Mg + Fe]; mol%), CaO, MnO and NiO wt% contents. We classified olivine into three major groups based on their Fo core compositions: (1) normally zoned crystals with high Fo content (Fo > 85), (2) crystals with intermediate Fo contents (84–81), and (3) reversely zoned crystals with lower Fo core content (<80). Olivine chemical zoning (diffusion) profiles were modelled in the context of a specific magmatic environment linked with changes in thermodynamic variables during storage (temperature, pressure, and oxygen fugacity). We propose that the normally zoned olivine crystals grew in one magmatic environment (ME1), which subsequently intruded into a more evolved (lower MgO) environment (ME2), where they interacted and were stored for up to 135 days before their eruption. During magma ascent to the surface, a second magma mixing event occurred between ME2 and magma within a third magmatic environment (ME3), forming reversely-zoned olivine crystals yielding notably shorter ascent times of approximately a few days. The rocks from the opening phreatomagmatic phase of the eruption show a larger range in olivine group types compared to the final magmatic phase, where those from the deeper ME1 are more abundant. The short time scales of magma transport obtained in our study, on the order of days to months, should be informative of the warning times that may be encountered between the onset of volcanic unrest and an eruption in the Auckland Volcanic Field.

Conjugate shear model for emplacement of mid-crustal plutons during dextral transpression, southern Nova Scotia, Canada: tectonic and petrological consequences

The South Mountain Batholith (SMB) is a syn-tectonic composite batholith emplaced in the upper crust within the Meguma terrane between 380 and 370 Ma during the later stages of the Neoacadian orogeny. Coeval plutons in southern Nova Scotia are surrounded by mid-crustal (∼4 kbar) andalusite–staurolite aureoles and are discordant to northeast-trending, regional Neoacadian folds. Detailed field studies, combined with published results, indicate emplacement within a dextral transpressional regime during the transition from distributed (D1) to focused heterogeneous strain (D2), which provided vertical conduits that facilitated magma ascent. The Port Mouton Pluton intruded along P-orientated crustal-scale fractures as a series of subvertical granitic sheets, which were progressively rotated and folded with ongoing dextral shearing. By contrast, the Barrington Passage Pluton intruded between crustal-scale, antithetic (sinistral) P-shear fractures and spread laterally between them as pulsed increments to form a layered, subhorizontal, and sill-like complex resembling the laccolithic structure of the SMB. The SMB was emplaced below the Meguma Supergroup, with magma derived from underthrust Avalon terrane and Silurian Rockville Notch Group. The lack of mantle components in the SMB suggests transpressional orogenesis facilitated conductive crustal heating without significant mantle addition, consistent with low p-wave velocities of Meguma lower crust. 400–355 Ma zircons, recorded either as inherited grains in granites or in felsic granulite xenoliths, imply the Neoacadian thermal anomaly extended for 45 Ma, but magmatism represented only ∼40% of that perturbation.

Three-dimensional density model of the mantle beneath the Ukrainian shield

Purpose. Mantle density models are key tools for understanding the fundamental geological and physical processes occurring within the Earth and are essential to our scientific and applied understanding of the planet. Methodology. The tasks were solved by a complex research method, including analysis and generalization of literary and patent sources, analytical, experimental studies, using computer and mathematical modelling methods. Findings. One-dimensional models simplify the mantle density distribution by assuming that it is uniform only in the vertical direction. This limitation does not allow for horizontal variations in mantle density, which may be important on a regional scale. 3D models are more complex and require more data and computational resources, so their use may be limited. In this study, we present a quasi-three-dimensional model of mantle density beneath the Ukrainian Shield. This 3D model is obtained using a basic set of one-dimensional seismic tomographic velocity models calculated for 21 mantle domains in the depth range from 50 to 2,600 km. The process of converting the P-wave velocity model into a density model includes the following stages: 1) determining seismic boundaries in the mantle based on P-wave velocity curves for each mantle domain; 2) creating a synthetic mantle model beneath the Ukrainian Shield for the P,S-wave velocity curves; 3) solving the Adams-Williamson equation for each domain, considering polynomial corrections to extract heterogeneities during its solution; 4) analysing existing models by comparing the calculated gravitational potential at the central point of the Ukrainian Shield as the standard reference for selecting one of 5 reference models. Here, we focus on the final stages of constructing the mantle density model by: 1) balancing the mass of the upper and lower mantle for each domain when determining density using the Adams-Williamson equation and introducing polynomial corrections; 2) calculating densities for each of the 21 mantle domains and their 3D integration. Originality. The obtained mantle-density model of the Ukrainian Shield aligns well with the division of the mantle into three main layers: lithosphere, upper mantle, and lower mantle. Each of the mantle’s structural layers has its representation pattern in density heterogeneities. Anomalies of decreased density in the lithosphere of the Ukrainian Shield correlate with thermal anomalies, whereas anomalies of increased density correspond to tectonic zones dividing its megablocks. Practical value. Regions of increased density gradient are associated with mantle thrust faults, which in some cases can be boundaries between different petrological formations and serve as channels for magma ascent into the Earth’s crust at certain stages of geological development of the Ukrainian shield and, in turn, be sources of minerals.

Insights Into Venus' Crustal Plateaus From Dyke Trajectories Below Craters

AbstractOn Venus, radar observations of the surface have highlighted two categories of craters: bright‐floored, interpreted as pristine, and dark‐floored, interpreted as being partially filled by lava. While volcanic resurfacing occurs within and outside craters in the plains, it seems mainly concentrated within the interior of dark‐floored craters in the crustal plateaus, suggesting that the magma is negatively buoyant there. Indeed, crater unloading may facilitate vertical ascent of a negatively buoyant magma by decompressing the underlying crust. However, the crater topography also generates a shear stress which would tend to horizontalize the vertical propagation of a dyke. We use numerical simulations of magma ascent in an axisymmetric crater stress field to demonstrate that, depending on the crust thickness and the magma‐crust density contrast, a negatively buoyant magma can indeed erupt only in the crater interior while remaining stored in the crust elsewhere. In particular, we identify four different behaviors depending on if and where a magma‐filled crack ascending below a crater reaches the surface. We draw a regime diagram as a function of two characteristic dimensionless numbers. For eruption to occur only in the crater interior requires a crust thinner than 45 km and a limited range of magma‐crust density contrasts, between 40 and 280 kg m−3 for crust thicknesses between 20 and 45 km, the permissible range decreasing for increasing crustal thicknesses. These results suggest that the crustal plateaus may not be particularly thick and could be slightly differentiated, but probably not very felsic.

The transition of the magma plumbing system of Tianchi shield-forming basalts, Changbaishan Volcanic Field, NE China: Constraints from dynamic Fe-Mg diffusion modelling in olivine

Abstract The depths of crustal reservoirs within volcanic systems may experience transitions over time. Here, we report the crystal and bulk rock compositions of the shield-forming basaltic lavas of the Tianchi composite volcano in the intraplate Changbaishan Volcanic Field, NE China to constrain the crustal magmatic evolution with time. We investigated samples covering the entire basaltic stratigraphic sequence, consisting of the Toudao (TD), Baishan (BS), and Laofangzixiaoshan (LFZ) units from bottom to top, respectively. The core compositions of olivine macrocrysts vary among the three units, i.e., the TD and BS olivine phenocrysts can both be divided into two populations: a high-Fo population (~Fo76-80) and a low-Fo population (~Fo72-74). The LFZ unit only exhibits a high-Fo population (~Fo77-80). Phase equilibria modelling using rhyolite-MELTS suggests that the high-Fo populations were stored at depths of ~20 km for the TD and BS units and ~15 km for the LFZ unit. The low-Fo populations crystallized at shallow depths, i.e., ≤15 km for the TD unit and ≤13 km for the BS unit. We employ a dynamic Fe-Mg interdiffusion modelling with constantly adapting boundary conditions in zoned olivine macrocrysts to constrain the magmatic environments and timescales during the pre-eruption and post-eruption, enabling clarify the magmatic histories recorded by two olivine populations. The dynamic Fe-Mg interdiffusion modelling considers the variable boundary condition caused by crystal growth and composition variation of melts during magma cooling. Calculated results suggest that the high-Fo populations from the TD and BS units recorded prolonged timescales ranging from six months to more than two years with lower cooling rates and slower crystal growth 37 rates. These characteristics reflect a relatively hot and slow-cooling magmatic environment; and the modelled timescales correspond to the sum time including shallow storage, magma ascent, and further cooling within the lava flows. Conversely, the high-Fo population from the LFZ unit and the low-Fo populations from the TD and BS units record shorter timescales (&amp;lt;140 days) with higher cooling rates and faster crystal growth rates. These results indicate relatively cold and highly undercooling magmatic environments; hence the timescales record magma ascent in the conduits and further cooling during lava emplacements. Our study demonstrates that the Tianchi basaltic plumbing system experienced a structural transition over time. In detail, the TD and BS magmas experienced multi-stage stalling and ascent, first accumulating in deep reservoirs and then transferring to shallow reservoirs for storage before the eruption. The LFZ magmas accumulated in a mid-crustal reservoir, followed by a direct ascent to the surface without additional residence.

Banding in the Margins of Basaltic Dykes Indicates Pulsatory Propagation During Emplacement

AbstractBasaltic fissure eruptions, which are the most common type of eruption on Earth, are fed by dykes which mediate magma transport through the crust. Dyke propagation processes are important because they determine the geometry of the transport pathway and the nature of any geophysical signals associated with magma ascent. Here, we investigate small‐scale (mm–cm wide) banding features at the margins of dykes in the Teno Massif (Tenerife, Spain) and the Columbia River Basalt Province (CRBP) (USA). Similar marginal bands have been reported for dykes in numerous localities around the world. Dyke margins record valuable information about propagation because they are the first material to solidify against the host rock at the propagating dyke tip. We find that the marginal bands are defined by cyclic variations in phenocryst concentration and vesicularity, and we infer that these cyclic variations in texture are a product of cyclic variations in magma flow rates and pressures within the dyke tip. This indicates that dyke emplacement occurs in pulses, with propagation repeatedly hindered by the rapid cooling and solidification of magma in the narrow dyke tip. Using a 1D conduction model, we estimate the time taken for each band to cool and solidify, which provides a timescale of several minutes to tens of minutes for the pulses. The occurrence of similar bands in various volcanic settings suggests that pulsatory propagation is a common, if not ubiquitous, process associated with dyke emplacement.

CO2 Flushing Triggers Paroxysmal Eruptions at Open Conduit Basaltic Volcanoes

AbstractOpen conduit volcanoes erupt with the highest frequency on Earth. Their activity is characterized by an outgassing flux that largely exceeds the gas that could be released by the erupted magma; and by frequent small explosions intercalated by larger events that pose a significant risk to locals, tourists, and scientists. Thus, identifying the signs of an impending larger explosion is of utmost importance for the mitigation of volcanic hazard. Larger explosive events have been associated with the sudden ascent of volatile rich magmas, however, where and why magma accumulates within the plumbing system remains unclear. Here we show that the interaction between CO2‐rich fluids and magma spontaneously leads to the accumulation of volatile‐rich, low density and gravitationally unstable magma at depth, without the requirement of permeability barriers. CO2‐flushing forces the exsolution of water and the increase of magma viscosity, which proceeds from the bottom of the magma column upward. This rheological configuration unavoidably leads to the progressive thickening of a gas‐rich and low density (i.e., gravitationally unstable) layer at the bottom of the feeding system. Our calculations account for observations, gas monitoring and petrological data; moreover, they provide a basis to trace the approach to deeply triggered large or paroxysmal eruptions and estimate their size from monitoring data. Our model is finally applied to Stromboli volcano, an emblematic example of open conduit volcano, but can be applied to any other open conduit volcano globally and offers a framework to anticipate the occurrence of unexpectedly large eruptions.

Evaluating the Role of Titanomagnetite in Bubble Nucleation: Rock Magnetic Detection and Characterization of Nanolites and Ultra‐Nanolites in Rhyolite Pumice and Obsidian From Glass Mountain, California

AbstractWe document the presence, composition, and number density (TND) of titanomagnetite nanolites and ultra‐nanolites in aphyric rhyolitic pumice, obsidian, and vesicular obsidian from the 1060 CE Glass Mountain volcanic eruption of Medicine Lake Volcano, California, using magnetic methods. Curie temperatures indicate compositions of Fe2.40Ti0.60O4 to Fe3O4. Rock‐magnetic parameters sensitive to domain state, which is dependent on grain volume, indicate a range of particle sizes spanning superparamagnetic (&lt;50–80 nm) to multidomain (&gt;10 μm) particles. Cylindrical cores drilled from the centers of individual pumice clasts display anisotropy of magnetic susceptibility with prolate fabrics, with the highest degree of anisotropy coinciding with the highest vesicularity. Fabrics within a pumice clast require particle alignment within a fluid, and are interpreted to result from the upward transport of magma driven by vesiculation, ensuing bubble growth, and shearing in the conduit. Titanomagnetite number density (TND) is calculated from titanomagnetite volume fraction, which is determined from ferromagnetic susceptibility. TND estimates for monospecific assemblages of 1,000 nm–10 nm cubes predict 1012 to 1020 m−3 of solid material, respectively. TND estimates derived using a power law distribution of grain sizes predict 1018 to 1019 m−3. These ranges agree well with TND determinations of 1018 to 1020 m−3 made by McCartney et al. (2024), and are several orders of magnitude larger than the number density of bubbles in these materials. These observations are consistent with the hypothesis that titanomagnetite crystals already existed in extremely high number‐abundance at the time of magma ascent and bubble nucleation.

Forecasting the onset of volcanic eruptions using the increase in seismicity during magma ascent

The Failure Forecast Method (FFM) has been used to forecast the onset of volcanic eruptions with varying degrees of success. The method involves fitting its empirical equation to precursory observables, e.g., seismic data. Current models explaining the empirical equation assume that the seismic observables used (e.g., seismic event-rate, seismic moment-rate or seismic energy-rate) are related to accelerated fracture growth. In these models, such fracture growth results in an apparent increase in seismic activity. In this study, however, we propose an alternative explanation for the increase in seismicity culminating in volcanic eruptions. We suggest that the increase in seismicity, particularly in the seismic amplitudes, can result from ascending seismic sources towards the surface and that the increase can be used for eruption forecasts using the FFM. Our argument is based on the need for magma to migrate towards the surface to produce an eruption.To test our argument, we used synthetic and field/real data and fitted them with the FFM model. The synthetic data is generated using the body wave attenuation equation, and the real data is from the Piton de la Fournaise eruption on January 2nd, 2010. We used seismic amplitudes and seismic amplitude ratios as the precursory observables and used a combination of Simulated Annealing and Bayesian Inversion algorithms to obtain the best-fit parameters of the FFM model and their posterior distributions. For the synthetic cases, we found that the method gave reasonable estimates of failure or eruption time, suggesting that the seismic amplitudes and the seismic amplitude ratios can be used for forecasting using the FFM. For the field/real cases, the results gave information about the timing of magma chamber failure, magma propagation, and the eruption onset. Our results highlight an example of modelling volcano monitoring data with the FFM that gives reasonable estimates of eruption time even if the data may not fully represent failure processes. We suggest that the result may still be reasonable owing to the empirical and asymptotic nature of the FFM combined with the underlying physics of magma and seismic migrations.