Dynamic modeling suggests terrace zone asymmetry in the Chicxulub crater is caused by target heterogeneity

We investigate the cause of terrace zone asymmetry in the Chicxulub impact crater using dynamic models of crater formation. Marine seismic data acquired across the crater show that the geometry of the crater's terrace zone, a series of sedimentary megablocks that slumped into the crater from the crater rim, varies significantly around the offshore half of the crater. The seismic data also reveal that, at the time of impact, both the water depth and sediment thickness varied with azimuth around the impact site. To test whether the observed heterogeneity in the pre-impact target might have affected terrace zone geometry we constructed two end-member models of upper-target structure at Chicxulub, based on the seismic data at different azimuths. One model, representing the northwest sector, had no water layer and a 3-km thick sediment layer; the other model, representing the northeast sector, had a 2-km water layer above a 4-km sediment layer. Numerical models of vertical impacts into these two targets produced final craters that differ substantially in terrace zone geometry, suggesting that the initial water depth and sediment thickness variations affected the structure of the terrace zone at Chicxulub. Moreover, the differences in terrace zone geometry between the two numerical models are consistent with the observed differences in the geometry of the terrace zone at different azimuths around the Chicxulub crater. We conclude that asymmetry in the pre-impact target rocks at Chicxulub is likely to be the primary cause of asymmetry in the terrace zone.     

Heating-induced changes in the anisotropy of magnetic susceptibility of impact breccias, Chicxulub Crater (Mexico)

Initial results of a thermal treatment study on the anisotropy of magnetic susceptibility (AMS) of impact breccias from Chicxulub crater are used to investigate the nature of the magnetic fabrics. Chicxulub impact breccias are heterogeneous materials, with carbonate, basement and melt clasts within carbonate-rich or melt-rich matrix. Samples studied come from the carbonate-rich basal unit Lower Suevite in the Yaxcopoil-1 borehole impactite sequence (core depth interval: 885?C895 m). The Lower Suevite is characterized by mixed prolate and oblate ellipsoids with shallow to steep principal susceptibility axes, which had been related to emplacement as an excavation flow with ground-surge components during the early cratering stages. Thermal treatment results in changes in the fabrics with a tendency to oblate fabrics. Stepwise thermal treatment up to 700°C reveals different behaviors for the oblate, neutral and prolate fabrics marked by changes in AMS parameters and principal susceptibility axis orientations. A sample with oblate fabrics and vertical minimum axes showed an increase of magnetic susceptibility at high temperatures, indicating formation of secondary magnetite and fabric enhancement. A sample with neutral ellipsoid showed heating-induced changes towards oblate fabrics and vertical minimum susceptibility axes. Samples characterized by prolate ellipsoids with horizontal maximum axes showed no directional changes. In a sample with apparent intermediate or inverse fabrics, vertical maximum axes showed changes to horizontal inclinations, with the intermediate and maximum axes switching positions. Changes induced by stepwise thermal treatment appear useful to characterize the fabrics of impact lithologies. Further investigation of heating-induced effects in mineralogy, grain size and textural changes is, however, required to relate the different behaviors observed after stepwise thermal treatment with the magnetic mineralogy and emplacement mode of the breccias.     

小行星撞击对地球的上地幔对流的影响

本文分别基于数值结果和地质学模型,在假定地球的上地幔存在稳定Rayleigh-Bénard对流的基础上,模拟了直径为10 km（陨石坑直径约180 km,以Chicxulub为例）和直径为100 km（陨石坑直径约1000 km）的小行星撞击对地球的上地幔对流格局的影响.本文将直径10 km小行星的撞击效果等效为热异常,将直径100 km小行星的撞击效果等效为热异常和速度异常（主要指陨石坑底部的回弹）的叠加.计算结果表明,当小行星的直径在10 km左右时,撞击对上地幔对流的影响十分微弱,热扰动时间仅2—3 Ma;而当小行星的直径达到100 km时,撞击就会对上地幔对流产生强烈影响.这时,对流从扰动到新的稳态有一定模式可循（依次为：调整、多个对流环、调整、稳定）,扰动的持续时间受黏度和撞击点位置影响,同时稳定后地幔热柱会向着撞击点的方向产生一定的位移.     

