Probabilistic seismic hazard assessment for South-Eastern France

The accurate evaluation and appropriate treatment of uncertainties is of primary importance in modern probabilistic seismic hazard assessment (PSHA). One of the objectives of the SIGMA project was to establish a framework to improve knowledge and data on two target regions characterized by low-to-moderate seismic activity. In this paper, for South-Eastern France, we present the final PSHA performed within the SIGMA project. A new earthquake catalogue for France covering instrumental and historical periods was used for the calculation of the magnitude-frequency distributions. The hazard model incorporates area sources, smoothed seismicity and a 3D faults model. A set of recently developed ground motion prediction equations (GMPEs) from global and regional data, evaluated as adequately representing the ground motion characteristics in the region, was used to calculate the hazard. The magnitude-frequency distributions, maximum magnitude, faults slip rate and style-of-faulting are considered as additional source of epistemic uncertainties. The hazard results for generic rock condition (Vs30 = 800 m/s) are displayed for 20 sites in terms of uniform hazard spectra at two return periods (475 years and 10,000 years). The contributions of the epistemic uncertainties in the ground motion characterizations and in the seismic source characterization to the total hazard uncertainties are analyzed. Finally, we compare the results with existing models developed at national scale in the framework of the first generation of models supporting the Eurocode 8 enforcement, (MEDD 2002 and AFPS06) and at the European scale (within the SHARE project), highlighting significant discrepancies at short return periods.     

Determination of GMPE functional form for an active region with limited strong motion data: application to the Himalayan region

Advancement in the seismic networks results in formulation of different functional forms for developing any new ground motion prediction equation (GMPE) for a region. Till date, various guidelines and tools are available for selecting a suitable GMPE for any seismic study area. However, these methods are efficient in quantifying the GMPE but not for determining a proper functional form and capturing the epistemic uncertainty associated with selection of GMPE. In this study, the compatibility of the recent available functional forms for the active region is tested for distance and magnitude scaling. Analysis is carried out by determining the residuals using the recorded and the predicted spectral acceleration values at different periods. Mixed effect regressions are performed on the calculated residuals for determining the intra- and interevent residuals. Additionally, spatial correlation is used in mixed effect regression by changing its likelihood function. Distance scaling and magnitude scaling are respectively examined by studying the trends of intraevent residuals with distance and the trend of the event term with magnitude. Further, these trends are statistically studied for a respective functional form of a ground motion. Additionally, genetic algorithm and Monte Carlo method are used respectively for calculating the hinge point and standard error for magnitude and distance scaling for a newly determined functional form. The whole procedure is applied and tested for the available strong motion data for the Himalayan region. The functional form used for testing are five Himalayan GMPEs, five GMPEs developed under NGA-West 2 project, two from Pan-European, and one from Japan region. It is observed that bilinear functional form with magnitude and distance hinged at 6.5 M _w and 300 km respectively is suitable for the Himalayan region. Finally, a new regression coefficient for peak ground acceleration for a suitable functional form that governs the attenuation characteristic of the Himalayan region is derived.     

Comparison of conventional and Monte Carlo simulation-based probabilistic seismic hazard analyses for Shiraz city,southern Iran

Estimation of ground-motion amplitudes of different hazard levels is of paramount importance in planning of urban development of any metropolis. Such estimation can be computed through a probabilistic seismic hazard analysis (PSHA). This paper concentrates on the PSHA of an area located in Shiraz city, southern Iran. The area includes whole of Shiraz city (i.e., one of the largest and most populous cities of Iran) and its outskirts. Conventional and Monte Carlo simulation-based approaches are utilized to perform the PSHA of the studied area. Two areal seismic source models are delineated, and thence seismicity parameters of all zones associated with their corresponding uncertainties are computed. Uncertainties in ground-motion prediction are accounted for via three ground-motion prediction equations (GMPEs) within the logic tree framework. These GMPEs are applied to estimate bedrock ground shaking (V_s30?=?760 m/s) for several return periods (i.e., 75, 475, 975, and 2475 years). In general, the results of the two abovementioned PSHA approaches show relatively similar results. However, the Monte Carlo simulation-based approach overpredicts bedrock spectral accelerations at periods of 0.4–2.5 s compared to the conventional PSHA approach for return periods of 475, 975, and 2475 years.     

