Search

In volcanic belts, magma influx into magma chambers generates excess pressure, amplifying the initial stress field to cause crustal deformation with significant ground displacements, which manifests in topographic relief. Quantifying such volcano-driven ground surface displacements is a fundamental requirement to embark on a criticality analysis of volcanotectonic events and associated hazard monitoring strategies. This study theoretically examines the underlying dynamics of surface displacements in a volcanic plumbing system comprising multiple magma chambers. The classical Mogi equation is extended to derive a set of analytical solutions to evaluate surface displacements as a function of separations between two off-axis chambers, measured along both horizontal and vertical directions. The resulting surface displacement plots, from the analytical solutions are compared with those calculated from a set of finite element (FE) model simulations run with the same parameters considered for the analytical formulations. Both the analytical and FE results suggest that horizontal (Sh) and vertical (Sv) separations of magma chambers largely control the vertical (Uz) and lateral (Ur) ground-displacement components. Spatially varying Uz attains its peak value at a specific location above the chambers, but increasing horizontal separation (Sh ~ 10 km) transforms the single-peak Uz pattern to a weakly developed double-peak Uz pattern, which eventually give way to two prominent high-amplitude peaks above the chambers when Sh ~ 25 km. Similarly, a large vertical separation (Sv ~ 6 km) yields double peaks in the Uz profile, which merge to form a single peak for small Sv (~ 1.5 km). The FE model results are used to map the stress fields around the two magma chambers to show that inter-chamber mechanical interaction can influence the deformation behaviour around the chambers, depending on Sh and Sv magnitudes. Finally, the model estimates are evaluated using available reports on the naturally occurring volcanoes: Teide volcano (Tenerife, Spain) and Long Valley Caldera (USA). [ABSTRACT FROM AUTHOR]

The studies are focused in the Mara–Kamenka–Uvat interfluve of the Biryusa Uplift of the basement of the Siberian Platform, where the Mara paleovolcano was recognized during the prospecting works for manganese more than 50 years ago. The specific volcano-sedimentary rocks are considered in the structure of the Late Riphean Karagas Group. Our mineralogical and petrographic studies allowed us to establish abundant high-K pyroclastic rocks, ignimbrites, and trachybasalts, which indicate a subaerial explosive character of the Mara volcano. The age of high-K volcanism of 640 Ma is based on U–Th–Pb dating of zircon. The Lu–Hf isotope systematics of zircon indicates a link of volcanism with mantle magmas. The composition and the period of the formation of rocks prevent their correlation (as was considered before) with sedimentary Late Riphean quartz and quartz–feldspar sandstones of the Karagas Group and dolerites of the Nersa intrusive complex. Specific mineralogical–petrographic features of rocks allow their use as a regional stratigraphic reference. [ABSTRACT FROM AUTHOR]

We report whole-rock chemistry, mineral chemistry, and volatile content from Villarrica volcano's major recent paroxysms and background activity. Composition of the volcanic products are basalt to basaltic andesite with whole-rock SiO2 content between 50 and 56 wt%, and a mineralogy dominated by olivine (Fo71-80), clinopyroxene (Mg# ~ 50) and plagioclase (An60–80). Volatile contents in melt inclusions are up to 1.5 wt% H2O, 500 ppm CO2, 1230 ppm sulfur and 580 ppm chlorine. Regardless of the type of activity, there are no substantial variations in whole-rock composition or the volatile content when the activity switches from background activity to a major paroxysm, strongly suggesting that this shift does not just depend on the arrival of new magma in the shallow magmatic system. Geothermobarometry constrains crystallization of the major mineral phases at various depths between 3 and 12.7 km, suggesting that degassing of a volatile-rich recharge magma occurs deeper than 12 km, producing efficient mixing throughout the whole system, and sustaining the lava lake activity in Villarrica's summit crater. The occurrence of a permanent lava lake also suggests that the magma recharge must be close to continuous and therefore sudden changes between background and paroxysmal volcanic activity are likely controlled by relatively small changes in the rate of recharge and/or the volatile release rate in the magmatic system. This has important implications for the understanding of eruption triggers and the forecasting of volcanic eruptions. [ABSTRACT FROM AUTHOR]

The Southern Volcanic Zone (SVZ) in Chile is an active continental arc with a complex history of volcanism, where a range of magmatic compositions have been erupted in a variety of styles. In the Central SVZ, both monogenetic and polygenetic volcanoes exist, in close proximity to the Liquiñe-Ofqui Fault System (LOFS), but with variable local stress states. Previous studies have inferred varying crustal storage timescales, controlled by the orientation of volcanic centres relative to the N-S striking LOFS and σHMax in this region. To assess the relationship between volcanism and crustal stress states affected by large-scale tectonic structures and edifice controls, we present whole rock geochemical data, to ensure consistency in source dynamics and crustal processing, mineral-specific compositional data, thermobarometry, and Fe–Mg diffusion modelling in olivine crystals from mafic lavas, to assess ascent timescales, from the stratovolcanic edifice of Puyehue-Cordón Caulle and proximal small eruptive centres. Textural observations highlight differences in crystal maturation timescales between centres in inferred compression, transpression, and extension, yet source melting dynamics remain constant. Only samples from the stratovolcanic edifice (in regional compression) preserve extensive zonation in olivine macrocrysts; these textures are generally absent from proximal small eruptive centres in transtension or extension. The zonation in olivines from stratovolcanic lavas yields timescales on the order of a few days to a few weeks, suggesting that even in environments which inhibit ascent, timescales between unrest and eruption of mafic magmas may be short. Significantly, high-resolution compositional profiles from olivine grains in the studied samples record evidence for post-eruptive growth and diffusion, highlighting the importance of careful interpretation of diffusion timescales from zoned minerals in more slowly cooled lavas when compared with tephra samples. [ABSTRACT FROM AUTHOR]

We investigate ten of the most recent tephra-fall deposits emplaced between ≤21 and 2 ka from the Pacheco stage of Misti volcano, Peru, to elucidate magma dynamics and explosive eruption triggers related to magma storage, recharge, and remobilization. Whole-rock, glass, and mineral textures and compositions indicate the presence of broadly felsic, intermediate, and mafic magmas in a chemically and thermally stratified magma storage system (Zones 1–3) that interact to differing extents prior to eruption. Intermediate magmas are defined by plagioclase + amphibole + two-pyroxenes + Fe-Ti oxides and phase equilibria indicate they formed at ~300 to 600 MPa and ~950°C to 1000°C. Intermediate magmas dominate the Pacheco stage and either erupted alone as hybridized magmas or mingled with minor volumes of cool felsic magmas (~800°C) in which only plagioclase + Fe-Ti oxides are stable. Felsic magmas do not exclusively comprise any tephra-fall deposit emplaced during the Pacheco stage but were remobilized by recharge and mixing with intermediate magmas in order to erupt. Furthermore, felsic-hosted amphibole cognate to the intermediate magmas are reacted despite the felsic magmas being water saturated, which suggests they are staged above the amphibole stability limit (≤200 MPa). The cryptic presence of mafic magmas is indicated by high-An plagioclase cores (An74–88), rare anhedral olivine (Fo77–80), and possibly high Mg# augite and amphibole (up to Mg# 84 and 77, respectively). The dearth of basalt to basaltic andesite melts recorded in erupted glasses and exclusivity of high-An plagioclase to crystal cores signals mafic magmas are staged deeper in the crust than the intermediate magmas. Periodic interactions between these magmas tracked via glass compositions and crystal exchange reveal an alternation between the production of mingled magmas and their eruption shortly after a recharge event, followed by a period of homogenization and eruption of hybridized magmas. As such, we identify magma recharge as a key mechanism by which half of the explosive eruptions were triggered in the Pacheco stage. A >100°C increase in Misti's fumarole temperatures from 1967 to 2018 coincident with changes in fumarolic gas compositions is consistent with degassing of a mafic recharge magma, signaling that Misti could produce similar explosive eruptions in the future. [ABSTRACT FROM AUTHOR]

A Bayesian linear regression to determine the bias in the Nafe–Drake relationship between compressional velocity and density provides an improved model for the density structure of Kīlauea volcano, Hawaiʻi. In previous work, we combined the results of seismic tomography with the Nafe–Drake relationship between compressional velocity and density to explain the large values of gravity disturbances overlying the summits and rift zones of the island's volcanoes. These results were used to determine mechanisms for gravitational instability of the island flanks. Here, we use laboratory measurements of the relationship of velocity and density for a wide range of Hawaiʻi island rocks as a prior in a Bayesian regression, with seismic tomography, to refine the 3-D density structure for Kīlauea volcano. This refined structure shows dense bodies (3220 kg m−3) between 5 and 8 km below sea level that underly regions of magma storage, found from geodetic and geophysical studies, beneath the summit and East Rift Zone of Kīlauea volcano. Above these bodies, density isosurfaces surround and cradle sources of pressure change determined from geodetic models, both at the summit and along the East Rift Zone. Continued subsidence of the summit following the 2018 eruption is aligned with a bowl-shaped density structure, formed primarily by density isosurfaces between 2800 and 2900 kg m−3 at 4–6 km depth. These surfaces underly the ∼3 km depth at which dyke injection initiates, are largely aseismic, and from their density values are inferred to contain high concentrations of olivine. Taken together, these density structures are consistent with an olivine-rich mush with variable porosity that increases in density with depth and provides a mechanism to form olivine cumulates both at the summit and along the rift zones. This structural framework for Kīlauea volcano is consistent with melt and mush transport occurring over a large range of depths to accommodate the growth and spreading of the volcano. [ABSTRACT FROM AUTHOR]