Numerical model of crater lake eruptions

We present results from a numerical investigation of subaqueous eruptions involving superheated steam released through a lake mimicking the volcanic setting at Mt. Ruapehu. The simulations were conducted using an adaptive mesh, multi-material, hydrodynamics code with thermal conduction SAGE, (Simple Adaptive Grid Eulerian). Parameters investigated include eruption pressure, lake level and mass of superheated vapor. The simulations produced a spectrum of eruption styles from vapor cavities to radial jets that resulted in hazards that ranged from small-scale waves to high amplitude surges that reached and cascaded over the edge of the crater rim. There was an overall tendency for lake surface activity to increase (including wave amplitude) with increasing mass of superheated vapor and eruption pressure. Surface waves were induced by the formation and collapse of a gas cavity. The collapse of the cavity is considered to play a major role in the characteristic features observed during a subaqueous eruption. The additional mass of superheated vapor produced a larger cavity that displaced a larger area of the lake surface resulting in fast moving surges upon the collapse of the cavity. High lake levels (>90 m) appear to suppress the development of explosive jetting activity when eruption pressures are <10 MPa. At very large eruption pressures (>10 MPa), vertical jets and radial ejections of steam and water can occur in water depths >90 m. Less explosive eruption styles can produce hazardous events such as lahars by the outward movement of surface waves over the crater rim.     

A preparation zone for volcanic explosions beneath Naka-dake crater,Aso volcano,as inferred from magnetotelluric surveys

The 1st crater of Naka-dake, Aso volcano, is one of the most active craters in Japan, and known to have a characteristic cycle of activity that consists of the formation of a crater lake, drying-up of the lake water, and finally a Strombolian-type eruption. Recent observations indicate an increase in eruptive activity including a decrease in the level of the lake water, mud eruptions, and red hot glows on the crater wall. Temporal variations in the geomagnetic field observed around the craters of Naka-dake also indicate that thermal demagnetization of the subsurface rocks has been occurring in shallow subsurface areas around the 1st crater. Volcanic explosions act to release the energy transferred from magma or volcanic fluids. Measurement of the subsurface electrical resistivity is a promising method in investigating the shallow structure of the volcanic edifices, where energy from various sources accumulates, and in investigating the behaviors of magma and volcanic fluids. We carried out audio-frequency magnetotelluric surveys around the craters of Naka-dake in 2004 and 2005 to determine the detailed electrical structure down to a depth of around 1 km. The main objective of this study is to identify the specific subsurface structure that acts to store energy as a preparation zone for volcanic eruption. Two-dimensional inversions were applied to four profiles across the craters, revealing a strongly conductive zone at several hundred meters depth beneath the 1st crater and surrounding area. In contrast, we found no such remarkable conductor at shallow depths beneath the 4th crater, which has been inactive for 70 years, finding instead a relatively resistive body. The distribution of the rotational invariant of the magnetotelluric impedance tensor is consistent with the inversion results. This unusual shallow structure probably reflects the existence of a supply path of high-temperature volcanic gases to the crater bottom. We propose that the upper part of the conductor identified beneath the 1st crater is mainly composed of hydrothermally altered zone that acts both as a cap to upwelling fluids supplied from deep-level magma and as a floor to infiltrating fluid from the crater lake. The relatively resistive body found beneath the 4th crater represents consolidated magma. These results suggest that the shallow conductor beneath the active crater is closely related to a component of the mechanism that controls volcanic activity within Naka-dake.     

Magnetic susceptibility logging of Chicxulub proximal impact breccias in the Santa Elena borehole: implications for emplacement mode

