Patterns of seismic swarm activity in the Corinth Rift in 2000–2005

Based on the data of the detailed earthquake catalog provided on the website of the Corinth Rift Laboratory, zones of swarm activity are revealed and the variations in the statistical parameters of seismic swarms that occurred in the western part of the Gulf of Corinth are calculated. The preliminary analysis of the catalogue is carried out; the magnitude of completeness and the accuracy of the location of the earth-quake are estimated; the changes in these parameters associated with the development of the observational network are assessed. The b-value (b-values) and the cluster dimension of the set of hypocenters are estimated, and time variations in these parameters in the course of the evolution of swarm activity are revealed. The style of changes in the parameters characterizing the seismic regime during intervals of swarm activity indicates that the process of failure exhibits scale redistribution over the course of time, changing from upscaling (progression from smaller to larger scales) at the stage of increasing seismicity to downscaling (progression from larger to lower scales) at the stage of decay. These particular features of enhancement and reduction of swarm seismicity are qualitatively similar to the scenarios of source preparation and aftershock relaxation of strong earthquakes. The pattern of variations of the swarm seismicity studied is similar to those identified in the previous laboratory and field modeling of various transient modes of seismicity. This fact confirms the relevancy of the retrieved results and conclusions based on the laboratory studies of transient modes, and suggests that the latter have a universal governing mechanism.     

Recent explosive eruptions and volcano hazards at Soputan volcano—a basalt stratovolcano in north Sulawesi, Indonesia

Soputan is a high-alumina basalt stratovolcano located in the active North Sulawesi-Sangihe Islands magmatic arc. Although immediately adjacent to the still geothermally active Quaternary Tondono Caldera, Soputan’s magmas are geochemically distinct from those of the caldera and from other magmas in the arc. Unusual for a basalt volcano, Soputan produces summit lava domes and explosive eruptions with high-altitude ash plumes and pyroclastic flows—eight explosive eruptions during the period 2003–2011. Our field observations, remote sensing, gas emission, seismic, and petrologic analyses indicate that Soputan is an open-vent-type volcano that taps basalt magma derived from the arc-mantle wedge, accumulated and fractionated in a deep-crustal reservoir and transported slowly or staged at shallow levels prior to eruption. A combination of high phenocryst content, extensive microlite crystallization and separation of a gas phase at shallow levels results in a highly viscous basalt magma and explosive eruptive style. The open-vent structure and frequent eruptions indicate that Soputan will likely erupt again in the next decade, perhaps repeatedly. Explosive eruptions in the Volcano Explosivity Index (VEI) 2–3 range and lava dome growth are most probable, with a small chance of larger VEI 4 eruptions. A rapid ramp up in seismicity preceding the recent eruptions suggests that future eruptions may have no more than a few days of seismic warning. Risk to population in the region is currently greatest for villages located on the southern and western flanks of the volcano where flow deposits are directed by topography. In addition, Soputan’s explosive eruptions produce high-altitude ash clouds that pose a risk to air traffic in the region.     

Geochemical model of a magmatic–hydrothermal system at the Lastarria volcano, northern Chile

Lastarria volcano (25°10′ S, 68°31′ W; 5,697 m above sea level), located in the Central Andes Volcanic Zone (northern Chile), is characterized by four distinct fumarolic fields with outlet temperatures ranging between 80°C and 408°C as measured between May 2006–March 2008 and April–June 2009. Fumarolic gasses contain significant concentrations of high temperature gas compounds (i.e., SO₂, HCl, HF, H₂, and CO), and isotopic ratios (³He/⁴He, δ¹³C–CO₂, δ¹⁸O–H₂O, and δD–H₂O) diagnostic of magmatic gas sources. Gas equilibria systematics, in both the H₂O-H₂-CO₂-CO-CH₄ and alkane–alkene C₃ system, suggest that Lastarria fumarolic gasses emanate from a superheated vapor that is later cooled and condensed at relatively shallow depths. This two-stage process inhibits the formation of a continuous aquifer (e.g., horizontal liquid layer) at relatively shallow depth. Recent developments in the magmatic gas system may have enhanced the transfer and release of heat causing shallow aquifer vaporization. The consequent pressure increase and aquifer vaporization likely triggered the inflation events beginning in 2003 at the Lastarria volcano.     

Quantification of the 1998–1999 explosion sequence at Popocatépetl volcano,Mexico