Imaging the intra-sediment magma chamber in the Damavand region, northern Iran, is beyond the resolution of the local seismic observations. Gravity anomalies can precisely image the lateral extension of magma reservoir. In order to provide vertical extension of the magma chamber, we apply the inversion of magnetic data with a higher sensitivity to shallow structure in comparison to the gravity data. More importantly, the knowledge of the magma chamber’s density allows the prediction of its mechanical behaviour, including the potential of eruption. As Damavand is estimated to be an active volcano, it is important to revisit the physical properties of the magma chamber to be able to evaluate the potential of eruption. Here, we apply the sparse norm inversion of the Bouguer gravity anomaly and magnetic data to model the uppermost crust beneath the Damavand volcano. The qualitative analysis of the Bouguer anomaly shows that the power of the spectrum remains almost unchanged by upward continuation using heights greater than 4 km. Thus, we conclude that the 4-km upward continued Bouguer anomaly represents the regional gravitational effects free from very shallow effects. The inversion of the magnetic anomaly, interestingly, shows a susceptibility structure, with susceptibility contrast of up to ~ +0.025 SI, in the same place as the density anomaly. This study proposes a 10- km wide magma chamber beneath Damavand from depth ~3 km to depth ~12 km. The resulting density structure is comparable with the obtained values from derived densities (using the thermodynamic mineral phase equilibrium) based on geochemical data and those from the conversion of seismic velocity to density. According to the geochemical data analysis, the lava is andesitic which is categorized among dense crustal rocks (2.8 g/cm³ ). But our modelling results shows a density contrast of maximum + 0.25 g/cm³ between the magma chamber and the surrounding sedimentary rocks (with a density of 2.45 g/cm³ ) above 5 km. Therefore, we can conclude that the shallow magma chamber, composed of dense andesite, is relatively warm and probably not completely consolidated. The high temperature of the magma chamber appears to be neutralized by the impact of the high density of andesite (naturally dense rocks) to result in a moderate negative anomaly in tomography (i.e., ∆Vs=~ -2 %). The magma chamber’s temperature might exceed 750-800 ºC which is still beyond the solidus-liquid transition temperature of 1100 ºC. Therefore, we can conclude that the magma is no liquid and is partially consolidated. [ABSTRACT FROM AUTHOR]

For active volcanoes, knowledge of the architecture of the plumbing system and the conditions of magma storage prior to an eruption are highly important, given their influence on the eruptive style and, thus, the management of future volcanic crises. Here, chlorine is used as a geobarometer for potassic alkaline magmas at the Campi Flegrei volcanic complex, revealing the shallowest depth of fluid-melt equilibration with respect to Cl. The results for representative fallout deposits of selected explosive eruptions show the existence of a multi-depth equilibration zone through time, including shallow magma storage. We describe evidence for the shallowest zone located at a depth equivalent to 65 MPa for the Agnano Monte Spina eruption (4482–4625 cal. yrs BP), at ~100 MPa for the Pomici Principali (11 915–12 158 cal. yrs BP), and the Astroni 6 (4098–4297 cal. yrs BP) eruptions, and close to 115 MPa for the last explosive eruption of Monte Nuovo (AD 1538). For comparison, the pressure estimated for a possible reservoir feeding the Cretaio eruption of Ischia island (AD 430), the only studied eruption on Ischia, is ~140 MPa. The pressure estimates for the two largest magnitude eruptions, the Campanian Ignimbrite (40 ka) and the Neapolitan Yellow Tuff (14.9 ka), are also discussed with respect to available magma withdrawal models. The pressures estimated using the Cl geobarometer for the magma leading to the fallout phases of these two eruptions provide evidence for a low-volume, shallow domain (~40 MPa) for the Plinian phase of the Campanian Ignimbrite eruption and a main, deeper reservoir (~130–165 MPa) for the Neapolitan Yellow Tuff eruption. The inferred shallowest equilibration pressures are interpreted here as corresponding to transitory, short-lived magma apophyses, whose eruption may have been facilitated by optimum tectonic stresses, rheological behavior of the crust, and efficiency of volatile exsolution. Alternatively, these magma apophyses may represent an evolved, crystal-rich ponded magma into which a volatile-rich magma ascending from depth was injected. The transient nature of such very shallow reservoirs is suggested by the short timescales inferred from diffusion modeling on crystals available in the literature for the studied Campi Flegrei eruptions. The influence of sulfur (S) on Cl solubility is assessed through Cl solubility modeling and applied to different eruptions. In addition, the pressure at which magmatic fluids and melts equilibrated with respect to Cl is shallower for the Campi Flegrei volcanic complex than the Somma-Vesuvio volcanic complex, erupting more homogeneous differentiated magma, of trachytic or phonolitic composition. This approach of using Cl to investigate the architecture of the plumbing system can be extended to all alkali-rich magma systems. [ABSTRACT FROM AUTHOR]

The size of deep-seated magma chambers is an important parameter for understanding pre-eruptive signals such as surface deformation. The constantly inflating Hekla volcano in Iceland has had relatively simple eruptive behaviour during the historical period. The eruptions start explosively with production of differentially evolved andesite magma to dacite, related to the length of the foregoing quiescence period, and ends with an emission of a basaltic andesite lava of uniform composition. The basaltic andesite is formed by fractional crystallisation from a deeper-seated basalt source in a steady-state manner. How fast such a differentiation mechanism operates is unknown. Measured Ra–Th radioactive disequilibrium in both the basalt and the basaltic andesite reveal a decrease from a 14% excess of 226Ra over 230Th to only 5% with magma differentiation. The decrease in 226Ra excess to 5% in the basaltic andesite of Hekla is shown to be controlled by plagioclase fractionation alone. Therefore, the magma differentiation time from basalt to intermediate magma beneath Mt. Hekla is significantly shorter than three centuries, the time needed to detect significant 226Ra-decay. Given the steady-state production of basaltic andesite magma and the estimated magma production rate, the volume of the basaltic andesite magma reservoir can be estimated as less than 2 km3. [ABSTRACT FROM AUTHOR]

Volcanic islands are often subject to flank instability, resulting from a combination of magmatic intrusions along rift zones and gravitational spreading causing extensional faulting at the surface. Here, we study the Koaʻe fault system (KFS), located south of the summit caldera of Kīlauea volcano in Hawaiʻi, one of the most active volcanoes on Earth, prone to active faulting, episodic dike intrusions, and flank instability. Two rift zones and the KFS are major structures controlling volcanic flank instability and magma propagation. Although several magmatic intrusions occurred over the KFS, the link between these faults, two nearby rift zones and the flank instability, is still poorly studied. To better characterize the KFS and its structural linkage with the surrounding fault and rift zones, we performed a detailed structural analysis of the extensional fault system, coupled with a helicopter photogrammetric survey, covering part of the south flank of Kīlauea. We generated a high-resolution DEM (~ 8 cm) and orthomosaic (~ 4 cm) to map the fracture field in detail. We also collected ~ 1000 ground structural measurements of extensional fractures during our three field missions (2019, 2022, and 2023). We observed many small, interconnected grabens, monoclines, rollover structures, and en-echelon fractures that were in part previously undocumented. We estimate the cumulative displacement rate across the KFS during the last 600 ~ 700 years and found a decrease toward the west of the horizontal component from 2 to 6 cm per year, consistent with GNSS data. Integrating morphology observations, fault mapping, and kinematic measurements, we propose a new kinematic model of the upper part of the Kīlauea's south flank, suggesting a clockwise rotation and a translation of a triangular wedge. This wedge is bordered by the extensional structures (ERZ, SWRZ, and the KFS), largely influenced by gravitational spreading. These findings illustrate a structural linkage between the two rift zones and the KFS, the latter being episodically affected by dike intrusions. [ABSTRACT FROM AUTHOR]

The Gulf of Guinea exhibits a continuous emission of narrow-band and long-period signals (16, 26 and 27 s) on teleseismic records, yet the underlying excitation mechanism remains unclear. This study establishes a connection between these tremors and the vibration of thin, decoupled crustal plates at unexplored volcanoes in the gulf. We first formulate the damped plate oscillation equation, by incorporating the vibration of the thin surface crustal plate and magma flow in the subsurface sill. The findings reveal that a fundamental-mode vibration with a period of several dozen seconds can be induced by a crustal plate that is less than 1.0 km thick but extends over tens of kilometres in both length and width, given a subsurface sill depth exceeding 10.0 cm. The thin plate hypothesis also allows for excitation of a few overtone modes, but such waves in higher frequencies diminish over long distances, leaving only the monotonous fundamental-mode vibration at teleseismic stations. The long duration of Guinea tremors at each recurrence is attributed to the presence of low viscosity basaltic magma, which influences the damping factor. Direct wave loads at the shallow gulf serve as the primary vibration source, accounting for seasonal variations and recurring patterns. Sporadic energy bursts may also occur due to large storms. Radiation patterns of Guinea tremors are linked to the geometric structure of the thin plate. Our theoretical estimates of tremor spectra closely align with observed data, confirming the model's accuracy in capturing reported Guinea tremor characteristics. This study provides valuable insights into the origins of very long-period tremors at continental volcanoes. [ABSTRACT FROM AUTHOR]

