Twelve available two-way time high-resolution seismic reflection profiles located in the central part of the middle Rhône valley are interpreted. In addition, one of the profiles was reprocessed to determine the P-wave velocities of the main geological units and to convert this profile into a depth cross section. The Lower and Upper Cretaceous units are clearly identifiable on all the profiles, along with the Messinian Erosion Surface (MES) carved out during the Messinian Salinity Crisis (MSC) by the paleo-Rhône and its western tributaries, the Ardèche and Cèze paleo-canyons. The Plio-Quaternary fill of these paleo-canyons shows at least 4 main units with an overall transgression. The combination of geological data from geological maps, geological field surveys and borehole data made it possible to model the MES in 3D at the scale of the region, and to produce depth/elevation model. From a geological point of view, the interpretation of the seismic profiles enabled us to reconstruct the stages in the sub-aquatic filling of the Messinian-Pliocene aggradation of the paleo-river. Several Mass Transport Deposits (MTDs) were identified both during the drop and during the rise in the Mediterranean Sea level. From a geomorphological point of view, this study provides new insights in the route and longitudinal profile of paleo-rivers and, in particular, it deepens the profile of the Paleo-Rhône at the latitude of the Tricastin region (up to −700 m b.s.l.) and significantly modifies the course and depth of the Ardèche proposed in previous studies. The N-Ardèche river, known to develop a karstic system during the MSC, is connected to a deep canyon, most likely through a karstic pocket valley, as suggested by the very steep longitudinal profile of the MES. Finally, from a structural point of view, our interpretation of the seismic profiles shows a broad ENE-trending anticline structure associated with a normal fault which apparently did not affect the Mio-Pliocene fill. In the southern part of the area, near the Uchaux anticline, the imaged structures suggest the presence of a recent (syn- to post-Pliocene) fault propagation fold. In addition to all the new information on the geology, morphology and methods of excavation and filling of the Messinian paleo-canyon, the proposed topographic model of the paleo-canyon is crucial for modelling seismic movement in the context of a basin with a complex geometry and, in particular, for the numerical assessment of site effects in a context of low seismicity.

sensitivity and resolution capability have an impact on the numerical amplification predicted in the basin. In particular, this ANSWT model lacks clear basin edges in order to efficiently trap seismic waves in the basin and to generate significant 3D wave propagation effects (diffractions, reflections, and generation of laterally propagating surface waves at the edges of the basin). As a result, the numerical amplification predicted in the ANSWT model remains dominated by a 1D response and does not reproduce the broadband character of the observed amplification at locations affected by significant 3D propagation effects. On the other hand, the numerical amplification predicted in the ANSWT model shows a good agreement with the observations at locations that seem less affected by 3D propagation effects, including in complex regions of the model where lateral variations must be taken into account. Our results therefore contribute to identify and better understand the potential and limitations of using ANSWT models for numerical site effect estimation. This study allows us to propose perspectives for future work to improve the approach, which remains promising for site effect assessment in low- to moderate-seismicity contexts.

SUMMARY The local ground motion amplification related to the geology at a specific site (i.e. the so-called site effects) may be classically quantified through the SSR (standard spectral ratio) technique applied on earthquake recordings. However, such a quantification might be challenging to carry out in low-to-moderate seismicity regions. Methods based on background ambient noise, such as noise-based standard spectral ratio (SSRn), might be of great interest in these areas. But noise-derived amplification is particularly sensitive to local anthropogenic sources, which may introduce biases in the evaluation of site effects by dominating the geological effects, especially for frequencies higher than 1 Hz. A hybrid approach (SSRh), developed to reduce biases in noise-based spectral ratios by combining classical earthquake-based spectral ratio (SSR) and SSRn, was recently introduced and relies on a site reference. We here investigate the applicability of SSRn and SSRh in a heavily industrialized environment in the Tricastin Valley (south-east France), where critical facilities are located. We continuously recorded ambient noise from 2020 February to March on a 400-sensor seismic array covering an area of about 10 km by 10 km. We demonstrate that SSRn and SSRh computed below 1 Hz are able to reproduce amplification factors provided by SSR. By contrast, at frequencies higher than 1 Hz, SSRn strongly deviates from SSR. SSRh shows closer results to SSR but presents a dependence to the choice of the local site reference, thereby questioning the possibility to use SSRh blindly to estimate the local amplification in our context. These discrepancies reflect the impact of local anthropogenic sources. We therefore introduced a two-step workflow to mitigate the influence of local sources. The first step is to define a characteristic time window to optimally isolate significant transient signals. The second step consists in selecting the time segments that do not contain these transients with a clustering-based approach. By applying this workflow, we were able to remove some strong anthropogenic transient signals likely to be generated by local sources at some sites and therefore to locally improve the amplification assessment through noise-based spectral ratios. However, stationary sources, whose impact cannot be removed through our procedure, remain a major issue. This study provides some insights into the application of SSRn and SSRh in noisy industrialized areas, especially regarding the impact of local noise sources. It illustrates the difficulty of having a procedure for mitigating the impact of these sources that is efficient everywhere inside such a complex anthropized environment, where different types of sources (including stationary sources) cohabit.