The sequence of large Vulcanian explosions occurring at the andesitic Popocatépetl volcano, Mexico during November 1998 to April 1999 was studied. The size of 26 largest explosions was estimated from broadband seismic records at the distance of 4 km from the crater. The sequence began with the largest explosion (E = 2.6 × 10¹² J) occurring on 25 November at 08:05, and following largest daily explosions were characterized by gradual decrease in the energy. The energy of 20 large (E ≥ 10¹¹ J) explosions was distributed as Student's t-distribution with a geometrical mean Log E = 11.81 (J).     

Relationship between tilt changes and effusive-explosive episodes at an andesitic volcano: the 2004–2005 eruption at Volcán de Colima,México

Continuous tilt changes during the 2004–2005 effusive-explosive episodes at Volcán de Colima (México) were recorded simultaneously by two tiltmeters installed on opposite sides of the volcano at elevations of 2200 m and 3060 m above sea level. Data indicate that the 2004 lava extrusion was preceded by an inflation accompanied by a deflation. The 2005 explosion sequences were associated with a deflationary–inflationary tilt. The period between the 2004 extrusion and the 2005 explosions was characterized by an inflationary tilt during a 3 month period. Two deformation sources were located. The first was situated at a depth between 300 m and 1800 m beneath the crater at the northern flank of the volcano and was responsible for volcano deformation during the preliminary September 2004 stage, the October 2004 extrusion, and the initial stage of the transition period and the March 2005 explosion sequence. The second source was located at a depth between 1800 m and 2800 m beneath the crater at the southern flank of the volcano and was responsible for volcano deformation during the final stage of the transition period and the May–June 2005 explosion sequence.     

The February–July 2007 eruption of the summit crater of Klyuchevskoi volcano, Kamchatka

Observations of the summit eruption of Klyuchevskoi volcano in the period from February 15, 2007 to July 9, 2007 are considered. This typical (for this volcano) summit eruption was explosive-effusive in character. The ejectamenta volume is estimated at 0.025 km³. Calculation of active phases of the volcano was carried out in accordance with V.A. Shirokov’s technique. The identified active phases agree well with the eruptive periods. The 2007 summit eruption corresponds to an active phase (May 2006 to May 2009) favorable for the volcano’s eruption. Geodetic observations carried out since 1979 along a radial profile have revealed uplifts and subsidences of the northeastern slope of the volcano. The maximum displacement of 23 cm was recorded in 2007 on the site closest to the volcano crater at a distance of 11 km from the summit crater center. In the course of two previous summit eruptions (2003–2004 and 2005) insignificant uplifts and subsidences of the slope were also noted, although the general ascent of the slope remained. This indicated possible repeated eruptions in the nearest future. Changes in the seismicity before, during and after the eruption are also discussed.     

Magma-hydrothermal system interaction inferred from volcanic gas measurements obtained during 2003–2008 at Meakandake volcano,Hokkaido, Japan

Phreatic eruptions occurred at the Meakandake volcano in 1988, 1996, 1998, 2006, and 2008. We conducted geochemical surveillance that included measurements of temperature, SO₂ emission rates, and volcanic gas composition from 2003 to 2008 at the Nakamachineshiri (NM), Northwest (NW), and Akanuma (AK) fumarolic areas, and the 96–1 vent, where historical eruptions had occurred. The elemental compositions of the gases discharged from the different areas are similar compared with the large variations observed in volcanic gases discharged from subduction zones. All the gases showed high apparent equilibrium temperatures, suggesting that all these gases originated from a common magmatic gas. The gases discharged from each area also exhibited different characteristics, which are probably the results of differences in the conditions of meteoric water mixing, quenching of chemical reactions, and vapor-liquid separation. The highest apparent equilibrium temperatures (about 500°C) were observed in the case of NW fumarolic gases, despite the low outlet temperature of about 100°C at these fumaroles. Since the NW fumaroles were formed as a result of the 2006 phreatic eruption, the high-temperature gas supply to the NW fumarole suggests that the phreatic eruption was caused by the ascent of high-temperature magmatic gases. The temperatures, compositions, and emission rates of the NM and 96–1 gases did not show any appreciable change after the 2006 eruption, indicating that each fumarolic system had a separate magmatic-hydrothermal system. The temperatures, compositions, and emission rates of the NM fumarolic gases were apparently constant, and these fumaroles are inferred to be formed by the evaporation of a hydrothermal system with a constant temperature of about 300°C. The 96–1 gas compositions showed large changes during continuous temperature decrease from 390° to 190°C occurred from 2003 to 2008, but the sulfur gas emission rates were almost constant at about four tons/day. At the 96–1 vent, the SO₂/H₂S ratio decreased, while the H₂/H₂O ratio remained almost constant; this was probably caused by the rock-buffer controlled chemical reaction during the temperature decrease.     