Determining the time spans of processes related to the assembly of eruptible magma at active volcanoes is fundamental to understand magma chamber dynamics and assess volcanic hazard. This information can be recorded in the chemical zoning of crystals. Nevertheless, this kind of study is still poorly employed for the active volcanoes of the Neapolitan area (Southern Italy), in particular, for Ischia island where the risk is extremely high and this information can provide the basis for probabilistic volcanic hazard assessment. For these reasons, we acquired chemical composition on clinopyroxene crystals erupted at Ischia during the Zaro eruption (6.6 ± 2.2 ka) and performed numerical simulations of the input of mafic magma into a trachytic reservoir, in order to investigate various aspects of pre-eruptive dynamics occurring at different timescales. This event emplaced a ~ 0.1 km3 lava complex, in which the main trachytic lava flows host abundant mafic to felsic enclaves. Previous petrological investigation suggested that mafic magma(s) mixed/mingled with a trachytic one, before the eruption. In this work, the clinopyroxene zoning patterns depict the growth of crystals in different magmatic environments, recording sequential changes occurred in the plumbing system before the eruption. The evolution of the plumbing system involved a hierarchy of timescales: a few hours for magma mingling caused by mafic recharge(s) and likely occurred multiple times over a decade during which a dominant magmatic environment was sustained before the eruption. Such timescales must be considered in volcanic hazard assessment at Ischia and similar active volcanoes in densely populated areas. [ABSTRACT FROM AUTHOR]

The thermal activity of a magmatic–hydrothermal system commonly changes at various stages of volcanic activity. Few studies have provided an entire picture of the thermal activity of such a system over an eruptive cycle, which is essential for understanding the subsurface heat transport process that culminates in an eruption. This study quantitatively evaluated a sequence of thermal activity associated with two phreatic eruptions in 2021 at Aso volcano. We estimated plume-laden heat discharge rates and corresponding H2O flux during 2020–2022 by using two simple methods. We then validated the estimated H2O flux by comparison with volcanic gas monitoring results. Our results showed that the heat discharge rate varied substantially throughout the eruptive cycle. During the pre-eruptive quiescent period (June 2020–May 2021), anomalously large heat discharge (300–800 MW) were observed that were likely due to enhanced magma convection degassing. During the run-up period (June–October 2021), there was no evident change in heat discharge (300–500 MW), but this was accompanied by simultaneous pressurization and heating of an underlying hydrothermal system. These signals imply progress of partial sealing of the hydrothermal system. In the co-eruptive period, the subsequent heat supply from a magmatic region resulted in additional pressurization, which led to the first eruption (October 14, 2021). The heat discharge rates peaked (2000–4000 MW) the day before the second eruption (October 19, 2021), which was accompanied by sustained pressurization of the magma chamber that eventually resulted in a more explosive eruption. In the post-eruptive period, enhanced heat discharge (~ 1000 MW) continued for four months, and finally returned to the background level of the quiescent period (< 300 MW) in early March 2022. Despite using simple models, we quantitatively tracked transient thermal activity and revealed the underlying heat transport processes throughout the Aso 2021 eruptive activity. [ABSTRACT FROM AUTHOR]

The geochemical monitoring of volcanic activity today relies largely on remote sensing, but the combination of this approach together with soil gas monitoring, using the appropriate parameters, is still not widely used. The main purpose of this study was to correlate data from crater gas emissions with flank emissions of soil gases at Mt. Etna volcano from June 2006 to December 2020. Crater SO2 fluxes were measured from fixed stations around the volcano using the DOAS technique and applying a modeled clear-sky spectrum. The SO2/HCl ratio in the crater plume was measured with the OP-FTIR technique from a transportable instrument, using the sun as an IR source. Soil CO2 efflux coupled with the 220Rn/222Rn activity ratio in soil gases (named SGDI) were measured at a fixed monitoring site on the east flank of Etna. All signals acquired were subject both to spectral analysis and to filtering of the periodic signals discovered. All filtered signals revealed changes that were nicely correlated both with other geophysical signals and with volcanic eruptions during the study period. Time lags between parameters were explained in terms of different modes of magma migration and storage inside the volcano before eruptions. A comprehensive dynamic degassing model is presented that allows for a better understanding of magma dynamics in an open-conduit volcano. [ABSTRACT FROM AUTHOR]

Open 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. Plain Language Summary: Open conduit volcanoes erupt with the highest frequency on Earth and release significantly more gas than magma. Their activity is characterized by almost continuous and small explosions that are intercalated by significantly larger explosive events. The relatively continuous explosive activity makes these volcanoes attractive both for scientific research and tourism, however, unexpectedly large explosions pose a significant risk to scientists and tourists. The excess gas these volcanoes release is CO2‐rich and its interaction with magma leads to the accumulation of gas‐rich and low‐density magma at depth, below higher density magma. We suggest that the gravitational destabilization of this deep and gas rich magma is ultimately responsible for the large explosions that punctuate the activity of open conduit volcanoes. The approach to this larger events can be traced with continuous gas monitoring, which therefore provides a unique opportunity to anticipate this otherwise unexpected explosions. Key Points: Continuous CO2 flushing can prime basaltic volcanic systems toward more energetic eruptionsAccumulation of gas‐rich and low‐density magma at depth results from CO2 flushing alone and no ad‐hoc physical barrier is requiredBuild‐up to larger explosions can be monitored through continuous gas monitoring at the surface [ABSTRACT FROM AUTHOR]

The nearly continuous volcanic eruption record at Mt. Etna dating back ~700 years provides an excellent opportunity to investigate the geochemical evolution of a highly active volcano. Of particular interest is elucidating the cause of selective enrichment in alkali elements (K and Rb) and 87Sr/86Sr observed in various episodes of past activity. More recently, this alkali enrichment trend started to manifest in the 17th century and accelerated after 1971, and was accompanied by an increase in the volume, frequency, and explosivity of eruptions. Explanations for this signature include recharge of alkali-enriched magmas and/or crustal contamination from the subvolcanic basement. This study quantitatively examines the role of crustal contamination in post-1971 Etnean magma compositions via hundreds of open-system phase equilibria and trace element calculations based upon whole-rock major oxides, trace elements, 87Sr/86Sr ratios, and mineral compositional data. Available pre-1971 petrochemical data are satisfactorily reproduced by fractional crystallization of a high whole-rock MgO (12–17 wt.%), Ni (135–285 ppm), and Cr (920–1330 ppm) parental magma composition that is documented in Etna's ~4-ka fall-stratified deposit. Observed post-1971 whole-rock and glass trends and phase equilibria are reproduced via modeled assimilation of a skarn and flysch mixture, lithologies that represent the uppermost 10 to 15 km of sedimentary rocks beneath Etna. Notably, models show that K2O (wt.%) and Rb (ppm) behave incompatibly during partial melting of skarn/flysch. Additionally, the observed elevation of 87Sr/86Sr in post-1971 samples is consistent with the addition of radiogenic Sr from wallrock partial melts. In best-fit models, which yield observed post-1971 K2O, Rb, and 87Sr/86Sr trends, ~17% anatectic melt is assimilated and there may be a subordinate stoped wallrock component of ≤2% (percentage is relative to the starting mass of pristine magma). Previous work has shown that metasomatized spinel lherzolite and garnet pyroxenite can be melted in different proportions to reproduce long- and short-term changes observed in Etna's geochemical products. We propose that the alkali enrichment signature observed after 1971 can be fully explained through the combination of mantle heterogeneity and crustal contamination. In particular, up to ~20% crustal input coupled with mantle heterogeneity of primitive melts explains the geochemical signals quite well. The influence of crustal contamination on post-1971 lavas is, in part, the result of frequent recharge of magmas that thermally primed the middle to upper crust and enhanced its partial melting. [ABSTRACT FROM AUTHOR]

The major- and trace-element compositions of amphiboles in andesite from Quaternary Yufu Volcano, northeastern Kyushu, Japan were analysed to investigate the generation processes of andesitic magma from Yufu Volcano. The amphiboles in andesite from Yufu volcano can be divided into two groups based on major-element composition: pargasite and magnesio-hornblende. To estimate temperature, pressure, and major- and trace-element compositions of melts in equilibrium with amphiboles, we used the recently proposed methods that can calculate temperature, pressure, major element compositions, and partition coefficients of trace-element between amphibole and melt using only the major-element compositions of amphibole. The estimated temperature, pressure, and major-element composition of melt in equilibrium with the amphibole phenocrysts indicate that each group crystallised under different conditions. These differences suggest that two magma chambers at different depths existed beneath Yufu Volcano and that the andesitic magma of Yufu Volcano was formed by mixing of the two magmas. The trace-element compositions of melts in equilibrium with the pargasite and magnesio-hornblende, estimated by applying the partition coefficients calculated from major-element compositions of amphibole to trace-element compositions of amphiboles, indicate magma derived from slab melt and the partial melting of crustal material, respectively. Because magma is a mixture of minerals and melt, we estimate the chemical compositional ranges of the two end-member magmas on the Y versus SiO2 diagram from the mixing relationship between amphibole and estimated melt, as well as phenocrysts of plagioclase, clinopyroxene, and orthopyroxene. The overlap of the estimated compositional range with the trend of whole-rock composition represents the chemical compositions of the end-members of magma mixing, yielding estimates of the mafic (SiO2 ≈ 45 wt %) and felsic (SiO2 ≈ 68 wt %) end-member magmas. Furthermore, we estimate the concentrations of other elements in the end-member magmas by substituting the estimated SiO2 concentrations of the magmas into linear regression equations between the whole-rock contents of other elements and SiO2. The trace-element compositions of the mafic and felsic end-member magmas, as estimated in this study, have similar features to those of gabbroids and Cretaceous granitic rocks, respectively, that are presumed to lie beneath Yufu Volcano. These similarities could be explained by the possibility that the compositions of the end-member magmas were influenced by basement rocks. [ABSTRACT FROM AUTHOR]

