versión On-line ISSN 0718-7106
AndGeo vol.39 no.3 Santiago set. 2012
Andean Geology 39 (3): 380-393. September, 2012
formerly Revista Geológica de Chile
On the likelihood of future eruptions in the Chilean Southern Volcanic Zone: interpreting the past century's eruption record based on statistical analyses
Probabilidades de futuras erupciones en la Zona Volcánica del Sur de Chile: interpretación estadística de la serie temporal de erupciones del siglo pasado
Yvonne Dzierma1, Heidi Wehrmann2
1 SFB 574, Christian-Albrechts-Universitat zu Kiel, Otto-Hahn-Platz 1, D-24118 Kiel, Germany. firstname.lastname@example.org
2 SFB 574, GEOMAR Helmholtz Centre for Ocean Research Kiel, Wischhofstr. 1-3, 24148 Kiel, Germany email@example.com
ABSTRACT. A sequence of 150 explosive eruptions recorded during the past century at the Chilean Southern Volcanic Zone (SVZ) is subjected to statistical time series analysis. The exponential, Weibull, and log-logistic distribution functions are fit to the eruption record, separately for literature-assigned volcanic exploslvlty indices (VEI) ≥ 2 and VEI ≥ 3. Since statistical tests confirm the adequacy of all the fits to describe the data, all models are used to estimate the likelihood of future eruptions. Only small differences are observed between the different distribution functions with regard to the eruption forecast, whereby the log-logistic distribution predicts the lowest probabilities. There is a 50% probability for VEI ≥ 2 eruptions to occur in the SVZ within less than a year, and 90% probability to occur within the next 2-3 years. For the larger VEI ≥ 3 eruptions, the 50% probability is reached in 3-4 years, while the 90% level is reached in 9-11 years.
Keywords: Chilean Southern Volcanic Zone, Eruption time series analysis, Probabilistic eruption forecasting, Volcanic hazards.
RESUMEN. Se presenta un análisis estadístico de la serie temporal de 150 erupciones volcánicas explosivas registradas durante el siglo pasado en la Zona Volcánica del Sur de Chile. Se modeló el conjunto de erupciones mediante la distribución exponencial, de Weibull y log-logística, restringiendo el análisis a erupciones de índice de explosividad volcánica (IEV) mayores a 2 y 3, respectivamente. Como los modelos pasan las pruebas estadísticas, los tres modelos se aplican para estimar la probabilidad de erupciones futuras. Se observan solo diferencias menores entre las predicciones mediante los distintos modelos, con la distribución log-logística dando las probabilidades más bajas. Para erupciones de IEV ≥ 2, la probabilidad de producirse una erupción dentro de un año es más del 50%, creciendo al 90% en 2-3 años. Para erupciones más grandes, de IEV ≥ 3, el 50% de probabilidad se alcanza dentro de 3-4 años, y el 90% dentro de 9-11 años.
Palabras clave: Zona Volcánica del Sur de Chile, Análisis de series de tiempo de erupciones volcánicas, Predicción estadística de erupciones, Peligros volcánicos.
The Chilean Southern Volcanic Zone (SVZ) is the volcanologically most active part of the Chilean Andes, stretching from Tupungatito (33.40°S, 69.80°W) in the wider Santiago region (Región Metropolitana) southwards to Cerro Hudson (45.90°S, 72.97°W) in central Patagonia. Several population centres used for living, agriculture, industry and commerce are in the vicinity of these volcanoes. Coastal lowlands to mountainous regions are touristically attractive and frequently visited for recreation throughout the year. Larger cities, in particular Santiago de Chile, while not in the immediate vicinity of the volcanic chain, are located inside drainage channels from the Andes and would, in the event of a volcanic eruption, be threatened by possible lahars and pyroclastic flows. The entire region is endangered by fallout tephra from explosive eruptions, which may cover large areas at long distances from the eruptive source, also reaching Argentina to the east. Estimating the risk posed by this volcanic arc segment is thus of paramount importance for the most densely populated region of Chile and Argentina.
Continuous monitoring, i.e. the surveillance of seismic activity, gas release and composition, ground deformation, scientific constraints on the regional tectonics and geochemistry are the most powerful tools for assessing the current state of a volcano. Recent efforts by the Observatorio Volcanológico de los Andes del Sur' (OVDAS, part of SERNAGEOMIN, the Servicio Nacional de Geología y Minería), in Chile have greatly improved the monitoring of these volcanoes by the installation of a large number of seismometers on the Chilean volcanoes considered as most active (Fernando Gil, personal communication, 2010). However, in-detail monitoring requires substantial funds for staff and instrumentation and can therefore only be performed for selected volcanoes. Monitoring and scientific investigation in the Andes are additionally complicated by the still poor accessibility of some volcanoes at high altitudes in rugged terrains. Moreover, the unexpected 2008 violent eruption of Chaitén volcano should be taken as an alert that not only those volcanoes with a recent (historical) eruptive record should be recognised in hazard considerations.
Since all these volcanoes are located in the same volcanic arc, formed by subduction of the Nazca Plate beneath the South American Plate, they share a collective susceptibility to the tectonic driving forces, which may potentially affect the system on a large scale. Although individual volcanoes may not be matter of statistical analysis because of too few eruptions, a compiled investigation will pay credit to their importance for relaxing the system by volumetric magma release and stress changes, and may reveal changes in the eruption rates of the entire arc. Also in the context of an interdisciplinary hazard and risk estimation, it is fundamental to evaluate the overall eruption hazard in addition to the individual analyses of single eruption sites. Furthermore, few cases of registered eruptions were revealed where an eruption was initially assigned to a wrong source volcano - such a problem becomes inconsequential when regarding compiled data sets.
An additional advantage of combining all SVZ eruptions is the possibility to regard larger eruptions (VEI ≥ 3), which are usually too infrequent to be assessed for each volcano individually. Although this collective approach does not allow us to specify where the next eruption is most likely to occur, the exact site becomes less important by the potential of a larger eruption to be more devastating for a much wider area.
This paper therefore aims at describing the typical temporal eruption behaviour of the SVZ as a coherent subduction system by means of several statistical distributions. These models are then used to infer an estimate for the probability of future eruptions. With these results we provide a probabilistic part of volcanic hazard assessment, to be combined with results of the various continuous monitoring techniques.
2. Tectonic setting and volcanic activity at the SVZ
Volcanism in the SVZ is driven by the subduc-tion of the Nazca Plate beneath the South American Plate, which takes place in an oblique direction of approximately 80-82°ENE at a rate of 70-90 mm/ year. The incoming plate dips uniformly at about 25° to a depth of 90-100 km (Stern, 2004). At the northern boundary, the dip angle at depth shallows to less than 5°, considered as the reason for the lack of volcanism in the North-adjacent Pampean flat slab segment. The southern end of the SVZ is defined by the subducting Chile Ridge.
In the active arc over a length of 1,400 km, several calderas and more than 60 large volcanic complexes tower over numerous smaller cones (Fig. 1; López-Escobar et al, 1995; Stern et al, 2007). Of those, twenty have erupted not only in the Pleistocene/ Holocene, but also produced a series of eruptions in historical times.
|FIG. 1. Overview SRTM image of the Chilean Southern Volcanic Zone (SVZ) displaying the location of the volcanoes considered in this study (red triangles), other major volcanic centres within the SVZ (yellow triangles), the plate boundary (dashed blue line), major cities (white circles). Inset shows the location of the SVZ on the South American continent. Image courtesy SRTM Team NASA/JPL/NIMA, modified.|
The historical records investigated in this study are taken from the Global Volcanism Program webpage (www.volcano.si.edu, Siebert and Simkin, 2002-) and literature compilations, in particular from Petit-Breuilh (2004). Only eruptions from a minimum volcanic explosivity index (VEI, defined by Newall and Self, 1982; Simkin and Siebert, 1994) upwards are considered in this study. The eruption chronologies from the different sources are largely in agreement with each other. In the cases where discrepancies exist on eruption occurrence and VEI-assignment, the eruptions were included in the analysis if at least one source states them to be large enough (VEI ≥ 2), and no other source explicitly discredited them. The eruption records were converted into time series by calculating the repose times between successive eruptions. Following the convention established by Klein (1982), the repose time is defined as the interval from the onset of one eruption to the onset of the next one, thereby neglecting the duration of the eruptions. The repose times are resolved on a yearly scale. This results from the precision of the available source data, which sometimes do not provide exact dates, and is sufficient for most purposes. For those years in which several eruptions took place, repose times of 'zero years' are accumulated and as such included in the statistical analysis. In several cases, volcanoes showed eruptive cycles with variably strong activity and variable interruptions therein over two or several years. For example, an eruption might commence with weak activity in one year, and then culminate to a degree of explosivity sufficient for being considered in this study in the following year/s. The decision on whether to treat such eruptive activity as one longer-lasting complex eruption, or as a sequence of individual events, is not within the scope of this study. We therefore adopted the eruption demarcation as given in the available literature. As such, they form the data base for the statistical processing, and are provided with the literature-assigned VEI in table 1.
|TABLE 1. ERUPTIONS OF ALL THE SVZ VOLCANOES COMBINED, TAKEN FROM THE SMITHSONIAN'S GLOBAL VOLCANISM PROGRAM WEBSITE AND FROM THE COMPILATION OF PETIT-BREUILH (2004). ONLY ERUPTIONS WITH A KNOWN VEI ≥ 2 ARE LISTED. IN SOME CASES, THE VEI-ASSIGNMENT IS NOT IDENTICAL IN THE DIFFERENT SOURCES, THIS IS COMMENTED IN THE TABLE. THOSE ERUPTIONS FOR WHICH AT LEAST ONE SOURCE CLAIMS THEM TO BE VEI ≥ 2 ARE CONSIDERED IN THE STATISTICAL ANALYSIS.|
Due to the limited historical record, our study will only include eruptions from the year 1900 onwards. Although the statistical method would be more reliable if based on a longer eruption record, in particular since short-term fluctuations in activity over time scales of centuries might bias the analysis, the older part of the record becomes progressively more unreliable. For those volcanoes analysed individually by Dzierma and Wehrmann (2010a), the year 1900 presents an acceptable threshold after which the eruption record can be reasonably assumed to be complete based on statistical tests.
The eruption filtering by a minimum magnitude is performed here because very small eruptions are more prone not to be registered in the eruption record. For eruptions of higher explosivity index (VEI ≥ 4), the eruption record is too scarce for statistical analysis, although their devastating energy is obviously much larger. Our measure here is eruption
occurrence; the events are not weighted by their magnitude. The analysis is twofold performed here, (1) including all eruptions with a VEI ≥ 2, and (2) for all eruptions with VEI ≥ 3. The collective eruption record is visualised as a cumulative frequency plot over time in figure 2.
|FIG. 2. Cumulative numbre of eruptions as a function of time for all volcanoes of the SVZ combined; top image VEI ≥ 2, low image: VEI ≥ 3). The green line corresponds to an overall constant eruption rate; wich presents a relatively good approximation to the data.|
The method employed here is a simple application of the standard statistical lifetime and failure analysis, presented elsewhere in detail (Dzierma and Wehrmann, 2010a, b). Several studies have used this technique or parts of it to characterise the eruption behaviour and quantify the eruption hazard of volcanoes around the world (e.g., Wickman, 1966; Klein, 1982; Ho, 1990; Bebbington and Lai, 1996a; De la Cruz-Reyna, 1996; Scarpa and Tilling, 1996; Connor et al., 2003; Mendoza-Rosas and De la Cruz-Reyna, 2008, 2009, 2010).
In brief, to determine whether the eruptions can be modelled as simple stochastic point processes, it is assessed to what extent successive eruptions are independent of each other by calculating the correlation between successive repose times. Whether or not the time series can be considered stationary is checked using a 5-point moving average of the repose times and assessing the deviations from the average.
Given the assumption of a stationary stochastic point process, several standard models are available which can be used to represent the repose time distribution function:
or the corresponding survival function:
The most widely-used statistical models are: a. the exponential distribution; b. the Weibull distribution and c. the log-logistic distribution function.
The exponential distribution mathematically represents the most simple scenario where the eruptions can be described by a Poisson process (for an introduction, see, e.g., Cox and Lewis, 1966; Cox and Oakes, 1984). In this case, the rate of eruptions λ is constant with time, and the survival function is given by the functional form:
As an extension of this model, the Weibull distribution takes into account a shape parameter d which facilitates to describe systems with a temporally changing hazard rate, usually attributed to one dominant process acting in the system:
Depending on whether the shape parameter d is larger or smaller than one, the hazard rate of the system is increasing or decreasing with time, respectively (Ho, 1991; Bebbington and Lai, 1996a, b; Watt et al., 2007).
As additional conceptual adjustment, the log-logistic model (Connor et al., 2003, 2006; Dirksen, 2006) accommodates the inclusion of competing processes, where some factors increase the probability of an eruption with time, whereas others counteract them, decreasing the eruption probability:
All three survival functions are fit to the data, optimised by a least-squares method, and the quality and adequateness of the fits evaluated using several statistical measures, the Kolmogorov-Smirnov-test (KS-test, Gibbons, 1976), χ2 test and corrected Akaike Information Criterion (AICc, Akaike, 1973; Sigiura, 1978; Burnham and Anderson, 1998).
The fit functions can then be used to estimate the probability of at least one VEI ≥ 2 or VEI ≥ 3 eruption within a given time span t in the future as:
where x is the time that has elapsed since the last eruption (Marshall and Olkin, 2007). In the present case, since the year 2011 accommodates the onset of the Puyehue eruption, the statistical forecast is calculated as seen from the beginning of 2012.
The first visual impression of the chronological eruptive sequence is obtained from the cumulative number of eruptions versus time plot (Fig. 2). For the VEI ≥ 3 limit, the cumulative number of eruptions increases approximately linearly with time. The plot displaying the VEI ≥ 2 eruptions shows an overall increase over time that can generally be described as narrow scattering bent over a linear increase, indicating that the series of eruption occurrence has remained relatively stable since 1900, with maybe a slight decrease in eruption rate over time. If the time series was to be described in much greater detail, a stepwise consideration could be invoked, distinguishing several short regimes of higher and lower activity levels. This is, however, beyond the scope of this study. When regarding the subduction system on a large scale and in simplified terms, the long-term driving forces involved in provoking volcanic activity have not undergone any drastic changes during the short eruption record considered here, which in part allows us to assume near-stationarity for the overall eruption rate.
Statistically testing the time series for stationarity, this first visual impression is confirmed by the outcome from the correlation and moving-average analyses (Figs. 3, 4), although sensu stricto the requirements of no simultaneously occurring events is not entirely fulfilled at the chosen time resolution, as the data set contains several years in which more than one eruption took place. The correlation coefficients of R=-0.111 and R=-0.181 for VEI ≥ 2 and VEI ≥ 3 eruptions, respectively, do not indicate a significant correlation at the 5% level of significance (p-values are 0.180 and 0.408, respectively). The stationarity checks are equally satisfied, with variations in the 5-point average repose times remaining below 2 standard deviations from the mean. While there appears to be a slight trend towards increasing repose times in the recent past, this has not yet reached the point of statistical significance. From a mathematical point of view it is therefore justified to treat both datasets as stationary for the following analyses.
|FIG. 3. Serial correlation scatter plot (top image: VEI ≥ 2, low image: VEI ≥ 3). Each repose time is plotted as a function of the preceding repose time. Since many repose times occur several times and dots would overlap, jitter is applied to the data before plotting. In this way, each data point is displaced randomly a small distance around its coordinate in the plot. This is done for better visual inspection of the plot and was not included in the correlation/regression analysis. Because of the yearly resolution, points appear as being clustered on a grid instead of showing an irregular distribution. The red line gives a linear fit of the data points, the green line would be expected for perfect correlation between successive repose times.|
|FIG. 4. Stationarity check by plotting the 5-point moving repose time averages as a function of time. The red line gives the overall mean 5-point repose time average, the green dashed line denotes two standard deviations from the mean.|
The life time distributions can be fit stably in all cases, and all fits pass the goodness-of-fit test (Table 2, Fig. 5). For the VEI ≥ 2 eruptions, the exponential and Weibull distributions give the best fit to the data according to the AICc. The Weibull model differs from the exponential distribution in that it can take into account a shape parameter, which for both data sets is slightly larger than one, indicating a weakly increasing hazard rate. For the VEI ≥ 3 eruptions, the AICc attributes the best fit quality to the log-logistic and exponential distributions. Both the VEI ≥ 2 and VEI ≥ 3 eruptions can be reasonably described by a Poisson process, which agrees well with the approximately linear shape of the cumulative eruption curve (Fig. 2). More complicated models (the Weibull distribution for VEI ≥ 2 eruptions and the log-logistic model for VEI ≥ 3 eruptions) are also in agreement with the data, but do not significantly improve the fit of the exponential distribution.
|TABLE 2. FIT PARAMETERS AND RESULTS OF THE GOODNESS-OF-FIT AND QUALITY TESTS.|
|FIG. 5. Fits of the repose time distribution, for VEI ≥ 2 eruptions (top image) and VEI ≥ 3 eruptions (low image), respectively.|
Since all distribution functions pass the goodness-of-fit test, the predictions for the future eruption probabilities are calculated for all distribution functions (Fig. 6), although it appears on the basis of the AICc that the Weibull or exponential model should be preferred for the VEI ≥ 2 eruptions and the exponential or log-logistic behaviour for the VEI ≥ 3 eruptions. The differences in the predictions for the three distribution functions are slight. All of them reach over 50% probability for at least one VEI ≥ 2 eruption in less than one year, and a probability of 90% in 2-3 years. These predictions stochastically underline the high hazard posed by the SVZ volcanoes, as also directly observed by the short recurrence intervals in the order of no more than a few years and sometimes several eruptions occurring in the same year.
|FIG. 6. Probability (in %) that at least one eruption with VEI ≥ 2 (top image) or with VEI ≥ 3 (low image) will occur in the near future, based on the three different distribution function fits.|
For the VEI ≥ 3 eruptions, the probabilities are somewhat smaller. The analysis gives at most 20% probability of at least one VEI ≥ 3 eruption within one year, and the 50% probability is reached after 3-4 years. With over 90% probability, it can be expected that at least one VEI ≥ 3 eruption will occur over the next 9-11 years, depending on the distribution function.
Based on time series analysis of 150 eruptions in the historical record since 1900, we provide the likelihood for future eruptions of VEI ≥ 3 and VEI ≥ 2 to occur in the SVZ. The models predict a statistical probability of 65-72% of a VEI ≥ 3 eruption within the next 5 years, and 96-99% of a VEI ≥ 2 eruption within this time. For obvious reasons, hazard evaluation should never build up on statistical analyses alone. Volcanic settings are known to be able to produce catastrophic eruptions at short time scales, which cannot always be reliably predicted based on the average behaviour. Statistical eruption analysis as presented here can only estimate the likelihood that eruptions will take place within a given time interval. As a main outcome with regard to the high probability for future eruptions, the estimate provided here underpins the necessity for continuous volcano surveillance in the SVZ. Monitoring of the volcano-seismic signals, gas emissions, inflation/deflation, and understanding magmatic and tectonic systems are indispensable for successful short-term disaster prevention, management and mitigation. A joint assessment of these physical parameters combined with insight from statistical analyses yields the potential to reliably indicate signs of an impending eruption at short notice.
Credit is due to Dr. H. Moreno for detailed comments on the eruption records, especially for providing additional information beyond published data and beyond the compilation of the GVP. We are grateful to Prof. Dr. Ch. Stern for his encouraging evaluation of this work, and to Dr. M. Suaréz for the editorial handling of the manuscript. Dr. M.E. Petit-Breuilh kindly facilitated access to her book on the eruptive history of hispanoamerican volcanoes. This paper is contribution number 200 of the Collaborative Research Centre (Sonderforschungsbereich) 574 'Volatiles and Fluids in Subduction Zones' at Kiel University, funded by the German Research Foundation (DFG).
Akaike, H. 1973. Information Theory and an Extension of the Maximum Likelihood Principle (Petrov, B.N.; Csaki, F.; editors). In International Symposium on Information Theory, No. 2. Akademiai Kiado: 267-281. Budapest. [ Links ]
Bebbington, M.S.; Lai, C.D. 1996a. Statistical analysis of New Zealand volcanic occurrence data. Journal Volcanology Geothermal Research 74: 101-110. [ Links ]
Bebbington, M.S.; Lai, C.D. 1996b. On nonhomogenous models for volcanic eruptions. Mathematical Geology 28 (5): 585-600. [ Links ]
Burnham, K.P.;Anderson, D.R. 1998. Model selection and inference: a practical information-theoretic approach. Springer-Verlag: 353 p. New York, USA. [ Links ]
Connor, C.B.; Sparks, R.S.J.; Mason, R.M.; Bonadonna, C. 2003. Exploring links between physical and probabilistic models ofvolcanic eruptions: The Soufrière Hills Volcano, Montserrat. Geophysical Research Letters 30 (13): 1701. doi:10.1029/2003GL017384. [ Links ]
Connor, C.B.; McBirney, A.R.; Furlan, C. 2006. What is the probability of explosive eruption at long-dormant volcano? In Statistics in Volcanology (Mader, H.M.; Coles, S.G.; Connor, C.B.; Connor, L.J.; editors). Special Publications of International Association of Volcanology and Chemistry of the Earth's Interior 1 (IAVCEI). Geological Society of London: 39-46. [ Links ]
Cox, D.R.; Oakes, D. 1984. Analysis of survival data. Chapman and Hall: 201 p. London. [ Links ]
Cox, D.R.; Lewis, P.A.W. 1966. The statistical Analysis of Series of Events, Methuen & Co., London: 285 p. [ Links ]
De la Cruz-Reyna, S. 1996. Long-term probabilistic analysis of future explosive eruptions. In Monitoring and mitigation of volcanic hazards (Scarpa, R.; Tilling, R.I.; editors). Springer-Verlag: 599-629. Berlin. [ Links ]
Dirksen, O.; Humphreys, M.C.S.; Pletchov, P.; Melnik, O.; Demyanchuk, Y.; Sparks, R.S.J.; Mahony, S. 2006. The 2001-2004 dome-forming eruption of Shiveluch volcano, Kamchatka: Observation, petrological investigation and numerical modelling. Journal Volcanology Geothermal Research 155 (3-4): 201-226. [ Links ]
Dzierma, Y.; Wehrmann, H. 2010a. Statistical eruption forecast for the Chilean Southern Volcanic Zone: typical probabilities of volcanic eruptions as baseline for possibly enhanced activity following the large 2010 Concepción earthquake. Natural Hazards and Earth System Sciences 10: 1-16. doi:10.5194/ nhess-10-1-2010. [ Links ]
Dzierma, Y.; Wehrmann, H. 2010b. Eruption time series statistically examined: Probabilities of future eruptions at Villarrica and Llaima Volcanoes, Southern Volcanic Zone, Chile. Journal Volcanology Geother-mal Research 193: 82-92. doi:10.1016/j.jvolgeores.2010.03.009. [ Links ]
Gibbons, J.P. 1976. Nonparametric method for quantitative analysis. Holt, Rinehart and Winston: 463 p. New York. [ Links ]
Hildreth, W.; Drake, R.E. 1992. Volcán Quizapu, Chilean Andes. Bulletin of Volcanology 54: 93-125. [ Links ]
Ho, C.-H. 1990. Bayesian analysis of volcanic eruptions. Journal Volcanology Geothermal Research 43: 91-98. [ Links ]
Ho, C.-H. 1991. Time trend analysis of basaltic volcanism for the Yukka Mountain site. Journal Volcanology Geothermal Research 46: 61-72. [ Links ]
Klein, F. W. 1982. Patterns of historical eruptions at Hawaiian volcanoes. Journal Volcanology Geothermal Research 12: 1-35. [ Links ]
López, L.; Moreno H. 1981. Erupción de 1979 del volcán Mirador, Andes del Sur, 40°21'S: características geoquímicas de las lavas y xenolitos graníticos. Revista Geológica de Chile 13-14: 17-33. [ Links ]
López-Escobar, L.; Cembrano, J.; Moreno, H. 1995. Geochemistry and tectonics of the Chilean Southern Andes basaltic Quaternary volcanism (37°-46°S). Revista Geológica de Chile 22 (2): 219-234. [ Links ]
Marshall, A.W.; Olkin, I. 2007. Life Distributions: Structure of Nonparametric, Semiparametric, and Parametric Families. Springer Series in Statistics. Springer, New York: 782 p. [ Links ]
Mendoza-Rosas, A.T.; De la Cruz-Reyna, S. 2008. A statistical method linking geological and historical eruption time series for volcanic hazard estimations: Application to active polygenetic volcanoes. Journal Volcanology Geothermal Research 176: 277-290. doi:10.1016/j.volgeores.2008.04.005. [ Links ]
Mendoza-Rosas, A.T.; De la Cruz-Reyna, S. 2009. A mixture of exponentials distribution for a simple and precise assessment of the volcanic hazard. Natural Hazards and Earth System Sciences 9: 425-431. [ Links ]
Mendoza-Rosas, A.T.; De la Cruz-Reyna, S. 2010. Hazard estimates for El Chichón volcano, Chiapas, México: a statistical approach for complex eruptive histories. Natural Hazards and Earth System Sciences 10: 1159-1170. [ Links ]
Moreno, H. 1980. La erupción del volcán Mirador en Abril-Mayo de 1979, Lago Ranco-Riñinahue, Andes del Sur. Universidad de Chile, Comunicaciones 28: 1-23. [ Links ]
Naranjo, J.A. 1991. Nueva erupción del volcán Hudson. Revista Geológica de Chile 18: 183-184. [ Links ]
Naranjo, J.A.; Moreno, H.; Banks, N. 1993. La erupción del volcán Hudson en 1991 (46°S), Región XI, Aisén, Chile. Servicio Nacional de Geología y Minería Boletín 44: 50 p. [ Links ]
Naranjo, J.A.; Polanco, E. 2004. The 2000 AD eruption of Copahue Volcano, Southern Andes. Revista Geológica de Chile 31 (2): 279-292. doi: 10.4067/ S0716-02082004000200007. [ Links ]
Naranjo, J.A.; Moreno, H. 2005. Geología del Volcán Llaima, Región de la Araucaria. Servicio Nacional de Geología y Minería, Carta Geológica de Chile, Serie Geología Basica 88: 33 p. escala 1:50.000. [ Links ]
Newhall, C.G.; Self S. 1982. The Volcanic Explosivity Index (VEI): an estimate of explosive magnitude for historical volcanism. Journal of Geophysical Research 87 (C2): 1231-1238. [ Links ]
Petit-Breuilh Sepúlveda, M.E. 2004. La historia eruptiva de los volcanes hispanoamericanos (siglos XVI al XX). Servicio de Publicaciones del Exmo. Cabildo Insular de Lanzarote-Casa de los volcanes: 431 p. [ Links ]
Scarpa, R.; Tilling, R.I. 1996. Monitoring and mitigation of volcanic hazards. Springer-Verlag: 841 p. Berlin. [ Links ]
Sigiura, N. 1978. Further analysis of the data by akaike's information criterion and the finite corrections. Comm. Statistics-Theory and Methods 7 (1): 13-26. doi:10.1080/03610927808827599. [ Links ]
Simkin, T.; Siebert, L. 1994: Volcanoes of the World (2a edition). Geosciences Press: 349 p. Tucson, Arizona. [ Links ]
Siebert, L.; Simkin, T. 2002. Volcanoes of the World: an Illustrated Catalog of Holocene Volcanoes and their Eruptions. Smithsonian Institution, Global Volcanism Program, Digital Information Series, GVP-3, http://www.volcano.si.edu/world/, access 2011. [ Links ]
Stern, C.R.; Moreno, H.; López-Escobar, L.; Clavero, J. E.; Lara, L.E.; Naranjo, J.A.; Parada, M.A.; Skewes, M.A. 2007. Chilean volcanoes. In The Geology of Chile (Moreno, T.; Gibbons, W.; editors). Geological Society of London: 147-178. [ Links ]
Stern, C.R. 2004. Active Andean volcanism: its geologic and tectonic setting. Revista Geológica de Chile 31 (2): 161-206. [ Links ]
Watt, S.F.L.; Mather, T.A.; Pyle, D.M. 2007. Vulcanian explosion cycles: Patterns and predictability. Geology: 35 (9): 839-842. doi:10.1130/G23562A.1. [ Links ]
Wickman, F.E. 1966. Repose period patterns of volcanoes. Arkiv fur Mineralogie Geologie 4 (7): 291-367. [ Links ]
Manuscript received: May 11, 2011; revised/accepted: March 23, 2012; available online: June 19, 2012.