Carbon dioxide diffuse emission from the soil: ten years of observations at Vesuvio and Campi Flegrei (Pozzuoli), and linkages with volcanic activity

Carbon dioxide flux from the soil is regularly monitored in selected areas of Vesuvio and Solfatara (Campi Flegrei, Pozzuoli) with the twofold aim of i) monitoring spatial and temporal variations of the degassing process and ii) investigating if the surface phenomena could provide information about the processes occurring at depth. At present, the surveyed areas include 15 fixed points around the rim of Vesuvio and 71 fixed points in the floor of Solfatara crater. Soil CO₂ flux has been measured since 1998, at least once a month, in both areas. In addition, two automatic permanent stations, located at Vesuvio and Solfatara, measure the CO₂ flux and some environmental parameters that can potentially influence the CO₂ diffuse degassing. Series acquired by continuous stations are characterized by an annual periodicity that is related to the typical periodicities of some meteorological parameters. Conversely, series of CO₂ flux data arising from periodic measurements over the arrays of Vesuvio and Solfatara are less dependent on external factors such as meteorological parameters, local soil properties (porosity, hydraulic conductivity) and topographic effects (high or low ground). Therefore we argue that the long-term trend of this signal contains the “best” possible representation of the endogenous signal related to the upflow of deep hydrothermal fluids.     

The determination of deep temperatures by means of the CO-CO2-H2-H2O geothermometer: an example using fumaroles in the Campi Flegrei,Italy

Chromatographic analyses of fumarolic gases, collected in sampling bottles containing an alkaline solution, have been carried out using a thermal conductivity detector and a flame ionization detector, after catalytic conversion of CO and CH₄. The latter method enables the concentration of carbon monoxide to be measured with sufficient accuracy for use in a CO-CO₂-H₂-H₂O geothermometer. Application of this geothermometer to fumaroles in the crater of Solfatara in the Campi Flegrei, Italy, indicates that they are fed from a steam reservoir at 250±15 °C and at 10^–36±2atm of oxygen. On the other hand, the CH₄-CO₂-H₂-H₂O geothermobarometer seems to re-equilibrate at superficial temperatures and cannot be used for infering thermodynamic conditions at depth. Regular sampling of these fumaroles together with a geothermometric interpretation of the gas analyses provides a means of monitoring, with comparative accuracy, the chemical and thermal evolution of the hydrothermal reservoir below the Solfatara crater. Such monitoring would probably detect an increase in temperature at depth and the injection of magmatic gas into the reservoir.     

CO2 and H2S concentrations in the atmosphere at the Solfatara of Pozzuoli

The CO₂ and H₂S concentration in the Solfatara atmosphere has been measured. The concentrations of both gases are higher neraby the more active areas and decrease away from them. A sharp horizontal and vertical gradient of the CO₂ content has been recognized.Such gradient is assumed to result from a diffusion of gas from the ground to the atmosphere.The total output of CO₂ has been computed based on a turbulent diffusion model. The obtained value is in good agreement with previously abserved values (Italiano et al., 1984).The feasibility of monitoring the atmosphere of Solfatara for either gas hazard and surveillance of volcanic activity has also been evaluated.     

Monitoring volcanic hazard using eddy covariance at Solfatara volcano, Naples, Italy

