Southern Volcanic Zone Research Articles

Targeting deeply-sourced seeps along the Central Volcanic Zone

At convergent margins, plates collide producing a subduction process. When an oceanic plate collides with a continental plate, the denser (i.e., oceanic) plate subducts beneath the less dense (continental) plate. This process results in the transportation of carbon and other volatiles into Earth’s deep interior and is counterbalanced by volcanic outgassing. Sampling deeply-sourced seeps and fumaroles throughout a convergent margin allows us to assess the processes that control the inventory of volatiles and their interaction with the deep subsurface microbial communities. The Andean Convergent Margin is volcanically active in four distinct zones: the Northern Volcanic Zone, the Central Volcanic Zone, the Southern Volcanic Zone and the Austral Volcanic Zone, which are each characterised by significantly different subduction parameters like crustal thickness, age of subduction and subduction angle. These differences can change subduction dynamics along the convergent margin, possibly influencing the recycling efficiency of carbon and volatiles and its interaction with the subsurface microbial communities. We carried out a scientific expedition, sampling along a ~800 km convergent margin segment of the Andean Convergent Margin in the Central Volcanic Zone of northern Chile, between 17 °S and 24 °S, sampling fluids, gases and sediments, in an effort to understand interactions between microbiology, deeply-sourced fluids, the crust, and tectonic parameters. We collected samples from 38 different sites, representing a wide diversity of seep types in different geologic contexts. Here we report the field protocols and the descriptions of the sites and samples collected.

Syn and post-glacial volcanic activity in the western retroarc area of the andean southern volcanic zone: Overo volcano

Located in a poorly studied zone of the rear arc in the Transitional Southern Volcanic Zone, the Overo volcano presents a complete record of syn-glacial, inter-glacial, and post-glacial activity, marked by significant changes in eruptive mechanisms over time. Following pre-glacial units, an effusive syn-glacial basaltic phase was succeeded by an extensive inter-glacial diking that fed andesitic monogenetic fields. Subsequently, an eruptive pulse, possibly coupled with the destabilization of the main edifice, resulted in a block and ash deposit channeled into a glacial valley towards the southwest, later covered by extensive basaltic lava flows. Then, different types of pyroclastic flows and glacial deposits covered these earlier units. Finally, these younger glacial deposits were overlain by a localized postglacial pyroclastic deposit and widespread lava flows, which ultimately reconstructed the deeply eroded glacial morphology of the Overo volcano. This volcanic record reveals the young and complex activity of this poorly known western rear-arc volcanic center in the Transitional Southern Volcanic Zone of the Andes, whose evolution is more akin to the arc front volcanoes to the west than to their eastern within-plate counterparts. We propose that this magmatic evolution was influenced by the Pleistocene-Holocene Last Glacial Maximum activity (LGM), beginning with an initial phase of volcanic activity encapsulated within glacial conditions. Glacial retreat then triggered isostatic uplift and crustal relaxation, resulting in dike formation, flank instability, and syn-eruptive collapse of the volcanic edifice. Subsequently, a lesser glacial pulse partially eroded the volcano morphology and produced the currently preserved glacial deposits. Finally, a latest and limited post-glacial pyroclastic eruption occurred simultaneously with the formation of an extensive monogenetic basaltic field, shaping the final morphology of the volcano.

Upper‐flow‐regime deposits related to glacio‐volcanic interactions in Patagonia: Insights from the Pleistocene record in Southern Andes

