Caldera Research Articles

Gravity Data Allow to Image the Shallow‐Medium Subsurface Below Mud Volcanoes

AbstractThe debate about the conceptual model of mud volcanoes functioning is still alive in the literature. A large part of the literature focuses on the characterization of the deep reservoir where expelled fluids are expected to originate. Another part of literature is focused on the study of the shallow system of mud volcanoes, which could influence the short‐term variations in mud volcanoes activity. We present and analyze a new data set of micro‐gravimetric data to study the area of the Nirano Salse, Italy. Unlike what is commonly assumed for the study area, our results suggest that the geomorphology of the Nirano Salse is not related to a caldera collapse above a shallow mud chamber, but to the surface expression of slip distribution of a fault termination along which the fluids ascended to the surface. We believe that gravimetric data can significantly improve the study of hydrocarbon seeps and mud volcanism.

A shifty Toba magma reservoir: Improved eruption chronology and petrochronological evidence for lateral growth of a giant magma body

Polycyclic caldera complexes hold clues to understanding why some magmatic systems develop into supersized magma bodies and how they can recover to produce several caldera-forming eruptions. However, the geologic records of the transitions between successive caldera events are very often inaccessible due to limited preservation of eruptive products of inter-caldera activity, prompting the search for alternative archives of magma evolution such as accessory minerals. Here we applied multiple geochemical tools to study one of the most active caldera centres of the Quaternary, the Toba caldera complex in Sumatra (Indonesia), which produced at least four caldera-forming eruptions in the last 1.6 My, including the iconic Youngest Toba Tuff at 74 ka. We combined feldspar 40Ar/39Ar and zircon U–Pb geochronology of proximal pyroclastic deposits with glass and mineral chemistry of both the tuffs and distal marine tephra to revise the eruption chronology of Toba, obtaining new eruption ages of 1417 -31/+14 ka (zircon) or 1339 ± 39/39 ka (plagioclase, internal/full external 2σ uncertainty) for the Haranggaol Dacite Tuff, 783.81 ± 0.85/1.32 ka (sanidine) for Oldest Toba Tuff, and 503.61 ± 1.36/1.50 ka (sanidine) for Middle Toba Tuff. Isotope dilution thermal ionisation mass spectrometry (ID-TIMS) U–Pb crystallisation ages, trace element contents and Hf isotopic ratios of zircons illuminate changes in the shallow magma reservoir which saw near-continuous zircon crystallisation over 1.6 My. Prolonged build-ups to each eruption with highly scattered zircon trace element compositions reflect a complex, heterogeneous character of the shallow reservoir, without a clear temporal trend or indications of the eruption trigger. In contrast, hafnium isotopes in zircon display a pronounced shift towards unradiogenic values immediately after the OTT caldera collapse, followed by a gradual recovery to a baseline value of εHf = -7 at the time of YTT eruption, interpreted as a reflection of the shift in magma reservoir position corresponding to change in the character of assimilated crust. We can show in unprecedented detail how a large caldera collapse affects magma geochemistry; however, identification of patterns in the behaviour of the Toba system and making geochemistry-based predictions about its future development remain a challenge.

The Colíder Paleoproterozoic felsic volcanism: New insights into stratigraphy and petrogenesis in the southern Amazonian Craton

