Ground-motion prediction equations (GMPEs) for inelastic response and structural behaviour factors

The objective of this paper is to present ground-motion prediction equations describing constant-ductility inelastic spectral ordinates and structural behaviour factors. These equations are intended for application within the framework of Eurocode 8. Most of the strong-motion data used in the present work is obtained from the ISESD (Internet Site for European Strong-motion Data) databank. Present analysis includes ground motion records from significant Icelandic earthquakes, which are augmented by records obtained from continental Europe and the Middle East. In all cases the selected ground motion records are generated during shallow earthquakes within a distance of 100 km from the recording station. The classification of site conditions in the present work is based on the Eurocode 8 definition.     

Comparing predicted and observed ground motions from subduction earthquakes in the Lesser Antilles

This brief article presents a quantitative analysis of the ability of eight published empirical ground-motion prediction equations (GMPEs) for subduction earthquakes (interface and intraslab) to estimate observed earthquake ground motions on the islands of the Lesser Antilles (specifically Guadeloupe, Martinique, Trinidad, and Dominica). In total, over 300 records from 22 earthquakes from various seismic networks are used within the analysis. It is found that most of the GMPEs tested perform poorly, which is mainly due to a larger variability in the observed ground motions than predicted by the GMPEs, although two recent GMPEs derived using Japanese strong-motion data provide reasonably good predictions. Analyzing separately the interface and intraslab events does not significant modify the results. Therefore, it is concluded that seismic hazard assessments for this region should use a variety of GMPEs in order to capture this large epistemic uncertainty in earthquake ground-motion prediction for the Lesser Antilles.     

A strong-motion database from the Peru–Chile subduction zone

Earthquake hazard along the Peru–Chile subduction zone is amongst the highest in the world. The development of a database of subduction-zone strong-motion recordings is, therefore, of great importance for ground-motion prediction in this region. Accelerograms recorded by the different networks operators in Peru and Chile have been compiled and processed in a uniform manner, and information on the source parameters of the causative earthquakes, fault-plane geometries and local site conditions at the recording stations has been collected and reviewed to obtain high-quality metadata. The compiled database consists of 98 triaxial ground-motion recordings from 15 subduction-type events with moment magnitudes ranging from 6.3 to 8.4, recorded at 59 different sites in Peru and Chile, between 1966 and 2007. While the database presented in this study is not sufficient for the derivation of a new predictive equation for ground motions from subduction events in the Peru–Chile region, it significantly expands the global database of strong-motion data and associated metadata that can be used in the derivation of predictive equations for subduction environments. Additionally, the compiled database will allow the assessment of existing predictive models for subduction-type events in terms of their suitability for the Peru–Chile region, which directly influences seismic hazard assessment in this region.     

Hazard-consistent response spectra in the Region of Murcia (Southeast Spain): comparison to earthquake-resistant provisions

Hazard-consistent ground-motion characterisations of three representative sites located in the Region of Murcia (southeast Spain) are presented. This is the area where the last three damaging events in Spain occurred and there is a significant amount of data for comparing them with seismic hazard estimates and earthquake-resistant provisions. Results of a probabilistic seismic hazard analysis are used to derive uniform hazard spectra (UHS) for the 475-year return period, on rock and soil conditions. Hazard deaggregation shows that the largest hazard contributions are due to small, local events for short-period target motions and to moderate, more distant events for long-period target motions. For each target motion and site considered, the associated specific response spectra (SRS) are obtained. It is shown that the combination of two SRS, for short- and long-period ground motions respectively, provides a good approximation to the UHS at each site. The UHS are compared to design response spectra contained in current Spanish and European seismic codes for the 475-year return period. For the three sites analysed, only the Eurocode 8 (EC8) type 2 spectrum captures the basic shape of the UHS (and not the EC8 type 1, as could be expected a priori). An alternative response spectrum, anchored at short- and long-period accelerations, is tested, providing a close match to the UHS spectra at the three sites. Results underline the important contribution of the frequent, low-to-moderate earthquakes that characterize the seismicity of this area to seismic hazard (at the 475-year return period).     

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.     

Consistency of ground-motion predictions from the past four decades

