首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Event tree structures constitute one of the most useful and necessary tools in modern volcanology for assessment of hazards from future volcanic scenarios (those that culminate in an eruptive event as well as those that do not). They are particularly relevant for evaluation of long- and short-term probabilities of occurrence of possible volcanic scenarios and their potential impacts on urbanized areas. In this paper, we introduce Hazard Assessment Event Tree (HASSET), a probability tool, built on an event tree structure that uses Bayesian inference to estimate the probability of occurrence of a future volcanic scenario and to evaluate the most relevant sources of uncertainty from the corresponding volcanic system. HASSET includes hazard assessment of noneruptive and nonmagmatic volcanic scenarios, that is, episodes of unrest that do not evolve into volcanic eruption but have an associated volcanic hazard (e.g., sector collapse and phreatic explosion), as well as unrest episodes triggered by external triggers rather than the magmatic system alone. Additionally, HASSET introduces the Delta method to assess precision of the probability estimates, by reporting a 1 standard deviation variability interval around the expected value for each scenario. HASSET is presented as a free software package in the form of a plug-in for the open source geographic information system Quantum Gis (QGIS), providing a graphically supported computation of the event tree structure in an interactive and user-friendly way. We also include further in-depth explanations for each node together with an application of HASSET to Teide-Pico Viejo volcanic complex (Spain).  相似文献   

2.
Long-term multi-hazard assessment for El Misti volcano (Peru)   总被引:1,自引:1,他引:0  
We propose a long-term probabilistic multi-hazard assessment for El Misti Volcano, a composite cone located <20 km from Arequipa. The second largest Peruvian city is a rapidly expanding economic centre and is classified by UNESCO as World Heritage. We apply the Bayesian Event Tree code for Volcanic Hazard (BET_VH) to produce probabilistic hazard maps for the predominant volcanic phenomena that may affect c.900,000 people living around the volcano. The methodology accounts for the natural variability displayed by volcanoes in their eruptive behaviour, such as different types/sizes of eruptions and possible vent locations. For this purpose, we treat probabilistically several model runs for some of the main hazardous phenomena (lahars, pyroclastic density currents (PDCs), tephra fall and ballistic ejecta) and data from past eruptions at El Misti (tephra fall, PDCs and lahars) and at other volcanoes (PDCs). The hazard maps, although neglecting possible interactions among phenomena or cascade effects, have been produced with a homogeneous method and refer to a common time window of 1 year. The probability maps reveal that only the north and east suburbs of Arequipa are exposed to all volcanic threats except for ballistic ejecta, which are limited to the uninhabited but touristic summit cone. The probability for pyroclastic density currents reaching recently expanding urban areas and the city along ravines is around 0.05 %/year, similar to the probability obtained for roof-critical tephra loading during the rainy season. Lahars represent by far the most probable threat (around 10 %/year) because at least four radial drainage channels can convey them approximately 20 km away from the volcano across the entire city area in heavy rain episodes, even without eruption. The Río Chili Valley represents the major concern to city safety owing to the probable cascading effect of combined threats: PDCs and rockslides, dammed lake break-outs and subsequent lahars or floods. Although this study does not intend to replace the current El Misti hazard map, the quantitative results of this probabilistic multi-hazard assessment can be incorporated into a multi-risk analysis, to support decision makers in any future improvement of the current hazard evaluation, such as further land-use planning and possible emergency management.  相似文献   