The eastern side of Alta Floresta Mineral Province (AFMP) (União do Norte region) comprises an effusive/subvolcanic phase, explosive phase, and epiclastic rocks. The effusive/subvolcanic phase comprises massive porphyritic rhyolite with different degrees of crystallinity, porphyritic rhyolite with flow structure, and microcrystalline rhyolite. The explosive phase is characterized by welded lapilli tuff, lapilli tuff with flow structure, stratified lapilli tuff, and tuff with accretionary lapilli. Sandstones and siltstones characterize epiclastic rocks. There are porphyritic rhyolite dikes that cut the volcanic sequence. The rocks of the Colíder Group are rhyolites (>70 % SiO2) with calcic-alkaline to alkaline-calcic affinity. A rhyolitic rock is predominantly ferroan and metaluminous, while a pyroclastic rock is predominantly magnesian to ferroan and meta- to peraluminous in composition. The rocks have a high LREE content, and slight fractionation of HREE suggests they have hybrid geochemical characteristics of high potassium calc-alkaline, shoshonitic, and suggestive affinities of A-type magmas, probably produced in a post-collisional orogenic environment. These rocks have a relatively high level of Ba and Rb and a negative anomaly for Sr, P, Nb, Ta, and Ti, suggesting an intra-plate mantle source with contributions from a metasomatized mantle edge. The rocks do not have adakitic characteristics nor fertility for Cu-Au (Sr/Y < 1.9 and LaN/YbN < 38.0) as in other regions of the AFMP. The welded lapilli tuff rocks show U-Pb age of 1838 ± 17 Ma and 1817 ± 2 Ma and are the older explosive phases. The effusive/subvolcanic phase has U-Pb ages of 1800 ± 3 Ma and 1792 ± 3 Ma obtained in the massive porphyritic rhyolite. The age interval of approximately 46 Ma (1838–1792 Ma) suggests the presence of two or more volcanic cycles or distinct volcanic events. Based on the interpretation of volcanic and epiclastic deposits and remote sensing products, the paleoenvironmental reconstruction includes subaerial with effusive/subvolcanic and explosive stages in a probable eroded volcanic caldera system with epiclastic rocks from a fluvial depositional system. This volcanic scenario provides more information about the Colíder–Teles Pires volcano-plutonic event, uniquely mapped in AFMP, and offers new insights into the tectonic frame of the Amazonian craton.

Numerical models of restless caldera volcanoes

Caldera volcanoes undergo cycles of inflation and deflation on both long and short timescales. This study uses the Distinct Element Method to model numerically a pressurized magma reservoir within a brittle host rock in order to investigate the effects on these cycles of: (i) host rock strength, (ii) strain localization, (iii) fault reactivation and (iv) surface deformation. Additionally, the critical over- or under-pressure within the magma reservoir - i.e. the pressure prior to reservoir roof failure leading to eruption or collapse - is determined. The model results indicate that with increasing number of inflation/deflation cycles, the host rock weakens and strain is increasingly localized onto faults that are reactivated during both inflation and deflation. Consequently, both the critical over- and under-pressure magnitudes for eruption or collapse decrease from cycle to cycle. The relationships of surface uplift and subsidence with magma reservoir over- and under-pressure, respectively, display hysteresis behavior. This means that the magnitude of elevation change as a function of the change of magma pressure differs for inflating and deflating magma reservoirs. These results provide new insight into the relationships between structural development, elevation change, pressure change, and the generation of potential eruptive pathways at caldera volcanoes.

Insights into the explosive eruption history of the Campanian volcanoes prior to the Campanian Ignimbrite eruption

The Campanian Volcanic Zone (CVZ) comprises multiple active volcanoes and includes the highly productive Campi Flegrei and Ischia caldera systems. These caldera volcanoes have produced probably the largest eruptions in Europe in the past 200 ka, such as the Monte Epomeo Green Tuff (MEGT; Ischia) at ca. 56 ka and the Campanian Ignimbrite (CI; Campi Flegrei) at ca. 40 ka, which form widespread isochrons across the Mediterranean region. These closely-spaced volcanic centres erupt phonolitic to trachytic glass compositions that are similar, and thus it can be challenging to correlate tephra deposits to specific volcanic sources. Here we present a detailed tephrostratigraphy for pre-CI eruption activity using the units preserved within a sequence at the coastal Acquamorta outcrop, on the western side of the CI caldera rim. Both the MEGT and CI units are present in the section, and they bracket twelve eruption units that were logged and sampled. New major and trace element glass chemistry data have been acquired for these Acquamorta tephra deposits. Three eruption deposits from Ischia and nine from Campi Flegrei are identified, which helps constrain the tempo of volcanic activity of these centres between the large caldera-forming eruptions. The three Ischia tephra deposits between the MEGT and the CI are indistinguishable based on both major and trace element glass chemistry and cannot be correlated to a specific or known eruption in this interval, such as the Schiappone tephra. The compositional variations between the Campi Flegrei eruptions reveal temporal shifts in the composition of the tephra deposits that reflect changes in the magmatic system prior to the CI eruption. These deposits indicate that there were at least nine eruptions at Campi Flegrei within 16 ka of the enormous CI eruption, and suggest that there was no significant period of repose before the caldera generating eruption.