Geodetic Strain Observations and Return Period of the Strongest Earthquakes of a Given Seismic Source Zone

We present the basis for a method for estimating the return period of large and medium earthquakes that is independent of current deterministic and probabilistic approaches. The two standard techniques of seismic hazard assessment??probabilistic seismic hazard assessment (PSHA) and deterministic seismic hazard assessment (DSHA)??suffer from limited knowledge of seismic prehistory. A further weakness of PSHA is its requirement of homogeneous seismic activity within a seismic zone. Moreover, PSHA and DSHA were developed for seismically active areas and, thus, cannot reliably be used in areas of medium and low activity. In this paper we propose the combined use of geodetic strain rate data and the seismic moment data set determined for past seismic events. This combination represents a new and independent approach to estimation of future seismic activity. Using a modified version of Kostrov??s (Phys Solid Earth 1:23?C40, 1974) equation and the catalogue of seismic moments, the minimum return period of the strongest earthquakes of a source area is estimated.     

Physics‐based probabilistic seismic hazard and loss assessment in large urban areas: A simplified application to Istanbul

A set of 3D physics‐based numerical simulations (PBS) of possible earthquakes scenarios in Istanbul along the North Anatolian Fault (Turkey) is considered in this article to provide a comprehensive example of application of PBS to probabilistic seismic hazard (PSHA) and loss assessment in a large urban area. To cope with the high‐frequency (HF) limitations of PBS, numerical results are first postprocessed by a recently introduced technique based on Artificial Neural Networks (ANN), providing broadband waveforms with a proper correlation of HF and low‐frequency (LF) portions of ground motion as well as a proper spatial correlation of peak values also at HF, that is a key feature for the seismic risk application at urban scale. Second, before application to PSHA, a statistical analysis of residuals is carried out to ensure that simulated results provide a set of realizations with a realistic within‐ and between‐event variability of ground motion. PBS results are then applied in a PSHA framework, adopting both the “generalized attenuation function” (GAF) approach, and a novel “footprint” (FP)‐based approach aiming at a convenient and direct application of PBS into PSHA. PSHA results from both approaches are then compared with those obtained from a more standard application of PSHA with empirical ground motion models. Finally, the probabilistic loss assessment of an extended simplified portfolio of buildings is investigated, comparing the results obtained adopting the different approaches: (i) GMPE, (ii) GAF, and (iii) FP. Only FP turned out to have the capability to account for the specific features of source and propagation path, while preserving the proper physically based spatial correlation characteristics, as required for a reliable loss estimate on a building portfolio spatially distributed over a large urban area.     

Site specific probabilistic seismic hazard analysis at Dubai Creek on the west coast of UAE

A probabilistic seismic hazard analysis (PSHA) was conducted to establish the hazard spectra for a site located at Dubai Creek on the west coast of the United Arab Emirates (UAE). The PSHA considered all the seismogenic sources that affect the site, including plate boundaries such as the Makran subduction zone, the Zagros fold-thrust region and the transition fault system between them; and local crustal faults in UAE. PSHA indicated that local faults dominate the hazard. The peak ground acceleration (PGA) for the 475-year return period spectrum is 0.17 g and 0.33 g for the 2,475-year return period spectrum. The hazard spectra are then employed to establish rock ground motions using the spectral matching technique.     

A probabilistic seismic hazard assessment for the Turkish territory: part II—fault source and background seismicity model

Over the years, several local and regional seismic hazard studies have been conducted for the estimation of the seismic hazard in Turkey using different statistical processing tools for instrumental and historical earthquake data and modeling the geologic and tectonic characteristics of the region. Recently developed techniques, increased knowledge and improved databases brought the necessity to review the national active fault database and the compiled earthquake catalogue for the development of a national earthquake hazard map. A national earthquake strategy and action plan were conceived and accordingly with the collaboration of the several institutions and expert researchers, the Revision of Turkish Seismic Hazard Map Project (UDAP-Ç-13-06) was initiated, and finalized at the end of 2014. The scope of the project was confined to the revision of current national seismic hazard map, using the state of the art technologies and knowledge of the active fault, earthquake database, and ground motion prediction equations. The following two seismic source zonation models are developed for the probabilistic earthquake hazard analysis: (1) Area source model, (2) Fault and spatial smoothing seismic source model (FSBCK). In this study, we focus on the development and the characterization of the Fault Source model, the background spatially smoothed seismicity model and intrinsic uncertainty on the earthquake occurrence-rates-estimation. Finally, PSHA results obtained from the fault and spatial smoothed seismic source model are presented for 43, 72, 475 and 2475 years return periods (corresponding to 69, 50, 10, and 2% probability of exceedance in 50 years) for PGA and 5% damped spectral accelerations at 0.2 and 1.0 s.     