Magnetic susceptibility logging is used to study the impact breccias in the Chicxulub crater. The basic premise is that the high contrasts in magnetic properties can be used to characterize the breccias. The Santa Elena borehole was drilled 110 km radial distance from crater center and sampled a 172 m thick sequence of impact breccias, between 332 and 504 m depth. Breccia units are distinguished from differences in composition, size, and relative contents of clasts, type of matrix and textural and lithological assemblages, which can be resolved in the susceptibility logs. The whole-core log shows characteristic variation patterns with high, intermediate and low susceptibilities. High resolution logging of matrix and clasts records the heterogeneous nature of impactites, with higher variability at smaller spatial scales. Measurements confirm that diamagnetic susceptibilities characterize the carbonate clasts, high susceptibilities the basement granitic clasts and intermediate values the silicate melt-rich and silicate-poor matrix. Intermediate variable susceptibilities characterize breccias rich in melt particles. Correlation of matrix and clast logs with whole-core log shows that signal is controlled by the matrix. Logs for clast shows a discrete distribution with peaks of intermediate to high values, which correlate with large clast distributions. The ejecta blanket includes the fallback suevites rich in silicate melt particles and shocked minerals, the high temperature vapor deposits from ejecta curtain collapse and high velocity basal flows, and the carbonate rich deposits from lateral basal flows and secondary cratering. Late fallback suevites record minor turbulent conditions resulting from progressive cooling of the ejecta plume.     

Groundwater transport of crater-lake brine at Poa´s Volcano, Costa Rica

Poa´s Volcano is an active stratovolcano in Costa Rica that has a lake in its active crater. The crater lake has high temperatures (50–90 °C), high acidity (pH ≈ 0.0), and a high dissolved-solids content (100 g/kg). The volcano has numerous freshwater springs on its flanks, but a few on the northwestern flank are highly acidic (pH = 1.6–2.5) and have high dissolved-solids concentrations (2–22 g/kg). This study analyzes the regional groundwater system at Poa´s and demonstrates the likelihood that the water discharging from the acidic springs in the Rio Agrio watershed originates at the acidic crater lake. Both heat and solute transport are analyzed on a regional scale through numerical simulations using the HST3D finite-difference model, which solves the coupled equations for fluid flow, heat transport, and solute transport. The code allows fluid viscosity and density to be functions of both temperature and solute concentration. The simulations use estimates for recharge to the mountain and a range of values and various distributions of permeability and porosity. Several sensitivity analyses are performed to test how the uncertainty in many of the model parameters affects the simulation results. These uncertainties yield an estimated range of travel times from the crater lake to the Rio Agrio springs of 1–30 years, which is in close agreement with the results of tritium analyses of the springs. Calculated groundwater fluxes into and out of the crater lake are both about several hundred kg/s. These fluxes must be accounted for in water budgets of the crater lake.     

Calculations of Asteroid Impacts into Deep and Shallow Water

Contrary to received opinion, ocean impacts of small (<500?m) asteroids do not produce tsunamis that lead to world-wide devastation. In fact the most dangerous features of ocean impacts, just as for land impacts, are the atmospheric effects. We present illustrative hydrodynamic calculations of impacts into both deep and shallow seas, and draw conclusions from a parameter study in which the size of the impactor and the depth of the sea are varied independently. For vertical impacts at 20?km/s, craters in the seafloor are produced when the water depth is less than about 5?C7 times the asteroid diameter. Both the depth and the diameter of the transient crater scale with the asteroid diameter, so the volume of water excavated scales with the asteroid volume. About a third of the crater volume is vaporised, because the kinetic energy per unit mass of the asteroid is much larger than the latent heat of vaporisation of water. The vaporised water carries away a considerable fraction of the impact energy in an explosively expanding blast wave which is responsible for devastating local effects and may affect worldwide climate. Of the remaining energy, a substantial portion is used in the crown splash and the rebound jet that forms as the transient crater collapses. The collapse and rebound cycle leads to a propagating wave with a wavelength considerably shorter than classical tsunamis, being only about twice the diameter of the transient crater. Propagation of this wave is hindered somewhat because its amplitude is so large that it breaks in deep water and is strongly affected by the blast wave??s perturbation of the atmosphere. Even if propagation were perfect, however, the volume of water delivered per metre of shoreline is less than was delivered by the Boxing Day 2004 tsunami for any impactor smaller than 500?m diameter in an ocean of 5?km depth or less. Near-field effects are dangerous for impactors of diameter 200?m or greater; hurricane-force winds can extend tens of kilometers from the impact point, and fallout from the initial splash can be extremely violent. There is some indication that near-field effects are more severe if the impact occurs in shallow water.     

A telemetering system for monitoring aqueous polythionates in an active crater lake