Due to the limited observational datasets available for the derivation of ground-motion prediction equations (GMPEs) there is always epistemic uncertainty in the estimated median ground motion. Because of the increasing quality and quantity of strong-motion datasets it would be expected that the epistemic uncertainty in ground-motion prediction (related to lack of knowledge and data) is decreasing. In this study the predicted median ground motions from over 200 GMPEs for various scenarios are plotted against date of publication to examine whether the scatter in the predictions (a measure of epistemic uncertainty) is decreasing with time. It is found that there are still considerable differences in predicted ground motions from the various GMPEs and that the variation between estimates is not reducing although the ground motion estimated by averaging median predictions is roughly constant. For western North America predictions for moderate earthquakes have show a high level of consistency since the 1980s as do, but to a lesser extent, predictions for moderate earthquakes in Europe, the Mediterranean and the Middle East. A good match is observed between the predictions from GMPEs and the median ground motions based on observations from similar scenarios. Variations in median ground motion predictions for stable continental regions and subduction zones from different GMPEs are large, even for moderate earthquakes. The large scatter in predictions of the median ground motion shows that epistemic uncertainty in ground-motion prediction is still large and that it is vital that this is accounted for in seismic hazard assessments.     

A procedure on ground motion selection and scaling for nonlinear response of simple structural systems

This study presents a ground-motion selection and scaling methodology that preserves the basic seismological features of the scaled records with reduced scatter in the nonlinear structural response. The methodology modifies each strong-motion recording with known fundamental seismological parameters using the estimations of ground-motion prediction equations for a given target hazard level. It provides robust estimations on target building response through scaled ground motions and calculates the dispersion about this target. This alternative procedure is not only useful for record scaling and selection but, upon its further refinement, can also be advantageous for the probabilistic methods that assess the engineering demand parameters for a given target hazard level. Case studies that compare the performance of the proposed procedure with some other record selection and scaling methods suggest its usefulness for building performance assessment and loss models. Copyright © 2012 John Wiley & Sons, Ltd.     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

Ground-Motion Prediction Equations (GMPEs) for inelastic displacement and ductility demands of constant-strength SDOF systems

The objective of this paper is to present ground-motion prediction equations for ductility demand and inelastic spectral displacement of constant-strength perfectly elasto-plastic single-degree-of-freedom (SDOF) oscillators. Empirical equations have been developed to compute the ductility demand as a function of two earthquake parameters; moment magnitude, and source-to-site distance; one site parameter, the ground type; and three oscillator parameters, an undamped natural period, critical damping ratio, and the mass-normalized yield strength. In addition, a comparative study of the proposed model with selected previous studies and recommendations of Eurocode 8 is presented. Proposed equations can easily be incorporated in existing probabilistic seismic hazard analysis (PSHA) software packages with the introduction of an additional parameter. This leads to hazard curves for inelastic spectral displacement, which can provide better estimates of target displacement for nonlinear static procedures and an efficient intensity measure for probabilistic seismic demand analysis (PSDA). Proposed equations will be useful in performance evaluation of existing structures.     

Strong motion uncertainty determined from observed records by dense network in Japan

The variation of ground motions at specific stations from events in six narrow areas was inspected by using K-NET and KiK-net records. A source-area factor for individual observation stations was calculated by averaging ratios between observed values for horizontal peak acceleration and velocity, as well as acceleration response spectra for 5% damping, and predicted values using a ground-motion model (usually known as an attenuation relation) by Kanno et al. (Bull Seismol Soc Am, 96:879–897, 2006). Standard deviations between observed and predicted amplitudes after the correction factor are less than 0.2 on the logarithmic scale and decrease down to around 0.15 in the short-period range. Intra-event standard deviation clearly increases with decreasing distance due to differing paths around near source area. Standard deviations may increase with amplitude or decrease with magnitude; however, both amplitude and magnitude of the data are strongly correlated with distance. The standard deviation calculated in this study is obviously much smaller than that of the original ground-motion model, as epistemic uncertainties are minimized by grouping ground motions at specific stations. This result indicates that the accuracy of strong ground motion prediction could be improved if ground-motion models for specified region are determined individually. For this to be possible, it is necessary to have dense strong-motion networks in high-seismicity regions, such as K-NET and KiK-net.     

A partially non-ergodic ground-motion prediction equation for Europe and the Middle East