This paper provides information on the volumes of ejecta discharged by lateral and summit eruptions during the period from 1932 to 1989. We show that the volumes of ejecta from lateral eruptions were larger than those of the summit eruptions for the period. We examine subsequent variations in the character of eruptions. Lateral and summit eruptions had been occurring on the volcano from 1932 to 1989, and summit eruptions alone from 1990 to 2016. Changes occurred in 2016 with the resumption of lateral eruptions on the volcano's slopes. We note one characteristic feature, namely, all summit eruptions terminated in collapse events in the central crater. The collapse structures are characterized. The formation of the collapse structures is related to changes in the density and volume of magma as it was ascending to the ground surface. Data are provided relating to the summit eruptions of 2019–2021 and on the 2021 lateral eruption named after G.S. Gorshkov. New data are presented showing changes in the morphology of the summit crater during eruptions and collapses in the crater after the eruptions. [ABSTRACT FROM AUTHOR]

Lake Taupō, New Zealand, conceals a large caldera volcano which last erupted in 232 AD. A surveying experiment in 1979 which utilised the surface of the lake for detecting small movements of the lakebed soon became a programme to monitor vertical deformation by making regular observations at 22 fixed points – survey stations – around the lake. The programme continues today with 4 surveys being made annually. In this paper, we review the methodology and present the 42 year-long dataset in graphical forms, and draw conclusions on its implications in terms of volcanic and tectonic deformation. We also make the dataset available for download. The dataset confirms that vertical deformation is occurring in the lakebed. Long periods of slow deformation, dominated by subsidence totalling 140 mm in the Taupō Fault Belt at the northern end of the lake, and to a lesser degree at the southern end, are interrupted by uplift episodes now reaching 160 mm in the northeastern sector containing the most recent active vents. We suggest that the subsidence is primarily tectonic while inflation episodes are driven by upward migration of magma to shallow levels during periods of unrest. [ABSTRACT FROM AUTHOR]

The application of whole-rock compositions to trace magma evolution or crystal-melt equilibrium may be called into question when foreign crystals are incorporated into host magmas. To address this challenge, establishing the origins (orthocrysts, antecrysts, xenocrysts, etc.) of minerals with variable textures in magmatic plumbing systems is necessary. In this paper, we describe complex oscillatory zoning patterns of clinopyroxenes (Cpx) from the Miocene Majuangou (MJG) alkali basalts from West Qinling, China. Our detailed petrographic, mineralogical, and geochemical analyses unravel the origins of various Cpx crystals with distinct textural patterns, thereby providing information about magma storage, recharging and/or mixing, and transportation, as well as the reactions between crystals and melts/fluids. Based on textural patterns, Cpx may be divided into four types: normal (Type-1 Cpx), simple oscillatory (Type-2 Cpx), complex oscillatory (Type-3 Cpx), and grains that lack zoning (Type-4 Cpx, suggested to be orthocrysts). Through the textural characterization of Cpx, the comparison between different types of Cpx, and the relationships between Cpx major compositions from different lithologies, we concluded that Type-1–3 Cpx cores are antecrysts or xenocrysts with diverse origins: primitive magma (Type-1 Cpx cores), magma mush (Type-2 Cpx cores), and crustal granulite (Type-3 Cpx cores). The zoning patterns and the compositions of these Cpx crystals indicate at least three batches of magmatic activity, i.e., the Batch-1 low-Mg# magma (Mg#: 47.4–53.3), the Batch-2 primitive magma (Mg#: 57.2–64.5), and the Batch-3 low-Mg# host alkali magma (Mg#: 47.2–54.6). Cpx-melt thermobarometry demonstrates that at least two crustal magma reservoirs existed in the magma plumbing system at depths of 30.1 and 40.9 km. The antecrystic/xenocrystic Cpx cores were captured by, continued to grow in, and subsequently reacted with ascending K-rich melt/fluid. The spongy textures in Cpx cores/mantles are attributed to this reaction, which may be expressed as: Melt 1 (primitive or evolved) + K-rich melts/fluids + Cpx (CaMgSi2O6) = K-feldspar (KAlSi3O8) + ilmenite (FeTiO3) + Melt 2 (derivative). The products of this reaction (K-feldspar and ilmenite) filled the sieves in the spongy zones of Type-1–3 Cpx. This detailed investigation of compositional and textural features of Cpx antecrysts/xenocrysts suggests that the interactions between various interconnected magma reservoirs are widespread beneath the magmatic plumbing system. Our study emphasizes the importance of the incorporation of foreign crystals and the Cpx-melt/fluid reaction in magmatic plumbing system, which can significantly modify the whole-rock compositions and lead to the formation of spongy textures without the need for fractures and cracks in minerals. [ABSTRACT FROM AUTHOR]

Unrest began in July 2021 at Askja volcano in the Northern Volcanic Zone (NVZ) of Iceland. Its most recent eruption, in 1961, was predominantly effusive and produced ∼0.1 km3 lava field. The last plinian eruption at Askja occurred in 1875. Geodetic measurements between 1983 and 2021 detail subsidence of Askja, decaying in an exponential manner. At the end of July 2021, inflation was detected at Askja volcano, from GNSS observations and Sentinel‐1 interferograms. The inflationary episode can be divided into two periods from the onset of inflation until September 2023. An initial period until 20 September 2021 when geodetic models suggest transfer of magma (or magmatic fluids) from within the shallowest part of the magmatic system (comprising an inflating and deflating source), potentially involving silicic magma. A following period when one source of pressure increase at shallow depth can explain the observations. Plain Language Summary: Askja volcano, situated in the Northern Volcanic Zone in Iceland, has been quiet since its last eruption in 1961, with surface deformation measurements from 1983 to 2021 displaying a decaying subsidence signal within the Askja caldera. However, at the end of July 2021, the volcano began to inflate. This was detected on both GNSS and satellite observations. As of September 2023, ∼65 cm of uplift had been measured at GNSS station OLAC. Modeling of surface deformation measurements indicates that the inflation was triggered by upward migration of melt (or magmatic fluids). Key Points: At the end of July 2021, Askja volcano began to inflate—detected on both GNSS and satellite observations, ending 1983–2021 subsidenceGeodetic modeling indicates upward migration of magma, feeding a magma body at an inferred depth of 2.5–3.1 km under the main Askja calderaStart of unrest was associated with magma transfer within the upper part of the system, followed by possible additional influx from depth [ABSTRACT FROM AUTHOR]

The Great Dyke of Atla Regio (GDAR) is traced for ~3700 km on Venus, as a surface graben (narrow trough) interpreted to overlie a continuous laterally-emplaced underlying mafic dyke (vertical magma-filled crack). The GDAR belongs to a giant radiating dyke swarm associated with Ozza Mons (volcano), Atla Regio plume, and was fed from a magma reservoir ~600 km south of the Ozza Mons centre. A 50-degree counter-clockwise swing of the GDAR at 1200 km from the centre is consistent with a 1200 km radius for the underlying Ozza Mons plume head, and a stress link to the 10,000 km long Parga Chasmata rift system. Our discovery of the GDAR, should spur the search for additional long continuous single dykes on Venus (and Earth), with implications for estimating plume head size, locating buffered magma reservoirs, mapping regional stress variation at a geological instant, and revealing relative ages (through cross-cutting relationships) over regional-scale distances. In this study, the authors trace lateral magma movement in a ~ 3700 km long fracture on Venus. This Great Dyke of Atla Regio is the longest so-far traced on Venus (and in the solar system) and belongs to a giant radiating dyke swarm of the Ozza Mons volcano of the Atla Regio plume. [ABSTRACT FROM AUTHOR]

Using a multidisciplinary approach to understand the subsurface processes behind the formation of maar‐diatreme volcanoes is of growing interest. While geophysical characterization can visualize the diatreme and the feeding dike system beneath the volcano at a reasonable scale, such data are rare and generally unavailable. Stratigraphic‐controlled sampling and geochemical analysis of pyroclasts within the ejecta ring can, however, provide substantial information on dike evolution and the influence of the magmatic plumbing system on the growth of these volcanoes. Such investigation is presented here for the Barombi Mbo Maar (BMM), a complex maar of the Cameroon Volcanic Line (CVL) composed of a pile of tephra units linked to multiple explosive phases that were grouped into three eruptive episodes. Major and trace element compositions of lavas collected from the different eruptive units indicate that the erupted magmas at BMM consist mainly of basalt, trachybasalt and basanite, with Oceanic Island Basalts (OIB) and high μ (μ = 238U/204Pb) (HIMU) signatures. Compositional modelling suggests that partial melting occurred at different degrees in the garnet‐to‐spinel transition zone from one episode to another. The repetition of eruptions with big gaps between them, the presence of another large adjacent old maar crater next to the 2.5 km crater of the BMM, and the overall similarity in geochemical compositions from one eruption to another suggest a deep high‐productive zone in the mantle beneath the BMM. The latter productive zone was capable of generating magma batches episodically to fuel several individual monogenetic eruptions at the same location. [ABSTRACT FROM AUTHOR]