We have developed an in-situ monitoring system for aqueous polythionates in the Yugama crater lake at Kusatsu-Shirane volcano. A commercially available nitrate-selective electrode has a Nernstian response to the concentration of polythionates in solution in the range of pH 2.5–11. However, many active crater lakes contain highly acidic water (pH ranges from 0.2 to 2). The monitoring system reduces the acidity of the lake water to pH 4 without changing the concentration of polythionates using a high efficiency cation-exchange membrane technique. The monitoring system was placed in a submersed station, because the Yugama surface lake water freezes in mid-winter. The system is described and preliminary data are reported.     

The summit hydrothermal system of Stromboli. New insights from self-potential,temperature, CO<Subscript>2</Subscript> and fumarolic fluid measurements,with structural and monitoring implications

Accurate and precisely located self-potential (SP), temperature (T) and CO₂ measurements were carried out in the summit area of Stromboli along 72 straight profiles. SP data were acquired every metre and T data every 2.5 m. CO₂ concentrations were acquired with the same density as T, but only along seven profiles. The high density of data and the diversity of the measured parameters allows us to study structures and phenomena at a scale rarely investigated. The shallow summit hydrothermal activity (Pizzo–Fossa area) is indicated by large positive SP, T and CO₂ anomalies. These anomalies are focused on crater faults, suggesting that the fracture zones are more permeable than surrounding rocks at Stromboli. The analysis of the distribution of these linear anomalies, coupled with the examination of the geologic, photographic and topographic data, has led us to propose a new structural interpretation of the summit of Stromboli. This newly defined structural framework comprises (1) a large Pizzo circular crater, about 350 m in diameter; (2) a complex of two concealed craters nested within the Pizzo crater (the Large and the Small Fossa craters), thought to have formed during the eruption of the Pizzo pyroclastites unit; the Small Fossa crater is filled with highly impermeable material that totally impedes the upward flow of hydrothermal fluids; and (3) The present complex of active craters. On the floor of the Fossa, short wavelength SP lows are organized in drainage-like networks diverging from the main thermal anomalies and converging toward the topographic low in the Fossa area, inside the Small Fossa crater. They are interpreted as the subsurface downhill flow of water condensed above the thermal anomalies. We suspect that water accumulates below the Small Fossa crater as a perched water body, representing a high threat of strong phreatic and phreatomagmatic paroxysms. T and CO₂ anomalies are highly correlated. The two types of anomalies have very similar shapes, but the sensitivity of CO₂ measurements seems higher for lowest hydrothermal flux. Above T anomalies, a pronounced high frequency SP signal is observed. Isotopic analyses of the fluids show similar compositions between the gases rising through the faults of the Pizzo and Large Fossa craters. This suggests a common origin for gases emerging along different structural paths within the summit of Stromboli. A site was found along the Large Fossa crater fault where high gas flux and low air contamination made gas monitoring possible near the active vents using the alkaline bottle sampling technique.     

Chicxulub impact ejecta from Albion Island, Belize

Impact ejecta from the Albion Formation are exposed in northern Belize. The ejecta come from the outer portion of the continuous ejecta blanket of the Chicxulub crater, which is located 360 km to the northwest. The basal unit of the Albion Formation is a 1-m-thick clay and dolomite spheroid bed composed of up to four discrete flows. The clay spheroids are altered impact glass, and the dolomite spheroids are accretionary lapilli. The upper unit is a 15-m-thick coarse diamictite bed containing altered glass, large accretionary blocks, striated, polished, and impacted cobbles, and rare shocked quartz. The abundance of accretionary clasts, evidence for atmospheric drag sorting, and the presence of multiple flows in the Albion Formation indicate that atmospheres play an important role in the formation of the outer portions of continuous ejecta blankets of large craters.     

美国Chesapeake湾撞击坑研究进展及其启示