AbstractSupercritical flows and their associated upper‐flow‐regime bedforms have been increasingly recognized in several contrasting depositional environments around the world, including volcano‐sedimentary settings. In recent years, there have been significant advances in flow modelling, architecture and depositional dynamics in upper‐flow‐regime bedforms related to glacio‐volcanic interaction, however most of these works come from specific volcanic landscapes. Examples of other regions with different tectonic and volcanic settings are needed in order to gain a more global perspective regarding the generation and preservation of upper‐flow‐regime glacio‐volcanic bedforms. This study presents a detailed analysis of Pleistocene volcaniclastic deposits accumulated in alluvial to lacustrine deltaic environments in the Southern Andes within the Southern Volcanic Zone. Six facies were defined and related to distinctive bedforms by examining sedimentary features, including textural characteristics, sedimentary structures, architecture and composition. Among these facies, five are distinguished by deposition from supercritical flows and upper‐flow‐regime bedforms, such as cyclic steps, chutes and pools, and antidunes. The stratigraphic succession reveals a predominance of repeated supercritical flows, which appear to be influenced by the availability of detritus for removal along with sudden water discharges. Compositional analyses indicate that explosive volcanism was the primary source of detritus, occurring synchronously with sedimentation. Additionally, the widespread occurrence of Pleistocene glacial conditions in this region suggests that volcanic‐induced melting ice or snow may have been an important source of water supply. Regarding the Southern Andes, the presence of local depressions, such as volcano‐tectonic calderas plays an essential role in the generation and preservation of the upper‐flow‐regime bedforms related to glacio‐volcanism. Furthermore, the importance of carrying out detailed and multidisciplinary studies seems to be critical for the recognition of these deposits in the record of the Andes.

Gas geothermometry, soil CO2 degassing, and heat release estimation to assess the geothermal potential of the Alpehue Hydrothermal Field (Sollipulli volcano, Southern Chile)

The Alpehue Hydrothermal Field (AHF) near the Sollipulli Volcano in the Southern Volcanic Zone of Chile shows promise as a significant geothermal resource. A comprehensive geothermal exploration survey was conducted, including the evaluation of hydrothermal gases, geothermometer calculations, and CO2 flux measurements, to assess the AHF's geothermal potential. Our results indicate that the hydrothermal gasses at the AHF primarily originate from primitive, mantle-derived sources, with some contribution from crustal sediments. Two different CO2 populations of fluxes were identified. One corresponds to the background emission related to the soil biological activity (mean ∼7.7 g·m−2·d−1), and the other, much more significant, emanates from an endogenous source related to the Alpehue hydrothermal reservoir (mean ∼461 g·m−2·d−1). Reservoir temperatures were calculated using gas geothermometry yielding average temperatures of 249 °C. The calculated heat flow rate of the AHF is approximately 3.3 MW and the heat flux corresponds to 156 thermal MW⋅km−2, which could be considered a medium geothermal potential comparable to other systems worldwide. Although further studies are needed to fully address its exploitability, this study presents favorable characteristics of the AHF that make it a promising avenue for further exploration.

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.

Effect of crustal stress state on magmatic stalling and ascent: case study from Puyehue-Cordón Caulle, Chile

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.

Glacial and volcanic landforms in the upper catchment of Claro River, Central Chile (35.5°S): A Late Quaternary geomorphological case study

AbstractThe upper catchment of the Claro river is located within the Transitional Southern Volcanic Zone in the Maule Region of Chile (35.5°S), at the foothills of Manantial Pelado volcano, a Pleistocene stratovolcano whose morphology suggests the occurrence of important glacial processes in the area. The observed glacial landforms can be sequenced and associated with episodes of volcanic activity. Different units of volcanic and glacial origin were identified based on field observations, detailed mapping, and digital elevation model and satellite image analysis. The contact relations between them made it possible to discern four events that can be linked to Marine Isotope Stages 8, 6, 4‐2, and to the Neoglacial. We propose seven stages in the geomorphological evolution from the late Pleistocene to the late Holocene, where there is coexistence between volcanic and glacial episodes, such as the emplacement of the Manantial Pelado volcano overlying the Abanico and Cola de Zorro formations; a subsequent debris avalanche; the enlargement of the Del Indio valley during the Penultimate Glaciation; the construction of Cerro Redondo (a minor eruptive center at the head of Del Indio valley); ongoing fluvial incision; and Neoglacial advances. The relative chronology proposed in this work contributes to clarifying the Pleistocene–Holocene geomorphological history in this catchment and to further understand the interplay between volcanic and glacial processes in the central Chilean Andes during the Quaternary.

