@article{PilzParolaiLeytonetal.2009, author = {Pilz, Marco and Parolai, Stefano and Leyton, Felipe and Campos, Jaime and Zschau, Jochen}, title = {A comparison of site response techniques using earthquake data and ambient seismic noise analysis in the large urban areas of Santiago de Chile}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2009.04195.x}, year = {2009}, abstract = {Situated in an active tectonic region, Santiago de Chile, the country's capital with more than six million inhabitants, faces tremendous earthquake risk. Macroseismic data for the 1985 Valparaiso event show large variations in the distribution of damage to buildings within short distances, indicating strong effects of local sediments on ground motion. Therefore, a temporary seismic network was installed in the urban area for recording earthquake activity and a study was carried out aiming to estimate site amplification derived from horizontal-to- vertical (H/V) spectral ratios from earthquake data (EHV) and ambient noise (NHV), as well as using the standard spectral ratio (SSR) technique with a nearby reference station located on igneous rock. The results lead to the following conclusions: The analysis of earthquake data shows significant dependence on the local geological structure with respect to amplitude and duration. An amplification of ground motion at frequencies higher than the fundamental one can be found. This amplification would not be found when looking at NHV ratios alone. The analysis of NHV spectral ratios shows that they can only provide a lower bound in amplitude for site amplification. P-wave site responses always show lower amplitudes than those derived by S waves, and sometimes even fail to provide some frequencies of amplification. No variability in terms of time and amplitude is observed in the analysis of the H/V ratio of noise. Due to the geological conditions in some parts of the investigated area, the fundamental resonance frequency of a site is difficult to estimate following standard criteria proposed by the SESAME consortium, suggesting that these are too restrictive under certain circumstances.}, language = {en} } @phdthesis{Pilz2010, author = {Pilz, Marco}, title = {A comparison of proxies for seismic site conditions and amplification for the large urban area of Santiago de Chile}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-52961}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Situated in an active tectonic region, Santiago de Chile, the country´s capital with more than six million inhabitants, faces tremendous earthquake hazard. Macroseismic data for the 1985 Valparaiso and the 2010 Maule events show large variations in the distribution of damage to buildings within short distances indicating strong influence of local sediments and the shape of the sediment-bedrock interface on ground motion. Therefore, a temporary seismic network was installed in the urban area for recording earthquake activity, and a study was carried out aiming to estimate site amplification derived from earthquake data and ambient noise. The analysis of earthquake data shows significant dependence on the local geological structure with regards to amplitude and duration. Moreover, the analysis of noise spectral ratios shows that they can provide a lower bound in amplitude for site amplification and, since no variability in terms of time and amplitude is observed, that it is possible to map the fundamental resonance frequency of the soil for a 26 km x 12 km area in the northern part of the Santiago de Chile basin. By inverting the noise spectral rations, local shear wave velocity profiles could be derived under the constraint of the thickness of the sedimentary cover which had previously been determined by gravimetric measurements. The resulting 3D model was derived by interpolation between the single shear wave velocity profiles and shows locally good agreement with the few existing velocity profile data, but allows the entire area, as well as deeper parts of the basin, to be represented in greater detail. The wealth of available data allowed further to check if any correlation between the shear wave velocity in the uppermost 30 m (vs30) and the slope of topography, a new technique recently proposed by Wald and Allen (2007), exists on a local scale. While one lithology might provide a greater scatter in the velocity values for the investigated area, almost no correlation between topographic gradient and calculated vs30 exists, whereas a better link is found between vs30 and the local geology. When comparing the vs30 distribution with the MSK intensities for the 1985 Valparaiso event it becomes clear that high intensities are found where the expected vs30 values are low and over a thick sedimentary cover. Although this evidence cannot be generalized for all possible earthquakes, it indicates the influence of site effects modifying the ground motion when earthquakes occur well outside of the Santiago basin. Using the attained knowledge on the basin characteristics, simulations of strong ground motion within the Santiago Metropolitan area were carried out by means of the spectral element technique. The simulation of a regional event, which has also been recorded by a dense network installed in the city of Santiago for recording aftershock activity following the 27 February 2010 Maule earthquake, shows that the model is capable to realistically calculate ground motion in terms of amplitude, duration, and frequency and, moreover, that the surface topography and the shape of the sediment bedrock interface strongly modify ground motion in the Santiago basin. An examination on the dependency of ground motion on the hypocenter location for a hypothetical event occurring along the active San Ram{\´o}n fault, which is crossing the eastern outskirts of the city, shows that the unfavorable interaction between fault rupture, radiation mechanism, and complex geological conditions in the near-field may give rise to large values of peak ground velocity and therefore considerably increase the level of seismic risk for Santiago de Chile.}, language = {en} } @article{PilzParolaiStupazzinietal.2011, author = {Pilz, Marco and Parolai, Stefano and Stupazzini, Marco and Paolucci, Roberto and Zschau, Jochen}, title = {Modelling basin effects on earthquake ground motion in the Santiago de Chile basin by a spectral element code}, series = {Geophysical journal international}, volume = {187}, journal = {Geophysical journal international}, number = {2}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2011.05183.x}, pages = {929 -- 945}, year = {2011}, abstract = {Simulations of strong ground motion within the Santiago de Chile Metropolitan area were carried out by means of 3-D deterministic wave propagation tool based on the spectral element method. The simulated events take into account the pronounced interface between the low-velocity sedimentary basin and the bedrock as well as topography of the area. To verify our model we simulated a regional earthquake recorded by a dense network installed in the city of Santiago for recording aftershock activity after the 2010 February 27 Maule main shock. The results proof the alluvial basin amplification effects and show a strong dependence of spectral amplification in the basin on the local site conditions. Moreover, we studied the seismic response due to a hypothetical M(w) = 6.0 event occurring along the active San Ramon Fault, which is crossing the eastern edge of the city. The scenario earthquakes exhibit that an unfavourable interaction between fault rupture, radiation mechanism and complex geological and topographic conditions in the near-field region may give rise to large values of peak ground velocity in the basin. Finally, 3-D numerical predictions of ground motion are compared with the one computed according to ground motion prediction equations selected among the next generation attenuation relationships, in terms of ground motion peak values and spectral acceleration. The comparison underlines that the 3-D scenario simulations predict a significantly higher level of ground motion in the Santiago basin, especially over deep alluvial deposits. Moreover, also the location of the rupture nucleation largely influences the observed shaking pattern.}, language = {en} } @article{PilzParolaiPicozzietal.2011, author = {Pilz, Marco and Parolai, Stefano and Picozzi, Matteo and Zschau, Jochen}, title = {Evaluation of proxies for seismic site conditions in large urban areas the example of Santiago de Chile}, series = {Physics and chemistry of the earth}, volume = {36}, journal = {Physics and chemistry of the earth}, number = {16}, publisher = {Elsevier}, address = {Oxford}, issn = {1474-7065}, doi = {10.1016/j.pce.2011.01.007}, pages = {1259 -- 1266}, year = {2011}, abstract = {Characterizing the local site response in large cities is an important step towards seismic hazard assessment. To this regard, single station seismic noise measurements were carried out at 146 sites in the northern part of Santiago de Chile. This extensive survey allowed the fundamental resonance frequency of the sedimentary cover, derived from horizontal-to-vertical (H/V) spectral ratios, to be mapped. By inverting the spectral ratios under the constraint of the thickness of the sedimentary cover, known from previous gravimetric measurements, local S-wave velocity profiles have been retrieved. After interpolation between the individual profiles, the resulting high resolution 3D S-wave velocity model allows the entire area, as well as deeper parts of the basin, to be represented in great detail. Since one lithology shows a great scatter in the velocity values only a very general correlation between S-wave velocity in the uppermost 30 m (v(s)(30)) and local geology is found. Local S-wave velocity profiles can serve as a key factor in seismic hazard assessment, since they allow an estimate of the amplification potential of the sedimentary cover. Mapping the intensity distribution of the 27 February 2010 Maule, Chile, event (Mw = 8.8) the results indicate that local amplification of the ground motion might partially explain the damage distribution and encourage the use of the low cost seismic noise techniques for the study of seismic site effects.}, language = {en} } @article{PilzCotton2019, author = {Pilz, Marco and Cotton, Fabrice Pierre}, title = {Does the One-Dimensional Assumption Hold for Site Response Analysis?}, series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, volume = {35}, journal = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, number = {2}, publisher = {Earthquake Engineering Research Institute}, address = {Oakland}, issn = {8755-2930}, doi = {10.1193/050718EQS113M}, pages = {883 -- 905}, year = {2019}, abstract = {The one-dimensional (1-D) approach is still the dominant method to incorporate site effects in engineering applications. To bridge the 1-D to multidimensional site response analysis, we develop quantitative criteria and a reproducible method to identify KiK-net sites with significant deviations from 1-D behavior. We found that 158 out of 354 show two-dimensional (2-D) and three-dimensional (3-D) effects, extending the resonance toward shorter periods at which 2-D or 3-D site effects exceed those of the classic 1-D configurations and imposing an additional amplification to that caused by the impedance contrast alone. Such 2-D and 3-D effects go along with a large within-station ground motion variability. Remarkably, these effects are found to be more pronounced for small impedance contrasts. While it is hardly possible to identify common features in ground motion behavior for stations with similar topography typologies, it is not over-conservative to apply a safety factor to account for 2-D and 3-D site effects in ground motion modeling.}, language = {en} } @article{PilzCottonZaccarellietal.2019, author = {Pilz, Marco and Cotton, Fabrice Pierre and Zaccarelli, Riccardo and Bindi, Dino}, title = {Capturing Regional Variations of Hard-Rock Attenuation in Europe}, series = {Bulletin of the Seismological Society of America}, volume = {109}, journal = {Bulletin of the Seismological Society of America}, number = {4}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120190023}, pages = {1401 -- 1418}, year = {2019}, abstract = {A proper assessment of seismic reference site conditions has important applications as they represent the basis on which ground motions and amplifications are generally computed. Besides accounting for the average S-wave velocity over the uppermost 30 m (V-S30), the parameterization of high-frequency ground motions beyond source-corner frequency received significant attention. kappa, an empirical parameter introduced by Anderson and Hough (1984), is often used to represent the spectral decay of the acceleration spectrum at high frequencies. The lack of hard-rock records and the poor understanding of the physics of kappa introduced significant epistemic uncertainty in the final seismic hazard of recent projects. Thus, determining precise and accurate regional hard-rock kappa(0) values is critical. We propose an alternative procedure for capturing the reference kappa(0) on regional scales by linking thewell-known high-frequency attenuation parameter kappa and the properties of multiple-scattered coda waves. Using geological and geophysical data around more than 1300 stations for separating reference and soft soil sites and based on more than 10,000 crustal earthquake recordings, we observe that kappa(0) from multiple-scattered coda waves seems to be independent of the soil type but correlated with the hard-rock kappa(0), showing significant regional variations across Europe. The values range between 0.004 s for northern Europe and 0.020 s for the southern and southeastern parts. On the other hand, measuring kappa (and correspondingly kappa(0)) on the S-wave window (as classically proposed), the results are strongly affected by transmitted (reflected, refracted, and scattered) waves included in the analyzed window biasing the proper assessment of kappa(0). This effect is more pronounced for soft soil sites. In this way, kappa(coda)(0) can serve as a proxy for the regional hard-rock kappa(0) at the reference sites.}, language = {en} } @article{ZhuCottonPilz2019, author = {Zhu, Chuanbin and Cotton, Fabrice Pierre and Pilz, Marco}, title = {Testing the Depths to 1.0 and 2.5 km/s Velocity Isosurfaces in a Velocity Model for Japan and Implications for Ground-Motion Modeling}, series = {Bulletin of the Seismological Society of America}, volume = {109}, journal = {Bulletin of the Seismological Society of America}, number = {6}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120190016}, pages = {2710 -- 2721}, year = {2019}, abstract = {In the Next Generation Attenuation West2 (NGA-West2) project, a 3D subsurface structure model (Japan Seismic Hazard Information Station [J-SHIS]) was queried to establish depths to 1.0 and 2.5 km/s velocity isosurfaces for sites without depth measurement in Japan. In this article, we evaluate the depth parameters in the J-SHIS velocity model by comparing them with their corresponding site-specific depth measurements derived from selected KiK-net velocity profiles. The comparison indicates that the J-SHIS model underestimates site depths at shallow sites and overestimates depths at deep sites. Similar issues were also identified in the southern California basin model. Our results also show that these underestimations and over-estimations have a potentially significant impact on ground-motion prediction using NGA-West2 ground-motion models (GMMs). Site resonant period may be considered as an alternative to depth parameter in the site term of a GMM.}, language = {en} } @article{ZhuPilzCotton2019, author = {Zhu, Chuanbin and Pilz, Marco and Cotton, Fabrice Pierre}, title = {Which is a better proxy, site period or depth to bedrock, in modelling linear site response in addition to the average shear-wave velocity?}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {18}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-019-00738-6}, pages = {797 -- 820}, year = {2019}, abstract = {This study aims to identify the best-performing site characterization proxy alternative and complementary to the conventional 30 m average shear-wave velocity V-S30, as well as the optimal combination of proxies in characterizing linear site response. Investigated proxies include T-0 (site fundamental period obtained from earthquake horizontal-to-vertical spectral ratios), V-Sz (measured average shear-wave velocities to depth z, z = 5, 10, 20 and 30 m), Z(0.8) and Z(1.0) (measured site depths to layers having shear-wave velocity 0.8 and 1.0 km/s, respectively), as well as Z(x-infer) (inferred site depths from a regional velocity model, x = 0.8 and 1.0, 1.5 and 2.5 km/s). To evaluate the performance of a site proxy or a combination, a total of 1840 surface-borehole recordings is selected from KiK-net database. Site amplifications are derived using surface-to-borehole response-, Fourier- and cross-spectral ratio techniques and then are compared across approaches. Next, the efficacies of 7 single-proxies and 11 proxy-pairs are quantified based on the site-to-site standard deviation of amplification residuals of observation about prediction using the proxy or the pair. Our results show that T-0 is the best-performing single-proxy among T-0, Z(0.8), Z(1.0) and V-Sz. Meanwhile, T-0 is also the best-performing proxy among T-0, Z(0.8), Z(1.0) and Z(x-infer) complementary to V-S30 in accounting for the residual amplification after V-S30-correction. Besides, T-0 alone can capture most of the site effects and should be utilized as the primary site indicator. Though (T-0, V-S30) is the best-performing proxy pair among (V-S30, T-0), (V-S30, Z(0.8)), (V-S30, Z(1.0)), (V-S30, Z(x-infer)) and (T-0, V-Sz), it is only slightly better than (T-0, V-S20). Considering both efficacy and engineering utility, the combination of T-0 (primary) and V-S20 (secondary) is recommended. Further study is needed to test the performances of various proxies on sites in deep sedimentary basins.}, language = {en} } @article{ZhuPilzCotton2020, author = {Zhu, Chuanbin and Pilz, Marco and Cotton, Fabrice Pierre}, title = {Evaluation of a novel application of earthquake HVSR in site-specific amplification estimation}, series = {Soil dynamics and earthquake engineering}, volume = {139}, journal = {Soil dynamics and earthquake engineering}, publisher = {Elsevier}, address = {Oxford}, issn = {0267-7261}, doi = {10.1016/j.soildyn.2020.106301}, pages = {14}, year = {2020}, abstract = {Ground response analyses (GRA) model the vertical propagations of SH waves through flat-layered media (1DSH) and are widely carried out to evaluate local site effects in practice. Horizontal-to-vertical spectral ratio (HVSR) technique is a cost-effective approach to extract certain site-specific information, e.g., site fundamental frequency (f(0)), but HVSR values cannot be directly used to approximate the levels of S-wave amplifications. Motivated by the work of Kawase et al. (2019), we propose a procedure to correct earthquake HVSR amplitudes for direct amplification estimations. The empirical correction compensates HVSR by generic vertical amplification spectra categorized by the vertical fundamental frequency (f(0v)) via kappa-means clustering. In this investigation, we evaluate the effectiveness of the corrected HVSR in approximating observed linear amplifications in comparison with 1DSH modellings. We select a total of 90 KiK-net (Kiban Kyoshin network) surface-downhole sites which are found to have no velocity contrasts below their boreholes and thus of which surface-to-borehole spectral ratios (SBSRs) can be taken as their empirical transfer functions (ETFs). 1DSH-based theoretical transfer functions (TTFs) are computed in the linear domain considering uncertainties in Vs profiles through randomizations. Five goodness-of-fit metrics are adopted to gauge the closeness between observed (ETF) and predicted (i.e., TTF and corrected HVSR) amplifications in both amplitude and spectral shape over frequencies from f(0) to 25 Hz. We find that the empirical correction to HVSR is highly effective and achieves a "good match" in both spectral shape and amplitude at the majority of the 90 KiK-net sites, as opposed to less than one-third for the 1DSH modelling. In addition, the empirical correction does not require a velocity model, which GRAs require, and thus has great potentials in seismic hazard assessments.}, language = {en} } @article{ZhuCottonPilz2020, author = {Zhu, Chuanbin and Cotton, Fabrice Pierre and Pilz, Marco}, title = {Detecting site resonant frequency using HVSR}, series = {Bulletin of the Seismological Society of America : BSSA}, volume = {110}, journal = {Bulletin of the Seismological Society of America : BSSA}, number = {2}, publisher = {Seismological Society of America}, address = {El Cerito, Calif.}, issn = {0037-1106}, doi = {10.1785/0120190186}, pages = {427 -- 440}, year = {2020}, abstract = {In this investigation, we examine the uncertainties using the horizontal-to-vertical spectral ratio (HVSR) technique on earthquake recordings to detect site resonant frequencies at 207 KiK-net sites. Our results show that the scenario dependence of response (pseudospectral acceleration) spectral ratio could bias the estimates of resonant frequencies for sites having multiple significant peaks with comparable amplitudes. Thus, the Fourier amplitude spectrum (FAS) should be preferred in computing HVSR. For more than 80\% of the investigated sites, the first peak (in the frequency domain) on the average HVSR curve over multiple sites coincides with the highest peak. However, for sites with multiple peaks, the highest peak frequency (f(p)) is less susceptible to the selection criteria of significant peaks and the extent of smoothing to spectrum than the first peak frequency (f(0)). Meanwhile, in comparison to the surface-to-borehole spectral ratio, f(0) tends to underestimate the predominant frequency (at which the largest amplification occurs) more than f(p). In addition, in terms of characterizing linear site response, f(p) shows a better overall performance than f(0). Based on these findings, we thus recommend that seismic network operators provide f(p) on the average HVSRFAS curve as a priority, ideally together with the average HVSRFAS curve in site characterization.}, language = {en} }