Atmosphere injection of sea salts during large explosive submarine volcanic eruptions

The 15 January 2022 submarine eruption at Hunga volcano was the most explosive volcanic eruption in 140 years. It involved exceptional magma and seawater interaction throughout the entire submarine caldera collapse. The submarine volcanic jet breached the sea surface and formed a subaerial eruptive plume that transported volcanic ash, gas, sea salts and seawater up to ~ 57 km, reaching into the mesosphere. We document high concentrations of sea salts in tephra (volcanic ash) collected shortly after deposition. We also discuss the potential climatic consequences of large-scale injection of salts into the upper atmosphere during submarine eruptions. Sodium chloride in these volcanic plumes can reach extreme concentrations, and dehalogenation of chlorides and bromides poses the risk of long-term atmospheric and weather impact. Salt content in rapidly collected tephra samples may also be used as a proxy to estimate the water:magma ratio during eruption, with implications for quantification of fragmentation efficiency in submarine breaching events. The balance between salt loading into the atmosphere versus deposition in ash aggregates is a key factor in understanding the atmospheric and climatic consequences of submarine eruptions.

NUMERICAL SIMULATION OF TSUNAMI IMPACT FROM THE 1/15/22 ERUPTION OF THE HUNGA TONGA - HUNGA HA’APAI VOLCANO

Tsunamis from volcanic ‘explosive’ eruptions are rare, with the last catastrophic event being Krakatau in 1883 (Verbeek, 1885), during which, tsunamis were generated in the far-field by pressure shock-waves and in the nearfield of the volcano, in the Sunda Straits, by several potential geological mechanisms including pyroclastic flows, ash column, and/or caldera collapse. On 1/22/55, at about 4:15 UTC, a one in 1,000 year eruption of the Hunga Tonga-Hunga Ha’a-pai Volcano (HTHHV), that had started on12/20/21, reached its paroxysm with a series of large underwater explosions, releasing enormous energy (4-18 Mt of TNT), and ejecting a large ash plume 58 km into the stratosphere. We simulate both the near- and far-field tsunami generation from the eruption, but in this paper we focus on analyzing and validating the near-field impact against field data.

Complex styles of phreatomagmatic explosions at Kīlauea Volcano, Hawaii, controlled by magma structure

Explosive eruptions at basaltic volcanoes remain poorly understood. Kīlauea Volcano is a type locality for basaltic eruptions and is well-known for effusive activity. However, more than 7 m of phreatomagmatic Keanakākoʻi Tephra unit D deposits from explosive eruptions crown the southern rim of the summit caldera and provide a stark reminder of Kīlauea’s explosive past and future potential. We used detailed field observations as well as granulometric and morphological analysis of 100 samples from two proximal sections to assess the eruption style and fragmentation mechanism. The deposits can be divided into four subunits, six different lithofacies, and contain three juvenile tephra components. Each juvenile component shows distinct shape variability resulting from molten fuel-coolant interaction (MFCI) explosions of magma of variable vesicularity. Fragmentation of dense glass generates olive-green ash, fragmentation of low to moderately vesicular magma generates a dark gray ash-lapilli component, and fragmentation of highly vesicular magma generates light-yellow pumice. Our work shows that magma structure impacts MFCI explosion efficiency. Small-scale planar bedding throughout most of the deposit points to a general eruption style of small, frequent explosions generating low plumes. Thicker beds of accretionary lapilli of fine-extremely fine ash are related to very efficient magma-water mixing. Pyroclastic density current (PDC) deposits in the upper part of the stratigraphy contain at least three flows but show no significant dune or cross-bedding structures. We suggest that this is a function of the vent being situated in a caldera that was then ∼600 m deep, where the caldera wall acted as a barrier and changed the flow dynamics to very dilute overspills and co-PDC plume falls over the wall. Deconvolution modeling of the polymodal grain size distributions is used to assess grain size changes of each juvenile component for this deposit, which greatly improves interpretation of lithofacies generation and eruption dynamics. Size-correlated shape parameters show that shape data across a wide size range are needed to accurately track grain shapes. This study demonstrates how careful examination of grain size and shape of juvenile tephra clasts can help volcanologists understand how effusive basaltic volcanoes can become violently explosive.