An eddy covariance (EC) station was deployed at Solfatara crater, Italy, June 8–25, 2001 to assess if EC could reliably monitor CO₂ fluxes continuously at this site. Deployment at six different locations within the crater allowed areas of focused gas venting to be variably included in the measured flux. Turbulent (EC) fluxes calculated in 30-min averages varied between 950 and 4460 g CO₂ m⁻² d⁻¹; the highest measurements were made downwind of degassing pools. Comparing turbulent fluxes with chamber measurements of surface fluxes using footprint models in diffuse degassing regions yielded an average difference of 0% (±4%), indicating that EC measurements are representative of surface fluxes at this volcanic site. Similar comparisons made downwind of degassing pools yielded emission rates from 12 to 27 t CO₂ d⁻¹ for these features. Reliable EC measurements (i.e. measurements with sufficient and stationary turbulence) were obtained primarily during daytime hours (08:00 and 20:00 local time) when the wind speed exceeded 2 m s⁻¹. Daily average EC fluxes varied by ±50% and variations were likely correlated to changes in atmospheric pressure. Variations in CO₂ emissions due to volcanic processes at depth would have to be on the same order of magnitude as the measured diurnal variability in order to be useful in predicting volcanic hazard. First-order models of magma emplacement suggest that emissions could exceed this rate for reasonable assumptions of magma movement. EC therefore provides a useful method of monitoring volcanic hazard at Solfatara. Further, EC can monitor significantly larger areas than can be monitored by previous methods.     

The gas/steam ratio as indicator of heat transfer at the Solfatara fumaroles,Phlegraean Fields (Italy)

A simple geochemical model of Solfatara, Phlegraean Fields (Italy), is proposed on the basis of gas composition and temperature at the surface.Data on the Solfatara fumaroles have been collected since 1979 within the framework of a geochemical monitoring for the surveillance of the Phlegraean volcanic system.Surface manifestations of Solfatara are likely to be fed through isoenthalpic expansion of dry steam, which separates from a geothermal liquid in an intensively fractured zone at about 236°C. This value is consistent both with gas composition and surface temperature.The gas/steam ratio appears to be the most effective parameter to detect changes of heat flow at depth.Actually a remarkable decrease in the gas/steam ratio has been observed since 1981, while the gas composition and the temperature did not change significantly. These facts suggest increased heat flow at depth.     

Geothermal energy release at the Solfatara of Pozzuoli (Phlegraean Fields): Phreatic and phreatomagmatic explosion risk implications

The H₂O, CO₂ and H₂S outputs at the Solfatara of Pozzuoli have been measured and a map of the exhaling areas has also been made. The energy released at the surface by the fluids has been estimated to be 10¹⁹ ergs/day.The presence of aquifers at Phlegraean Fields increases the phreatic and phreatomagmatic explosion risk.Our results suggest that even if an uprising magma may interact with water at depth, an explosion could occur only at the shallow levels of a few hundred meters. Since the transfer of energy toward the surface is favoured by the presence of fractures, a detailed analysis of the deep fracture network would help to evaluate the risk levels of the various areas of Phlegraean Fields.     

Composition of hydrothermal fluids during the bradyseismic crisis which commenced at Phlegraean Fields in 1982

A systematic geochemical surveillance on the fumaroles of Solfatara and the boiling pools of Pisciarelli was carried out by discontinuous monitoring of the chemical composition of the emitted fluids during the Phlegraean Fields bradyseismic crisis which has begun in 1982. The fluids are considered to be produced by the ebullition of shallow aquifers receiving a convective gaseous inflow from the underlying magma chamber.Increased water vapor concentrations at a constant temperature of about 155 °C throughout the investigated period, along with the occurrence of ground deformations and seismic phenomena, are interpreted as resulting from an increased heat supply to the boiling water bodies.Dissolution processes and reactions with the confining rocks can alter the chemical composition of fluids escaping from magma to a large extent. Therefore it does not appear correct to consider the absolute values of any chemical constituent for geochemical surveillance without taking this modifying factor into account. Acid gases will be preferably absorbed by the above mentioned aquifers, while other species like H₂, N₂, O₂, CH₄, will instead increase their relative concentrations. Because of this, water vapour concentrations and the ratios H₂S/CO₂ and H₂/CH₄ in surface thermal manifestations appear to reflect better the varying extent of the observed phenomenon.On the basis of these parameters, and of both the upheaval rate and the intensity of seismic events, maximum values in the convective input of magmatic origin are estimated to have occurred at the beginning of the crisis and in September–October 1983.As long as water bodies at shallow depth are able to buffer the convective flow both thermally and chemically, no important volcanic activity can develop. When the absorbing capacity of these aquifers is exhausted, the increasing temperature and the changing characteristics of fluids towards a magmatic composition will indicate a higher probability of eruptive phenomena.CNR —Centro di studio per la mineralogia e la geochimica dei sedimenti.     

