# Site effects observations and mapping on the weathered volcanic formations of Mayotte Island

Agathe Roulle, Marie Baillet, Didier Bertil, Cécile Cornou

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

## To cite this version:

Agathe Roulle, Marie Baillet, Didier Bertil, Cécile Cornou. Site effects observations and mapping on the weathered volcanic formations of Mayotte Island. Comptes Rendus. Géoscience, 2022, 354 (S2), pp.317-341. 10.5802/crgeos.151. hal-04041185

HAL Id: hal-04041185

https://brgm.hal.science/hal-04041185

Submitted on 31 Mar 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Comptes Rendus
## Géoscience
### Sciences de la Planète

Agathe Roullé, Marie Baillet, Didier Bertil and Cécile Cornou

**Site effects observations and mapping on the weathered volcanic formations of Mayotte Island**

Volume 354, Special Issue S2 (2022), p. 317-341

Published online: 9 December 2022

Issue date: 17 January 2023

https://doi.org/10.5802/crgeos.151

**Part of Special Issue:** The Mayotte seismo-volcanic crisis of 2018-2021 in the Comoros archipelago (Mozambique channel)

**Guest editors:** Jérôme Van der Woerd (Institut Terre Environnement de Strasbourg, UMR 7063 CNRS / Université de Strasbourg, 67084 Strasbourg, France), Vincent Famin (Laboratoire Géosciences Réunion, Université de La Réunion - IPGP, 97744 Saint-Denis, France) and Eric Humler (Professeur Université de Nantes, Laboratoire de Planétologie et Géosciences, UMR 6112, Faculté des Sciences et Techniques, Nantes Université, 44322 Nante, France)

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

This article is licensed under the
CREATIVE COMMONS ATTRIBUTION 4.0 INTERNATIONAL LICENSE.
http://creativecommons.org/licenses/by/4.0/

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

Les Comptes Rendus. Géoscience — Sciences de la Planète sont membres du
Centre Mersenne pour l'édition scientifique ouverte

www.centre-mersenne.org

e-ISSN : 1778-7025
régionale pour permettre le calcul en temps réel de cartes de mouvements forts et de dommages à des fins de gestion de crise. Enfin, pour les 8 sites de stations sismologiques étudiés dans cet article, nous proposons une première estimation des paramètres de sites clefs, à savoir : la géologie superficielle, la fréquence de résonance, un facteur d'amplification du PGA, une valeur du paramètre $V_{S,30}$, si disponible, et une estimation de la classe de sol EC8 au site.

**Keywords.** Site effects, Weathering, Mayotte, Seismic hazard, Seismic risk, Volcanic formations.

**Mots-clés.** Effets de site, Altération, Mayotte, Risque sismique, Aléa sismique, Formations volcaniques.

Published online: 9 December 2022, Issue date: 17 January 2023

# 1. Introduction

After the seismic sequence, initiated in May 2018, and characterized by hundreds of felt earthquakes and consecutive building damages [Sira et al., 2018], the evaluation of the seismic hazard in the whole of Mayotte territory became an imperative need for both crisis management and risk mitigation.

Mayotte is a cluster of volcanic islands on the southern part of the Somalia basin between Africa and Madagascar. It is part of the Comoros archipelago and comprises two main islands: Grande Terre on the west and Petite Terre on the east. Since 2011, the Mayotte territory must comply with the French Building Code regulations, which classifies it as a moderate seismicity level area (decrees 2010-1254 and 2010-1255). In May 2018, an exceptional telluric activity started ~50 km east of Mayotte, and remains active after more than two years, with 29 events of magnitude greater than 5, the highest magnitude reaching Mw 5.9 [Cesca et al., 2020, Lemoine et al., 2020, Bertil et al., 2021, Feuillet et al., 2021]. Until June 2018, only three accelerometric stations were operating on the island and could record the earthquakes: the two RESIF-RAP stations YTMZ and MDZA [RESIF, 1995] and the BRGM station MILA. During the main shock of May 15, 2018, the peak ground acceleration (PGA) reached 180 mg at the MILA station, located at Iloni on the eastern coast of Mayotte, near the top of a hill with outcropping clayey weathered formations, whereas the PGA recorded at the YTMZ station, situated on a rock site, at the same epicentral distance, was about 50 mg. These observations remind us that local geology plays an important role in modifying the characteristics of strong ground motion, both in terms of amplification and duration of seismic waves. The particularity of the superficial geological formations of Mayotte lies in their ancient volcanic origin [Nehlig et al., 2013] and the tropical conditions, which induces strong weathering of outcropping formations. The

presence of those particular formations is a complex subject, which is still difficult to take into account, especially in a tropical volcanic context. While many papers show examples of weathering profile imaging [for example Anbazhagan and Sitharam, 2009 in India; Wang et al., 2019, Keifer et al., 2019 in Wyoming, United States] with some of them being applied in geological contexts similar to Mayotte [Von Voigtlander et al., 2018, Nelson and McBride, 2019, in Hawaii], a few of them intended to address site effects. Nevertheless, there are some precedents such as Davis [1995] who modelled the effect of the presence of a weathered layer on the seismic response of a site or Spudich et al. [1996] and Graizer [2009] who studied the directional amplification effects related to topography and the amplification effects related to the soft surface layer (weathered shale) at the Tarzana site in California. Narayan and Kumar [2015] and Kumar et al. [2017] also worked on the effects of weathering on seismic motion on a hillside via 2D modelling. In a similar manner, Wang et al. [2018] performed 3D modelling of lithological and topographic effects under weathering conditions in Hong Kong. Some authors also treated site effects in weathered areas for the estimation of induced ground motion, as is the case of Havenith et al. [2002] on the Ananevo site in Kyrgyzstan on weathered granite or of Ma et al. [2019] in Tokushima, Japan. In both cases, the authors used H/V ratios to analyze the directivity of the seismic response according to geomorphological criteria. Di Naccio et al. [2017] analyzed the amplification effects related to the shape of the topography and the presence of fractures on the surface via, in particular, the polarization analysis of H/V signals. This work on the San Gregorio site in Italy, located on fractured rock and which suffered significant damage during the L'Aquila earthquake, shows the importance of weathering in the directionality of seismic motion. All these works show the importance of weathering in seismic strong motions, both in terms of amplification, modification of the spectral content and
directivity of the seismic energy. Since the 2000s, site effects have been studied on the island of Mayotte through specific site studies prior to the construction of public buildings [mainly scholarly buildings, e.g., Rey et al., 2012, Roullé et al., 2019]. Informative maps, including lithological and topographical site effect maps, were also produced in 2004 for the atlas of natural hazards on behalf of the local authorities [Audru et al., 2010]. Those maps were based on a “solid geology” low-resolution version of the geological map [Stieltjes, 1988] and are now regarded as obsolete in terms of description and mapping of superficial formations. Despite this limitation, they have shown high susceptibility to site effects of the volcanic formations of Mayotte. The seismic sequence initiated in 2018 led us (i) to update the site effects mapping to account for the amplification of strong motion related to the surficial geology in crisis management tools such as shake maps and damage scenarios and (ii) to characterize site effects at the seismological station sites in order to improve the metadata of the stations used for seismic monitoring of the territory.

## 2. Geological context

