Single-station cross-correlation analysis of ambient seismic noise: application to stations in the surroundings of the 2008 Iwate-Miyagi Nairiku earthquake

We introduce the single-station cross-correlation (SC) technique of processing ambient seismic noise and compare its results with the established cross-correlation (CC) and autocorrelation (AC) techniques. While CC is the correlation of the signals of two seismic stations with each other and AC is the correlation of a signal with itself, SC is the correlation of two different components of a single three-component seismic sensor. The comparison of the three different correlation techniques shows that CCs give the best results at frequencies below 0.5 Hz and that SCs give the best results at higher frequencies. In all three processing techniques, ambient seismic noise is correlated in order to reconstruct the Green's function describing the wave propagation between the first and the second sensor. By relating the coda parts of the daily Green's functions with the long-term reference Green's functions, shear wave velocity changes are determined. Here, we apply this technique to the data of 20 seismic stations in the surroundings of the fault zone of the Iwate-Miyagi Nairiku earthquake (MW=6.9), which occurred on 2008 June 13, UTC (2008 June 14, Japan Standard Time) in the northern part of the Japanese island Honshu. The data range from 2008 January to 2011 June and therefore include the Tohoku earthquake (MW=9.0), which occurred on 2011 March 11, off the coast of northern Honshu. The data are analysed in five different frequency ranges between 0.125 and 4.0 Hz. The data show coseismic velocity changes for both earthquakes followed by a post-seismic velocity recovery. In general, the coseismic velocity changes increase with frequency. For the Iwate-Miyagi Nairiku earthquake, the strongest velocity changes occur close to the fault zone. Quickly recovering coseismic velocity changes can be separated from changes not recovering during the study period. For the Tohoku earthquake, the complete area is affected by coseismic velocity changes. A modelling of the depth of the coseismic velocity changes indicates that the Iwate-Miyagi Nairiku earthquake can be explained either by large shallow velocity changes or by small, but deep changes. For one station, the observations can only be explained by assuming deeper changes. For the Tohoku earthquake, the modelling shows that different parts of the study area are affected in different ways, some showing shallow changes, others deeper changes. Furthermore, seasonal velocity variations occur, which are compatible for the different stations above 0.5 Hz, with velocity maxima in autum

Geophysical surveys for the characterization of the seismic local response at instrumented sites : a case study from a station of the Swiss strong motion network

Site characterization surveying is one of the key efforts for the understanding of local amplification effects in earthquake engineering. In this perspective, geophysical measurements for the estimation of the physical properties of the subsurface at instrumented sites are particularly important. For instance, the comparison between simulated site amplification based on the retrieved Vs profile and empirical local response observed at the station allows assessing whether the site is characterized by a purely 1D response or additional 2D-3D effects are present. We present a case study involving a site characterization measurement at a station of the Swiss strong-motion network. The survey includes active seismic data acquisition along two geophone lines, of 15 and 115 m length, for the imaging of the near-surface and deeper layers, respectively. The acquired data were interpreted in terms of P-wave refraction and Rayleigh-wave multi-modal propagation analysis. The reconstructed Vs-Vp profiles achieve a high level of accuracy over a wide depth range and allow modelling the SH-transfer function of the local soil column; the latter is in excellent agreement with the local amplification function estimated by means of empirical spectral modelling of the station’s recordings. The agreement confirms a purely 1D site response for the surveyed station

Current status of the IAG working group 4.3.7 on geodetic GNSS-R

Presentación realizada online en el Scientific Assembly of the International Association of Geodesy (2021) celebrado del 28 de junio al 2 de julio en Beijing

Polarisation des ondes de surface : caractérisation, inversion et application à l'étude de l'aléa sismique