Multiple volcano deformation sources in a post-rifting period: 1989–2005 behaviour of Krafla,Iceland constrained by levelling,tilt and GPS observations

The Krafla rifting episode, which occurred in North Iceland in 1975–1984, was followed by inflation of a shallow magma chamber until 1989. At that time, gradual subsidence began above the magma chamber and has continued to the present at a declining rate. Pressure decrease in a shallow magma chamber is not the only source of deformation at Krafla, as other deformation processes are driven by exploitation of two geothermal fields, together with plate spreading. In addition, deep-seated magma accumulation appears to take place, with its centre ∼ 10 km north of the Krafla caldera. The relative strength of these sources has varied with time. New results from a levelling survey and GPS measurements in 2005 allow an updated view on the deformation field. Deformation rates spanning 2000–2005 are the lowest recorded in the 30-year history of geodetic studies at the volcano. The inferred rate of 2000–2005 subsidence related to processes in the shallow magma chamber is less than 0.3 cm/yr whereas it was ∼ 5 cm/yr in 1989–1992. Currently, the highest rate of subsidence takes place in the Leirbotnar area, within the Krafla caldera, and appears to be a result of geothermal exploitation.     

The 2005 eruption of Sierra Negra volcano,Galápagos,Ecuador

Sierra Negra volcano began erupting on 22 October 2005, after a repose of 26 years. A plume of ash and steam more than 13 km high accompanied the initial phase of the eruption and was quickly followed by a ~2-km-long curtain of lava fountains. The eruptive fissure opened inside the north rim of the caldera, on the opposite side of the caldera from an active fault system that experienced an m_b 4.6 earthquake and ~84 cm of uplift on 16 April 2005. The main products of the eruption were an `a`a flow that ponded in the caldera and clastigenic lavas that flowed down the north flank. The `a`a flow grew in an unusual way. Once it had established most of its aerial extent, the interior of the flow was fed via a perched lava pond, causing inflation of the `a`a. This pressurized fluid interior then fed pahoehoe breakouts along the margins of the flow, many of which were subsequently overridden by `a`a, as the crust slowly spread from the center of the pond and tumbled over the pahoehoe. The curtain of lava fountains coalesced with time, and by day 4, only one vent was erupting. The effusion rate slowed from day 7 until the eruption’s end two days later on 30 October. Although the caldera floor had inflated by ~5 m since 1992, and the rate of inflation had accelerated since 2003, there was no transient deformation in the hours or days before the eruption. During the 8 days of the eruption, GPS and InSAR data show that the caldera floor deflated ~5 m, and the volcano contracted horizontally ~6 m. The total eruptive volume is estimated as being ~150×10⁶ m³. The opening-phase tephra is more evolved than the eruptive products that followed. The compositional variation of tephra and lava sampled over the course of the eruption is attributed to eruption from a zoned sill that lies 2.1 km beneath the caldera floor.     

The rhyolitic–andesitic eruptive history of Cotopaxi volcano,Ecuador

