We examine seismic velocity changes due to the Earth tide by conducting cross-correlation function (CCF) analyses of ambient seismic noise recorded at a small array composed of seven seismometers in northeastern Japan. We calculate CCFs for the dilatational and contractional episodes that are predicted from theoretical tidal volumetric strains. CCFs of the two episodes are highly correlated, but tiny differences are found in their phases. The phase differences are explained by seismic velocity changes of À0.19 ± 0.06% at 1-2 Hz, which are interpreted to be caused by opening/closure of cracks or pores in the shallow subsurface due to the tidal strain. Strain sensitivities of the seismic velocity changes are estimated to be 6.9 × 10 4 strain À1 , which are almost consistent with those reported in previous studies using artificial sources.

We investigate the stress sensitivity of velocity changes at Izu‐Oshima in Japan using seismic interferometry method. We calculate cross correlation functions (CCFs) of ambient noises recorded by four seismic stations on the active volcano from 1 January 2012 to 31 December 2015 at 0.5–1 Hz, 1–2 Hz, and 2–4 Hz. Applying moving time windows to calculate cross spectrum between daily CCFs and reference CCF, which is the average for the all observation period, we compute daily velocity changes. The obtained velocity changes vary from −1% to 3% with a dominant period of about 1 year. We calculate areal strains by using Global Navigation Satellite Systems data at the volcano. It is evident that the velocity changes are well correlated with the areal strain changes whose linear trend is removed. The stress sensitivity of velocity changes, which are obtained from the observed velocity changes and areal strains, is lower at 0.5–1 Hz than that at higher frequencies: (7.1 ± 1.3) × 10−8 Pa−1 at 0.5–1 Hz, (1.4 ± 0.1) × 10−7 Pa−1 at 1–2 Hz, and (1.3 ± 0.1) × 10−7 Pa−1 at 2–4 Hz. Modeling the velocity changes of Rayleigh wave propagating in layered structures, we find that the observed velocity changes are concentrated in the upper 1 km of the structure. Since similar frequency dependency is recognized among the stress sensitivity of velocity changes reported in previous studies, we conclude that the velocity changes are localized in the shallow depth because of the increase of confining pressure.

We investigate seismic velocity changes in response to the tidal strain at Izu‐Oshima volcano, Japan, by analyzing the data of permanent seismic stations and a small seismic array to evaluate the characteristics of strain sensitivity of velocity changes. We estimate the seismic velocity changes by phase differences between cross‐correlations functions of ambient noises at the frequency of 2–4 Hz stacked for time periods with different tidal strain amplitudes. The seismic velocity changes decrease and increase during dilatation and contraction periods, respectively, when analyzing the cross‐correlations functions at early lapse times ranging from 2 to 7 s. The strain sensitivity of seismic velocity changes is estimated to be ()−2.1±0.2×104 at the early lapse times. However, we find that strain sensitivity of the seismic velocity changes decreases when analyzing the cross‐correlation functions at later lapse times from 7 s to 35 s. Applying an array analysis to the cross‐correlation functions, we observe apparent velocities of about 1 km/s at the early lapse times and those of higher than 1 km/s at the late lapse times. Since the group velocity of Rayleigh waves is 1.1 km/s at Izu‐Oshima volcano, the apparent velocities at the late lapse times may indicate the scattered or reflected body waves incident from a deeper region. Decrease of strain sensitivity with the lapse times therefore results from the emergence of body waves on the late lapse times. These results highlight the need to pay attention to wave types of cross‐correlation functions and their paths to interpret seismic velocity changes.

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.