Extending the mid-18th century eruptive record of the Cordón Caulle Volcanic Complex (40,5° S) through the study of three lake sedimentary sequences

The research focuses on the characterization of a tephra level extracted from sedimentary sequences obtained from three lacustrine environments from the northern Patagonia Andean range to define the volcano of origin. The dating of the tephra level corresponds to the mid-18th century but there is no scientific record nor geological evidence of this event. The sediment sequences were collected from lakes Nahuel Huapi, Moreno Oeste, and Tonček, in northern Patagonia, located downwind from the Cordón Caulle Volcanic Complex (CCVC), the potential source, at 45, 86, and 95 km respectively. These environments are frequently impacted by tephra falls from different sources in the Southern Volcanic Zone. Particle morphology, glass and mineral chemistry, together with textural analysis, provided a detailed characterization of the tephra layer recovered at each site. This characterization allowed the association of the tephra material with the CCVC predominantly, with minor Osorno volcano contributions. Historical records of both volcanic systems included eruptions in the mid-18th century (1759 CE and 1752–1796 CE, respectively). Hence, the characterization and the comparison with historical records of volcanic eruptions let us conclude that the tephra is composed of materials dispersed by CCVC in the nearby lakes, whereas contributions from Osorno volcano are observed in environments to the south of the study site, during the mid-18th century. The study demonstrates the potential of this type of detailed characterization in lake sequences to develop and support increasing research on historical volcanic records, helping to improve the knowledge of the volcanic hazard of a region.

Two-stage evolution of a bimodal reservoir: The case of Holocene lavas of the Lanín composite volcano, Southern Volcanic Zone, Chile

Lanín is a composite volcano (∼39.38°S; 71.30°W) located in the Central Southern Volcanic Zone of the Andes. It is the southernmost member of the NW-SE volcanic alignment with the Villarrica and Quetrupillán composite volcanoes. During the Holocene, Lanín volcano erupted many lava flows, which we classify into three compositional groups: BAT (basalt, andesite, and basaltic trachyandesite), intermediate, and trachyte samples. The main recognised phenocrysts in the Lanín lavas are plagioclase, clinopyroxene, orthopyroxene, Fe–Ti oxides, and apatite. These mineral phases allow us to constrain the pre-eruptive (from the magma reservoir) conditions: temperature (985–1089 °C), pressure (1.8–4.1 kbar; equivalent to 6.5–15 km depth), and oxygen fugacity patch (∼QFM+0.4). In addition, textural features combined with pre-eruptive conditions suggest the existence of a crystal mush within the reservoir, which developed in two stages. The first stage considers the emplacement, crystal mush development, and eruption related to a melt of BAT composition. The residual melt evolved from BAT to trachytic primarily through fractional crystallisation. The second stage involves an evolved trachytic melt in the reservoir and the intrusion of a new, more mafic magma (presumably of BAT composition), which generates heating of the trachytic residual melt and crystal mush disaggregation. The fractional crystallisation also gives rise to scarce products of intermediate composition. This observation is consistent with the chemical constituents plots (TAS and Harker) and thermodynamic modelling.

Origin of the compositionally zoned Paso Puyehue Tephra, Antillanca Volcanic Complex, Chile

