Empirical ground-motion relations using moderate earthquakes recorded by Medellín–Aburrá Valley (Colombia) strong-motion networks

We tested attenuation relations obtained for different regions of the world to verify their suitability to predict strong-motion data recorded by Medellín and Aburrá Valley Accelerographic Networks. We used as comparison criteria, the average of the difference between the observed and the predicted data as a function of epicenter distance and its standard deviation. We also used the approach developed by Sherbaum et al. (Bull Seism Soc Am 94:2164–2185, 2004) that provides a method to evaluate the overall goodness-of-fit of ground-motion prediction equations. The predictive models selected use a generic focal depth. We found that this parameter has an important influence in the ground-motion predictions and must be taken into account as an independent variable. We also found important to characterize the local soil amplification to improve the attenuation relations. We found empirical relations for peak horizontal acceleration PGA and velocity PGV based on the Kamiyama and Yanagisawa (Soils Found 26:16–32, 1986) approach. $$\begin{aligned} \log _{10} (PGA)=0.5886M_L -1.0902\log _{10}(R)-0.0035H+C_{st}\pm 0.\text{29} \end{aligned}$$ $$\begin{aligned} \log _{10} (PGV)=0.7255M_L -1.8812\log _{10}(R)-0.0016H+C_{st}\pm 0.36 \end{aligned}$$ where PGA is measured in cm/s $^{2}$ and PGV in cm/s, $M_{L}$ is local magnitude in the range 2.8–6.5, $R$ is epicentral distance up to 290 km, $H$ is focal depth in km and $C_{st}$ is a coefficient that accounts for the site response due to soil conditions of each recording station. The introduction of focal depth and local site conditions as independent variables, minimize the residuals and the dispersion of the predicted data. We conclude that $H$ and $C_{st}$ are sensitive parameters, having a strong influence on the strong-motion predictions. Using the same functional form, we also propose an empirical relation for the root mean square acceleration a $_\mathrm{rms}$ : $$\begin{aligned} \log _{10} \left( {a_{rms} } \right)=0.4797M_L -1.1665\log _{10} (R)-0.00201H+C_{st}\pm 0.40 \end{aligned}$$ where a $_\mathrm{rms}$ is measured in cm/s $^{2}$ , from the S-wave arrival and using a window length equal to the rupture duration. The other variables are the same as those for PGA and PGV. The site correction coefficients $C_{st}$ found for PGA, PGV and a $_\mathrm{rms}$ show a similar trend indicating a good correlation with the soil conditions of the recording sites.     

Relevant features of the valley seismic response: the case study of Tuscan Northern Apennine sector

The so called “valley effect” relates to the typical seismic response of basin shaped bedrock filled by quaternary sediments. It is an aspect of the renown “local seismic effect” that shall be taken into account when dealing with microzoning studies. Several experimental surveys and numerical simulations performed worldwide over the last 40 years, confirmed that valley responses under seismic excitations show common features in various geological contexts as far as the sedimentary valleys (e.g. alluvial and lacustrine plains), the intermountain valleys (e.g. alpine valleys) and graben shaped basins. Such features mainly depend on the basin geometry, referred to as the shape ratio SR, and the sediment and basin impedance contrast IC. Although researchers agree on the prominent role of local seismic effects for interpreting erratic damages caused by seismic shaking in urbanized areas, no fully shared strategies have been identified for taking into account valley effect within microzoning studies. In this paper, a numerical simulations on three models of trapezoidal shaped basins have been performed. These valley models relate to sediments and basins detected within the Tuscany Region territory during the VEL project. Results, in terms of the amplification index $\text{ F }_{\mathrm{A}}$ F A have been provided. Three “valley effect charts” for various SR and IC values have been propose for taking into account the local seismic effects due to the basin amplifications within microzoning maps.     

A Crustal Model for the Eastern Alps Region and a New Moho Map in Southeastern Europe

