On the wobbles of phase-velocity dispersion curves

To calculate phase-velocity dispersion curves,we introduce amethodwhich reflects both structural and dynamic effects of wave propagation and interference. Rayleigh-wave fundamentalmode surface waves from the South Atlantic Ocean earthquake of 19 August 2016, M = 7.4, observed at the AlpArray network in Europe are strongly influenced by the upper-mantle lowvelocity zone under the Cameroon Volcanic Line in Central Africa. Predicting phase-delay times affected by diffraction from this heterogeneity for each station gives phase velocities as they would be determined using the classical two-station method as well as the advanced array-beamforming method. Synthetics from these two methods are thus compared with measurements. We show how the dynamic phase velocity differs from the structural phase velocity, howthese differences evolve in space and howtwo-station and arraymeasurements are affected. In principle, arrays are affected with the same uncertainty as the two-station measurements. The dynamic effects can be several times larger than the error caused by the unknown arrival angle in case of the two-station method. The non-planarity of the waves and its relation to the arrival angle and dynamic phase-velocity deviations is discussed. Our study is complemented by extensive review of literature related to the surface wave phase-velocity measurement of the last 120 years

Azimuthal anisotropy in the wider Vienna basin region: a proxy for the present-day stress field and deformation

International audienc

Coda-Q in the 2.5-20 s period band from seismic noise : application to the greater Alpine area

Coda-Q is used to estimate the attenuation and scattering properties of the Earth. So far focus has been on earthquake data at frequencies above 1 Hz, as the high noise level in the first and second microseismic peak, and possibly lower scattering coefficient, hinder stable measurements at lower frequencies. In this work, we measure and map coda-Q in the period bands 2.5-5 s, 5-10 s and 10-20 s in the greater Alpine region using noise cross-correlations between station pairs, based on data from permanent seismic stations and from the temporary AlpArray experiment. The observed coda-Q for short interstation distances is independent of azimuth so there is no indication of influence of the directivity of the incoming noise field on our measurements. In the 2.5-5 s and 5-10 s period bands, our measurements are self-consistent, and we observe stable geographic patterns of low and high coda-Q in the period bands 2.5-5 s and 5-10 s. In the period band 10-20 s, the dispersion of our measurements increases and geographic patterns become speculative. The coda-Q maps show that major features are observed with high resolution, with a very good geographical resolution of for example low coda-Q in the Po Plain. There is a sharp contrast between the Po Plain and the Alps and Apennines where coda-Q is high, with the exception a small area in the Swiss Alps which may be contaminated by the low coda-Q of the Po Plain. The coda of the correlations is too short to make independent measurements at different times within the coda, so we cannot distinguish between intrinsic and scattering Q. Measurements on more severely selected data sets and longer time-series result in identical geographical patterns but lower numerical values. Therefore, high coda-Q values may be overestimated, but the geographic distribution between high and low coda-Q areas is respected. Our results demonstrate that noise correlations are a promising tool for extending coda-Q measurements to frequencies lower than those analysed with earthquake data

Imaging seismic wave-fields with AlpArray and neighboring European networks

The AlpArray experiment and the deployment of Swath-D together with the dense permanent network in Italy allow for detailed imaging of the spatio-temporal imaging complexity of seismic wave-fields within the greater Alpine region. The distance of any point within the area to the nearest station is less than 30 km, resulting in an average inter-station distance of about 45 km. With a much denser deployment in a smaller region of the Alps (320 km in length and 140 km wide), the Swath-D network possesses an average inter-station distance of about 15 km. We show that seismogram sections with a spatial sampling of less than 5 km can be obtained using recordings of these regional arrays for just a single event. Multiply reflected body waves can be observed for up to 2 h after source time. In addition, we provide and describe animations of long-period seismic wave-fields using recordings of about 1300-1600 broadband stations for six representative earthquakes. These illustrate the considerable spatio-temporal variability of the wave-field's properties at a high lateral resolution. Within denser station distributions like those provided by Swath-D, even shorter period body and surface wave features can be recovered. The decrease of the horizontal wavelength from P to S to surface waves, deviations from spherically symmetric wavefronts, and the capability to detect multi-orbit arrivals are demonstrated qualitatively by the presented wave-field animations, which are a valuable tool for educational, quality control, and research purposes. We note that the information content of the acquired datasets can only be adequately explored by application of appropriate quantitative methods accounting for the considerable complexity of the seismic wave-fields as revealed by the now available station configuration

Evidence for radial anisotropy in the lower crust of the Apennines from Bayesian ambient noise tomography in Europe

Probing seismic anisotropy of the lithosphere provides valuable clues on the fabric of rocks. We present a 3-D probabilistic model of shear wave velocity and radial anisotropy of the crust and uppermost mantle of Europe, focusing on the mountain belts of the Alps and Apennines. The model is built from Love and Rayleigh dispersion curves in the period range 5–149 s. Data are extracted from seismic ambient noise recorded at 1521 broad-band stations, including the AlpArray network. The dispersion curves are first combined in a linearized least squares inversion to obtain 2-D maps of group velocity at each period. Love and Rayleigh maps are then jointly inverted at depth for shear wave velocity and radial anisotropy using a Bayesian Monte Carlo scheme that accounts for the trade-off between radial anisotropy and horizontal layering. The isotropic part of our model is consistent with previous studies. However, our anisotropy maps differ from previous large scale studies that suggested the presence of significant radial anisotropy everywhere in the European crust and shallow upper mantle. We observe instead that radial anisotropy is mostly localized beneath the Apennines while most of the remaining European crust and shallow upper mantle is isotropic. We attribute this difference to trade-offs between radial anisotropy and thin (hectometric) layering in previous studies based on least-squares inversions and long period data (>30 s). In contrast, our approach involves a massive data set of short period measurements and a Bayesian inversion that accounts for thin layering. The positive radial anisotropy (VSH > VSV) observed in the lower crust of the Apennines cannot result from thin layering. We rather attribute it to ductile horizontal flow in response to the recent and present-day extension in the region