The origin of gaps or zoning in the composition of erupted products is critical to understanding how sub-volcanic reservoirs operate. We characterize the compositionally zoned magma that produced the 2053 ± 50 cal. yr BP Paso Puyehue Tephra from the Antillanca Volcanic Complex in the Andean Southern Volcanic Zone (SVZ). The 3.7 km3 Paso Puyehue Tephra is zoned from dacite (69 wt% SiO2) lapilli and ash comprising the lowermost 80% of the deposit that abruptly transitions upward into basaltic andesite scoria (54 wt% SiO2) making up the remaining ∼20%. Variations in whole-rock, matrix glass, and mineral compositions through the deposit allow us to estimate pre-eruptive magma storage conditions and to develop a model of how this magma body was generated.Our findings suggest that amphibole-bearing basaltic andesitic magma stored at ∼8.0 ± 1.3 km depth fractionally crystallized and cooled from 1048 ± 1.1 to 811 ± 28.6 °C under highly oxidizing conditions to produce silicic a melt that upon extraction and rise, pooled at ∼6.4 ± 1.2 km depth at temperatures as low as 810 °C before eruption. MELTS models suggest that crystallization of a basaltic andesite parent magma with 4 wt% dissolved H2O can produce the dacite under conditions predicted by mineral thermobarometers with phase compositions comparable to those measured in minerals. Pervasive normal zoning at the rims of plagioclase crystals—most pronounced at the transition between dacite and basaltic andesite, and compatible vs. incompatible trace element concentrations, suggest that magma mixing was limited and likely occurred at the interface between the dacitic and basaltic andesitic magmas during ascent within the conduit upon eruption. Compositionally bimodal tephras are increasingly recognized throughout the SVZ with several interpreted to reflect basaltic recharge and mixing into extant rhyolitic reservoirs. In further contrast to other SVZ rhyolitic products, e.g., from the nearby Cordón Callue and Mocho Choshuenco volcanoes, the Paso Puyehue magma was highly oxidized. This may reflect enhanced delivery of H2O from the subducting plate into the mantle wedge, which in turn may facilitate efficient extraction and separation of buoyant, low-viscosity rhyolitic melt from crystal-rich basaltic andesitic parent magmas and the co-eruption of both end members.

Decoding the state of stress and fluid pathways along the Andean Southern Volcanic Zone

Decoding means decrypting a hidden message. Here, the encrypted messages are the state of stress, fluid pathways, and volcano tectonic processes occurring in volcanoes of the Andean Southern Volcanic Zone (SVZ). To decode these messages, we use earthquake focal mechanisms, fault slip data, and a Monte Carlo simulation that predicts potential pathways for magmatic and hydrothermal fluids. From this analysis, we propose that SVZ volcanoes have three end-member stress patterns: (i) Stress-A, a strike-slip regime coupled with the regional far-field tectonic stress; (ii) Stress-B, an extensional regime that may be promoted by volcanic edifice loading and upward pressure due to magma inflation occurring within the upper brittle-crust; and (iii) Stress-C, a local and transient fluid-driven stress rotated ~90 degrees from Stress-A. Notoriously, Stress-C pattern was observed in most volcanoes with historical eruptions. We propose that volcanoes presenting Stress-B are attractive geothermal targets, while Stress-C could be used as a predicting signal for impending eruptions.

Plio-Pleistocene rear-arc volcanism in the Southern Volcanic Zone: Eruptive styles of the Varvarco Volcanic Field

The Varvarco Volcanic Field (VVF) is located in the southern part of the Las Loicas Trough, as part of the Late Pliocene-Early Pleistocene rear-arc volcanic belt in the Transitional Southern Volcanic Zone (34.5–37°S). Its volcanic products show an elliptical distribution, elongated parallel to NW-SE main structures that regionally controlled the Las Loicas Trough. A detailed field and petrographic study was carried out to identify main lithofacies and establish its eruptive styles.The VVF magmatic evolution is initially characterized by a voluminous explosive stage represented in the area by dense and dilute pyroclastic density currents (PDC) deposits (massive lapilli tuffs, cross-stratified lapilli tuffs and diffuse-stratified tuffs). Afterwards, it evolved into an effusive stage represented by basaltic lava flows (coherent basalts), associated with Hawaiian to Strombolian eruptive style, which constitutes most of the VVF volume. The final stage of the VVF history was linked to a stratovolcano-type activity where both effusive (coherent basalts and andesites, and rhyolitic coulees) and explosive lithofacies (such as massive lithic breccias and massive lapilli tuffs) are described. Within this stage, the uppermost effusive levels were intruded by dacitic and rhyolitic domes and basaltic dykes. Available ages allow to conclude that the VVF emplacement was developed during Plio-Pleistocene times, linked to the re-steepening of the Nazca plate, after the Late Miocene Payenia shallow subduction regime.