At Cotopaxi volcano, Ecuador, rhyolitic and andesitic bimodal magmatism has occurred periodically during the past 0.5 Ma. The sequential eruption of rhyolitic (70–75% SiO₂) and andesitic (56–62% SiO₂) magmas from the same volcanic vent over short time spans and without significant intermingling is characteristic of Cotopaxi’s Holocene behavior. This study documents the eruptive history of Cotopaxi volcano, presenting its stratigraphy and geologic field relations, along with the relevant mineralogical and chemical nature of the eruptive products, in order to determine the temporal and spatial relations of this bimodal alternation. Cotopaxi’s history begins with the Barrancas rhyolite series, dominated by pumiceous ash flows and regional ash falls between 0.4 and 0.5 Ma, which was followed by occasional andesitic activity, the most important being the ample andesitic lava flows (∼4.1 km³) that descended the N and NW sides of the edifice. Following a ∼400 ka long repose without silicic activity, Cotopaxi began a new eruptive phase about 13 ka ago that consisted of seven rhyolitic episodes belonging to the Holocene F and Colorado Canyon series; the onset of each episode occurred at intervals of 300–3,600 years and each produced ash flows and regional tephra falls with DRE volumes of 0.2–3.6 km³. Andesitic tephras and lavas are interbedded in the rhyolite sequence. The Colorado Canyon episode (4,500 years BP) also witnessed dome and sector collapses on Cotopaxi’s NE flank which, with associated ash flows, generated one of the largest cohesive debris flows on record, the Chillos Valley lahar. A thin pumice lapilli fall represents the final rhyolitic outburst which occurred at 2,100 years BP. The pumices of these Holocene rhyolitic eruptions are chemically similar to those of older rhyolites of the Barrancas series, with the exception of the initial eruptive products of the Colorado Canyon series whose chemistry is similar to that of the 211 ka ignimbrite of neighboring Chalupas volcano. Since the Colorado Canyon episode, andesitic magmatism has dominated Cotopaxi’s last 4,400 years, characterized by scoria bomb and lithic-rich pyroclastic flows, infrequent lava flows that reached the base of the cone, andesitic lapilli and ash falls that were carried chiefly to the W, and large debris flows. Andesitic magma emission rates are estimated at 1.65 km³ (DRE)/ka for the period from 4,200 to 2,100 years BP and 1.85 km³ (DRE)/ka for the past 2,100 years, resulting in the present large stratocone.     

Regimes of magma recharge and their control on the eruptive behaviour during the period 2001–2005 at Mt. Etna volcano

Mount Etna volcano (Italy) during the period 2001–2005 has undergone a period of intense eruptive activity marked by three large eruptions (2001, 2002–2003 and 2004–2005). These eruptions encompassed diverse eruptive styles and regimes: from intensely explosive, during 2001 and 2002–2003 eruptions, to exclusively effusive in the 2004–2005 event. In this work, we put forward the idea that these three eruptions are the response of the progressive arrival into the uppermost segment of the open-conduit system of a new magma, which was geochemically distinct in terms of trace element and Sr–Nd–Pb isotope signature from the products previously emitted by the Etnean volcano. The magma migrated upwards mainly through a peripheral tectonic system, which can be considered as eccentric in spite of its relative proximity to the main system. The ingress of the new magma and its gradual displacement from the eccentric system into the uppermost sector of the open-conduit gave rise to different eruptive behaviours. At the beginning, the ascent of the undegassed magma, able to exsolve a gas phase at depth, and its interaction with closed-system magma reservoirs less than 10 km deep gave rise to the explosive events of 2001 and 2002–2003. Later, when the same magma entered into the open-conduit system, it took part in the steady-state degassing and partially lost its volatile load, leading to a totally effusive eruption during the 2004–2005 event. One further consideration highlighted here is that in 2001–2005, migration of the feeding axis from an eccentric and peripheral position towards the main open-conduit has led to the development of a new vent (South East Crater 2) located at the eastern base of the South East Crater through which most of the subsequent Etnean activity occurred.     

Combined electromagnetic, geochemical and thermal surveys of Taal volcano (Philippines) during the period 2005–2006

Since 1572, 33 phreatic to phreatomagmatic eruptions have occurred on Taal volcano (Philippines), some of them causing several hundred casualties. Considering the time delay between two consecutive eruptions, there is an 88% probability that Taal volcano should have already erupted. Since 1992, several phases of seismic activity have been recorded accompanied by ground deformation, opening of fissures, and surface activity. The volcanic activity of Taal appears to be controlled by dike injections and magma supply, buffered by a hydrothermal system that releases fluids and heat through boiling and subsequent steaming. In early 2005, a multidisciplinary project was launched for studying the hydrothermal activity. To map the hydrothermal system, combined surveys were carried out to investigate self-potential, total magnetic field, ground temperature and carbon dioxide soil degassing, along with satellite thermal imaging of the Main Crater Lake. The elevated temperatures and high concentrations of carbon dioxide, as well as electromagnetic anomalies, indicate large-scale hydrothermal degassing. This process is enhanced along the tectonic features (e.g., crater rim and faults) of the volcano, while active fissures opened along the E–W northern flank during the 1992–1994 seismic activity. Heat and fluids from the hydrothermal system are essentially released in the northern part of the crater, which is bounded to the South by a suspected NW–SE fault along which seismicity seems to take place, and dikes are thought to be intruded. During the January 2005 surveys, a new seismic crisis started, and the felt earthquakes prompted spontaneous evacuation of hundreds of inhabitants living on the volcano. Repeated surveys show changes of self-potential, total magnetic field, and ground temperature with time, without any noticeable spatial enlargement. These observations suggest that the northern flank located between the crater rim and the 1992–1994 fissures is connected with a deep thermal source in Main crater and is reactivated during seismic crises. This sector could be subjected to flank failure.     