Geophysical data provide the chance to investigate a volcano's dynamics; considerable information can especially be gleaned on the stress and strain patterns accompanying the internal processes and the effect of magma ascent on the main structures triggering earthquakes. Here, we analysed in detail the seismicity recorded over the last two decades on Etna volcano (southern Italy), focusing on earthquakes distribution and focal mechanism clustering; the ground deformation pattern affecting the volcanic edifice with the inflation and deflation phases was also examined. Analysed data were compared in order to shed light on possible relationships with the volcanic activity and to better understand the internal dynamics of the volcano over time. Significant steps during or shortly before major eruptions in the seismic strain release and ground deformation temporal series highlight a straightforward relationship between seismicity occurring at shallow level, inflation/deflation and volcanism. Furthermore, at depths greater than 5–7 km, down to about 20 km, the orientation of the P- and T-axes clearly indicate the existence of a pressure source in the central part of the volcano. All the results underline that the stress field related to the volcano plumbing system interferes with the regional field, partly overriding it. [ABSTRACT FROM AUTHOR]

Lava fountains are a common manifestation of basaltic volcanism. While magma degassing plays a clear key role in their generation, the controls on their duration and intermittency are only partially understood, not least due to the challenges of measuring the most abundant gases, H2O and CO2. The 2021 Fagradalsfjall eruption in Iceland included a six-week episode of uncommonly periodic lava fountaining, featuring ~ 100–400 m high fountains lasting a few minutes followed by repose intervals of comparable duration. Exceptional conditions on 5 May 2021 permitted close-range (~300 m), highly time-resolved (every ~ 2 s) spectroscopic measurement of emitted gases during 16 fountain-repose cycles. The observed proportions of major and minor gas molecular species (including H2O, CO2, SO2, HCl, HF and CO) reveal a stage of CO2 degassing in the upper crust during magma ascent, followed by further gas-liquid separation at very shallow depths (~100 m). We explain the pulsatory lava fountaining as the result of pressure cycles within a shallow magma-filled cavity. The degassing at Fagradalsfjall and our explanatory model throw light on the wide spectrum of terrestrial lava fountaining and the subsurface cavities associated with basaltic vents. This study of volcanic gas chemistry during pulsatory lava fountaining at Fagradalsfjall volcano in Iceland reveals that the intermittency stems from pressure cycles and gas-melt separation within a shallow magma-filled cavity. [ABSTRACT FROM AUTHOR]

Changbaishan Volcano is considered one of the most hazardous active volcanic fields in northeastern Asia, and it has been the subject of increasing concern due to its unrest from 2002 to 2006. The physical conditions of magma chambers in the crust, particularly temperature and pressure, are crucial factors that determine the tempo and magnitude of volcanic eruptions, which are closely linked to potential hazards. However, the lithospheric thermal structure, which strongly influences rheological behaviour and melting, has been poorly studied. We conducted direct geothermal measurements and numerical forward modelling to address this issue to confirm whether a crustal magma chamber lies beneath the caldera as in earlier geophysical interpretations. We first reported four heat-flow data sets in Changbaishan Volcano. The research findings indicate that the heat flow value near the Tianchi caldera is remarkably high, at 270±16 mW m−2. As the distance from the caldera increases, the heat flow gradually decreases to a normal continental heat flow value of 75 mW m−2. 3-D transient heat simulations with a magma cooling model and a continuous magma supply model demonstrate the thermal effect of the magma chamber. The best-fitting model for the Tianchi magma system is an ellipsoidal magma chamber with a depth of 8–14 km, a 20 km east–west axis and a 70 km north–south axis, supplied with magma from the asthenosphere for 2 Myr. The high surface heat flow and crustal temperatures suggest that magma is active beneath Changbaishan Volcano within the middle-upper crust, and volcanic reactivation could occur. Thus, further research on the lithospheric thermal structure is necessary to understand the potential volcanic hazards associated with Changbaishan Volcano. [ABSTRACT FROM AUTHOR]

Copyright of Rudarsko-Geolosko-Naftni Zbornik is the property of Faculty of Mining, Geology & Petroleum Engineering and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Late Miocene-Quaternary adakitic rocks are widely distributed in northwestern Iran and are a key component of magmatism in the Turkish–Iranian Plateau of the collision zone between Eurasia and Arabia. Igneous rocks with adakitic affinity are distributed over large parts of northwestern Iran, eastern Turkey, and the Lesser Caucasus, and post-date a presumed slab break-off event at c. 10 Ma. Here, we present whole-rock Sr–Nd isotopic data for 4.5–0.149 Myr adakitic rocks of Sabalan volcano, NW Iran. These rocks are characterized by near-identical 87Sr/86Sr = 0.7044 to 0.7045 and ɛNd = + 2.24 to + 2.76. This is interpreted to indicate that Sabalan magmas were primarily generated by melting juvenile intrusions at the base of thickened lower crust and that assimilation of upper continental crustal rocks only played a minor role in their evolution relative to fractional crystallization. The delay between slab break-off beneath the Neo-Tethyan suture zone and magmatic activity at Sabalan is indicative for asthenospheric melting being triggered by small-scale mantle convection underneath the Turkish-Iranian Plateau. A preferred scenario is that the detached slab heated up in the mantle causing dehydration, and ascending fluids consequently lowered the viscosity of the mantle. Subsequently, the delaminated lower lithosphere was replaced by asthenospheric mantle which heated and partially melted lower crustal rocks to generate adakitic magmas. Collectively, these events occurred after a significant temporal hiatus that followed the earlier slab break-off event when remnants of subducted Neo-Tethyan lithosphere became detached. [ABSTRACT FROM AUTHOR]

We would like to thank Geist et al. (Contrib Miner Petrol, 2023) for the opportunity to further discuss the arguments presented in our paper "Long or short silicic magma residence time beneath Hekla volcano, Iceland?" (Sigmarsson et al. Contrib Miner Petrol 177:13, 2022). The disagreement centers around the origin of the silicic magmas at Hekla, namely whether it is by (1) fractional crystallization and a long crustal residence time before eruption or (2) partial melting of altered basaltic crust and short transfer time to the surface. We disagree with the arguments presented by Geist et al. (Contrib Miner Petrol, 2023) against the model for the origin of dacite at Hekla from dehydration melting of amphibolite, a model that still explains most if not all results obtained so far on the Hekla magma suite. [ABSTRACT FROM AUTHOR]

Recently reported major element, trace element, and volatile compositions of dacites erupted from Hekla Volcano, Iceland are inconsistent with an anatectic origin. In contrast, all of the data are consistent with dacitic melts being derived by modest amounts of crystallization of an anhydrous assemblage from an andesitic melt. Models that apply straightforward Gd, Yb, Th, and U concentrations, instead of ratios, to the problem show that the crystallization hypothesis remains valid, especially when uncertainties are propagated. [ABSTRACT FROM AUTHOR]

Kirishima volcano consists of more than 20 eruptive centers. Among them, Shinmoe-dake had magmatic eruptions in October 2017 and March 2018. Subsequently, another active cone, Iwo-yama, had phreatic eruptions in April 2018. These events were unique in that the 2018 eruption was the first effusion-dominated eruption of Shinmoe-dake and the first simultaneous activity of two cones of the Kirishima volcanic group ever documented. We report the detailed sequence of the events by combining areal photos, satellite images, and seismo-acoustic data analyses with the other published information. The seismo-acoustic data clarify the eruption onset and the transitions of the behaviors in three stages for each of the 2017 and 2018 eruptions. For both eruptions, we present regularly repeated tremors or 'drumbeat' earthquakes in the second stage, which interpret as gas separation from magma, leading to the ash-poor plume in the 2017 eruption or the effusive eruption in the 2018 event. We also propose that the 2017 and 2018 eruptions of Shinmoe-dake and the 2018 eruption of Iwo-yama are sequential events linked by the degassing of magma beneath Shinmoe-dake. An eruption like the 2017–2018 eruptions of Shinmoe-dake would leave few geological records and could be captured only by modern techniques. Although Shinmoe-dake has been believed to be an example of less-frequent eruptions, effusive eruptions like the 2018 case might have occurred more frequently in the past , but the following eruptions had obscured their records. The timelines summarized in this study will be useful in future studies of Kirishima volcanoes and world equivalences. [ABSTRACT FROM AUTHOR]

We obtain new images of the 3D resistivity structure of the Wudalianchi volcanic field (WVF). No low‐resistivity structure is found beneath Laoheishan and Huoshaoshan volcanoes in the WVF, which challenges the hypothesis of one or more magma reservoirs in the upper crust of this region. We observe a low‐resistivity body beneath Weishan volcano; however, estimated magma emplacement rates conflict with the observed volume of that body. Due to the relatively low geothermal gradient and surface heat flow in the area, magma would cool rapidly and could not be sustained for a long time unless new magma was regularly introduced into the system, but our observations do not support this scenario. We suggest that the magma in the WVF mainly originates from the mantle or the base of the crust with little crustal residence, and that the low‐resistivity body in the upper crust is likely due to saline aqueous free fluids. Plain Language Summary: The origin of magma, whether it is derived directly from the deep mantle or has previously resided in the crust, plays a crucial role in determining volcanic eruption mechanisms. The presence of a magma reservoir in the shallow crust often signifies an increased risk of volcanic eruptions. The Wudalianchi volcanic field (WVF) is situated in northeastern China. Earlier research suggested the existence of a magma reservoir beneath the Weishan volcano within the WVF. However, our study reveals an absence of magma reservoirs beneath the more recently erupted Laoheishan and Huoshaoshan volcanoes. This raises the question of why a magma reservoir is present beneath Weishan, which experienced an eruption 0.5 Ma ago, rather than beneath Laoheshan and Huoshaoshan, which erupted a mere 300 years ago. Furthermore, the WVF is characterized by an extremely low geothermal gradient, causing magma in the shallow crust to lose substantial heat and cool rapidly, unless replenished by a continuous influx of high‐temperature magma. Nevertheless, the volume of magma injection inferred from the low‐resistivity body beneath Weishan is inconsistent with numerical simulation results. Consequently, we suggest that the low‐resistivity body beneath Weishan volcano is instead a structure containing saline aqueous free fluids, not magma. Key Points: Electrical structure of the Wudalianchi volcanic field (WVF) is imagedSubsurface beneath Laoheishan and Huoshaoshan volcanoes is highly resistiveLow‐resistivity body beneath Weishan volcano is due to saline aqueous free fluids [ABSTRACT FROM AUTHOR]