Chironomidae responses to volcanic ash fall: a case study from Northern Patagonia for consideration in paleolimnological reconstructions

ABSTRACT Lacustrine environments from Northern Patagonia are exposed to volcanic events from the Southern Volcanic Zone. Chironomid larvae are an important paleo-proxy because, in addition to their environmental sensitivity, their chitinous cephalic capsules are well preserved in lacustrine sediments, allowing their recovery and identification. We analysed the effects of a recent volcanic eruption (Puyehue-Cordón Caulle Volcanic Complex 2011) on the chironomid community based on their cephalic remains accumulated following the ash deposition process. Two pulses of ash deposition were distinguished: a first stage comprising the immediate coarse ash fall and a delayed finer volcanic deposition. An abrupt decline in chironomid taxa and their densities were noted during the first coarse pulse, followed by an attempt of community recovery through the second finer pulse, characterised by main taxa turnover, observing an alternative community functional structure. Based on these results and comparing with past volcanic depositions, we observed different chironomid responses according to their particular ecological habits, and we noted possible percolation effects, mainly in coarse tephras, and probably a delayed record of very low abundance taxa. Our results suggest caution in analysing chironomid assemblages within and adjacent to tephra deposition as local volcanic eruption could affect chironomid community with possible consequences on paleoreconstructions.

Magma storage and transfer in the Villarrica volcanic chain, South Chile: MT insights into volcano-tectonic interactions

The active volcanic arc of the Andes is controlled by subduction of the Nazca and Antarctic Plates beneath the South American Plate. The Southern Volcanic Zone is the segment with the most active volcanoes in the Andes and the one where volcano-tectonic interactions are more evident. Magmatic fluids are there channeled through the brittle crust by faults and fractures, becoming an important region for understanding the mechanisms associated with melt migration and storage in subduction zones. The Liquiñe–Ofqui Fault System and Andean Transverse Faults, control the overall architecture of the volcanic arc and play an important role on magmatic transport and compositional partitioning. Evidence of this is the oblique chain composed of Villarrica, Quetrupillán and Lanín stratovolcanoes and a number of minor eruptive centers, where a wide range of magma compositions have been erupted with basalts and basaltic andesites dominating the suite. We used long–period and broad–band magnetotelluric stations deployed surrounding of Villarrica volcano to investigate the implications of crustal faulting and volcanism and their consequence on crustal reservoirs. Inversion of the data was used to generate three–dimensional electrical resistivity models. The resistivity distribution shows the upper crust as highly resistive, but below and east of Villarrica volcano, the model suggests the presence of a magmatic reservoir at shallow crustal levels (between 1.5 and 3 km b.s.l.), possibly corresponding to a transient magma storage. Meanwhile, the middle crust contains several intermediate to low conductive features interpreted as fluid pathways and/or melt storage regions, respectively, revealing the important role of fault systems. The lower crust also contains low resistivity zones indicating the presence of partial melt and/or fluids, associated with deep reservoirs (8–20 km), with a significant proportion of them likely non eruptible parts. This would suggest that melt is accumulated as highly crystallized mush or disconnected melt pockets and highlight the complex vertical extent of the structurally-controlled plumbing systems even in thin crust settings.

Hydrothermal Alteration in the Nevados de Chillán Geothermal System, Southern Andes: Multidisciplinary Analysis of a Fractured Reservoir

