@article{PilzCottonZhu2021, author = {Pilz, Marco and Cotton, Fabrice Pierre and Zhu, Chuanbin}, title = {How much are sites affected by 2-D and 3-D site effects?}, series = {Geophysical journal international}, volume = {228}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab454}, pages = {1992 -- 2004}, year = {2021}, abstract = {1-D site response analysis dominates earthquake engineering practice, while local 2-D/3-D models are often required at sites where the site response is complex. For such sites, the 1-D representation of the soil column can account neither for topographic effects or dipping layers nor for locally generated horizontally propagating surface waves. It then remains a crucial task to identify whether the site response can be modelled sufficiently precisely by 1-D analysis. In this study we develop a method to classify sites according to their 1-D or 2-D/3-D nature. This classification scheme is based on the analysis of surface earthquake recordings and the evaluation of the variability and similarity of the horizontal Fourier spectra. The taxonomy is focused on capturing significant directional dependencies and interevent variabilities indicating a more probable 2-D/3-D structure around the site causing the ground motion to be more variable. While no significant correlation of the 1-D/3-D site index with environmental parameters and site proxies seems to exist, a reduction in the within-site (single-station) variability is found. The reduction is largest (up to 20 per cent) for purely 1-D sites. Although the taxonomy system is developed using surface stations of the KiK-net network in Japan as considerable additional information is available, it can also be applied to any (non-downhole array) site.}, language = {en} } @article{WatheletGuillierRouxetal.2018, author = {Wathelet, Marc and Guillier, B. and Roux, P. and Cornou, C. and Ohrnberger, Matthias}, title = {Rayleigh wave three-component beamforming}, series = {Geophysical journal international}, volume = {215}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggy286}, pages = {507 -- 523}, year = {2018}, abstract = {The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.}, language = {en} } @article{LontsiGarciaJerezCamiloMolinaVillegasetal.2019, author = {Lontsi, Agostiny Marrios and Garcia-Jerez, Antonio and Camilo Molina-Villegas, Juan and Jose Sanchez-Sesma, Francisco and Molkenthin, Christian and Ohrnberger, Matthias and Kr{\"u}ger, Frank and Wang, Rongjiang and Fah, Donat}, title = {A generalized theory for full microtremor horizontal-to-vertical [H/V(z,f)] spectral ratio interpretation in offshore and onshore environments}, series = {Geophysical journal international}, volume = {218}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggz223}, pages = {1276 -- 1297}, year = {2019}, abstract = {Advances in the field of seismic interferometry have provided a basic theoretical interpretation to the full spectrum of the microtremor horizontal-to-vertical spectral ratio [H/V(f)]. The interpretation has been applied to ambient seismic noise data recorded both at the surface and at depth. The new algorithm, based on the diffuse wavefield assumption, has been used in inversion schemes to estimate seismic wave velocity profiles that are useful input information for engineering and exploration seismology both for earthquake hazard estimation and to characterize surficial sediments. However, until now, the developed algorithms are only suitable for on land environments with no offshore consideration. Here, the microtremor H/V(z, f) modelling is extended for applications to marine sedimentary environments for a 1-D layered medium. The layer propagator matrix formulation is used for the computation of the required Green's functions. Therefore, in the presence of a water layer on top, the propagator matrix for the uppermost layer is defined to account for the properties of the water column. As an application example we analyse eight simple canonical layered earth models. Frequencies ranging from 0.2 to 50 Hz are considered as they cover a broad wavelength interval and aid in practice to investigate subsurface structures in the depth range from a few meters to a few hundreds of meters. Results show a marginal variation of 8 per cent at most for the fundamental frequency when a water layer is present. The water layer leads to variations in H/V peak amplitude of up to 50 per cent atop the solid layers.}, language = {en} } @article{CristianoMeierKruegeretal.2016, author = {Cristiano, L. and Meier, T. and Kr{\"u}ger, F. and Keers, H. and Weidle, C.}, title = {Teleseismic P-wave polarization analysis at the Grafenberg array}, series = {Geophysical journal international}, volume = {207}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggw339}, pages = {1456 -- 1471}, year = {2016}, abstract = {P-wave polarization at the Grafenberg array (GRF) in southern Germany is analysed in terms of azimuthal deviations and deviations in the vertical polarization using 20 yr of broad-band recordings. An automated procedure for estimating P-wave polarization parameters is suggested, based on the definition of a characteristic function, which evaluates the polarization angles and their time variability as well as the amplitude, linearity and the signal-to-noise ratio of the P wave. P-wave polarization at the GRF array is shown to depend mainly on frequency and backazimuth and only slightly on epicentral distance indicating depth-dependent local anisotropy and lateral heterogeneity. A harmonic analysis is applied to the azimuthal anomalies to analyse their periodicity as a function of backazimuth. The dominant periods are 180A degrees A and 360A degrees. At low frequencies, between 0.03 and 0.1 Hz, the observed fast directions of azimuthal anisotropy inferred from the 180A degrees A periodicity are similar across the array. The average fast direction of azimuthal anisotropy at these frequencies is N20A degrees E with an uncertainty of about 8A degrees A and is consistent with fast directions of Pn-wave propagation. Lateral velocity gradients determined for the low-frequency band are compatible with the Moho topography of the area. A more complex pattern in the horizontal fast axis orientation beneath the GRF array is observed in the high-frequency band between 0.1 and 0.5 Hz, and is attributed to anisotropy in the upper crust. A remarkable rotation of the horizontal fast axis orientation across the suture between the geological units Moldanubicum and Saxothuringicum is observed. In contrast, the 360A degrees A periodicity at high frequencies is rather consistent across the array and may either point to lower velocities in the upper crust towards the Bohemian Massif and/or to anisotropy dipping predominantly in the NE-SW direction. Altogether, P-wave polarization analysis indicates the presence of layered lithospheric anisotropy in the area of the GRF array. Seismic anisotropy is more variable in the brittle upper crust compared to lower crustal and subcrustal depths.}, language = {en} } @article{PaloTilmannKruegeretal.2014, author = {Palo, Mauro and Tilmann, Frederik and Kr{\"u}ger, Frank and Ehlert, Lutz and Lange, Dietrich}, title = {High-frequency seismic radiation from Maule earthquake (M-w 8.8, 2010 February 27) inferred from high-resolution backprojection analysis}, series = {Geophysical journal international}, volume = {199}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggu311}, pages = {1058 -- 1077}, year = {2014}, abstract = {We track a bilateral rupture propagation lasting similar to 160 s, with its dominant branch rupturing northeastwards at about 3 kms(-1). The area of maximum energy emission is offset from the maximum coseismic slip but matches the zone where most plate interface aftershocks occur. Along dip, energy is preferentially released from two disconnected interface belts, and a distinct jump from the shallower belt to the deeper one is visible after about 20 s from the onset. However, both belts keep on being active until the end of the rupture. These belts approximately match the position of the interface aftershocks, which are split into two clusters of events at different depths, thus suggesting the existence of a repeated transition from stick-slip to creeping frictional regime.}, language = {en} } @article{TranThanhTuanScherbaumMalischewsky2011, author = {Tran Thanh Tuan, and Scherbaum, Frank and Malischewsky, Peter G.}, title = {On the relationship of peaks and troughs of the ellipticity (H/V) of Rayleigh waves and the transmission response of single layer over half-space models}, series = {Geophysical journal international}, volume = {184}, journal = {Geophysical journal international}, number = {2}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2010.04863.x}, pages = {793 -- 800}, year = {2011}, abstract = {One of the key challenges in the context of local site effect studies is the determination of frequencies where the shakeability of the ground is enhanced. In this context, the H/V technique has become increasingly popular and peak frequencies of H/V spectral ratio are sometimes interpreted as resonance frequencies of the transmission response. In the present study, assuming that Rayleigh surface wave is dominant in H/V spectral ratio, we analyse theoretically under which conditions this may be justified and when not. We focus on 'layer over half-space' models which, although seemingly simple, capture many aspects of local site effects in real sedimentary structures. Our starting point is the ellipticity of Rayleigh waves. We use the exact formula of the H/V-ratio presented by Malischewsky \& Scherbaum (2004) to investigate the main characteristics of peak and trough frequencies. We present a simple formula illustrating if and where H/V-ratio curves have sharp peaks in dependence of model parameters. In addition, we have constructed a map, which demonstrates the relation between the H/V-peak frequency and the peak frequency of the transmission response in the domain of the layer's Poisson ratio and the impedance contrast. Finally, we have derived maps showing the relationship between the H/V-peak and trough frequency and key parameters of the model such as impedance contrast. These maps are seen as diagnostic tools, which can help to guide the interpretation of H/V spectral ratio diagrams in the context of site effect studies.}, language = {en} } @article{DiGiacomoBindiParolaietal.2011, author = {Di Giacomo, Domenico and Bindi, Dino and Parolai, Stefano and Oth, Adrien}, title = {Residual analysis of teleseismic P-wave energy magnitude estimates: inter- and intrastation variability}, series = {Geophysical journal international}, volume = {185}, journal = {Geophysical journal international}, number = {3}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2011.05019.x}, pages = {1444 -- 1454}, year = {2011}, abstract = {P>Computing the magnitude of an earthquake requires correcting for the propagation effects from the source to the receivers. This is often accomplished by performing numerical simulations using a suitable Earth model. In this work, the energy magnitude M(e) is considered and its determination is performed using theoretical spectral amplitude decay functions over teleseismic distances based on the global Earth model AK135Q. Since the high frequency part (above the corner frequency) of the source spectrum has to be considered in computing M(e), the influence of propagation and site effects may not be negligible and they could bias the single station M(e) estimations. Therefore, in this study we assess the inter- and intrastation distributions of errors by considering the M(e) residuals computed for a large data set of earthquakes recorded at teleseismic distances by seismic stations deployed worldwide. To separate the inter- and intrastation contribution of errors, we apply a maximum likelihood approach to the M(e) residuals. We show that the interstation errors (describing a sort of site effect for a station) are within +/- 0.2 magnitude units for most stations and their spatial distribution reflects the expected lateral variation affecting the velocity and attenuation of the Earth's structure in the uppermost layers, not accounted for by the 1-D AK135Q model. The variance of the intrastation error distribution (describing the record-to-record component of variability) is larger than the interstation one (0.240 against 0.159), and the spatial distribution of the errors is not random but shows specific patterns depending on the source-to-station paths. The set of coefficients empirically determined may be used in the future to account for the heterogeneities of the real Earth not considered in the theoretical calculations of the spectral amplitude decay functions used to correct the recorded data for propagation effects.}, 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{HobigerCornouWatheletetal.2013, author = {Hobiger, M. and Cornou, C. and Wathelet, M. and Di Giulio, G. and Knapmeyer-Endrun, B. and Renalier, F. and Bard, Pierre-Yves and Savvaidis, Alexandros and Hailemikael, S. and Le Bihan, N. and Ohrnberger, Matthias and Theodoulidis, N.}, title = {Ground structure imaging by inversions of Rayleigh wave ellipticity sensitivity analysis and application to European strong-motion sites}, series = {Geophysical journal international}, volume = {192}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggs005}, pages = {207 -- 229}, year = {2013}, abstract = {The knowledge of the local soil structure is important for the assessment of seismic hazards. A widespread, but time-consuming technique to retrieve the parameters of the local underground is the drilling of boreholes. Another way to obtain the shear wave velocity profile at a given location is the inversion of surface wave dispersion curves. To ensure a good resolution for both superficial and deeper layers, the used dispersion curves need to cover a wide frequency range. This wide frequency range can be obtained using several arrays of seismic sensors or a single array comprising a large number of sensors. Consequently, these measurements are time-consuming. A simpler alternative is provided by the use of the ellipticity of Rayleigh waves. The frequency dependence of the ellipticity is tightly linked to the shear wave velocity profile. Furthermore, it can be measured using a single seismic sensor. As soil structures obtained by scaling of a given model exhibit the same ellipticity curve, any inversion of the ellipticity curve alone will be ambiguous. Therefore, additional measurements which fix the absolute value of the shear wave velocity profile at some points have to be included in the inversion process. Small-scale spatial autocorrelation measurements or MASW measurements can provide the needed data. Using a theoretical soil structure, we show which parts of the ellipticity curve have to be included in the inversion process to get a reliable result and which parts can be omitted. Furthermore, the use of autocorrelation or high-frequency dispersion curves will be highlighted. The resulting guidelines for inversions including ellipticity data are then applied to real data measurements collected at 14 different sites during the European NERIES project. It is found that the results are in good agreement with dispersion curve measurements. Furthermore, the method can help in identifying the mode of Rayleigh waves in dispersion curve measurements.}, language = {en} }