Collapse caldera walls: Mechanical controls on slope failure and morphology

Volcanic calderas are enclosed topographic depressions that are commonly delimited by steep annular cliffs (caldera walls). The walls represent degraded scarps of ring faults that accommodated caldera collapse. Here, we assess the mechanical controls on caldera wall morphology and failure by: (i) analysing slope properties of recently formed calderas; (ii) examining analytical solutions for slope failure; and (iii) simulating progressive slope failure via Distinct Element Method (DEM) numerical modelling. Average slope angles of the caldera walls in nature range from 20° to 65° and slope heights range from 99 m to 1080 m. Smaller slope heights are weakly tied to greater slope inclinations. Almost 60% of sampled caldera wall profiles show a steeper upper section and a gentler lower section, which in numerical models are reproduced as an upper slope comprising intact wall rock and a lower slope comprising failed material inclined at the angle of repose. Moreover, modelled material with a low angle of internal friction and low cohesion exhibits shear failure, whereas material with high friction and high cohesion shows tension cracking, buckling and toppling. The slope height and inclination data from the numerical models agree well with those in nature. The numerical model data correspond more closely to analytically-predicted failure planes than critical (i.e. limit equilibrium) slopes because they are short-lived and therefore unlikely to be sampled in either model or nature. If rock mass cohesions or friction angles constrained analytically from natural slopes are determined by assuming the slope is critical, they are likely to be underestimates. The homogeneous numerical models here give an estimated natural rock mass cohesion of 0.5 MPa or less. This is much lower than typical laboratory-scale cohesion values, and it highlights the importance of strength reduction when upscaling mechanical parameters from laboratory experiments for modelling slope stability and volcano deformation.

Mineralization age of the Xiangshan uranium ore field, South China redefined by hydrothermal apatite U-Pb geochronology