In the last two decades, south-central Europe and the Eastern Alps have been widely explored by many seismic refraction experiments (e.g., CELEBRATION 2000, ALP 2002, SUDETES 2003). Although quite detailed images are available along linear profiles, a comprehensive, three-dimensional crustal model of the region is still missing. This limitation makes this region a weak spot in continental-wide comprehensive representations of crustal structure. To improve on this situation, we select and collect 37 published active-source seismic lines in this region. After geo-referencing each line, we sample them along vertical profiles—every 50?km or less along the line—and derive P-wave velocities in a stack of homogeneous layers (separated by discontinuities: depth of crystalline basement, top of lower crust, and Moho). We finally merge the information using geostatistical methods, and infer S-wave velocity and density using empirical scaling relations. We present here the resulting crustal model for a region encompassing the Eastern Alps, Dinarides, Pannonian basin, Western Carpathians and Bohemian Massif, covering the region within $45^{\circ}\text{--}51^{\circ}\hbox{N}$ and $11^{\circ} \text{--} 22^{\circ}\hbox{E}$ with a resolution of $0.2^{\circ} \times 0.2^{\circ}.$ We are also able to extend and update the map of Moho depth in a wider region within $35^{\circ}\text{--}51^{\circ}\hbox{N}$ and $12^{\circ}\text{--}45^{\circ}\hbox{E},$ gathering Moho values from the collected seismic lines, other published dataset and using the European plate reference EPcrust as a background. All the digitized profiles and the resulting model are available online.     

Mantle Transition Zone Structure Beneath Southeastern China and its Implications for Stagnant Slab and Water Transportation in the Mantle

We determined depth variation of the 410- and 660-km discontinuities beneath southeastern China by common-converted-point stacking of \(\rm P\) -wave receiver functions of 121 permanent Chinese seismic stations. We then combined the results with seismic velocity variation to estimate temperature and water content variations in the mantle transition zone of the region. Previous tomographic studies have shown a stagnant slab in the mantle transition zone in eastern Asia that is connected to subduction of the western Pacific. Temperature variations obtained clearly outline the shape of the stagnant slab, with its western edge at 113.5 \(^\circ\) E and the southern edge at 28.5 \(^\circ\) N. The correlation between the location of the stagnant slab and surface tectonics suggests that the Cenozoic extension in eastern China is closely associated with the subduction of the western Pacific and its eastward migration. The water content of the stagnant slab is lower than in surrounding slabs, suggesting that the water has already been released from the subducting slab into the upper mantle.     

2013年甘肃岷县MS6.6地震区域地震波传播模拟与地面运动特征研究

2013年7月22日,甘肃省定西市岷县发生M_S6.6地震,造成了严重的人员伤亡和经济损失.地震发生于青藏高原东北部边缘与甘肃东南部地区的交界处,是该区域一百多年以来发生的最大地震.分析本次地震观测数据显示,岷县M_S6.6地震产生的地震地面运动呈现出明显的区域变化特征.为了研究这种区域变化特征,本文使用全球CRUST1.0和SRTM30模型数据建立了甘肃岷县及周边地区的地下三维传播介质模型,使用并行的三维有限差分方法进行了岷县地震的区域地震波传播模拟,并与研究区内数字地震台网记录的地震观测资料进行对比.通过分析地面运动的峰值速度和持续时间的区域分布特征,发现：在青藏高原东北部边缘,较大的地形起伏对地面运动的峰值速度分布起主要影响作用;在青藏高原外围地区,地形起伏较小,而沉积层主要影响地面运动的峰值速度和持续时间,从而导致了明显的盆地效应,如地面运动的振幅放大和持续时间加长.因此,研究区剧烈的地形起伏和表层沉积层是影响地面运动的重要因素.     

Compatible ground-motion prediction equations for damping scaling factors and vertical-to-horizontal spectral amplitude ratios for the broader Europe region

In a companion article Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4, 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) present a new ground-motion prediction equation (GMPE) for estimating 5 %-damped horizontal pseudo-acceleration spectral (PSA) ordinates for shallow active crustal regions in Europe and the Middle East. This study provides a supplementary viscous damping model to modify 5 %-damped horizontal spectral ordinates of Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) for damping ratios ranging from 1 to 50 %. The paper also presents another damping model for scaling 5 %-damped vertical spectral ordinates that can be estimated from the vertical-to-horizontal (V/H) spectral ratio GMPE that is also developed within the context of this study. For consistency in engineering applications, the horizontal and vertical damping models cover the same damping ratios as noted above. The article concludes by introducing period-dependent correlation coefficients to compute horizontal and vertical conditional mean spectra (Baker in J Struct Eng 137:322–331, 2011). The applicability range of the presented models is the same as of the horizontal GMPE proposed by Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b): as for spectral periods $0.01 \hbox { s}\le \,\hbox {T}\le \,4\hbox { s}$ as well as PGA and PGV for V/H model; and in terms of seismological estimator parameters $4\le \hbox {M}_\mathrm{w} \le 8, \hbox { R} \le 200 \hbox { km}, 150\hbox { m/s}\le \hbox { V}_\mathrm{S30}\le $ 1,200 m/s, for reverse, normal and strike-slip faults. The source-to-site distance measures that can be used in the computations are epicentral $(\hbox {R}_\mathrm{epi})$ , hypocentral $(\hbox {R}_\mathrm{hyp})$ and Joyner–Boore $(\hbox {R}_\mathrm{JB})$ distances. The implementation of the proposed GMPEs will facilitate site-specific adjustments of the spectral amplitudes predicted from probabilistic seismic hazard assessment in Europe and the Middle East region. They can also help expressing the site-specific design ground motion in several formats. The consistency of the proposed models together with the Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) GMPE may be advantageous for future modifications in the ground-motion definition in Eurocode 8 (CEN in Eurocode 8, Design of structures for earthquake resistance—part 1: general rules, seismic actions and rules for buildings. European Standard NF EN 1998-1, Brussels, 2004).     

