Site effects at Euroseistest—I. Determination of the valley structure and confrontation of observations with 1D analysis

This paper describes the process of construction of the 2D model of Volvi's geological structure and results of empirical and theoretical approaches to the evaluation of site response at Euroseistest. The construction of the 2D model is based on a re-interpretation of the available geophysical and geotechnical data in an effort to improve the definition of the subsoil structure at Euroseistest in terms of the most important parameters needed to model site response. The results of this re-interpretation are compared with a previous published 2D model of the same alluvial valley. Different analysis of the measurements and different criteria in the synthesis of data have led to a different model, even if both studies had access to the same field measurements. This underscores the fact that a model results of an interpretation and is not uniquely determined by the data, no matter how detailed they are. The well known subsoil structure opened the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from a detailed analysis of two events in frequency and time domains and 1D numerical modeling. The study of site response shows the important effect of the lateral variations on the ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In the case of Volvi's graben, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects in this basin.     

Complex Site Effects in Thessaloniki (Greece): II. 2D SH Modelling and Engineering Insights

This paper presents results of numerical modelling of site response for Thessaloniki, obtained with two different 2D methods; a finite difference and a finite element method. Ground motion across a 2D model of the subsoil of the city has been simulated for vertically incident SH waves. The predominance of locally generated surface waves is very clear in the synthetic seismograms of a weak event and of stronger ones. These results are then compared with the observations in time domain and frequency domain. The role of the soil formations with high attenuation in the lateral propagation and the effect of the differential motion close to the lateral variations are also pinpointed. The stronger events were finally used to compute strong ground motion in order to reveal and to discuss practical engineering aspects such as peak ground acceleration value, the most familiar indicator in seismic norms, the soil to rock spectral coefficients for the period bandwidth of interest, and the aggravation factor in terms of 2D to 1D response spectra as a useful ruler to account for complex site effects.     

Parametric analysis of the seismic response of a 2D sedimentary valley: implications for code implementations of complex site effects

A detailed 2D model has been constructed and validated for Euroseistest valley, in northern Greece. We take advantage of this model to investigate what parameters, in addition to surface soil conditions (obviously the most important parameter), can be used to correctly characterize site response in a 2D structure. Through a parametric analysis using 2D numerical simulations for SH waves, we explore the differences between the computed ground motion for different simplifications of the valley's structure. We consider variations in the velocity structure within the sediments, and variations of the shape between sediments and bedrock. We also compare the results from different 1D models reflecting current approaches to the determination of site response. Our results show clearly that, in the case of Euroseistest, site response owes fundamentally to its closed basin shape because it is largely controlled by locally generated surface waves. Thus, in terms of predicting site response, a rough idea of its shape ratio and of the average mechanical properties of the sediments are better than a very detailed 1D profile at the central site. Although the details of ground motion may vary significantly between the models, the relative amount of surface waves generated in the 2D models seems to be relatively constant. Moreover, if we quantify the additional amplification caused by the lateral heterogeneity in terms of the ‘aggravation factor’ introduced by Chávez-García & Faccioli [7], a roughly constant factor between 2 and 3 seems to appropriately take into account the effects of lateral heterogeneity. Of course, a correct estimate of the overall impedance contrast is necessary to correctly predict the maximum amplification, a caveat that also applies to 1D models. In this sense, Euroseistest rings an alarm bell. In this valley the more significant impedance contrast lies at about 200 m depth, and it is missed both by consideration of the average shear wave velocity of the first 30 m (the V_s30 criterion) or using the detailed velocity profile down to a depth where a shear wave velocity larger than 750 m/s is found. Our conclusions indicate that, in order to improve current schemes to take into account site effects in building codes, the more to be gained comes from consideration of lateral heterogeneity, at least in the case of shallow alluvial valleys, where locally generated surface waves are likely to be important.     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Site effects of the 2002 Molise earthquake, Italy: analysis of strong motion, ambient noise, and synthetic data from 2D modelling in San Giuliano di Puglia

