Due to the deep socioeconomic implications, induced seismicity is a timely and increasingly relevant topic of interest for the general public. Cases of induced seismicity have a global distribution and involve a large number of industrial operations, with many documented cases from as far back to the beginning of the twentieth century. However, the sparse and fragmented documentation available makes it difficult to have a clear picture on our understanding of the physical phenomenon and consequently in our ability to mitigate the risk associated with induced seismicity. This review presents a unified and concise summary of the still open questions related to monitoring, discrimination, and management of induced seismicity in the European context and, when possible, provides potential answers. We further discuss selected critical European cases of induced seismicity, which led to the suspension or reduction of the related industrial activities.

Few problems in elastodynamics have a closed-form analytical solution. The others can be investigated with semianalytical methods, but in general one is not sure whether these methods give reliable solutions. The same happens with numerical techniques: for instance, finite difference methods solve, in principle, any complex problem, including those with arbitrary inhomogeneities and boundary conditions. However, there is no way to verify the quantitative correctness of the solutions. The major problems are stability with respect to material properties, numerical dispersion, and the treatment of boundary conditions. In practice, these problems may produce inaccurate solutions. In this paper, the study of complex problems with two different numerical grid techniques in order to cross-check the solutions is proposed. Interface waves, in particular, are emphasized, since they pose the major difficulties due to the need to implement boundary conditions. The first method is based on global differential operators where the solution is expanded in terms of the Fourier basis and Chebyshev polynomials, while the second is the spectral element method, an extension of the finite element method that uses Chebyshev polynomials as interpolating functions. Both methods have spectral accuracy up to approximately the Nyquist wave number of the grid. Moreover, both methods implement the boundary conditions in a natural way, particularly the spectral element algorithm. We first solve Lamb’s problem and compare numerical and analytical solutions; then, the problem of dispersed Rayleigh waves, and finally, the two-quarter space problem. We show that the modeling algorithms correctly reproduce the analytical solutions and yield a perfect matching when these solutions do not exist. The combined modeling techniques provide a powerful tool for solving complex problems in elastodynamics.

International audiencen a low‐seismicity context, the use of numerical simulations becomes essential due to the lack of representative earthquakes for empirical approaches. The goals of the EUROSEISTEST Verification and Validation Project (E2VP) are to provide (1) a quantitative analysis of accuracy of the current, most advanced numerical methods applied to realistic 3D models of sedimentary basins (verification) and (2) a quantitative comparison of the recorded ground motions with their numerical predictions (validation). The target is the EUROSEISTEST site located within the Mygdonian basin, Greece. The site is instrumented with surface and borehole accelerometers, and a 3D model of the medium is available. The simulations are performed up to 4 Hz, beyond the 0.7 Hz fundamental frequency, thus covering a frequency range at which ground motion undergoes significant amplification. The discrete representation of material heterogeneities, the attenuation model, the approximation of the free surface, and nonreflecting boundaries are identified as the main sources of differences among the numerical predictions. The predictions well reproduce some, but not all, features of the actual site effect. The differences between real and predicted ground motions have multiple origins: the accuracy of source parameters (location, hypocentral depth, and focal mechanism), the uncertainties in the description of the geological medium (damping, internal sediment layering structure, and shape of the sediment‐basement interface). Overall, the agreement reached among synthetics up to 4 Hz despite the complexity of the basin model, with code‐to‐code differences much smaller than predictions‐to‐observations differences, makes it possible to include the numerical simulations in site‐specific analysis in the 3D linear case and low‐to‐intermediate frequency range

Differences between 3-D numerical predictions of earthquake ground motion in the Mygdonian basin near Thessaloniki, Greece, led us to define four canonical stringent models derived from the complex realistic 3-D model of the Mygdonian basin. Sediments atop an elastic bedrock are modelled in the 1D-sharp and 1D-smooth models using three homogeneous layers and smooth velocity distribution, respectively. The 2D-sharp and 2D-smooth models are extensions of the 1-D models to an asymmetric sedimentary valley. In all cases, 3-D wavefields include strongly dispersive surface waves in the sediments. We compared simulations by the Fourier pseudo-spectral method (FPSM), the Legendre spectral-element method (SEM) and two formulations of the finite-difference method (FDM-S and FDM-C) up to 4 Hz. The accuracy of individual solutions and level of agreement between solutions vary with type of seismic waves and depend on the smoothness of the velocity model. The level of accuracy is high for the body waves in all solutions. However, it strongly depends on the discrete representation of the material interfaces (at which material parameters change discontinuously) for the surface waves in the sharp models. An improper discrete representation of the interfaces can cause inaccurate numerical modelling of surface waves. For all the numerical methods considered, except SEM with mesh of elements following the interfaces, a proper implementation of interfaces requires definition of an effective medium consistent with the interface boundary conditions. An orthorhombic effective medium is shown to significantly improve accuracy and preserve the computational efficiency of modelling. The conclusions drawn from the analysis of the results of the canonical cases greatly help to explain differences between numerical predictions of ground motion in realistic models of the Mygdonian basin. We recommend that any numerical method and code that is intended for numerical prediction of earthquake ground motion should be verified through stringent models that would make it possible to test the most important aspects of accuracy.

Accurate and automated locations of microseismic events are desirable for many seismological and industrial applications. The analysis of microseismicity is particularly challenging because of weak seismic signals with low signal-to-noise ratio. Traditional location approaches rely on automated picking, based on individual seismograms, and make no use of the coherency information between signals at different stations. This strong limitation has been overcome by full-waveform location methods, which exploit the coherency of waveforms at different stations and improve the location robustness even in presence of noise. However, the performance of these methods strongly depend on the accuracy of the adopted velocity model, which is often quite rough; inaccurate models result in large location errors. We present an improved waveform stacking location method based on source-specific station corrections. Our method inherits the advantages of full-waveform location methods while strongly mitigating the dependency on the accuracy of the velocity model. With this approach the influence of an inaccurate velocity model on the results is restricted to the estimation of travel times solely within the seismogenic volume, but not for the entire source-receiver path. We finally successfully applied our new method to a realistic synthetic dataset as well as real data.