Shallow sedimentary layers have a strong impact on seismic motion. These so-called site effects may be responsible for dramatic ground motion amplification and increase the duration of shaking when an earthquake occurs. The quantification of such amplification effects for specific sites might be challenging to carry out in low-to-moderate seismicity regions where moderate to large earthquakes have long return periods. Therefore, methods based on background ambient noise might be of great interest for these areas. In this study, we investigate the potential of ambient noise in ground motion amplification assessment through SSRn (noise-based Standard Spectral Ratio) and SSRh (hybrid Standard Spectral Ratio, Perron et al., 2018) computation. We continuously recorded ambient noise from February to March 2020 on a 400-sensor seismic array covering an area of about 10 x 10 km in the Tricastin industrial region (French Rhone Valley) where critical facilities are located. This area is located on a very elongated valley, filled with Pliocene sediments (sands and clays), that was dug during the Messinian Salinity Crisis in Cretaceous sandstones and limestones. The strong lithological contrast between the sedimentary filling and the bedrock, as well as the valley's incised geometry, is prone to generate strong and complicated site effects. Previous studies have shown that SSRn is not able to reproduce earthquake-based SSR amplification factor for frequencies higher than 1 Hz. This disagreement may be explained by the influence of local noise sources. Here, we introduce an approach to mitigate the influence of strong local sources in SSRn and SSRh. Our workflow relies on a clustering algorithm to select the Fourier Amplitude Spectrum (FAS) used in the SSRn and SSRh computation. By applying this method, we were able to remove strong anthropic transient signals at some sites and therefore improve the amplification assessment above 1Hz through the SSRn and SSRh. However, half part of the array is located nearby permanent anthropic sources that remain a major issue in quantifying the amplification at the scale of the valley. This study provides some insights into the conditions of applications of SSRn and SSRh in noisy industrialized environments.

. On November 11, 2019, a M w 4.9 earthquake hit the region close to Montelimar (lower Rhône Valley, France), on the eastern margin of the Massif Central close to the external part of the Alps. Occuring in a moderate seismicity area, this earthquake is remarkable for its very shallow focal depth (between 1 and 3 km), its magnitude, and the moderate to large damages it produced in several villages. InSAR interferograms indicated a shallow rupture about 4 km long reaching the surface and the reactivation of the ancient NE–SW La Rouvière normal fault in reverse faulting in agreement with the present-day E–W compressional tectonics. The peculiarity of this earthquake together with a poor coverage of the epicentral region by permanent seismological and geodetic stations triggered the mobilisation of the French post-seismic unit and the broad French scientific community from various institutions, with the deployment of geophysical instruments (seismological and geodesic stations), geological field surveys, and field evaluation of the intensity of the earthquake. Within 7 days after the mainshock, 47 seismological stations were deployed in the epicentral area to improve the Le Teil aftershocks locations relative to the French permanent seismological network (RESIF), monitorthetemporalandspatialevolutionofmicroearthquakesclosetothefaultplaneandtemporal evolutionoftheseismicresponseof3damagedhistoricalbuildings,andtostudysuspectedsitee ff ects and their influence in the distribution of seismic damage. This seismological dataset, completed by data owned by di ff erent institutions, was integrated in a homogeneous archive and distributed through FDSN web services by the RESIF data center. This dataset, together with observations of surface rupture evidences, geologic, geodetic and satellite data, will help to unravel the causes and rupture mechanism of this earthquake, and contribute to account in seismic hazard assessment for earthquakesalongthemajorregionalCévennefaultsysteminacontextofpresent-daycompressional tectonics.