3.
The Earth's surface is characterized by the presence of many active volcanoes, most of which are surrounded by ancient villages. High-valued historical sites are often so exposed that it becomes imperative to perform volcanic risk assessment including cultural heritage. For the safeguard of the historical property in volcanic areas, two major problems are definition of (a) criteria for diagnosis and evaluation of hazard and vulnerability, and (b) methods for risk prevention and mitigation. In this paper, we first review the state-of-the-art and most outstanding geophysical prospecting and modeling methods currently on the use, which contribute to the solution of the problems mentioned above. We then show the results of an application on the most alarming volcano in Italy, Mount Vesuvius in the Neapolitan area. The imaged configuration of the feeding and plumbing systems induces to consider Vesuvius a high-risk volcano with a high probability of pyroclastic flow in case of reactivation. Finally, we show the results from a modeling approach of a pyroclastic flow simulating the eruptive scenario of Vesuvius compatible with its internal structure and dynamics. The simulation shows that the emplacement of artificial barriers close to the eruptive vent is a practical solution to reduce the local radial momentum of the pyroclastic flow and to transfer the related energy to the vertical buoyant cloud. The Vesuvius case history allows us to conclude that the integrated geophysical surveying and modeling approach can notably contribute to make decisions and also for the protection of the historical heritage in active volcanic areas.  相似文献   

4.
We propose a long-term volcanic hazards event tree for Teide-Pico Viejo stratovolcanoes, two complex alkaline composite volcanoes that have erupted 1.8–3 km3 of mafic and felsic magmas from different vent sites during the last 35 ka. This is the maximum period that can be investigated from surface geology and also represents an upper time limit for the appearance of the first phonolites on that volcano. The whole process of the event tree construction was divided into three stages. The first stage included the determination of the spatial probability of vent opening for basaltic and phonolitic eruptions, based on the available geological and geophysical data. The second, involved the analysis of the different eruption types that have characterised the volcanic activity from Teide during this period. The third stage focussed on the generation of the event tree from the information obtained in the two previous steps and from the application of a probabilistic analysis on the occurrence of each possible eruption type. As for other volcanoes, the structure of the Teide-Pico Viejo Event Tree was subdivided into several steps of eruptive progression from general to more specific events. The precursory phase was assumed as an unrest episode of any geologic origin (magmatic, hydrothermal or tectonic), which could be responsible for a clear increase of volcanic activity revealed by geophysical and geochemical monitoring. According to the present characteristics of Teide-Pico Viejo and their past history, we started by considering whether the unrest episode would lead to a sector collapse or not. If the sector collapse does not occur but an eruption is expected, this could be either from the central vents or from any of the volcanoes' flanks. In any of these cases, there are several possibilities according to what has been observed in the period considered in our study. In the case that a sector collapse occurs and is followed by an eruption we considered it as a flank eruption. We conducted an experts elicitation judgement to assign probabilities to the different possibilities indicated in the event tree. We assumed long term estimations based on existing geological and historical data for the last 35 Ka, which gave us a minimum estimate as the geological record for such a long period is incomplete. However, to estimate probabilities for a short term forecast, for example during an unrest episode, we would need to include in the event tree additional information from the monitoring networks, as any possible precursors that may be identified could tell us in which direction the system will evolve. Therefore, we propose to develop future versions of the event tree to include also the precursors that might be expected on each path during the initial stages of a new eruptive event.  相似文献   

5.
Citlaltépetl or Pico de Orizaba is the highest active volcano in the North American continent. Although Citlaltépetl is at present in repose, its eruptive history reveals repetitive explosive eruptions in the past. Its relatively low eruption rate has favored significant population growth in areas that may be affected by a potential eruptive activity. The need of some criteria for hazards assessment and land-use planning has motivated the use of statistical methods to estimate the time and space distribution of volcanic hazards around this volcano. The analysis of past activity, from late Pleistocene to historic times, and the extent of some well-identified deposits are used to calculate the recurrence probabilities of eruptions of various size during time periods useful for land-use planning.  相似文献   