Probabilistic seismic hazard analysis: Early history

Probabilistic seismic hazard analysis (PSHA) is the evaluation of annual frequencies of exceedence of ground motion levels (typically designated by peak ground acceleration or by spectral accelerations) at a site. The result of a PSHA is a seismic hazard curve (annual frequency of exceedence vs ground motion amplitude) or a uniform hazard spectrum (spectral amplitude vs structural period, for a fixed annual frequency of exceedence). Analyses of this type were first conceived in the 1960s and have become the basis for the seismic design of engineered facilities ranging from common buildings designed according to building codes to critical facilities such as nuclear power plants. This Historical Note traces the early history of PSHA. Copyright © 2007 John Wiley & Sons, Ltd.     

Ground motion prediction equation considering combined dataset of recorded and simulated ground motions

Himalayan region is one of the most active seismic regions in the world and many researchers have highlighted the possibility of great seismic event in the near future due to seismic gap. Seismic hazard analysis and microzonation of highly populated places in the region are mandatory in a regional scale. Region specific Ground Motion Predictive Equation (GMPE) is an important input in the seismic hazard analysis for macro- and micro-zonation studies. Few GMPEs developed in India are based on the recorded data and are applicable for a particular range of magnitudes and distances. This paper focuses on the development of a new GMPE for the Himalayan region considering both the recorded and simulated earthquakes of moment magnitude 5.3–8.7. The Finite Fault simulation model has been used for the ground motion simulation considering region specific seismotectonic parameters from the past earthquakes and source models. Simulated acceleration time histories and response spectra are compared with available records. In the absence of a large number of recorded data, simulations have been performed at unavailable locations by adopting Apparent Stations concept. Earthquakes recorded up to 2007 have been used for the development of new GMPE and earthquakes records after 2007 are used to validate new GMPE. Proposed GMPE matched very well with recorded data and also with other highly ranked GMPEs developed elsewhere and applicable for the region. Comparison of response spectra also have shown good agreement with recorded earthquake data. Quantitative analysis of residuals for the proposed GMPE and region specific GMPEs to predict Nepal–India 2011 earthquake of Mw of 5.7 records values shows that the proposed GMPE predicts Peak ground acceleration and spectral acceleration for entire distance and period range with lower percent residual when compared to exiting region specific GMPEs.     

Probabilistic seismic hazard assessment of Italy using kernel estimation methods

A representation of seismic hazard is proposed for Italy based on the zone-free approach developed by Woo (BSSA 86(2):353–362, 1996a), which is based on a kernel estimation method governed by concepts of fractal geometry and self-organized seismicity, not requiring the definition of seismogenic zoning. The purpose is to assess the influence of seismogenic zoning on the results obtained for the probabilistic seismic hazard analysis (PSHA) of Italy using the standard Cornell’s method. The hazard has been estimated for outcropping rock site conditions in terms of maps and uniform hazard spectra for a selected site, with 10 % probability of exceedance in 50 years. Both spectral acceleration and spectral displacement have been considered as ground motion parameters. Differences in the results of PSHA between the two methods are compared and discussed. The analysis shows that, in areas such as Italy, characterized by a reliable earthquake catalog and in which faults are generally not easily identifiable, a zone-free approach can be considered a valuable tool to address epistemic uncertainty within a logic tree framework.     

Probabilistic seismic hazard assessment of the Pyrenean region