Superficial geological layers can strongly modify the surface ground motion induced by an earthquake. These so-called site effects are highly variable from one site to another and still difficult to quantify for complex geological configurations. That is why site-specific studies can greatly contribute to improve the hazard prediction at a specific site. However, site-specific studies have historically been considered difficult to carry out in low-to-moderate seismicity regions. We present here seismological datasets acquired in the framework of the French-German dense array for seismic site effect estimation project in the heavily industrialized area surrounding the French Tricastin Nuclear Site (TNS). TNS is located above an ancient canyon dug by the Rhone River during the Messinian period. The strong lithological contrast between the sedimentary fill of the canyon and the substratum, as well as its expected confined geometry make this canyon a good candidate for generating site effects that are variable on short spatial scales. To investigate the impact of this geological structure on the seismic motion, we conducted complementary seismic campaigns in the area. The first main campaign consisted of deploying 400 nodes over a 10 x 10 km area for one month and aimed at recording the seismic ambient noise. A second seismic campaign involved the deployment of 49 broadband stations over the same area for more than eight months. This complementary campaign aimed at recording the seismicity (including local, regional, and teleseismic events). These different designs allowed us to target a variety of seismic data at different spatial and temporal scales. Beyond the interest for local operational seismic hazard applications, these datasets may be valuable for studying seismic wave propagation within complex kilometer-scale sedimentary structures. In this article, we present the deployment designs as well as initial analyses to provide information on the characteristics and the overall quality of the data acquired to future users.

The so-called site effects caused by superficial geological layers may be responsible for strong ground motion amplification in certain configurations. We focus here on the industrialized Tricastin area, in the French Rhône valley, where a nuclear site is located. This area lies above an ancient Rhône Canyon whose lithology and geometry make it prone to site effects. This study presents preliminary measurements to investigate the local seismic amplification. We deployed three seismic stations in the area for several months: two stations were located above the canyon, the third one was located on a nearby reference rock site. The recorded seismicity was analysed using the Standard Spectral Ratio technique (SSR). The estimated amplification from weak motions reaches a value of 6 for some frequencies. These first results confirm the possibility of estimating seismic amplification using earthquakes recorded for less than one year, in this highly anthropogenic and industrialized environment, despite the local low-to-moderate level of seismicity. Noise-based SSR, that presents an obvious interest in such seismic context, shows also promising results in the area. To complement this empirical approach, we estimated the amplification using 1D wave propagation modelling. This numerical estimate is based on shear wave velocity profiles resulting from geophysical characterization campaigns. Comparison of the two approaches at low frequency, where numerical estimate is considered as the most representative, tends to suggest that edge-generated surface waves may have a strong influence in the local seismic response. This interpretation will be further investigated in the future.