6.
The estimation of site intensity occurrence probabilities in low seismic activity regions has been studied from different points of view. However, no method has been definitively established because several problems arise when macroseismic historical data are incomplete and the active zones are not well determined. The purpose of this paper is to present a method that estimates site occurrence probabilities and at the same time measures the uncertainties inherent in these probabilities in low activity regions. The region to be studied is divided into very broad seismic zones. An exponential intensity probability law is adjusted for each zone and the degree of uncertainty in the assumed incompleteness of the catalogue is evaluated for each intensity. These probabilities are used to establish what may be termed ‘prior site occurrence models’. A Bayesian method is used to improve ‘prior models’ and to obtain the ‘posterior site occurrence models’. Epicentre locations are used to recover spatial information lost in the prior broad zoning. This Bayesian correction permits the use of specific attenuation for different events and may take into account, by means of conservative criteria, epicentre location errors. Following Bayesian methods, probabilities are assumed to be random variables and their distribution may be used to estimate the degree of uncertainty arising from (a) the statistical variance of estimators, (b) catalogue incompleteness and (c) mismatch of data to prior assumptions such as Poisson distribution for events and exponential distribution for intensities. The results are maps of probability and uncertainty for each intensity. These maps exhibit better spatial definition than those obtained by means of simple, broad zones. Some results for Catalonia (NE of Iberian Peninsula) are shown.  相似文献   

7.
The expert system approach enables geochemical evidence to be integrated with geological, petrological and mineralogical evidence in identifying the eruptive setting of ancient volcanic rocks. This paper explains the development of ESCORT, an Expert System for Characterization of Rock Types. ESCORT uses as its knowledge base a set of dispersion matrices derived from a geochemical data bank of some 8000 immobile element analyses, together with tables of magma-type membership probabilities based on a-priori knowledge, the dominant flow type, the facies of intercalated sediment, the basement to the volcanic sequence, the presence or absence of an iron-enrichment trend and the phenocryst phase assemblage. Inferencing is based on the Bayes Decision Rule, adapted to take into account possible uncertainties in the geological evidence. Upper and lower probability thresholds are used to decide whether an interpretation is likely or unlikely. Evaluation of the system shows that it is able to overcome many of the well-publicized ambiquities attached to geochemical discrimination diagrams, as well as providing a unified method of decision making and a measure of the correctness of the interpretation.  相似文献   

8.
Among the schemes for earthquake forecasting, the search for semi-periodicity during large earthquakes in a given seismogenic region plays an important role. When considering earthquake forecasts based on semi-periodic sequence identification, the Bayesian formalism is a useful tool for: (1) assessing how well a given earthquake satisfies a previously made forecast; (2) re-evaluating the semi-periodic sequence probability; and (3) testing other prior estimations of the sequence probability. A comparison of Bayesian estimates with updated estimates of semi-periodic sequences that incorporate new data not used in the original estimates shows extremely good agreement, indicating that: (1) the probability that a semi-periodic sequence is not due to chance is an appropriate estimate for the prior sequence probability estimate; and (2) the Bayesian formalism does a very good job of estimating corrected semi-periodicity probabilities, using slightly less data than that used for updated estimates. The Bayesian approach is exemplified explicitly by its application to the Parkfield semi-periodic forecast, and results are given for its application to other forecasts in Japan and Venezuela.  相似文献   

9.
Ash produced by a volcanic eruption on Iceland can be hazardous for both the transatlantic flight paths and European airports and airspace. In order to begin to quantify the risk to aircraft, this study explored the probability of ash from a short explosive eruption of Hekla Volcano (63.98°N, 19.7°W) reaching European airspace. Transport, dispersion and deposition of the ash cloud from a three hour ‘explosive’ eruption with an initial plume height of 12 km was simulated using the Met Office's Numerical Atmospheric-dispersion Modelling Environment, NAME, the model used operationally by the London Volcanic Ash Advisory Centre. Eruptions were simulated over a six year period, from 2003 until 2008, and ash clouds were tracked for four days following each eruption.Results showed that a rapid spread of volcanic ash is possible, with all countries in Europe facing the possibility of an airborne ash concentration exceeding International Civil Aviation Organization (ICAO) limits within 24 h of an eruption. An additional high impact, low probability event which could occur is the southward spread of the ash cloud which would block transatlantic flights approaching and leaving Europe. Probabilities of significant concentrations of ash are highest to the east of Iceland, with probabilities exceeding 20% in most countries north of 50°N. Deposition probabilities were highest at Scottish and Scandinavian airports. There is some seasonal variability in the probabilities; ash is more likely to reach southern Europe in winter when the mean winds across the continent are northerly. Ash concentrations usually remain higher for longer during summer when the mean wind speeds are lower.  相似文献   