The seismic hazard of a given site can be largely influenced by site effects. In order to estimate these effects, the local soil structure and the wave field properties have to be investigated. This can be done using surface waves (Love and Rayleigh waves) since their properties (dispersion curves, Rayleigh wave ellipticity) are closely linked to the soil structure. The key parameter for the correct wave type identification is the polarization. In a first part of this thesis, three new methods to estimate the polarization parameters of surface waves have been developed. Two methods, DELFI and RayDec, estimate the Rayleigh wave ellipticity using the recordings of a single seismic station. The third method, MUSIQUE, is a further developed version of the MUSIC algorithm, and uses seismic array recordings to discriminate between Love and Rayleigh waves and estimate their respective properties. In the second part of the thesis, a theoretical investigation of the inversion of ellipticity curves shows which parts of an ellipticity curve carry the important information on the soil structure and by which means the inversion can be disambiguated. The developed inversion scheme is then tested on real data measurements obtained at 14 different European sites. Finally, the seismological wave field in the Californian Santa Clara valley is investigated applying MUSIQUE to seismic array recordings for 22 earthquakes. The azimuthal energy repartition between the different wave types is investigated and the respective dispersion curves as well as the Rayleigh wave ellipticity curve are estimated. This study outlines the importance of diffracted surface waves in the recorded wave field. STARL'aléa sismique d'un site donné peut être largement influencé par des effets de site. Afin d'évaluer ces effets, la structure locale du sous-sol ainsi que les propriétés du champ d'onde doivent être étudiées. Les ondes de surface (ondes de Love et de Rayleigh) s'avèrent utiles, leurs propriétés (courbes de dispersion, ellipticité des ondes de Rayleigh) étant directement liées à la structure du sous-sol. Le paramètre clé pour l'identification du type d'onde est la polarisation. Dans la première partie de la thèse, de nouvelles méthodes pour l'estimation de paramètres de polarisation d'ondes de surface sont développées. Deux méthodes, DELFI et RayDec, estiment l'ellipticité des ondes de Rayleigh à partir d'enregistrements d'un seul capteur sismique. La troisième méthode, MUSIQUE, est basée sur la méthode MUSIC et utilise les enregistrements multi-composantes de réseaux sismiques afin de distinguer ondes de Love et ondes de Rayleigh et d'estimer leurs propriétés. Dans la deuxième partie de la thèse, une étude théorique de l'inversion de courbes d'ellipticité montre quelles parties de ces courbes véhiculent les informations importantes sur la structure du sol et comment l'inversion peut être améliorée. Le schéma d'inversion résultant est alors testé en l'appliquant à des données réelles mesurées pour 14 sites européens. Finalement, 22 séismes enregistrés par un réseau de capteurs dans la vallée de Santa Clara en Californie sont analysés par MUSIQUE. La répartition azimutale, les courbes de dispersion, la courbe d'ellipticité et les énergies des différents types d'ondes sont analysées et soulignent l'importance des ondes de surface diffractées dans le champ d'ondes enregistré

Polarization of surface waves : characterization, inversion and application to seismic hazard assessment

L'aléa sismique d'un site donné peut être largement influencé par des effets de site. Afin d'évaluer ces effets, la structure locale du sous-sol ainsi que les propriétés du champ d'onde doivent être étudiées. Les ondes de surface (ondes de Love et de Rayleigh) s'avèrent utiles, leurs propriétés (courbes de dispersion, ellipticité des ondes de Rayleigh) étant directement liées à la structure du sous-sol. Le paramètre clé pour l'identification du type d'onde est la polarisation. Dans la première partie de la thèse, de nouvelles méthodes pour l'estimation de paramètres de polarisation d'ondes de surface sont développées. Deux méthodes, DELFI et RayDec, estiment l'ellipticité des ondes de Rayleigh à partir d'enregistrements d'un seul capteur sismique. La troisième méthode, MUSIQUE, est basée sur la méthode MUSIC et utilise les enregistrements multi-composantes de réseaux sismiques afin de distinguer ondes de Love et ondes de Rayleigh et d'estimer leurs propriétés. Dans la deuxième partie de la thèse, une étude théorique de l'inversion de courbes d'ellipticité montre quelles parties de ces courbes véhiculent les informations importantes sur la structure du sol et comment l'inversion peut être améliorée. Le schéma d'inversion résultant est alors testé en l'appliquant à des données réelles mesurées pour 14 sites européens. Finalement, 22 séismes enregistrés par un réseau de capteurs dans la vallée de Santa Clara en Californie sont analysés par MUSIQUE. La répartition azimutale, les courbes de dispersion, la courbe d'ellipticité et les énergies des différents types d'ondes sont analysées et soulignent l'importance des ondes de surface diffractées dans le champ d'ondes enregistré.The seismic hazard of a given site can be largely influenced by site effects. In order to estimate these effects, the local soil structure and the wave field properties have to be investigated. This can be done using surface waves (Love and Rayleigh waves) since their properties (dispersion curves, Rayleigh wave ellipticity) are closely linked to the soil structure. The key parameter for the correct wave type identification is the polarization. In a first part of this thesis, three new methods to estimate the polarization parameters of surface waves have been developed. Two methods, DELFI and RayDec, estimate the Rayleigh wave ellipticity using the recordings of a single seismic station. The third method, MUSIQUE, is a further developed version of the MUSIC algorithm, and uses seismic array recordings to discriminate between Love and Rayleigh waves and estimate their respective properties. In the second part of the thesis, a theoretical investigation of the inversion of ellipticity curves shows which parts of an ellipticity curve carry the important information on the soil structure and by which means the inversion can be disambiguated. The developed inversion scheme is then tested on real data measurements obtained at 14 different European sites. Finally, the seismological wave field in the Californian Santa Clara valley is investigated applying MUSIQUE to seismic array recordings for 22 earthquakes. The azimuthal energy repartition between the different wave types is investigated and the respective dispersion curves as well as the Rayleigh wave ellipticity curve are estimated. This study outlines the importance of diffracted surface waves in the recorded wave field. STA

