# Characterization of site conditions (soil class, VS30, velocity profiles) for 33 stations from the French permanent accelerometric network (RAP) using surface-wave methods

Fabrice Hollender, Cécile Cornou, Aline Dechamp, Kaveh Oghalaei, Florence Renalier, Emeline Maufroy, Clément Burnouf, Sylvette Thomassin, Marc Wathelet, Pierre-Yves Bard, et al.

► To cite this version:

Fabrice Hollender, Cécile Cornou, Aline Dechamp, Kaveh Oghalaei, Florence Renalier, et al. Characterization of site conditions (soil class, VS30, velocity profiles) for 33 stations from the French permanent accelerometric network (RAP) using surface-wave methods. Bulletin of Earthquake Engineering, 2018, 16, pp.2337-2365. 10.1007/s10518-017-0135-5. hal-01693200

HAL Id: hal-01693200

https://hal.science/hal-01693200v1

Submitted on 21 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Characterization of site conditions (soil class, V<sub>S30</sub>,
velocity profiles) for 33 stations from the French
permanent accelerometric network (RAP) using surface-
wave methods

Fabrice Hollender<sup>1,2</sup> · Cécile Cornou<sup>2</sup> · Aline Dechamp<sup>3</sup> ·

Kaveh Oghalaei<fnref n="2" /> · Florence Renalier<fnref n="4" /> · Emeline Maufroy<fnref n="2" /> ·

Clément Burnouf<fnref n="5" /> · Sylvette Thomassin<fnref n="6" /> · Marc Wathelet<fnref n="2" /> ·

Pierre-Yves Bard<fnref n="2" /> · Vincent Boutin<fnref n="3" /> · Clément Desbordes<fnref n="1" /> ·

Isabelle Douste-Bacqué<fnref n="2" /> · Laetitia Foundotos<fnref n="1" /> ·

Cédric Guyonnet-Benaize<fnref n="1" /> · Vincent Perron<fnref n="1" /> ·

Julie Régnier<fnref n="7" /> · Agathe Roullé<fnref n="8" /> · Mickael Langlais<fnref n="2" /> ·

Deborah Sicilia<fnref n="9" />

Received: 14 October 2016/Accepted: 9 April 2017/Published online: 19 April 2017
© The Author(s) 2017. This article is an open access publication

**Abstract** Data provided by accelerometric networks are important for seismic hazard assessment. The correct use of accelerometric signals is conditioned by the station site metadata quality (i.e., soil class, V<sub>S30</sub>, velocity profiles, and other relevant information that can help to quantify site effects). In France, the permanent accelerometric network consists of about 150 stations. Thirty-three of these stations in the southern half of France have been characterized, using surface-wave-based methods that allow derivation of velocity profiles from dispersion curves of surface waves. The computation of dispersion curves and their subsequent inversion in terms of shear-wave velocity profiles has allowed estimation of V<sub>S30</sub> values and designation of soil classes, which include the corresponding uncertainties. From a methodological point of view, this survey leads to the following recommendations: (1) perform both active (multi-analysis surface waves) and passive (ambient

<figure><img src="image_3.png" /></figure>

<footnote>Fabrice Hollender
fabrice.hollender@cea.fr</footnote>

<footnote n="1">1 CEA, DEN, 13108 St Paul lez Durance, France</footnote>

<footnote n="2">2 ISTerre, CNRS, IRD, IFSTTAR, Univ. Grenoble Alpes, 38041 Grenoble, France</footnote>

<footnote n="3">³ CEA, DAM, DIF, 91297 Arpajon, France</footnote>

<footnote n="4">⁴ GeophyConsult, 159, quai des Allobroges, 73000 Chambéry, France</footnote>

<footnote n="5">⁵ SIXENSE Geophysics, 38300 Nivolas Vermelle, France</footnote>

<footnote n="6">⁶ Résonance Ingénieurs-Conseils SA, 1227 Carouge, Switzerland</footnote>

<footnote n="7">⁷ CEREMA, DTer Méditerranée, 06300 Nice, France</footnote>

<footnote n="8">⁸ BRGM, DRP/RSV, 45100 Orléans, France</footnote>

<footnote n="9">⁹ EdF, DPI-CEIDRE, 13097 Aix-en-Provence, France</footnote>
vibration arrays) measurements to derive dispersion curves in a broadband frequency range; (2) perform active acquisitions for both vertical (Rayleigh wave) and horizontal (Love wave) polarities. Even when the logistic contexts are sometimes difficult, the use of surface-wave-based methods is suitable for station-site characterization, even on rock sites. In comparison with previous studies that have mainly estimated $V_{S30}$ indirectly, the new values here are globally lower, but the EC8-A class sites remain numerous. However, even on rock sites, high frequency amplifications may affect accelerometric records, due to the shallow relatively softer layers.

**Keywords** Accelerometric network · Surface-wave methods · Velocity profile · $V_{S30}$

# 1 Introduction

Data provided by accelerometric networks are of great importance for seismic hazard assessment, especially for the derivation of ground-motion prediction equations (GMPEs), which are the basis of most deterministic or probabilistic seismic hazard analysis. However, the optimal use of accelerometric waveforms is closely linked to the quality of the station ‘site-metadata’, that needs to provide reliable information about local soil conditions, such as for the EC8 soil class (EUROCODE8 1998), $V_{S30}$, velocity profiles, and all of the relevant information that can help to quantify site effects. Unfortunately, these metadata are often only poorly known (Sandikkaya et al. 2013; Seyhan et al. 2014; Akkar et al. 2014). Soil classes or $V_{S30}$ parameters are indeed most often indirectly derived, such as through the use of geological maps, where a given lithology is associated to a given shear-wave velocity. However, this approach does not take into account possible weathering of geological material, nor the presence of shallow quaternary deposits that are not always indicated on geological maps. Another indirect approach is to correlate $V_{S30}$ and the slope of the topography (Allen and Wald 2009), which has been shown to require further calibration through additional local constraints to reach reliable prediction capabilities (Lemoine et al. 2012; Thompson et al. 2014; Yong 2016). Nevertheless, indirect approaches can also often lead to $V_{S30}$ estimates of high uncertainty. To tackle these issues, important efforts have been made in recent years to perform in situ site characterization of accelerometric stations: e.g., in Italy (Foti et al. 2011; Pileggi et al. 2011), Switzerland (Michel et al. 2014), and the USA (Yong et al. 2013).

In France, the ‘Réseau Accélérométrique Permanent’ (RAP; permanent accelerometric network) consists of approximately 150 stations (Guéguen et al. 2007; Pequegnat et al. 2008). Previous studies have already proposed soil classes or $V_{S30}$ for some or all of the RAP stations. For example, Régnier et al. (2010) gathered a wide range of pre-existing geophysical and geotechnical information in the vicinity of the RAP sites. Other studies (Drouet et al. 2008, 2010) have been aimed at estimation of $V_{S30}$ using generalized inversion methods.

The present paper addresses recent efforts that have led to the characterization of 33 of the RAP station sites that are located in the Pyrenees, Alps, Auvergne and Provence-Côte d’Azur regions (Fig. 1), using in situ systematic measurements of $V_S$ profiles. One key issue for this measurement campaign is the choice of methods to be used, which obviously impacts on the survey budget. The use of classical geotechnical methods such as cross-hole and down-hole measurements leads to high costs and high logistic constraints. Moreover, these methods can be limited in terms of the investigation depth. Alternatively, methods
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Map of the southern half of France with the locations of the RAP network stations and stations that benefited from the characterization presented here. Right upper frame zoom in on the Annecy area. Right lower frame zoom in on the Grenoble area</figcaption></figure>

based on analysis of Rayleigh-wave and/or Love-wave surface-wave dispersion characteristics can allow the derivation of shear-wave velocity profiles. Even if the capabilities of these methods have been debated for a long time, they provide good estimates when implemented with care, as demonstrated within the InterPacific project (Garofalo et al. 2016a, b; Foti et al. 2017). In particular, Garofalo et al. (2016b) outlined that $V_{S30}$ can be equally well estimated using ‘noninvasive’ approaches (e.g., surface-wave-based methods) and ‘invasive’ approaches (e.g., cross-holes, down-holes). However, Garofalo et al. (2016a) also showed that noninvasive methods are less accurate for locating a given lithological interface, low velocity layers or bedrock depth. They attribute this to surface wave depth resolution issues, non-uniqueness of the inversion procedures and/or the energy content of the signals used. The consequences of the uncertainties related to non-uniqueness of site responses have been a matter of discussion in recent years (Cornou et al. 2006; Foti et al. 2009; Socco et al. 2012; Boaga et al. 2013; Jakka et al. 2014; Comina and Foti 2015; Griffiths et al. 2016; Cox and Teague 2016), with some authors claiming that the great uncertainty in Vs profiles leads to high variability of the site responses, although others do not agree. However, considering that the bedrock velocity is well known (e.g., through large arrays or other indirect information) and the velocity profiles are coherent with other indicators (e.g., the ambient vibration H/V frequency peak), the impact of non-uniqueness on site responses remains acceptable. Therefore, with the main objective of the survey being the estimation of site effects, noninvasive surface-wave-based methods represent good candidates, and were chosen for the present study.

This study starts with a brief description of the chosen accelerometric stations, the methods used, and the related acquisition layout. Next, the different steps of processing are presented, which include data quality control, derivation of dispersion curves, and shear-wave velocity profile inversions. Finally, the results are shown, and we draw implications from these results, as well as for the overall methodological feedback.
## 2 Choice of rap sites

The choice of the accelerometric sites to be characterized was a compromise between the wish to characterize sites that produce the largest number of accelerograms, and the wish to characterize sites that were used as ‘reference sites’ in Drouet et al. (2010). Finally, some sites were ruled out due to logistical constraints (e.g., difficulties in site access).

The acquisitions were carried out in different phases, with each new phase taking into account the feedback of the previous phase, in terms of the survey plan, geophysical material, logistic organization in the field, and data processing. This explains the successive changes in the standard acquisition layout, as well as the overall level of successfully recorded data for each site (considering that we usually dedicated one working field day to one site, with a team of 6 to 8 people).