10.
Probabilistic characterizations of possible future eruptive scenarios at Vesuvius volcano are elaborated and organized within a risk-based framework. In the EXPLORIS project, a wide variety of topics relating to this basic problem have been pursued: updates of historical data, reinterpretation of previous geological field data and the collection of new fieldwork results, the development of novel numerical modelling codes and of risk assessment techniques have all been completed. To achieve coherence, many diverse strands of evidence had to be unified within a formalised structure, and linked together by expert knowledge. For this purpose, a Vesuvius ‘Event Tree’ (ET) was created to summarise in a numerical-graphical form, at different levels of detail, all the relative likelihoods relating to the genesis and style of eruption, development and nature of volcanic hazards, and the probabilities of occurrence of different volcanic risks in the next eruption crisis. The Event Tree formulation provides a logical pathway connecting generic probabilistic hazard assessment to quantitative risk evaluation. In order to achieve a complete parameterization for this all-inclusive approach, exhaustive hazard and risk models were needed, quantified with comprehensive uncertainty distributions for all factors involved, rather than simple ‘best-estimate’ or nominal values. Thus, a structured expert elicitation procedure was implemented to complement more traditional data analysis and interpretative approaches. The structure of the Vesuvius Event Tree is presented, and some of the data analysis findings and elicitation outcomes that have provided initial indicative probability distributions to be associated with each of its branches are summarized. The Event Tree extends from initiating volcanic eruption events and hazards right through to human impact and infrastructure consequences, with the complete tree and its parameterisation forming a quantitative synoptic framework for comprehensive hazard evaluation and mapping of risk impacts. The organization of the Event Tree allows easy updating, as and when new information becomes available.  相似文献   

11.
Analysis of the patterns of eruption occurrences may improve our understanding of volcanic processes. In this paper, the available historical data of an individual volcano, Colima in México, are used to classify its eruptions by size using the Volcanic Explosivity Index (VEI). The data shows that, if eruptions are only taken into account above a certain VEI level, the stochastic process associated with the explosive volcanic events can be represented by a non-stationary Poisson point process, which can be reduced to a homogeneous Poisson process through a transformation of the time axis. When eruptions are separated by VEI values, the occurrence patterns of each magnitude category can also be represented by a Poisson distribution. Analysis of the rate of occurrence of all eruptions with VEI greater than 1 permits the recognition of three distinct regimes or rates of volcanic activity during the last 430 years. A double stochastic Poisson model is suggested to describe this non-stationary eruptive pattern of Colima volcano and a Bayesian approach permits an estimation the present hazard.  相似文献   

12.
Bayes estimate of the probability of exceedance of annual floods   总被引:1,自引:1,他引:1  
In this paper Lindley's Bayesian approximation procedure is used to obtain the Bayes estimate of the probability of exceedence of a flood discharge. The Bayes estimates of the probability of exceedence has been shown by S.K. Sinha to be equivalent to the estimate of the probability of exceedence from the predictive or Bayesian disribution, of a future flood discharge. The evaluation of complex ratios of multiple integrals common in a Bayesian analysis is not necessary using Lindley's procedure. The Bayes estimates are compared to those obtained by the method of maximum likelihood and the method of moments. The results show that Bayes estimates of the probability of exceedence are larger as expected, but have smaller posterior standard deviations.  相似文献   

13.
The island of Tenerife is volcanically complex, and its eruptive history predominantly reflects the processes and products of two different eruptive styles: (1) non-explosive effusions of basaltic lavas from fissure vents mostly aligned along two ridges; and (2) less frequent but explosive salic eruptions from central vents associated with the Las Cañadas volcanic edifice and associated summit caldera. We have taken into account this fundamental distinction to develop a volcanic-hazards zonation (for lava flows and ash fall only) that includes: definition of the principal hazards; identification of the areas that have higher probability of containing emission centres; and numerical modelling of the vulnerable areas to be affected by volcanic hazards. Not only does the volcanic-hazards zonation map provide emergency-management officials with an updated assessment of the volcanic hazards, but it also represents a starting point for the preparation of a volcanic risk map for Tenerife. Finally, the hazards-zonation map also furnishes the basis for the design of a proposed volcano surveillance network.  相似文献   