The Xiangshan volcanic-related uranium ore field, located in the Ganzhou-Hangzhou Uranium-Polymetallic Metallogenic Belt (GHMB) in South China, has been the country’s largest historic uranium producer. Although this ore field has been intensely studied for decades, the absolute age of uranium ore formation was still debated. Conventional U-Pb dates on U-minerals (i.e., pitchblende, brannerite) and Ar-Ar dates on associated alteration minerals (i.e., illite) from different uranium deposits yield a wide range of ages (122 to 85 Ma). The uncertainty regarding the exact timing of uranium mineralization at Xiangshan has limited the development of accurate ore genesis models and precluded a better understanding of the links between metallogenesis and regional-scale geodynamics during the Cretaceous-Tertiary period in South China.In this work, we focus on applying in-situ LA-ICP-MS U-Pb geochronology to magmatic and hydrothermal apatite from two typical uranium ore deposits at Xiangshan, intending to indirectly constrain the exact timing of uranium mineralization and reconstruct the magmatic-hydrothermal history. U-Pb age determinations on magmatic zircon and apatite from the three types of ore-hosting volcanic rocks yielded precise crystallization ages of 133 to 136 Ma and 135 to 140 Ma, respectively. These results are identical within uncertainty to previously reported U-Pb dates obtained in the study area, and thereby tightly constrain the emplacement age of the Xiangshan volcanic complex that hosts uranium deposits. Moreover, the potassium feldspar-biotite Ar-Ar dating for the major volcanic rocks demonstrates that these host rocks were cooled and completely solidified through the biotite/potassium feldspar closure temperature (below 300 °C) at 133–131 Ma. Subhedral or anhedral aggregates of hydrothermal apatite ubiquitously formed in the high-grade uranium ores at the Shannan and Zoujiashan deposits and are commonly intergrown with fluorite, illite and uranium minerals. Relative to magmatic apatite from the host rocks, the hydrothermal apatite is typically characterized by much higher U and Th contents, depletion of LREEs, enrichment in HREEs and a much more pronounced negative Eu anomaly. As such, the hydrothermal apatite is unequivocal of hydrothermal origin and formed synchronously with uranium mineralization. LA-ICP-MS U-Pb dating of these hydrothermal apatite crystals yielded an absolute 206Pb/238U age of 126.7 ± 0.6 Ma for the Shannan deposit and 131.3 ± 7.2 Ma for the Zoujiashan deposit. The new ages are much older (>10 Ma) than those previously reported ages on ore or gangue minerals such as pitchblende and illite. Thus, one unidentified mineralization event likely occurred at ca. 131–127 Ma, soon after solidification and subsequent brittle deformation or faulting within the volcanic caldera. We, therefore, proposed that volcanic magmatism was highly likely to be the ‘heat engine’ driving hydrothermal fluid convection and possibly provide the essential conditions (i.e., magmatic fluids, U sources) for the formation of the U mineralization. Another implication of this study is that dating hydrothermal apatite with high U concentrations associated with uranium mineralization could be considered a complementary geochronological tool in the study of analogous hydrothermal uranium deposits worldwide.

The 2018 Eruption of Kīlauea: Insights, Puzzles, and Opportunities for Volcano Science

The science of volcanology advances disproportionately during exceptionally large or well-observed eruptions. The 2018 eruption of Kīlauea Volcano (Hawai‘i) was its most impactful in centuries, involving an outpouring of more than one cubic kilometer of basalt, a magnitude 7 flank earthquake, and the volcano's largest summit collapse since at least the nineteenth century. Eruptive activity was documented in detail, yielding new insights into large caldera-rift eruptions; the geometry of a shallow magma storage-transport system and its interaction with rift zone tectonics; mechanisms of basaltic tephra-producing explosions; caldera collapse mechanics; and the dynamics of fissure eruptions and high-volume lava flows. Insights are broadly applicable to a range of volcanic systems and should reduce risk from future eruptions. Multidisciplinary collaboration will be required to fully leverage the diversity of monitoring data to address many of the most important outstanding questions. ▪Unprecedented observations of a caldera collapse and coupled rift zone eruption yield new opportunities for advancing volcano science.▪Magma flow to a low-elevation rift zone vent triggered quasi-periodic step-like collapse of a summit caldera, which pressurized the magma system and sustained the eruption.▪Kīlauea's magmatic-tectonic system is tightly interconnected over tens of kilometers, with complex feedback mechanisms and interrelated hazards over widely varying timescales.▪The eruption revealed magma stored in diverse locations, volumes, and compositions, not only beneath the summit but also within the volcano's most active rift zone.

Frictional melting mechanisms of rocks during earthquake fault slip

Rapid slip, at rates in the order of 1 m/s or more, may induce frictional melting in rocks during earthquakes. The short-lived melting has been thought to be a disequilibrium process, for decades. We conducted frictional melting experiments on acidic, basic, and ultrabasic silicate rocks at a slip rate of 1.3 m/s. The experiments and microstructural observations reveal that all minerals in the rocks are melted at temperatures below their known melting temperatures (Tm); e.g., quartz is melted at ~ 1000–1200 °C, not ~ 1720 °C, while olivine at ~ 1300 °C, rather than ~ 1700 °C. The low-temperature melting is incompatible with the conventional disequilibrium melting, and may be caused predominantly by grain size reduction and phase boundary reactions during the early and later stages of slip, respectively. The newly estimated Tm and the melting mechanisms should be considered for understanding the mechanics of earthquakes, landslides, and caldera collapses.