Local site effect microzonation of Lorca town (SE Spain)

Local site effect assessment based on subsurface ground conditions is often the key to evaluate urban seismic hazard. The site effect evaluation in Lorca town (south-eastern Spain) started with a classification of urban geology through the geological mapping at scale 1:10,000 and the use of geotechnical data and geophysical surveys. The 17 geological formations identified were classified into 5 geological/seismic formations according to their seismic amplification capacity obtained from ambient vibration measurements as well as from simultaneous strong motion records. The shear-wave velocity structure of each geological/seismic formation was evaluated by means of inversion of Rayleigh wave dispersion data obtained from vertical-component array records of ambient noise. Nakamura’s method was applied to determine a predominant period distribution map. The spectral amplification factors were fourfold the values recorded in a reference hard-rock site. Finally, the capability of this study for explaining the damage distribution caused by the May 11th, 2011 Lorca destructive earthquake (Mw \(=\) 5.2) was examined. The methods used in this work are of assistance to evaluate ground amplification phenomena in urban areas of complex geology as Lorca town due to future earthquakes with applicability on urban seismic risk management.     

Gradual Fault Weakening with Seismic Slip: Inferences from the Seismic Sequences of L’Aquila, 2009, and Northridge, 1994

We estimate seismological fracture energies from two subsets of events selected from the seismic sequences of L’Aquila (2009), and Northridge (1994): 57 and 16 selected events, respectively, including the main shocks. Following Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), we postulate that fracture energy (G) represents the post-failure integral of the dynamic weakening curve, which is described by the evolution of shear traction as a function of slip. Following a direct-wave approach, we compute mainshock-/aftershock-source spectral ratios, and analyze them using the approach proposed by Malagnini et al. (Pure Appl. Geophys., this issue, 2014) to infer corner frequencies and seismic moment. Our estimates of source parameters (including fracture energies) are based on best-fit grid-searches performed over empirical source spectral ratios. We quantify the source scaling of spectra from small and large earthquakes by using the MDAC formulation of Walter and Taylor (A revised Magnitude and Distance Amplitude Correction (MDAC2) procedure for regional seismic discriminants, 2001). The source parameters presented in this paper must be considered as point-source estimates representing averages calculated over specific ruptured portions of the fault area. In order to constrain the scaling of fracture energy with coseismic slip, we investigate two different slip-weakening functions to model the shear traction as a function of slip: (i) a power law, as suggested by Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), and (ii) an exponential decay. Our results show that the exponential decay of stress on the fault allows a good fit between measured and predicted fracture energies, both for the main events and for their aftershocks, regardless of the significant differences in the energy budgets between the large (main) and small earthquakes (aftershocks). Using the power-law slip-weakening function would lead us to a very different situation: in our two investigated sequences, if the aftershock scaling is extrapolated to events with large slips, a power law (a la Abercrombie and Rice) would predict unrealistically large stress drops for large, main earthquakes. We conclude that the exponential stress evolution law has the advantage of avoiding unrealistic stress drops and unbounded fracture energies at large slip values, while still describing the abrupt shear-stress degradation observed in high-velocity laboratory experiments (e.g., Di Toro et al., Fault lubrication during earthquakes, Nature 2011).     

Exploration of S-wave velocity profiles at strong motion stations in Eskisehir,Turkey, using microtremor phase velocity and S-wave amplification

