SUMMARY The Groningen gas field is one of the largest gas fields in Europe. The continuous gas extraction led to an induced seismic activity in the area. In order to monitor the seismic activity and study the gas field many permanent and temporary seismic arrays were deployed. In particular, the extraction of the shear wave velocity model is crucial in seismic hazard assessment. Local S-wave velocity-depth profiles allow us the estimation of a potential amplification due to soft sediments. Ambient seismic noise tomography is an interesting alternative to traditional methods that were used in modelling the S-wave velocity. The ambient noise field consists mostly of surface waves, which are sensitive to the Swave and if inverted, they reveal the corresponding S-wave structures. In this study, we present results of a depth inversion of surface waves obtained from the cross-correlation of 1 month of ambient noise data from four flexible networks located in the Groningen area. Each block consisted of about 400 3-C stations. We compute group velocity maps of Rayleigh and Love waves using a straight-ray surface wave tomography. We also extract clear higher modes of Love and Rayleigh waves. The S-wave velocity model is obtained with a joint inversion of Love and Rayleigh waves using the Neighbourhood Algorithm. In order to improve the depth inversion, we use the mean phase velocity curves and the higher modes of Rayleigh and Love waves. Moreover, we use the depth of the base of the North Sea formation as a hard constraint. This information provides an additional constraint for depth inversion, which reduces the S-wave velocity uncertainties. The final S-wave velocity models reflect the geological structures up to 1 km depth and in perspective can be used in seismic risk modelling.

No abstract

SUMMARY Unveiling the mechanisms of earthquake and volcanic eruption preparation requires improving our ability to monitor the rock mass response to transient stress perturbations at depth. The standard passive monitoring seismic interferometry technique based on coda waves is robust but recovering accurate and properly localized P- and S-wave velocity temporal anomalies at depth is intrinsically limited by the complexity of scattered, diffracted waves. In order to mitigate this limitation, we propose a complementary, novel, passive seismic monitoring approach based on detecting weak temporal changes of velocities of ballistic waves recovered from seismic noise correlations. This new technique requires dense arrays of seismic sensors in order to circumvent the bias linked to the intrinsic high sensitivity of ballistic waves recovered from noise correlations to changes in the noise source properties. In this work we use a dense network of 417 seismometers in the Groningen area of the Netherlands, one of Europe's largest gas fields. Over the course of 1 month our results show a 1.5 per cent apparent velocity increase of the P wave refracted at the basement of the 700-m-thick sedimentary cover. We interpret this unexpected high value of velocity increase for the refracted wave as being induced by a loading effect associated with rainfall activity and possibly canal drainage at surface. We also observe a 0.25 per cent velocity decrease for the direct P-wave travelling in the near-surface sediments and conclude that it might be partially biased by changes in time in the noise source properties even though it appears to be consistent with complementary results based on ballistic surface waves presented in a companion paper and interpreted as a pore pressure diffusion effect following a strong rainfall episode. The perspective of applying this new technique to detect continuous localized variations of seismic velocity perturbations at a few kilometres depth paves the way for improved in situ earthquake, volcano and producing reservoir monitoring.