Continuous monitoring of distal gas emanations at Vulcano,southern Italy

The increasing activity of Vulcano Island (Italy) since 1985 led to the initiation of continuous geochemical monitoring of the lateral soil gas emissions. On the basis both of their relative geochemical characteristics and of local considerations, three gaseous components were selected for monitoring, namely CO₂, He and ²²²Rn. Monitoring has been performed by means of specific analysers. Gases extracted from a water well located at the foot of the active cone were selected for monitoring, on the basis of their geochemical and isotopic characters that indicate their genetic link with central high temperature fumarolic gases emitted at the crater. Very strong variations of gas composition can be observed within one day (from 1 to about 94% for CO₂). Some variations display a daily character and can be correlated with that of atmospheric pressure. The three monitored gases are highly correlated, suggesting very high kinetics of gas transfer in the system. Because of these considerable variations of chemical composition, bulk concentrations obviously are not suitable for monitoring at Vulcano. However, the evolution with time of ratios such as ²²²Rn/CO₂ and He/CO₂ (the latter being corrected for atmospheric contamination) supplies numerical parameters that the expected to characterize the intensity of the degassing process. A new input of magmatic gases, that would lead to an increase in the ²²²Rn/CO₂ and He/CO₂ ratios, should therefore be detected by such a monitoring station.     

The 1975–1977 crisis of la Soufriere de Guadeloupe (F.W.I): A still-born magmatic eruption

On July 8, 1976, eruptive activity broke out at la Soufrière de Guadeloupe (F.W.I) after about one year of increasing seismic activity. Seismic activity continued to increase until August 1976, reaching more than 1500 events (a 200-fold increase over the preceding quiet period of a few years) and an energy output of about 10¹⁷ ergs in a day. A total of 26 major phreatic eruptions similar to the July 8 outburst took place during an eight-months period. The steam blasts that characterized the eruptions gave rise to particle- and sometimes block-charged plumes that deposited an estimated 10⁶ m³ of solids. The H₂O-rich gases emitted during the blasts presumably contained other gases (H₂S, SO₂, CO₂...) that were partly adsorbed on solid particles. All material was erupted at temperatures of the order of 100° to 200°C.The observation of vertical migration of earthquake foci in less than a few hours and over about 6 km depth, and of abnormal variations of the geomagnetic field, indicate a deep energy source for the phreatic eruptions. A small proportion of the gases adsorbed on solid particles had a magmatic origin. However, most of the steam and the tephra seemed to originate from superficial levels of a hydrothermal system. Similar phreatic eruptions have occurred several times in recorded history. In the case of la Soufrière, the origin of the phreatic eruptions is best described by an abnormal energy input (versus steady-state) from a crustal magma chamber. The occurrence of truly magmatic eruptions is presumably inhibited by an extensive hydrothermal system. The abrupt release of more power from the magma chamber could have resulted in an explosive pyroclastic eruption.Substantial improvement of the Guadeloupe volcano observatory has followed the 1975–1977 crisis. Permanent telemetered geophysical networks and regular geochemical observations have provided a five year data base of the volcano behavior in its noneruptive state which can be compared to crisis situations.     

Evolution of Br/Cl ratios in fumarolic salammoniac from Vulcano (Aeolian Islands,Italy)

 Bromine/chlorine ratios in salammoniac samples, collected from 1983 to 1994 at the Fossa crater fumaroles of Vulcano (Aeolian Islands), were determined. The data obtained in this work have been compared to those obtained by Coradossi et al. (1985) in order to carry out a systematic study of the Br/Cl ratio variations since 1977. Br/Cl ratios range from 0.04 to 0.29. The Br/Cl ratio variations observed have been interpreted as the result of different contributions among deep and shallow components feeding the fumarolic fluids of Vulcano. The Br/Cl increases, which occurred from 1977 to 1979, in 1983–1984 and in 1992–1993, have been connected to an increased input of magmatic gases. Received: 10 November 1995 / Accepted: 12 July 1996     