K–Ar geochronology of a Quaternary monogenetic volcano group in Ojika Jima District, Southwest Japan

 An unspiked K–Ar dating method using a mass-fractionation correction procedure was applied to a Quaternary independent group of monogenetic volcanoes, Ojikajima Volcano Group, located in northwestern Kyushu in Southwest Japan, in order to clarify in detail secular variations in eruptive volume, locations of eruptive vents, and magma compositions in a single monogenetic volcano group. The major results were as follows: (a) K–Ar ages of monogenetic volcanoes distribute from 1.08 to 0.30 Ma, with voluminous peaks at approximately 1.0 and 0.6 Ma. (b) The volcanic activity commenced in the central part of the field, expanded to the whole field at approximately 0.6 Ma, and then shrank to the central area. (3) Concentrations of incompatible elements, such as Ba, K, and Nb, increase with decreasing age, whereas P, Y, and Zr concentrations remain constant. These concentrations suggest gradual decrease in the degree of partial melting from an identical mantle source with residual garnet. Received: 15 December 1997 / Accepted: 23 May 1998     

Dike emplacement near Parícutin volcano,Mexico in 2006

A major seismic swarm occurred near Parícutin volcano between the end of May and early July 2006. More than 700 earthquakes with magnitude (M _L ) exceeding 2.4 were located. Parícutin, located in the Michoacán–Guanajuato volcanic field in western Mexico, is well known as the site of the 1943 eruption in which a new 400 m cinder cone was constructed in what had been farmland. The 2006 swarm exhibits all of the characteristics typically associated with swarms of volcanic origins. The earthquake rate showed the typical ramp up and ramp down over the course of several days. Magnitudes were evenly distributed in time with a notably high b-value of 2.45. The earthquake locations cluster around a northeast-striking trend extending approximately 6 km. Over the first two weeks, hypocenters migrated steadily a few hundred meters per day, rising from 9 to 5 km depth and moving northeast about 5 km. On approximately June 7, the ascent of hypocenters stalled. For the next three weeks, hypocenters held their depth while migrating laterally back to the southwest. Focal mechanisms during the first part of the swarm reflected the increased stress caused by dike inflation. Following June 7, the stress orientation changed and became more consistent with the inflation of horizontal sill-like structures. Though only limited information is available from the seismic swarm preceding the 1943 eruption, several features, including the swarm duration and magnitude relationships, were comparable to those of the 2006 episode. The strong indicators of a magmatic origin to the 2006 swarm suggest that at this location there are few, if any, traditional seismic discriminants that could be used to distinguish which seismic swarms and dike emplacement events might culminate in eruption.     

Giant gas bubbles in a rheomorphic vent fill at the Las Cañadas caldera,Tenerife (Canary Islands)

During rheomorphism subsequent to fallout deposition, a portion of the densely welded fallout of the La Grieta Member flowed back into the vent from where it was erupted, while the rest of it flowed down the outer slopes of the Las Cañadas caldera in Tenerife. The welded fallout and conduit-vent structure are physically connected and constitute a rare example of this type of deposits rooted to its feeder conduit and exposed in the caldera wall. The lower part of the vent-filling rheomorphic rocks shows gas bubbles and cavities that increase in size (up to 4 m) down vent. Bubbles are deformed against other bubbles, against the steep vent walls, flattened parallel to the flow foliation planes, and elongated parallel to the flow lineation and flow fold axes. The preservation of such giant bubbles, rather than their formation, seems to be a pretty unique feature of the phonolitic products investigated here and it is likely the result of the combination of factors that acted to preserve, in the surrounding of the glass transition interval, the sealing and the late stage cooling of a pressurized system. In addition, strain drop at the base of the vent-filling rheomorphic flow caused by flow stopping against vertical vent walls may have promoted rapid gas exsolution and the formation of large bubbles.     

Seismicity of Ruapehu volcano,New Zealand, 1971–1996: a review

