SUMMARY About a decade ago, noise-based monitoring became a key tool in seismology. One of the tools is passive image interferometry (PII), which uses noise correlation functions (NCF) to retrieve seismic velocity variations. Most studies apply PII to vertical components recording oceanic low-frequent ambient noise ( < 1 Hz). In this work, PII is applied to high-frequent urban ambient noise ( > 1 Hz) on three three-component sensors. With environmental sensors inside the subsurface and in the air, we are able to connect observed velocity variations with environmental parameters. Temperatures below 0 °C correlate well with strong shear wave velocity increases. The temperature sensors inside the ground suggest that a frozen layer of less than 5 cm thickness causes apparent velocity increases above 2 % , depending on the channel pair. The observations indicate that the different velocity variation retrieved from the different channel pairs are due to different surface wave responses inherent in the channel pairs. With dispersion curve modelling in a 1-D medium we can verify that surfaces waves of several tens of metres wavelength experience a velocity increase of several percent due to a centimetres thick frozen layer. Moreover, the model verifies that Love waves show larger velocity increases than Rayleigh waves. The findings of this study provide new insights for monitoring with PII. A few days with temperature below 0 °C can already mask other potential targets (e.g. faults or storage sites). Here, we suggest to use vertical components, which is less sensitive to the frozen layer at the surface. If the target is the seasonal freezing, like in permafrost studies, we suggest to use three-component sensors in order to retrieve the Love wave response. This opens the possibility to study other small-scale processes at the shallow subsurface with surface wave responses.
Many applications in seismology rely on the accurate absolute timing of seismograms. However, both seismological land stations and ocean bottom seismometers (OBSs) can be affected by clock errors, which cause the absolute timing of seismograms to deviate from a highly accurate reference time signal, usually provided by GPS satellites. Timing problems can occur in land stations when synchronization with a GPS signal is temporarily or permanently lost. This can give rise to complicated, time-dependent clock drifts relative to GPS time, due to varying environmental conditions. Seismometers at the ocean bottom cannot receive GPS satellite signals, but operate in more stable ambient conditions than land stations. The standard protocol is to synchronize an OBS with a GPS signal immediately before deployment and after recovery. The measured timing deviation, called 'skew', is assumed to have accumulated linearly over the deployment interval, an assumption that is plausible but usually not verifiable. In recent years, cross-correlations of ambient microseismic noise have been put to use for correcting timing errors, but have been limited to interstation distances of at most a few tens of kilometres without reducing the temporal resolution. We apply noise cross-correlations to the evaluation of clock errors in four broad-band land stations and 53 wideband and broad-band OBSs, which were installed on and around the island of La Réunion in the western Indian Ocean during the RHUM-RUM (Réunion Hotspot and Upper Mantle-Réunions Unterer Mantel) experiment. We correlate all three seismic components, plus a hydrophone channel in OBS stations. Daily cross-correlation functions are derived for intermediate distances (∼20 km) for land-to-land station pairs; stable, 10 d stacks are obtained for very large interstation distances up to >300 km for land-to-OBS and OBS-to-OBS configurations. Averaging over multiple station pairs, and up to 16 component pairs per station, improves the accuracy of the method by a factor of four compared to the single-channel approaches of prior studies. The timing accuracy of our method is estimated to be ∼20 ms standard deviation or one sample at a sampling rate of 50 Hz. In land stations, nonlinear clock drifts and clock jumps of up to 6 min are detected and successfully corrected. For 52 out of 53 OBSs, we successfully obtain drift functions over time, which validate the common assumption of linear clock drift. Skew values that were available for 29 of these OBSs are consistent with our independent estimates within their observational error bars. For 23 OBSs that lacked skew measurements, linear OBS clock drifts range between 0.2 and 8.8 ms d−1. In addition to linear drift, three OBSs are affected by clock jumps of ∼1 s, probably indicating a missing sample problem that would otherwise have gone undetected. Thus we demonstrate the routine feasibility of high-accuracy clock corrections in land and OBSs over a wide range of interstation distances.
Abstract Using closely located seismographs at Piñon Flat (PFO), California, for 1 year long record (2015), we estimated the Rayleigh‐to‐Love wave energy ratio in the secondary microseism (0.1–0.35 Hz) in four seasons. Rayleigh wave energy was estimated from a vertical component seismograph. Love wave energy was estimated from rotation seismograms that were derived from a small array at PFO. Derived ratios are 2–2.5, meaning that there is 2–2.5 times more Rayleigh wave energy than Love wave energy at PFO. In our previous study at Wettzell, Germany, this ratio was 0.9–1.0, indicating comparable energy between Rayleigh waves and Love waves. This difference suggests that the Rayleigh‐to‐Love wave ratios in the secondary microseism may differ greatly from region to region. It also implies that an assumption of the diffuse wavefield is not likely to be valid for this low frequency range as the equipartition of energy should make this ratio much closer.
Coda Wave Interferometry analyzes changes in diffuse waveforms when source and receiver are fixed. Changes are often due to changing acoustic wave speed, perhaps a change of temperature. Alternatively, changes can be studied in noise correlations. The advantage is in not having to wait for, or generate, a precisely identical source. But interpretation is problematic; changes in the correlation of noise fields due to an uncontrolled set of sources could be due to the sources, or to the environment. Recent studies of ambient seismic noise correlations have shown apparent stretches (potentially corresponding to wave speed slowing by parts in 10 000) over periods of weeks. Changes seen at La Reunion correlate with subsequent volcanic eruptions. Changes seen at Parkfield, California occur quickly after earthquakes and then relax back over many months. Potential for seismic monitoring motivates a study of the relative effects of source and medium changes on apparent stretching or compression. Here it is shown that long duration coda like character in a waveform permits a highly accurate estimate of stretch, an estimate that is relatively independent of changes in sources. The theory is supported by laboratory measurements and seismic data.
In previous studies, the contribution of Love waves to the primary microseismic noise field was found to be comparable to those of Rayleigh waves. However, so far only few studies analysed both wave types present in this microseismic noise band, which is known to be generated in shallow water and the theoretical understanding has mainly evolved for Rayleigh waves only. Here, we study the relevance of different source region parameters on the observed primary microseismic noise levels of Love and Rayleigh waves simultaneously. By means of beamforming and correlation of seismic noise amplitudes with ocean wave heights in the period band between 12 and 15 s, we analysed how source areas of both wave types compare with each other around Europe. The generation effectivity in different source regions was compared to ocean wave heights, peak ocean gravity wave propagation direction and bathymetry. Observed Love wave noise amplitudes correlate comparably well with near coastal ocean wave parameters as Rayleigh waves. Some coastal regions serve as especially effective sources for one or the other wave type. These coincide not only with locations of high wave heights but also with complex bathymetry. Further, Rayleigh and Love wave noise amplitudes seem to depend equally on the local ocean wave heights, which is an indication for a coupled variation with swell height during the generation of both wave types. However, the wave-type ratio varies directionally. This observation likely hints towards a spatially varying importance of different source mechanisms or structural influences. Further, the wave-type ratio is modulated depending on peak ocean wave propagation directions which could indicate a variation of different source mechanism strengths but also hints towards an imprint of an effective source radiation pattern. This emphasizes that the inclusion of both wave types may provide more constraints for the understanding of acting generation mechanisms.
Use is free for scientific purposes, provided the aforementioned reference is appropriately cited. The dataset contains raw seismic data in mini-seed format which was used in the publication: doi will follow soon. The zip-file contains the raw data of three broadband stations and their inventories as xml-files. The recording time ranges from 2018-01-01 until 2018-04-19. The data presented here is partly funded through the German Research Foundation (DFG) under Germany‘s Excellence Strategy – EXC 2037 'CLICCS - Climate, Climatic Change, and Society' – Project Number: 390683824, and the Cluster of Excellence 'CliSAP' (EXC177), contribution to the Center for Earth System Research and Sustainability (CEN) of Universität Hamburg.
Seismic velocity changes and nonvolcanic tremor activity in the Parkfield area in California reveal that large earthquakes induce long-term perturbations of crustal properties in the San Andreas fault zone. The 2003 San Simeon and 2004 Parkfield earthquakes both reduced seismic velocities that were measured from correlations of the ambient seismic noise and induced an increased nonvolcanic tremor activity along the San Andreas fault. After the Parkfield earthquake, velocity reduction and nonvolcanic tremor activity remained elevated for more than 3 years and decayed over time, similarly to afterslip derived from GPS (Global Positioning System) measurements. These observations suggest that the seismic velocity changes are related to co-seismic damage in the shallow layers and to deep co-seismic stress change and postseismic stress relaxation within the San Andreas fault zone.
Abstract Our knowledge of the origin of Love waves in the ambient seismic noise is extremely limited. This applies in particular to constraints on source locations and source mechanisms for Love waves in the secondary microseism. Here three‐component beamforming is used to distinguish between the differently polarized wave types in the primary and secondary microseismic noise fields, recorded at several arrays across Europe. We compare characteristics of Love and Rayleigh wave noise, such as source directions and frequency content, measure Love to Rayleigh wave ratios for different back azimuths, and look at the seasonal behavior of our measurements by using a full year of data in 2013. The beamforming results confirm previous observations that back azimuths for Rayleigh and Love waves in both microseismic bands mainly coincide. However, we observe differences in relative directional noise strength between both wave types for the primary microseism. At those frequencies, Love waves dominate on average, with kinetic Love‐to‐Rayleigh energy ratios ranging from 0.6 to 2.0. In the secondary microseism, the ratios are lower, between 0.4 and 1.2. The wave type ratio is directionally homogeneous, except for locations far from the coast. In the primary microseism, our results support the existence of different generation mechanisms. The contribution of a shear traction‐type source mechanism is likely.
Continuous excitation of isolated noise sources leads to repeating wave arrivals in cross correlations of ambient seismic noise, including throughout their coda. These waves propagate from the isolated sources. We observe this effect on correlation wavefields computed from two years of field data recorded at the Gräfenberg array in Germany and two master stations in Europe. Beamforming the correlation functions in the secondary microseism frequency band reveals repeating waves incoming from distinct directions to the West, which correspond to well-known dominant microseism source locations in the Northeastern Atlantic Ocean. These emerge in addition to the expected acausal and causal correlation wavefield contributions by boundary sources, which are converging onto and diverging from the master station, respectively. Numerical simulations reproduce this observation. We first model a source repeatedly exciting a wavelet, which helps illustrate the fundamental mechanism behind repeated wave generation. Second, we model continuously acting secondary microseism sources and find good agreement with our observations. Our observations and modelling have potentially significant implications for the understanding of correlation wavefields and monitoring of relative velocity changes in particular. Velocity monitoring commonly assumes that only multiply scattered waves, originating from the master station, are present in the coda of the correlation wavefield. We show that repeating waves propagating from isolated noise sources may dominate instead, including the very late coda. Our results imply that in the presence of continuously acting noise sources, which we show is the case for ordinary recordings of ocean microseisms, velocity monitoring assuming scattered waves may be adversely affected with regard to measurement technique, spatial resolution, as well as temporal resolution. We further demonstrate that the very late coda of correlation functions contains useful signal, contrary to the common sentiment that it is dominated by instrument noise.