We have explored 1D S-wave velocity profiles of shallow and deep soil layers over a basement at strong motion stations in Eskisehir Province, Turkey. Microtremor array explorations were conducted at eight strong motion stations in the area to know shallow 1D S-wave velocity models. Rayleigh wave phase velocity at a frequency range from 3 to 30 Hz was estimated with the spatial autocorrelation analysis of array records of vertical microtremors at each station. Individual phase velocity was inverted to a shallow S-wave velocity profile. Low-velocity layers were identified at the stations in the basin. Site amplification factors from S-wave parts of earthquake records that had been estimated at the strong motion stations by Yamanaka et al. (2017) were inverted to the S-wave velocities and Q-values of the sedimentary layers. The depths to the basement with an S-wave velocity of 2.2 km/s are about 1 km in the central part of the basin, while the basement becomes shallow as 0.3 km in the marginal part of the basin. We finally discussed the effects of the shallow and deep sedimentary layers on the 1D S-wave amplification characteristics using the revealed profiles. It is found that the shallow soil layers have no significant effects in the amplification at a frequency range lower than 3 Hz in the area.     

A Numerical Experiment on the Horizontal to Vertical Spectral Ratio in Flat Sedimentary Basins

—?In this paper we study the seismic response of flat sedimentary basins and carry out numerical experiments to determine the extent to which we could go using the Horizontal to Vertical Spectral Ratio (HVSR) for a given site. The HVSR has been used by many researchers to characterize local conditions in terms of the dynamic response of the soil, and one of its variants, that proposed by NAKAMURA (1989) in which records of microtremors are used, is one of the most applied in recent years. We study the response of different configurations under incident waves coming from an explosive source using the Indirect Boundary Element Method (IBEM), and we investigate two cases: low- and high-velocity contrast. We compute the seismic response using the HVSR technique at various locations in the free surface of the basins, and compare it with the response calculated with the horizontal Sediment to Bedrock Spectral Ratio (SBSR) and with the Horizontal Component (HC) of the transfer function for the displacement at the same locations. The comparison shows that, in general, HVSR cannot provide the predominant period of a site due to the fact that this technique cannot predict accurately the Spectral amplification levels. On the other hand, the HVSR provides an erroneous response in the sedimentary basins which have a low-impedance contrast, with respect to bedrock, and with shape ratios like the one studied here, whereas it can reasonably well predict the fundamental local frequency when there is a high-impedance contrast, except in the center of the basin.     

Near real-time automatic moment magnitude estimation

In this paper we describe a stable automatic method to estimate in real time the seismic moment, moment magnitude and corner frequency of events recorded by a network comprising broad-band and accelerometer sensors. The procedure produces reliable results even for small-magnitude events $\hbox {M}_{\mathrm{W}}\approx 3$ . The real-time data arise from both the Transfrontier network at the Alps-Dinarides junction and from the Italian National Accelerometric Network (RAN). The data is pre-processed and the S-wave train identified through the application of an automatic method, which estimates the arrival times based on the hypocenter location, recording site and regional velocity model. The transverse component of motion is used to minimize conversion effects. The source spectrum is obtained by correcting the signals for geometrical spreading and intrinsic attenuation. Source spectra for both velocity and displacement are computed and, following Andrews (1986), the seismic moment and the first estimate of the corner frequency, $f_{0}$ , derived. The procedure is validated using the recordings of some recent moderate earthquakes (Carnia 2002; Bovec 2004; Parma 2008; Aquila 2009; Macerata 2009; Emilia 2012) and the recordings of some minor events in the SE Alps area for which independent seismic moment and moment magnitude estimates are available. The results obtained with a dataset of 843 events recorded by the Transfrontier and RAN networks show that the procedure is reliable and robust for events with $\hbox {M}_{\mathrm{W}}\ge 3$ . The estimates of $f_{0}$ are less reliable. The results show a scatter, principally for small events with $\hbox {M}_{\mathrm{W}}\le 3$ , probably due to site effects and inaccurate locations.     

New code site classification, amplification factors and normalized response spectra based on a worldwide ground-motion database