The first site that was characterized was NALS station, in the city of Nice (April, 2012), followed by three sites in Auvergne (OCLD, OCOR, OCOL) in late April 2012, and nine sites in the Pyrenees (EPF, PYAS, PYAT, PYBB, PYLI, PYLL, PYLO, PYLU, PYOR) in September 2012. In July and August 2014, five sites were characterized in the Provence-Côte d’Azur area (IRPV, OGCA, IRVG, CALF, NBOR) and 15 in the French Alps (OGDH, OGPC, OGCH, OGBL, SURF, OGFO, OGMU, OGIM, OGLE, GRN, OGAN, OGAP, OGME, OGMA, OGSI). Some necessary additional measurements were performed from November 2014 to January 2015 to complete the datasets for NBOR, OGBL, OGCA and IRVG stations. Table 1 gives the main geological features of the 33 studied stations, and Fig. 1 shows their locations. The former OCKE station was located very close to the OCLD one. The results obtained for OCLD are also valid for OCKE. The station OGFO corresponds to the Montbonnot borehole site. It involves three stations: a first one at the surface (also known as OGFH), a second one at 40 m depth (also known as OGFM) and a third one at 556 m depth (also known as OGFB).

## 3 Acquisition layout

The acquisition layout basically consisted of both active surface-wave recording through the implementation of the multi-analysis surface-wave (MASW) method, and passive surface-wave recording through the deployment of seismic ambient vibration arrays (AVAs).

For the MASW approach, we used a 24-channel acquisition device (Geometrics) connected to either 24 4.5-Hz vertical geophones, or 24 horizontal geophones. We usually performed the acquisition with both vertical polarization (for Rayleigh-wave analysis) and horizontal polarization (for Love-wave analysis). We tried to systematically record one profile with a geophone spaced every 2 m (i.e., a 46 m profile length), but according to the site specificity, this geophone spacing could be reduced to 1 m (i.e., a 23 m profile length). In some sites, several profiles were deployed, with inter-traces of up to 4 m (i.e., a 92 m profile length). For Rayleigh-wave acquisition, an aluminum plate on the ground was hit vertically with a 4-kg hammer, on both sides of the profile and in the middle of the profile. For Love-wave acquisition, the same hammer was used to hit the ground on both sides of the profile using an oak beam weighed down by a car or a few volunteers.