In-situ monitoring of 3He/4He in summit gases of Kilauea Volcano (Hawaii) prior to the 2020 eruption

We present He isotope (3He/4He) data from a fumarole and near-ground gases measured in-situ at the Sulfur Banks solfatara field at the summit of Kilauea Volcano, Hawaii. We used a field-deployable mass-spectrometer-based system: the Helium Isotope Monitor (HIM) previously described in McMurtry et al. (2019a, b). The in-situ instrument was deployed using solar power for the first time and results were ground-truthed against data determined using conventional gas analytical and noble gas mass spectrometry techniques. The HIM instrumentation, associated Vent Gas Purification System (VGPS), and related sampling equipment and strategy are described. Cloudy and rainy weather conditions hampered the deployment, which was reorganized to reduce power loads and resulted in less sampling than planned. Nevertheless, we obtained daily sampling of the volcanic vent gas. Results from the Old Well fumarole indicate a ~ 2 RA increase in 3He/4He on the day of the December 20th, 2020 eruption of nearby Halema‘uma‘u Crater, reaching 17.0 RA using the in-situ instrument and 16.0 ± 0.67 RA using conventional techniques. This finding suggests that a new 3He-enriched magma source is driving the current, ongoing eruption phase of Kilauea and, if so, confirms that the deep summit caldera fault system that hosts the Sulfur Banks field is connected to the Halema‘uma‘u Crater magmatic system. Overall, these findings illustrate how time-series helium isotope data, which are well established by ongoing discrete monitoring at low temporal resolution, can help forecast forthcoming eruptive events that may not be foreseen by other volcanic monitoring methods.

Migration of Mechanical Perturbations Estimated by Seismic Coda Wave Interferometry During the 2018 Pre‐Eruptive Period at Kīlauea Volcano, Hawaii

AbstractWe use seismic ambient noise correlation and coda wave interferometry to estimate velocity variations at high temporal resolution, during the pre‐eruptive period and the onset of the 2018 eruption of Kīlauea volcano. A progressive velocity increase is observed from March to the end of April. It is followed by rapid decrease starting a few days before the onset of the East Rift Zone (ERZ) eruption and then by sharp velocity drop when the eruption started. The change of trend from velocity increase to decrease is progressively delayed by a few days from the summit caldera toward the ERZ. The location of the velocity perturbations shows a migration of the sources of velocity changes from the summit caldera toward the ERZ before the eruption. Using a model of pressure source, we show that the simultaneous caldera inflation and velocity increase probably result from an anisotropic distribution of fault and crack orientations. The velocity decrease could be due to damaging processes above the shallow reservoir and to plastic deformations around the caldera. We introduce a forward model of rock damage associated with the volcano‐tectonic seismicity to calculate the velocity decrease. The good agreement between the calculated and the observed velocity variations shows that a large part of the velocity decrease results from damage of the medium. The delayed onsets of velocity decrease and the source migration of velocity perturbations are probably related to progressive fault openings in the Southern and Eastern parts of the caldera and to magma transfer toward the ERZ.

The Taverna San Felice Dike (NE of Roccamonfina Volcano): Unraveling Magmatic Intrusion Processes and Volcano‐Tectonics in the Tyrrhenian Margin of the Southern Apennines

AbstractThe Roccamonfina volcano is located within the Garigliano Graben (southern Apennines, Italy) and has been active throughout the Middle‐Late Pleistocene. Along its polyphase volcanic history (630–55 ka), including several caldera‐forming eruptions (385–230 ka), several effusive/mildly explosive monogenetic events occurred along the volcano slopes, within the summit caldera, and along the graben‐bounding carbonate reliefs. In this paper, we present a multidisciplinary study of a mafic magmatic feeder dike intruded within the Meso‐Tertiary carbonates and overlying Lower Pleistocene breccias of Mt Cesima, northeast of the Roccamonfina volcano. We performed a stratigraphic and structural survey of the area and petrographic analyses on several samples of the dike. Results indicate that a ∼1 km long fissure fed an eruption that also emplaced a Strombolian pyroclastic sequence. Petrological data show that an open‐system mafic recharge fueled the tephritic magma that fed the eruption, whereas no evidence of significant pre/syn‐eruptive assimilation of carbonate has been identified. Stratigraphic and petrological data do not allow to firmly constrain the timing of the eruption, which could belong both to the pre‐Brown Leucitic Tuff (&gt;354 ka) and to the post‐White Trachytic Tuffs (&lt;230 ka) epochs of activity of the Roccamonfina volcano. Structural data show that the dike is broadly oriented E‐W and changes direction toward NE‐SW in correspondence with a pre‐existing fault damage zone. We suggest that magma was intruded during an N‐S trending extensional event in the Middle Pleistocene, whose prolonged activity resulted in regional uplift and exhumation of regional significance.