On November 11, 2019, a M$_{\mathrm{w}}$ 4.9 earthquake hit the region close to Montelimar (lower Rhône Valley, France), on the eastern margin of the Massif Central close to the external part of the Alps. Occuring in a moderate seismicity area, this earthquake is remarkable for its very shallow focal depth (between 1 and 3 km), its magnitude, and the moderate to large damages it produced in several villages. InSAR interferograms indicated a shallow rupture about 4 km long reaching the surface and the reactivation of the ancient NE–SW La Rouvière normal fault in reverse faulting in agreement with the present-day E–W compressional tectonics. The peculiarity of this earthquake together with a poor coverage of the epicentral region by permanent seismological and geodetic stations triggered the mobilisation of the French post-seismic unit and the broad French scientific community from various institutions, with the deployment of geophysical instruments (seismological and geodesic stations), geological field surveys, and field evaluation of the intensity of the earthquake. Within 7 days after the mainshock, 47 seismological stations were deployed in the epicentral area to improve the Le Teil aftershocks locations relative to the French permanent seismological network (RESIF), monitor the temporal and spatial evolution of microearthquakes close to the fault plane and temporal evolution of the seismic response of 3 damaged historical buildings, and to study suspected site effects and their influence in the distribution of seismic damage. This seismological dataset, completed by data owned by different institutions, was integrated in a homogeneous archive and distributed through FDSN web services by the RESIF data center. This dataset, together with observations of surface rupture evidences, geologic, geodetic and satellite data, will help to unravel the causes and rupture mechanism of this earthquake, and contribute to account in seismic hazard assessment for earthquakes along the major regional Cévenne fault system in a context of present-day compressional tectonics.

The horizontal to vertical spectral ratio (HVSR) of seismic ambient noise has been proven to be a fast and efficient method for characterizing the 1-D resonance frequency of the local subsurface in a practical framework. Over the last decades, theories have been developed in order to extend the exploitation of HVSR beside the frequency of its first peak, notably the diffuse field assumption (DFA) which links the HVSR to the Green's function of the local medium assuming the diffuseness of the seismic ambient noise wavefield. However, the underlying assumption of the seismic ambient noise being a diffuse, equipartitioned field may not be satisfied under certain circumstances. In order to exploit the contribution of scattering in forging diffuse wave fields, we leverage the advantages of coda waves and present a novel procedure for computing the HVSR, using the coda part of ambient noise correlations. We applied this technique to data gathered at the plio-quaternary sedimentary basin of Argostoli, Greece. Results on this data set show the potential of the method to improve the temporal stability of the HVSR measurements compared to the classical computation, and the fit with the theoretical HVSR curve derived from the DFA theory. These results suggest that this procedure could help in extracting physical information from the HVSR and thus could lead to an extended use of these measurements to characterize the mechanical properties of the medium.

Site-effect assessments performed through earthquake-based approaches, such as the standard spectral ratio (SSR), require good quality records of numerous earthquakes. In contrast, the use of ambient noise appears to be an attractive solution for ease and rapid computation of site responses with sufficient spatial resolution (microzonation), especially in low seismicity areas. Two main approaches are tested here: the horizontal-to-vertical spectral ratio (HVSR) and the noise-based SSR (SSRn). The HVSR uses the relative amplitude of the horizontal and vertical components of the ambient noise. Instead, the SSRn defines the spectral ratio between the seismic noise recorded simultaneously at a site and at a rock reference station, similar to earthquake-based SSR. While the HVSR is currently used in hundreds of site-specific studies, the SSRn approach has been gradually abandoned since the 1990s. In this study, we compare the results obtain from these two approaches with those of earthquake-based SSR. This comparison is carried out for two sedimentary basins, in Provence (southeastern France) and in Argostoli (western Greece). In agreement with the literature, the HVSR does not provide more than the fundamental resonance frequency of the site (f(0)). The SSRn leads to overestimation of the SSR amplification factors for frequencies higher than the minimal f(0) of the basin (f(0min)). This discrepancy between SSRn and SSR is discussed, and appears to be mainly dependent on the local geological configuration. We thus introduce the hybrid standard spectral ratio (SSRh) approach, which aims to improve upon the SSRn by adding an intermediate station inside the basin for which the SSR is known. This station is used in turn as a local reference inside the basin for the SSRn computation. The SSRh provides site transfer functions very similar to those of the SSR, in a broad frequency range. Based on these results, the SSRn (or SSRh) should be further tested and should receive renewed attention for microzonation inside sedimentary basins.