A high‐precision earthquake catalog is made from continuous waveform data recorded at 360 portable seismographs deployed for one month in the Changbaishan‐Tianchi volcanic area based on an automatic earthquake detection and location workflow. Then we investigate detailed 3‐D structures of P and S wave velocities (Vp and Vs) and Vp/Vs ratio of the crust beneath the study region. Our results reveal a high‐Vp/Vs zone accompanied with active volcanic earthquakes beneath the Tianchi caldera at ∼5 km depth, which may reflect a magma chamber. The intense seismic activities within 2 km depth beneath the Tianchi caldera from 2002 to 2005 occurred directly above the magma chamber revealed by this study, which could be related to upward migration of high‐pressure volatiles originating from the magma chamber. A potential middle‐crustal magma chamber might provide continuous magmatic materials feeding the upper‐crust magma chamber and lead to ruptures of surrounding rocks, resulting in two swarms of volcanic seismicity during the observation period. We deem that the diverse eruption styles of the Changbaishan volcano are closely associated with the multilevel structure of magma chambers in the upper and middle crust. Plain Language Summary: It is difficult for the earthquake monitoring system to identify and locate microearthquakes by using traditional seismic detection methods due to their weak signals, leading to a relatively low quantity and quality of seismic data in a study region, which cannot meet the demands of high‐resolution seismic tomography. In this study, we detect and locate microseismic events beneath the Changbaishan‐Tianchi volcanic area using 1‐month waveform data recorded at our 360 portable seismographs deployed in the volcanic area. Then we investigate high‐resolution 3‐D crustal structures of P and S wave velocities (Vp, Vs) and Vp/Vs ratio beneath the study area. Our results reveal a possible crustal magma reservoir at ∼5 km depth beneath the Tianchi caldera, which is accompanied with volcanic earthquakes during our observation period. In addition, small earthquakes occurred actively at a depth of ∼2 km during 2002–2005, which are located right above the magma reservoir. We deem that the joint effect of magma reservoirs in the upper and middle crust could explain the active volcanic seismicity and distinctive eruption behaviors of the Changbaishan volcano. Key Points: A magma chamber at ∼5 km depth beneath the Tianchi caldera of the Changbaishan volcano (CBV) is revealedThe CBV unrest during 2002–2005 was related to reactivation of the crustal magma chamberDistinctive eruption behaviors of the volcano are related to multilevel magma chambers in the crust [ABSTRACT FROM AUTHOR]

Plagioclase phenocrysts in pre-historic andesites provide insight to the dynamics of magma formation and eruption at Taranaki volcano, New Zealand. The phenocryst population has a diversity of relic cores and a total in situ 87Sr/86Sr range of 0.70440–0.70486. Within-sample 87Sr/86Sr variations of 0.00018 to 0.00043 indicate that many phenocrysts are antecrysts and/or xenocrysts, derived from multiple crystal mush bodies. The Sr-isotopic differences in the phenocrysts of consecutive eruptions indicate that different magmas were tapped or formed on a centennial timescale. Most phenocrysts have multiple resorption/calcic regrowth zone(s) with elevated FeO* but invariant MgO zonation profiles. They likely record mafic melt inputs, and subsequent storage at elevated temperature caused re-equilibration of the Mg gradient. However, distinct rim types record different final pre/syn-eruptive magmatic conditions. Those in magmas erupted at 1030–1157 CE, 1290–1399 CE and 1780–1800 CE are characterised by resorption and calcic regrowth with sharp MgO and FeO* gradients. They record the entry of mafic melt into the system a few days or less before eruption based on Mg diffusion chronometry. In contrast, most phenocrysts erupted at 1755 CE, 1655 CE, and a few pre-1 ka events, have texturally uniform rims, compositionally consistent with closed-system crystallisation. This suggests alternating external and internal eruption triggers. Alternatively, the rate of magma reactivation via intrusion may dictate whether there was sufficient time for a mineralogical response to be recorded in part or all of the system. With respect to anticipating future eruptions, the plagioclase phenocrysts suggest multi-stage magma priming but rapid onset of the final trigger. [ABSTRACT FROM AUTHOR]

Between 38.5 ka cal BP and 32.4 ka cal BP, a dacitic Volcanic Explosivity Index 5 eruption at Misti volcano emplaced the Sacarosa tephra-fall deposit. Its biotite phenocrysts, fine grain size, scarce lithics, and abundant loose crystals characterize the deposit at locations sampled. The eruption's ~ 800 °C magma rose rapidly from ~ 10 km depth, culminating in a Plinian eruption which reached a mass eruption rate of 7.7 × 106–4.1 × 107 kg/s and emplaced about 3 km3 of tephra within tens of hours. The unit comprises two layers of subequal thickness separated by a diffuse contact with the upper distinguished by being slightly coarser and less well sorted than the lower. The deposit's coarser upper layer indicates either climactic conditions or a lesser degree of fragmentation during the latter half of the eruption. Strong winds distributed the deposit southwest of Misti, where it crops out over at least 800 km2 and drapes the present site of Arequipa with up to 100 cm of tephra. The Sacarosa deposit is the first among the Cayma stage deposits, a distinctive group of felsic, biotite-bearing units, to be carefully described and its eruption characterized. Several Cayma stage deposits were emplaced by voluminous explosive eruptions similar to the Sacarosa eruption, representing a ~ 8.9–15.5 ky interval of powerful eruptions. Such an explosive eruption today would threaten Arequipa's over 1,100,000 residents, many of whom live within the Sacarosa deposit's distribution. [ABSTRACT FROM AUTHOR]

Texture and particularly compositional variations of zoned minerals are important to understand trans-crustal magma plumbing systems of volcanoes. The Sipisupisu volcano in Sunda arc, northern Sumatra, occurs as a stratocone on the northwestern tip of the Toba Caldera and belongs to the large Toba volcanic system, which has produced two super-eruptions in the Quaternary. Compositionally zoned olivine, plagioclase and clinopyroxene grains, together with orthopyroxene reaction rims around olivine from basalts and basaltic andesites of the Sipisupisu volcano indicate open-system processes characterized by multiple episodes of magma recharge and mixing. Based on compositional variations of a variety of zoned olivine and clinopyroxene crystals, melts ascended from the large basaltic magma reservoir at the Moho have resided in at least three staging magma reservoirs within the crust. One magma reservoir is Si-saturated and located < 13 km deep in the upper crust, underlain by two Si-undersaturated ones at a depth from ~ 13 to 25 km. Fe–Mg diffusion modeling shows that olivine residence time from mixing of compositionally more evolved magmas to eruption varies from ~ 7 days to 9 months, whereas that from recharge of more primitive magmas to eruption is ~ 6 months to 5 years. This timescale is, on average, longer than olivine residence time of the majority of mafic arc volcanoes worldwide and may favor prolonged crystal fractionation. The multiple staging magma reservoirs and prolonged crystal residence time of the Sipisupisu volcano could potentially be linked to the large basaltic magma reservoir at depth. Finally, mixing of Si-saturated and Si-undersaturated magmas probably triggered Sipisupisu eruption. [ABSTRACT FROM AUTHOR]

The lithosphere--asthenosphere boundary (LAB), which can be seismically detected, stabilizes plate tectonics. Several conflicting hypotheses have been proposed as the causes of LAB discontinuity, such as the contribution of hydrated minerals, mineral anisotropy, and partial melts. The petit-spot melts ascending from the asthenosphere, owing to subducting plate flexures, support the partial melting at the LAB. Here, we observed the lava outcrops of six monogenetic volcanoes formed by petit-spot volcanism in the western Pacific. Thereafter, we determined the 40Ar/39Ar ages, major and trace element compositions, and Sr, Nd, and Pb isotopic ratios of the petit-spot basalts. The 40Ar/39Ar ages of two monogenetic volcanoes were ca. 2.6 Ma (million years ago) and ca. 0 Ma, respectively. The isotopic compositions of the western Pacific petit-spot basalts suggest their geochemically similar melting sources. They were likely derived from a mixture of high-μ (HIMU) mantle-like and enriched mantle (EM) -1-like components related to carbonatitic/carbonated materials and recycled crustal components. A mass balance-based melting model implied that the characteristic trace element composition (i.e., Zr, Hf, and Ti depletions) of the western Pacific petit-spot magmas could be explained by the partial melting of garnet lherzolite with a small degree of carbonatite melt flux with crustal components. This result confirms the involvement of carbonatite melt and recycled crust in the source of petit-spot melts and provides an implication for the genesis of tectonic-induced volcanism with similar geochemical signatures to those of petit-spots. [ABSTRACT FROM AUTHOR]