A partially non-ergodic ground-motion prediction equation is estimated for Europe and the Middle East. Therefore, a hierarchical model is presented that accounts for regional differences. For this purpose, the scaling of ground-motion intensity measures is assumed to be similar, but not identical in different regions. This is achieved by assuming a hierarchical model, where some coefficients are treated as random variables which are sampled from an underlying global distribution. The coefficients are estimated by Bayesian inference. This allows one to estimate the epistemic uncertainty in the coefficients, and consequently in model predictions, in a rigorous way. The model is estimated based on peak ground acceleration data from nine different European/Middle Eastern regions. There are large differences in the amount of earthquakes and records in the different regions. However, due to the hierarchical nature of the model, regions with only few data points borrow strength from other regions with more data. This makes it possible to estimate a separate set of coefficients for all regions. Different regionalized models are compared, for which different coefficients are assumed to be regionally dependent. Results show that regionalizing the coefficients for magnitude and distance scaling leads to better performance of the models. The models for all regions are physically sound, even if only very few earthquakes comprise one region.     

Treatment of ground-motion predictive relationships for the reference seismic hazard map of Italy

In the framework of the 2004 reference seismic hazard map of Italy the amplitude of the strong-motion (expressed in terms of Peak Horizontal Acceleration with 10% probability of non-exceedence in 50 years, referred to average hard ground conditions) was computed using different predictive relationships. Equations derived in Italy and in Europe from strong-motion data, as well as a set of weak and strong-motion based empirical predictive relationships were employed in a logic tree procedure, in order to capture the epistemic uncertainty affecting ground-motion attenuation. This article describes the adjustments and conversions required to eliminate the incompatibilities amongst the relations. Particularly significant are distance conversions and style-of-faulting adjustments, as well as the problems related to the use of regional relations, such as the selection of a reference depth, the quantification of random variability and the strong-motion prediction. Moreover, a regional attenuation relationship specific for volcanic areas was also employed, allowing a more realistic evaluation of seismic hazard, as confirmed by the attenuation of macroseismic intensities.     

Criteria for Selecting and Adjusting Ground-Motion Models for Specific Target Regions: Application to Central Europe and Rock Sites

A vital component of any seismic hazard analysis is a model for predicting the expected distribution of ground motions at a site due to possible earthquake scenarios. The limited nature of the datasets from which such models are derived gives rise to epistemic uncertainty in both the median estimates and the associated aleatory variability of these predictive equations. In order to capture this epistemic uncertainty in a seismic hazard analysis, more than one ground-motion prediction equation must be used, and the tool that is currently employed to combine multiple models is the logic tree. Candidate ground-motion models for a logic tree should be selected in order to obtain the smallest possible suite of equations that can capture the expected range of possible ground motions in the target region. This is achieved by starting from a comprehensive list of available equations and then applying criteria for rejecting those considered inappropriate in terms of quality, derivation or applicability. Once the final list of candidate models is established, adjustments must be applied to achieve parameter compatibility. Additional adjustments can also be applied to remove the effect of systematic differences between host and target regions. These procedures are applied to select and adjust ground-motion models for the analysis of seismic hazard at rock sites in West Central Europe. This region is chosen for illustrative purposes particularly because it highlights the issue of using ground-motion models derived from small magnitude earthquakes in the analysis of hazard due to much larger events. Some of the pitfalls of extrapolating ground-motion models from small to large magnitude earthquakes in low seismicity regions are discussed for the selected target region.     

Peak ground motion predictions in India: an appraisal for rock sites

Proper selection and ranking of Ground Motion Prediction Equations (GMPEs) is critical for successful logic-tree implementation in probabilistic seismic hazard analysis. The present study explores this issue in predicting peak ground accelerations at the rock sites in India. Macroseismic intensity data complemented with limited strong ground-motion recordings are used for the purpose. The findings corroborate the possible conformity between the GMPEs developed for tectonically active shallow crust across the globe. On the other hand, the relevant GMPEs in the intraplate regions cluster into two different groups with the equations of lower ranks catering to higher ground motions. The earthquakes in the subduction zones have significant regional implications. However, affinity in the ground-motion attenuations between the major interface events (M _W > 7.4) in Andaman-Nicobar, Japan and Cascadia, respectively, is noted. This can be also observed for the intraslab events in the Hindukush and Taiwan respectively. Overall, we do not observe any significant advantage with the equations developed using the regional data. These findings are expected to be useful in probabilistic seismic hazard analysis across the study region.     