位于美国弗吉尼亚东海岸直径85 km的Chesapeake湾撞击坑,是十几年前发现的由一颗陨星撞击形成的一个复杂撞击坑. 该坑的研究经历四个阶段：地下水调查、撞击坑的发现、美国多学科多部门的综合研究和即将进行的国际钻探取心项目. 钻井岩心中的角砾成份和微体化石,提示撞击坑的存在,并确定撞击发生在35 Ma前,即始新始晚期. 地震反射剖面资料帮助寻找到撞击坑的具体位置,确定撞击坑的结构和形态特征. Chesapeake湾撞击坑埋藏在新生界沉积层之下,是全球已知最大的、保存最好的撞击坑之一. Chesapeake湾撞击坑主要形态像一顶倒置的宽边大草帽,包括外缘、环状洼地、峰环（内缘）、内盆和中央峰. 撞击坑的形成破坏了原来的含水层,撞击坑当时即被富含咸水的抛射角砾岩和海啸角砾岩充填,再被后来的沉积层覆盖. Chesapeake湾撞击坑导致地面沉降、河流变向、海岸含水层的中断、内陆咸水楔的出现、地震,决定Chesapeake湾本身的位置,至今仍然影响当地居民的生活. 了解Chesapeake湾撞击坑对我国撞击坑研究具有借鉴作用.     

Crater formation in soils by raindrop impact

The process of crater formation by the impact of water drops on soil, sand and various other target material was studied. Craters of various shapes and sizes were observed on different target materials or conditions, ranging from circumferential depression to completely hemispherical shape. Crater shape was dependent upon target material, its ?ow stress or shear strength and the presence and thickness of water on the surface. Between 5 and 22 per cent of impact energy was spent on cratering, but the relationship between crater volume and kinetic energy of a raindrop was curvilinear, indicating a lower ef?ciency of impact energy in removing target material as the energy increases. Impact impulse, on the other hand, showed a more linear relationship with crater volume, and the ratio of impulse over crater volume (I/V) remained constant for the entire range of drop sizes, impact velocities, and surface conditions used in this study. Surface shear strength, represented by the penetration depth of fall‐cone penetrometer, appeared to be a key factor involved in this process. An equation was developed which related crater volume to cone penetration depth and impact impulse. Crater volume, which appeared to be a better indicator of the total amount of material dislodged by a raindrop than splash amount, can thus be predicted using this equation. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrochemical dynamics of the “lake–spring” system in the crater of El Chichón volcano (Chiapas,Mexico)