The lithosphere–asthenosphere boundary (LAB), which can be seismically detected, stabilizes plate tectonics. Several conflicting hypotheses have been proposed as the causes of LAB discontinuity, such as the contribution of hydrated minerals, mineral anisotropy, and partial melts. The petit-spot melts ascending from the asthenosphere, owing to subducting plate flexures, support the partial melting at the LAB. Here, we observed the lava outcrops of six monogenetic volcanoes formed by petit-spot volcanism in the western Pacific. Thereafter, we determined the 40Ar/39Ar ages, major and trace element compositions, and Sr, Nd, and Pb isotopic ratios of the petit-spot basalts. The 40Ar/39Ar ages of two monogenetic volcanoes were ca. 2.6 Ma (million years ago) and ca. 0 Ma, respectively. The isotopic compositions of the western Pacific petit-spot basalts suggest their geochemically similar melting sources. They were likely derived from a mixture of high-μ (HIMU) mantle-like and enriched mantle (EM) -1-like components related to carbonatitic/carbonated materials and recycled crustal components. A mass balance-based melting model implied that the characteristic trace element composition (i.e., Zr, Hf, and Ti depletions) of the western Pacific petit-spot magmas could be explained by the partial melting of garnet lherzolite with a small degree of carbonatite melt flux with crustal components. This result confirms the involvement of carbonatite melt and recycled crust in the source of petit-spot melts and provides an implication for the genesis of tectonic-induced volcanism with similar geochemical signatures to those of petit-spots. [ABSTRACT FROM AUTHOR]

The 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. Plain Language Summary: To improve eruption forecasting, it is important to put limits of the volume of eruptible magma stored beneath a volcano. There is evidence that regions of magma storage can be approximated by multiple "reservoirs" of magma beneath a single volcano. If magma can be transferred between these reservoirs, more magma may be available when an eruption occurs. We use models of ground displacements measured by satellite over 3 years to identify two reservoirs of magma at Ambrym volcano, located in Vanuatu, from 2019 to 2022. We then estimate how much magma entered or left these respective reservoirs during this time. Finally, we conclude that magma cannot be transferred between the two reservoirs efficiently, which controls the total amount of magma that can erupt at Ambrym. Key Points: At least two deformation sources were active at Ambrym volcano between the 2018 and 2022 eruptionsLava lake drainage caused closed‐system activity in 2019, resulting in uplift from ≥0.122 km3/year of magma inflow into a tilted reservoirParameter space exploration of a theoretical model indicates that the two distinct reservoirs are not efficiently hydraulically connected [ABSTRACT FROM AUTHOR]

SUMMARY: Unzen volcano, located on Shimabara Peninsula, Nagasaki, Japan, is an active volcano that has been intensively monitored since 1989, one year before the most recent eruption in 1990–1995. Previous earthquake and surface deformation studies have revealed that magma is transported obliquely from a magma reservoir beneath Tachibana Bay, to the west of Shimabara Peninsula. Here, we conduct broad-band magnetotelluric (MT) surveys at 99 sites around Shimabara Peninsula to investigate the crustal structure beneath Unzen volcano that is related to magma migration. A 3-D resistivity model that is constructed from 25 broad-band MT sites and 45 telluric sites shows a broad high-resistivity zone beneath Shimabara Peninsula and low-resistivity zones to the west and east of the peninsula. An unexpected observation is the spatial alignment of the high-resistivity zone with a seismic low-velocity zone (LVZ) at 3–15 km depth.  Quantitative analysis indicates this high-resistivity zone contains < 0.7 per cent melt under the assumption that the melt is stored in a good porosity network, while < 11 per cent melt in relatively poor pore network. We infer this high-resistivity, LVZ to be a highly crystallized mush zone (HCMZ) with low permeability. The hypocentres and pressure sources of the 1990–1995 eruption are distributed along the boundary between the high- and low-resistivity zones beneath the western part of the peninsula. We therefore conclude that the magma migrated along a structural boundary that possessed a relatively high permeability. Previous studies have suggested that eruptible magma is usually transported vertically upward through the centre of the mush zone, whereas the present results reveal that magma can be transported along the upper boundary of an HCMZ. [ABSTRACT FROM AUTHOR]

The Azufre volcano (21°47′S; 68°14′W) is emplaced above the western boundary of the Altiplano Puna Magma Body. The youngest lava of the Azufre volcano (50–331 ka) was selected for detailed studies because of the relevant petrological information that could emerge to understand its reservoir as a potential heat source of the neighboring Cerro Pabellón geothermal system. The studied lava corresponds to an andesite-dacite (61–63 SiO2 wt%) with phenocrysts of plagioclase, amphibole (Group 1), biotite, pyroxenes, quartz, and olivine. The lava also contains aphanitic enclaves (58–60 SiO2 wt%), whose groundmass have the same mineralogy of lava samples groundmass consisting of amphibole microphenocryts (Group 2) and microlites of plagioclase, pyroxenes, and Fe–Ti oxides. Disequilibrium textures are commonly observed in the studied samples such as partially resorbed plagioclase phenocrysts, amphibole breakdown, and reverse zoning in pyroxene and plagioclase phenocrysts. Thermobarometry calculations indicated pressures of ~ 2 kbar for Group 1 amphiboles, temperatures of 905–1097 °C for Mg-rich pyroxene phenocrysts and Ca-rich plagioclase (An ≥ 66), and near-solidus temperatures of 712–788 °C for Group 1 amphibole-plagioclase pairs. Group 1 amphiboles also indicate crystallization from an evolved liquid (63–79 wt% of SiO2). Oxidation conditions of QFM + 0.9–2.5 log units were recorded in amphibole and Fe–Ti oxides. Rhyolite-MELTS models reproduce the composition of the high-temperature phases from a melt composition (andesitic enclave) at similar P–T-fO2. The arrival of new hot andesite magma into a crystal-rich shallow reservoir, thermally and compositionally zoned, would have triggered the studied eruption. Diffusion models in Fe–Ti oxides microlites indicated cooling temperatures of 742–866 °C during the sub-aerial emplacement of the lava. [ABSTRACT FROM AUTHOR]

The basaltic caldera volcanoes in the Galápagos Islands characteristically erupt lavas via summit circumferential and radial flank dike intrusions, but the underlying magma plumbing systems remain enigmatic. Here, we document surface displacements of Wolf volcano using interferometric synthetic aperture radar (InSAR) data from 2015 to 2022. We show that Wolf volcano experienced 6‐years of continuous inflation after the 2015 eruption, followed by a shallow flank eruption in January 2022. The deformation is modeled with a vertical caldera ring‐fault and a radial dike on the southeast flank. The ring‐fault underwent opening and reverse faulting during the inflation period, and closure and normal faulting during and after the eruption. Stress interactions between the ring‐fault and the flank dike suggest the asymmetric opening of the ring‐fault promotes flank eruptions. The best‐fit deformation model differs from the previous models and offers an alternative view of how magma is fed into radial dikes during flank eruptions. Plain Language Summary: The 2022 eruption at Wolf volcano in the Galápagos Islands started from a radial fissure on the SE flank, which was different from the 2015 eruption that mainly occurred inside the caldera. We track the sequence of ground deformation before, during and after the 2022 eruption using InSAR data. We find that Wolf volcano experienced continuous inflation before the 2022 eruption, followed by deflation during and after the eruption, resulting in complex ground deformation within and beyond the caldera and on the SE flank. We propose a model with ring fault slip and opening/closure to explain the observed ground deformation within the caldera. We find that opening and reverse slip of the ring fault can explain the inflation before the 2022 eruption, while closure and normal slip caused the deflation period during and after the eruption. The asymmetric opening of the ring fault changed the stresses in the volcano and promoted the occurrence of the flank eruption. Our model is different from the dike rotation model that was proposed for Fernandina volcano and may allow for improved forecasts of the location and style of future eruptions on the Galapagos Islands. Key Points: Interferometric synthetic aperture radar data are used to study ground deformation before, during and after the 2022 Wolf eruptionA ring fault opening/closure model fits ground deformation of the caldera wellAsymmetric opening of the ring fault promotes flank eruptions at Wolf volcano [ABSTRACT FROM AUTHOR]

The Central Series of the Paleocene mafic to ultramafic Rum complex is the youngest of the three main cumulate series that make up the layered igneous complex of this famed Scottish island. The Central Series lies along the Long Loch Fault and provides insights into the feeder system to the Rum intrusion at an erosion level of about two kilometres below the former land surface. Much of the Central Series consists of a mélange of steep sided bodies of magmatic breccias that stretch along the Long Loch Fault (LLF) in a relatively narrow zone and is composed of blocks and clasts of all sizes derived largely from break‐up of the former conduit walls. Repeated movements of the LLF are thought to have been responsible for opening and closing of the magma conduit, resulting in repeated replenishment events, each of which gave rise to new cumulate formation within the Central Series and the bordering Eastern and Western Layered Series, which crystallized under relatively tranquil conditions. The Central Series probably acted as the feeder zone supplying the neighbouring layered series. The more complete of these is the Eastern Layered Series in which 16 conformable units can be distinguished. Others are presumed either to lie unseen at depth or to have been stripped by erosion. The Central Series, although often neglected because of its relative inaccessibility and complexity, formed from successive magma replenishments alternating with large‐volume side‐wall collapses of previously deposited cumulate material. It could thus be thought of as representing the 'pulsing heart' to the Rum volcano and deserves to be regarded as a site of major volcanological and petrological importance. Here we present a summary of some 60 years of investigation into the cumulate rocks of the Central Series. [ABSTRACT FROM AUTHOR]