Seismic hazard assessment of the Province of Murcia (SE Spain): analysis of source contribution to hazard

A probabilistic seismic hazard assessment of the Province of Murcia in terms of peak ground acceleration (PGA) and spectral accelerations [SA(T)] is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA making use of intensity-to-PGA relationships, hazard is here calculated in terms of magnitude and using European spectral ground-motion models. Moreover, we have considered the most important faults in the region as specific seismic sources, and also comprehensively reviewed the earthquake catalogue. Hazard calculations are performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for three different seismic source zonings and three different ground-motion models. Hazard maps in terms of PGA and SA(0.1, 0.2, 0.5, 1.0 and 2.0 s) and coefficient of variation (COV) for the 475-year return period are shown. Subsequent analysis is focused on three sites of the province, namely, the cities of Murcia, Lorca and Cartagena, which are important industrial and tourism centres. Results at these sites have been analysed to evaluate the influence of the different input options. The most important factor affecting the results is the choice of the attenuation relationship, whereas the influence of the selected seismic source zonings appears strongly site dependant. Finally, we have performed an analysis of source contribution to hazard at each of these cities to provide preliminary guidance in devising specific risk scenarios. We have found that local source zones control the hazard for PGA and SA(T ≤ 1.0 s), although contribution from specific fault sources and long-distance north Algerian sources becomes significant from SA(0.5 s) onwards.     

An Atlas of ShakeMaps and population exposure catalog for earthquake loss modeling

We present an Atlas of ShakeMaps and a catalog of human population exposures to moderate-to-strong ground shaking (EXPO-CAT) for recent historical earthquakes (1973–2007). The common purpose of the Atlas and exposure catalog is to calibrate earthquake loss models to be used in the US Geological Survey’s Prompt Assessment of Global Earthquakes for Response (PAGER). The full ShakeMap Atlas currently comprises over 5,600 earthquakes from January 1973 through December 2007, with almost 500 of these maps constrained—to varying degrees—by instrumental ground motions, macroseismic intensity data, community internet intensity observations, and published earthquake rupture models. The catalog of human exposures is derived using current PAGER methodologies. Exposure to discrete levels of shaking intensity is obtained by correlating Atlas ShakeMaps with a global population database. Combining this population exposure dataset with historical earthquake loss data, such as PAGER-CAT, provides a useful resource for calibrating loss methodologies against a systematically-derived set of ShakeMap hazard outputs. We illustrate two example uses for EXPO-CAT; (1) simple objective ranking of country vulnerability to earthquakes, and; (2) the influence of time-of-day on earthquake mortality. In general, we observe that countries in similar geographic regions with similar construction practices tend to cluster spatially in terms of relative vulnerability. We also find little quantitative evidence to suggest that time-of-day is a significant factor in earthquake mortality. Moreover, earthquake mortality appears to be more systematically linked to the population exposed to severe ground shaking (Modified Mercalli Intensity VIII+). Finally, equipped with the full Atlas of ShakeMaps, we merge each of these maps and find the maximum estimated peak ground acceleration at any grid point in the world for the past 35 years. We subsequently compare this “composite ShakeMap” with existing global hazard models, calculating the spatial area of the existing hazard maps exceeded by the combined ShakeMap ground motions. In general, these analyses suggest that existing global, and regional, hazard maps tend to overestimate hazard. Both the Atlas of ShakeMaps and EXPO-CAT have many potential uses for examining earthquake risk and epidemiology. All of the datasets discussed herein are available for download on the PAGER Web page (). T. I. Allen and M. G. Hearne—contracted through Synergetics Incorporated.     

Probabilistic seismic hazard assessment for Thailand