The interplay between a heat source, primary plus secondary permeability, and hydrothermal fluids makes geothermal systems a highly dynamic environment where evolving physico-chemical conditions are recorded in alteration mineralogy. A comprehensive characterization of hydrothermal alteration is therefore essential to decipher the major processes associated with geothermal system development. In this study, we defined the hydrothermal mineralogical evolution of the Nevados de Chillán Geothermal System (NChGS), located in the Southern Volcanic Zone (SVZ) of the central Andes, where the regional framework of the system is formed by a direct association with a currently active volcanic complex, a favorable structural control, and vertically inhibited fluid circulation. To characterize the secondary mineralogy present in the NChGS, we integrated optical petrography, Scanning Electron Microscopy (SEM) observations, X-ray Diffraction (XRD) analysis, and microthermometric measurements along a drill core with a depth of 1000 m at the Nieblas-1 well. These mineralogical approaches were combined with a structural field analysis to highlight the relevance of multidisciplinary study in understanding active geothermal systems. The results indicated that the evolution of the system involved four paragenetic stages, with the main processes in each phase being the heating, boiling, and mixing of fluids and re-equilibration to new physico-chemical conditions. Additionally, three hydrothermal zones were recognized: an upper argillic section, an intermediate sub-propylitic zone, and a deep propylitic domain. Sampled thermal springs are characterized by pH values of 2.4–5.9 and high SO4= concentrations (>290 ppm). These acid-sulfate steam-heated waters suggest the contribution of primary magmatic volatiles to the hydrothermal system. Alunite recorded in the alteration halos of veinlets presents at depths of 170–230 m denote the circulation of acidic fluids at these levels which were favored by reverse faults. These findings indicate that, at this depth range, the condensation of magmatic volatiles into shallow aquifers controls the recharge area of the superficial thermal manifestations. Conversely, deep-seated hydrothermal fluids correspond to near-neutral chloride fluids, with salinities ranging from 0.1 to 6.9 wt.% NaCl eq. The distribution of illite/smectite and chlorite/smectite mixed-layered minerals outline the presence of a significant clay cap, which, in this system, separates the steam-heated domain from the deep hydrothermal realm and restricts fluid circulation to existing permeable channels. Our mineralogical and structural study provides critical data for the interpretation of heat–fluid–rock interaction processes in the NChGS. The interplay between hydrothermal fluids and active faults is also discussed in the context of the complex of geological processes in active geothermal systems along the Chilean Southern Volcanic Zone.

Identifying analogues for data-limited volcanoes using hierarchical clustering and expert knowledge: a case study of Melimoyu (Chile)

Determining the eruption frequency-Magnitude (f-M) relationship for data-limited volcanoes is challenging since it requires a comprehensive eruption record of the past eruptive activity. This is the case for Melimoyu, a long-dormant and data-limited volcano in the Southern Volcanic Zone (SVZ) in Chile with only two confirmed Holocene eruptions (VEI 5). To supplement the eruption records, we identified analogue volcanoes for Melimoyu (i.e., volcanoes that behave similarly and are identified through shared characteristics) using a quantitative and objective approach. Firstly, we compiled a global database containing 181 variables describing the eruptive history, tectonic setting, rock composition, and morphology of 1,428 volcanoes. This database was filtered primarily based on data availability into an input dataset comprising 37 numerical variables for 438 subduction zone volcanoes. Then, we applied Agglomerative Nesting, a bottom-up hierarchical clustering algorithm on three datasets derived from the input dataset: 1) raw data, 2) output from a Principal Component Analysis, and 3) weighted data tuned to minimise the dispersion in the absolute probability per VEI. Lastly, we identified the best set of analogues by analysing the dispersion in the absolute probability per VEI and applying a set of criteria deemed important by the local geological service, SERNAGEOMIN, and VB. Our analysis shows that the raw data generate a low dispersion and the highest number of analogues (n = 20). More than half of these analogues are in the SVZ, suggesting that the tectonic setting plays a key role in the clustering analysis. The eruption f-M relationship modelled from the analogue’s eruption data shows that if Melimoyu has an eruption, there is a 49% probability (50th percentile) of it being VEI≥4. Meanwhile, the annual absolute probability of a VEI≤1, VEI 2, VEI 3, VEI 4, and VEI≥5 eruption at Melimoyu is 4.82 × 10−4, 1.2 × 10−3, 1.45 × 10−4, 9.77 × 10−4, and 8.3 × 10−4 (50th percentile), respectively. Our work shows the importance of using numerical variables to capture the variability across volcanoes and combining quantitative approaches with expert knowledge to assess the suitability of potential analogues. Additionally, this approach allows identifying groups of analogues and can be easily applied to other cases using numerical variables from the global database. Future work will use the analogues to populate an event tree and define eruption source parameters for modelling volcanic hazards at Melimoyu.