Modeling time dependence in a tidal boundary layer

Although the rate of acceleration and deceleration in a depth-limited tidal flow can be quite large, it has been found that quasi-steady logarithmic models or steady eddy diffusivitics can explain the turbulent field quite well (Gross andNowell, 1983;Long, 1981). A simple, time-dependent model was devised based on kinetic energy closure (Bradshaw et al., 1967) to simulate a depth-limited tidal boundary layer without the assumption of a time-independent eddy diffusivity. Computer solutions revealed that the parameterS = u_s/δΩ (u_s is a turbulent velocity scale, δ is depth, Ω is tidal frequency) governs the phasing of the velocity to the forcing pressure gradient. However, for values ofS reasonable for depth-limited estuarine tidal flow the ratio of the turbulent energy production-dissipation balance to rate of change of turbulent energy is always quite large, indicating that the turbulent adjustment rate is fast enough for the flow to be modeled with quasi-steady techniques.     

Measuring Concentrations of Dissolved Methane and Ethane and the 13C of Methane in Shale and Till

Baseline characterization of concentrations and isotopic values of dissolved natural gases is needed to identify contamination caused by the leakage of fugitive gases from oil and gas activities. Methods to collect and analyze baseline concentration‐depth profiles of dissolved CH₄ and C₂H₆ and δ¹³C‐CH₄ in shales and Quaternary clayey tills were assessed at two sites in the Williston Basin, Canada. Core and cuttings samples were stored in Isojars^® in a low O₂ headspace prior to analysis. Measurements and multiphase diffusion modeling show that the gas concentrations in core samples yield well‐defined and reproducible depth profiles after 31‐d equilibration. No measurable oxidative loss or production during core sample storage was observed. Concentrations from cuttings and mud gas logging (including IsoTubes^®) were much lower than from cores, but correlated well. Simulations suggest the lower concentrations from cuttings can be attributed to drilling time, and therefore their use to define gas concentration profiles may have inherent limitations. Calculations based on mud gas logging show the method can provide estimates of core concentrations if operational parameters for the mud gas capture cylinder are quantified. The δ¹³C‐CH₄ measured from mud gas, IsoTubes^®, cuttings, and core samples are consistent, exhibiting slight variations that should not alter the implications of the results in identifying the sources of the gases. This study shows core and mud gas techniques and, to a lesser extent, cuttings, can generate high‐resolution depth profiles of dissolved hydrocarbon gas concentrations and their isotopes.     

Restoration of new volcanic gas analyses from basalts of the afar region: Further evidence of CO2 -degassing trends

The restored compositions for approximately 70 new analyses reported recently for Erta'Alelava lake (LeGuern et al., 1979) are in good agreement with restored compositions (Gerlach, 1980a) based on previously published data. The results confirm earlier indications that gas collections taken at different times from the lava lake are related principally by variations in CO₂ content. Restored compositions for gas samples collected in the final stages of a November 1978 Ardoukoba eruption along the Asal Rift spreading axis resemble the Erta'Ale gases except for a much lower CO₂ content. The Ardoukoba gases fall close to a CO₂-decreasing control line for gases with initial compositions similar to the 1971–1973 Erta'Ale gases. These results suggest that gases released from basaltic lava along zones of crustal spreading follow compositional trends dominated by changes in CO₂ content.     

Modeling hydrothermal fluid circulation and gravity signals at the Phlegraean Fields caldera