The Gofar transform fault (GTF), 4°S on the East Pacific Rise, can generate Mw 5.5–6 earthquakes quasiperiodically on some specific patches that are separated by stationary rupture barriers. Small earthquakes along strike show a clear spatial and temporal evolution. To better understand the cause of the observed behaviors of large and small earthquakes, we have determined high-resolution earthquake locations within a period of one year covering the 2008 Mw 6.0 (M6) earthquake, as well as Vp, Vs, and Vp/Vs models along the westernmost segment of the GTF, using a well recorded ocean bottom seismograph dataset and a new Vp/Vs model consistency-constrained double-difference tomography method. Compared to the previous P-wave tomography study in this area, the use of a new automatic arrival picking algorithm significantly improves the accuracy of S-wave arrival times, thereby allowing for the inversion of Vs and Vp/Vs models in addition to Vp model. High-precision waveform cross-correlation differential times are also used. The tomographic Vp/Vs model reveals strong structural variations at multiple scales along the fault, which likely control the behaviors of large and small earthquakes. The M6 mainshock is generated within a specific ∼8-km-long fault patch composed of intact rocks. By contrast, multiple fluid-filled damaged zones on both sides of this asperity are imaged and have varying size which is suggested to be critical in their ability of stopping ∼M6 ruptures. High-resolution earthquake relocations and velocity models also indicate that the occurrence of small earthquakes is also correlated with structural variations. Combined with previous studies, our results further suggest that strong structural variations control the fault mechanics and earthquake behavior along the GTF.

Earthquakes deep in the continental lithosphere are rare and hard to interpret in our current understanding of temperature control on brittle failure. The recent lithospheric mantle earthquake with a moment magnitude of 4.8 at a depth of ~75 km in the Wyoming Craton was exceptionally well recorded and thus enabled us to probe the cause of these unusual earthquakes. On the basis of complete earthquake energy balance estimates using broadband waveforms and temperature estimates using surface heat flow and shear wave velocities, we argue that this earthquake occurred in response to ductile deformation at temperatures above 750°C. The high stress drop, low rupture velocity, and low radiation efficiency are all consistent with a dissipative mechanism. Our results imply that earthquake nucleation in the lithospheric mantle is not exclusively limited to the brittle regime; weakening mechanisms in the ductile regime can allow earthquakes to initiate and propagate. This finding has significant implications for understanding deep earthquake rupture mechanics and rheology of the continental lithosphere.

The aim of this study is to improve the temporal resolution of seismic wave velocity variations measured using ambient noise correlations. We first reproduce the result obtained by Chen et al. using a network of 21 broad-band stations ideally located around the fault system activated during the Wenchuan earthquake.We measure a velocity drop of 0.07 per cent that was associated with the main shock, with a temporal resolution of 30 days. To determine whether this velocity drop is co-seismic or post-seismic, we attempt to increase the temporal resolution of our observations. By taking advantage of the properties of the curvelet transform, we increase the signal-to-noise ratio of the daily correlations computed between each station pair. It is then possible to measure the velocity drop associated with the Wenchuan earthquake with a temporal resolution of 1 day. This shows that the velocity drop started on 2008 May 12, which was the day of the earthquake, and the velocity reached its lowest value 2 days after the main shock. Moreover, there was a second velocity drop on 2008 May 27, which might relate to strong aftershocks.

We locate temporal changes of seismic wave speed and scattering properties of the crust associated with the 2008 Mw7.9 Wenchuan earthquake. To that end, we analyze ambient seismic noise records from 2007 until the end of 2008 in the 1–3 s period band and in the 12–20 s period band measured in a region that covers the southern two thirds of the fault activated during the earthquake. To locate the changes, we use a refined imaging procedure based on the sensitivity of scattered waves to weak perturbations. This inverse method uses the radiative transfer approximation to describe the intensity of the noise correlation coda. Our results show that the largest structural changes are observed northeast of the Wenchuan epicenter, around the fault zone where aftershocks are distributed. The spatial and temporal characteristics of the behavior of the crust at depth around the Wenchuan earthquake suggest a postseismic and geology‐dependent signature of the middle crust. We also observe a clear seasonal signature within the Sichuan basin at depth that we relate to the heavy rainfalls during the monsoon season. This seasonality can also be observed at the surface in form of a waveform decoherence. Our study highlights that seismic velocity changes, and waveform decoherence are independent measurements that are sensitive to different parameters and can thus give complementary information. In our study we also point out the value of studying the evolution of the changes at different times in the coda.