A global dataset of more than 3,000 ground motion records from 536 sites from Greece, Italy, Turkey, USA and Japan is compiled and used to propose code-oriented elastic acceleration response spectra and soil amplification factors for a new site classification system, which, besides the classical geotechnical parameters $N_{SPT}, S_{u}$ and PI, uses also the fundamental period of the site, the thickness of soil deposits and the average shear wave velocity to the seismic bedrock, instead of $V_{s,30}$ . We propose a new classification system with the associated amplification factors and normalized response spectra for two seismicity levels, i.e. $M_{s}\le 5.5$ and $M_{s}>5.5$ . Uncertainties in the estimation of soil amplification factors are captured using a logic-tree approach, which allows the efficient use of alternative models and methods. The aim of this work is to improve the present EC8 soil classification. The effectiveness of the proposed classification system is compared to that of EC8 classification system using an error term, which represents the average dispersion of data within all categories of a given classification scheme. Error terms for the new classification system are lower than the error terms for EC8 classification system at all periods.     

An Evaluation of Earthquake Hazard Potential for Different Regions in Western Anatolia Using the Historical and Instrumental Earthquake Data

We applied the maximum likelihood method produced by Kijko and Sellevoll (Bull Seismol Soc Am 79:645–654, 1989; Bull Seismol Soc Am 82:120–134, 1992) to study the spatial distributions of seismicity and earthquake hazard parameters for the different regions in western Anatolia (WA). Since the historical earthquake data are very important for examining regional earthquake hazard parameters, a procedure that allows the use of either historical or instrumental data, or even a combination of the two has been applied in this study. By using this method, we estimated the earthquake hazard parameters, which include the maximum regional magnitude $ \hat{M}_{\max } , $ the activity rate of seismic events and the well-known $ \hat{b} $ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. The whole examined area is divided into 15 different seismic regions based on their tectonic and seismotectonic regimes. The probabilities, return periods of earthquakes with a magnitude M?≥?m and the relative earthquake hazard level (defined as the index K) are also evaluated for each seismic region. Each of the computed earthquake hazard parameters is mapped on the different seismic regions to represent regional variation of these parameters. Furthermore, the investigated regions are classified into different seismic hazard level groups considering the K index. According to these maps and the classification of seismic hazard, the most seismically active regions in WA are 1, 8, 10 and 12 related to the Alia?a Fault and the Büyük Menderes Graben, Aegean Arc and Aegean Islands.     

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.     

New Seismic Noise Models Obtained Using Very Broadband Stations

It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.     

On the seismic response of under-designed caisson foundations

