scholarly journals A unified probabilistic framework for volcanic hazard and eruption forecasting

2021 ◽  
Vol 21 (11) ◽  
pp. 3509-3517
Author(s):  
Warner Marzocchi ◽  
Jacopo Selva ◽  
Thomas H. Jordan

Abstract. The main purpose of this article is to emphasize the importance of clarifying the probabilistic framework adopted for volcanic hazard and eruption forecasting. Eruption forecasting and volcanic hazard analysis seek to quantify the deep uncertainties that pervade the modeling of pre-, sin-, and post-eruptive processes. These uncertainties can be differentiated into three fundamental types: (1) the natural variability of volcanic systems, usually represented as stochastic processes with parameterized distributions (aleatory variability); (2) the uncertainty in our knowledge of how volcanic systems operate and evolve, often represented as subjective probabilities based on expert opinion (epistemic uncertainty); and (3) the possibility that our forecasts are wrong owing to behaviors of volcanic processes about which we are completely ignorant and, hence, cannot quantify in terms of probabilities (ontological error). Here we put forward a probabilistic framework for hazard analysis recently proposed by Marzocchi and Jordan (2014), which unifies the treatment of all three types of uncertainty. Within this framework, an eruption forecasting or a volcanic hazard model is said to be complete only if it (a) fully characterizes the epistemic uncertainties in the model's representation of aleatory variability and (b) can be unconditionally tested (in principle) against observations to identify ontological errors. Unconditional testability, which is the key to model validation, hinges on an experimental concept that characterizes hazard events in terms of exchangeable data sequences with well-defined frequencies. We illustrate the application of this unified probabilistic framework by describing experimental concepts for the forecasting of tephra fall from Campi Flegrei. Eventually, this example may serve as a guide for the application of the same probabilistic framework to other natural hazards.

2021 ◽  
Author(s):  
Warner Marzocchi ◽  
Jacopo Selva ◽  
Thomas H. Jordan

Abstract. The main purpose of this article is to emphasize the importance of clarifying the probabilistic framework adopted for volcanic hazard and eruption forecasting. Eruption forecasting and volcanic hazard analysis seeks to quantify the deep uncertainties that pervade the modeling of pre-, sin- and post-eruptive processes. These uncertainties can be differentiated into three fundamental types: (1) the natural variability of volcanic systems, usually represented as stochastic processes with parameterized distributions (aleatory variability); (2) the uncertainty in our knowledge of how volcanic systems operate and evolve, often represented as subjective probabilities based on expert opinion (epistemic uncertainty); and (3) the possibility that our forecasts are wrong owing to behaviors of volcanic processes about which we are completely ignorant and, hence, cannot quantify in terms of probabilities (ontological error). Here we put forward a probabilistic framework for hazard analysis recently proposed by Marzocchi & Jordan (2014), which unifies the treatment of all three types of uncertainty. Within this framework, an eruption forecasting or a volcanic hazard model is said to be complete only if it (a) fully characterizes the epistemic uncertainties in the model's representation of aleatory variability and (b) can be unconditionally tested (in principle) against observations to identify ontological errors. Unconditional testability, which is the key to model validation, hinges on an experimental concept that characterizes hazard events in terms of exchangeable data sequences with well-defined frequencies. We illustrate the application of this unified probabilistic framework by describing experimental concepts for the forecasting of tephra fall from Campi Flegrei. Eventually, this example may serve as a guide for the application of the same probabilistic framework to other natural hazards.


2020 ◽  
Author(s):  
Mark Novakovic ◽  
Emrah Yenier ◽  
Andrew Hovey ◽  
Joseph Quinn ◽  
Benjamin Witten

<p>A seismic hazard analysis was conducted for a site in Papua New Guinea which is located in a seismically-active region that experiences frequent large earthquakes generated by crustal and subduction sources.  A suite of ground motion prediction equations (GMPEs) was developed for each source type (crustal, interface and in-slab) using the scaled-backbone approach.  To this end, a ground-motion database consisting of events of 4.0<Mw<8.0 was compiled from available local and regional monitoring stations.  Ground motions were classified based on the source type and converted to a common reference site condition.  The site-corrected motions were compared against alternative GMPEs to examine residual trends between observed and predicted amplitudes.  A backbone model that represents the best estimate of the median ground motions for each source type was selected.  The backbone models were then adjusted to the median of the ground motions observed at the study site.</p><p>The epistemic uncertainty in median predictions was modeled using a logic-tree approach, where the distribution of potential median predictions is approximated by a lower, central and upper model.  The central model is represented by the site-adjusted backbone model; it was scaled to define the lower and upper branches.  The scaling factor was determined considering: (i) the standard deviation in median prediction of alternative GMPEs; and (ii) epistemic uncertainties recommended in other studies.  The available data were insufficient to model aleatory variability with confidence; therefore, the standard deviation of observed motions in data-rich regions is used for guidance.  Two alternative aleatory variability models (ergodic and single-station sigma) adopted from other studies are recommended.</p>


2020 ◽  
Vol 18 (14) ◽  
pp. 6119-6148
Author(s):  
Graeme Weatherill ◽  
Fabrice Cotton

Abstract Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.