The Phlegraean Fields caldera is an active volcanic system where episodes of ground deformation are accompanied by significant changes in geochemical and geophysical parameters monitored at the surface. These changes derive from a complex interaction between magmatic system and hydrothermal fluid circulation. We calculate the gravity changes associated with the variable density of hydrothermal fluids. We simulate the multi-phase and multi-component fluid circulation triggered by a pulsating magma degassing, periodically increasing the discharge of CO₂-enriched fluids into the shallow hydrothermal system. The simulated evolution of the hydrothermal system successfully reproduces the observed composition of gas discharged at the surface. At the same time, results indicate that changes in average fluid density generate a detectable gravity signal that is of the same order of magnitude of the observed changes. This contribution to gravity changes can explain the peculiar behavior of gravity data collected at Solfatara, where surface hydrothermal phenomena are present. Simultaneous fitting of two independent sets of monitoring data (gas composition and gravity changes) confirms the conceptual model proposed for the hydrothermal system at Solfatara, and it provides new insights for the interpretation of gravity data.     

Spatial and temporal variations of nutrient concentration in the groundwater of a floodplain: effect of hydrology,vegetation and substrate

Spatio‐temporal variations in nitrogen and phosphorus concentrations in groundwater were analysed and related to the variations in hydrological conditions, vegetation type and substrate in an alluvial ecosystem. This study was conducted in the Illwald forest in the Rhine Plain (eastern France) to assess the removal of nutrients from groundwater in a regularly flooded area. We compared both forest and meadow ecosystems on clayey‐silty soils with an anoxic horizon (pseudogley) at 1·5–2 m depth (eutric gley soil) and a forest ecosystem on a clayey‐silty fluviosoil rich in organic matter with a gley at 0·5 m depth (calcaric gley soil). Piezometers were used to measure the nutrient concentrations in the groundwater at 2 m depth in the root layer and at 4·5 m depth, below the root layer. Lower concentrations of nitrate and phosphate in groundwater were observed under forest than under meadow, which could be explained by more efficient plant uptake by woody species than herbaceous plants. Thus NO₃‐N inputs by river floods were reduced by 73% in the shallow groundwater of the forested ecosystem, and only by 37% in the meadow. Compared with the superficial groundwater layer, the lowest level of nitrate nitrogen (NO₃‐N) and the highest level of ammonium nitrogen (NH₄‐N) were measured in the deep layer (under the gley horizon at 2·5 m depth), which suggests that the reducing potential of the anoxic horizon in the gley soils contributes to the reduction of nitrate. Nitrate concentrations were higher in the groundwater of the parcel rich in organic matter than in the one poorer in organic matter. Phosphate (PO₄‐P) concentrations in both shallow and deep groundwater are less than 62 to 76% of those found in surface water which can be related to the retention capacity of the clay colloids of these soils. Moreover, the temporal variations in nutrient concentrations in groundwater are directly related to variations in groundwater level during an annual hydrological cycle. Our results suggest that variations in groundwater level regulate spatio‐temporal variations in nutrient concentrations in groundwater as a result of the oxidation–reduction status of soil, which creates favourable or unfavourable conditions for nutrient bioavailability. The hydrological variations are much more important than those concerning substrate and type of vegetation. Copyright © 1999 John Wiley & Sons, Ltd.     

Geochemistry of the Submarine Gaseous Emissions of Panarea (Aeolian Islands, Southern Italy): Magmatic vs. Hydrothermal Origin and Implications for Volcanic Surveillance

The marine sector surrounding Panarea Island (Aeolian Islands, South Italy) is affected by widespread submarine emissions of CO₂ -rich gases and thermal water discharges which have been known since the Roman Age. On November 3^rd, 2002 an anomalous degassing event affected the area, probably in response to a submarine explosion. The concentrations of minor reactive gases (CO, CH₄ and H₂) of samples collected in November and December, 2002 show drastic compositional changes when compared to previous samples collected from the same area in the 1980s. In particular the samples collected after the November 3^rd phenomenon display relative increases in H₂ and CO and a strong decrease in the CH₄ contents, while other gas species show no significant change. The interaction of the original gas with seawater explains the variable contents of CO₂, H₂S, N₂, Ar and He which characterize the different samples, but cannot explain the large variations of CO, CH₄ and H₂ which are instead compatible with changes in the redox, temperature and pressure conditions of the system. Two models, both implying an increasing input of magmatic fluids are compatible with the observed variations of minor reactive species. In the first one, the input of magmatic fluids drives the hydrothermal system towards atypical (more oxidizing) redox conditions, slowly pressurizing the system up to a critical state. In the second one, the hydrothermal system is flashed by the rising high-T volcanic fluid, suddenly released by a magmatic body at depth. The two models have different implications for volcanic surveillance and risk assessment: In the first case, the November 3^rd event may represent both the culmination of a relatively slow process which caused the overpressurization of the hydrothermal system and the beginning of a new phase of quiescence. The possible evolution of the second model is unforeseeable because it is mainly related to the thermal, baric and compositional state of the deep magmatic system that is poorly known.     