The seismic behaviour of caisson foundations supporting typical bridge piers is analysed with 3D finite elements, with due consideration to soil and interface nonlinearities. Single-degree-of freedom oscillators of varying mass and height, simulating heavily and lightly loaded bridge piers, founded on similar caissons are studied. Four different combinations of the static ( $\text{ FS }_\mathrm{V}$ FS V ) and seismic ( $\text{ FS }_\mathrm{E}$ FS E ) factors of safety are examined: (1) a lightly loaded ( $\text{ FS }_\mathrm{V}= 5$ FS V = 5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson, (2) a lightly loaded seismically over-designed ( $\text{ FS }_\mathrm{E} >1$ FS E > 1 ) caisson, (3) a heavily loaded ( $\text{ FS }_\mathrm{V} = 2.5$ FS V = 2.5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson and (4) a heavily loaded seismically over-designed caisson. The analysis is performed with use of seismic records appropriately modified so that the effective response periods (due to soil-structure-interaction effects) of the studied systems correspond to the same spectral acceleration, thus allowing their inelastic seismic performance to be compared on a fair basis. Key performance measures of the systems are then contrasted, such as: accelerations, displacements, rotations and settlements. It is shown that the performance of the lightly loaded seismically under-designed caisson is advantageous: not only does it reduce significantly the seismic load to the superstructure, but it also produces minimal residual displacements of the foundation. For heavily loaded foundations, however ( $\text{ FS }_{V} = 2.5$ FS V = 2.5 ), the performance of the two systems (over and under designed) is similar.     

Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East

This article presents the latest generation of ground-motion models for the prediction of elastic response (pseudo-) spectral accelerations, as well as peak ground acceleration and velocity, derived using pan-European databases. The models present a number of novelties with respect to previous generations of models (Ambraseys et al. in Earthq Eng Struct Dyn 25:371–400, 1996, Bull Earthq Eng 3:1–53, 2005; Bommer et al. in Bull Earthq Eng 1:171–203, 2003; Akkar and Bommer in Seismol Res Lett 81:195–206, 2010), namely: inclusion of a nonlinear site amplification function that is a function of $\text{ V }_\mathrm{S30}$ and reference peak ground acceleration on rock; extension of the magnitude range of applicability of the model down to $\text{ M }_\mathrm{w}$ 4; extension of the distance range of applicability out to 200 km; extension to shorter and longer periods (down to 0.01 s and up to 4 s); and consistent models for both point-source (epicentral, $\text{ R }_\mathrm{epi}$ , and hypocentral distance, $\text{ R }_\mathrm{hyp}$ ) and finite-fault (distance to the surface projection of the rupture, $\text{ R }_\mathrm{JB}$ ) distance metrics. In addition, data from more than 1.5 times as many earthquakes, compared to previous pan-European models, have been used, leading to regressions based on approximately twice as many records in total. The metadata of these records have been carefully compiled and reappraised in recent European projects. These improvements lead to more robust ground-motion prediction equations than have previously been published for shallow (focal depths less than 30 km) crustal earthquakes in Europe and the Middle East. We conclude with suggestions for the application of the equations to seismic hazard assessments in Europe and the Middle East within a logic-tree framework to capture epistemic uncertainty.     

Recent evolution and challenges in the Seismic Hazard Analysis of the Po Plain region, Northern Italy

SH evaluation in Po Plain region is revisited with a view to: ( $a)$ updating available seismic source and ground motion attenuation models with recent data and testing the sensitivity of the results, ( $b)$ exploiting available non-stationary models to check their predictive performance at particular sites and, ( $c)$ handling the ground motion hazard generated by composite fault systems. Update of area source models and GMPEs does not strongly affect SH evaluated at representative sites, while usefulness of predictions derived from non-stationary models needs to be assessed on a site-by-site basis. A hybrid (deterministic-probabilistic) approach seems to be a promising tool for assessing site-specific hazard contributed by composite sources.     

Modeling Seismic Attributes of Pn Waves using the Spectral-Element Method

To investigate the nature of Pn propagation, we have implemented the spectral-element method (SEM) for vertically and laterally varying media with and without attenuation. As a practical measure, essential features of the Pn waves are distilled into seismic attributes including arrival times, amplitudes and pulse frequencies. To validate the SEM simulations, we first compare the SEM results with reflectivity calculations of Braile and Smith (Geophys. J.R. Astr. Soc. 40, 145–176, 1975) and then to the asymptotic results of ?erveny and Ravindra (Theory of Seismic Headwaves, University of Toronto Press, pp. 235–250, 1971). Models with random, laterally varying Moho structures are then simulated, where the amplitude and pulse frequency characteristics are found to be stable to small Moho interface perturbations. SEM calculations for models with different upper-mantle velocity gradients are next performed where it is found that interference effects can strongly influence the Pn amplitudes and pulses frequencies. For larger-scale, laterally varying structures, SEM models similar to that found along the Hi-CLIMB array in Tibet are then performed. It is observed that large-scale structures, along with small-scale structures, upper-mantle velocity gradients and attenuation, can all significantly affect the Pn attributes. Ambiguities between upper-mantle velocity gradients and attenuation are also found when using Pn amplitudes and pulse frequency attributes. These ambiguities may be resolved, to some degree, by using the curvature of the travel times at longer regional distance, however, this would also be complicated by lateral variability.     

Implications of the 2011 M9.0 Tohoku Japan earthquake for the treatment of site effects in large earthquakes

Site response in Japan is characterized using thousands of surface and borehole recordings from events of moment magnitude $(\mathbf{M}) > 5.5$ collected by the KiK-net network, including the 2011 M9.0 Tohoku earthquake. Site amplification is defined by the ratio of motions at the surface to those at depth (within the borehole), corrected for the depth effect due to destructive interference using a technique based on cross-spectral ratios between surface and down-hole motions. Site effects were particularly strong at high frequencies, despite the expectation that high-frequency response may be damped by nonlinear effects. In part, the large amplitudes at high frequencies are due to the prevalence of shallow soil conditions in Japan. We searched for typical symptoms for soil nonlinearity, such as a decrease in the predominant frequency and/or amplification, using spectral ratios of weak to strong ground motions. Localized nonlinearity occurred at some recording sites, but was not pervasive. We developed a general empirical model to express site amplification for the KiK-net sites as a function of common site variables, such as the average shear-wave velocity in the uppermost 30 m ( $\text{ V}_\mathrm{S30})$ and the horizontal-to-vertical (H/V) spectral ratio. We use the model to estimate site-corrected ground-motions for the Tohoku mainshock for a reference site condition; these motions are in reasonable agreement with the predictions of some of the published ground motion prediction equations for subduction zones.