14.
A long-standing problem in operational seismology is that of reliable focal depth estimation. Standard analyst practice is to pick and identify a ‘phase’ in the P-coda. This picking will always produce a depth estimate but without any validation it cannot be trusted. In this article we ‘hunt’ for standard depth phases like pP, sP and/or PmP but unlike the analyst we use Bayes statistics for classifying the probability that polarization characteristics of pickings belong to one of the mentioned depth phases given preliminary epicenter information. In this regard we describe a general-purpose PC implementation of the Bayesian methodology that can deal with complex nonlinear models in a flexible way. The models are represented by a data-flow diagram that may be manipulated by the analyst through a graphical-programming environment. An analytic signal representation is used with the imaginary part being the Hilbert transform of the signal itself. The pickings are in terms of a plot of posterior probabilities as a function of time for pP, Sp or PmP being within the presumed azimuth and incident angle sectors for given preliminary epicenter locations. We have tested this novel focal depth estimation procedure on explosion and earthquake recordings from Cossack Ranger II stations in Karelia, NW Russia, and with encouraging results. For example, pickings deviating more than 5° off ‘true’ azimuth are rejected while Pn-incident angle estimate exhibit considerable scatter. A comprehensive test of our approach is not quite easy as recordings from so-called Ground Truth events are elusive.  相似文献   

15.
Compositional Bayesian indicator estimation   总被引:1,自引:1,他引:0  
Indicator kriging is widely used for mapping spatial binary variables and for estimating the global and local spatial distributions of variables in geosciences. For continuous random variables, indicator kriging gives an estimate of the cumulative distribution function, for a given threshold, which is then the estimate of a probability. Like any other kriging procedure, indicator kriging provides an estimation variance that, although not often used in applications, should be taken into account as it assesses the uncertainty of the estimate. An alternative approach to indicator estimation is proposed in this paper. In this alternative approach the complete probability density function of the indicator estimate is evaluated. The procedure is described in a Bayesian framework, using a multivariate Gaussian likelihood and an a priori distribution which are both combined according to Bayes theorem in order to obtain a posterior distribution for the indicator estimate. From this posterior distribution, point estimates, interval estimates and uncertainty measures can be obtained. Among the point estimates, the median of the posterior distribution is the maximum entropy estimate because there is a fifty-fifty chance of the unknown value of the estimate being larger or smaller than the median; that is, there is maximum uncertainty in the choice between two alternatives. Thus in some sense, the latter is an indicator estimator, alternative to the kriging estimator, that includes its own uncertainty. On the other hand, the mode of the posterior distribution estimator, assuming a uniform prior, is coincidental with the simple kriging estimator. Additionally, because the indicator estimate can be considered as a two-part composition which domain of definition is the simplex, the method is extended to compositional Bayesian indicator estimation. Bayesian indicator estimation and compositional Bayesian indicator estimation are illustrated with an environmental case study in which the probability of the content of a geochemical element in soil being over a particular threshold is of interest. The computer codes and its user guides are public domain and freely available.  相似文献   

16.
17.
A singularly complete record exists for the eruptive activity of Etna volcano. The time series of occurrence of flank eruptions in the period 1600–1980, in which the record is presumably complete, is found to follow a stationary Poisson process. A revision of the available data shows that eruption durations are rather well correlated with the estimates of the volume of lava flows. This implies that the magnitude of an eruption can be defined directly by its duration. Extreme value statistics are then applied to the time series, using duration as a dependent variable. The probability of occurrence of a very long (300 days) eruption is greater than 50% only in time intervals of the order of 50 years. The correlation found between duration and total output also allows estimation of the probability of occurrence of a major event which exceeds a given duration and total flow of lava. The composite probabilities do not differ considerably from the pure ones. Paralleling a well established application to seismic events, extreme value theory can be profitably used in volcanic risk estimates, provided that appropriate account is also taken of all other variables.  相似文献   