Tsunami damage to coastal defences and buildings in the March 11th 2011 M w 9.0 Great East Japan earthquake and tsunami

On March 11th 2011 a M _w 9.0 mega-thrust interface subduction earthquake, the Great East Japan Earthquake, occurred 130 km off the northeast coast of Japan in the Pacific Ocean at the Japan Trench, triggering tsunami which caused damage along 600 km of coastline. Observations of damage to buildings (including vertical evacuation facilities) and coastal defences in Tōhoku are presented following investigation by the Earthquake Engineering Field Investigation Team (EEFIT) at 10 locations in Iwate and Miyagi Prefectures. Observations are presented in the context of the coastal setting and tsunami characteristics experienced at each location. Damage surveys were carried out in Kamaishi City and Kesennuma City using a damage scale for reinforced concrete (RC), timber and steel frame buildings adapted from an earlier EEFIT tsunami damage scale. Observations show that many sea walls and breakwaters were overtopped, overturned, or broken up, but provided some degree of protection. We show the extreme variability of damage in a local area due to inundation depth, flow direction, velocity variations and sheltering. Survival of many RC shear wall structures shows their high potential to withstand local earthquake and significant tsunami inundation but further research is required into mitigation of scour, liquefaction, debris impact, and the prevention of overturning failure. Damage to steel and timber buildings are also discussed. These observations are intended to contribute to mitigation of future earthquake and tsunami damage by highlighting the key features which influence damage level and local variability of damage sustained by urban coastal infrastructure when subjected to extreme tsunami inundation depths.     

The determination of deep temperatures by means of the CO-CO2-H2-H2O geothermometer: an example using fumaroles in the Campi Flegrei, Italy. A comment

Three points raised in the paper by Tedesco and Sabroux (1987) are dealt with. (1) The inconsistency between the water partial pressure calculated by Tedesco and Sabroux (1987) and saturation pressure is due to the improper use of the water-gas-shift reaction as a geothermometer. In fact Tedesco and Sabroux (1987) do not take into account the distribution of gas species between the coexisting vapour and liquid phases. (2) The depth of the steam reservoir is evaluated by Tedesco and Sabroux (1987) in too simplistic a way. This matter should be treated with greater care owing to the high social impact of any consideration on the Phlegraean Fields system. (3) The reliability of carbon monoxide determination at the concentration level encountered at Solfatara depends on the collection method rather than on the gas-chromatographic technique.     

The determination of deep temperatures by means of the CO-CO2-H2-H2O geothermometer: an example using fumaroles in the Campi Flegrei,Italy. A comment

Three points raised in the paper by Tedesco and Sabroux (1987) are dealt with. (1) The inconsistency between the water partial pressure calculated by Tedesco and Sabroux (1987) and saturation pressure is due to the improper use of the water-gas-shift reaction as a geothermometer. In fact Tedesco and Sabroux (1987) do not take into account the distribution of gas species between the coexisting vapour and liquid phases. (2) The depth of the “steam reservoir” is evaluated by Tedesco and Sabroux (1987) in too simplistic a way. This matter should be treated with greater care owing to the high social impact of any consideration on the Phlegraean Fields system. (3) The reliability of carbon monoxide determination at the concentration level encountered at Solfatara depends on the collection method rather than on the gas-chromatographic technique.     

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.