A set of probabilistic seismic hazard maps for Thailand has been derived using procedures developed for the latest US National Seismic Hazard Maps. In contrast to earlier hazard maps for this region, which are mostly computed using seismic source zone delineations, the presented maps are based on the combination of smoothed gridded seismicity, crustal-fault, and subduction source models. Thailand’s composite earthquake catalogue is revisited and expanded, covering a study area limited by 0°–30°N Latitude and 88°–110°E Longitude and the instrumental period from 1912 to 2007. The long-term slip rates and estimates of earthquake size from paleoseismological studies are incorporated through a crustal fault source model. Furthermore, the subduction source model is used to model the megathrust Sunda subduction zones, with variable characteristics along the strike of the faults. Epistemic uncertainty is taken into consideration by the logic tree framework incorporating basic quantities, such as different source modelling, maximum cut-off magnitudes and ground motion prediction equations. The ground motion hazard map is presented over a 10 km grid in terms of peak ground acceleration and spectral acceleration at 0.2, 1.0, and 2.0 undamped natural periods and a 5% critical damping ratio for 10 and 2% probabilities of exceedance in 50 years. The presented maps give expected ground motions that are based on more extensive data sources than applied in the development of previous maps. The main findings are that northern and western Thailand are subjected to the highest hazard. The largest contributors to short- and long-period ground motion hazard in the Bangkok region are from the nearby active faults and Sunda subduction zones, respectively.     

A new seismic hazard assessment in the region of Andalusia (Southern Spain)

A probabilistic seismic hazard assessment of Andalusia (Southern Spain) in terms of peak ground acceleration, PGA, and spectral accelerations, SA(T), is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA, making use of Intensity-to-PGA relationships, hazard was here calculated in terms of magnitude, using published spectral ground-motion models. Moreover, we considered different ground-motion models for the Atlantic sources, since the attenuation of those motions seems to be slower, as evidenced in the case of the extensive macroseismic areas of earthquakes like those occurred in the years 1755, 1969 and 2007. A comprehensive review of the seismic catalogue and of the seismogenic models proposed for the region was carried out, including those for Northern Africa, which is part of the influence area. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for six different seismic source zonings and five different ground-motion attenuation relationships. Hazard maps in terms of PGA and SA (0.2 s) and SA (1 s) and coefficient of variation (COV) maps, for the 475-year return period were first obtained in rock sites. A geotechnical classification and amplification factors were proposed and new hazard maps including local effects were represented, showing PGA values ranging from 24 to 370 cm/s² for the whole Andalusian territory, with the highest expected values (PGA > 300 cm/s²) in some parts of the Granada Province and in the town of Vélez Málaga. Lowest values (PGA < 50 cm/s²) correspond to some towns of the Huelva and Córdoba provinces. The inclusion of soil effects provides a more detailed picture of the actual hazard the region is subjected to.     

Investigating strong ground-motion variability using analysis of variance and two-way-fit plots

A statistical method to quantitatively assess the relative importance of unmodelled site and source effects on the observed variability (σ) in ground motions is presented. The method consists of analysis of variance (ANOVA) using the computed residuals with respect to an empirical ground-motion model for strong-motion records of various earthquakes recorded at a common set of stations. ANOVA divides the overall variance (σ ²) into the components due to site and source effects (respectively σ _S ² and σ _E ²) not modelled by the ground-motion model plus the residual variance not explained by these effects (σ _R ²). To test this procedure, four sets of observed strong-motion records: two from Italy (Umbria-Marche and Molise), one from the French Antilles and one from Turkey, are used. It is found that for the data from Italy, the vast majority of the observed variance is attributable to unmodelled site effects. In contrast, the variation in ground motions in the French Antilles and Turkey data is largely attributable, especially at short periods, to source effects not modelled by the ground-motion estimation equations used.     

Italian strong motion records in ITACA: overview and record processing

The published version 1.0 of the new Italian strong-motion database ITACA (Italian ACcelerometric Archive, ) includes to date (December 2010) about 4,000 three-component waveforms up to M 6.9, from more than 1,800 earthquakes up to 6.9, recorded by about 400 stations in the period 1972–2009. The uncorrected and corrected strong motion data are archived and can be retrieved with their metadata, concerning events, stations and waveforms. The aim of this paper is to present the procedures for processing the records included in ITACA, accounting for the heterogeneity of this data set, both in terms of quality and amplitude of records as well as illustrating the main features of the ITACA strong-motion dataset. Later, we focus on the “exceptional” ground-motion records, that we, conventionally, denote as those having peak acceleration and peak velocity larger than 300 cm/s² and 15 cm/s, respectively. These records are less than 2% of the whole ITACA dataset but they are the most relevant for the seismic hazard and engineering implications. Such large peak values, recorded at distances up to 30 km, are related not only to the strongest Italian earthquakes, but also to events with magnitude down to 4. Furthermore, we investigate the dependence of the largest peak values on horizontal and vertical directions and on source-to-site distance.