18.
The Auckland Volcanic Field (AVF) is a young basaltic field that lies beneath the urban area of Auckland, New Zealand’s largest city. Over the past 250,000 years the AVF has produced at least 49 basaltic centers; the last eruption was only 600 years ago. In recognition of the high risk associated with a possible future eruption in Auckland, the New Zealand government ran Exercise Ruaumoko in March 2008, a test of New Zealand’s nation-wide preparedness for responding to a major disaster resulting from a volcanic eruption in Auckland City. The exercise scenario was developed in secret, and covered the period of precursory activity up until the eruption. During Exercise Ruaumoko we adapted a recently developed statistical code for eruption forecasting, namely BET_EF (Bayesian Event Tree for Eruption Forecasting), to independently track the unrest evolution and to forecast the most likely onset time, location and style of the initial phase of the simulated eruption. The code was set up before the start of the exercise by entering reliable information on the past history of the AVF as well as the monitoring signals expected in the event of magmatic unrest and an impending eruption. The average probabilities calculated by BET_EF during Exercise Ruaumoko corresponded well to the probabilities subjectively (and independently) estimated by the advising scientists (differences of few percentage units), and provided a sound forecast of the timing (before the event, the eruption probability reached 90%) and location of the eruption. This application of BET_EF to a volcanic field that has experienced no historical activity and for which otherwise limited prior information is available shows its versatility and potential usefulness as a tool to aid decision-making for a wide range of volcano types. Our near real-time application of BET_EF during Exercise Ruaumoko highlighted its potential to clarify and possibly optimize decision-making procedures in a future AVF eruption crisis, and as a rational starting point for discussions in a scientific advisory group. It also stimulated valuable scientific discussion around how a future AVF eruption might progress, and highlighted areas of future volcanological research that would reduce epistemic uncertainties through the development of better input models.  相似文献   

19.
A compilation of paleomagnetic data from volcanic deposits of Mount St. Helens is presented in this report. The database is used to determine signature paleomagnetic directions of products from its Holocene eruptive events, to assign sampled units to their proper eruptive period, and to begin the assembly of a much larger database of paleomagnetic directions from Holocene volcanic rocks in western North America. The paleomagnetic results from Mount St. Helens are mostly of high quality, and generally agree with the division of its volcanic deposits into eruptive episodes based on previous geologic mapping and radiocarbon dates. The Muddy River andesite's paleomagnetic direction, however, indicates that it is more likely part of the Pine Creek eruptive period rather than the Castle Creek period. In addition, the Two-Fingers andesite flow is more likely part of the Middle Kalama eruptive period and not part of the Goat Rocks period. The paleomagnetic data from Mount St. Helens and Mount Hood document variation in the geomagnetic field's pole position over the last ~2,500 years. A distinct feature of the new paleosecular variation (PSV) record, similar to the Fish Lake record (Oregon), indicates a sudden change from rapid clockwise movement of the pole about the Earth's spin axis to relatively slow counterclockwise movement at ~800 to 900 years B.P.  相似文献   

20.
利用遥感影像的光谱分析和纹理分析技术,研究了长白山天池火山地区IKONOS影像的光谱特征以及纹理特征。对照野外实际调查结果,选取了8种训练样本类别,首先使用波谱角技术(SAM)和最大似然法进行分类,并通过混淆矩阵评价分类结果;然后,引入了ENVI的二阶概率矩阵法的纹理分析,利用差异性(Dissimilarity)纹理特征图对地物建立了纹理解译标志;最后,通过光谱和纹理结合的方法得出了最终的分类结果,使得分类精度得到了很大提高。以研究区广泛分布的浮岩为例,确定了浮岩的分布范围和相对厚度分布,分析了不同厚度浮岩的地形特征  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号