Caldera-Like Delineation and Geological Resources Identification of Majenang Area, Indonesia by Gravity and Magnetic Data Inversions

Indonesia’s archipelago formed by the tectonic evolution proceeded with subduction, which was accompanied by volcanism. The systematic subduction zones produce the magmatic arcs with different periods since the Permian up to the Tertiary. However, only the recent Quaternary volcanic arc is recognized and lacked of information about the ancient volcanic environment. Calderas are a crucial feature in any volcanic environment due to the prospect site of geological resources. The Gravity and Magnetic methods are commonly used for preliminary study in almost any cases due to their light-weight, low-cost, and ability to map a wide area rapidly. The inverse-modeling scheme was invoked to estimate the sub-surface situation during the interpretation processes. The study intended to show the ability of both gravity and magnetic method for delineating of caldera-like environment including geological resources prospective site identification. According to the research, a northwest-southeast dextral strike-slip fault found in the area and belongs to the Pamanukan-Cilacap Fault Zone (PCFZ). A circular caldera-like anomaly delineated and interpreted as the ring fault of an ancient volcanic caldera in the study area. Several high gravity anomalies found within the caldera rims are interpreted as lava domes or intrusion rocks, while the high-density found following the outer part of the ring-fault terrain is interpreted as the buried lava. The ancient eruption point inferred around Majenang city, thus the study area proposed as the Majenang Caldera. A mineralization zone identified around the study area, which is comparable with the Cihonje people’s gold mining site as the proven prospective area.

Object-Oriented Remote Sensing Approaches for the Detection of Terrestrial Impact Craters as a Reconnaissance Survey

The purpose of this study is to employ a remote sensing reconnaissance survey based on optimal segmentation parameters and an object-oriented random forest approach to the identification of possible terrestrial impact craters from the global 30-m resolution SRTM DEM. A dataset consisting of 94 confirmed and well-preserved terrestrial impact craters, 104 volcanic calderas, and 124 valleys were extracted from real-world surface features. For craters with different sizes, eight optimal scale parameters from 80 to 3000 have been identified using multi-resolution segmentation, where the scale parameters have a positive correlation (R2 = 0.78) with the diameters of craters. The object-oriented random forest approach classified the tested impact craters, volcanic calderas, and valleys with an overall accuracy of 88.4% and a Kappa coefficient of 0.8. The investigated terrestrial impact craters, in general, have relatively lower rim circularity, higher length-to-width ratio, and lower relief, slope, and elevation than volcanic calderas. The topographic characteristics can be explained by geological processes associated with the formation and post-deformation of impact craters. The excavation and ejection by initial impact and rebound of excavated materials contribute to low elevation. The post-impact deformation, including inward collapse and slump of unstable rims, weathering, erosion, and sediment deposition, further reduces elevation and relief and modifies shapes resulting in lower circularity and higher length-to-width ratio. Due to the resolution limitation of the source DEM data and the number of real-world samples, the model has only been validated for craters of 0.88 to 100 km in diameter, which can be generalized to explore undiscovered terrestrial impact craters using cloud computing with global datasets provided by platforms such as Google Earth Engine and Microsoft Planetary Computer.

The 2022 Eruption of Wolf Volcano, Galápagos: The Role of Caldera Ring‐Faults During Magma Transfer From InSAR Deformation Data