On 31 October and 1 November 2002, the Basso Molise area (Southern Italy) was struck by two earthquakes of moderate magnitude (M _L = 5.4 and 5.3). The epicentral area showed a high level of damage, attributable both to the high vulnerability of existing buildings and to site effects caused by the geological and geomorphological settings. Specifically, the intensity inside the town of San Giuliano di Puglia was two degrees higher than in neighbouring towns. Also, within San Giuliano di Puglia, the damage varied notably. The site response in the city was initially evaluated from horizontal-to-vertical spectral ratios (HVSR) from a limited number of strong motion recordings of the most severe aftershocks. Several microtremor measurements were also available. Both data sets indicated the simultaneous presence of two amplification peaks: one around 6 Hz, attributed in previous studies to the strong, shallow impedance contrast among landfill/clay and calcarenites, and one at 2 Hz related to the first S-wave arrivals and predominantly seen only on one receiver component. Further studies performed on weak-motion recordings also showed strong amplification on the vertical receiver component, thus indicating an underestimation of the amplification by the HVSR technique. Additionally, a 2D-model of the geology of the sub-surface was developed, reproducing the flower-shaped structure generated during the late orogenic transpressive regime. The numerical (finite-difference hybrid) simulation reproduced the two peaks of the observed data at slightly higher frequencies. The model also confirmed that the borders of the flower structure define a boundary between amplification levels, with higher amplification inside.     

Seismic site effects by an optimized 2D BE/FE method II. Quantification of site effects in two-dimensional sedimentary valleys

This paper deals with the evaluation of seismic site effects due to the local topographical and geotechnical characteristics. The amplification of surface motions is calculated by a numerical method combining finite elements in the near field and boundary elements in the far field (FEM/BEM). The numerical technique is improved by time truncation. In the first part of this article, the accuracy and the relevance of this optimized method are presented. Moreover, parametric studies are done on slopes, ridges and canyons to characterize topographical site effects. The second part deals with sedimentary valleys. The complexity of the combination of geometrical and sedimentary effects is underlined. Extensive parametrical studies are done to discriminate the topographical and geotechnical effects on seismic ground movement amplifications in two-dimensional irregular configurations. Characteristic coefficients are defined to predict the amplifications of horizontal displacements. The accuracy of this quantitative evaluation technique is tested and discussed.     

Site response in Tecoman,Colima, Mexico—II. Determination of subsoil structure and comparison with observations

In a companion paper local transfer functions were estimated at Tecoman using earthquake and microtremor data. In this paper, the subsoil structure at this city is investigated using seismic refraction and cross-correlation of noise records as a case study. P- and S-wave refraction profiles were measured at five sites within the city. Standard analysis constrained only very shallow layers. The P-wave refraction deployment was also used to record ambient vibration. These data were processed using an extension of the SPAC (SPatial AutoCorrelation (Aki, 1957) [1]) method; cross-correlation is computed between station pairs and the results are inverted to obtain a phase velocity dispersion curve. Penetration depth was larger than that from the refraction experiments but the shear-wave velocity of the basement could not be determined. For this reason, additional microtremor measurements were made using broad band seismometers with a larger spacing between stations. The results allowed to constrain the shear-wave velocity of the basement. Site amplification computed for the final profiles compare well with observed ground motion amplification at Tecoman. The case of Tecoman illustrates that even a simple subsoil structure may require crossing data from different experiments to correctly constrain site effects.     

Complex Site Effects in Thessaloniki (Greece): I. Soil Structure and Comparison of Observations with 1D Analysis

This paper presents a 2D model of the geological structure of Thessaloniki city and results of empirical and theoretical approaches for the evaluation of site response due to complex site effects. The construction of the 2D model is based on the available geophysical and geotechnical data in terms of the most important parameters needed to model site response. The well-known subsoil structure, despite the existence of some local uncertainties, gave the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from the analysis of one event in frequency and time domains and 1D numerical modelling. The study of site response shows the effect of the lateral variations on ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In this case history, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects.     

Seismic site effects by an optimized 2D BE/FE method I. Theory, numerical optimization and application to topographical irregularities

This paper deals with the evaluation of seismic site effects due to the local topographical and geotechnical characteristics. The amplification of surface motions is calculated by a numerical method combining finite elements in the near field and boundary elements in the far field (FEM/BEM). The numerical technique is improved by time truncation. In the first part of this article, the accuracy and the relevance of this optimized method are presented. Moreover, parametric studies are done on slopes, ridges and canyons to characterize topographical site effects. The second part deals with sedimentary valleys. The complexity of the combination of geometrical and sedimentary effects is underlined. Extensive parametrical studies are done to discriminate the topographical and geotechnical effects on seismic ground movement amplifications in two-dimensional irregular configurations. Characteristic coefficients are defined to predict the amplifications of horizontal displacements. The accuracy of this quantitative evaluation technique is tested and discussed.