3D crustal structure of the Ligurian Basin revealed by surface wave tomography using ocean bottom seismometer data
- 1GEOMAR Helmholtz Centre for Ocean Research Kiel, 24148 Kiel, Germany
- 2Institute of Geosciences, Kiel University, 24118 Kiel, Germany
- 3Laboratoire de Géosciences Marines, Institut de Physique du Globe de Paris, Paris, 75238 CEDEX 5, France
- A full list of authors appears at the end of the paper.
Correspondence: Felix N. Wolf (email@example.com)
The Liguro-Provençal basin was formed as a back-arc basin of the retreating Calabrian–Apennines subduction zone during the Oligocene and Miocene. The resulting rotation of the Corsica–Sardinia block is associated with rifting, shaping the Ligurian Basin. It is still debated whether oceanic or atypical oceanic crust was formed or if the crust is continental and experienced extreme thinning during the opening of the basin. We perform ambient noise tomography, also taking into account teleseismic events, using an amphibious network of seismic stations, including 22 broadband ocean bottom seismometers (OBSs), to investigate the lithospheric structure of the Ligurian Basin. The instruments were installed in the Ligurian Basin for 8 months between June 2017 and February 2018 as part of the AlpArray seismic network. Because of additional noise sources in the ocean, OBS data are rarely used for ambient noise studies. However, we carefully pre-process the data, including corrections for instrument tilt and seafloor compliance and excluding higher modes of the ambient-noise Rayleigh waves. We calculate daily cross-correlation functions for the AlpArray OBS array and surrounding land stations. We also correlate short time windows that include teleseismic earthquakes, allowing us to derive surface wave group velocities for longer periods than using ambient noise only. We obtain group velocity maps by inverting Green's functions derived from the cross-correlation of ambient noise and teleseismic events, respectively. We then used the resulting 3D group velocity information to calculate 1D depth inversions for S-wave velocities. The group velocity and shear-wave velocity results compare well to existing large-scale studies that partly include the study area. In onshore France, we observe a high-velocity area beneath the Argentera Massif, roughly 10 km below sea level. We interpret this as the root of the Argentera Massif. Our results add spatial resolution to known seismic velocities in the Ligurian Basin, thereby augmenting existing seismic profiles. In agreement with existing seismic studies, our shear-wave velocity maps indicate a deepening of the Moho from 12 km at the south-western basin centre to 20–25 km at the Ligurian coast in the north-east and over 30 km at the Provençal coast. The maps also indicate that the south-western and north-eastern Ligurian Basin are structurally separate. The lack of high crustal ratios beneath the south-western part of the Ligurian Basin preclude mantle serpentinisation there.
The Ligurian Basin is a marginal basin located in the north-western Mediterranean Sea at the transition from the Alpine orogen to the Apennine system (Fig. 1). It formed as a back-arc basin to the south-eastward trench retreat of the Apennines–Calabrian subduction zone during the late Oligocene and Miocene (Gueguen et al., 1998; Rollet et al., 2002). Rifting in the Liguro-Provençal basin initiated about 32 Myr ago (Jolivet et al., 2015). From 21 Ma, the rifting was followed by a counter-clockwise rotation of the Corsica–Sardinia block by approximately 30∘ (e.g. Vigliotti and Langenheim, 1995; Jolivet and Faccenna, 2000; Rollet et al., 2002; Speranza et al., 2002; Schettino and Turco, 2006). Gattacceca et al. (2007) estimate a rotation of 45∘, based on geochronological investigations of Miocene volcanic sequences in Sardinia. Le Breton et al. (2017) describe a total amount of counter-clockwise rotation of the Corsica–Sardinia block by at least 53∘ during the last 35 Myr. At the end of the Burdigalian Age (about 16–18 Ma), the Corsica–Sardinia block was stranded between the Apennines and the European margin in southern France (Rosenbaum et al., 2002). The opening of the Ligurian Basin terminated, while the roll-back of the Calabrian subduction zone continued and initiated the opening of the Tyrrhenian Sea (e.g. Faccenna et al., 2001). Today, the Ligurian Basin is 150–225 km wide (Dannowski et al., 2020), broadening from the north-east to the south-west. The continental margin is narrow (10–20 km) and steep along the Ligurian coast (Finetti et al., 2005) and broader (20–50 km) on the Corsican side (e.g. Rollet et al., 2002). The marine bedrock is covered by a sedimentary layer (e.g. Schettino and Turco, 2006) of variable thickness: less than 3 km thick near the Tuscany coast, increasing towards the south-west to a thickness of up to 8 km in offshore Marseille. Rollet et al. (2002) identify several areas of magmatic intrusions related to three phases of calcalkaline and alkaline volcanism. The first is linked to the opening of the basin, the second links to the transition of the Calabrian subduction zone to the Tyrrhenian Sea, and the third mainly occurred north of Corsica and in the Gulf of Genova (12–11 Ma).
The crust–mantle boundary is well defined along several seismic lines (detailed overview in Dannowski et al., 2020) but otherwise only estimated in parts from large-scale surface wave studies (e.g. Molinari et al., 2015b; Kästle et al., 2018; Lu et al., 2018). Parallel to the basin, Dannowski et al. (2020) explain the satellite-derived free-air anomaly (Sandwell et al., 2014) by gravity modelling along their refraction seismic line (Fig. 1). Dannowski et al. (2020) also include the directly connecting wide-angle reflection seismic line by Makris et al. (1999, Fig. 1). Both seismic and gravity modelling reveal similar values for the Moho depth, showing a gradual thickening of continental crust towards the north-east. At the south-western end of the seismic refraction profile, the Moho is about 12 km deep. It gradually deepens towards the north-east, reaching a depth of 22 km close to the Italian coast. Contrucci et al. (2001) estimated the Moho depth along the multichannel seismic line LISA01 (Fig. 1). They observe a decrease in Moho depth from 18 km at the Corsican margin to 13 km in the basin centre and an increase to over 30 km towards the Provençal coast. This variability is supported by the surface-wave-derived Moho map of Kästle et al. (2018), showing an increasing Moho depth from the Ligurian Basin (<20 km) towards the coast (>25 km).
Many studies addressed whether continental crust was extremely thinned during the rifting or if oceanic spreading occurred and formed oceanic crust in the basin centre. Several authors (Rollet et al., 2002; Gailler et al., 2009; Jolivet et al., 2015) propose an area of atypical oceanic crust, characterised as being very thin (<4 km) and showing complex magnetic anomalies that cannot be correlated to isochrons (e.g. Rollet et al., 2002; Schettino and Turco, 2006), in the basin centre. This area is located next to a broad transition zone towards continental crust at the basin's edges. Based on a recent refraction seismic study, Dannowski et al. (2020) propose that seafloor spreading did not occur during the formation of the Ligurian Basin. They show that beneath the south-western part of the basin, the continental crust thins and possibly gives way to partly serpentinised mantle lying directly beneath an up to 7 km thick sediment cover. Schettino and Turco (2006) find a similar sediment thickness based on a joint interpretation of magnetic and seismic data.
Another open question relates to the location of the prolongation of the Alpine front. It is well defined in onshore France and Corsica, but it remains unclear if and where the connection of these parts of the Alpine front is preserved offshore. At the scale of the entire Alpine belt region, land-data-based ambient noise tomographies have been performed by Molinari et al. (2015b), Kästle et al. (2018), and Lu et al. (2018). These studies revealed the onshore geometry but did not cover the Ligurian Basin. Guerin et al. (2020) conducted an ambient noise surface-wave tomography study along the south-western Alps and a small part of the Ligurian margin using five ocean bottom seismometers (OBSs) and two offshore cabled seismometers close to Nice. They identify a low-velocity zone offshore Nice that is linked to salt and evaporite deposits.
To better understand the present-day crustal velocity structure and its implications on the evolution of the Ligurian Basin, we use a unique amphibious seismic network covering the entire Ligurian Basin and adjacent coastal areas, providing high-resolution group velocity maps and a three-dimensional shear velocity model.
A network of 22 broadband ocean-bottom seismometers (OBSs) was installed jointly by the Institut de physique du globe de Paris (IPGP, Paris, France), the Institut des Sciences de la Terre (ISTerre, Grenoble, France) and GEOMAR Helmholtz Centre for Ocean Research Kiel (Kiel, Germany) (Fig. 1) to investigate the velocity structure of the crust and upper mantle beneath the Ligurian Basin. The AlpArray OBS array is the offshore component of the AlpArray seismic network (Hetényi et al., 2018). The instruments were deployed from the RV Pourquoi Pas? in June 2017 and were recovered in February 2018 by RV Maria S. Merian.
The AlpArray OBS network consisted of 6 French OBSs (BBOBSs) and 16 OBSs (LOBSTER type: “Longterm OBS for Tsunami and Earthquake Research”) provided by the German instrument pool (DEPAS, Schmidt-Ausch and Haberland, 2017). The BBOBSs were equipped with three-component Nanometrics Trillium 240 broadband seismometers with a lower corner period of 240 s and a differential pressure gauge (DPG) designed by the Scripps Institution of Oceanography (Cox et al., 1984). The sampling rate of the installed LCHEAPO recorder was 62.5 Hz. The DEPAS OBSs were equipped with Trillium compact seismometers by Nanometrics Inc. with a lower corner period of 120 s and HTI-01-PCA hydrophones from High Tech Inc. The sampling rate of the K.U.M. 6D6 recorder was 250 Hz (K.U.M. Umwelt and Meerestechnik Kiel GmbH).
The instrument clocks were synchronised with GPS time before deployment and after recovery to reveal any internal clock drift and apply a linear clock drift correction. We calculated every station's probabilistic power spectral densities (PPSDs) (McNamara and Buland, 2004). The lowest spectral levels on the vertical seismometer components fall in between the mean minimum and maximum noise levels for land stations (Peterson, 1993) for both the German (Fig. S1a–b in the Supplement) and French OBSs (Fig. S1c–d). Regarding the pressure component, the French DPGs yield high-quality data (−38 to 40 dB) while the HTI hydrophones have a range of −20 to 40 dB with a lesser resolution for periods larger than 10 s (Fig. S1a, c). These results are comparable to similar instrument set-ups (Stähler et al., 2016) used during the RHUM-RUM OBS experiment in the Indian Ocean. To resolve the entire structure of the Ligurian Basin and the surrounding onshore areas, we incorporated 16 temporary and 42 permanent land stations in our analysis (Table S1 in the Supplement).
The ambient noise technique was developed during the last 20 years (e.g. Lobkis and Weaver, 2001; Wapenaar et al., 2010a, b) and is based on the concept of Aki (1957) regarding the spectra of stationary stochastic waves. Ambient noise techniques exploit the “noise” of long-term recordings as the desired signal. This part of the measured signal includes, for example, anthropogenic noise, microseismic signals from ocean-coast interactions, and highly scattered waves of teleseismic origin (Campillo and Paul, 2003; Campillo and Roux, 2014). Given a continuous measurement and uniformly distributed noise sources, the cross-correlation of recordings of two stations is used as the empirical Green's function representing the subsurface response to a wave propagating from one station to the other. These empirical Green's functions from different station pairs are used to invert for two-dimensional (2D) group velocity maps, 1D velocity–depth profiles or 3D velocity distribution maps.
Although the technique is well established for land data (e.g. Barmin et al., 2001; Campillo and Paul, 2003; Shapiro et al., 2005; Prieto et al., 2009; Goutorbe et al., 2015; Kästle et al., 2019), it is not yet used regularly for ambient noise analysis on ocean bottom seismometer data. Previous studies show that ambient noise can be calculated using OBS data (e.g. Harmon et al., 2007, 2012; Takeo et al., 2014; Dewangan et al., 2018). However, compared to land stations, seismic recordings on OBSs contain less anthropogenic noise but other additional noise sources like tilt and compliance noise (Crawford et al., 1998; Webb, 1998; Bell et al., 2015).
3.1 Pre-processing – tilt and compliance correction
Adimah and Padhy (2020) showed that reducing tilt and compliance noise before running the cross-correlation proves beneficial, as tilt and compliance noise are not part of the useful signal. Therefore, we pre-process the OBS data as proposed by Crawford and Webb (2000). First, we cut the continuous OBS recordings into daily files and resample the data at 1 Hz. Next, we remove tilt and compliance noise. Tilt noise is introduced by a slight inclination of the instrument, causing horizontal movement to appear on the vertical component (Crawford and Webb, 2000). Although the instruments level themselves to an accuracy of ±0.5∘ (LOBSTER type) and ±5∘ (BBOBS), respectively, the remaining tilt is sufficient to create tilt noise. The tilt of the instrument can be caused by processes such as ocean bottom currents. On the other hand, compliance is a signal generated by ocean infragravity waves introducing pressure fluctuations that cause micrometre-scale deformation of the seafloor (Webb and Crawford, 1999). The variations of the gravitational forces of the water column, the deformation of the seafloor itself, and the variation caused in the distance of the OBS to the Earth's gravitational centre all introduce changes to the measured acceleration (Crawford et al., 1998). Thus, compliance increases vertical acceleration noise level by 10 to 25 dB for 30–100 s (Webb and Crawford, 1999).
To correct for tilt and compliance noise, we applied the procedure described in Crawford and Webb (2000) and Bell et al. (2015). First, we calculate a transfer function between the vertical seismometer component and the hydrophone component. Next, we subtract the coherent part of the signal (in this case: compliance) from the vertical seismometer component. We also corrected both horizontal components for compliance before removing tilt noise (Crawford and Webb, 2000). Subsequently, the same routine is used to remove the coherent signal between the vertical component and each horizontal component to remove tilt noise. Thus, we calculate the transfer functions between the vertical component and each of the horizontal components. Finally, we obtain a vertical component corrected for tilt and compliance. The order in which the components are corrected is interchangeable. The land station recordings were not corrected for tilt and compliance noise but are also resampled to 1 Hz.
3.2 Ambient noise technique – cross-correlation and mode identification
We use the tilt- and compliance-corrected daily files to estimate cross-correlation functions (CCFs) for every vertical component OBS–OBS and OBS–land station pair (Bensen et al., 2007). Additionally, we calculate CCFS for all land–land pairs for the land stations A317A, ARBF, and DIX (see Fig. 1) and CCFs for all combinations of 20 selected land stations (namely AJAC, BLAF, BOB, BSTF, CALF, CARD, EILF, ENR, GBOS, IMI, ISO, MSSA, PCP, PLMA, ROTM, SAOF SMPL, TRBF, TURF, and VLC) to increase the ray coverage onshore (Fig. 2). The cross-correlation is calculated day-wise for every station pair. Afterwards, we stack the single-day CCFs to estimate one CCF per station pair.
In addition to ambient noise cross-correlations, we correlate time windows (45 min long) that include strong teleseismic events using the two-station method (e.g. Meier et al., 2004; Boschi et al., 2013; Tonegawa et al., 2020). We only use station pairs for which the stations' azimuth equals the great circle from the event to within ±7∘. The correlation window, starting at the origin time of the event, is dominated by the earthquake signals. The further processing is identical to correlating ambient noise day files but is performed for longer periods (20 to 90 s).
3.2.2 Group velocity dispersion curves – fundamental mode and higher modes
To estimate group velocity dispersion curves, we apply the multiple filter technique (MFT) (Dziewonski et al., 1969). A narrow bandpass filter is applied to the CCFs to derive the velocity for a distinct period from the maximum correlation (e.g. Meier et al., 2004).
Extra care has to be taken when picking the dispersion curves since for some station pairs the first higher mode has stronger amplitudes than the fundamental mode. Different modes have different sensitivity kernels (e.g. Harmon et al., 2007), and, unfortunately, our tomography program cannot process input data from more than one mode at a time. Therefore, we picked manually by comparing each ray path to the theoretical dispersion curves predicted from a 2D model of the research area (Fig. S2) that includes results from Makris et al. (1999), Gailler et al. (2009), and Dannowski et al. (2020).
First, we picked the maximum signal on all dispersion curves. Station pairs showing no detectable maximum were excluded. After comparing the group velocities with synthetic dispersion curves, we excluded about 100 station pairs that showed velocities likely related to higher modes (Fig. 3b). Higher modes were mainly observed for ray paths in the southern part of the Ligurian Basin and parallel to the basin axis (Fig. 4). The origin could be layers in which the first higher-mode couples more strongly than the fundamental mode, as previously observed by Takeo et al. (2014) for CCFs from OBSs in the NW Pacific. During the MFT revision, we did not observe a degradation of the signal depending on the station distance.
Identifying and rejecting the higher-mode dispersion curves resulted in 1342 dispersion curves for the fundamental mode from ambient noise CCFs and 1963 additional dispersion curves from teleseismic CCFs (Fig. 5) used for further analysis steps. We use the ambient noise CCFs to derive dispersion curves for periods from 4–15 s. The CCFs from the correlation of teleseismic events were used to derive dispersion curves from 20–90 s (Fig. 5). The dispersion curves' frequency bands are complementary and provide a bandwidth ranging from 4–90 s.
3.3 Surface wave tomography for group velocities of ambient noise data and teleseismic data
We use the fast marching surface tomography method (FMST, Rawlinson and Sambridge, 2004, 2005) to derive 2D Rayleigh group velocity maps (Figs. 6 and S11) from the picked dispersion curves. FMST inverts for 2D map slices of group velocities for a given period. The forward prediction of travel times is achieved using the fast marching method (Sethian, 1996; Sethian and Popovici, 1999), a finite-difference solution of the eikonal equation. The inversion scheme is non-linear and repeated iteratively. Prior to the inversion, we deleted data outside the allowed velocity range: 0.5–3.6 km s−1 for periods <10 s, 1.0–4.0 km s−1 for periods between 10–20 s, and 2.0–6.0 km s−1 for periods 20 s and larger. We derived these thresholds based on the seismic velocity model of an active seismic refraction profile in the centre of the Ligurian Basin (Dannowski et al., 2020). The damping parameter for every period was estimated from L curves (e.g. Hansen, 1992; Fig. S2). The smoothing parameter was chosen visually depending on the resolution of the inversion (see Table S3 for inversion parameters) and the result of the checkerboard tests (Fig. 7). The input error is based on the picking error and linearly increases with the increasing periods from ± 0.75–2.0 s. We use homogeneous starting models with period-dependent velocities (Table S2). These are based on a group velocity model calculated from the seismic refraction line by Dannowski et al. (2020). The inversion grid consists of 28 × 35 nodes, resulting in one node every 18 km for both N–S and W–E directions. Due to grid refinements, the output grids consist of 406 × 511 grid points with a spacing of 0.0157∘ × 0.0136∘ or 1.5 km × 1.23 km. We calculate 2D group velocity maps for 5, 6, 7, 8, 9, 10, 12, and 15 s from ambient noise CCFs and 20 to 90 s in 10 s steps from teleseismic CCFs.
In the initial group velocity maps, we observed a low-velocity area west of Marseille associated with station ARBF. We ran a tomography with all ray paths from ARBF excluded, and the result did not show the low-velocity area. Since the station is positioned on sediments in the outer Rhône delta, we assume the low-velocity zone to be caused by a locally “slower” subsurface. We decided to exclude the station from our dataset to prevent the smearing of local low velocities into the group velocity maps. A similar low-velocity zone was observed close to station A430A, which was excluded as well.
Next, we use these initial group velocity maps to calculate residuals between the model input and the tomography output. We evaluate the residuals and keep those station pairs corresponding to 1.28 standard deviations (σ; 80 % of all pairs) for periods of 5 to 15 s. For longer periods, we observe smaller residuals and therefore keep 90 % of the station pairs (1.64σ). Then, we recalculate the 2D tomographies with the updated dataset (see Table S2 for final numbers) to create group velocity maps from ambient noise CCFs.
3.3.1 1D depth inversion
To remove effects of the highly variable topography and bathymetry, we invert for 1D shear-velocity–depth profiles using the iterative, weighted inversion code from Herrmann (2013). We produce one 1D vS–depth profile for every 10th grid point, corresponding to one profile every 12.3 km. To account for the non-uniqueness of the solution (Foti et al., 2018), we set up a starting model with fixed layers (Table S4) based on the vP velocities from Dannowski et al. (2020) for up to 16 km depth and on PREM (Dziewonski and Anderson, 1981) below. Since the dispersion curves for FMST represent a cumulated velocity profile for the subsurface between two stations, it is crucial to correctly parameterise the topography and water column prior to the velocity–depth inversions. To consider the effect of topography and bathymetry, we set up the two uppermost model layers independently: on land, the top layer reaches from the local elevation to the sea level. The second layer reaches from sea level to a depth of 4 km. Offshore, the uppermost layer represents the water column with fixed velocities of vP=1.52 km s−1 and vS=0 km s−1, reaching from the sea surface down to the seafloor, followed by a second layer below reaching from the seafloor down to 4 km depth. Therefore, below 4 km depth, all input models are identical. The layer thicknesses are not varied during the inversion, and the velocity uncertainty is estimated as 2 % of the input group velocity. After parameterisation, we perform iterative 1D depth inversions for vS (Herrmann, 2013). We obtain 1D velocity–depth profiles from the surface to a depth of 30 km.
3.4 Data quality
In general, the OBS stations have noise characteristics comparable to land data (Fig. S1). However, roughly 50 % of all possible CCF combinations do not show a clear correlation of the group velocities and hence were not considered further. Each OBS is part of combinations resulting in high-quality and low-quality CCFs. Similar effects have been observed by Harmon et al. (2012) and Adimah and Padhy (2020). One reason for this may be the variability in station sites. For the OBS, the water depth is highly variable (1133 to 2773 m). Also, the seafloor characteristics and the coupling to the subsurface are most likely very variable due to the varying sediment thickness beneath each site (Schettino and Turco, 2006). Overall, the essential difference between the AlpArray OBS stations and most previous studies (e.g. Harmon et al., 2007; Lin et al., 2016) is the shallow location of AlpArray OBSs in the marginal Ligurian Basin. To our knowledge, this study provides the shallowest OBS water depths used for CCFs, and the shallow water might not prove beneficial for the correlation quality. Five stations are at water depths of 2 km or less (the shallowest station A434A is at 1.1 km depth); none is deeper than 2.8 km. Harmon et al. (2012) estimated CCFs of similar quality using OBS stations at 2.5–3.5 km water depth in offshore Sumatra. Adimah and Padhy (2020) use OBSs in deeper water (14 OBSs in 4.3–5.1 km depth and only 1 OBS at 2.7 km). They observe variations in CCF quality as well, but their overall quality of CCFs is better than for our dataset.
Other reasons for our comparably low CCF quality include the form of the basin itself, for which noise sources are not uniformly distributed, and probably also the highly variable weather conditions in our research area. The Mediterranean Sea lies in a westerly wind system, but especially during winter, mistral events change the flow pattern of regional ocean currents (e.g. Millot and Wald, 1980; André et al., 2005). Moreover, mistral winds might create significant wave heights of 4 m and more (e.g. Pasi et al., 2011). Those temporary changes of the water column and currents alter the pressure on the OBS and the ocean floor and might therefore introduce highly variable noise. Additionally, the land station locations vary in topography and geological settings ranging from sediment basins to Alpine mountains. Nevertheless, overall we estimated more than 3300 high-quality dispersion curves.
3.4.1 Resolution tests
To estimate the resolution of the group velocity maps, we calculated two checkerboard tests for every period (5–90 s) with tiles of 0.4∘ × 0.4∘ and 0.8∘ × 0.8∘, respectively (Figs. 7 and S3–S6). The tiles' deviation from the input velocity was set to ±25 %. Synthetic data are calculated and inverted, using the same set-up as for the picked data. Overall, the resolution is good in the Ligurian Basin and along the northern coast. We defined one polygon where vG is reasonably well defined for all investigated periods. This was necessary to use the group velocity maps as input for the 1D shear-wave velocity inversion. Additionally, we performed a restoration test based on a synthetic 2D model of the research area (Fig. S7a). Different parts of the research area were assigned to distinct group velocity profiles (Fig. S7b). The synthetic group velocity maps for distinct periods are shown in Fig. S8. Based on this model, we calculated synthetic lag times for all station pairs used in the real dataset. This was done by projecting the ray path and estimating a total lag time, based on the proportion that the ray travels through the different areas. We then used this synthetic dataset to calculate group velocity tomographies, using the same settings as for the real dataset (Fig. S9). Comparing these to the synthetic group velocities (Fig. S8) supports our checkerboard test results. The Ligurian Basin itself and the Liguro-Provençal coast are well resolved. We observe some artefacts caused by the ray coverage (e.g. finger-like high-velocity areas in the north and from Corsica to the Italian mainland) that lie outside the interpreted area.
To evaluate the resolution of the 1D-shear-wave inversion, we used the group velocity maps (Fig. S9) to run a synthetic 1D-shear-wave inversion based on the restoration test. For this, we also use the same set-up as for the real dataset. The resulting shear-velocity–depth layers are a good hint of resolvable areas of the lithospheric structure (Fig. S10). Also, we estimate a root-mean-square (RMS) error for every 1D inversion.
4.1 Rayleigh wave group velocity
We show 2D group velocity maps for periods of 5, 8, 12, 20, and 40 s (Fig. 6b–f) accompanied by the tomography input as coloured ray path plots for 8 s (Fig. 6a). The resolvable area, marked by the red polygon in Fig. 6b–f, is determined from the checkerboard tests (Fig. 7); poorly resolved parts are transparent in the group velocity maps. Group velocity maps for all other periods used are shown in Fig. S11. The ray coverage for ambient noise tomography and the cross-correlation of teleseismic events differs (Fig. 2). Still, the resolved area of both datasets covers the Ligurian Basin and adjacent coastal areas (Figs. 7 and S3–S6), whereby the Liguro-Provençal coast is better resolved than the Corsican margin.
Along the Liguro-Provençal coastline, we observe a clear velocity change for periods of 5–12 s (Fig. 6b–d): vG ≅ 1–1.5 km s−1 offshore and vG ≅ 2.5–3 km s−1 onshore for 5 and 8 s, vG ≅ 2–2.5 km s−1 offshore and vG≥2.8 km s−1 onshore for 12 s. For longer periods, this distinction becomes less sharp, and the velocity gradient changes direction. For 20 and 40 s (Fig. 6e, f), vG is approximately 0.5 km s−1 slower onshore compared to the Ligurian Basin. The group velocity maps for periods 20 and 40 s appear more homogenous than those for shorter periods. For 20 s we observe vG = 3–3.5 km s−1; for 40 s it is vG = 3.5–4 km s−1.
The Ligurian Basin appears to be separated into a south-western (labelled SW in Fig. 6c) and a north-eastern part (labelled NE in Fig. 6c) of the Ligurian Basin. The onshore–offshore separation appears less distinct in the north-eastern basin, where the group velocity increases gradually towards the coast. In short periods, the NE part of the basin is faster (NE: 1.5–2.5 km s−1 at 5 s, SW: ∼ 1 km s−1) than the south-western part. At 12 s, the velocity gradient is smaller (NE: ∼ 2.5 km s−1, SW: ∼ 2 km s−1), and at 20 s the gradient vanishes.
Overall, the group velocity increases with increasing period. The velocity gradient is strongest (5 s period: vG≅1 km s−1; 12 s period: vG = 2–3 km s−1) beneath the south-western basin, less strong (vG≅1.75 to 2.5 km s−1) beneath the north-eastern basin, and least strong (vG≅2.5 to 3–3.25 km s−1) beneath the mainland.
4.2 1D shear-wave velocity inversion
We calculated 1D depth inversions for vS based on the group velocity maps (5–90 s) described above. It was a crucial step to remove the topographical effects that result from the amphibious nature of our study area. The average RMS of the 1D inversions is 0.15 km s−1 (Fig. 8i).
4.2.1 Liguro-Provençal coast
At shallow depth, the onshore velocity structure is heterogeneous. At a depth of 6–9 km below sea level (Fig. 8c), we see vS ≅ 2.75–3 km s−1 for the Po Plain, vS≥3.5 km s−1 along the Alpine belt, vS ≅ 2.7–3 km s−1 west of Nice, lower vS directly at the coast, and an increase in vS (vS≅3 km s−1) towards the Maures–Esterel Massif (Fig. 1). Just on the shore of Liguria (Fig. 8c), our results also indicate a narrow band of vS≅3.5 km s−1 at 6–9 km depth accompanied by lower vS≅3.2 km s−1 offshore.
At 9–12 km depth, we observe a high-velocity area north of Nice (dashed circle in Fig. 8d), showing vS≅4.2 km s−1. In other depth layers, this area does not show large velocity differences compared to the surrounding area. At depths of up to 12–15 km, we observe high S-wave velocities beneath the Alpine belt that decrease towards the Rhône delta in the south-west and towards the Po Plain in the north-east (Fig. 8a–e). At larger depths, the high-velocity anomaly disappears, and the velocity field along the coastline gets smoother. At 18–21 km depth (Fig. 8g), we observe large areas of vS=3.5 km s−1. At 21–25 km depth (Fig. 8h), the velocity reaches vS=4 km s−1 locally.
4.2.2 South-western and central Ligurian Basin
In the south-western and central basin (labelled SW in Fig. 8a), the shear-wave velocities in the uppermost 4 km (Fig. 8a) are ∼ 1.5 km s−1. The velocity increases towards the Provençal coast and the Gulf of Lion. The 4–6 km layer (Fig. 8b) shows vS≅2 km s−1 with areas of higher velocity (vS≥3 km s−1) just off the coast of Marseille and north-west of Corsica. Throughout the basin, but mainly along the basin axis, we observe areas of higher S-wave velocity of up to 3.5 km s−1 at 6–9 km depth (Fig. 8c).
With increasing depth, the S-wave velocity increases and at 12 km depth, vS≥4.3 km s−1 is reached locally in the basin centre. These fast areas broaden in the 12–15 km depth slice (Fig. 8e). The S-wave velocity is slower towards the Provençal coast. At 15–18 km depth (Fig. 8f), we observe vS≥4.3 km s−1 along the basin axis of the whole south-western and central Ligurian Basin. However, vS is slower (3.7–4 km s−1) south of Marseille and in the outer Gulf of Lion. At a depth of approximately 21 km (Fig. 8g), vS≥4.3 km s−1 applies to most of the south-western and central basin, except for the aforementioned areas. Close to the Provençal coast, vS≥4.3 km s−1 is reached only in the 21–25 km depth layer (Fig. 8h).
The “fingers” of high vS leading from the basin axis towards the coast east and west of Nice (e.g. Fig. 8g) are probably caused by an insufficient ray coverage of the group velocity tomography in that area. The ray coverage is better just off the shore of Nice. Therefore, we expect a similar vS to that just off the coast of Nice (vS=3.5 km s−1).
4.2.3 North-eastern basin
In the north-eastern basin (labelled NE in Fig. 8a), the shear-wave velocity is higher than in the south-western basin for shallow depths. North of Corsica, vS≅2 km s−1 at up to 4 km depth (Fig. 8a) with higher velocity vS≅2.5 km s−1 close to the Italian coast. The offshore velocity increases to vS ≅ 2.5–3 km s−1 at 4–6 km depth (Fig. 8b) and vS≅3 km s−1 at 6–9 km below the sea surface (Fig. 8c). In both layers, we identify an area of higher velocity north-east of Corsica. This patch shows vS>3 km s−1 in the 4–6 km layer (Fig. 8b) and vS>3.5 km s−1 at 6–9 km depth (Fig. 8c). From 9 km to up to 21 km depth (Fig. 8d–g), the offshore vS increases slowly from approximately 3.5 to 3.8 km s−1. Close to the Italian coast, the velocity gradient direction switches at approximately 12–15 km depth (Fig. 8e). For deeper layers, vS is lower near the Italian coast than towards the south-eastern basin. At 21–25 km depth (Fig. 8h), vS≥4.3 km s−1 accounts for the whole basin, except for a narrow band at the Ligurian coast that shows lower velocities of vS≅4 km s−1.
In the following, we discuss three regions that show differing characteristics in the velocity maps: the Liguro-Provençal coast, the south-western and central Ligurian Basin, and the north-eastern basin. Also, we discuss the proposed offshore prolongation of the Alpine front.
5.1 Liguro-Provençal coast
Along the Liguro-Provençal coast, we can compare our results to existing larger-scale ambient noise studies from Molinari et al. (2015b) and Kästle et al. (2018), as well as a local ambient noise study by Guerin et al. (2020). Guerin et al. (2020) conducted an ambient noise study covering the Provençal coast from Marseille to the Argentera Massif north of Nice (Fig. 1). Guerin et al. (2020) show Rayleigh wave group velocities as coloured ray coverage maps, as we do for the 8 s period (Fig. 6a). In the 8 s period, we observe vG = 2.75–3.2 km s−1 in the coastal area. Guerin et al. (2020) find approximately vG=3 km s−1 (their Fig. 6). In shear-wave velocity maps, Guerin et al. (2020) observe vS ≅ 3–3.5 km s−1 at a depth of 6.4 km (their Fig. 12) along the Provençal coast. This fits our results nicely (Fig. 8c). For shallower depths, Guerin et al. (2020) found that the S-wave velocity increases with depth faster than in our dataset, a feature that is probably controlled by a denser station spacing compared to our study.
At up to 9 km depth (Fig. 8a–c), we observe laterally varying shear-wave velocities on land that we assume to be caused by variations in the geology. At the Rhône delta (Fig. 1), where the sedimentary cover is up to 12 km thick (Le Pichon et al., 2010), we observe vS≅2.7 km s−1 in the layer at 4–6 km depth (Fig. 8b) and vS≅3 km s−1 in the 6–9 km depth range (Fig. 8c). Similarly, the Po Basin has an average sedimentary thickness of 4–5 km (Molinari et al., 2015a) with a shear-wave velocity increasing from vS≅2.5 km s−1 to vS≅3.1 km s−1 at 4–9 km depth. In contrast to the sedimentary basins, we observe higher vS ≅ 3–3.5 km s−1 (4–9 km depth) beneath the Alpine belt, composed of crystalline and metamorphic rocks (e.g. Molinari et al., 2015b). This S-wave variation is most probably caused by the different rock types and structure of the Alpine belt and the sedimentary basins. West of Nice, we observe vS ≅ 2.7–3 km s−1 at 6–9 km depth (Fig. 8c), lower vS directly at the coast, and an increase in vS (vS≅3 km s−1) towards the Maures–Esterel Massif (Fig. 1). These results all compare well to Molinari et al. (2015b, compare their Fig. 8). The large-scale structures also compare nicely to the results of Kästle et al. (2018, their Fig. 9) and Lu et al. (2018, their Fig. 7). In contrast to Molinari et al. (2015b), they observe slightly higher velocities at 10 km depth that are closer to those in our 9–12 km depth slice (Fig. 8d). Both our results and Kästle et al. (2018) indicate a narrow band of vS≅3.5 km s−1 at 6–9 km depth just on the shore of Liguria (Fig. 8c), accompanied by lower vS≅3.2 km s−1 offshore. This observation of a local high-velocity area in onshore Liguria is supported by seismic refraction profiles evaluated by Laubscher et al. (1992). They identified several high-velocity bodies (vP>6 km s−1) linked to ophiolites.
Our results show a high-velocity area (vS≅4.2 km s−1) north of Nice at 9–12 km depth (dashed circle in Fig. 8d), coinciding with a small area of higher velocity in the 10 km depth map of Kästle et al. (2018). This is probably linked to the Argentera Massif (Fig. 1). The massif is composed of crystalline rocks and was identified as a high-velocity area in more shallow depths (vS≅3.4 km s−1 at 6 km depth) by Guerin et al. (2020). The high-velocity anomaly cannot be tracked at greater depth (Fig. 8e). Instead, the velocity field along the coastline gets smoother. At 18–21 km depth (Fig. 8g), we observe large areas of vS=3.5 km s−1. This velocity is similar to Molinari et al. (2015b) and Kästle et al. (2018). We still observe crustal velocities (vS≅4 km s−1) at 21–25 km depth (Fig. 8h), in line with the Moho depth of ∼ 35 km beneath the Liguro-Provençal coast (Kästle et al., 2018).
5.2 South-western and central Ligurian Basin
At shallow depths (Fig. 8a–b), the S-wave velocity is mainly dominated by sediments. We observe vS = 1–1.5 km s−1 in up to 4 km depth (below the sea surface) and vS = 2–2.5 km s−1 at 4–6 km depth. Studies by Schettino and Turco (2006) revealed thick sediment layers in the south-western Ligurian Basin. Off the coast of western Corsica, the sediments are 6–7 km thick, with the maximum thickness of 8 km occurring to the south-west of Marseille (e.g. Schettino and Turco, 2006). This is supported by the findings of Moulin et al. (2015). Their wide-angle reflection seismic data show up to 7.6 km of sediment in the south-eastern Gulf of Lion, thinning to 6.3 km in the Ligurian Basin. Throughout the basin, but mainly along the basin axis, we observe areas of higher S-wave velocity of up to 3.5 km s−1 at 6–9 km depth (Fig. 8c). These are also the areas with the highest RMS error (Fig. 8i). Some of these, e.g. north of Corsica, are in locations where Rollet et al. (2002) observe magmatic anomalies related to magmatic intrusions. We deduce that the velocity gradient is stronger for fast areas along the basin axis than away from the basin axis. The changing gradient is probably caused by the observed thinning of continental crust (Dannowski et al., 2020) and possible exhumation of denser lower crust and upper mantle rock (Gailler et al., 2009; Jolivet et al., 2015) observed further to the south-west. Both would lead to a higher S-wave velocity near the basin axis.
Dannowski et al. (2020) observe a Moho depth of about 12 km in the basin centre, where we observe patches of vS≥4.3 km s−1 along the basin axis (Fig. 8e). Comparing the P-wave velocity of Dannowski et al. (2020) to our S-wave velocity, we calculate a ratio of = 1.74 at the south-western end of their profile. Following Carlson and Miller (2003), this does not indicate mantle serpentinisation. This interpretation is also supported by the local seismicity study of Thorwart et al. (2021). They observe vP=8.1 km s−1 and vS=4.7 km s−1 () in the basin centre roughly 3–4 km below the Moho. The fast area along the axis broadens in the 12–15 km depth slice (Fig. 8f). Linking this observation to the seismic lines (e.g. Jolivet et al., 2015; and Dannowski et al., 2020), our results indicate that the Moho depth increases towards the Provençal and Corsican margins. At larger depths, the velocity maps get more homogeneous. This hints at fewer heterogeneities in the mantle but might also be caused by the decreasing sensitivity of group velocities with increasing periods (e.g. Adimah and Padhy, 2020). For a 5 s period, the group velocity is sensitive to a narrow depth range that peaks at ∼ 5 km. For 20 s, the overall sensitivity is lower and has a much broader range of approximately 10–25 km depth.
In the central Ligurian Basin, Contrucci et al. (2001) investigated a multichannel seismic profile (LISA01) from Antibes, close to Nice, to L'Île Rousse on Corsica (Figs. 1 and 8b–c). They find that the transition from sediments to crust (at vP ∼ 4.8–5 km s−1) is shallow at the Provençal coast (3 km below sea level), deepens towards the basin centre (8 km below sea level), and rapidly shallows again at the Corsican margin (from 5 to 1.5 km below sea level). Also, the salt (Messinian) and sediment (Miocene) layers (vP=3.8 to 5 km s−1) thicken towards the basin centre (Contrucci et al., 2001), where the Moho is ∼ 12–13 km deep. Our vS maps for 4–9 km depth (Fig. 8b–c) show a local velocity high with increasing vS (2.75–3 km s−1) directly off the coast of Nice. Further south-east along the LISA01 profile (Fig. 8b–c, solid black line), the velocity decreases to vS≅2.1 km s−1 at 4–6 km depth and vS ≅ 2.5–2.7 km s−1 at 6–9 km depth. The resolution is poor at the Corsican margin, but vS increases to 3 km s−1 towards Corsica (Fig. 8c). The observed velocity structure fits the findings of Contrucci et al. (2001) nicely, supporting their finding of thicker sediment and salt layers near the basin axis. Comparing vP of the LISA01 profile to our vS gives a ratio of 1.75 for the high-velocity area in offshore Nice and ≅ 2.1 for the sediment layers at the basin axis. Shillington et al. (2007) found similar values for sediments up to 1 km below the seafloor.
Dannowski et al. (2020) suggest that continental crust was (extremely) thinned along their profile, but that no spreading occurred. This is in-line with our results. A possible spreading centre must be located to the south-west. At the Gulf of Lion margin along the south-western edge of our research area, Gailler et al. (2009) interpreted their results as oceanic crust, also observing a transition zone made up of “lower crustal material or mixture of serpentinised upper mantle material with lower crustal material” (Gailler et al., 2009). Later, Jolivet et al. (2015) attributed the shallow high velocities to exhumed lower crust and possibly also partially serpentinised mantle. Therefore, a spreading centre may have been located south-west of our research area, possibly as close as the Gulf of Lion margin.
5.3 North-eastern basin
North-east of the LISA01 profile, the north-eastern basin (labelled NE in Fig. 8a) exhibits different characteristics than the south-western and central Ligurian Basin. For shallow depths of up to 12–15 km, the S-wave velocity is higher in the north-eastern basin compared to the south-western basin. For greater depth, this switches and the north-east is slower. The velocity increase with depth is smaller in the north-east compared to the south-west. Overall, the north-eastern basin is more homogenous than the south-west, and the transition from the basin to onshore Italy is not as sharp as at the Provençal coastline. These observations are supported by large-scale ambient noise studies by Molinari et al. (2015b) and Kästle et al. (2018), who observed a similar velocity distribution.
The sediment thickness map by Schettino and Turco (2006) shows a sediment thickness of 3–4 km in the north-east, increasing to the south-west. The north-eastward thinning of the sediment layer explains the higher vS at shallow depths (Fig. 8b) compared to the south-west. Additionally, Makris et al. (1999) suggest that the sediments' compactness increases to the north-east. Increasing compactness would add to the velocity increase.
The crust–mantle boundary is well defined along seismic profiles across the Ligurian Basin. Our shear-velocity model adds spatial information to these studies, allowing for a broader understanding of the Moho. For the north-eastern basin, we observe vS<4 km s−1 in the 18–21 km layer and mantle-like vS≥4.3 km s−1 in the 21–25 km depth layer. This compares well to the Moho depth of 22 km near the Italian coast, observed by Dannowski et al. (2020). They observe an increasing Moho depth towards the north-east. Linking our 3D shear-wave velocity data to the seismic observations indicates that the Moho depth gradually increases towards the north-east and from the basin axis towards the Provençal coast. Close to the Ligurian coast, we observe mantle-like vS≥4.3 km s−1 in most of the 21–25 km layer, except for a slim band of lower vS at the coast. At the coastline, Kästle et al. (2018) predict a Moho depth of 30–40 km. The apparent thickening of the continental crust towards the north-east is likely related to the position of the rotational pole of the opening of the Ligurian Basin during the Oligocene–Miocene. According to Speranza et al. (2002) and Gattacceca et al. (2007), the rotational pole was located in the north-eastern Ligurian Basin at 43.5∘ N, 9.5∘ E. Therefore, the south-western basin was more extensively opened, and the continental crust was thinned further than in the north-east.
5.4 Alpine front
Rollet et al. (2002) raised the question of whether an offshore prolongation of the Alpine front can be observed in onshore France and onshore Corsica. These authors suggested that the south-western and north-eastern parts of the Ligurian Basin form, respectively, the footwall and hanging wall of the Alpine front. Thus, the Alpine front would be located approximately at the boundary between the north-eastern and south-western crustal domains distinguished in our data (illustrated by the dashed line in Fig. 8e). Dannowski et al. (2020) observe a gradual thickening of the continental crust towards the north-eastern part of the Ligurian Basin. To explain the free-air anomaly derived by Sandwell et al. (2014), they did not need the sharp step that Makris et al. (1999) introduced between Corsica and the Liguro-Provençal coast. In keeping with Dannowski et al. (2020), our spatial shear-wave velocity data do not show a sharp lateral boundary, but a gradual change of the velocity layers. Detection of an offshore Alpine Front is therefore not feasible with the current resolution.
Applying ambient noise techniques and the correlation of teleseismic events to amphibious data results in the first 3D high-resolution seismic group and shear velocity models for the Ligurian Basin. Data processing of the OBS data included correction for tilt and compliance. The dataset differs from most previous ambient noise studies using OBS data. Our stations are comparably shallow, and the fundamental mode is not always the most prominent signal in the marine ray paths. Higher modes are primarily observed in the south-east. Onshore, our results compare well with existing larger-scale ambient noise studies. We reveal a high-velocity area at the Argentera Massif, approximately 10 km below sea level. Offshore, the lithospheric structure in the Ligurian Basin mostly mimics the geometry of the basin. Shear-wave velocity maps indicate a gradual deepening of the Moho from 12–15 km in the south-western basin centre towards 20–25 km in the north-eastern basin and a more rapid deepening from the basin axis to the Provençal coast (>30 km). Based on the low ratios of 1.74, we exclude mantle serpentinisation in the basin centre. Overall, the off-shore region north of Corsica is faster than the south-western basin at shallow depths (<12 km) and slower at greater depth. This is linked to the varying sediment cover and the crustal thickness. In the south-western part, the opening of the basin is more developed, but we do not observe oceanic crust in our study area. The change between these domains appears gradual.
Data from the temporary AlpArray land stations as well as the OBSs are available through the AlpArray seismic network (Z3, https://doi.org/10.12686/alparray/z3_2015, AlpArray Seismic Network, 2015). Additional permanent land stations are available through the RESIF-RLBP French broadband network (https://doi.org/10.15778/RESIF.FR, RESIF, 1995), the Regional Seismic Network of North-Western Italy (https://doi.org/10.7914/SN/GU, University of Genoa, 1967), the Mediterranean Very Broadband Seismographic Network (MedNet, https://doi.org/10.13127/SD/FBBBTDTD6Q, MedNet Project Partner Institutions, 1990), and the Italian National Seismic Network (https://doi.org/10.13127/SD/X0FXnH7QfY, INGV Seismological Data Centre, 2006). Data from temporary AlpArray stations (including the OBSs) are accessible to AlpArray members. They will be freely accessible after March 2022.
The Supplement contains the lag-time input for the 2D group velocity tomography, the resulting grids, the codes we use to create the 1D inversion input (from group velocity maps), and the shear-wave velocity *.xyz files resulting from the 1D shear-wave inversion.
The supplement related to this article is available online at: https://doi.org/10.5194/se-12-2597-2021-supplement.
Complete list of the AlpArray Working Group: György Hetényi, Rafael Abreu, Ivo Allegretti, Maria-Theresia Apoloner, Coralie Aubert, Simon Besançon, Maxime Bès de Berc, Götz Bokelmann, Didier Brunel, Marco Capello, Martina Čarman, Adriano Cavaliere, Jérôme Chèze, Claudio Chiarabba, John Clinton, Glenn Cougoulat, Wayne C. Crawford, Luigia Cristiano, Tibor Czifra, Ezio D'Alema, Stefania Danesi, Romuald Daniel, Anke Dannowski, Iva Dasović, Anne Deschamps, Jean-Xavier Dessa, Cécile Doubre, Sven Egdorf, ETHZ-SED Electronics Lab, Tomislav Fiket, Kasper Fischer, Wolfgang Friederich, Florian Fuchs, Sigward Funke, Domenico Giardini, Aladino Govoni, Zoltán Gráczer, Gidera Gröschl, Stefan Heimers, Ben Heit, Davorka Herak, Marijan Herak, Johann Huber, Dejan Jarić, Petr Jedlička, Yan Jia, Hélène Jund, Edi Kissling, Stefan Klingen, Bernhard Klotz, Petr Kolínský, Heidrun Kopp, Michael Korn, Josef Kotek, Lothar Kühne, Krešo Kuk, Dietrich Lange, Jürgen Loos, Sara Lovati, Deny Malengros, Lucia Margheriti, Christophe Maron, Xavier Martin, Marco Massa, Francesco Mazzarini, Thomas Meier, Laurent Métral, Irene Molinari, Milena Moretti, Anna Nardi, Jurij Pahor, Anne Paul, Catherine Péquegnat, Daniel Petersen, Damiano Pesaresi, Davide Piccinini, Claudia Piromallo, Thomas Plenefisch, Jaroslava Plomerová, Silvia Pondrelli, Snježan Prevolnik, Roman Racine, Marc Régnier, Miriam Reiss, Joachim Ritter, Georg Rümpker, Simone Salimbeni, Marco Santulin, Werner Scherer, Sven Schippkus, Detlef Schulte-Kortnack, Vesna Šipka, Stefano Solarino, Daniele Spallarossa, Kathrin Spieker, Josip Stipčević, Angelo Strollo, Bálint Süle, Gyöngyvér Szanyi, Eszter Szűcs, Christine Thomas, Martin Thorwart, Frederik Tilmann, Stefan Ueding, Massimiliano Vallocchia, Luděk Vecsey, René Voigt, Joachim Wassermann, Zoltán Wéber, Christian Weidle, Viktor Wesztergom, Gauthier Weyland, Stefan Wiemer, Felix Noah Wolf, David Wolyniec, Thomas Zieke, Mladen Živčić, Helena Žlebčíková. For further information regarding the team, please visit http://www.alparray.ethz.ch/ (last access: 26 April 2021).
DL, MT, IG, and HK were responsible for the conception of this study. DL, WC, AD, and HK were responsible for designing the OBS network. FNW, DL, AD, MT, WC, and HK acquired OBS data. FNW analysed the data with the support of all co-authors. WC and LW provided software and expertise for the compliance correction and multiple filter technique. All authors interpreted the data. FNW prepared the paper, and all authors critically reviewed it.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article is part of the special issue “New insights into the tectonic evolution of the Alps and the adjacent orogens”. It is not associated with a conference.
This contribution is part of the German priority program “Mountain Building Processes in Four Dimensions (MB-4D)” SPP 2017 (http://www.spp-mountainbuilding.de/index.html, last access: 26 April 2021) and of the international research initiative AlpArray. We thank the captains and crews of RV Pourquoi Pas? and RV Maria S. Merian for their effort during deployment and recovery of the OBS network. We also thank the participating scientific crews. We thank Anne Paul for attending the cruise and her support with the OBSs from the Institut de physique du globe de Paris. The deployment of the French component of the AlpArray Seismic Network was funded by the AlpArray-FR project of the Agence Nationale de la Recherche (contract ANR-15-CE31-0015). For calculations, we use the Python-based tool ObsPy (Beyreuther et al., 2010; Krischer et al., 2015) and MATLAB (https://de.mathworks.com/, last access: 26 April 2021). Figures were created using Generic Mapping Tools version 6 (Wessel et al., 2019), MATLAB, and Inkscape (https://www.inkscape.org, last access: 26 April 2021). We thank Thomas Meier and his working group at Kiel University for their help with the multiple filter technique calculations and feedback on the dispersion curves. Further, we thank the AlpArray Seismic Network Team; a complete list of members can be found here: http://www.alparray.ethz.ch/en/seismic_network/backbone/data-policy-and-citation/ (last access: 26 April 2021). We thank two anonymous reviewers and editor Mark Handy for their fruitful comments. Also, we thank Louisa Murray-Bergquist for checking the language of our manuscript.
This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. LA 2970/3-1).
The article processing charges for this open-access publication were covered by the GEOMAR Helmholtz Centre for Ocean Research Kiel.
This paper was edited by Mark R. Handy and reviewed by two anonymous referees.
Adimah, N. I. and Padhy, S.: Ambient noise Rayleigh wave tomography across the Madagascar island, Geophys. J. Int., 220, 1657–1676, https://doi.org/10.1093/gji/ggz542, 2020.
Aki, K.: Space and Time Spectra of Stationary Stochastic Waves, with Special Reference to Microtremors, Bull. Earthq. Res. Inst. Univ. Tokyo, 415–456, 1957.
AlpArray Seismic Network: AlpArray Seismic Network (AASN) temporary component, AlpArray Working Group, Other/Seismic Network [data set], https://doi.org/10.12686/alparray/z3_2015, 2015.
André, G., Garreau, P., Garnier, V., and Fraunié, P.: Modelled variability of the sea surface circulation in the North-western Mediterranean Sea and in the Gulf of Lions, Ocean Dynam., 55, 294–308, https://doi.org/10.1007/s10236-005-0013-6, 2005.
Barmin, M. P., Ritzwoller, M. H., and Levshin, A. L.: A Fast and Reliable Method for Surface Wave Tomography, in: Monitoring the Comprehensive Nuclear-Test-Ban Treaty: Surface Waves, Birkhäuser Basel, Basel, 1351–1375, https://doi.org/10.1007/978-3-0348-8264-4_3, 2001.
Bell, S. W., Forsyth, D. W., and Ruan, Y.: Removing Noise from the Vertical Component Records of Ocean-Bottom Seismometers: Results from Year One of the Cascadia Initiative, B. Seismol. Soc. Am., 105, 300–313, https://doi.org/10.1785/0120140054, 2015.
Bensen, G. D., Ritzwoller, M. H., Barmin, M. P., Levshin, A. L., Lin, F., Moschetti, M. P., Shapiro, N. M., and Yang, Y.: Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements, Geophys. J. Int., 169, 1239–1260, https://doi.org/10.1111/j.1365-246X.2007.03374.x, 2007.
Beyreuther, M., Barsch, R., Krischer, L., Megies, T., Behr, Y., and Wassermann, J.: ObsPy: A Python Toolbox for Seismology, Seismol. Res. Lett., 81, 530–533, https://doi.org/10.1785/gssrl.81.3.530, 2010.
Boschi, L., Weemstra, C., Verbeke, J., Ekström, G., Zunino, A., and Giardini, D.: On measuring surface wave phase velocity from station–station cross-correlation of ambient signal, Geophys. J. Int., 192, 346–358, https://doi.org/10.1093/gji/ggs023, 2013.
Campillo, M. and Paul, A.: Long-Range Correlations in the Diffuse Seismic Coda, Science, 299, 547–549, https://doi.org/10.1126/science.1078551, 2003.
Campillo, M. and Roux, P.: Seismic Imaging and Monitoring with Ambient Noise Correlations, in: Treatise on Geophysics, Elservier, Amsterdam, 2014.
Carlson, R. L. and Miller, D. J.: Mantle wedge water contents estimated from seismic velocities in partially serpentinized peridotites, Geophys. Res. Lett., 30, 12–15, https://doi.org/10.1029/2002gl016600, 2003.
Contrucci, I., Nercessian, A., Béthoux, N., Mauffret, A., and Pascal, G.: A Ligurian (Western Mediterranean Sea) geophysical transect revisited, Geophys. J. Int., 146, 74–97, https://doi.org/10.1046/j.0956-540x.2001.01418.x, 2001.
Cox, C., Deaton, T., and Webb, S.: A Deep-Sea Differential Pressure Gauge, J. Atmospheric Ocean. Technol., 1, 237–246, https://doi.org/10.1175/1520-0426(1984)001<0237:ADSDPG>2.0.CO;2, 1984.
Crawford, W. C. and Webb, S. C.: Identifying and Removing Tilt Noise from Low-Frequency (<0.1 Hz) Seafloor Vertical Seismic Data, B. Seismol. Soc. Am., 90, 952–963, https://doi.org/10.1785/0119990121, 2000.
Crawford, W. C., Webb, S. C., and Hildebrand, J. A.: Estimating shear velocities in the oceanic crust from compliance measurements by two-dimensional finite difference modeling, J. Geophys. Res.-Sol. Ea., 103, 9895–9916, https://doi.org/10.1029/97JB03532, 1998.
Dannowski, A., Kopp, H., Grevemeyer, I., Lange, D., Thorwart, M., Bialas, J., and Wollatz-Vogt, M.: Seismic evidence for failed rifting in the Ligurian Basin, Western Alpine domain, Solid Earth, 11, 873–887, https://doi.org/10.5194/se-11-873-2020, 2020.
Dewangan, P., Reddy, R., Kamesh Raju, K. A., Singha, P., Aswini, K. K., Yatheesh, V., Samudrala, K., and Shuhail, M.: Nature of the Ambient Noise, Site Response, and Orientation of Ocean-Bottom Seismometers (OBSs): Scientific Results of a Passive Seismic Experiment in the Andaman Sea, B. Seismol. Soc. Am., 108, 248–259, https://doi.org/10.1785/0120170163, 2018.
Dziewonski, A. M. and Anderson, D. L.: Preliminary reference Earth model, Phys. Earth Planet. In., 25, 297–356, https://doi.org/10.1016/0031-9201(81)90046-7, 1981.
Dziewonski, A. M., Bloch, S., and Landisman, M.: A Technique for the Analysis of Transient Seismic Signals, B. Seismol. Soc. Am., 59, 427–444, 1969.
Faccenna, C., Funiciello, F., Giardini, D., and Lucente, P.: Episodic back-arc extension during restricted mantle convection in the Central Mediterranean, Earth Planet. Sc. Lett., 187, 105–116, https://doi.org/10.1016/S0012-821X(01)00280-1, 2001.
Finetti, I. R., Boccaletti, M., Bonini, M., Del Ben, A., Pipan, M., Prizzon, A., and Sani, F.: Lithospheric tectono-stratigraphic setting of the Ligurian Sea–Northern Apennines–Adriatic Foreland from integrated CROP seismic data, in: Deep Seismic Exploration of the Central Mediterranean and Italy, CROP PROJECT, Vol. 8, Trieste, Italy, 119–158, 2005.
Foti, S., Hollender, F., Garofalo, F., Albarello, D., Asten, M., Bard, P. Y., Comina, C., Cornou, C., Cox, B., Di Giulio, G., Forbriger, T., Hayashi, K., Lunedei, E., Martin, A., Mercerat, D., Ohrnberger, M., Poggi, V., Renalier, F., Sicilia, D., and Socco, V.: Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project, B. Earthq. Eng., 16, 2367–2420, https://doi.org/10.1007/s10518-017-0206-7, 2018.
Gailler, A., Klingelhoefer, F., Olivet, J.-L. L., and Aslanian, D.: Crustal structure of a young margin pair: New results across the Liguro–Provencal Basin from wide-angle seismic tomography, Earth Planet. Sc. Lett., 286, 333–345, https://doi.org/10.1016/j.epsl.2009.07.001, 2009.
Gattacceca, J., Deino, A., Rizzo, R., Jones, D. S., Henry, B., Beaudoin, B., and Vadeboin, F.: Miocene rotation of Sardinia: New paleomagnetic and geochronological constraints and geodynamic implications, Earth Planet. Sc. Lett., 258, 359–377, https://doi.org/10.1016/j.epsl.2007.02.003, 2007.
Goutorbe, B., de Oliveira Coelho, D. L., and Drouet, S.: Rayleigh wave group velocities at periods of 6–23 s across Brazil from ambient noise tomography, Geophys. J. Int., 203, 869–882, https://doi.org/10.1093/gji/ggv343, 2015.
Gueguen, E., Doglioni, C., and Fernandez, M.: On the post-25 Ma geodynamic evolution of the western Mediterranean, Tectonophysics, 298, 259–269, https://doi.org/10.1016/S0040-1951(98)00189-9, 1998.
Guerin, G., Rivet, D., Deschamps, A., Larroque, C., Mordret, A., Dessa, J.-X., and Martin, X.: High resolution ambient noise tomography of the Southwestern Alps and the Ligurian margin, Geophys. J. Int., 220, 806–820, https://doi.org/10.1093/gji/ggz477, 2020.
Hansen, P. C.: Analysis of Discrete Ill-Posed Problems by Means of the L-Curve, SIAM Rev., 34, 561–580, https://doi.org/10.1137/1034115, 1992.
Harmon, N., Forsyth, D., and Webb, S.: Using Ambient Seismic Noise to Determine Short-Period Phase Velocities and Shallow Shear Velocities in Young Oceanic Lithosphere, B. Seismol. Soc. Am., 97, 2009–2023, https://doi.org/10.1785/0120070050, 2007.
Harmon, N., Henstock, T., Tilmann, F., Rietbrock, A., and Barton, P.: Shear velocity structure across the Sumatran Forearc-Arc, Geophys. J. Int., 189, 1306–1314, https://doi.org/10.1111/j.1365-246X.2012.05446.x, 2012.
Herrmann, R. B.: Computer Programs in Seismology: An Evolving Tool for Instruction and Research, Seismol. Res. Lett., 84, 1081–1088, https://doi.org/10.1785/0220110096, 2013.
Hetényi, G., Molinari, I., Clinton, J., Bokelmann, G., Bondár, I., Crawford, W. C., Dessa, J.-X., Doubre, C., Friederich, W., Fuchs, F., Giardini, D., Gráczer, Z., Handy, M. R., Herak, M., Jia, Y., Kissling, E. E., Kopp, H., Korn, M., Margheriti, L., Meier, T., Mucciarelli, M., Paul, A., Pesaresi, D., Piromallo, C., Plenefisch, T., Plomerová, J., Ritter, J., Rümpker, G., Šipka, V., Spallarossa, D., Thomas, C., Tilmann, F., Wassermann, J., Weber, M., Wéber, Z., Wesztergom, V., and Živčić, M.: The AlpArray Seismic Network: A Large-Scale European Experiment to Image the Alpine Orogen, Surv. Geophys., 39, 1009–1033, https://doi.org/10.1007/s10712-018-9472-4, 2018.
INGV Seismological Data Centre: Rete Sismica Nazionale (RSN), Istituto Nazionale di Geofisica e Vulcanologia (INGV), Italy [data set], https://doi.org/10.13127/SD/X0FXNH7QFY, 2006.
Jolivet, L. and Faccenna, C.: Mediterranean extension and the Africa-Eurasia collision, Tectonics, 19, 1095–1106, https://doi.org/10.1029/2000TC900018, 2000.
Jolivet, L., Gorini, C., Smit, J., and Leroy, S.: Continental breakup and the dynamics of rifting in back-arc basins: The Gulf of Lion margin, Tectonics, 34, 662–679, https://doi.org/10.1002/2014TC003570, 2015.
Kästle, E. D., El-Sharkawy, A., Boschi, L., Meier, T., Rosenberg, C., Bellahsen, N., Cristiano, L., and Weidle, C.: Surface Wave Tomography of the Alps Using Ambient-Noise and Earthquake Phase Velocity Measurements, J. Geophys. Res.-Sol. Ea., 123, 1770–1792, https://doi.org/10.1002/2017JB014698, 2018.
Kästle, E. D., Rosenberg, C., Boschi, L., Bellahsen, N., Meier, T., and El-Sharkawy, A.: Slab Break-offs in the Alpine Subduction Zone, Solid Earth Discuss. [preprint], https://doi.org/10.5194/se-2019-17, 2019.
Krischer, L., Megies, T., Barsch, R., Beyreuther, M., Lecocq, T., Caudron, C., and Wassermann, J.: ObsPy: a bridge for seismology into the scientific Python ecosystem, Comput. Sci. Discov., 8, 014003–014003, https://doi.org/10.1088/1749-4699/8/1/014003, 2015.
Laubscher, H., Biella, G. C., Cassinis, R., Gelati, R., Lozej, A., Scarascia, S., and Tabacco, I.: The collisional knot in Liguria, Geol. Rundsch., 81, 275–289, https://doi.org/10.1007/BF01828598, 1992.
Le Breton, E., Handy, M. R., Molli, G., and Ustaszewski, K.: Post-20 Ma Motion of the Adriatic Plate: New Constraints From Surrounding Orogens and Implications for Crust-Mantle Decoupling, Tectonics, 36, 3135–3154, https://doi.org/10.1002/2016TC004443, 2017.
Le Pichon, X., Rangin, C., Hamon, Y., Loget, N., Lin, J. Y., Andreani, L., and Flotte, N.: Geodynamics of the France Southeast Basin, B. Soc. Géol. Fr., 181, 477–501, 2010.
Lin, P.-Y. P., Gaherty, J. B., Jin, G., Collins, J. A., Lizarralde, D., Evans, R. L., and Hirth, G.: High-resolution seismic constraints on flow dynamics in the oceanic asthenosphere, Nature, 535, 538–541, https://doi.org/10.1038/nature18012, 2016.
Lobkis, O. I. and Weaver, R. L.: On the emergence of the Green's function in the correlations of a diffuse field, J. Acoust. Soc. Am., 110, 3011–3017, https://doi.org/10.1121/1.1417528, 2001.
Lu, Y., Stehly, L., and Paul, A.: High-resolution surface wave tomography of the European crust and uppermost mantle from ambient seismic noise, Geophys. J. Int., 214, 1136–1150, https://doi.org/10.1093/gji/ggy188, 2018.
Makris, J., Egloff, F., Nicolich, R., and Rihm, R.: Crustal structure from the Ligurian Sea to the Northern Apennines – a wide angle seismic transect, Tectonophysics, 301, 305–319, https://doi.org/10.1016/S0040-1951(98)00225-X, 1999.
McNamara, D. E. and Buland, R. P.: Ambient Noise Levels in the Continental United States, B. Seismol. Soc. Am., 94, 1517–1527, https://doi.org/10.1785/012003001, 2004.
MedNet Project Partner Institutions: Mediterranean Very Broadband Seismographic Network (MedNet), Istituto Nazionale di Geofisica e Vulcanologia (INGV) [data set], https://doi.org/10.13127/SD/FBBBTDTD6Q, 1990.
Meier, T., Dietrich, K., Stöckhert, B., and Harjes, H. P.: One-dimensional models of shear wave velocity for the eastern Mediterranean obtained from the inversion of Rayleigh wave phase velocities and tectonic implications, Geophys. J. Int., 156, 45–58, https://doi.org/10.1111/j.1365-246X.2004.02121.x, 2004.
Millot, C. and Wald, L.: The effects of Mistral wind on the Ligurian current near Provence, Oceanol. Acta, 3, 399–402, 1980.
Molinari, I., Argnani, A., Morelli, A., and Basini, P.: Development and Testing of a 3D Seismic Velocity Model of the Po Plain Sedimentary Basin, Italy, B. Seismol. Soc. Am., 105, 753–764, https://doi.org/10.1785/0120140204, 2015a.
Molinari, I., Verbeke, J., Boschi, L., Kissling, E., and Morelli, A.: Italian and Alpine three-dimensional crustal structure imaged by ambient-noise surface-wave dispersion, Geochem. Geophy. Geosy., 16, 4405–4421, https://doi.org/10.1002/2015GC006176, 2015b.
Moulin, M., Klingelhoefer, F., Afilhado, A., Aslanian, D., Schnurle, P., Nouzé, H., Rabineau, M., Beslier, M.-O., and Feld, A.: Deep crustal structure across a young passive margin from wide-angle and reflection seismic data (The SARDINIA Experiment) – I. Gulf of Lion's margin, Bull. Société Géologique Fr., 186, 309–330, https://doi.org/10.2113/gssgfbull.186.4-5.309, 2015.
Pasi, F., Orlandi, A., Onorato, L. F., and Gallino, S.: A study of the 1 and 2 January 2010 sea-storm in the Ligurian Sea, Adv. Sci. Res., 109–115, https://doi.org/10.5194/asr-6-109-2011, 2011.
Prieto, G. A., Lawrence, J. F., and Beroza, G. C.: Anelastic Earth structure from the coherency of the ambient seismic field, J. Geophys. Res., 114, B07303, https://doi.org/10.1029/2008JB006067, 2009.
Rawlinson, N. and Sambridge, M.: Wave front evolution in strongly heterogeneous layered media using the fast marching method, Geophys. J. Int., 156, 631–647, https://doi.org/10.1111/j.1365-246X.2004.02153.x, 2004.
Rawlinson, N. and Sambridge, M.: The fast marching method: an effective tool for tomographic imaging and tracking multiple phases in complex layered media, Explor. Geophys., 36, 341–350, https://doi.org/10.1071/EG05341, 2005.
RESIF: RESIF-RLBP French Broad-band network, RESIF-RAP strong motion network and other seismic stations in metropolitan France, RESIF – Réseau Sismologique et géodésique Français [data set], https://doi.org/10.15778/resif.fr, 1995.
Rollet, N., Déverchère, J., Beslier, M.-O., Guennoc, P., Réhault, J.-P., Sosson, M., and Truffert, C.: Back arc extension, tectonic inheritance, and volcanism in the Ligurian Sea, Western Mediterranean, Tectonics, 21, 6-1–6-23, https://doi.org/10.1029/2001TC900027, 2002.
Rosenbaum, G., Lister, G. S., and Duboz, C.: Reconstruction of the tectonic evolution of the western Mediterranean since the Oligocene, J. Virtual Explor., 08, 107–126, https://doi.org/10.3809/jvirtex.2002.00053, 2002.
Ryan, W. B. F., Carbotte, S. M., Coplan, J. O., O'Hara, S., Melkonian, A., Arko, R., Weissel, R. A., Ferrini, V., Goodwillie, A., Nitsche, F., Bonczkowski, J., and Zemsky, R.: Global Multi-Resolution Topography synthesis, Geochem. Geophy. Geosy., 10, Q03014, https://doi.org/10.1029/2008GC002332, 2009.
Sandwell, D. T., Müller, R. D., Smith, W. H. F., Garcia, E., and Francis, R.: New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure, Science, 346, 65–67, https://doi.org/10.1126/science.1258213, 2014.
Schettino, A. and Turco, E.: Plate kinematics of the Western Mediterranean region during the Oligocene and Early Miocene, Geophys. J. Int., 166, 1398–1423, https://doi.org/10.1111/j.1365-246X.2006.02997.x, 2006.
Schmid, S. M., Fügenschuh, B., Kissling, E., and Schuster, R.: Tectonic map and overall architecture of the Alpine orogen, Eclogae Geol. Helv., 97, 93–117, https://doi.org/10.1007/s00015-004-1113-x, 2004.
Schmidt-Ausch, M. C. and Haberland, C.: DEPAS (Deutscher Geräte-Pool für amphibische Seismologie): German Instrument Pool for Amphibian Seismology, J. Large-Scale Res. Facil. JLSRF, 3, A122, https://doi.org/10.17815/jlsrf-3-165, 2017.
Sethian, J. A.: A fast marching level set method for monotonically advancing fronts, P. Natl. Acad. Sci., 93, 1591–1595, https://doi.org/10.1073/pnas.93.4.1591, 1996.
Sethian, J. A. and Popovici, A. M.: 3-D traveltime computation using the fast marching method, Geophysics, 64, 516–523, https://doi.org/10.1190/1.1444558, 1999.
Shapiro, N. M., Campillo, M., Stehly, L., and Ritzwoller, M. H.: High-Resolution Surface-Wave Tomography from Ambient Seismic Noise, Science, 307, 1615–1618, https://doi.org/10.1126/science.1108339, 2005.
Shillington, D. J., Minshull, T. A., Peirce, C., and O'Sullivan, J. M.: P- and S-wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland, Geophys. Prospect., 56, 197–211, https://doi.org/10.1111/j.1365-2478.2007.00669.x, 2007.
Speranza, F., Villa, I. M., Sagnotti, L., Florindo, F., Cosentino, D., Cipollari, P., and Mattei, M.: Age of the Corsica–Sardinia rotation and Liguro–Provençal Basin spreading: new paleomagnetic and Ar/Ar evidence, Tectonophysics, 347, 231–251, https://doi.org/10.1016/S0040-1951(02)00031-8, 2002.
Stähler, S. C., Sigloch, K., Hosseini, K., Crawford, W. C., Barruol, G., Schmidt-Aursch, M. C., Tsekhmistrenko, M., Scholz, J.-R., Mazzullo, A., and Deen, M.: Performance report of the RHUM-RUM ocean bottom seismometer network around La Réunion, western Indian Ocean, Adv. Geosci., 41, 43–63, https://doi.org/10.5194/adgeo-41-43-2016, 2016.
Takeo, A., Forsyth, D. W., Weeraratne, D. S., and Nishida, K.: Estimation of azimuthal anisotropy in the NW Pacific from seismic ambient noise in seafloor records, Geophys. J. Int., 199, 11–22, https://doi.org/10.1093/gji/ggu240, 2014.
Thorwart, M., Dannowski, A., Grevemeyer, I., Lange, D., Kopp, H., Petersen, F., Crawford, W., Paul, A., and the AlpArray Working Group: Basin inversion: Reactivated rift structures in the Ligurian Sea revealed by OBS, Solid Earth Discuss. [preprint], https://doi.org/10.5194/se-2021-9, in review, 2021.
Tonegawa, T., Yamashita, Y., Takahashi, T., Shinohara, M., Ishihara, Y., Kodaira, S., and Kaneda, Y.: Spatial relationship between shallow very low frequency earthquakes and the subducted Kyushu-Palau Ridge in the Hyuga-nada region of the Nankai subduction zone, Geophys. J. Int., 222, 1542–1554, https://doi.org/10.1093/gji/ggaa264, 2020.
University of Genoa: Regional Seismic Network of North Western Italy, International Federation of Digital Seismograph Networks [data set], https://doi.org/10.7914/SN/GU, 1967.
Vigliotti, L. and Langenheim, V. E.: When did Sardinia stop rotating? New palaeomagnetic results, Terra Nova, 7, 424–435, https://doi.org/10.1111/j.1365-3121.1995.tb00538.x, 1995.
Wapenaar, K., Draganov, D., Snieder, R., Campman, X., and Verdel, A.: Tutorial on seismic interferometry: Part 1 – Basic principles and applications, Geophysics, 75, 75A195-75A209, https://doi.org/10.1190/1.3457445, 2010a.
Wapenaar, K., Slob, E., Snieder, R., and Curtis, A.: Tutorial on seismic interferometry: Part 2 – Underlying theory and new advances, Geophysics, 75, 75A211–75A227, https://doi.org/10.1190/1.3463440, 2010b.
Webb, S. C.: Broadband seismology and noise under the ocean, Rev. Geophys., 36, 105–142, https://doi.org/10.1029/97RG02287, 1998.
Webb, S. C. and Crawford, W. C.: Long-period seafloor seismology and deformation under ocean waves, B. Seismol. Soc. Am., 89, 1535–1542, 1999.
Wessel, P., Luis, J. F., Uieda, L., Scharroo, R., Wobbe, F., Smith, W. H. F., and Tian, D.: The Generic Mapping Tools Version 6, Geochem. Geophy. Geosy., 20, 5556–5564, https://doi.org/10.1029/2019GC008515, 2019.