2019 ◽  
Vol 19 (10) ◽  
pp. 2097-2115 ◽  
Author(s):  
Panjamani Anbazhagan ◽  
Ketan Bajaj ◽  
Karanpreet Matharu ◽  
Sayed S. R. Moustafa ◽  
Nassir S. N. Al-Arifi

Abstract. Peak ground acceleration (PGA) and study area (SA) distribution for the Patna district are presented considering both the classical and zoneless approaches through a logic tree framework to capture the epistemic uncertainty. Seismicity parameters are calculated by considering completed and mixed earthquake data. Maximum magnitude is calculated using three methods, namely the incremental method, Kijko method, and regional rupture characteristics approach. The best suitable ground motion prediction equations (GMPEs) are selected by carrying out an “efficacy test” using log likelihood. Uniform hazard response spectra have been compared with Indian standard BIS 1893. PGA varies from 0.38 to 0.30 g from the southern to northern periphery considering 2 % probability of exceedance in 50 years.


2004 ◽  
Vol 31 (3-4) ◽  
pp. 345-360 ◽  
Author(s):  
Sreten Mastilovic ◽  
Dusan Krajcinovic

The present review focuses on the plane strain problem of high strain rate expansion of a cylindrical cavity within an infinite brittle material with random microstructure. The material is represented by an ensemble of "continuum particles" forming a two-dimensional geometrically and structurally disordered lattice. The proposed model includes the aleatory variability and epistemic uncertainty of the process. The dynamic particle simulations are performed at seven different cavity expansion rates. The resulting damage evolution process is non-stationary, non-local, and non-equilibrium. This problem, therefore, belongs to the class of phenomena for which the traditional continuum models are not well suited, and detailed experimental data are either difficult to get or not available at all. The present study explores the potential role of the particle dynamics in addressing these problems. .


2015 ◽  
Vol 31 (2) ◽  
pp. 661-698 ◽  
Author(s):  
Julian J. Bommer ◽  
Kevin J. Coppersmith ◽  
Ryan T. Coppersmith ◽  
Kathryn L. Hanson ◽  
Azangi Mangongolo ◽  
...  

A probabilistic seismic hazard analysis has been conducted for a potential nuclear power plant site on the coast of South Africa, a country of low-to-moderate seismicity. The hazard study was conducted as a SSHAC Level 3 process, the first application of this approach outside North America. Extensive geological investigations identified five fault sources with a non-zero probability of being seismogenic. Five area sources were defined for distributed seismicity, the least active being the host zone for which the low recurrence rates for earthquakes were substantiated through investigations of historical seismicity. Empirical ground-motion prediction equations were adjusted to a horizon within the bedrock at the site using kappa values inferred from weak-motion analyses. These adjusted models were then scaled to create new equations capturing the range of epistemic uncertainty in this region with no strong motion recordings. Surface motions were obtained by convolving the bedrock motions with site amplification functions calculated using measured shear-wave velocity profiles.


2021 ◽  
pp. 875529302110560
Author(s):  
Yousef Bozorgnia ◽  
Norman A Abrahamson ◽  
Sean K Ahdi ◽  
Timothy D Ancheta ◽  
Linda Al Atik ◽  
...  

This article summarizes the Next Generation Attenuation (NGA) Subduction (NGA-Sub) project, a major research program to develop a database and ground motion models (GMMs) for subduction regions. A comprehensive database of subduction earthquakes recorded worldwide was developed. The database includes a total of 214,020 individual records from 1,880 subduction events, which is by far the largest database of all the NGA programs. As part of the NGA-Sub program, four GMMs were developed. Three of them are global subduction GMMs with adjustment factors for up to seven worldwide regions: Alaska, Cascadia, Central America and Mexico, Japan, New Zealand, South America, and Taiwan. The fourth GMM is a new Japan-specific model. The GMMs provide median predictions, and the associated aleatory variability, of RotD50 horizontal components of peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral acceleration (PSA) at oscillator periods ranging from 0.01 to 10 s. Three GMMs also quantified “within-model” epistemic uncertainty of the median prediction, which is important in regions with sparse ground motion data, such as Cascadia. In addition, a damping scaling model was developed to scale the predicted 5%-damped PSA of horizontal components to other damping ratios ranging from 0.5% to 30%. The NGA-Sub flatfile, which was used for the development of the NGA-Sub GMMs, and the NGA-Sub GMMs coded on various software platforms, have been posted for public use.


2019 ◽  
Author(s):  
Panjamani Anbazhagan ◽  
Ketan Bajaj ◽  
Karanpreet Matharu ◽  
Sayed S. R. Moustafa ◽  
Nassir S. N. Al-Arifi

Abstract. PGA and SA distribution for Patna district is presented considering both classical and zoneless approach through logic tree frame work to capture the epistemic uncertainty. Seismicity parameters are calculated by considering completed and mixed earthquake data. Maximum magnitude was calculated using three methods namely incremental method, Kijko method and regional rupture characteristics approach. Best suitable GMPE was selected by carrying out efficacy test using log likelihood. Uniform hazard response spectra have been compared with Indian standard BIS 1893. PGA varies from 0.38 g to 0.30 g from southern to northern periphery considering 2 % probability of exceedence in 50 years.


Sign in / Sign up

Export Citation Format

Share Document