From 1971 until 1995, the style of seismicity at Ruapehu changed little, reflecting a period of relatively low eruptive activity and consequent long-term stability within the vent system. Volcanic earthquakes and volcanic tremor were both dominated by a frequency of about 2 Hz. Volcanic earthquakes accompanied all phreatic and phreatomagmatic eruptions, but not small hydrothermal eruptions that originated within Crater Lake. Furthermore, more than half of the M_L>3 volcanic earthquakes and changes in the reduced displacement of 2 Hz volcanic tremor by as much as a factor of 20 occurred without any accompanying eruptive activity. Three and 7 Hz volcanic tremor were also recorded, although never at lower-elevation seismometers. At times, this tremor was stronger at the summit seismometer than the 2 Hz tremor. Their source regions were independent of the 2 Hz source, and located at shallower depths. Volcano-tectonic earthquakes were generally unrelated to eruptive activity. The seismicity accompanying the 1995–1996 eruptive activity was significantly different from that of the period 1971 to 1995, and included volcanic tremor with a frequency of less than 1 Hz, simultaneous changes in the amplitude of the previously independent 2 Hz and 7 Hz volcanic tremor, and finally a change in the frequency content of volcanic earthquakes and volcanic tremor from 2 Hz to wideband. Path transmission effects play an important role in determining the characteristics of seismograms at Ruapehu. The presence of Crater Lake affects both the style of eruptions and the accompanying seismicity.     

The behavior of seasonal variations in induced seismicity in the Koyna–Warna region,western India

Based on the earthquake catalog data for the Koyna–Warna region of induced seismicity in western India, the seasonal variations in seismic activity associated with annual fluctuations in the reservoir water level are analyzed over the time span of the entire history of seismological observations in this region. The regularities in the time changes in the structure of seasonal variations are revealed. The seasonal seismic activity is minimal in May–June when the reservoir level is lowest. During the remaining part of the year, the activity has three peaks: the fall peak in September, winter peak in November–December, and spring peak in February–March. The first mentioned peak, which falls in the phase of the water level reaching its maximal seasonal value is considered as the immediate response of the fluid saturated medium to the additional loading under the weight of reservoir water. The two subsequent maxima concur with the decline phase in the reservoir level and are interpreted as the delayed response associated with the changes in the properties of the medium due to water diffusion. It is shown that the intensities of the immediate and delayed responses to the seasonal water level variations both vary with time as does their ratio. The probable factors affecting the variations in the intensity of the seasonal components of the reservoir-induced seismicity are discussed.     

Aqueous sulfur speciation possibly linked to sublimnic volcanic gas-water interaction during a quiescent period at Yugama crater lake,Kusatsu–Shirane volcano,Central Japan

After the phreatic eruption in 1982–83, volcanic activities at Kusatsu–Shirane volcano had been decreasing and reached a minimum in 1990, had turned to a temporal rise in activity up to 1994 and then decreased again at least up to 1997. During this low-activity period we observed a relatively short (≤ 1 y) cyclic variation in polythionates (PT) in the Yugama lake water. Spectral power density analysis of the PT time-series by an autoregressive (maximum entropy spectral estimation, MESE) method indicates that the major frequency in the PT variations is 1.0 y^− 1 and the minor is 2–3 y^− 1 (1.0 and 0.3–0.5 y in periodicity, respectively). Annual variations in the lake temperature are ruled out for explaining these periodicities. We attribute these cyclic variations to a cyclic magnification-reduction in meteoric-water injection into a hydrothermal regime where volcanic gases from cooling magma bodies at depth and cooler oxidized groundwater come into contact with each other. This interaction may result in a periodical change in the composition and flux of SO₂ and H₂S gases being discharged into the lake and forming PT. From a phase deviation (2–3 months) in the cycles between the annual precipitation, including snowmelt, and the PT time-series, we estimated the maximal depth of a hydrothermal reservoir beneath the lake assuming a vertical hydraulic conductivity (5 × 10^− 3 cm/s) of the volcanic detritus around the summit hydrothermal system. Chloride in the lake water reached a maximum 1.5 years faster than PT. This is most likely due to a gradual elevation of the potentiometric front of a concentrated sublimnic solution in the hydrothermal reservoir. Variations of dissolved SO₂ and H₂S in the lake water were not consistent with those of the fumarolic gases on the north flank of the volcano. This fact together with additional observations strongly suggests that these fumaroles may have the same origin but are chemically modified by a subsurface aquifer. The PT monitoring at active crater lakes during a quiescent period can provide insight into the annual expansions and reductions of a volcano-hosted hydrothermal reservoir.