A unified probabilistic seismic hazard assessment (PSHA) for the Pyrenean region has been performed by an international team composed of experts from Spain and France during the Interreg IIIA ISARD project. It is motivated by incoherencies between the seismic hazard zonations of the design codes of France and Spain and by the need for input data to be used to define earthquake scenarios. A great effort was invested in the homogenisation of the input data. All existing seismic data are collected in a database and lead to a unified catalogue using a local magnitude scale. PSHA has been performed using logic trees combined with Monte Carlo simulations to account for both epistemic and aleatory uncertainties. As an alternative to hazard calculation based on seismic sources zone models, a zoneless method is also used to produce a hazard map less dependant on zone boundaries. Two seismogenic source models were defined to take into account the different interpretations existing among specialists. A new regional ground-motion prediction equation based on regional data has been proposed. It was used in combination with published ground-motion prediction equations derived using European and Mediterranean data. The application of this methodology leads to the definition of seismic hazard maps for 475- and 1,975-year return periods for spectral accelerations at periods of 0 (corresponding to peak ground acceleration), 0.1, 0.3, 0.6, 1 and 2 s. Median and percentiles 15% and 85% acceleration contour lines are represented. Finally, the seismic catalogue is used to produce a map of the maximum acceleration expected for comparison with the probabilistic hazard maps. The hazard maps are produced using a grid of 0.1°. The results obtained may be useful for civil protection and risk prevention purposes in France, Spain and Andorra.     

Regional adjustment factors for three NGA-West2 ground-motion prediction equations to be applicable in northern Iran

Ground-motion prediction equations (GMPEs) are essential tools in seismic hazard studies to estimate ground motions generated by potential seismic sources. Global GMPEs which are based on well-compiled global strong-motion databanks, have certain advantages over local GMPEs, including more sophisticated parameters in terms of distance, faulting style, and site classification but cannot guarantee the local/region-specific propagation characteristics of shear wave (e.g., geometric spreading behavior, quality factor) for different seismic regions at larger distances (beyond about 80 km). Here, strong-motion records of northern Iran have been used to estimate the propagation characteristics of shear wave and determine the region-specific adjustment parameters for three of the NGA-West2 GMPEs to be applicable in northern Iran. The dataset consists of 260 three-component records from 28 earthquakes, recorded at 139 stations, with moment magnitudes between 4.9 and 7.4, horizontal distance to the surface projection of the rupture (R _JB) less than 200 km, and average shear-wave velocity over the top 30 m of the subsurface (V _S30) between 155 and 1500 m/s. The paper also presents the ranking results for three of the NGA-West2 GMPEs against strong motions recorded in northern Iran, before and after adjustment for region-dependent attenuation characteristics. The ranking is based on the likelihood and log-likelihood methods (LH and LLH) proposed by Scherbaum et al. (Bull Seismol Soc Am 94: 2164–2185, 2004, Bull Seismol Soc Am 99, 3234–3247, 2009, respectively), the Nash–Sutcliffe model efficiency coefficient (Nash and Sutcliffe, J Hydrol 10:282–290, 1970), and the EDR method of Kale and Akkar (Bull Seismol Soc Am 103:1069–1084, 2012). The best-fitting models over the whole frequency range are the ASK14 and BSSA14 models. Taking into account that the models’ performances were boosted after applying the adjustment factors, at least moderate regional variation of ground motions is highlighted. The regional adjustment based on the Iranian database reveals an upward trend (indicated as high Q factor) for the selected database. Further investigation to determine adjustment factors based on a much richer database of the Iranian strong-motion records is of utmost important for seismic hazard and risk analysis studies in northern Iran, containing major cities including the capital city of Tehran.     

Seismic hazard assessment of Tehran,Iran with emphasis on near-fault rupture directivity effects

Many destructive earthquakes happened in Tehran, Iran in the last centuries. The existence of active faults like the North Tehran is the main cause of seismicity in this city. According to past investigations, it is estimated that in the scenario of activation of the North Tehran fault, many structures in Tehran will collapse. Therefore, it is necessary to incorporate the near field rupture directivity effects of this fault into the seismic hazard assessment of important sites in Tehran. In this study, using calculations coded in MATLAB,Probabilistic Seismic Hazard Analysis(PSHA) is conducted for an important site in Tehran. Following that, deaggregation technique is performed on PSHA and the contribution of seismic scenarios to hazard is obtained in the range of distance and magnitude. After identifying the North Tehran fault as the most hazardous source affecting the site in 10000-year return period, rupture directivity effects of this fault is incorporated into the seismic hazard assessment using Somerville et al.(1997) model with broadband approach and Shahi and Baker(2011) model with narrowband approach. The results show that the narrowband approach caused a 27% increase in the peak of response spectrum in 10000-year return period compared with the conventional PSHA. Therefore, it is necessary to incorporate the near fault rupture directivity effects into the higher levels of seismic hazard assessment attributed to important sites.     