Mayotte Island is made up of a succession of volcanic deposits: the oldest goes back to 10 My in the south and the north-west, the most recent, found in Petite Terre, has been proposed to be of Holocene age [Zinke et al., 2013], but these deposits have not been dated. The succession of eruptive phases and quiescence, which allowed weathering and erosion, have shaped the present geology of Mayotte [for detailed information on the geological history of Mayotte, see Stieltjes, 1988, Debeuf, 2004, Audru et al., 2006, 2010, Nehlig et al., 2013]. The tropical climate of Mayotte is characterized by high rainfall levels, with an average annual rainfall exceeding 1500 mm over the whole island (http://www.meteofrance.yt/climat/description-du-climat, last access 14/12/2021). Coupled with the age of the rocks, this is largely responsible for the strong weathering of Mayotte's superficial formations.

In 2013, a 1:30,000 geological map of Mayotte [Lacquement et al., 2013, Nehlig et al., 2013], with focus on superficial geology and natural hazards, was produced with an attempt to better characterize weathering levels. A simplified version of this map is given in Figure 1. The main geological units that

outcrop on land are basaltic, phonolitic or trachytic lavas (later referred to as lava formations), volcaniclastic deposits and large thicknesses of weathered rocks (later referred to as superficial formations). It distinguishes, among the superficial formations, the anthropic deposits, the autochthonous formations resulting from supergene weathering processes (e.g., alloterites and isalterites) and the allochthonous formations which bring together all the sedimentary deposits from erosion and transport process (e.g., alluvium, beach sands and slope formations as colluvium or scree). A conceptual diagram to explain the typical pattern of the weathering profiles of Mayotte and the relation between the geological formations of Figure 1 is given in Figure 2.

## 3. Data description

### 3.1. *Earthquake records*

The earthquake data used in this paper come from the records of the YTMZ, MCHI, MILA, MTSB, KNKL, PMZI and R1EE2 stations described in Bertil et al. [2021] and Saurel et al. [2022] which provide continuous data in the framework of the Mayotte seismo-volcanic monitoring network (Réseau de surveillance volcanologique et sismologique de Mayotte—REVOSIMA). The geographic location of these stations is shown in Figure 3 and these are described in terms of name, network and type in Table 1. A particular case is the TBAD station, which corresponds to a test site for the future installation of a permanent seismological station on Petite Terre to replace the PMZI temporary station. The analyzed seismic events occurred from 08/03/2018 to 09/02/2020 and include the strongest events of the seismo-volcanic crisis initiated in May 2018 (Figure 4). Their locations are extracted from the work of Lemoine et al. [2020] until May 2019 and Bertil et al. [2021] from June 2019 up to February 2020. All their magnitudes come from Bertil et al. [2021] and correspond to local $M_{lv}$ magnitudes estimated from the regional localizations, generally stronger than the magnitudes calculated from the global networks. In the catalogue of Bertil et al. [2021], the comparison between the $M_w$ value given by G-CMT [Dziewonski et al., 1981, Ekström et al., 2012] and $M_{lv}$ (for $M_w$ between 4.9 and 5.9) and between the $m_b$ value given by the International Seismological Centre [2016] and $M_{lv}$ (between $m_b = 4.0$ and 5.9) shows that $M_{lv}$ overestimates
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Simplified geology of Mayotte derived from Lacquement et al. [2013].</figcaption></figure>

$M_w$ by 0.3 on average and $m_b$ by 0.4 on average. From this main dataset, regrouping more than 4900

events, we extracted a subset of 84 events in order to calculate H/V earthquake spectral ratios (i.e., the
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Conceptual diagram representing the typical weathering profile in Mayotte [modified from Nehlig et al., 2013]. Alloterites thicknesses are generally about few meters, whereas isalterites thicknesses can reach 20–70 m depending on the area.</figcaption></figure>

ratio between the Fourier amplitude spectra of the horizontal and the vertical components of seismic motion, Lermo and Chávez-García [1993]; see Supplementary Table S1). This subset of data includes (i) the 53 highest magnitude earthquakes ($M_{Iv} > 5$) recorded between 08/03/2018 and 27/06/2018 for stations YTMZ and MILA (the only stations operating during this period) and (ii) 31 additional earthquakes with magnitudes between 2.8 and 5.0 recorded between 25/05/2019 and 08/02/2020 for stations YTMZ, MCHI, MILA, MTSB, KNKL, PMZI and R1EE2.

## 3.2. *H/V noise spectral ratios and MASW profiles*

For 20 years, BRGM has been carrying out seismic hazard studies on sites in Mayotte Island. These studies generally include a lithological site effects analysis and have led us to acquire numerous geophysical measurements to characterize (i) the resonance frequency of the sites via the calculation of the H/V noise spectral ratio [i.e., the ratio between the Fourier spectra of the horizontal and the vertical components of the seismic noise recorded at the surface, Nakamura, 1989] and (ii) the S-wave velocity profile with depth inferred from the MASW active surface wave method (Multichannel Analysis of Surface Wave). The compilation of all these data allowed us to obtain 557 H/V measurements and 117 MASW profiles distributed over the whole island (Figure 3). An exhaustive list of the public reports from which the H/V and MASW data were collected is given in Supplementary Table S2.

# 4. Method of analysis

## 4.1. *Peak ground accelerations*

The PGA, defined as the highest amplitude of the absolute value of a signal in acceleration between the two horizontal components of the ground motion, is a common parameter to describe strong ground motion for seismic damage evaluation and mitigation purposes. PGAs were extracted directly from the raw signals through the scwfparam module of the Seiscomp monitoring tool [Cauzzi et al., 2013]. This module uses distance and magnitude criteria to select the signals to be processed and the size of the time window. The PGA is calculated on a filtered signal cut at 40 Hz for signals with a signal-to-noise ratio superior to a fixed threshold.

## 4.2. *H/V noise spectral ratios*

We calculated the H/V noise spectral ratio following the method proposed by Nakamura [1989] to characterize the site resonance frequency. The method assumes that the spectral ratio between the horizontal and vertical component of seismic noise is a good indicator for site effect evaluation in 1D subsurface conditions. In case of high impedance contrast between soft filling and stiff basement and under the 1D assumption, the H/V peak frequency can be associated to the fundamental resonance frequency of
<figure><img src="image_7.png" /><figcaption>Figure 3. Location of the H/V noise spectral ratios data (black dot), the MASW profiles (blue dot), and the seismological stations (red dot) used in this paper.</figcaption></figure>

the soft soil [e.g., Field and Jacob, 1995, Bonnefoy-Claudet et al., 2006] following the equation:

$$f_0 = V_s/4H, \qquad (1)$$

where *H* corresponds to the thickness of the soft layer and *V*<sub>S</sub> is the average shear-wave velocity in the soft layer.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Location of the seismic events used in this paper, from Lemoine et al. [2020] and Bertil et al. [2021]. The full set of earthquakes considered are represented in grey while the black dots indicate the subset of the data of 84 earthquakes used to calculate the H/V earthquake spectral ratios.</figcaption></figure>

<table><caption>**Table 1.** Description of seismic stations</caption><thead><tr><th>Stations</th><th>YTMZ</th><th>MILA</th><th>MTSB</th><th>MCHI</th><th>KNKL</th><th>TBAD</th><th>PMZI</th><th>R1EE2</th></tr></thead><tbody><tr><td>Network</td><td>RA</td><td>RA</td><td>1T</td><td>ED</td><td>QM</td><td>None</td><td>1T</td><td>AM</td></tr><tr><td>Site name</td><td>Mamoudzou</td><td>Iloni</td><td>M'Tsamboro</td><td>Chiconi</td><td>Kani-Keli</td><td>Badamiers</td><td>Pamandzi</td><td>Coconi</td></tr><tr><td>Latitude</td><td>-12.7577</td><td>-12.8481</td><td>-12.6804</td><td>-12.8329</td><td>-12.9571</td><td>-12.770231</td><td>-12.7993</td><td>-12.8354</td></tr><tr><td>Longitude</td><td>45.2307</td><td>45.1928</td><td>45.0847</td><td>45.1237</td><td>45.1042</td><td>45.277658</td><td>45.2743</td><td>45.1365</td></tr><tr><td>Sensor type</td><td>Acc</td><td>Acc</td><td>BB (HH)</td><td>BB (BH)</td><td>BB (HH)</td><td>BB (HH)</td><td>BB (HH)</td><td>RaspB</td></tr></tbody></table>

Station types are: Acc = accelerometer; BB (HH) = broadband 0–100 Hz; BB (BH) = broadband 0–50 Hz; RaspB = Raspberry shakes. Networks references are: AM [Raspberry Shake Community, 2016]; ED = http://www.edusismo.org/; RA [RESIF, 1995]; 1T [Feuillet et al., 2022] and QM: CNRS-INSU Tellus SISMAYOTTE project.

For characterization of seismological stations, we calculated the H/V noise spectral ratios by extracting 1 h of signal from the continuous data stream of the station during the night (quiet period). The calculations were made following the guidelines from the

SESAME project [SESAME, 2004] using the Geopsy software [Wathelet et al., 2020]. We used the following parameters for calculations of the H/V spectral ratios: data were filtered with a bandpass filter between 0.1 and 20 Hz, the selected windows are 50 s long, and
<table><caption>**Table 2.** Minimum, maximum and average S-wave window lengths in seconds used for the H/V earthquake spectral ratio calculations</caption><thead><tr><th></th><th>MCHI</th><th>MILA</th><th>MTSB</th><th>KNKL</th><th>YTMZ</th></tr></thead><tbody><tr><th>Minimum</th><td>62.3209</td><td>44.3768</td><td>25.6447</td><td>13.5924</td><td>26.6476</td></tr><tr><th>Maximum</th><td>147.9226</td><td>136.0315</td><td>114.9714</td><td>57.8797</td><td>87.2135</td></tr><tr><th>Average</th><td>101.0179</td><td>92.7994</td><td>50.2577</td><td>35.2106</td><td>56.2333</td></tr></tbody></table>

we applied a 40% Kono-Homachi smoothing on the Fourier spectra.

### 4.3. *H/V earthquake spectral ratios*

For seismological stations YTMZ, MCHI, MILA, MTSB and KNKL (see Figure 3 for stations location), the H/V noise spectral ratios were compared with the H/V earthquake spectral ratios from the subset of data comprising the 84 earthquakes listed in Supplementary Table S1. To calculate them, we followed the procedure described by Lermo and Chávez-García [1993], keeping only the S-wave part of our records. This work has been done manually. Selected window lengths vary between 13 s and 136 s, with an average length greater than 50 s for four of the five studied stations (Table 2). Discrepancies between window lengths are due to the available signal for each event at each station. The shortest windows occur when events happened in bursts. The filtering and smoothing parameters are the same as those used for the H/V noise spectral ratio calculation to obtain comparable spectral ratios.

For stations PMZI and R1EE2, we only calculated the H/V noise spectral ratio, mainly because of data availability problems due to the late installation of these stations. However, we chose to include our results to provide preliminary indications of site effects at these stations.

### 4.4. *Azimuthal variation of H/V noise spectral ratios*

We computed the azimuthal variation of H/V curves. The direction of the maximum H/V ratio in the horizontal plane is used as an estimate of the predominant polarization [Di Giulio et al., 2015]. This tool allows to study the directivity of the seismic response as a function of the morphology and geological structure of a site [Theodoulidis et al., 2018], for example as a function of slope [Del Gaudio et al., 2008], rock

fracturing [Di Naccio et al., 2017], or more broadly as a function of weathering and the presence of a soft layer on the surface [Havenith et al., 2002, Ma et al., 2019]. In this paper, we only show those analyses indicating a strong directivity effect, i.e., for which the amplitude of the H/V peaks varies with azimuth.

### 4.5. *Damping*

We finally used the damping tool available in Geopsy package [Wathelet et al., 2020] to check the natural origin of the resonance frequency [Dunand et al., 2002]. Following Dunand et al. [2002] and the SESAME guidelines [2004], we consider that the damping of a natural origin peak will present very low values (<1%). The damping test has been done in a systematic way but only the particular case of the temporary station TBAB has been presented in this article.

## 5. Site effects analysis at seismological station sites

### 5.1. *Peak ground acceleration comparisons*

In this paragraph, we compared the PGAs extracted from the earthquake signals listed in the Lemoine et al. [2020] catalog for the MILA, KNKL, MCHI, MTSB, MPZI and R1EE2 stations with those of the YTMZ station, which is considered as the reference station in Mayotte and which recorded the whole seismo-volcanic crisis initiated in 2018 (Figure 5). For each station and each PGA value, a PGA correction was applied to account for the difference in hypocentral distance between the considered station and the YTMZ station. For this, we applied the empirical relationship between PGA, local magnitude and hypocentral distance proposed by Bertil and Hoste-Colomer [2020] for Mayotte:

$$ \log(\text{PGA}) = 0.757 \times \text{M}_{\text{lv}} - 5.79 \times 10^{-5} \times \text{R}_{\text{hypo}} \\ - \log \text{R}_{\text{hypo}} - 4.433 \pm 0.29, \quad (2) $$
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Comparison of PGAs between the reference station YTMZ and the considered station. Black dots give the value of PGA for each recorded event, the grey line represents the reference line for which $PGA_{YTMZ} = PGA_{station}$, the dotted grey line represents the median amplification coefficient line. The median amplification coefficient corresponds to the ratio between $PGA_{station}$ and $PGA_{YTMZ}$ and is indicated in the legend.</figcaption></figure>
where $R_{hypo}$ represents the hypocentral distance in km, $M_{Iv}$ is the local magnitude on the vertical component, PGA unit is g.

We have a good sampling of PGA values between 0.1 and 10 mg with little data above 10 mg except for the MILA station, which has recorded the strongest ground motions of the crisis. In general, we observe that the PGA values at the stations studied are mostly higher than those observed at the YTMZ station, with median amplification values ranging from 1.4 for MCHI to 4.9 for MTSB. Except for station MCHI, these coefficients are superior to the soil coefficients given by the French building code (between 1.35 for soil class B and 1.8 for soil class E). We also observe a strong variability of PGA ratios between the considered station and YTMZ, with amplifications that can go from a factor of 1 to factors greater than 8 on some earthquakes (e.g., on MILA).

## 5.2. Earthquake and noise spectral ratios

We then calculated the H/V noise spectral ratios for each station of the study and compared them to the H/V earthquake spectral ratio if available (Figures 6 and 7). For all the stations, H/V noise spectral ratios and H/V earthquake spectral ratios show similar shapes, particularly in terms of resonance frequency. The amplitude of the H/V earthquake spectra ratio is higher than that of the noise one as observed before [see for example Haghshenas et al., 2008, for a detailed comparisons between noise and earthquake H/V spectral ratios], especially at frequencies lower than the resonance frequency.

The YTMZ station, located on recent volcaniclastic formations with little alteration, is considered since the beginning of the seismo-volcanic crisis as the reference station. It is a key station for monitoring the crisis since it is the only one that recorded the entire crisis. The H/V noise spectral ratio shows a value close to 1 for most of the frequency range and a slight amplification around 2.5 Hz. This amplification is not visible on the H/V earthquake spectral ratio, which shows a slight amplification at low frequencies. Although these spectral ratios are not perfectly flat and equal to 1 over the entire frequency range considered (the ideal case of a rocky site with no site effect), the observed amplifications have much lower amplitudes than those observed at the other sites. We can therefore consider that this station does not present

a significant site effect and is a good reference station for Mayotte.

The MILA station, located on a hill of about 50 m height made of isalterites, presents two close peaks at 1.6 and 2.2 Hz on the H/V ratios, both for the noise and earthquake spectral ratios. MASW measurements made at the site yield an S-wave velocity of approximately 290 m/s between 0 and 15 m depth as shown in Figure 8. The $V_{S,30}$ value could be estimated by extrapolating the $V_S(z)$ profile using Boore's [2004] formula and is about 350 m/s, which would correspond to a minimum thickness of soft formation at the surface of 40 m if (1) is considered, consistent with the geological knowledge of the site. For this station, the azimuthal analysis of the H/V noise spectral ratio (Figure 9) shows a variation in the amplitude of the H/V peaks but without the disappearance of any peak, with a maximum amplitude between N060 and N170, indicating a potential geometric effect. These observations suggest a complex combination of effects related to the nature of the involved rocks and the topography as observed by Spudich et al. [1996] on the Tarzana hill during the 1994 Northridge earthquake, or on the Rognes hill which suffered severe damages during the 1909 Provence earthquake [Glinsky et al., 2019], or in a fractured rock site subjected to significant damages during the L'Aquila earthquake in Italy [Di Naccio et al., 2017], or even in landslide areas as studied for example by Del Gaudio et al. [2014] in the Terano landslide site, Italy. In all those cases, both the geomechanical characteristics of the geological layers (fracturing, presence of a fault, presence of a low velocity layer in surface), combined with a geometric effect linked to the surrounding topography, produced important amplifications of ground motion associated with strong directivity patterns. The MTSB and KNKL stations, both located at the bottom of the slope on colluvium, show a single resonance centered around 5 Hz (Figure 6), corresponding to soft material of lesser thickness than on MILA, in coherence with the alteration model proposed in Figure 2 and the known geology. For these two stations, the study of the directivity of the H/V noise spectral ratio shows a strong variability of the peaks according to the azimuth (Figure 9). This is of particular interest for the MTSB station for which we observe a H/V resonance peak of lower amplitude than the other stations and spread between 4.5 and 6.5 Hz, and for which the directivity analysis of
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Comparison of H/V spectral ratios performed on earthquake and on seismic noise for stations YTMZ, MILA, MTSB, KNKL and MCHI.</figcaption></figure>

the H/V ratio indicates a predominance of the lowest resonance frequency (around 4.5 Hz) at N080 and a predominance of the highest resonance frequency (around 6.5 Hz) at N180. This could correspond to a topography-related surface wave separation [Roten et al., 2006] with the predominance of Love waves at low frequency and N080 and Rayleigh waves at higher

frequency and N180. For the KNKL station, the observed resonance seems simpler with a very strong peak at 5 Hz and the appearance of a secondary peak around 3 Hz for azimuths N040 to N080. Here again, geometric effects are possible.

The MCHI station shows a lower frequency peak around 2.9 Hz, with a very good correspondence
<figure><img src="image_11.png" /><figcaption>**Figure 7.** H/V noise spectral ratio for stations R1EE2 and PMZI.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 8.** MASW profiles for the seismological station MILA. The black profile was acquired on the station site whereas the grey profiles were acquired at 200 m of the station site. Each grey $V_S(z)$ profile corresponds to a different shot of the same MASW profile.</figcaption></figure>

between the H/V noise and earthquake spectral ratios both in terms of frequency and amplitude (Figure 6). This station is located on a high plateau area filled by a large thickness of isalterites (Combani pa-leosurface), consistent with the observed resonance frequency, which corresponds to a filling of several tens of meters thickness (of the order of 25 m thickness according to (1) if we consider a $V_S$ value of the

order of 305 m/s as indicated in Table 4, in the next section). For this station, the H/V directivity analysis (not shown here) indicates a stable peak in amplitude, independent of the azimuth considered, in coherence with the flat topography of the site.

The Station R1EE2, located on alloterites at the top of the alteration profile, has two distinct resonance frequencies (Figure 7): one, not very pronounced, around 1 Hz and the second, of higher amplitude, at 3.7 Hz. This weak low frequency amplification is also found on the MCHI station located about 1.5 km east of R1EE2 and on the H/V noise spectral ratios acquired north of the area, in Combani, in a similar geological context. The response of this station would therefore correspond to the presence of two distinct interfaces: one quite deep, around 75 m deep, and the other, more superficial, around 20 m deep, assuming a tabular environment with a surface layer of $V_S = 300$ m/s. These hypotheses will need to be validated by a more detailed study of the 3D structure underlying the station. The PMZI station, located on the volcanic formations of Petite Terre, presents a unique amplification at 2.6 Hz, probably corresponding to a rather thick soft layer (several tens of meters). The absence of MASW data in the vicinity of the area and the high spatial variability of the volcanic formations do not allow us to go further in the interpretation of the data. The last station analyzed here is the Badamiers site, TBAD. It is in fact a site located to the north of Petite Terre, on the volcanic formations, and potentially intended to host a permanent seismological station as part of the seismological monitoring of the Mayotte territory to
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Azimuthal variation of H/V spectral ratio for stations MILA, MTSB and KNKL.</figcaption></figure>

replace the PMZI temporary station currently used. Although this station was installed for a short time in 2020 and therefore does not have a sufficient earthquake record dataset for a complete analysis of site effects on its site, we have chosen to show the preliminary results obtained on this site to demonstrate the importance of analyzing site effects before any permanent installation. The H/V ratio of the TBAD

station (Figure 10) shows a rather peculiar pattern with peaks at 1.5 Hz, 8.3 Hz, and 12.4 Hz. The analysis of the damping on each of these frequencies confirms the natural origin of the peak at 1.5 Hz with an attenuation value of 4% and shows us that the peak at 8 Hz is of anthropic origin with a damping lower than 1% [Dunand et al., 2002, Guillier et al., 2007]. We observe the same result at 12 Hz. Looking at the H/V
<figure><img src="image_14.png" /><figcaption>**Figure 10.** H/V noise spectral ratio (top left), azimuthal variation of H/V spectral ratio (top right) and damping at 1.5 Hz (bottom left) and 8.3 Hz (bottom right) at the station site TBAD.</figcaption></figure>

ratio directivity, we can see a strong H/V peak directivity at 8 Hz, with a maximum amplitude between N020 and N080. This could correspond to the signature of the thermal power plant located about 500 m north of this site, in the N040 direction.

**6. Preliminary geology-based site effects mapping at regional scale**

The strong ground motion observations made during the seismo-volcanic crisis confirmed that the superficial geology of the island could lead to significant site effects, including on weathered volcanic formations [Bertil et al., 2019, 2021, Roullé et al., 2022]. Since the beginning of the 2000s, numerous local seismic hazard studies have shown that the Mayotte territory presents superficial geological formations likely to undergo lithological site effects [e.g., Audru et al., 2002, 2010]. This led us to produce, in the first months of the crisis, a new lithological site effects map at regional scale dedicated to the calculation of strong motion maps and damage scenarios for crisis management purposes [Taillefer et al., 2019, 2022].

This map has been upgraded since then including the analysis done in the framework of the revision of the local seismic hazard mapping of Mayotte, a study in progress on behalf of the local authorities. This paper describes the upgraded version of the preliminary map described in Taillefer et al. [2019, 2022].

For practical reasons of constrained time delays and costs, we decided to use only the available data (i.e., without carrying out additional measurement) and to base our site effects classification on the soil classes A–E derived from the European building code EuroCode 8 (EC8) [NF EN 1998-5, 2005] and described in Table 3. Our site effects map is based on the 1:30,000 geological map of Mayotte Island [Nehlig et al., 2013, Lacquement et al., 2013, see Figure 1]. To complete the geological data, we compiled the geophysical data of the local seismic hazard studies realized by BRGM since 2000 (Figure 3). The H/V measurements and MASW profiles have been analyzed respectively in terms of resonance frequency $f_0$, corresponding to the frequency of the maximum peak of the H/V spectral ratio, and $V_{S,30}$ value, a site effect proxy from EC8 soil classification. This is calculated
<figure><img src="image_15.png" /><figcaption>**Figure 11.** Distribution of the H/V peak values (dominant resonance frequency) regarding simplified surface geology given in Figure 1. For each box, the central line represents the median, the edges of the box represent the 25th and 75th percentiles, the whiskers (dotted lines) extend to 1.5 times the interquartile ranges, and the outliers (data with values beyond the ends of the whiskers) are plotted individually (crosses). Corresponding values are indicated in Table 4.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 12.** Distribution of $V_{S,30}$ values regarding simplified surface geology. The box description is the same as for Figure 11. Corresponding Q25, Q50 and Q75 values are indicated in Table 4.</figcaption></figure>

following the equation [Borcherdt, 1994]:

$$V_{S,30} = \frac{30}{\sum_{i=1,N} \frac{h_i}{V_i}}, \quad (3)$$

where $h_i$ and $V_i$ represent respectively the thickness and shear wave velocity of the $i$th layer and $N$ corresponds to the number of layers identified in the upper 30 m of the ground. The $f_0$ and $V_{S,30}$ values were directly extracted from existing reports with no extra calculation. For the geological formations for which geophysical data were available, we calculated the distribution of the $f_0$ and $V_{S,30}$ parameters (Figures 11 and 12).

The first observation that can be made is that all the considered geological formations present a significant site effect with resonance frequencies globally between 1 and 5 Hz (Figure 11 and Table 4). If we consider a single layer medium characterized by

shear-wave velocities between 200 and 300 m/s according to the available MASW profiles [see for example Rey et al., 2012, Roullé et al., 2019], this would correspond to thicknesses of superficial formations up to several tens of meters. These values are consistent with the geological outcrop data acquired for the achievement of the 2013 geological map where weathering profile thicknesses from 20 m to 70 m were observed [see Nehlig et al., 2013]. These observations corroborate the importance of site effects on Mayotte observed in Figure 5 both in terms of amplification and spatial extent. The spatial variability of the observed resonance frequencies can occur over distances of the order of a few hundred meters as shown by the example of the Dembeni site (Figure 13). On this example, we observe a strong variability of $f_0$ values on the isalterite formations located on the south-western zone of the map with values ranging from 3.7 to 13 Hz in less than 200 m of distance. On the contrary, the $f_0$ values observed on the south-eastern hill with outcropping isalterites (values ranging from 2.2 to 3.5 over a distance of 200 m) are rather stable. A detailed knowledge of the geometry and geomechanical characteristics of those weathered formations would be necessary to better understand such variability patterns. The analysis of
<table><caption>**Table 3.** Description of the EC8 soil classes used in this paper</caption><thead><tr><th>Soil class</th><th>Description of soil profile</th><th>V<sub>S,30</sub> parameter (m/s)</th></tr></thead><tbody><tr><td>A</td><td>Rock or other rock-like geological formation, including at most 5 m of weaker material at surface</td><td>&gt;800</td></tr><tr><td>B</td><td>Deposits of very dense gravel, or very stiff clay, at least several tens of m in thickness, characterized by a gradual increase of mechanical properties with depth</td><td>360–800</td></tr><tr><td>C</td><td>Deep deposits of dense or medium-dense sand, gravel or stiff clay with thickness from several tens to many hundreds of m</td><td>180–360</td></tr><tr><td>D</td><td>Deposits of loose-to-medium cohesionless soil (with or without some soft cohesive layers), or of predominantly soft-to-firm cohesive soil</td><td>&lt;180</td></tr><tr><td>E</td><td>A soil profile consisting of a surface alluvium layer with V<sub>S,30</sub> values of type C or D and thickness varying between about 5 m to 20 m, underlain by stiffer material with V<sub>S,30</sub> &gt; 800 m/s</td><td></td></tr></tbody></table>

<table><caption>**Table 4.** Statistical values for $f_0$ and $V_{S,30}$ parameters distribution (data count, Q25, Q50 and Q75) and EC8 soil class according to the simplified geology</caption><thead><tr><th>Simplified geology</th><th>f<sub>0</sub> (count)</th><th>f<sub>0</sub> (Q25)</th><th>f<sub>0</sub> (Q50)</th><th>f<sub>0</sub> (Q75)</th><th>V<sub>S,30</sub> (count)</th><th>V<sub>S,30</sub> (Q25)</th><th>V<sub>S,30</sub> (Q50)</th><th>V<sub>S,30</sub> (Q75)</th><th>EC8 class</th></tr></thead><tbody><tr><td>Anthropic fills</td><td>23</td><td>1.5</td><td>2.0</td><td>2.4</td><td>1</td><td>—</td><td>—</td><td>—</td><td>As surrounding formation</td></tr><tr><td>Alloterites</td><td>25</td><td>2.4</td><td>2.9</td><td>3.1</td><td>3</td><td>267</td><td>287</td><td>297</td><td>C</td></tr><tr><td>Isalterites</td><td>111</td><td>2.5</td><td>3.2</td><td>4.3</td><td>23</td><td>285</td><td>305</td><td>322</td><td>C</td></tr><tr><td>Slope formations (colluvium, screes, breccias)</td><td>153</td><td>3</td><td>4.2</td><td>5.3</td><td>42</td><td>296</td><td>336</td><td>380</td><td>B or C or E</td></tr><tr><td>Alluvium, beach sands</td><td>176</td><td>1.8</td><td>2.6</td><td>4.2</td><td>46</td><td>235</td><td>261</td><td>304</td><td>C</td></tr><tr><td>Volcaniclastic formations</td><td>66</td><td>1.9</td><td>2.5</td><td>3.6</td><td>0</td><td>—</td><td>—</td><td>—</td><td>A on Grande-Terre island, B or C on Petite Terre island</td></tr><tr><td>Lava formations</td><td>3</td><td>—</td><td>—</td><td>—</td><td>1</td><td>—</td><td>—</td><td>—</td><td>A</td></tr></tbody></table>

the $V_{S,30}$ parameter distribution (Figure 12) shows a quite similar trend for the four analyzed geological formations with most of the $V_{S,30}$ distribution (i.e., the interval between Q25 and Q75) ranging from 230 to 380 m/s. The autochthonous altered volcanic formations (alloterites and isalterites), which form the upper weathering profile on top of the fractured lava formations (see Figure 2 for a conceptual presentation of a characteristic weathering profile in Mayotte Island), present values around 270–320 m/s. Both their lithological characteristics, geometry, and $V_{S,30}$ values lead us to consider them as a C soil class. The

allochthonous formations classified as slope formations (mainly colluvium) present higher $V_{S,30}$ values between 300 and 380 m/s. These formations drape the sides of the reliefs: their thickness is prone to increase downstream and can reach 10 m at particular areas of accumulation. They present different facies from fine colluvium to boulder colluvium with strong lateral variations of facies and thickness. At the bottom of the reliefs, these formations tend to rest directly on bedrock with a clear contact between the two formations. Those formations are thus difficult to classify following EC8 criteria since they can
<figure><img src="image_17.png" /></figure>

**Figure 13.** Dominant resonance frequencies $f_0$ observed at the Dembeni site on a simplified geological map and topographic background. The location of the site is indicated by a rectangle in the map at the top right. The colours correspond to the geological formations shown in Figure 1. The shading used as a base map corresponds to the IGN DEM at a 5 m resolution.

be considered either as B, C or even E soil class in
certain configurations. At present, neither geologi-
cal nor geomorphological information permits us to
distinguish the different possible configurations: so
we decided to keep a B, or C or E classification in
our preliminary map but this point should be clari-
fied in future works. The alluvium and beach sands,
on the other hand, present $V_{S,30}$ values mainly be-
tween 235 and 304 m/s (Q25-Q75 interval) with a
median value of 261 m/s, which corresponds clearly
to a C soil class. The visible spread towards higher
values of $V_{S,30}$ (above 300 m/s) corresponds to ar-
eas at the edge of the basins where the alluvial de-
posit is thinner. For geological formations without
sufficient $V_{S,30}$ measurements to lead to a simple
EC8 classification (anthropic fills and lava or volcani-
clastic formations), we proceeded by expert analy-
sis. Thus, as the anthropic fills present variable me-

chanical characteristics and relatively small thick-
nesses and a very small area (less than 1% of the en-
tire area of the island), we preferred to assign them
an EC8 soil class similar to the surrounding geo-
logical formation. As for the lava formations, cor-
responding to the least altered volcanic deposits of
the island and essentially present on the coasts, the
tops of crests and the bottoms of gullies, we con-
sidered them as bedrock, that is to say, as an A
soil class.

Finally, the case of volcaniclastic formations
is more complex. They consist mainly in the
phreatomagmatic projections of Petite Terre and
Mamoudzou, which are characterized by ash
deposits that are often poorly consolidated but
very heterogeneous both in terms of granulometry
and lithology. The $f_0$ resonance frequency values
issued from the H/V measurements available to date
indicate the presence of lithological site effects on
these formations (Figure 11 and Table 4), but the spa-
tial distribution of the H/V measurements, mostly
around the airport area on Petite Terre, does not
allow us to confirm that observation on a more ex-
panded area. In addition, the reference station YTMZ
used to monitor the seismo-volcanic crisis and sit-
uated in Grande Terre, is also located on these for-
mations and does not present a significant site effect
(Figure 6). This station being preponderant in the
calculation of the strong motion maps, we decided to
assign the A soil class to all the volcaniclastic forma-
tions located on Grande Terre Island and a B or C soil
class for the specific deposits of phreatomagmatic
projections of Petite Terre. This decision ensures
that strong motions are not overestimated during
the real-time shake map calculation process, based
on the acceleration values observed at terrestrial
seismological stations and which takes into account
station site conditions. Work in progress for the local
authorities of Mayotte, with a forthcoming H/V and
MASW survey on Petite Terre, should make it possi-
ble to modulate this classification in the future. The
EC8 soil classes assigned to each surficial geological
formation defined in Figure 1 allowed us to build a
preliminary lithological site effects map (Figure 14)
based on a simplified geology of the island and us-
able for real-time calculation of strong motion maps
(shakemaps) and damage maps (SEISAID releases) in
case of significant seismic event [Auclair et al., 2015,
Guérin-Marthe et al., 2021]. The results show that
only 6% of the Mayotte territory can be considered as
EC8 site class A (rock site). Consequently, the urban
areas are mainly located on soils prone to site effects,
which could imply severe consequences in terms of
seismic risk. This map was built from existing and
partial data, especially in terms of spatial distribu-
tion. It therefore needs to be improved to better take
into account the geometric complexity and spatial
variability of the geological formations, for example
for volcaniclastic formations. Work is currently un-
derway as part of a project to revise the local seismic
hazard on behalf of the local authorities. The idea is
to integrate information on the weathering thickness
and geometry of surface formations via additional
in-situ geophysical measurements and analysis of
airborne EM data. This study could lead us shortly
to revise the preliminary lithological site effects map
presented in this paper.

7. Discussion and conclusion

The seismo-volcanic crisis initiated in May 2018 in
Mayotte led us to assess the strong motions on the
territory for crisis management and seismic monitor-
ing purposes via the calculation of strong motion and
damage maps (for risk scenarios, shake maps, seis-
mological bulletins, etc ...). In this context, an esti-
mation of the site effects, potentially responsible for
a local amplification of the seismic motion, is essen-
tial.

We first reanalyzed all the H/V and MASW data
acquired in the last 20 years by comparing them to
the recent geological map of the surface formations
of Mayotte, which takes into account the levels of
weathering of the various outcropping geological for-
mations [Lacquement et al., 2013, Nehlig et al., 2013].
The results show that the alteration of volcanic rocks
in Mayotte induces generalized site effects since all
the geological formations for which we have data
show a clear resonance (Figure 4, Table 4) and low
values of V<sub>S,30</sub> (Figure 5, Table 4).

These results are based on weak ground motion
(seismic noise) and need to be validated by compar-
ison with strong ground motion observations. The
seismo-volcanic crisis initiated in 2018 has allowed
the constitution of an unparalleled set of strong
ground motion data for the study of the response
of highly weathered volcanic formations under seis-
mic load. The first analyses show the importance of
the amplifications observed on the analyzed seismo-
logical stations, with median amplification factors of
1.4 to 4.9 for the PGA compared to the reference sta-
tion YTMZ (Table 5). In addition, the analysis of the
H/V spectral ratios showed the possible combina-
tion of lithological effects linked to the presence of
soft layers on the surface and topographical effects
linked to the geometry of the site (relief in particu-
lar), as is the case, for example, at the MILA station,
which recorded the strongest motions (180 mg) dur-
ing the main earthquake of 15/05/2018 (Mw = 5.9).
However, the work done here presents several limi-
tations since it only covers PGA values between 0.1
and 10 mg with little data above 10 mg (except for the
MILA station) and present a strong bias in azimuthal
covering since most of the analyzed event epicen-
ters are located full east. It would be interesting to
pursue the data collection and analysis to complete
those preliminary results. Further works on this ex-
<figure><img src="image_18.png" /><figcaption>**Figure 14.** Geology-based lithological site effects map in terms of EC8 soil classes. The corresponding surface areas for each soil class are indicated in the legend.</figcaption></figure>

ceptional dataset should advance our understanding of site effects in a highly weathered volcanic environment with direct operational application in Mayotte but also in the West Indies, the highest seismic hazard zone in France.

Beyond the analysis of this exceptional strong ground motion dataset, we also proposed a prelimi-

nary lithologic site effects map based on the 1:30,000 geological map produced in 2013 and EC8-derived soil classes estimated from the existing geophysical database (Figure 7). It shows that 94% of the Mayotte territory is likely to present a site effect linked to its superficial geology. This map is intended to be used in the framework of crisis management and seismic
<table><caption>**Table 5.** Summary of the information characterising the sites of the seismological stations studied in this paper</caption><thead><tr><th>Stations</th><th>YTMZ</th><th>MILA</th><th>MTSB</th><th>MCHI</th><th>KNKL</th><th>TBAD</th><th>PMZI</th><th>R1EE2</th></tr></thead><tbody><tr><td>Network</td><td>RA</td><td>RA</td><td>1T</td><td>ED</td><td>QM</td><td>None</td><td>1T</td><td>AM</td></tr><tr><td>Site name</td><td>Mamoudzou</td><td>Iloni</td><td>M'Tsamboro</td><td>Chiconi</td><td>Kani-Keli</td><td>Badamiers</td><td>Pamandzi</td><td>Coconi</td></tr><tr><td>Superficial geology</td><td>Volcaniclastic formation</td><td>Isalterites</td><td>Colluvium</td><td>Isalterites</td><td>Colluvium</td><td>Volcaniclastic formation</td><td>Volcaniclastic formation</td><td>Alloterites</td></tr><tr><td>Lithological effect</td><td>No</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td></tr><tr><td>Suspected topographical site effect</td><td></td><td>Yes</td><td>Yes</td><td></td><td>Yes</td><td></td><td></td><td></td></tr><tr><td>Resonance frequencies (Hz)</td><td></td><td>1.6 and 2.2</td><td>5</td><td>1 and 2.8</td><td>4.8</td><td>1.5</td><td>2.6</td><td>1 and 3.7</td></tr><tr><td>Median PGA ratio between the station and YTMZ</td><td>1</td><td>3.3</td><td>4.9</td><td>1.4</td><td>3.8</td><td>Not calculated</td><td>3.0</td><td>2.2</td></tr><tr><td>V<sub>S,30</sub> (m/s)</td><td></td><td>350</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Estimated EC8 class</td><td>A</td><td>C</td><td>B-C-E</td><td>C</td><td>B-C-E</td><td>B-C</td><td>B-C</td><td>C</td></tr></tbody></table>

For each station, we indicate: the station code, its corresponding network, the name of the site, the simplified surface geology (Figure 1) from the 1:30,000 geological map of Mayotte [Lacquement et al., 2013], the presence or absence of a lithological and topographical site effect on the station, the associated resonance frequency, the median PGA ratio between the site and the reference station YTMZ on the analyzed earthquakes, and the measured $V_{S,30}$ value, if available. In case of double H/V peaks, both frequencies are indicated and the dominant one (with the maximum amplitude) is underlined.

monitoring. Nevertheless, it presents several limitations:

* a poor spatial distribution of collected geophysical data, including a lack of data inland and in the south of the island, outside the main urban areas (Figure 2);

* a lack of knowledge of the 3D geometry of the interfaces of interest, which does not allow us, for example, to distinguish between soil classes B, C and E at the colluvium level (Figures 6 and 7);

* a lack of knowledge of the spatial variability of the geological formations concerned, which does not allow us, for example, to objectively characterize the volcaniclastic formations present especially on Petite Terre (Figure 7).

The mapping carried out will therefore evolve as knowledge of the Mayotte subsoil improves, particularly in terms of the thickness of surface layers and their geomechanical characteristics. Finally, the importance of the site effects observed on the strong recordings analyzed in this paper leads us to recommend

a systematic characterisation of the site parameters of the seismological stations in place within the REVOSIMA framework, following the recommendations of the European SERA project [Cultrera et al., 2021, Di Giulio et al., 2021] on the key parameters to be estimated as a priority and the protocols to be used to obtain good quality information. The most relevant indicators identified by Cultrera et al. [2021] are: the fundamental resonance frequency, the shear-wave velocity profile, the time-averaged shear-wave velocity over the first 30 m, the depth of both seismological and engineering bedrock, surface geology and soil class. A systematic calculation of the Frequency-Scaled Curvature proxy developed by Maufroy et al. [2015] should also provide information on whether a site is prone to such effect or not. In our case, the synthesis of the site characteristics from our work (Table 5) shows that a dedicated geophysical measurement campaign would provide the main missing proxy, namely the velocity profile $V_S(z)$ and the resulting parameter $V_{S,30}$, and thus complete the site metadata of each station to obtain a standardised characterisation according to European quality criteria.
## Conflicts of interest

Authors have no conflict of interest to declare.

## Acknowledgments

We thank all the technicians, engineers and researchers who participated in the installation and operation of the Mayotte seismological stations. RE-SIF data centre provides data of RA and 1T networks stations (YTMZ, MILA, MTSB and PMZI). We thank IPGP for having shared the data of QM network station KNKL, EDUSISMO (www.edusismo.org) for having installed and maintained the station MCHI, and IRIS data centre and Raspberry Shake SA data centre for having provided R1EE2 data.

We also thank J.-M. Saurel for providing the TBAD data.

This work was funded by the Ministry of Environment (Ministère de la Transition Ecologique—MTE), the Department of the Environment, Planning and Housing (DEAL) of Mayotte, and the Ministry of Research (Ministère de l'Enseignement Supérieur et de la Recherche et de l'Innovation) through the BRGM own funds.

Detailed comments provided by Julie Régnier, Vincent Famin and one anonymous reviewer greatly improved the manuscript.

## Supplementary data

Supporting information for this article is available on the journal's website under https://doi.org/10.5802/crgeos.151 or from the author.

## References

Anbazhagan, P. and Sitharam, T. G. (2009). Spatial variability of the depth of weathered and engineering bedrock using multichannel analysis of surface wave method. *Pure Appl. Geophys.*, 166, 409–428.

Auclair, S., Monfort, D., Colas, B., Langer, T., and Perrier, P. (2015). Evaluation rapide des bilans matériels et humains: une aide essentielle à la gestion opérationnelle des crises sismiques. In *9ème Colloque National AFPS*.

Auclair, S., Vanoudheusden, E., Rey, J., Roullé, A., Audru, J.-C., and Eucher, G. (2006). Projet de

piste longue de l'aéroport de Pamandzi (Mayotte) - Étude de l'aléa sismique et détermination des mouvements sismiques. Rapport BRGM/RP-55151-FR, 55 p.

Audru, J.-C., Bitri, A., Desprats, J.-F., Dominique, P., Eucher, G., Hachim, S., Jossot, O., Mathon, C., Nédellec, J.-L., Sabourault, P., Sedan, O., Stollsteiner, P., and Terrier-Sedan, M. (2010). Major natural hazards in a tropical volcanic island: a review for Mayotte Island, Comoros archipelago, Indian Ocean. *Eng. Geol.*, 114, 264–381.

Audru, J.-C., Bitri, A., Desprats, J.-F., Mathon, C., Maurillon, N., Sabourault, P., Terrier-Sedan, M., and Sedan, O. (2002). Projets risques naturels et érosion à Mayotte : résultats d'année 1. Rapport BRGM/RP-51738-FR, 49 p.

Audru, J. C., Guennoc, P., Thinon, I., and Abellard, O. (2006). BATHYMAY: underwater structure of Mayotte Island revealed by multibeam bathymetry. *C. R. Géosci.*, 338, 1240–1249.

Bertil, D. and Hoste-Colomer, R. (2020). Loi Magnitude-Intensité pour Mayotte. Rapport final. Rapport BRGM/RP-70297-FR, 54 p.

Bertil, D., Lemoine, A., and Roullé, A. (2019). Evaluation d'un niveau d'aléa sismique probabiliste à Mayotte et à la Réunion. Rapport final. Rapport BRGM/RP69481-FR, 93 p.

Bertil, D., Mercury, N., Doubre, C., Lemoine, A., and Van der Woerd, J. (2021). The unexpected Mayotte 2018-2020 seismic sequence: a reappraisal of the regional seismicity of the Comoros. *C. R. Géosci.*, 353, 211–235.

Bonnefoy-Claudet, S., Cotton, F., and Bard, P. Y. (2006). The nature of noise wavefield and its applications for site effects studies: a literature review. *Earth-Sci. Rev.*, 79(3–4), 205–227.

Boore, D. M. (2004). Estimating Vs30 (or NEHRP site classes) from shallow velocity models (depths <30 m). *Bull. Seismol. Soc. Am.*, 94(2), 591–597.

Borcherdt, R. D. (1994). Estimates of site-dependent response spectra for design (methodology and justification). *Earthq. Spectra*, 10(4), 617–653.

Cauzzi, C., Clinton, J., Becker, J., and Kästli, P. (2013). Scwfparam: a tool for rapid parameterisation of ground motions and input to ShakeMap in Seis-ComP3. In *SSA2013 Annual Meeting, 17–19 April, Salt Lake City—Abstract no. 13-347*.

Cesca, S., Letort, J., Razafindrakoto, H. N. T., Heimann, S., Rivalta, E., Isken, M. P., Passarelli, L.,
Petersen, G. M., Cotton, F., and Dahm, T. (2020). Drainage of a deep magma reservoir near Mayotte inferred from seismicity and deformation. *Nat. Geosci.*, 13, 87–93.

Cochery, R., Lauret, G., Rey, J., and Forestier, L. (2014). Caractérisation des aléas mouvements de terrain, inondation et sismiques sur le futur site d'implantation du futur lycée de la Kwalé. Rapport BRGM/RP-63836-FR, 78 p.

Cultrera, G., Cornou, C., Di Giulio, G., and Bard, P.-Y. (2021). Indicators for site characterization at seismic station: recommendation from a dedicated survey. *Bull. Earthq. Eng.*, 19, 4171–4195.

Davis, R. (1995). Effects of weathering on site response. *Earthq. Eng. Struct. Dyn.*, 24, 301–309.

Debeuf, D. (2004). *Étude de l'évolution volcano-structurale et magmatique de Mayotte, Archipel des Comores, Océan Indien: approches structurale, pétrographique, géochimique et géochronologique*. Doctoral dissertation, Université de la Réunion, https://tel.archives-ouvertes.fr/tel-00646477.

Del Gaudio, V., Coccia, S., Wasowski, J., Gallipoli, M. R., and Mucciarelli, M. (2008). Detection of directivity in seismic site response from microtremor spectral analysis. *Nat. Hazards Earth Syst. Sci.*, 8, 751–762.

Del Gaudio, V., Muscillo, S., and Wasowski, J. (2014). What we can learn about slope response to earthquakes from ambient noise analysis: an overview. *Eng. Geol.*, 182, 182–200.

Di Giulio, G., Amoroso, S., Di Naccio, D., Falcucci, E., Gori, S., Hailemikael, S., Vasallo, M., Ciampaglia, A., De Luca, G., Del Grosso, A., Di Marcantonio, P., D'Onofrio, K., Porreca, M., Trotta, D., Villani, F., and Milana, G. (2015). The seismic microzonation of San Gregorio through a multidisciplinary approach. Seismic amplification in a stiff site. *Eng. Geol. Soc. Territ.*, 5, 1137–1141.

Di Giulio, G., Cultrera, G., Cornou, C., Bard, P.-Y., and Al-TFaily, B. (2021). Quality assessment for site characterization at seismic stations. *Bull. Earthq. Eng.*, 19, 4643–4691.

Di Naccio, D., Vassallo, M., Di Giulio, G., Amoroso, S., Cantore, L., Hailemikael, S., Falcucci, E., Gori, S., and Milana, G. (2017). Seismic amplification in a fractured rock site. The case study of San Gregorio (L'Aquila, Italy). *Phys. Chem. Earth, Parts A/B/C*, 98, 90–106.

Dunand, F., Bard, P. Y., Chatelain, J. L., Guéguen,

P., Vassail, T., and Farsi, M. N. (2002). Damping and frequency from Randomdec method applied to in situ measurements of ambient vibrations. Evidence for effective soil structure interaction. In *12th European Conference on Earthquake Engineering, London*. Paper (Vol. 869).

Dziewonski, A. M., Chou, T. A., and Woodhouse, J. H. (1981). Determination of earthquake source parameters from waveform data for studies of global and regional seismicity. *J. Geophys. Res.*, 86, 2825–2852.

Ekström, G., Nettles, M., and Dziewonski, A. M. (2012). The global CMT project 2004–2010: centroid-moment tensors for 13,017 earthquakes. *Phys. Earth Planet. Inter.*, 200–201, 1–9.

Feuillet, N., Jorry, S., Crawford, W., Deplus, C., Thimon, I., Jacques, E., Saurel, J. M., Lemoine, A., Pauquet, F., Satriano, C., Aiken, C., Foix, O., Kowalski, P., Laurent, A., Rinnert, E., Cathalot, C., Donval, J. P., Guyader, V., Gaillot, A., Scalabrin, C., Moreira, M., Peltier, A., Beauducel, F., Grandin, R., Ballu, V., Daniel, R., Pelleau, P., Besancon, S., Geli, L., Bernard, P., Bachelery, P., Fouquet, Y., Bertil, D., Lemarchand, A., and Van der Woerd, J. (2021). Birth of a large volcano offshore Mayotte through lithosphere-scale rifting. *Nat. Geosci.*, 14, 787–795.

Feuillet, N., Van der Woerd, J., and RESIF (2022). Seismic sequence monitoring on land and at sea in Mayotte: SISMAYOTTE (RESIF-SISMOB) [Data set]. RESIF - Réseau Sismologique et géodésique Français.

Field, E. H. and Jacob, K. H. (1995). A comparison and test of various site-response estimation techniques, including that are not reference-site dependent. *Bull. Seismol. Soc. Am.*, 85, 1127–1143.

Glinsky, N., Bertrand, E., and Régnier, J. (2019). Numerical simulation of topographical and geological site effects. Applications to canonical topographies and Rognes hill, South East France. *Soil Dyn. Earthq. Eng.*, 116, 620–636.

Graizer, V. (2009). Low-velocity zone and topography as a source of site amplification effect on Tarzana hill, California. *Soil Dyn. Earthq. Eng.*, 29(2), 324–332.

Guérin-Marthe, S., Gehl, P., Fayjaloun, R., Caterina, N., and Auclair, S. (2021). Rapid earthquake response: the state-of-the art and recommendations with a focus on European systems. *Int. J. Disaster
*Risk Reduct.*, 52, article no. 101958.

Guillier, B., Chatelain, J.-L., Bonnefoy-Claudet, S., and Haghshenas, E. (2007). Use of ambient noise: from spectral amplitude variability to H/V stability. *J. Earthq. Eng.*, 11(6), 925–942.

Haghshenas, E., Bard, P. Y., and Theodulidis, N. (2008). Empirical evaluation of microtremor H/V spectral ratio. *Bull. Earthq. Eng.*, 6(1), 75–108.

Havenith, H.-B., Jongmans, D., Faccioli, E., Abdrakhmatov, K., and Bard, P.-Y. (2002). Site effect analysis around the seismically induced Ananevo Rockslide. *Kyrgyzstan Bull. Seismol. Soc. Am.*, 92, 3190–3209.

International Seismological Centre (2016). On-line Bulletin, http://www.isc.ac.uk, Internatl. Seismol. Cent., Thatcham, United Kingdom.

Keifer, I., Dueker, K., and Chen, P. (2019). Ambient Rayleigh wave field imaging of the critical zone in a weathered granite terrane. *Earth Planet. Sci. Lett.*, 510, 198–208.

Kumar, V., Narayan, J. P., and Kumar, M. (2017). Numerical study of effect of ridge-weathering thickness effect on the ground motion amplification. *Earth Sci. India*, 10, 167–178.

Lacquement, F., Nehlig, P., and Bernard, J. (2013). *Carte géol. France (1/30,000), feuille Mayotte (1179)*. BRGM, Orléans.

Lemoine, A., Briole, P., Bertil, D., Roullé, A., Foumelis, M., Thinon, I., Raucoules, D., de Michele, M., Valty, P., and Hoste Colomer, R. (2020). The 2018–2019 seismo-volcanic crisis east of Mayotte, Comoros islands: seismicity and ground deformation markers of an exceptional submarine eruption. *Geophys. J. Int.*, 223, 22–44.

Lermo, J. and Chávez-García, F. J. (1993). Site effect evaluation using spectral ratios with only one station. *Bull. Seismol. Soc. Am.*, 83(5), 1574–1594.

Ma, N., Wang, G., Kamai, T., Doi, I., and Chigira, M. (2019). Amplification of seismic response of a large deep-seated landslide in Tokushima, Japan. *Eng. Geol.*, 249, 218–234.

Maufroy, E., Cruz-Atienza, V. M., Cotton, F., and Gaffet, S. (2015). Frequency-scaled curvature as a proxy for topographic site-effect amplification and ground-motion variability. *Bull. Seismol. Soc. Am.*, 105(1), 354–367.

Nakamura, Y. (1989). A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. *Railw. Tech. Res.*

*Inst. Q. Rep.*, 30(1), 25–33.

Narayan, J. P. and Kumar, V. (2015). A numerical study of effects of ridge-weathering and ridge-shape-ratio on the ground motion characteristics. *J. Seismol.*, 19(1), 83–104.

Nehlig, P., Lacquement, F., Bernard, J., Caroff, M., Deparis, J., Jaouen, T., Pelleter, A.-A., Perrin, J., Prognon, C., and Vittecoq, B. (2013). *Notice explicative, Carte géol. France (1/30,000), feuille Mayotte (1179)*. BRGM, Orléans.

Nelson, S. and McBride, J. (2019). Application of HVSR to estimating thickness of laterite weathering profiles in basalt. *Earth Surf. Process. Landf.*, 44(7), 1365–1376.

NF EN 1998-5 (2005). *Eurocode 8 - Calcul des structures pour leur résistance aux séismes - Partie 5 : Fondations, ouvrages de soutènement et aspects géotechniques (P 06-035-1)*. AFNOR.

Oppermann, A., Auclair, S., Bitri, A., Nédellec, J. L., Rey, J., Roullé, A., and Vanoudheusden, E. (2009). Evaluation de l'aléa sismique sur les sites de l'extension du collège de Labattoir, de la construction du collège de Bouéni, et du lycée de Chirongui (Mayotte). Rapport BRGM/RP-57715-FR, 126 p.

Raspberry Shake Community (2016). OSOP, S.A.; Gempa GmbH.

RESIF (1995). RESIF-RAP French Accelerometric Network. RESIF – Réseau Sismologique et géodésique Français.

Rey, J., Auclair, S., Vanoudheusden, E., Oppermann, A., François, B., and Bitri, A. (2010). Evaluation de l'aléa sismique sur le site de l'extension de la maison d'arrêt de Majicavo-Lamir (Mayotte). Rapport BRGM/RP-58025-FR, 73 p.

Rey, J., Bastone, V., Bitri, A., Cochery, R., François, B., Roullé, A., and Tardy, D. (2012). Evaluation de l'aléa sismique pour 8 sites prévus pour la construction d'établissements scolaires à Mayotte. Rapport BRGM/RP-61170-FR, 176 p.

Rey, J., Bastone, V., Bitri, A., François, B., Lauret, G., and Tardy, D. (2013a). Evaluation de l'aléa sismique pour 3 sites d'extensions d'établissements scolaires à Mayotte. Rapport BRGM/RP-62972-FR, 67 p.

Rey, J., Bastone, V., Bitri, A., François, B., Lauret, G., and Tardy, D. (2013b). Evaluation de l'aléa sismique pour le site du projet d'extension de la Préfecture de Mayotte. Rapport BRGM/RP-62965-FR, 38 p.
Rey, J. and Cochery, R. (2012). Evaluation de l'aléa sismique selon la nouvelle réglementation pour la construction d'un collège à Bambo Ouest à Mayotte. Rapport BRGM/RP-61593-FR, 32 p.

Rey, J., Roullé, A., De Martin, F., Audru, J.-C., and Eucher, G. (2006a). Étude de l'aléa sismique et détermination des mouvements sismiques. Lycée de Pamandzi (Mayotte). Rapport BRGM/RP-54471-FR, 31 p., 7 Illust., 6 tabl., 1 annexe.

Rey, J., Roullé, A., De Martin, F., Audru, J.-C., Miehe, J.-M., Bitri, A., and Nédellec, J.-L. (2006b). Étude de l'aléa sismique et détermination des mouvements sismiques. LEP de Kaouéni (Mayotte). Rapport BRGM/RP-54470-FR, 49 p.

Rey, J., Roullé, A., De Martin, F., Audru, J.-C., Miehe, J.-M., Nédellec, J.-L., Eucher, G., and Bitri, A. (2006c). Étude de l'aléa sismique et détermination des mouvements sismiques. Futur lycée de Dembéni (Mayotte). Rapport BRGM/RP-54467-FR, 45 p.

Rey, J., Roullé, A., Vanoudheusden, E., and Audru, J. C. (2006d). Étude de l'aléa sismique et détermination des mouvements sismiques sur le site de l'IFM de Dembéni (Mayotte). Rapport BRGM/RP-54960-FR, 31 p.

Rey, J., Vanoudheusden, E., and Lauret, G. (2013c). Evaluation de l'aléa sismique pour le site du projet d'extension du Centre Hospitalier de Mamoudzou (CHM). Rapport BRGM/RP-63175-FR, 34 p.

Roten, D., Fäh, D., Cornou, C., and Giardini, D. (2006). Two-dimensional resonances in Alpine valleys identified from ambient vibration wavefields. *Geophys. J. Int.*, 165, 889-905.

Roullé, A., Audru, J.-C., Bernardie, S., and Nédellec, J.-L. (2005a). Étude de l'aléa sismique et détermination des mouvements sismiques sur le site du lycée de Chirongui (Mayotte). Rapport BRGM/RP-53832-FR, 61 p.

Roullé, A., Belvaux, M., Abad, J., François, B., Colombain, A., and Noury, G. (2019). Evaluation de l'aléa sismique pour 10 sites de construction/extension d'établissements scolaires à Mayotte. Rapport BRGM/RP-68294-FR, 236 p.

Roullé, A., Bernardie, S., Bitri, A., Miehe, J. M., Audru, J. C., and Eucher, G. (2005b). Étude de l'aléa sismique et détermination des mouvements sismiques sur le site de l'hôpital de Mamoudzou (Mayotte). Rapport BRGM/RC-53975-FR, 45 p., 11 fig., 11 tabl., 2 annexes.

Roullé, A., Raingeard, A., Cochery, R., Baptiste,

J., Vanoudheusden, E., Reninger, P. A., Bitri, A., Sadeski, L., Beaubois, F., and Dectot, G. (2022). Méthodologie de révision de l'aléa sismique local et de la liquéfaction : application aux zones pilotes de Dembeni, Longoni et Combani. Rapport final. BRGM/RP-71365-FR.

Sabourault, P., Sedan, O., Audru, J. C., Mathon, C., Nédellec, J. L., Bitri, A., Bour, M., and Lebrun, B. (2002). Microzonage sismique de Mamoudzou centre, Mayotte. Rapport BRGM/RP-51839-FF, 107 p.

Saurel, J.-M., Jacques, E., Aiken, C., Lemoine, A., Retailleau, L., Lavayssière, A., Foix, O., Dofal, A., Laurent, A., Mercury, N., Crawford, W., Lemarchand, A., Daniel, R., Pelleau, P., Bès de Berc, M., Dectot, G., Bertil, D., Roullé, A., Broucke, C., Colombain, A., Jund, H., Besançon, S., Guyavarch, P., Kowalski, P., Roudaut, M., Apprioual, R., Battaglia, J., Bodihar, S., Boissier, P., Bouin, M. P., Brunet, C., Jamale, K., Catherine, P., Desfete, N., Doubre, C., Dretzen, R., Dumouche, T., Fernagu, P., Ferrazzini, V., Fontaine, F. R., Gaillot, A., Géli, L., Griot, C., Grunberg, M., Can Guzel, E., Hoste-Colomer, R., Lambotte, S., Lauret, F., Léger, F., Maros, E., Peltier, A., Vergne, J., Satriano, C., Tronel, F., Van der Woerd, J., Fouquet, Y., Jorry, S. J., Rinnert, E., Thinon, I., and Feuillet, N. (2022). Mayotte seismic crisis: building knowledge in near real-time by combining land and ocean-bottom seismometers, first results. *Geophys. J. Int.*, 228(2), 1281-1293.

SESAME (2004). Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation. Available at http://sesame-fp5 OBS. ujf-grenoble.fr/Papers/HV_User_Guidelines.pdf.

Sira, C., Schlupp, A., Bontemps, M., Regis, E., and Van der Woerd, J. (2018). Essaim sismique à l'est de Mayotte. Analyse pour la période du 10 mai au 15 juin 2018. Note préliminaire du BCSF-RENASS, BCSF-RENASS2018-R4, 62 p.

Spudich, P., Hellweg, M., and Lee, W. H. K. (1996). Directional topographic site response at Tarzana observed in aftershocks of the 1994 Northridge, California, earthquake: implications for mainshock motions. *Bull. Seismol. Soc. Am.*, 86(1B), S193-S208.

Stieltjes, L. (1988). *Geological map of Mayotte at 1/50,000*. BRGM, Orléans.

Taillefer, N., Bertil, D., Negulescu, C., Roullé, A., Iasio,
C., Auclair, S., and Quique, R. (2022). Seismic damage scenarios for Mayotte: a tool for disaster management. *C. R. Géosci.*, 354(S2), 343-359.

Taillefer, N., Negulescu, C., Iasio, C., Bertil, D., and Roullé, A. (2019). Scénarios de dommages sismiques pour Mayotte. Rapport final. BRGM/RP-69518 -FR, 63 p.

Theodoulidis, N., Cultrera, G., Cornou, C., Bard, P. Y., Boxberger, T., DiGiulio, G., Imtiaz, A., Kementzetzidou, D., Makra, K., and the Argostoli NERA Team (2018). Basin effects on ground motion: the case of a high-resolution experiment in Cephalonia (Greece). *Bull. Earthq. Eng.*, 16(2), 529-560.

Vermeersch, F., Auclair, S., Vanoudheusden, E., Rey, J., Roullé, A., Bitri, A., and Cluzet, C. (2008). Étude de l'aléa sismique et détermination des mouvements sismiques. Futur lycée de Dembéni, site de Tsararano (Mayotte). Rapport BRGM/RP-56431-FR, 52 p.

Von Voigtlander, J., Clark, M. K., Zekkos, D., Greenwood, W. W., Anderson, S. P., Anderson, R. S., and Godt, J. W. (2018). Strong variation in weathering of layered rock maintains hillslope-scale strength un-

der high precipitation. *Earth Surf. Process. Landforms*, 43(6), 1183-1194.

Wang, G., Du, C., Huang, D., Jin, F., Koo, R. C., and Kwan, J. S. (2018). Parametric models for ground motion amplification considering 3D topography and subsurface soils. *Soil Dyn. Earthq. Eng.*, 115, 41-54.

Wang, W., Chen, P., Keifer, I., Dueker, K., Lee, E. J., Mu, D., Jiao, J., Zhang, Y., and Carr, B. (2019). Weathering front under a granite ridge revealed through full-3D seismic ambient-noise tomography. *Earth Planet. Sci. Lett.*, 509, 66-77.

Wathelet, M., Chatelain, J. L., Cornou, C., Giulio, G. D., Guillier, B., Ohrnberger, M., and Savvaidis, A. (2020). Geopsy: a user-friendly open-source tool set for ambient vibration processing. *Seismol. Res. Lett.*, 91(3), 1878-1889.

Zinke, J., Reijmer, J. J. G., Thomassin, B. A., Dullo, W.-C., Grootes, P. M., and Erlenkeuser, H. (2013). Postglacial flooding history of Mayotte Lagoon (Comoro Archipelago, southwest Indian Ocean). *Mar. Geol.*, 194, 181-196.