Retrieval of Rayleigh wave ellipticity from ambient vibration recordings

ISSN:0956-540XISSN:1365-246

Ambient vibration analysis on seismic arrays to investigate the properties of the upper crust: an example from Herdern in Switzerland

The difficulty and the high cost to assess the subsurface properties led to the development of several geophysical techniques. Generally, the focus of a site study is the reconstruction of the S-wave velocity profile down to few tens to hundreds of metres (e.g. 30–300 m), but not the investigation of deeper structures, such as the transition to the crystalline basement. However, such deeper structures are of interest when seismic hazard products have to relate to a reference rock-velocity profile, for example in regional seismic hazard assessment and microzonation studies. To estimate the S-wave velocity profiles down to several kilometres, we study the potential of Rayleigh and Love waves at low (down to 0.1 Hz) and high (up to 20 Hz) frequencies using two seismic arrays of increasing size. The small array, with a maximum inter-station distance of 900 m and a recording time of 3 hr, was aimed at constraining the shallow subsurface down to about 350–400 m, while the big one, with a maximum inter-station distance of more than 29 km and 23 hr of recording had the goal to constrain the deeper structure. The arrays were deployed in northern Switzerland (east of the village of Herdern) within the Swiss Molasse basin, a sedimentary basin north of the Alps stretching from Lake Constance to Lake Geneva; its thickness increases from 800 to 900 m in the northeast to more than 5 km in the southwest. The seismic data recorded by the two arrays were analysed using the techniques developed for the analysis of small-aperture arrays. The results were inverted for the S-wave velocity profile in two steps: first, the Rayleigh and Love wave phase dispersion curves were inverted together. Secondly, the previous dispersion curves were jointly inverted with the measured Rayleigh wave ellipticity angle. The resulting S-wave velocity profiles are similar and show agreement with the available geological and geophysical data, confirming the potential of surface waves to investigate deep structures. Moreover, our analysis proves the feasibility of site characterization techniques to large arrays and the possibility to estimate the P- and S-wave velocity profiles down to 5 km, deeper than the contrast between Molasse basin and crystalline rock at around 2.1 km.ISSN:0956-540XISSN:1365-246

Polarisation des ondes de surface (caractérisation, inversion et application à l'étude de l'aléa sismique)

L'aléa sismique d'un site donné peut être largement influencé par des effets de site. Afin d'évaluer ces effets, la structure locale du sous-sol ainsi que les propriétés du champ d'onde doivent être étudiées. Les ondes de surface (ondes de Love et de Rayleigh) s'avèrent utiles, leurs propriétés (courbes de dispersion, ellipticité des ondes de Rayleigh) étant directement liées à la structure du sous-sol. Le paramètre clé pour l'identification du type d'onde est la polarisation. Dans la première partie de la thèse, de nouvelles méthodes pour l'estimation de paramètres de polarisation d'ondes de surface sont développées. Deux méthodes, DELFI et RayDec, estiment l'ellipticité des ondes de Rayleigh à partir d'enregistrements d'un seul capteur sismique. La troisième méthode, MUSIQUE, est basée sur la méthode MUSIC et utilise les enregistrements multi-composantes de réseaux sismiques afin de distinguer ondes de Love et ondes de Rayleigh et d'estimer leurs propriétés. Dans la deuxième partie de la thèse, une étude théorique de l'inversion de courbes d'ellipticité montre quelles parties de ces courbes véhiculent les informations importantes sur la structure du sol et comment l'inversion peut être améliorée. Le schéma d'inversion résultant est alors testé en l'appliquant à des données réelles mesurées pour 14 sites européens. Finalement, 22 séismes enregistrés par un réseau de capteurs dans la vallée de Santa Clara en Californie sont analysés par MUSIQUE. La répartition azimutale, les courbes de dispersion, la courbe d'ellipticité et les énergies des différents types d'ondes sont analysées et soulignent l'importance des ondes de surface diffractées dans le champ d'ondes enregistré.The seismic hazard of a given site can be largely influenced by site effects. In order to estimate these effects, the local soil structure and the wave field properties have to be investigated. This can be done using surface waves (Love and Rayleigh waves) since their properties (dispersion curves, Rayleigh wave ellipticity) are closely linked to the soil structure. The key parameter for the correct wave type identification is the polarization. In a first part of this thesis, three new methods to estimate the polarization parameters of surface waves have been developed. Two methods, DELFI and RayDec, estimate the Rayleigh wave ellipticity using the recordings of a single seismic station. The third method, MUSIQUE, is a further developed version of the MUSIC algorithm, and uses seismic array recordings to discriminate between Love and Rayleigh waves and estimate their respective properties. In the second part of the thesis, a theoretical investigation of the inversion of ellipticity curves shows which parts of an ellipticity curve carry the important information on the soil structure and by which means the inversion can be disambiguated. The developed inversion scheme is then tested on real data measurements obtained at 14 different European sites. Finally, the seismological wave field in the Californian Santa Clara valley is investigated applying MUSIQUE to seismic array recordings for 22 earthquakes. The azimuthal energy repartition between the different wave types is investigated and the respective dispersion curves as well as the Rayleigh wave ellipticity curve are estimated. This study outlines the importance of diffracted surface waves in the recorded wave field. STARSAVOIE-SCD - Bib.électronique (730659901) / SudocGRENOBLE1/INP-Bib.électronique (384210012) / SudocGRENOBLE2/3-Bib.électronique (384219901) / SudocSudocFranceF