Seismic Hazard Assessment: Issues and Alternatives

Seismic hazard and risk are two very important concepts in engineering design and other policy considerations. Although seismic hazard and risk have often been used interchangeably, they are fundamentally different. Furthermore, seismic risk is more important in engineering design and other policy considerations. Seismic hazard assessment is an effort by earth scientists to quantify seismic hazard and its associated uncertainty in time and space and to provide seismic hazard estimates for seismic risk assessment and other applications. Although seismic hazard assessment is more a scientific issue, it deserves special attention because of its significant implication to society. Two approaches, probabilistic seismic hazard analysis (PSHA) and deterministic seismic hazard analysis (DSHA), are commonly used for seismic hazard assessment. Although PSHA has been proclaimed as the best approach for seismic hazard assessment, it is scientifically flawed (i.e., the physics and mathematics that PSHA is based on are not valid). Use of PSHA could lead to either unsafe or overly conservative engineering design or public policy, each of which has dire consequences to society. On the other hand, DSHA is a viable approach for seismic hazard assessment even though it has been labeled as unreliable. The biggest drawback of DSHA is that the temporal characteristics (i.e., earthquake frequency of occurrence and the associated uncertainty) are often neglected. An alternative, seismic hazard analysis (SHA), utilizes earthquake science and statistics directly and provides a seismic hazard estimate that can be readily used for seismic risk assessment and other applications.     

Assessment of ground motion variability and its effects on seismic hazard analysis: a case study for iceland

Probabilistic seismic hazard analysis (PSHA) generally relies on the basic assumption that ground motion prediction equations (GMPEs) developed for other similar tectonic regions can be adopted in the considered area. This implies that observed ground motion and its variability at considered sites could be modelled by the selected GMPEs. Until now ground-motion variability has been taken into account in PSHA by integrating over the standard deviation reported in GMPEs, which significantly affects estimated ground motions, especially at very low probabilities of exceedance. To provide insight on this issue, ground-motion variability in the South Iceland Seismic Zone (SISZ), where many ground-motion records are available, is assessed. Three statistical methods are applied to separate the aleatory variability into source (inter-event), site (inter-site) and residual (intra-event and intra-site) components. Furthermore, the current PSHA procedure that makes the ergodic assumption of equality between spatially and temporal variability is examined. In contrast to the ergodic assumption, several recent studies show that the observed ground-motion variability at an individual location is lower than that implied by the standard deviation of a GMPE. This could imply a mishandling of aleatory uncertainty in PSHA by ignoring spatial variability and by mixing aleatory and epistemic uncertainties in the computation of sigma. Station correction coefficients are introduced in order to capture site effects at different stations. The introduction of the non-ergodic assumption in PSHA leads to larger epistemic uncertainty, although this is not the same as traditional epistemic uncertainty modelled using different GMPEs. The epistemic uncertainty due to the site correction coefficients (i.e. mean residuals) could be better constrained for future events if more information regarding the characteristics of these seismic sources and path dependence could be obtained.     

The effect of causal parameter bounds in PSHA‐based ground motion selection

In this paper the effect of causal parameter bounds (e.g. magnitude, source‐to‐site distance, and site condition) on ground motion selection, based on probabilistic seismic hazard analysis (PSHA) results, is investigated. Despite the prevalent application of causal parameter bounds in ground motion selection, present literature on the topic is cast in the context of a scenario earthquake of interest, and thus specific bounds for use in ground motion selection based on PSHA, and the implications of such bounds, is yet to be examined. Thirty‐six PSHA cases, which cover a wide range of causal rupture deaggregation distributions and site conditions, are considered to empirically investigate the effects of various causal parameter bounds on the characteristics of selected ground motions based on the generalized conditional intensity measure (GCIM) approach. It is demonstrated that the application of relatively ‘wide’ bounds on causal parameters effectively removes ground motions with drastically different characteristics with respect to the target seismic hazard and results in an improved representation of the target causal parameters. In contrast, the use of excessively ‘narrow’ bounds can lead to ground motion ensembles with a poor representation of the target intensity measure distributions, typically as a result of an insufficient number of prospective ground motions. Quantitative criteria for specifying bounds for general PSHA cases are provided, which are expected to be sufficient in the majority of problems encountered in ground motion selection for seismic demand analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Ground-motion prediction equations for interface earthquakes of M7 to M9 based on empirical data from Japan