A 2200-year record of Andean Condor diet and nest site usage reflects natural and anthropogenic stressors.

Understanding how animals respond to large-scale environmental changes is difficult to achieve because monitoring data are rarely available for more than the past few decades, if at all. Here, we demonstrate how a variety of palaeoecological proxies (e.g. isotopes, geochemistry and DNA) from an Andean Condor (Vultur gryphus) guano deposit from Argentina can be used to explore breeding site fidelity and the impacts of environmental changes on avian behaviour. We found that condors used the nesting site since at least approximately 2200 years ago, with an approximately 1000-year nesting frequency slowdown from ca 1650 to 650 years before the present (yr BP). We provide evidence that the nesting slowdown coincided with a period of increased volcanic activity in the nearby Southern Volcanic Zone, which resulted in decreased availability of carrion and deterred scavenging birds. After returning to the nest site ca 650 yr BP, condor diet shifted from the carrion of native species and beached marine animals to the carrion of livestock (e.g. sheep and cattle) and exotic herbivores (e.g. red deer and European hare) introduced by European settlers. Currently, Andean Condors have elevated lead concentrations in their guano compared to the past, which is associated with human persecution linked to the shift in diet.

Editorial: Volcanism in the Central Volcanic Zone of the Andes

This Research Topic was born out of the 1 st Symposium on the Central Volcanic Zone of the Andes 21 (CVZA; carried out in virtual mode during 2021) which counted for the presence of geoscientists 22 working in the volcanological, petrological and geochemical fields, among the others. Our objective 23 was to bring together contributions focused on the evolution and impact of volcanism in the CVZA 24 from its formation (Late Oligocene) up to present times, including specific and multidisciplinary 25 studies, with both detailed and regional approaches. The nine articles published in this Research Topic 26 cover a spectrum of topics applied in different areas of the CVZA, including poly-and monogenetic 27 volcanism, the influence of the crust in magmatic processes, fluid geochemistry, volcano seismicity, 28 petrology, geochronology, hazards assessment, relationships between volcanoes and 29 archeology/anthropology/sociology. 30 Barry et al. (2022) compare the degassing processes between the Central and Southern volcanic zones 31 of the Andes using He and carbon isotopes, and CO2/ 3 He ratios obtained from gas and water samples 32 related to fumaroles and thermal springs. They conclude that crustal thickness exerts a primary control 33 on the extent of fluid-crust interaction, as helium and other volatiles rise through the upper plate in the 34 Andean convergent margin, with more significant crustal inputs in the CVZA due its thicker crust. 35

The Magmatic Evolution and the Regional Context of the 1835 AD Osorno Volcano Products (41°06’S, Southern Chile)