For the AVA acquisition, 30-s broadband seismometers (Güralp CMG6TD) were used for most of the sites. The survey conducted for nine sites in the Pyrenees used 5-s velocimeters (Lennartz LE-3D) connected to a wireless array analysis system (Ohrnberger
<table><caption>**Table 1** Brief description of the geological formations underlying the accelerometric stations studied here</caption><thead><tr><th>Region</th><th>Station</th><th>Geological description</th></tr></thead><tbody><tr><td rowspan="3">Auvergne</td><td>OCLD (OCKE)</td><td>Quaternary deposits overlying Oligocene and volcanic deposits</td></tr><tr><td>OCOL</td><td>Hercynian granite</td></tr><tr><td>OCOR</td><td>Quaternary colluvium overlying anatexites</td></tr><tr><td rowspan="10">Pyrenees</td><td>EPF</td><td>Mesozoic limestones</td></tr><tr><td>PYAS</td><td>Shallow quaternary deposits overlying Mesozoic limestones</td></tr><tr><td>PYAT</td><td>Lower Cretaceous (Albian) marls</td></tr><tr><td>PYBB</td><td>Cretaceous (Albian/Cenomanian) limestones and marls alternations</td></tr><tr><td>PYLI</td><td>Shallow Quaternary alluvium overlying Lower Cretaceous (Aptian) limestone</td></tr><tr><td>PYLL</td><td>Precambrian gneiss</td></tr><tr><td>PYLO</td><td>Lower Cretaceous limestones</td></tr><tr><td>PYLU</td><td>Quaternary alluvial deposits overlying Devonian mudstones</td></tr><tr><td>PYOR</td><td>Quaternary deposits (or weathered formations) overlying magmatic rocks (migmatites)</td></tr><tr><td rowspan="10">Provence and Côte d'Azur</td><td>CALF</td><td>Upper Jurassic (Oxfordian) karstified limestones</td></tr><tr><td>IRPV</td><td>Sand and molassic Miocene deposits overlying Oligocene marls</td></tr><tr><td>IRVG</td><td>Upper Jurassic (Tithonian) limestones</td></tr><tr><td>NALS</td><td>Quaternary alluvium deposits overlying deep Jurassic limestones</td></tr><tr><td>NBOR</td><td>Jurassic (Tithonian and Kimmeridgian) limestones</td></tr><tr><td>OGCA</td><td>Cretaceous (Valanginian) limestone</td></tr><tr><td>GRN</td><td>Upper Jurassic limestone and marls alternations</td></tr><tr><td>OGAN</td><td>Cretaceous (Urgonian) limestone, weathered in surface</td></tr><tr><td>OGAP</td><td>Würmian fluvio-glacial deposits overlying deep Cretaceous limestones</td></tr><tr><td>OGBL</td><td>Recent alluvial and lacustrine deposits overlying Cretaceous limestones</td></tr><tr><td rowspan="10">Alps</td><td>OGCH</td><td>Upper Jurassic limestones</td></tr><tr><td>OGDH</td><td>Post-Würmian alluvial and lacustrine deposits, overlying deep Jurassic limestones</td></tr><tr><td>OGFO</td><td>Post-Würmian alluvial and lacustrine deposits, overlying deep Jurassic limestones</td></tr><tr><td>OGIM</td><td>Recent torrential cones/alluvial deposits, overlying deep Jurassic limestones and marls</td></tr><tr><td>OGLE</td><td>Fractured granites and schists</td></tr><tr><td>OGMA</td><td>Würmian glacial deposits overlying Jurassic limestones and marls</td></tr><tr><td>OGME</td><td>Recent alluvial deposits overlying a terrace, overlying deep Cretaceous limestones</td></tr><tr><td>OGMU</td><td>Upper Jurassic limestone-marl alternations</td></tr><tr><td>OGPC</td><td>Recent alluvial and lacustrine deposits valley overlying deep Jurassic limestones</td></tr><tr><td>OGSI</td><td>Lower Cretaceous limestone</td></tr><tr><td>SURF</td><td>Morainic Würmian deposits overlying Albian-Cenomanian flyschs and black shales</td></tr></tbody></table>
et al. 2006). The first experiments were performed in 2012, and consisted of the successive deployment of several single-circle array geometries at each site, using either 10 sensors (one in the center; nine equally spaced in a circle) or eight sensors (one in the center; seven equally spaced in a circle). The number of circles and the circle diameters were chosen site by site, according to the site conditions, the remaining time during the survey, and other factors. For the acquisitions performed in 2014, a more systematic layout was applied to most sites, due to the larger number of seismological instruments available and the successful initial feedback of the InterPacific project (Foti et al. 2017). Indeed, we usually used 15 sensors, which allowed double-circle geometry (one in the center; 2× seven equally spaced around two circles). The azimuth of the sensors for the larger circle was shifted by ~26° (360°/14) with respect to the azimuth of the smaller circle sensors, to optimize the azimuthal coverage of the whole array. A radius ratio of 3.0 was chosen for the increase between two consecutive array radii, starting from a circle with a 5-m radius, up to a circle with a 405-m radius. For consecutive acquisitions, the inner/smaller circle was moved to form the next larger circle. Hence, in addition to a center station, the four different sets of array layouts used the increasing paired radii of: 5 and 15 m (denoted as array #1); 15 and 45 m (array #2); 45 and 135 m (array #3); and 135 and 405 m (array #4). This standard layout was adjusted to take into account local constraints. For (a priori) rock sites, array #1 was usually omitted (due to lack of energy at high frequency in the seismic ambient noise wavefield), and instead two MASW profiles were performed (the standard 46-m profile length, and a 92-m profile length). For sites that had very difficult access, array #4 was also sometimes omitted. Figure 2 shows the example of the acquisition layout of the OGIM station survey.

Only MASW measurements were performed for four sites (OGLE, OGMU, GRN, OGSI), as their topography was too pronounced to allow AVA acquisition (for both logistic and theoretical considerations). For the OGFO site, we benefited from previously measured AVAs (Bettig et al. 2001; Cornou et al. 2008), and so only the new MASW data were recorded in the framework of the present study. Figure 3 summarizes the different acquisition geometries that were performed at each site.

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Acquisition layout used for the OGIM station measurements (presented at two scales). The ‘ideal’ geometry is adapted to the local logistic constrains (e.g., buildings, roads...). All of the array use the center sensor; ‘array #1’ is composed of stations with radii ‘R1’ and ‘R2’; ‘array #2’ is composed of stations with radii ‘R2’ and ‘R3’; ‘array #3’ is composed of stations with radii ‘R3’ and ‘R4’; ‘array #4’ is composed of stations with radii ‘R4’ and ‘R5’</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Summary of MASW and AVA layouts used for each station studied here. For MASW, ‘V’ indicates vertically polarized acquisition (for further Rayleigh-wave processing); ‘H’ indicates horizontally polarized acquisition (for further Love-wave processing); the number following the *letter* indicates the geophone inter-distance. As all acquisitions used a 24-channel system, ‘V2’ indicates vertically polarized acquisition with 2 m geophone inter-distance, for a 46-m length profile. For AVA, *dotted lines* indicated single-circle arrays, *solid lines* indicate double-circle arrays. The *number* indicates total number of sensors used for further processing (sensors that failed are not counted here). *Gray lines* indicate AVA that failed to produce dispersion curves (whatever the processing method: FK, HRFK or MSPAC). *Other colors* indicate AVA that produced a dispersion curve with at least one of the processing methods. The *different bright colors* are just used for ease of reading</figcaption></figure>
# 4 Processing

## 4.1 Quality monitoring and H/V

The overall processing procedure from the ‘raw data’ to the velocity profiles and the other outcomes was consistent with the current state-of-the-art practices. A description of this procedure can be found in the recommendations of the InterPacific project guidelines (Foti et al. 2017), which focus on fundamental Rayleigh-wave processing and inversion. In the present study, when possible, we also use multi-mode inversion, and Rayleigh-wave and Love-wave joint inversion. The data processing was performed using the Geopsy software (Wathelet et al. 2004, 2008; Wathelet 2005, 2008). The first step of the analysis consisted of examining the Fourier amplitude spectra to estimate the energy content of the seismic ambient noise wavefield, as well as to identify stations that produced spurious estimates, which were then removed from the analysis (e.g., Fig. 4). The data quality monitoring then included evaluation of the classical ambient vibration H/V ratio (Nakamura 1989; SESAME team 2004; Bard et al. 2010) (Fig. 4). As well as providing the fundamental resonance frequency of the site that could be used later in the dispersion curve inversion, the H/V curves provided estimates of the spatial variation of the sediment-to-bedrock depth interface, and allowed checking whether one-dimensional (1D) wave propagation can be assumed throughout the array.

## 4.2 Dispersion curves

The seismic ambient vibrations recorded by the AVA arrays were systematically processed by using: (1) the frequency-wavenumber (FK) (Neidell and Taner 1971; Douze and Laster

<figure><img src="image_7.png" /><figcaption>**Fig. 4** H/V ambient vibration ratio for all of the sensors of the four arrays of the OGIM station measurements. Array #1 and array #2 show good site homogeneity at array scale. Array #3 still shows relatively good site homogeneity, while homogeneity is degraded for array #4</figcaption></figure>
1979; McMechan and Yedlin 1981); (2) the high-resolution frequency-wavenumber (HRFK) (Capon 1969); and (3) the modified spatial autocorrelation (MSPAC) method (Bettig et al. 2001; Köhler et al. 2007). Dispersion curves were then extracted for each array and each processing method. For the FK and HRFK methods, only phase-velocity estimates within the array resolution capabilities were considered. At low frequencies, the minimum resolvable wavenumber was controlled by the array aperture and was easily derived from the array response (Henstridge 1979; Di Giulio et al. 2006; Wathelet et al. 2008). While the array resolution for the FK method was fixed to the width of the main lobe of the array response at its mid-height ($k_{min}$), the HRFK technique usually provided improved resolution capabilities. Although numerous authors (e.g., Asten and Henstridge 1984; Tokimatsu 1997; Satoh et al. 2001) have reported resolution improvements of a factor of three to six, we chose to limit the minimum resolvable wavenumber to $k_{min}/2$ (Cornou et al. 2006; Wathelet et al. 2008), which allowed for the retrieval of larger wavelengths compared to the FK approach. At high frequencies, all of the dispersion curves were picked as high as the aliasing features allowed.

For MSPAC processing, the phase-velocity estimates were derived from the autocorrelation curves (Wathelet 2005). MSPAC based methods are usually reported to resolve wavelengths that are typically five to seven times larger than those predicted by $k_{min}$ (e.g., Horike 1985; Asten et al. 2004).

When FK (or HRFK) and MSPAC provided different dispersion curves within the same frequency band, we monitored the azimuth of the dominant seismic-noise waves. Indeed, the MSPAC approach assumes that the seismic ambient-noise sources are randomly distributed in space. When this assumption was not fulfilled (i.e., the seismic noise wavefield was dominated by a single dominant source direction), we disregarded the MSPAC estimates. Conversely, when MSPAC provided lower phase-velocity estimates than FK-based methods in the low frequency range, we favored the MSPAC estimates, on the assumption of a lack of resolution capability of the FK-based methods.

Not all of the arrays produced usable dispersion curves. The passive experiments that failed to provide phase-velocity estimates were mainly those that corresponded to very ‘quiet’ sites, in connection with a too limited duration of the recording, and a reduced number of stations. Figure 3 (gray) shows the arrays that did not produce usable information in terms of dispersion curves from passive experiments. Globally, over the whole dataset, the HRFK method provided most of the dispersion curves. One example is given in Fig. 5, which shows the FK, HRFK, and MPSAC dispersion images, together with the extracted dispersion curves, for the OGIM station. In this example, all of the processing methods succeeded in providing dispersion estimates, except the MSPAC analysis for the largest array, due to the very strong dominating source between 1 and 4 Hz.

Active MASW measurements were processed for both horizontal and vertical polarization using the FK method. Here, the dispersion curves were also carefully picked, especially at low frequencies, to avoid near-field effects (e.g., cylindric propagation, near-field). For all of the sites, for the shots that were performed with a relatively short distance between the source and the first geophones, the closest geophones from the source were excluded from the processing (typically, 5–10 m far from the source). Next, for the sites where AVA processing led to dispersion curves up to high frequencies (indeed, in most cases), the ‘O’Neil criterion’ (O’Neil 2003) low-frequency limit was used (maximum interpretable wavelength = 0.4 × profile length; e.g., 18.4 m for a line of 24 geophones with a 2-m interstation interval). This criterion was rather constraining, especially for near-source offset; however, using large source offsets helped to relax this rule for larger interpretable wavelengths (O’Neil 2003). For other sites, where only the MASW
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Systematic AVA processing of the OGIM station survey. Left FK, centre HRFK, right MSPAC. Dashed and continuous black lines indicate the theoretical array resolution; namely $k_{min}$ and $k_{min}/2$. All frames have the same scale. All processing allowed the picking of a dispersion curve, except for MSPAC for array #4</figcaption></figure>

measurements were performed or where small aperture AVA failed to produce dispersion curves at high frequency, the dispersion curves beyond the O’Neil limit were picked with care, although wavelengths that were larger than the profile length were never picking, for the sake of safety. Figure 6 shows some examples of the processing and dispersion-curve picking for active MASW measurements performed for the OGIM site.

All of the dispersion curves were then combined to build a broadband dispersion curve for each site. Most of the sites showed several dispersion curves: Rayleigh-wave fundamental mode and higher modes, and Love-wave fundamental mode and higher modes. Note also that in some cases, the dispersion curve corresponding to one given mode was not necessarily continuous from low to high frequencies as frequency gaps can exist. At this step in the whole procedure, the mode attribution of each dispersion curve was indicative and could be further monitored and revised during the inversion tests. Figure 7 shows an example of the individual dispersion curves and combined dispersion curves for

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 6** Example of MASW processing of the OGIM station survey. Left vertically polarized acquisition leading to Rayleigh-wave dispersion curve. Right horizontally polarized acquisition leading to Love-wave dispersion curve. Black solid line indicates wavelength corresponding to the O’Neil criterion (0.4× profile length)</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 7** Gathering of all dispersion curves (DCs) derived from AVA and MASW analysis of the OGIM station survey. Dispersion curves are shown as frequency–velocity plot (left) and wavelength–velocity plot (right). Colored lines shows dispersion curves for each type of processing. Black line shows the final ‘broadband composite’ dispersion curve that is used for inversion. Here, the MSPAC dispersion curve and the lower frequency ($<$40 Hz) segment of the MASW dispersion curve are rejected due to a wavefield dominated by a single source of vibration</figcaption></figure>

the OGIM site. Here, the fundamental mode of Rayleigh waves (denoted as R0) was identified, with no higher modes. The combined R0 dispersion curve for the OGIM station covered a wide frequency range, from 1.3 to 70 Hz. FK and HRFK provided dispersion estimates up to 30 Hz, while MASW allowed the extension of the dispersion curve up to 70 Hz. The dispersion curve derived from MSPAC was rejected here: it was not coherent with the other dispersion curves, and a complementary analysis showed that the wavefield was dominated by a single dominant source of vibration, when MSPAC processing needs a relatively good azimuthal distribution of sources. Figures 8 and 9 show all of the broadband dispersion curves for 30 of the 33 sites, with the phase velocity either expressed as a function of the frequency or as a function of the wavelength. OGBL was the site that was characterized by a dispersion curve that started at the lowest frequency, at 0.4 Hz, inferred from the MSPAC processing on a very wide passive-array aperture (1200 m) leading to the measurement of a 5-km wavelength. At high frequencies, the higher picked frequency (95 Hz) was found for the PYLO site. Several sites (PYOR, IRPV and OGAP) had
<figure><img src="image_12.png" /><figcaption>**Fig. 8** Frequency–phase velocity plots of all of the broadband composite dispersion curves for the 30 sites that led to a successful processing. The *different colors* indicate the dispersion-curve mode (e.g., *red* Rayleigh fundamental mode; *orange* Rayleigh first higher mode; etc.). The *green line* indicates the wavelength of 40 m that is a proxy to roughly infer $V_{S30}$ from the R0 dispersion curve. All frames have the same scale</figcaption></figure>

dispersion curves down to 2 m wavelengths. Statistically, and not surprisingly, sites that were investigated with the most extensive surveys (several 15-sensor double-circle AVA and MASW measurements) had the most continuous dispersion curves over the wider frequency band.

Three sites did not allow derivation of significant dispersion curves. For two of these (EPF, PYLL), AVA did not succeed and the MASW results were very heterogeneous, and only the information at high frequencies (i.e., associated to very shallow layers) was ‘extractable’ from the data. For OGSI, the MASW was indeed performed within a tunnel, which did not allow the generation of a surface wave at a sufficiently large wavelength.
<figure><img src="image_13.png" /><figcaption>Fig. 9 As for legend to Fig. 8, but for the phase velocity-wavelength plots</figcaption></figure>

## 5 Inversion and derivation of velocity profiles

The inversions were performed with the inversion tool of the Geopsy software package, which uses a global search approach with a neighborhood algorithm (Wathelet 2008). The broad-band dispersion curves previously derived were inverted, which in most cases involved joint inversion of both the Rayleigh and Love dispersion curves, with sometimes several various modes. The inversion scheme used here required the identification of surface wave modes. Such identification could be done iteratively starting by the ‘simplest’ solution (fundamental mode for both Rayleigh and Love dispersion curve) and checking whether this assumption provided satisfactory ground profiles in terms of inversion misfit. The Geopsy software package also allows letting free the mode identification and computes different mode attribution possibilities, which helps attributing propagation modes.
Other information than dispersion curves can be used for inversion. For example, for sites where a clear fundamental frequency can be deduced from the H/V analysis, the $f_0$ value was also used in a joint inversion. In general, the usable ‘a-priori’ available geological information was poor, although it could be used in some cases to select between different ground structure parameterization hypotheses. For example, in the case of rock sites that showed visible hard rock outcrops with no evidence of large weathering, the parameterizations that resulted in profiles with shallow layers associated with too low velocity values were rejected in favor of the parameterizations that produced more realistic results. At the OGFO site, the bedrock depth was also set as a priori information, due to the availability of a deep borehole that reached the bedrock at 535 m in depth (Nicoud et al. 2002). For the OGDH and OGFO sites, the availability of nearby geotechnical boreholes also allowed the inclusion of a low-velocity layer in the ground-model parameterization. A low velocity layer was also introduced for the OGBL site. Table 2 summarizes the information that was used for the inversion of each site.

For sites where horizontally polarized MASW acquisition was available and usable, a standard shear-wave seismic refraction analysis was also performed, to provide an alternative evaluation of the ‘first layer’ shear-wave velocity and related thickness. Even if the concept of the first layer is not the same between surface-wave inversion and shear-wave seismic refraction, the order of magnitude of these values were compared and helped to identify sites that showed inconsistencies. These sites were then revised in the inversion procedure. An important additional ‘control’ was included using the phase velocity of the fundamental Rayleigh-wave dispersion curve (when available) at the wavelengths of 40 and 45 m. Some studies have indeed shown that these values provide a good proxy of the $V_{S30}$ (Martin and Diehl 2004; Albarello and Gargani 2010). Here also, when the difference between $V_{S30}$ inferred from the velocity profiles and $V_{\lambda,40}$ or $V_{\lambda,45}$ were too high, the inversion procedures were revised.

In order to limit the depth of provided profiles to account for the real investigation depth, we followed the usual rule of thumb that relates the maximum resolution depth to half the largest wavelength, while also being aware that the actual maximum resolvable depth could be much smaller, especially for sites that show large impedance contrast in the near surface. As for the thickness of the first (shallower) layer, the usual rule of thumb suggests that the smallest thickness should not be smaller than half to a third of the smallest wavelength. While this is suitable when only the fundamental Rayleigh-wave mode is inverted, this criterion can be relaxed with joint inversion of Love-wave mode(s) and Rayleigh-wave mode(s) (Eslick et al. 2008; Xia et al. 2012; Yin et al. 2014).

In order to provide an optimized evaluation of uncertainties associated to velocity profiles inferred from inversion and to $V_{S30}$ values, the following procedures was used. Once the mode numbering determined, different inversion parametrizations were tested (e.g., 2 layers, 3 layers, etc. up to 6 layers, using uniform velocity layers or gradient velocity layers). Only parametrization that produced satisfactory results and misfit were retained. For those parametrization, we used the ‘acceptable misfit’ approach (Lomax and Snieder 1995; Souriau et al. 2011), which led to a set of equivalent $V_S$ profiles that explain the dispersion estimates within their uncertainty bounds. In such a case, the search algorithm does not try to find the best misfit $V_S$ profile, but rather a large number of $V_S$ profiles (50,000 or more) that explain the data equally well. For each site, at least three different parametrizations were used. This way helps in the identification of the $V_S$ profile uncertainties that were associated to the non-uniqueness of the inversion, and the effects of the dispersion-curve uncertainty on the inversion data. Once the inversions done, at least 6000 profiles were randomly extracted from the obtained tens of thousands profiles. These
<table><caption>**Table 2** Summary of information used in the inversion procedures for each site studied here</caption><thead><tr><th>Region</th><th>Station</th><th>Information used in inversion</th></tr></thead><tbody><tr><td rowspan="3">Auvergne</td><td>OCLD (OCKE)</td><td>R0 + L0</td></tr><tr><td>OCOL</td><td>R0 + L0</td></tr><tr><td>OCOR</td><td>R0 + L0</td></tr><tr><td rowspan="8">Pyrenees</td><td>PYAS</td><td>R0 + L0</td></tr><tr><td>PYAT</td><td>R0 + L0</td></tr><tr><td>PYBB</td><td>R0</td></tr><tr><td>PYLI</td><td>R0</td></tr><tr><td>PYLO</td><td>R0 + R1 + L1</td></tr><tr><td>PYLU</td><td>R0 + L1</td></tr><tr><td>PYOR</td><td>R0 + L0</td></tr><tr><td>Provence and Côte d'Azur</td><td>CALF</td><td>R0 + L0</td></tr><tr><td rowspan="7">Alps</td><td>IRPV</td><td>R0 + R1 + L0</td></tr><tr><td>IRVG</td><td>R0</td></tr><tr><td>NALS</td><td>R0 + R1 + L1 + f<sub>0</sub></td></tr><tr><td>NBOR</td><td>R0 + R1</td></tr><tr><td>OGCA</td><td>R0</td></tr><tr><td>GRN</td><td>R0 + R1 + L0</td></tr><tr><td>OGAN</td><td>R0</td></tr><tr><td></td><td></td><td></td></tr><tr><td></td><td>OGAP</td><td>R0 + R1 + L0 + f<sub>0</sub></td></tr><tr><td></td><td>OGBL</td><td>R0 + L1 + f<sub>0</sub> + LVZ</td></tr><tr><td></td><td>OGCH</td><td>R0 + L0</td></tr><tr><td></td><td>OGDH</td><td>R0 + R1 + L1 + L2 + L3 + LVZ</td></tr><tr><td></td><td>OGFO</td><td>R0 + R1 + L0 + H<sub>bed</sub> + LVZ</td></tr><tr><td></td><td>OGIM</td><td>R0 + L1</td></tr><tr><td></td><td>OGLE</td><td>R0 + L0</td></tr><tr><td></td><td>OGMA</td><td>R0 + R1 + L0 + f<sub>0</sub></td></tr><tr><td></td><td>OGME</td><td>R0 + L0 + f<sub>0</sub></td></tr><tr><td></td><td>OGMU</td><td>R0 + R1 + L0</td></tr><tr><td></td><td>OGPC</td><td>R0 + R1 + L1</td></tr><tr><td></td><td>SURF</td><td>R0 + L0</td></tr></tbody></table>

R0 = fundamental mode of Rayleigh-wave dispersion curve; R1 = first higher mode of Rayleigh-wave dispersion curve; R2 = etc.; L0 = fundamental mode of Love-wave dispersion curve; L1 = first higher mode of Love-wave dispersion curve; L2 = etc.; f₀ = fundamental resonance frequency measured with the H/V method (used in inversions when very clear and close to the minimum measured phase velocity). H<sub>bed</sub> = bedrock depth (used when known by boreholes). LVZ = low velocity zone (depths used in inversions when known by boreholes)

profiles were then used to compute average V<sub>S30</sub> value and associated uncertainty, σ. In order to deal with a reasonable set of profiles (for example to compute 1D soil response), we then randomly extracted 33 profiles from the former large set, such as this set of profiles reproduce the V<sub>S30</sub> average and uncertainty σ. We have thus extracted:

* 7 profiles that led to V<sub>S30</sub> values between V̄<sub>S30</sub> – 0.25σ and V̄<sub>S30</sub> + 0.25σ,

* 6 profiles that led to V<sub>S30</sub> values between V̄<sub>S30</sub> – 0.75σ and V̄<sub>S30</sub> – 0.25σ,

* 6 profiles that led to V<sub>S30</sub> values between V̄<sub>S30</sub> + 0.25σ and V̄<sub>S30</sub> + 0.75σ,
* 4 profiles that led to $V_{S30}$ values between $\bar{V}_{S30} - 1.25\sigma$ and $\bar{V}_{S30} - 0.75\sigma$,

* 4 profiles that led to $V_{S30}$ values between $\bar{V}_{S30} + 0.75\sigma$ and $\bar{V}_{S30} + 1.25\sigma$,

* 2 profiles that led to $V_{S30}$ values between $\bar{V}_{S30} - 1.25\sigma$ and $\bar{V}_{S30} - 1.75\sigma$,

* 2 profiles that led to $V_{S30}$ values between by $\bar{V}_{S30} + 1.25\sigma$ and $\bar{V}_{S30} + 1.75\sigma$,

* 1 profile that led to $V_{S30}$ value between by $\bar{V}_{S30} - 1.75\sigma$ and $\bar{V}_{S30} - 2.25\sigma$,

* 1 profile that led to $V_{S30}$ value between by $\bar{V}_{S30} + 1.75\sigma$ and $\bar{V}_{S30} + 2.25\sigma$,

$\bar{V}_{S30}$ being the average $V_{S30}$ value and $\sigma$ being the standard deviation.

Figure 10 shows the inversion results (subset of the ‘acceptable’ $V_S$ profiles) for the OGIM site by using a 3-layers-over-half-space parametrization. To visualize the robustness of these results, the forward-modeled dispersion curves were compared to the measured dispersion curves (i.e., Rayleigh and Love, fundamental and first higher modes). For the OGIM site, the most consistent results were obtained by attributing the fundamental mode (R0) to the Rayleigh dispersion curve, while the first higher mode was chosen for the Love dispersion curve (L1). Figure 11 shows the $V_S$ profiles derived for all of the 30 stations (here, for each station, only the 33 selected $V_S$ profiles are represented, together with the envelope of all ‘acceptable’ $V_S$ profiles).

The derived $V_S$ profiles have to be considered as profiles that satisfyingly explain the experimental dispersion curves(s). The $V_S$ profiles do not necessary reflect the exact lithological layering. For example, the most general ‘staircase shape’ of the ground profiles may reflect the way the inversion accommodated a ‘gradient-like’ real $V_S$ profile (i.e., continuous increase in $V_S$ with depth).

When the seismic bedrock velocity was well constrained, we could also compute the local 1D site responses, even though the $V_S$ profiles themselves were non-unique as a result of the inversion process, or were ‘simplified’ by the inversion ground-model parameterization choices. The determination of the bedrock velocity could appear rather difficult, but nonetheless achievable as shown in this study. This estimation was made possible in this study in most cases through the use of rather large arrays. Until recently, it was considered that passive surface wave-based methods were not suitable for rock sites, especially due to the low level of ambient vibration characterizing such sites. However, as far as careful signal analysis is performed (seismic noise level, instrument setup, recording durations), consistent results can be obtained. Another example of hard rock site characterization could be found in Garofalo et al. (2016a, b).

<figure><img src="image_14.png" /><figcaption>**Fig. 10** Inversion results for the OGIM station survey for one given parametrization. Left set of statistically acceptable $V_S$ profiles explaining the observed dispersion data within their uncertainty bounds. Center and right theoretical dispersion curves inferred from the full inverted profile set for Rayleigh phase (center) and Love phase (right), compared to the measured dispersion curves and estimated standard deviation</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Fig. 11** ‘Acceptable’ $V_s$ profiles (blue lines) obtained by inversion for the 30 sites that led to successful processing. For each site, a set of 33 profiles are shown, randomly chosen to explore the uncertainty within the full set of ‘acceptable’ profiles. Red lines show the envelopes of all computed ‘acceptable’ profiles. All frames have the same scale, but note the unusual log-scale on the depth axis, used here to show all of the stations results in the same way</figcaption></figure>

# 6 Contribution of Love-wave active measurements

In most common surveys, only Rayleigh waves are recorded and considered for inversion. In the present study, for the MASW acquisition, we recorded both Rayleigh waves (with vertical geophones and excitation) and Love waves (with horizontal geophones and excitation). In several cases, and especially on rock sites, the Love-wave MASW acquisition allowed us to obtain clearer dispersion curves than the Rayleigh-wave MASW acquisition. This observation has also been outlined through numerical and experimental studies (Safani et al. 2005; Dal Moro and Ferigo 2011; Xia et al. 2012; Boaga et al. 2013).
<figure><img src="image_16.png" /><figcaption>**Fig. 12** Example of the contribution of Love-wave dispersion curves (DC) in addition to Rayleigh-wave dispersion curves within the inversion, for the OGCH station survey. On the top result of the inversion when only the Rayleigh dispersion curve is used to constrain the inversion (the empirical Love dispersion curve is shown in blue but is not used for inversion). On the bottom both Rayleigh and Love dispersion curves are used, leading to a more constrained velocity profile (the empirical Love dispersion curve is shown in black, all direct dispersion curves computed with inverted profiles fit now the empirical Love dispersion curve). Left ensemble of statistically acceptable V~S~ profiles. Center and right theoretical dispersion curves inferred from the full inverted profile set for Rayleigh phase (center) and Love phase (right), compared to the measured dispersion curves and estimated standard deviations</figcaption></figure>

Measuring Love waves thus greatly helps in the identification of the modes of propagation, especially at those sites that show weathered rock material at the near surface overlying much stiffer geological units, or at sites with low-velocity zones.

When possible, we jointly inverted Rayleigh and Love wave dispersion curves, thus allowing for better-constrained inversions and a more robust set of V~S~ profiles. Figure 12 shows the example of the OGCA site with only Rayleigh R0 mode inversion, and with both Rayleigh and Love inversion (R0 and L0 modes). For this site, the use of the Love dispersion curve allowed the V~S~ model to be better constrained within the first meters below the surface, and especially for the weathered near-surface layers. This information in the near surface is also important to understand the high frequency content of the accelerograms recorded at these stations.

# 7 V$_{S30}$ and soil class results

Table 3 gathers all the results in terms of V~S30~ and the associated uncertainties and EC8 soil classes. Most of the V~S30~ values in Table 3 (i.e., 30 out of 33) were computed from velocity profiles derived from surface-wave dispersion-curve inversion. These sites
<table><caption>**Table 3** Summary of the V<sub>S30</sub> and soil class estimations (medians and standard deviation) for each studied site</caption><thead><tr><th>Region</th><th>Station</th><th>Site V<sub>S30</sub> (m/s)</th><th>Station V<sub>S30</sub> (m/s)</th><th>σ (m/s)</th><th>COV (%)</th><th>Station EC8 class</th><th>Soil class (preferred)</th><th>V<sub>S30</sub> source</th><th>σ source</th></tr></thead><tbody><tr><td rowspan="2">Auvergne</td><td>OCLD</td><td>406</td><td>406</td><td>52</td><td>13</td><td>B or E</td><td>B</td><td>1</td><td>4</td></tr><tr><td>OCOL</td><td>≥1184</td><td>≥1484</td><td>297</td><td>20</td><td>A</td><td>A</td><td>1</td><td>3</td></tr><tr><td rowspan="10">Pyrenees</td><td>OCOR</td><td>1072</td><td>1072</td><td>214</td><td>20</td><td>A</td><td>A</td><td>1</td><td>3</td></tr><tr><td>EPF</td><td>1100</td><td>1100</td><td>400</td><td>36</td><td>likely A, maybe B</td><td>A</td><td>2</td><td>5</td></tr><tr><td>PYAS</td><td>1089</td><td>1089</td><td>425</td><td>39</td><td>likely A, maybe B</td><td>A</td><td>1</td><td>4</td></tr><tr><td>PYAT</td><td>1164</td><td>1164</td><td>233</td><td>20</td><td>A</td><td>A</td><td>1</td><td>3</td></tr><tr><td>PYBB</td><td>894</td><td>894</td><td>99</td><td>11</td><td>likely A, maybe B</td><td>A</td><td>1</td><td>4</td></tr><tr><td>PYLI</td><td>1257</td><td>1257</td><td>251</td><td>20</td><td>A</td><td>A</td><td>1</td><td>3</td></tr><tr><td>PYLL</td><td>1100</td><td>1100</td><td>400</td><td>36</td><td>likely A, maybe B</td><td>A</td><td>2</td><td>5</td></tr><tr><td>PYLO</td><td>964</td><td>1639</td><td>328</td><td>20</td><td>A</td><td>A</td><td>1</td><td>3</td></tr><tr><td>PYLU</td><td>388</td><td>388</td><td>39</td><td>10</td><td>likely B, maybe C</td><td>B</td><td>1</td><td>2</td></tr><tr><td>PYOR</td><td>426</td><td>426</td><td>85</td><td>20</td><td>likely B, maybe C</td><td>B</td><td>1</td><td>3</td></tr><tr><td rowspan="5">Provence and Côte d'Azur</td><td>CALF</td><td>1746</td><td>1746</td><td>349</td><td>20</td><td>A</td><td>A</td><td>1</td><td>3</td></tr><tr><td>IRPV</td><td>516</td><td>605</td><td>61</td><td>10</td><td>B</td><td>B</td><td>1</td><td>2</td></tr><tr><td>IRVG</td><td>2090</td><td>2090</td><td>209</td><td>10</td><td>A</td><td>A</td><td>1</td><td>2</td></tr><tr><td>NALS</td><td>171</td><td>171</td><td>17</td><td>10</td><td>C or D</td><td>D</td><td>1</td><td>2</td></tr><tr><td>NBOR</td><td>927</td><td>927</td><td>202</td><td>22</td><td>likely A, maybe B</td><td>A</td><td>1</td><td>4</td></tr><tr><td rowspan="7">Alps</td><td>OGCA</td><td>1383</td><td>1383</td><td>138</td><td>10</td><td>A</td><td>A</td><td>1</td><td>2</td></tr><tr><td>GRN</td><td>819</td><td>1062</td><td>106</td><td>10</td><td>A</td><td>A</td><td>1</td><td>2</td></tr><tr><td>OGAN</td><td>1003</td><td>1003</td><td>143</td><td>14</td><td>A</td><td>A</td><td>1</td><td>4</td></tr><tr><td>OGAP</td><td>240</td><td>240</td><td>24</td><td>10</td><td>C</td><td>C</td><td>1</td><td>2</td></tr><tr><td>OGBL</td><td>196</td><td>196</td><td>30</td><td>15</td><td>likely C, maybe D</td><td>C</td><td>1</td><td>4</td></tr><tr><td>OGCH</td><td>1454</td><td>1454</td><td>211</td><td>14</td><td>A</td><td>A</td><td>1</td><td>4</td></tr><tr><td>OGDH</td><td>194</td><td>194</td><td>19</td><td>10</td><td>likely C, maybe D</td><td>C</td><td>1</td><td>2</td></tr><tr><td></td><td>OGEO</td><td>259</td><td>259</td><td>26</td><td>10</td><td>C</td><td>C</td><td>1</td><td>2</td></tr><tr><td></td><td>OGIM</td><td>544</td><td>599</td><td>64</td><td>11</td><td>B</td><td>B</td><td>1</td><td>4</td></tr></tbody></table>
<table><caption>**Table 3** continued</caption><thead><tr><th>Region</th><th>Station</th><th>Site V<sub>S30</sub> (m/s)</th><th>Station V<sub>S30</sub> (m/s)</th><th>σ (m/s)</th><th>COV (%)</th><th>Station EC8 class</th><th>Soil class (preferred)</th><th>V<sub>S30</sub> source</th><th>σ source</th></tr></thead><tbody><tr><td></td><td>OGLE</td><td>≥763</td><td>≥763</td><td>162</td><td>21</td><td>A or B</td><td>B</td><td>1</td><td>4</td></tr><tr><td></td><td>OGMA</td><td>845</td><td>845</td><td>93</td><td>11</td><td>likely A, maybe B</td><td>A</td><td>1</td><td>4</td></tr><tr><td></td><td>OGME</td><td>389</td><td>389</td><td>52</td><td>13</td><td>likely B, maybe C</td><td>B</td><td>1</td><td>4</td></tr><tr><td></td><td>OGMI</td><td>≥913</td><td>≥1338</td><td>268</td><td>20</td><td>A</td><td>A</td><td>1</td><td>1</td></tr><tr><td></td><td>OGPC</td><td>507</td><td>620</td><td>62</td><td>10</td><td>B</td><td>B</td><td>1</td><td>2</td></tr><tr><td></td><td>OCSI</td><td></td><td>1500</td><td>500</td><td>33</td><td>A</td><td>A</td><td>2</td><td>5</td></tr><tr><td></td><td>SURF</td><td>411</td><td>518</td><td>59</td><td>11</td><td>B</td><td>B</td><td>1</td><td>4</td></tr></tbody></table>

‘Site V<sub>S30</sub>’ gives the values computed from the surface, inferred directly from the survey; ‘Station V<sub>S30</sub>’ gives the values computed for the accelerometric station, taking into account its depth when not on the surface. V<sub>S30</sub> source: 1 = values inferred from surface-wave dispersion-curve inversions (and possibly other information: f<sub>0</sub>, bedrock depth...); 2 = values roughly estimated with refraction seismic information and in geological analogy with other sites. σ source: 2 = minimum postulated uncertainty of 10% (when all other uncertainties sources are <10%, including R0 dispersion-curve uncertainty at λ = 40 m); 3 = minimum postulated uncertainties of 20% (when all other uncertainties sources are <20%, but the R0 dispersion curve is not available at λ = 40 m); 4 = uncertainties estimated with the R0 dispersion-curve uncertainty at λ = 40 m; 5 = very large postulated uncertainties when V<sub>S30</sub> source = 2. COV = coefficient of variation

<figure><img src="image_17.png" /></figure>
correspond to the indication of ‘1’ in the ‘$V_{S30}$ source’ column of Table 3. For three of the 33 sites (EPF, PYLL, OGSI), it was not possible to build up a broadband dispersion curve, and $V_S$ profile inversion was thus not possible (or simply for the few meters of the uppermost layers). In this case, $V_{S30}$ was roughly evaluated on the basis of the velocity within the shallowest layers determined from the $V_S$ seismic refraction analysis and the phase velocity at very high frequencies, and by analogy with the other sites with the same velocity within the shallower layers. These sites correspond to the indication ‘2’ in the ‘$V_{S30}$ source’ column of Table 3. For three other sites (OCOL, OGLE, OGMU), the $V_S$ profiles do not reach the depth of 30 m as far as we respect the rule of thumb that limits the maximum resolution depth to half the largest wavelength of the measured dispersion curve. For these sites, we used the deepest velocity value of the calculated velocity profile (which is for all these cases rather high and corresponds to unweathered rock) and used it to complete the profile up to 30 m. From a general point of view, as the velocity profile usually increases with depth, these extrapolated values may underestimate the real $V_{S30}$ value.

We attempted to estimate the uncertainties in three ways. The first was from the inversion using the ‘acceptable misfit’ concept applied on multi-parametrization results. This option probably underestimates the real uncertainty. The second option consisted of the use of the Rayleigh experimental dispersion-curve standard deviations at the wavelength of 40 m (i.e., when the dispersion curve was ‘available’ at this wavelength, and when the dispersion curve corresponded to the fundamental Rayleigh mode). Indeed, some studies have shown that this value is a good proxy of the $V_{S30}$ (Martin and Diehl 2004; Albarello and Gargani 2010). For some sites where the dispersion curves were not strictly available at 40 m, we slightly extrapolated (or interpolated) dispersion curves to get the standard deviation at 40 m. The third option was to apply a ‘minimum’ fixed coefficient of variation based on the InterPacific project feedback (Garofalo et al. 2016b) for the variability between the various analysts on the ‘best’ $V_S$ estimates. This minimum coefficient of variation was set at 10% when the R0 dispersion curve was available at $\lambda = 40$ m, and was set at 20% when it was not available (this was often the case when only the MASW experiments were usable). Obviously, these ‘uncertainty estimates’ were not comparable to one another. However, due to a lack of scientific consensus and any well-established method to quantify the uncertainty of inverted $V_S$ profiles (Cornou et al. 2006; Di Giulio et al. 2012; Michel et al. 2014; Cox and Teague 2016), we finally chose the ‘worst’ uncertainty estimation. In Table 3, the ‘$\sigma$ source’ column gives the chosen uncertainty estimates:

* 1: based on the acceptable misfit, mutli-parametrization approach;

* 2: based on a fixed minimum 10% uncertainty;

* 3: based on a fixed minimum 20% uncertainty;

* 4: based on the dispersion-curve uncertainty at $\lambda = 40$ m;

* 5: large uncertainty in the case of $V_{S30}$ roughly estimated for the $V_S$ refraction data.

In Table 3, the ‘site $V_{S30}$’ and ‘station $V_{S30}$’ were distinguished. The site $V_{S30}$ corresponds to values computed from the surface to 30 m, inferred directly from the survey. In some cases, the accelerometric station was not installed at the surface, but at a certain depth (e.g., when the station was installed in the basement of a building). In this case, a second value was computed, from $h$ to $h + 30$ m, where $h$ is the sensor depth. This value corresponds to the ‘station $V_{S30}$’ value. In our study, the maximum station depth is 5 m. When the ‘station $V_{S30}$’ values are higher than the ‘site $V_{S30}$’ values, this means that the shallower (and often softer) layers are not accounted for in the calculation. In Table 3, the
uncertainties values are associated to the “station $V_{S30}$”. When this uncertainty comes from the experimental dispersion-curve standard deviation at the wavelength of 40 m, no distinction is made between site and station $V_{S30}$ since there are no significant differences between 40 and 45 m on the dispersion curves in terms of in standard deviation.

A comparison between the present study and the previous studies can be carried out. For example, Régnier et al. (2010) gathered a wide range of pre-existing geophysical and geotechnical information in the vicinity of the RAP sites. However, the sources of this information were very heterogeneous and the soil conditions of numerous stations were inferred from very indirect inputs, such as geological maps. The comparison between the present work and the previous one by Régnier et al. (2010) can be carried out for 28 sites. Among these stations:

* 3 previously assumed ‘A class’ stations now ‘downgraded’ to B class;

* 1 previously assumed ‘B class’ station now ‘downgraded’ to C class;

* 2 previously assumed ‘C class’ stations now ‘downgraded’ to D class.

* 1 previously assumed ‘E class’ station may be now ‘upgraded’ to B class (but an uncertainty remains for the OCLD site).

In terms of $V_{S30}$ values:

* 9 stations now have lower $V_{S30}$;

* 15 stations remain unchanged (within ±25%).

* 4 stations now have higher $V_{S30}$.

Another previous study was done by Drouet et al. (2010) that have been aimed at estimation of $V_{S30}$ using generalized inversion methods and by comparing site-transfer functions from the inversion to generic amplification curves derived from generic $V_s$ profiles (Boore and Joyner 1997; Cotton et al. 2006). This study was focused on 21 presumed rock sites. Here, the comparison can be made for 14 stations. In terms of $V_{S30}$ values:

* 5 stations have now lower $V_{S30}$;

* 7 stations remain unchanged (within ±25%);

* 2 stations now have higher $V_{S30}$.

Of course, these previous studies did not use extensive in situ measurements, and the station soil classes and $V_{S30}$ were often estimated indirectly. However, this overall observation shows that $V_{S30}$ of the network stations is usually and statistically overestimated when classification is not based on in situ measurements (this has already been observed for other networks, such as in Italy, see e.g., Pileggi et al. 2011).

## 8 Reference stations

Beyond the discussion of the $V_{S30}$ values, one important issue is the concept of the reference stations. These reference stations can be defined as stations that were not affected by local amplification due to site effects. To comment on this concept, we chose the stations from our set that were clearly within the A soil class, and for which it was possible to derive a shear-wave velocity profile (OCOL, OCOR, PYAT, PYLI, PYLO, CALF, IRVG, OGCA, GRN, OGAN, OGCH, OGMU). For each of these, the theoretical 1D transfer function (1DTF) for vertically incident S-waves was then computed, using the Haskell-Thomson computation. All 33 representative profiles were used; mean and
<figure><img src="image_18.png" /><figcaption>**Fig. 13** Example of 1D transfer functions (1DTF) computed using the shear-wave velocity profiles inferred from surface-wave inversion. We chose stations that are clearly within the EC8-A soil class and for which the ambient vibration H/V ratio frequency peaks are coherent with 1DTFs. For each station, 33 1DTF were computed using 33 different ‘acceptable’ profiles (grey lines), as well as the means and standard deviation of these (red line). All stations show high-frequency amplification due to shallow ‘softer’ layers. The frequency identified by the green vertical line is the one above which amplification >1.5</figcaption></figure>

standard deviation were computed for each site. These curves are divided by 2 to correct for the free surface effects (i.e., the ‘1’ value means no amplification). The simulations were carried out with damping defined as $Q_s = V_s/20$. As an example, the amplification of three of these sites (OGAN, OGCH, PYLI), are given in Fig. 13 (these stations also have an ambient vibration H/V ratios characterize by a peak frequency that was coherent with the 1DTF). Most of the stations showed important amplification at high frequency. This amplification can have an impact on the evaluation of the ‘high-frequency parameters’, such as $κ$ (Ktenidou et al. 2014). It would then be interesting to define a parameter that can determine for a given station the frequency up to which the station site response function can be assumed to remain ‘flat’. For example, we identify for each station the frequency above which the amplification is >1.5. This frequency is relatively low for some of the stations (e.g., OGAN: 6.6 Hz), whereas it is high for others (e.g., OGCH: 14.6 Hz).

## 9 Topographic amplification

Although this study is not intended to provide a systematic assessment of the topographic site effect, the magnitude of this effect could be assessed in relation to the lithologic effect. In order to perform this evaluation, the ‘frequency-scaled curvature’ (FSC) proxy proposed by Maufroy et al. (2015) was computed using available digital terrain models. This approach was applied on four selected stations (IRVG, OGCA, GRN and IRPV). IRVG and GRN are interesting for this study since they are likely the stations that are located on the steepest relief among our station set. IRVG is also the station that has the highest $V_{S30}$ value (2090 m/s) and the only one that shows almost no amplification (even at high frequency) according to the previous lithological amplification analysis. In contrast, OGCA is a station that has a high $V_{S30}$ (1383 m/s), while being located on a gentle slope. IRPV has a lower $V_{S30}$ (605 m/s) and is located on a sand and molassic Miocene deposit basin. IRVG, OGCA and IRPV FSC proxies were computed using digital terrain models (DTM) with a high resolution of 1 m, whereas GRN was computed with a 25 m DTM. FSC proxy is expressed as a function of shear wave wavelength. In order to convert as function of frequency, the velocity of 3000 m/s was used for all of the four considered stations. This value corresponds to the order of magnitude of velocity at depth for all stations, the one to be considered for the large wavelengths impacted by topographic effect. FSC proxies were
<figure><img src="image_19.png" /><figcaption>**Fig. 14** Estimation of lithological amplification and topographic amplification on four sites (IRVG, OGCA, GRN and IRPV). Lithological amplification is computed with 1D numerical simulation using the 33 shear wave velocity profiles for each site obtained in the present paper. Topographic amplification is predicted using the frequency-scaled curvature proxy proposed by Maufroy et al. (2015) using digital terrain models (with a resolution of 1 m for IRVG, OGCA and IRPV; and 25 m for GRN—causing the evaluation limitation at high frequency)</figcaption></figure>

computed down to the minimum wavelength of 100 m for IRVG, OGCA and IRPV, and down to 300 m for GRN. Results are shown on Fig. 14 where topographic amplification (computed with Maufroy et al. 2015, FSC proxy) and lithological 1D amplification (computed following the procedure described above) are compared.

Among the tested sites, two of them show no topographic amplification, i.e., OGCA and IRPV. While IRPV is located close to a small ridge, its closeness to that topographic feature is not enough for that station to be affected. The station itself is located over a very weak and regular slope that doesn’t produce amplification according to the FSC proxy. Similarly, OGCA is located over a gentle slope and no topographic amplification is awaited. Only a weak deamplification occurs at OGCA below 3 Hz, but that effect is too weak to be observable. Given such results, any amplification effect observed at OGCA or IRPV should be attributed to lithological origins. On the contrary, IRVG and GRN are both located over sharp topographies, and thus the FSC proxy predicts a noticeable topographic amplification at both sites. IRVG appears affected over a large frequency bandwidth, with the level of amplification being remarkably stable among frequencies, decreasing only towards 1 Hz. GRN is even more affected, but the low-resolution DTM used for the proxy calculation doesn’t allow observing if the effect remains in the higher frequencies.

This figure shows that the topographic amplification is a broadband smooth function, compared to the lithological effect that produces strong peaks of amplification in narrow bandwidths. The topographic site effect may affect rock sites typically used as reference sites, which underlines the importance to carefully analyse the amplitude of the ground
motion at those sites even if no clear peak of amplification appears, keeping in mind that the topography may smoothly amplify the ground motion over a large bandwidth.

One remaining issue is the combination of both lithological and topographic site effects at the same site. Predictions are available for separated effects, but there is currently no method to conjointly predict both effects, or on how to combine both predictions. As illustrated by our results at IRVG and GRN, both effects can overlap in the same frequency band, and thus this issue should be tackled in order to better understand the origins of the amplification measured at topographic sites.

Note also the strong deamplification around 25 Hz for the GRN station. In fact, the station is not located at surface, but in a small cave, at a depth of around 5 m. This deamplification is due to the destructive interference of the downgoing wave. This situation should also be taken into account for high frequency analysis of accelerometer records since it can also widely affect signal spectrum.

# 10 Conclusions

We carried out site characterization surveys on 33 sites of the French RAP network. This campaign was performed using surface-wave-based methods that allowed the derivation of the shear-wave velocity profiles from the dispersion characteristics of the surface waves. Even if the logistic context was sometimes difficult, the use of surface-wave-based methods is suitable for accelerometric station characterization, and also on rock sites (where the applicability of these methods has sometimes been disputed). It was possible to retrieve broadband dispersion curves for 30 sites to derive shear-wave velocity profiles, and then to propose $V_{S30}$ values and EC8 soil classes, with their associated uncertainties. For three sites, the EC8 soil classes could be deduced from S-wave refraction seismic processing, although the surface-wave-based processing failed to produce significant results, probably due to a too-low level of ambient vibration and/or a too-short duration of recording and/or a too-low number of sensors. To maximize the chances of obtaining usable data, we recommend the use of a relatively high number of sensors for AVA (15 sensors in this study) with the proposed double-circle geometry and with successive increasing apertures. We also strongly recommend the joint use of active (MASW) and passive (AVA) methods, to obtain ‘broadband’ dispersion curves. The use of Love waves (for MASW in this case) is also very valuable in the inversion process, to reduce the probability of misidentification of surface-wave modes and to increase the inverted $V_S$ profile resolution within the shallower layers.

In comparison with previous studies, which have mainly estimated $V_{S30}$ indirectly, these new values are globally lower, and a few of the stations were ‘downgraded’ in terms of the EC8 soil classes, although there remain many EC8-A class sites. Interestingly here, the seismic bedrock velocity could be measured at some rock sites due to the large-aperture arrays. Some high frequency amplification is suggested at most of the rock sites due to the near-surface weathering. This feature should definitely be accounted for in terms of the good use of accelerometric data in the case of high-frequency studies.

As mentioned by Michel et al. (2014), future comparisons of predicted 1D seismic responses with those obtained on earthquake recordings, as well as inclusion of other dispersion features (e.g., polarization, ellipticity, use of Love-wave dispersion curves from passive recordings, use of various ground model parameterization) will definitely help in
the constraining of $V_s$ profiles and the related uncertainties, to reach robust and very reliable site characterization, and hereafter, better estimates of local site effects.

## 11 Data and resources

The data from field measurement can be obtained by simple request to the authors. The results presented in this paper will be available on the RESIF webportal (http://www.resif.fr/).

**Acknowledgements** This work was conducted within the framework of the Cadarache Seismic Hazard Integrated Multidisciplinary Assessment (CASHIMA; funded by CEA, ILL and ITER Organization) and Seismic Ground-Motion Assessment (SIGMA; funded by EdF, CEA, Areva and Enel) research programs. This work is also supported by EU-Horizon2020 project EPOS-IP, Grant Number 676564. We would like to thank all of the contributors to this work (especially for help in the field): Alice Renault (Soldata), Anne Loevenbruck (CEA/LDG), Anthony Bonjour (BRGM), Aurore Laurendeau (CEA/LDG), Catherine Pequegnat (ISTerre), Cécile Cretin (ISTerre), Claire Labonne (CEA/LDG), Corinne Lacave (Résonance), David Wolyniec (ISTerre), Didier Brunel (GeoAzur), Diego Mercerat (CEREMA), Emmanuel Chaljub (ISTerre), Etienne Bertrand (CEREMA), Franck Grimaud (OMP), Franck Tilloloy (BRGM), Jean Letort (ISTerre), Jean-Pierre Deslandes (ISTerre), Mathieu Causse (ISTerre), Michel Pernoud (CEREMA), Paul Calou (CEA/LDG), Robin Barbier (BRGM), Roser Hoste-Colomber (CEA/LDG), Sadrac Saint-Fleur (GeoAzur), Stéphane Netchschein (IRSN). C. Berrie provided valuable improvement of the manuscript. We also thank two anonymous reviewers for the constructive comments which greatly helped us to improve our manuscript.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## References

Akkar S, Sandıkkaya MA, Şenyurt M, Azari Sisi A, Ay BÖ, Traversa P, Douglas J, Cotton F, Luzi L, Hernandez B, Godey S (2014) Reference database for seismic ground-motion in Europe (RESOURCE). Bull Earthq Eng 12:311–339. doi:10.1007/s10518-013-9506-8

Albarello D, Gargani G (2010) Providing NEHRP soil classification from the direct interpretation of effective Rayleigh-wave dispersion curves. Bull Seismol Soc Am 100:3284–3294. doi:10.1785/0120100052

Allen TI, Wald DJ (2009) On the use of high-resolution topographic data as a proxy for seismic site conditions (VS30). Bull Seismol Soc Am 99:935–943. doi:10.1785/0120080255

Asten M, Henstridge J (1984) Array estimators and the use of microseisms for reconnaissance of sedimentary basins. Geophysics 49:1828–1837. doi:10.1190/1.1441596

Asten M, Dhu T, Lam N (2004) Optimised array design for microtremor array studies applied to site classification; observations, results and future use. In: Proceedings of the 13th world conference of earthquake engineering, Vancouver

Bard P-Y, Cadet H, Endrun B, Hobiger M, Renalier F, Theodulidis N, Ohrnberger M, Fäh D, Sabetta F, Teves-Costa P, Duval A-M, Cornou C, Guillier B, Wathelet M, Savvaidis A, Köhler A, Burjanek J, Poggi V, Gassner-Stamm G, Havenith HB, Hailemikael S, Almeida J, Rodrigues I, Veludo I, Lacave C, Thomassin S, Kristekova M (2010) From non-invasive site characterization to site amplification: recent advances in the use of ambient vibration measurements. In: Garevski M, Ansal A (eds) Earthquake engineering in Europe. Springer, Dordrecht, pp 105–123

Bettig B, Bard PY, Scherbaum F, Riepl J, Cotton F, Cornou C, Hatzfeld D (2001) Analysis of dense array noise measurements using the modified spatial auto-correlation method (SPAC): application to the Grenoble area. Boll Geofis Teor Ed Appl 42:281–304
Boaga J, Cassiani G, Strobbia CL, Vignoli G (2013) Mode misidentification in Rayleigh waves: ellipticity as a cause and a cure. Geophysics 78:EN17–EN28. doi:10.1190/geo2012-0194.1

Boore DM, Joyner WB (1997) Site amplifications for generic rock sites. Bull Seismol Soc Am 87:327–341

Capon J (1969) High-resolution frequency-wavenumber spectrum analysis. Proc IEEE 57:1408–1418

Comina C, Foti S (2015) Discussion on “Implications of surface wave data measurement uncertainty on seismic ground response analysis” by Jakka et al. Soil Dyn Earthq Eng 74:89–91. doi:10.1016/j.soildyn.2014.10.001

Cornou C, Ohrnberger M, Boore DM, Kudo K, Bard P-Y (2006) Derivation of structural models from ambient vibration array recordings: Results from an international blind test. In: Proceedings of the third international symposium on the effects of surface geology on seismic motion, Grenoble, France, pp 1127–1219

Cornou C, Chaljub E, Tsuno S (2008) Real and synthetic ambient noise recordings in the Grenoble basin: reliability of the 3D numerical model. In: Proceedings of 14th world conference on earthquake engineering. Paper 03-03-0070, Beijing, China

Cotton F, Scherbaum F, Bommer JJ, Bungum H (2006) Criteria for selecting and adjusting ground-motion models for specific target regions: application to central Europe and rock sites. J Seismol 10:137–156. doi:10.1007/s10950-005-9006-7

Cox BR, Teague DP (2016) Layering ratios: a systematic approach to the inversion of surface wave data in the absence of a priori information. Geophys J Int 207:422–438. doi:10.1093/gji/ggw282

Dal Moro G, Ferigo F (2011) Joint analysis of Rayleigh- and Love-wave dispersion: issues, criteria and improvements. J Appl Geophys 75:573–589. doi:10.1016/j.jappgeo.2011.09.008

Di Giulio G, Cornou C, Ohrnberger M, Wathelet M, Rovelli A (2006) Deriving wavefield characteristics and shear-velocity profiles from two-dimensional small-aperture arrays analysis of ambient vibrations in a small-size alluvial basin, Colfiorito, Italy. Bull Seismol Soc Am 96:1915–1933. doi:10.1785/0120060119

Di Giulio G, Savvaidis A, Ohrnberger M, Wathelet M, Cornou C, Knapmeyer-Endrun B, Renalier F, Theodoulidis N, Bard P-Y (2012) Exploring the model space and ranking a best class of models in surface-wave dispersion inversion: application at European strong-motion sites. Geophysics 77:B147–B166. doi:10.1190/geo2011-0116.1

Douze EJ, Laster SJ (1979) Statistics of semblance. Geophysics 44:1999–2003. doi:10.1190/1.1440953

Drouet S, Chevrot S, Cotton F, Souriau A (2008) Simultaneous inversion of source spectra, attenuation parameters, and site responses: application to the data of the french accelerometric network. Bull Seismol Soc Am 98:198–219. doi:10.1785/0120060215

Drouet S, Cotton F, Guéguen P (2010) Vs30, κ, regional attenuation and Mw from accelerograms: application to magnitude 3-5 French earthquakes. Geophys J Int 182:880–898. doi:10.1111/j.1365-246X.2010.04626.x

Eslick R, Tsoflias G, Steeples D (2008) Field investigation of Love waves in near-surface seismology. Geophysics 73:G1–G6. doi:10.1190/1.2901215

EUROCODE8 (1998) Design of structures for earthquake resistance—Part 1: general rules, seismic actions and rules for buildings. European Committee for Standardization, Brussels

Foti S, Comina C, Boiero D, Socco LV (2009) Non-uniqueness in surface-wave inversion and consequences on seismic site response analyses. Soil Dyn Earthq Eng 29:982–993. doi:10.1016/j.soildyn.2008.11.004

Foti S, Parolai S, Bergamo P, Di Giulio G, Maraschini M, Milana G, Picozzi M, Puglia R (2011) Surface wave surveys for seismic site characterization of accelerometric stations in ITACA. Bull Earthq Eng 9:1797–1820. doi:10.1007/s10518-011-9306-y

Foti S, Hollender F, Garofalo F, Albarello D, Asten M, Bard PY, Comina C, Cornou C, Cox BR, Di Giulio G, Forbriger T, Hayashi K, Lunedei E, Martin A, Mercerat D, Ohrnberger M, Poggi V, Renalier F, Sicilia D, Socco LV (2017) Guidelines for the good practice of surface wave analysis—a product of the InterPACIFIC project. Bull Earthq Eng 0:0–0 (submitted to)

Garofalo F, Foti S, Hollender F, Bard PY, Cornou C, Cox BR, Dechamp A, Ohrnberger M, Perron V, Sicilia D, Teague D, Vergniault C (2016a) InterPACIFIC project: comparison of invasive and non-invasive methods for seismic site characterization. Part II: inter-comparison between surface-wave and borehole methods. Soil Dyn Earthq Eng 82:241–254. doi:10.1016/j.soildyn.2015.12.009

Garofalo F, Foti S, Hollender F, Bard PY, Cornou C, Cox BR, Ohrnberger M, Sicilia D, Asten M, Di Giulio G, Forbriger T, Guillier B, Hayashi K, Martin A, Matsushima S, Mercerat D, Poggi V, Yamanaka H (2016b) InterPACIFIC project: comparison of invasive and non-invasive methods for seismic site characterization. Part I: intra-comparison of surface wave methods. Soil Dyn Earthq Eng 82:222–240. doi:10.1016/j.soildyn.2015.12.010
Griffiths SC, Cox BR, Rathje EM, Teague DP (2016) Mapping dispersion misfit and uncertainty in Vs profiles to variability in site response estimates. J Geotech Geoenviron Eng. doi:10.1061/(ASCE)GT. 1943-5606.0001553

Guéguen P, Péquegnat C, Souriau A, Bard P-Y, Cara M, Dominique P, Régnier M (2007) The French Accelerometric Network (RAP): current state in 2007. In: Proceedings of the international symposium on seismic risk reduction, Bucarest, Romania

Henstridge J (1979) A signal processing method for circular arrays. Geophysics 44:179–184. doi:10.1190/1.1440959

Horike M (1985) Inversion of phase velocity of long-period microtremors to the S-wave-velocity structure down to the basement in urbanized areas. J Phys Earth 33:59–96. doi:10.4294/jpe1952.33.59

Jakka RS, Roy N, Wason HR (2014) Implications of surface wave data measurement uncertainty on seismic ground response analysis. Soil Dyn Earthq Eng 61–62:239–245. doi:10.1016/j.soildyn.2014.02.004

Köhler A, Ohrnberger M, Scherbaum F, Wathelet M, Cornou C (2007) Assessing the reliability of the modified three-component spatial autocorrelation technique. Geophys J Int 168:779–796. doi:10.1111/j.1365-246X.2006.03253.x

Ktenidou O-J, Cotton F, Abrahamson NA, Anderson JG (2014) Taxonomy of kappa: a review of definitions and estimation approaches targeted to applications. Seismol Res Lett 85:135–146. doi:10.1785/0220130027

Lemoine A, Douglas J, Cotton F (2012) Testing the applicability of correlations between topographic slope and Vs30 for Europe. Bull Seismol Soc Am 102:2585–2599. doi:10.1785/0120110240

Lomax A, Snieder R (1995) Identifying sets of acceptable solutions to non-linear, geophysical inverse problems which have complicated misfit functions. Nonlinear Process Geophys 2:222–227

Martin A, Diehl J (2004) Practical experience using a simplified procedure to measure average shear-wave velocity to a depth of 30 meters (VS30). In: Proceedings of the 13th world conference in earthquake engineering, Vancouver Canada, p 952

Maufroy E, Cruz-Atienza VM, Cotton F, Gaffet S (2015) Frequency-scaled curvature as a proxy for topographic site-effect amplification and ground-motion variability. Bull Seismol Soc Am 105:354–367. doi:10.1785/0120140089

McMechan GA, Yedlin MJ (1981) Analysis of dispersive waves by wave field transformation. Geophysics 46:869–874. doi:10.1190/1.1441225

Michel C, Edwards B, Poggi V, Burjanek J, Roten D, Cauzzi C, Fah D (2014) Assessment of site effects in alpine regions through systematic site characterization of seismic stations. Bull Seismol Soc Am 104:2809–2826. doi:10.1785/0120140097

Nakamura Y (1989) A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Railway Technical Research Institute, Quarterly Reports

Neidell NS, Taner MT (1971) Semblance and other coherency measures for multichannel data. Geophysics 36:482–497. doi:10.1190/1.1440186

Nicoud G, Royer G, Corbin J-C, Lemeille F, Paillet A (2002) Creusement et remplissage de la vallée de l'Isère au Quaternaire récent. Géologie Fr 4

O'Neil A (2003) Full-waveform reflectivity for modelling, inversion and appraisal of seismic surface wave dispersion in shallow site investigations. PhD, The University of Western Australia and School of Earth and Geographical Sciences

Ohrnberger M, Vollmer D, Scherbaum F (2006) WARAN—a mobile wireless array analysis system for in field ambient vibration dispersion curve estimation. In: Proceedings of the first european conference on earthquake engineering and seismology, Geneva, Switzerland

Pequegnat C, Gueguen P, Hatzfeld D, Langlais M (2008) The French accelerometric network (RAP) and national data centre (RAP-NDC). Seismol Res Lett 79:79–89. doi:10.1785/gssrl.79.1.79

Pileggi D, Rossi D, Lunedei E, Albarello D (2011) Seismic characterization of rigid sites in the ITACA database by ambient vibration monitoring and geological surveys. Bull Earthq Eng 9:1839–1854. doi:10.1007/s10518-011-9292-0

Régnier J, Laurendeau A, Duval AM, Gueguen P (2010) From heterogeneous set of soil data to Vs profile: application on the French permanent accelerometric network (RAP) sites. In: Proceedings of the fourteenth ECEE—European conference of earthquake engineering, Ohrid, Republic of Macedonia

Safari J, O'Neill A, Matsuoka T, Sanada Y (2005) Applications of love wave dispersion for improved shear-wave velocity imaging. J Environ Eng Geophys 10:135–150. doi:10.2113/JEEG10.2.135

Sandikkaya MA, Akkar S, Bard P-Y (2013) A nonlinear site-amplification model for the next pan-European ground-motion prediction equations. Bull Seismol Soc Am 103:19–32. doi:10.1785/0120120008

Satoh T, Kawase H, Matsushima SI (2001) Differences between site characteristics obtained from microtremors, S-waves, P-waves, and codas. Bull Seismol Soc Am 91:313–334
SESAME team (2004) Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations: measurements, processing and interpretation. SESAME European research project

Seyhan E, Stewart JP, Ancheta TD, Darragh RB, Graves RW (2014) NGA-West2 site database. Earthq Spectra 30:1007–1024

Socco LV, Foti S, Comina C, Boiero D (2012) Comment on “Shear wave profiles from surface wave inversion: the impact of uncertainty on seismic site response analysis”. J Geophys Eng 9:241–243. doi:10.1088/1742-2132/9/2/241

Souriau A, Chaljub E, Cornou C, Margerin L, Calvet M, Maury J, Wathelet M, Grimaud F, Ponsolles C, Pequegnat C, Langlais M, Gueguen P (2011) Multimethod characterization of the French-Pyrenean Valley of Bagneres-de-Bigorre for seismic-hazard evaluation: observations and models. Bull Seismol Soc Am 101:1912–1937. doi:10.1785/0120100293

Thompson EM, Wald DJ, Worden CB (2014) A Vs30 map for California with geologic and topographic constraints. Bull Seismol Soc Am 104:2313–2321. doi:10.1785/0120130312

Tokimatsu K (1997) Geotechnical site characterization using surface waves. In: Proceedings of the first international conference on earthquake geotechnical engineering, Balkema, Ishihara, pp 1333–1368

Wathelet M (2005) Direct inversion of spatial autocorrelation curves with the neighborhood algorithm. Bull Seismol Soc Am 95:1787–1800. doi:10.1785/0120040220

Wathelet M (2008) An improved neighborhood algorithm: parameter conditions and dynamic scaling. Geophys Res Lett. doi:10.1029/2008GL033256

Wathelet M, Jongmans D, Ohrnberger M (2004) Surface-wave inversion using a direct search algorithm and its application to ambient vibration measurements. Surf Geophys 2:211–221

Wathelet M, Jongmans D, Ohrnberger M, Bonnefoy-Claudet S (2008) Array performances for ambient vibrations on a shallow structure and consequences over Vs inversion. J Seismol 12:1–19. doi:10.1007/s10950-007-9067-x

Xia J, Xu Y, Luo Y, Miller RD, Cakir R, Zeng C (2012) Advantages of using multichannel analysis of Love waves (MALW) to estimate near-surface shear-wave velocity. Surv Geophys 33:841–860. doi:10.1007/s10712-012-9174-2

Yin X, Xia J, Shen C, Xu H (2014) Comparative analysis on penetrating depth of high-frequency Rayleigh and Love waves. J Appl Geophys 111:86–94. doi:10.1016/j.jappgeo.2014.09.022

Yong A (2016) Comparison of measured and proxy-based Vs30 values in California. Earthq Spectra 32:171–192. doi:10.1193/013114EQS025M

Yong A, Martin A, Stokoe K, Diehl J (2013) ARRA-funded Vs30 measurements using multi-technique approach at strong-motion stations in California and central-eastern United States, U.S. Geological Survey Open-File Report. USGS