Ground motion prediction equations (GMPEs) have a major impact on seismic hazard estimates, because they control the predicted amplitudes of ground shaking. The prediction of ground-motion amplitudes due to mega-thrust earthquakes in subduction zones has been hampered by a paucity of empirical ground-motion data for the very large magnitudes (moment magnitude (M) $>$ 7) of most interest to hazard analysis. Recent data from Tohoku M9.0 2011 earthquake are important in this regard, as this is the largest well-recorded subduction event, and the only such event with sufficient data to enable a clear separation of the overall source, path and site effects. In this study, we use strong-ground-motion records from the M9 Tohoku event to derive an event-specific GMPE. We then extend this M9 GMPE to represent the shaking from other M $>$ 7 interface events in Japan by adjusting the source term. We focus on events in Japan to reduce ambiguity that results when combining data in different regions having different source, path and site effect attributes. Source levels (adjustment factors) for other Japanese events are determined as the average residuals of ground-motions with respect to the Tohoku GMPE, keeping all other coefficients fixed. The mean residuals (source terms) scale most steeply with magnitude at the lower frequencies; this is in accord with expectations based on overall source-scaling concepts. Interpolating source terms over the magnitude range of 7.0–9.0, we produce a GMPE for large interface events of M7–M9, for NEHRP B/C boundary site conditions (time-averaged shear-wave velocity of 760 m/s over the top 30 m) in both fore-arc and back-arc regions of Japan. We show how these equations may be adjusted to account for the deeper soil profiles (for the same value of $\hbox {V}_\mathrm{S30})$ in western North America. The proposed GMPE predicts lower motions at very long periods, higher motions at short periods, and similar motions at intermediate periods, relative to the simulation-based GMPE model of Atkinson and Macias (2009) for the Cascadia subduction zone.     

Seismic hazard assessment for Iceland in terms of macroseismic intensity using a site approach

We present the results of probabilistic seismic hazard assessment for Iceland in the framework of the EU project UPStrat-MAFA using the so-called site approach implemented in the SASHA computational code. This approach estimates seismic hazard in terms of macroseismic intensity by basically relying on local information about documented effects of past seismic events in the framework of a formally coherent and complete treatment of intensity data. In the case of Iceland, due to the lack of observed intensities for past earthquakes, local seismic histories were built using indirect macroseismic estimates deduced from epicentral information through an empirical attenuation relationship in probabilistic form. Seismic hazard was computed for four exceedance probabilities for an exposure time of 50 years, equivalent to average return periods of 50, 200, 475 and 975 years. For some localities, further return periods were examined and deaggregation analysis was performed. Results appear significantly different from previous seismic hazard maps, though just a semi-qualitative comparison is possible because of the different shaking measure considered (peak ground acceleration versus intensity), and the different computational methodology and input data used in these studies.     

Uncertainty Analysis and Expert Judgment in Seismic Hazard Analysis

The large uncertainty associated with the prediction of future earthquakes is usually regarded as the main reason for increased hazard estimates which have resulted from some recent large scale probabilistic seismic hazard analysis studies (e.g. the PEGASOS study in Switzerland and the Yucca Mountain study in the USA). It is frequently overlooked that such increased hazard estimates are characteristic for a single specific method of probabilistic seismic hazard analysis (PSHA): the traditional (Cornell?CMcGuire) PSHA method which has found its highest level of sophistication in the SSHAC probability method. Based on a review of the SSHAC probability model and its application in the PEGASOS project, it is shown that the surprising results of recent PSHA studies can be explained to a large extent by the uncertainty model used in traditional PSHA, which deviates from the state of the art in mathematics and risk analysis. This uncertainty model, the Ang?CTang uncertainty model, mixes concepts of decision theory with probabilistic hazard assessment methods leading to an overestimation of uncertainty in comparison to empirical evidence. Although expert knowledge can be a valuable source of scientific information, its incorporation into the SSHAC probability method does not resolve the issue of inflating uncertainties in PSHA results. Other, more data driven, PSHA approaches in use in some European countries are less vulnerable to this effect. The most valuable alternative to traditional PSHA is the direct probabilistic scenario-based approach, which is closely linked with emerging neo-deterministic methods based on waveform modelling.