AbstractThe 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.

Modelling the post-caldera plumbing system of Changbaishan volcano (China) from integrated geochemical, isotopic, geobarometric, and geophysical data

The Changbaishan Tianchi volcano is the largest active volcano of China and the 946 CE VEI 6 “Millennium eruption” (ME) is responsible for the formation of its summit caldera. After the ME, the volcanic activity is characterized by smaller eruptions most of which are recorded only in historical documents. Here we analyze the deposits of the ME and 1403 CE eruption outcropping in the Heishigou section, which is located east of the Chagbaishan Tianchi caldera. We present whole-rock, glass, mineral, isotopic, and geobarometric data to obtain information on the vertical extent of the caldera and post-caldera plumbing system of the volcano. Results show that the plumbing system of the 1403 CE eruption is located between 1.9 and 5.4 km depth, which is also the depth range of the shallower reservoir of the ME. Therefore, the magma reservoir of the ME was active also after the caldera collapse. The 1403 CE magma stored in the ME reservoir evolved by fractional crystallization of the residual ME trachytic magma and underwent the contamination of upper crust material. The 1403 CE event was triggered by the injection of less evolved, mafic magma from depth. The inferred depth of the shallower ME reservoir and that of the 1403 CE eruption is the same of the present-day main geophysical anomalies (low seismic velocity and low resistivity) recorded below the Changbaishan Tianchi caldera. In addition, the source of the deformations related to the uplift during unrest episodes between 2002 and 2020 suggests a refilling of the 1403 CE plumbing system. We conclude that the plumbing system of the 1403 CE, post-caldera activity at Changbaishan Tianchi volcano is the same of the ME eruption and it is still active. This system may be recharged from the arrival of less evolved magma from depth. The proposed conceptual model of the post-caldera volcanism at Changbaishan Tianchi represents a reference framework to interpret the signals of unrest and decipher the dynamics of the present-day plumbing system.

New geological and geochronological constraints on the evolution of the Cotacachi - Cuicocha volcanic complex (Ecuador)

Extensive fieldwork at the Cotacachi-Cuicocha Volcanic Complex (CCVC, North of Ecuador) resulted in a new collection of geological data including cartography, chronology, petrography, geochemistry, and morphology. This volcanic complex is formed by a central volcano (Cotacachi: 4939 m asl, current bulk volume of 56 ± 4 km3), several peripheral domes, and a 3 km wide volcanic caldera (Cuicocha: 4.2 ± 0.1 km3). CCVC comprises three stratigraphic members: The first, Cotacachi Basal, represents the initial phase of construction, which started at 173 ± 4 ka with a basal andesitic lava flow succession (∼500 m-thick) including isolated basaltic-andesitic lavas (Verde Tola unit; NE: 113 ± 6 ka, SE: 133 ± 9 ka), the construction of some peripheral amphibole-bearing andesitic domes, such as Muyurcu and Loma Negra (138 ± 4 ka and <108 ka, respectively), and a debris-avalanche deposit to the north-west (0.5–1.8 km3, older than 108 ka). The second member, Upper Cotacachi, consists of an andesitic lava flow succession (∼300 m-thick), younger than 108 ± 6 ka. A gap of activity occurred afterwards from 100 to 70 ka, during which a second debris-avalanche (0.2–1.1 km3, 108 to 65 ka) occurred to the NE, followed by the effusion of the dacitic Piribuela dome (65 ± 2 ka). Afterwards, several superimposed andesitic lava flows were emplaced at the summit, possibly around 15–10 ka since they lack glacial erosion. The third member includes the extrusion of the andesitic Cuicocha pre-caldera domes, which marks the beginning of a new eruptive stage of activity of CCVC during the Holocene, resulting in a violent eruption (3525 ± 35 to 2980 ± 30 a BP; VEI = 5) that partially destroyed the young dome and formed a funnel-shaped caldera (Cuicocha Caldera-Lake), ending with the emplacement of the Wolf and Yerovi post-caldera domes.