We investigated magma storage conditions prior to the 2011 and 2018 eruptions at Kirishima Shinmoe-dake volcano in Japan based on the relationship between geodetic volume change of magma chamber and erupted magma volume. We derived an analytical expression for the ratio of the erupted magma volume to the geodetic volume change ("volume ratio"), which was formulated as a function of parameters related to the magma storage conditions. This expression shows that the volume ratio is strongly dependent on the effective compressibility of the magma chamber, which in turn depends on the rigidity of surrounding host rocks and shape of the chamber. For the Shinmoe-dake eruptions, the magnitude of the volume change (i.e., deflation) of a spherical magma chamber associated with lava effusion was estimated based on geodetic observations. The erupted magma volume was estimated from a SAR image analysis of the lava accumulation inside the summit crater. Based on these observations, we estimated that the volume ratio in 2011 and 2018 was 2.69 and 2.33, respectively. Substituting the estimated volume ratio into the analytical expression revealed that the observed geodetic data and volume ratio can be explained only when the magma chamber, which was assumed to be spherical, is filled with bubble-free magma. This result suggests that efficient gas segregation from the chamber occurred prior to the eruptions. Our results indicate that combining multi-observation data based on the volume ratio provides valuable information about the magma storage process, such as the behavior of the gas phase in the magma chamber. [ABSTRACT FROM AUTHOR]

Shallow magma transfer is difficult to detect at poorly monitored volcanoes. Magma transfer before the last 1538 eruption at Campi Flegrei caldera (Italy) was exceptionally tracked using historical, archeological, and geological data. Here, we extend that data set to 1650 to uncover any magma transfer during post‐eruptive subsidence. Results show two post‐eruptive subsidence phases, separated by a previously undocumented uplift during 1540–1582. Uplift highlights the pressurization of the central (∼3.5 km depth) and peripheral (∼1 km depth) pre‐eruptive sources, suggesting an aborted eruption. The subsidence events mainly require the depressurization of the central source and pressurization of a deeper magmatic layer (∼8 km depth). Therefore, despite the overall post‐eruptive deflation, after 1538 the deeper reservoir experienced continuous magma supply, with magma almost erupting between 1540 and 1582, challenging the common assumption of post‐eruptive deflation. This underlies the importance of monitoring the deeper magmatic systems, also after eruptions, to properly assess their eruptive potential. Plain Language Summary: Today, volcanic activity is monitored by satellite and ground networks. However, very little is known about the pre‐ and post‐eruptive behavior of volcanoes before the instrumental era. Here, we present a unique set of archeological records related to the elevation changes at Campi Flegrei caldera (Italy) during 1515–1650, mainly focusing on its behavior after the last 1538 eruption at Monte Nuovo. After the eruption, subsidence occurred from 1538 to 1540, followed, from 1540 to 1582, by a previously unreported uplift; the latter was followed by renewed subsidence until 1650 at least. Modeling the sources responsible for the deformation, we find that, despite the depressurization of shallow sources, after 1538, a deeper magmatic layer (∼8 km depth) experienced continuous magma supply, with magma almost erupting between 1540 and 1582. This questions the common notion of post‐1538 deflation at Campi Flegrei, where the depressurization of the shallower sources masks the pressurization of the deeper magmatic system for more than a century. Key Points: 20 archeological sites at Campi Flegrei show two post‐1538 eruption subsidence phases, separated by an undocumented uplift in 1540–1582During 1540–1582, a central sill‐like source (∼3.5 km depth) transfers magma below Monte Nuovo, representing an aborted eruptionFrom 1538 to 1650 a deeper magmatic layer (∼8 km depth) experienced continuous magma supply, also during caldera subsidence [ABSTRACT FROM AUTHOR]

Volcanic seismicity provides essential insights into the behavior of an active volcano across multiple time scales. However, to understand how magma moves as the eruption cycle develops, better knowledge of the geometry and physical properties of the magma plumbing system is required. In this study, using full‐wave ambient noise tomography, we image the three‐dimensional (3‐D) crustal shear‐wave velocity structure below Great Sitkin Volcano in the central Aleutian Arc. The velocity model reveals two low‐velocity anomalies correlating with the migration of volcanic seismicity. With a bulk melt fraction of about 2.5%–9%, these low‐velocity anomalies are interpreted as mushy magma reservoirs. We propose a six‐stage eruption cycle to explain the migration of seismicity and the alternating eruption of the two reservoirs with different recharging histories. These findings have broad implications for the dynamics of magma plumbing systems and the structural control of eruption behaviors. Plain Language Summary: Understanding magma accumulation and transport systems below active volcanoes is essential for predicting eruption behavior and assessing the potential hazards. The distribution of earthquakes can partly be used to infer the development of magmatic activity at different times. However, to understand how magma moves at different stages of an eruption cycle, better knowledge of what the magma plumbing system looks like is necessary. In this study, we use an advanced seismic imaging method to construct the 3‐D crustal shear‐wave velocity structure below Great Sitkin Volcano in the central Aleutian Arc. The velocity model reveals two crustal magma reservoirs, which correlate with the migration of seismicity. We propose a six‐stage eruption cycle to explain the evolution of seismicity in space and time across the island and the alternating eruption of two reservoirs. The findings in this study help to understand better the control of eruption behaviors by the underlying magma plumbing system at active volcanoes. Key Points: The pre‐ and co‐eruptive seismicity below Great Sitkin Volcano, Alaska, shows a spatiotemporal migrationA new 3‐D shear‐wave velocity model reveals two crustal low‐velocity anomalies that correlate with the migrating seismicityWe propose a six‐stage eruption cycle involving two magma reservoirs to explain the long‐term and short‐term seismicity patterns [ABSTRACT FROM AUTHOR]

Lava domes exhibit highly unpredictable and hazardous behavior, which is why imaging their morphological evolution to decipher the underlying governing mechanisms remains a major challenge. Using high-resolution satellite radar imagery enhanced with deep-learning, we image the repetitive dome construction-subsidence cycles at Popocatépetl volcano (Mexico) with very high temporal and spatial resolution. We show that these cycles resemble gas-driven rise and fall of the upper magma column, where buoyant bubble-rich magma is extruded from the conduit (in ~hours-days), and successively drained back (in ~days-months) as magma degasses and crystallizes. These cycles are superimposed on a progressive decadal crater deepening, accompanied by heat and gas flux decrease, which could be partially explained by gas depletion within the magma plumbing system. Results reinforce the idea that gas retention and escape from the magma column play a key role in the short- and long-term morphological evolution of low-viscosity lava domes and their associated hazards. Satellite imagery enhanced with deep learning sheds light on the mechanisms driving lava dome construction-destruction cycles. Results suggest that gas retention and escape from the magma system control the dome and crater morphological evolution. [ABSTRACT FROM AUTHOR]

Kharchinsky and Zarechny volcanoes and the Kharchinsky Lake zone of monogenic cones are unique eruptive centers of magnesian lavas located above the northern margin of the Pacific plate subducting beneath Kamchatka. This paper presents new geochemical data on the composition of rocks (55 samples) and minerals (over 900 analyses of olivine, pyroxenes, amphibole, and plagioclase) of these centers analyzed by XRF and LA-ICP-MS (rocks) and electron microprobe (minerals). Most of the studied rocks are magnesian (Mg# = 60–75 mol %) medium-K basalts and basaltic andesites. Moderate-magnesian (Mg# = 52–59 mol %) basaltic andesites are present among the monogenic cones of Kharchinsky Lake. The rare rock varieties include the high-K basalts–basaltic andesites of dikes in the center of Kharchinsky volcano and the magnesian andesites (Mg# = 58–61 mol %) of the extrusions of Zarechnу volcano. The distribution of trace-element contents in these samples demonstrates enrichment in large-ion lithophile elements and light REEs at depletion in high field strength elements and heavy REEs, as is typical of arc rocks. The high-K basalts and basaltic andesites show anomalous enrichment in Ba (>1000 ppm), Th (>3.8 ppm), U (>1.8 ppm), Sr (> 800 ppm, Sr/Y > 50), and light REE (La > 20 ppm), and their compositions are close to those of low-Si adakites. The basalts and basaltic andesites contain phenocrysts of high-Mg olivine (up to Fo92.6) and clinopyroxene (Mg # up to 91 mol %). The rocks show petrographic and geochemical evidence of fractional crystallization, along with the processes of mineral accumulation and magma mixing. Some of the olivine phenocrysts show high NiO contents (up to 5000 ppm) and an elevated Fe/Mn ratio (up to 80), which were interpreted as evidence of the participation of a pyroxenite source in the magma generation processes. The use of the Ca/Fe and Ni/Mg ratios allowed us to distinguish the composition fields and evolution trends of olivine associated with different sources: peridotite and pyroxenite, which were formed by a reaction between mantle-wedge peridotites and high-Si melts of the subducted oceanic crust. The new data are consistent with other lines of evidence of melting of the subducted Pacific plate edge beneath the northern part of the Central Kamchatka Depression at the Kurile–Kamchatka and Aleutian subduction zone junction and testify to a significant heterogeneity of the mantle in this area. [ABSTRACT FROM AUTHOR]