Vol. 48, 01, 05ok.qxd 33 ANNALS OF GEOPHYSICS, VOL. 48, N. 1, February 2005 Key words Poas volcano – geochemistry – radon – precursor 1. Introduction The historical activity of the Poas volcano (Costa Rica) was observed and reported throughout the last century, showing periods of strombolian activity in 1910 and 1953-1955. Phreatic eruptions of geyser type, with plumes of water, gravel, mud and blocks originating from the principal crater lake were also reported (Barquero et al., 1981; Casertano et al., 1983). Since radon-in-soil fluctuations and chemical variations observed in water bodies associated with active volcanoes may give information con- cerning the activity changes (Matsuo et al., 1982; De la Cruz Reyna et al., 1985, 1989; Not- su et al., 1991; Brondi and Dall’Aglio, 1991; D’Alessandro and Parello, 1994; Dall’Aglio et al., 1994; Tedesco and Pece, 1996), a radon-in- soil survey was initiated in 1981 and a systemat- ic monitoring of the chemical composition of the hot water lake has also been implemented since 1983. In the present paper, we report the radon- in-soil behaviour together with the geochemical data obtained from the water sampling as a func- tion of the volcanic activity and of its environ- mental effects during the 1981-1994 period. 2. Poas volcanic activity The Poas volcano (10.20°N, 84.23°W) is one of the main active volcanoes of the Central Range of Costa-Rica, rising up to 2708 m a.s.l. (fig. 1a). It is located in the N.E. part of the Cen- tral Valley where the most populated cities of the country are also found. The morphology of the summit presents two calderic structures with several cones such as Von Frantzius to the north and Botos to the south, where a cold freshwater lake is found. Prominent features of the active crater include a plain of eroded pyroclastic material to the south and a 30 m high cone formed during the 1953- 1954 eruptions. The core of the cone consists of Mailing address: Dr. Michel Monnin, Hydrosciences, Université Montpellier 2 (UM2) - CNRS, 34095 Montpel- lier Cedex 05, France; e-mail: monnin@msem.univ- montp2.fr Water chemistry and soil radon survey at the Poas volcano (Costa Rica) Jorge Barquero (1), Erick Fernandez (1), Michel Monnin (2), Jean Luc Seidel (2) and Nuria Segovia (3) (1) Observatorio Vulcanológico y Sismológico de Costa Rica, Universidad Nacional (OVSICORI-UNA), Heredia, Costa Rica (2) Hydrosciences, Université Montpellier 2 (UM2) - CNRS, Montpellier, France (3) Instituto Nacional de Investigaciones Nucleares (ININ), México D.F., México Abstract Radon-in-soil monitoring at the Poas volcano (Costa Rica) has been performed together with water chemistry from the hot crater lake since 1981 and 1983 respectively. The results are discussed as a function of the eruptive evolution of the volcano over a 13 years period (1981-1994). It is shown that no definitely clear precursory radon signals have been recorded. On the contrary, ionic species concentrations are likely to be considered good pre- cursors, together with the temperature variations of the crater lake water. 34 Jorge Barquero, Erick Fernandez, Michel Monnin, Jean Luc Seidel and Nuria Segovia fractured lava of basaltic-andesitic composition; the northern portion of the cone is truncated by a pit crater which contains the acidic crater lake (Rowe et al., 1992a,b). The activity of the volcano during the in-soil radon monitoring period showed a quiescence phase since 1979 which lasted for several years. In June 1987, phreatic eruptions started and con- tinued during the remaining part of the year, char- acterising a new activity phase of the volcano. The eruptions originated from the central and western part of the crater lake, ejecting water, sed- iments and gas columns. The water temperature rose from 58°C in January 1988 to 70°C in June. Fig. 1a,b. Location map of Poas volcano; radon network. a The lake level decreased by 5 m during 1987. In April 1988, the ejecta column reached more than one kilometre height and the ejected materials af- fected particularly the southern flank. The lake level decreased by almost 40 m and the tempera- ture fluctuated between 65 and 85°C. Fumaroles were also observed at the lake bottom with tem- peratures ranging from 75 to 98°C. Gas emission increased and the acidic rain affected areas of lo- cal economic importance such as coffee, grass, and other plantations. The sporadic phreatic activ- ity continued during the year at a lower degree. From April 1989, the activity intensified and the formation of a cone in the crater took b 35 Water chemistry and soil radon survey at the Poas Volcano (Costa Rica) place with the ejected material. The crater lake disappeared at the end of the same month. On April 22th an eruption projected solid material 600 m away from the crater while on April 28th the sediments ejected from the lake bottom reached 2.5 km along the southern flank. On April 30th ash eruptions started and by May 1st the ash column reached one kilometre height. The crater lake recuperated during the rainy season. A phreatic eruption occurred at the end of May destroying the lake cone. From June 1989 to early 1990, the gas emanations were permanent with sporadic phreatic eruptions. The crater lake disappeared again in April 1990. Eruptive activity was less intense than dur- ing the previous year, although ash falls reached 11 km on the southern flank around May 1990, affecting the cities of San Miguel, San Luis, Gre- cia and Trojas. Gas venting was observed at dif- ferent locations of the crater with temperatures ranging from 90°C to 550°C. From 1990 to 1994 the activity continued at a lower rate, with the hot lake almost disappearing during the dry seasons. 3. Experimental 3.1. Radon measurements Originally a radon-in-soil network consisting of 7 monitoring stations was set up in 1981, while two new stations were added in 1986 (fig. 1b). The monitoring has been carried out using Solid State Nuclear Tracks Detectors, LR-115 type 2 from Dosirad Co., France. The technique for exposure and detectors processing has been previously described (Seidel et al., 1988). How- ever, in brief, the set-up consists of a 1m long PVC pipe inserted in the ground which received a shorter 30 cm tube with an open bottom. Its up- per bottom is closed and holds inside itself the polymeric Rn detector facing down. This shorter tube plays the role of a diffusion chamber. The remaining 70 cm upper part of the PVC tube is thermally insulated with rockwool and closed. By this means one measures actually the local radon flux accross the Earth, which results after a while in a specific concentration within the dif- fusion chamber. Hence, data are expressed in kBqm–3. The detectors were recovered essential- Table I. Average radon-in-soil concentration values (X+2σ/n) and uranium content in soil at the 9 moni- toring stations. Station Radon Uranium Rn/U (kBqm–3) (ppm) a.u. S1 7.8 ± 0.4 7 1.11 S2 1.5 ± 0.1 6 0.25 S3 0.96 ± 0.1 10 0.10 S4 1.9 ± 0.2 10 0.19 S5 1.0 ± 0.2 10 0.10 S6 1.4 ± 0.1 10 0.14 S7 1.5 ± 0.1 7 0.21 S8 1.6 ± 0.3 10 0.16 S9 3.1 ± 0.3 – – ly on a fortnightly basis. However, due to logis- tic problems, e.g., temporary inaccessibility, some of the detectors were recovered after dif- ferent periods of exposure. The detectors were processed by batches in a NaOH solution and read by means of a jumping spark counter. 3.2. Chemical and physical data Since 1983, temperature and pH measure- ments of the hot water lake have been performed monthly, together with water samples collec- tion. Water samples were analysed for sulphates, chlorides and fluorides. Chlorides were quanti- fied with a combination electrode (model 93-17, ORION) with a liquid membrane. The fluoride determination was performed with a specific electrode (model 96-07, ORION) and a Tissab buffer solution. Sulphates were determined by a modified turbidimetric technique from the American Public Health Association. The con- centrations of the specific ions were quantified using analytical standard curves and the results are reported in ppm. The pH was measured with a Beckman 39820 combination electrode. 4. Results and discussions Table I presents the average radon value and the 95% confidence limit (2σ/n), during the mon- 36 Jorge Barquero, Erick Fernandez, Michel Monnin, Jean Luc Seidel and Nuria Segovia itoring period, and the uranium content of the soil at each station. Regarding stations 2, 3, 4, 5, 6, 7, and 8, the average radon concentrations show small variations and their ratio to the local urani- um contents (expressed in arbitrary units) keep within a rather narrow range (0.10 to 0.25), indi- cating that the measured radon is either support- ed by the uranium in soil content or that the ad- ditional Rn flux is, on average, similar in intensi- ty from point to point. This fact does not mean that Rn variations are not correlated with the ac- tivity of the volcano since it is well known that even when Rn signals are due to some deep phe- nomenon they show the modulation of in-soil Rn which originates from a small volume around the measuring site and is not from deep origin. The average Rn concentration as measured on these 7 stations is lower than 2 kBqm–3. Fig. 2. Raw radon time series and rainfalls from 1981 to 1994 for the whole stations. The two vertical lines in- dicate the periods of hot lake disappearence. 37 Water chemistry and soil radon survey at the Poas Volcano (Costa Rica) On the contrary, on stations 1 and 9 the av- erage is higher than 3 kBqm–3. In addition, for station 1 located at the edge of the crater, the Rn/U ratio is substantially higher (1.11), indi- cating that an additional outgassing occurred on this site. U content on station 9, which is locat- ed in an altered zone, has not been measured. The crude data recorded on the nine stations are exhibited in fig. 2, together with monthly rainfall recorded during the same period of time. At first site it seems that a strong peak appears particularly on station 1 in correlation with the 1989 eruption. However, before being able to make any data analysis some corrections must be made. It has been observed by several authors (e.g., Fleischer, 1997) that the efficiency of the type of detectors used in this survey may vary, es- pecially with the exposure time (ageing of the de- tectors and moisture). In order to account for this behaviour, on every station, the detectors were separated in several groups according to exposure time: 8 to 14 days (the most abundant) taken as reference; 1 to 7 days, 15 to 21 days and longer than 21 days. The average radon concentrations measured on every site for the above mentioned periods of times during 13 years (which legiti- mates the procedure) were calculated for a given period of time, namely 7 days. They have then be compared and related to those measured during the 8-14 exposure times. The corresponding mul- tiplicating correction factors are indicated in table II. The data have been corrected accordingly, keeping in mind that the effect of ventilation (when changing the detector) is negligible with respect to the loss of detector efficiency. Only relative variations as a function of time have significance as possible precursory signals of an eruption. Accordingly, all the data were transformed, for the sake of comparison, under the form of centred and reduced values according to the well established formula (e.g., Bry, 1995) zz z zi= - v ) ^ h where < z> is the mean value and σz the standard deviation of a given set of zi values of the Rn concentration. The same treatment was applied to rainfall data. The resulting Rn concentration and rainfall curves are exhibited in fig. 3. The importance of the corrections carried out is clearly observed by comparing the two sets of data (figs. 2 and 3). On each set of curves, the 2σ horizontal line (which is commonly accepted as the criterion for an «anomalous» deviation) has been drawn. The two vertical lines show the time of the two lake disappearances in 1989 and 1990. For instance, if one focuses on the period of qui- escence before 1989, several Rn peaks above the 2σ line are present. Among them, one is particu- larly «strong» around the end of 85-beginning of 86, on stations 3 and 4. After performing the cor- rection, the «anomaly» on station 3 has complete- ly disappeared. Another possible correlation can be observed in December 1984 at stations 2, 3, 4, 5, 6, 7. On all stations the anomalous peaks are still present afer corrrection. Yet, they all occur at a time when the volcano was particularly quiet ac- cording to the observations of the Volcanological and Seismological Observatory of Costa Rica (Barquero Ed., 1983/1987) and our own tempera- ture and chemistry measurements (see below). Over most of the sites atmospheric temperature and pressure variations do influence the in-soil Rn concentration. Before going further in the data analysis and discussion, it must be pointed out that it is not likely in this case since the environmental temperature (19.2 ± 0.8°C) and the atmospheric pressure (779 ± 0.66 mm Hg) are quite constant during the year, with very small night-to-day vari- ations. A Fast Fourier Transform analysis of the Table II. Correction coefficient calculated for each station, in function of the detectors exposure time. Exposure S1 S2 S3 S4 S5 S6 S7 S8 S9 1-7 days 0.31 0.58 0.34 0.34 0.23 0.38 0.38 0.49 0.43 8-14 days 1 1 1 1 1 1 1 1 1 15-21 days 1.22 1.20 1.32 1.41 1.20 1.20 1.13 1.06 1.19 38 Jorge Barquero, Erick Fernandez, Michel Monnin, Jean Luc Seidel and Nuria Segovia radon time series shows a 365 days periodicity peak that corresponds to the precipitation pattern of the zone (rainy season from May to November and dry season from January to April). In general the rainy season exhibits the lowest radon-in-soil concentrations due to the water presence in the pores of the upper part of the soil. A radon increase occurs during the dry season for all the stations. The conclusions that can be drawn from this data analysis are: i) Most probably, none of the «anomalous peaks» recorded during the present survey have any relation with volcanic activity. ii) Rn data from LR-115 type SSNTD detec- tors must be corrected for detection efficiency. iii) Rn concentrations measured on a long- term basis are mostly governed by rainfall. iv) If Rn precursory signals are to be found on such volcanic edifices, they must be sought us- ing a short term Rn measuring probe such as the Fig. 3. Corrected radon time series and rainfalls from 1981 to 1994 for the whole stations. 39 Water chemistry and soil radon survey at the Poas Volcano (Costa Rica) one devised by Monnin et al. (1998). As a matter of fact such a device takes care of secondary ef- fects such as external atmospheric parameters and it also detects short-term Rn outbursts that are likely to occur on volcanoes (Garcia et al., 2000). Nevertheless, if one focuses on the period when the volcano entered an active phase until it went into a quieter period, strong anomalous peaks can be observed at stations 1, 2, 3, 4, 5, 6, 7, 8 and 9. These peaks could be governed by rain- fall but we would like to point out that all stations responded to the increase in activity of the volcano and also that the amplitude of the peaks are sub- stantially higher than «dry-season» induced peaks observed previously. It could be said that both phenomena, the drying of the soil and the volcanic activity both contributed to the excess Rn release. Figure 4 shows the evolution of the hot lake water temperature during the 1983-1993 period. From January 1983 to June 1986, the temperature was quite stable (around 49°C). From July 1986, the temperature increased, reaching a maximum (94°C) in July 1990. Seismicity also correlates with the eruptive pattern, increasing and reaching a maximum during the 1989-1991 period. The variations of SO42–, Cl– and F– concen- trations of the hot crater lake water samples are shown in fig. 5. The results obtained by other Fig. 5. Variations of SO42 –, Cl – and F – concentrations. Fig. 4. Temporal variations of crater lake tempera- ture (°C) and local seismicity. 40 Jorge Barquero, Erick Fernandez, Michel Monnin, Jean Luc Seidel and Nuria Segovia investigators (Rowe et al., 1992a,b) are also shown in the same figure. Having such a double set of data covering the same period (November 1984-September 1990) and extending sideways for our data to January 1983 and to March 1993 allows us to draw some conclusions. From the end of 1984 to the beginning of 1988, the two sets of data are remarkably similar regarding in both concentration values and the behaviour of the ions. The two surveys confirm each other. SO42–, Cl–, F– and temperature measurements exhibit a stable behaviour from January 1983 to the end of 1986, with a lag time of 6 months with respect to the temperature increase which started in March 1987. A sharp increase then oc- curred by mid-1986 for all the ionic species. This behaviour is equally observed in Rowe and co-workers’ report and in the present work. The highest concentrations were measured during the first lake disappearance in April 1989. Then if the general behaviour is similar to Rowe and co-workers’ report and ours, large discrepan- cies are observed in concentrations values. In oth- er words, our chemical analysis data are affected by a large dispersion. These differences can be at- tributed to the mere fact that at this time the dis- tribution of the ionic species in the lake was per- Fig. 6. SO42 –/Cl – ratio and F –/Cl – ratio. turbed and also because sampling was not per- formed exactly at the same place, or at the same times by the teams of investigators. What remains nevertheless and, as a common finding, is that: – Regarding sulphates: there is an increase that reaches its maximum before the April 1989 erup- tion and lake disappearance; followed by a gener- al decrease till the beginning of 1991 i.e. 6 months after the second eruption and lake disappearance. Our data continue and show another increase starting in March 1991, reaching a maximum in February 1992 to decrease to a very low value at the end of the survey at the beginning of 1993. – Regarding chlorides: there is an increase, more marked at our samples than at those of Rowe and his co-workers, but at the same time, with its maximum by mid-1988, followed in both cases by a decrease reaching its lowest val- ues right before the 1989 eruption. Then again the chloride concentration increases, all the way after the second eruption. Then, our data show a slow decline of the Cl concentration till 1993. – Regarding fluorides: we again observe val- ues that are higher than those of Rowe and his co-workers but that exhibit the same behaviour i.e. an increase reaching maximum values be- tween the two eruptions, followed by an erratic 41 Water chemistry and soil radon survey at the Poas Volcano (Costa Rica) decrease, then, as in the case of chlorides, a max- imum but finally ending in an erratic distribution up to 1993. Out of these results it could be said that the eruptive phase of the Poas volcano started in depth, when the temperature of the lake water started to increase 3 years before the actual erup- tion. The ionic activities in the lake confirm the start of the eruptive phase, the ion species being brought up to the level of the lake, as was sug- gested by Rowe et al. (1992, 1995), by the cir- culation of a hydrothermal fluid convection. However, when one looks at the chloride based element ratios, it seems that it cannot be concluded that a SO42 – peak occurred before the 1989 eruption. Figure 6 displays the SO42 –/Cl– ratio values obtained by the two teams are dis- played. The peak observed by Rowe and his co- workers does not show on our data, or at least is barely marked (one point that looks more as be- ing out of the statistics that having a genuine meaning). Regarding the F–/Cl– ratio (fig. 6), the behaviour is more in agreement in the two surveys, yet with a difference of two months at least for the time of maximum value. However, we observe a comparably high ratio long after the second eruption, in 1992 i.e., 2 years later after a time when the volcano was going into a period of quiescence. In the specific case of the Poas volcano and for this very eruptive phase, it is tempting to con- clude that the precursory signal of the volcanic activity is the increased concentration of the ion- ic species, which is enhanced by the decrease of the lake volume, as Rowe et al. pointed out. This decrease of the lake volume is due itself (Brown et al., 1989) to the evaporation from the lake sur- face caused by an increase in the power output attributed to the emplacement of a small magma intrusion inferred from microgravity measure- ments (Rymer and Brown, 1989). However, be- cause of the large dispersion of the data, one can- not really rely on the chloride based ratios to es- timate the imminence of the eruption. 5. Conclusions The radon survey reported in this paper shows that radon data should be taken with great care. After appropriate corrections, this survey neither shows nor allows us to conclude posi- tively that any precursory signal was observed at the measuring stations prior to the eruptive phase 1989-1990. It could have been expected that some Rn excess could show at the stations locat- ed within the area of acidic rain since radon and HCl have similar properties regarding their solu- bility and water/air fraction ratios. It is not the case. Rn activities are similar in all the stations including station 6 which was located out of this area and which does not show any Rn depletion with respect to the other stations that were in- cluded within the area of acidic rain. It seems that long exposure Rn measurement by means of SSNTD is not suitable for this type of survey and that automatic short term Rn measurements should be preferred. The hydrochemical survey shows that an in- crease in all the species measured is definitely related to the volcanic activity but that the chlo- ride-based ratios should be taken with great care because of the large dispersion of the data. Acknowledgements Part of this project was performed under IAEA and CNRS (France)-CONACYT (Mexi- co) scientific agreements. REFERENCES BARQUERO, J., E. MALAVASSI and A. ARAYA (1981): The Poas volcano activity during the 1978-1980, Bol. Vul- canol. (Escuela de Ciencias Geograficas, Univ. Na- cional), 11,10-12 (in Spanish). BRONDI, M. and M. DALL’AGLIO (1991): Evolution of mer- cury, arsenic, antimony, radon and helium contents in ground waters and fumaroles since 1983 through 1989 at volcano Island (Southern Italy), Volcanol. Acta, 1, 233-241. BROWN, G.C., H. RYMER, J. DOWDEN, P. KAPADIA, D. STEVEN- SON, J. BARQUERO and L.D. MORALES (1989): Energy budget analysis for Poas crater lake: implications for pre- dicting volcanic activity, Nature, 339, 370-373. BRY, X. (1995): Analyse Factorielle Simple (Economica Ed., Paris), 1-112. CASERTANO, L., A. BORGIA and C. CIGOLINI (1983): The Poas volcano, Costa Rica: chronology and activity, Ge- ofis. Int., 22 (3), 215-236 (in Spanish). D’ALESSANDRO, W. and F. PARELLO (1994): Soil degassing from the lower flanks of Mt. Etna, Environ. Geochem. Health, 16 (suppl.), 189-201. 42 Jorge Barquero, Erick Fernandez, Michel Monnin, Jean Luc Seidel and Nuria Segovia DALL’AGLIO, M., F. PAGANI and F. QUATTROCCHI (1994): Geochemistry of the groundwaters in the Etna region before and after the paroxysmal phase of the eruption of December 1991: implications for the geochemical surveillance of Mt. Etna, Acta Vulcanol., 4, 149-156. DE LA CRUZ REYNA, S., M. MENA, N. SEGOVIA, J.F. CHALOT, J-L. SEIDEL and M. MONNIN (1985): Radon emanometry in soil gases and activity in ashes from El Chichon vol- cano (Mexico), Pure Appl. Geophys., 123, 410-421. DE LA CRUZ REYNA, S., M.A. ARMIENTA, V. ZAMORA and F. JUAREZ (1989): Chemical changes in spring waters at Tacana volcano, Chiapas, Mexico: a possible precursor of the May 1986 seismic crisis and phreatic explosion, J. Volcanol. Geotherm. Res., 38, 345-353. FLEISCHER, R.L. (1997): Radon and earthquake prediction, in Radon measurements by Etched Track Detectors, (World Scientific Press, Singapore), 285-300. GARCIA, R., G. NATALE, M. MONNIN and J.L. SEIDEL (2000): Shock wave radon surface signals associat- ed with the upsurge of T-P solitons in volcanic sys- tems, J. Volcanol. Geotherm. Res., 96, 15-24. MATSUO, S., J. OSSAKA, J. HIRABAYASHI, T. OZAWA and K. KI- MISHIMA (1982): Chemical nature of volcanic gases of the Usu Volcano in Japan, Bull. Volcanol., 45 (3), 261-264. MONNIN, M.M. and J.L. SEIDEL (1998): An automatic radon probe for Earth Sciences studies, J. Appl. Geophys., 39, 209-220. NOTSU, K., H. WAKITA, G. IGARASHI and T. SATO (1991): Hydrological and geochemical changes related to the 1989 seismic and volcanic activities of the Izu Penin- sula, J. Phys. Earth, 39, 245-254. RYMER, H. and G.C. BROWN (1989): Gravity changes as a precursor to volcanic eruptions at Poas volcano, Costa Rica, Nature, 342, 902-905. ROWE, G.L, S.L BRANTLEY, M. FERNANDEZ, J.F. FERNAN- DEZ, A. BORGIA and J. BARQUERO (1992a): Fluid-vol- cano interaction in an active stratovolcano: the crater lake system of Poas volcano, Costa Rica, J. Volcanol. Geotherm. Res., 49, 23-51. ROWE, G.L., S. OHSAWA, B. TAKANO, S.L. BRANTLEY, J. F. FERNANDEZ and J. BARQUERO (1992b): Using crater lake chemistry to predict volcanic activity at Poas vol- cano, Costa Rica, Bull. Volcanol., 54, 494-503. ROWE, G.L., S.L. BRANTLEY, J.F. FERNANDEZ and A. BORGIA (1995): Fluid flow and mass transfer involving an acid crater lake at Poas Volcano, Costa Rica, J. Volcanol. Geotherm. Res., 64, 233-267. SEIDEL, J.L., M. MONNIN, A. BONNEVILLE, H. DELORME, J- F. LENAT and N. SEGOVIA (1988): Radon in soil meas- urements related to the activity of the Piton de la Four- naise volcano (Reunion), Nucl. Tracks Rad. Meas., 15, 567-570. TEDESCO, D., R. PECE and R. AVINO (1996): Radon, pH and temperature monitoring in water wells at Campi Flegrei caldera (Southern Italy), Geochem. J., 30 (2), 131-138.