Abstract Osorno volcano (41°06’S, 72°20’W) is a composite stratovolcano of the Central Southern Volcanic Zone of the Chilean Andes. It is the southernmost member of a NE–SW trending alignment of volcanic edifices including La Picada and Puntiagudo volcanoes and the Cordón Cenizos chain. According to contemporary descriptions recorded by Charles Darwin in 1835, two eruptive events occurred: the first during January–February, and the second during November–December 1835 and January 1836. The volcano erupted basaltic andesite lavas and tephra fall deposits (52.4 to 52.9 SiO2 wt. %), which contain phenocrysts of olivine, plagioclase, clinopyroxene, and spinel. The compositions of these phenocryst phases, together with those of olivine-hosted melt inclusions, allowed us to constrain intensive parameters for the pre-eruptive magmas. These varied from 1060°C to 1140°C, with an oxygen fugacity buffer of ~ΔQFM +1.1, dissolved water concentrations of up to 5.6 wt. % (average of ~4.2 wt. %) and maximum pressures equivalent to ~7-km depth. Textural relations, such as crystal accumulations and clots, zoning in crystals and other indications of disequilibrium, lead us to infer the involvement of a crystal mush, rich in individual crystals and clots of crystals, which underwent a degree of disaggregation and entrainment into the transiting magma prior to eruption. Comparison of trace element abundances, including rare earth elements, fluid-mobile elements, and relatively fluid-immobile elements, combined with 87Sr/86Sr and 143Nd/144Nd isotope ratios, allows us to consider variations in slab-derived fluid input and the minor role of crustal contamination on the Osorno eruptive products and those from neighboring volcanic systems. Our results suggest both a greater contribution from slab-derived fluid and a higher degree of partial melting in the systems supplying stratovolcanoes (Osorno, Calbuco, and La Picada) relative to those supplying small eruptive centers built over the major regional Liquiñe-Ofqui Fault Zone.

Petrology of the magmatic system beneath Osorno volcano (Central Southern Volcanic Zone, Chile)

Petrological studies of volcanic systems and magmatic processes provide small-scale information that are inaccessible to the current resolution of geophysical tools. In the Southern Volcanic Zone of the Chilean arc, Osorno ranks sixth in the official risk classification in a region of increasing vulnerability due to soaring human activities coupled with lahar hazard and the threat of a tsunami induced by a flank collapse. However, its magmatic system remained poorly investigated. Here, we present and integrate to existing geophysical models a comprehensive set of petrographic and geochemical data including whole-rock and mineral major and trace element analyses, associated with detailed numerical modelling to constrain the storage conditions below Osorno. To capture the full complexity of the system, over 154 samples of all known major units of the volcano have been collected. They cover a large range from tholeiitic primary basalt (Mg# = 0.72 and ≥ 50 wt% SiO 2 ) to rhyodacite (Mg# = 0.18 and ≥ 70 wt% SiO2), displaying a compositional gap within andesite (59–63 wt% SiO 2 ). This gap results either from andesitic melt thermal instability or from the interstitial melt extraction of the crystal mush. Basalt and basaltic andesite lava with plagioclase and olivine (± clinopyroxene) are dominant. Their crystallinity largely varies from 1 to 53 area % with the lowest values in the most evolved basaltic andesites. Dacites are limited to three small domes with low crystallinities between 7 and 13 area %. The presence of diktytaxitic enclaves within the dacites indicate minor mingling with a less differentiated melt. Water-bearing phases are generally absent, except for one dacite sample where few small amphibole crystals occur. Petrology, chemical data and thermobarometric results imply shallow fractional crystallization of troctolitic, gabbroic and gabbronoritic cumulates. Differentiation dominantly takes place between 2 and 3 kbar (6–10 km) and results from a water-poor (≈1 wt% H 2 O), tholeiitic parental melt. No evidence of high-pressure fractionation was observed. We interpret this differentiation depth as Osorno's main storage zone. It correlates with the depth of the intracrustal discontinuity and seismic reports below the volcano. Only the upper half of the storage zone, imaged with geophysical methods, was erupting. We suggest a comparable behavior for a potential future event. • Basalt evolves to dacite at ≈2 kbar near the depth of the intracrustal discontinuity. • Fractional crystallization dominates under water undersaturated conditions. • Petrographic observations and numerical models were integrated with geophysical data.