El Chichón volcano (Chiapas, Mexico) erupted violently in March–April 1982, breaching through the former volcano–hydrothermal system. Since then, the 1982 crater has hosted a shallow (1–3.3 m, acidic (pH ∼ 2.2) and warm (∼ 30 °C) crater lake with a strongly varying chemistry (Cl/SO₄ = 0–79 molar ratio). The changes in crater lake chemistry and volume are not systematically related to the seasonal variation of rainfall, but rather to the activity of near-neutral geyser-like springs in the crater (Soap Pool). These Soap Pool springs are the only sources of Cl for the lake. Their geyser-like behaviour with a long-term (months to years) periodicity is due to a specific geometry of the shallow boiling aquifer beneath the lake, which is the remnant of the 1983 Cl-rich (24,000 mg/l) crater lake water. The Soap Pool springs decreased in Cl content over time. The zero-time extrapolation (1982, year of the eruption) approaches the Cl content in the initial crater lake, meanwhile the extrapolation towards the future indicates a zero-Cl content by 2009 ± 1. This particular situation offers the opportunity to calculate mass balance and Cl budget to quantify the lake–spring system in the El Chichón crater. These calculations show that the water balance without the input of SP springs is negative, implying that the lake should disappear during the dry season. The isotopic composition of lake waters (δD and δ¹⁸O) coincide with this crater lake-SP dynamics, reflecting evaporation processes and mixing with SP geyser and meteoric water. Future dome growth, not observed yet in the post-1982 El Chichón crater, may be anticipated by changes in lake chemistry and dynamics.     

The time-dependent variation of hydrogeochemical characteristics, thermal regime, and biocenoses in the fresh and thermal waters of the Lake Karymskii basin following the catastrophic underwater eruption of 1996 in the Akademii Nauk Caldera, Kamchatka

The results of biohydrogeochemical monitoring are used to study time-dependent variations in the hydrogeochemical characteristics of the Lake Karymskii water mass, the state and characteristics of underwater discharge zones in the Tokarev crater (formed in 1996), hydrogeochemical characteristics of thermal springs around the lake, and the biota succession in the lake for the period 1996–2006 (2007). We detected a stratification in the chemical composition of water over depth and the presence of persistent zones of increased concentrations of dissolved oxygen. We found an alkalization of lake water and a decrease in its total salinity. The new thermal springs and underwater discharges of thermal water and gases were found to be continuing. The first data were obtained on the concentration of microelements in the thermal springs of the Karymskii basin. The biodiversity of algae in Lake Karymskii was largely increased by the species diversity of benthic Bacillariophyta. The plankton phytocomponent of the precatastrophic period was found to have been regenerated in the lake as of April 2007.     

Brownish discoloration of the summit crater lake of Mt. Shinmoe-dake,Kirishima Volcano,Japan: volcanic–microbial coupled origin

A drastic change in lake water color from blue-green to brown was observed in the summit crater lake of Mt. Shinmoe-dake, Kirishima Volcano about 8 months after its 2008 eruption. The color change lasted for about 2 months (April–June 2009). The discoloration was attributed to a brownish color suspension that had formed in the lake water. X-ray fluorescence and Fourier transform infrared analyses of a sample of the suspension identified schwertmannite (Fe₈O₈(OH)₆(SO₄)). A cultivation test of iron-oxidizing bacteria for the sampled lake water with lakebed sediment revealed that the crater lake hosts iron-oxidizing bacteria, which likely participated in schwertmannite formation. We suggest that pyrite (FeS₂) provided an energy source for the iron-oxidizing bacteria since the mineral was identified in hydrothermally altered tephra ejected by the August 2008 eruption. From consideration of these and other factors, the brownish discoloration of the summit crater lake of Mt. Shinmoe-dake was inferred to have resulted from a combined volcanic–microbial process.     

Evolution and hydrological conditions of a maar volcano (Atexcac crater,Eastern Mexico)

The Atexcac maar is located in the central part of the Serdán–Oriental lacustrine/playa basin in the eastern Mexican Volcanic Belt. It is part of a dispersed and isolated monogenetic field consisting of maar volcanoes, basaltic cinder cones and rhyolitic domes. Atexac is a maar volcano excavated into pyroclastic deposits, basaltic lava flows and the flanks of a cinder cone cluster, which itself was built on a topographic high consisting of limestone. It has an ENE-trending elliptical shape with beds, mostly unconsolidated deposits that dip outward at 16–22°. The Atexcac crater was formed from vigorous phreatomagmatic explosions in which fluctuations in the availability of external water, temporal migration of the locus of the explosion, and periodic injection of new magma were important controls on the evolution of the maar crater. Variations in grain sizes and component proportions of correlated deposits from the different sections suggest a migration of the locus of explosions, producing different eruptive conditions with fluctuating water–magma interactions. Deposits rich in large intrusive and limestone blocks are associated with a matrix enriched in small andesitic lapilli. This could suggest differential degrees of fragmentation due to inherited (previously acquired) fragmentation and/or relative distance to the locus of explosions. Initial short-lived phreatic explosions started at the southwest part of the crater and were followed by an ephemeral vertical column and the influx of external water that led to relatively shallow explosive interactions with the ascending basaltic magma. Drier explosions progressed downward and/or laterally northward, sampling subsurface rock types, particularly intrusive, limestone and andesitic zones as well as localized altered zones (N-NE), caused by repetitive injection of basaltic magma. A final explosive phase involved a new injection of magma and a new influx of external water producing wetter conditions at the end of the maar formation. We infer the aquifer was formed by fractured rocks, predominantly andesitic lava flows and limestone rocks. Andesitic accessory clasts dominate in all stratigraphic levels but these rocks are not exposed in the nearby area. These local hydrogeological conditions contrast with those at nearby maar volcanoes, where the water for the magma/water interactions apparently mostly came from a dominantly unconsolidated tuffaceous aquifer, producing tuff rings with a much lower profile than Atexcac.     

Disappearance of a crater lake: implications for potential explosivity at Soufrière volcano,St Vincent,Lesser Antilles

Soufrière volcano in St Vincent, West Indies, is one of the most active volcanoes in the Eastern Caribbean with at least six eruptions since 1718 AD, the latest of which occurred in 1979. Prior to the 1979 eruption, the active crater hosted deep-water lakes during periods of repose, which were always replenished within a few years after the eruptions. In 1979, the crater was filled with 10⁸ m³ of fragmental material and, despite constant precipitation, has remained virtually dry ever since, with the exception of a small shallow pond. A resistivity survey was conducted in July 2006 to investigate groundwater occurrence in the crater. Results from the resistivity data inversion on several 2-D profiles show a shallow horizontal conductor across the crater floor, consistent with a water-saturated aquifer. They also show that the post-1979 pond, currently present in the crater lake is in fact an outcropping part of the groundwater water reservoir. The reservoir water table is ∼28 m above the pre-1979 lake level and reflects mass equilibrium in the system where constant seepage underground balances the meteoric recharge. We suggest that the groundwater body extends at depth to the bottom of the pre-1979 crater lake, either due to a significant structural discontinuity or because of a reduction of permeability at depth. The estimated maximum volume of water stored underground is 10–30 × 10⁶ m³ and energy considerations indicate that 2.4–7.3 × 10¹⁰ kg of magma would potentially be sufficient to vaporise the whole groundwater body. This amount of magma represents only 13–41% of the mass erupted during the last eruption in 1979 which was the smallest of the past 3 eruptions (1902, 1971–72, 1979). Since explosive phreatic or phreatomagmatic eruptions at Soufrière seem to be linked to magma-water interaction within confined space, the results from this survey suggests that phreatic or phreatomagmatic activity is a distinct possibility during future magma intrusion in the summit area, despite the apparent disappearance of water in the summit crater.     

Modeling Groundwater Inflow to the New Crater Lake at Kīlauea Volcano,Hawai'i

During the 2018 eruption of Kīlauea Volcano, Hawai'i, scientists relied heavily on a conceptual model of explosive eruptions triggered when lava-lake levels drop below the water table. Numerical modeling of multiphase groundwater flow and heat transport revealed that, contrary to expectations, liquid water inflow to the drained magma conduit would likely be delayed by months to years, owing to the inability of liquid water to transit a zone of very hot rock. The summit of Kīlauea subsequently experienced an ∼2-month period of consistent repeated collapses, and the crater now extends below the equilibrium position of the water table. Liquid water first emerged into the deepened crater in late July 2019. The timing of first appearance of liquid water (about 14 months postcollapse) and the rate of crater lake filling (currently ∼27 kg/s) were well-predicted by the numerical modeling done in late spring 2018, which forecast liquid inflow after 3 to 24 months at rates of 10 to 100 kg/s. A second-generation groundwater model, reflecting the current crater geometry, forecasts lake filling over the next several years. The successful 2018 to present forecasts with both models are based on unadjusted in situ permeability estimates (1 to 6 × 10⁻¹⁴ m²) and water-table elevations (600 to 800 m) from a nearby research drillhole and geophysical surveys. Important unknowns that affect the reliability of longer-term forecasts include the equilibrium water-table geometry, the rate of evaporation from the hot and growing crater lake (currently ∼29,000 m² at 70-80 °C), and heterogenous permeability changes caused by the 2018 collapse.     

Subaqueous geothermal activity revealed by lacustrine sediments of the acidic Nakadake crater lake, Aso Volcano, Japan

Lacustrine sediments were sampled from the inaccessible acidic (pH = 0.43) Nakadake crater lake of Aso Volcano, Japan by a simple method. The sediments contain an extremely high content (74 wt.%) of sulfur, which exits as elemental sulfur, gypsum and anhydrite. The abundant elemental sulfur is likely formed by the reaction of SO₂ and H₂S gases and by the SO₂ disproportionation reaction in magmatic hydrothermal system below the crater lake. Based on the sulfur content of sediments and measurements of elevation change of the crater bottom, the sulfur accumulation rate at the Nakadake crater lake was calculated as 250 tonne/day, which is comparable with the SO₂ emission rate (200–600 tonne/day) from the Nakadake crater. The sediments include a small amount (9%) of clear glass shards that are apparently not altered in spite of the high reactivity of hyperacid lake water. This finding suggests that the clear glass shards are fragments of recently emitted magmas from fumaroles on the bottom of the crater lake and the magma emissions continuously occur even in quiescent periods.