Ambient Vibration Analysis on Large Scale Arrays When Lateral Variations Occur in the Subsurface: A Study Case in Switzerland

The ambient vibration analysis is a non-invasive and low-cost technique used in site characterization studies to reconstruct the subsurface velocity structure. Depending on the goal of the research, the investigated depth ranges from tens to hundreds of meters. In this work, we aimed at investigating the deeper contrasts within the crust and in particular down to the sedimentary-rock basement transition located at thousands of meters of depth. To achieve this goal, three seismic arrays with minimum and maximum interstation distances of 7.9 m and 26.8 km were deployed around the village of Schafisheim. Schafisheim is located in the Swiss Molasse Basin, a sedimentary basin stretching from Lake Constance to Lake Geneva with a thickness ranging from 800 to 900 m in the north to 5 km in the south. To compute the multimodal dispersion curves for Rayleigh and Love waves and the Rayleigh wave ellipticity angles, the data were processed using two single-station and three array processing techniques. A preliminary analysis of the inversion results pointed out a good agreement with the fundamental modes of Rayleigh and Love waves used in the inversion and a quite strong disagreement with the higher modes. The impossibility to explain at the same time most of the dispersion curves was interpreted as the co-existence, within the investigated area, of portions of the subsurface with different geophysical properties. The hypothesis was confirmed by the Horizontal-to-Vertical spectral analysis (H/V) which indicated the presence of two distinguished areas. The observation allowed a new interpretation and the identification of the Rayleigh and Love wave fundamental modes and of the S-wave velocity profiles to be reconstructed for each investigated zone. It results in two S-wave velocity profiles with similar velocities down to 15 km deferring only in their shallow portions due to the occurrence of a low velocity zone at a depth of 50–150 m at the centre of the investigated area.ISSN:0033-4553ISSN:1420-9136ISSN:1557-736

Analysis of Rayleigh waves with circular wavefront: a maximum likelihood approach

Analysis of Rayleigh waves is an important task in seismology and geotechnical investigations. In fact, properties of Rayleigh waves such as velocity and polarization are important observables that carry information about the structure of the subsoil. Applications analysing Rayleigh waves include active and passive seismic surveys. In active surveys, there is a controlled source of seismic energy and the sensors are typically placed near the source. In passive surveys, there is not a controlled source, rather, seismic waves from ambient vibrations are analysed and the sources are assumed to be far outside the array, simplifying the analysis by the assumption of plane waves. Whenever the source is in the proximity of the array of sensors or even within the array it is necessary to model the wave propagation accounting for the circular wavefront. In addition, it is also necessary to model the amplitude decay due to geometrical spreading. This is the case of active seismic surveys in which sensors are located near the seismic source. In this work, we propose a maximum likelihood (ML) approach for the analysis of Rayleigh waves generated at a near source. Our statistical model accounts for the curvature of the wavefront and amplitude decay due to geometrical spreading. Using our method, we show applications on real data of the retrieval of Rayleigh wave dispersion and ellipticity. We employ arrays with arbitrary geometry. Furthermore, we show how it is possible to combine active and passive surveys. This enables us to enlarge the analysable frequency range and therefore the depths investigated. We retrieve properties of Rayleigh waves from both active and passive surveys and show the excellent agreement of the results from the two surveys. In our approach we use the same array of sensors for both the passive and the active survey. This greatly simplifies the logistics necessary to perform a survey.ISSN:0956-540XISSN:1365-246