S U M M A R YSeismic interferometry, also known as Green's function retrieval by crosscorrelation, has a wide range of applications, ranging from surface-wave tomography using ambient noise, to creating virtual sources for improved reflection seismology. Despite its successful applications, the crosscorrelation approach also has its limitations. The main underlying assumptions are that the medium is lossless and that the wavefield is equipartitioned. These assumptions are in practice often violated: the medium of interest is often illuminated from one side only, the sources may be irregularly distributed, and losses may be significant. These limitations may partly be overcome by reformulating seismic interferometry as a multidimensional deconvolution (MDD) process. We present a systematic analysis of seismic interferometry by crosscorrelation and by MDD. We show that for the non-ideal situations mentioned above, the correlation function is proportional to a Green's function with a blurred source. The source blurring is quantified by a so-called interferometric point-spread function which, like the correlation function, can be derived from the observed data (i.e. without the need to know the sources and the medium). The source of the Green's function obtained by the correlation method can be deblurred by deconvolving the correlation function for the point-spread function. This is the essence of seismic interferometry by MDD. We illustrate the crosscorrelation and MDD methods for controlled-source and passive-data applications with numerical examples and discuss the advantages and limitations of both methods.

The increase in number and strength of shallow induced seismicity connected to the Groningen gas field since 2003 and the occurrence of a M L 3.6 event in 2012 started the development of a full probabilistic seismic hazard assessment (PSHA) for Groningen, required by the regulator. Densification of the monitoring network resulted in a decrease of the location threshold and magnitude of completeness down to ∼ M L = 0.5. Combined with a detailed local velocity model, epicentre accuracy could be reduced from 0.5-1 km to 0.1-0.3 km and a vertical resolution ∼0.3 km. Time-dependent seismic activity is observed and taken into account into PSHA calculations. Development of the Ground Motion Model for Groningen resulted in a significant reduction of the hazard. Comparison of different implementations of the PSHA, using different source models, based on either a compaction model and production scenarios or on extrapolation of past seismicity, and methods of calculation, shows similar results.

We introduce seismic interferometry of passive data by multidimensional deconvolution ͑MDD͒ as an alternative to the crosscorrelation method. Interferometry by MDD has the potential to correct for the effects of source irregularity, assuming the first arrival can be separated from the full response. MDD applications can range from reservoir imaging using microseismicity to crustal imaging with teleseismic data.

The last few years there has been a growing number of body-wave observations in noise records. In 1973 Vinnik conjectured that P-waves would even be the dominant wavemode, at epicentral distances of about 40 degrees and onwards from an oceanic source. At arrays far from offshore storms, surface waves induced by nearby storms would not mask the body-wave signal and hence primarily P-waves would be recorded. We measured at such an array in Egypt and indeed found a large proportion of P-waves. At the same time, a new methodology is under development to characterize the lithosphere below an array of receivers, without active sources or local earthquakes. Instead, transmitted waves are used which are caused by distant sources. These sources may either be transient or more stationary. With this new methodology, called seismic interferometry, reflection responses are extracted from the coda of transmissions. Combining the two developments it is clear that there is a large potential for obtaining reflection responses from low-frequency noise. A potential practical advantage of using noise instead of earthquake responses would be that an array only needs to be deployed for a few days or weeks instead of months, to gather enough illumination. We used a few days of continuous noise, recorded with an array in the Abu Gharadig basin, Egypt. We split up the record in three distinct frequency bands and in many small time windows. Using array techniques and taking advantage of all three-component recordings we could unravel the dominant wavemodes arriving in each time window and in each frequency band. The recorded wavemodes, and hence the noise sources, varied significantly per frequency band, and -to a lesser extent-per time window. Primarily P-waves were detected on the vertical component for two of the three frequency bands. For these frequency bands, we only selected the time windows with a favorable illumination. By subsequently applying seismic interferometry, we retrieved P-wave reflection responses and delineated reflectors in the crust, the Moho and possibly the Lehmann discontinuity.

[1] A number of seismic methods exist to image the lithosphere below a collection of receivers, using distant earthquakes. In the current practice, especially mode-conversions in teleseismic phases are utilized. We present a new method that takes advantage of the availability of global phases. This method is called global-phase seismic interferometry (GloPSI). With GloPSI, zero-offset reflections are extracted from reverberations near the array caused by global seismicity. We exemplify GloPSI with data from the Hi-CLIMB experiment (2002)(2003)(2004)(2005) and migrate the obtained reflection responses. This results in a 800 km long reflectivity profile through the Himalayas and a large part of the Tibetan Plateau.Citation: Ruigrok, E., and K. Wapenaar (2012), Global-phase seismic interferometry unveils P-wave reflectivity below the

S U M M A R YIn recent years, there has been an increase in the deployment of relatively dense arrays of seismic stations. The availability of spatially densely sampled global and regional seismic data has stimulated the adoption of industry-style imaging algorithms applied to converted-and scattered-wave energy from distant earthquakes, leading to relatively high-resolution images of the lower crust and upper mantle. We use seismic interferometry to extract reflection responses from the coda of transmitted energy from distant earthquakes. In theory, higher-resolution images can be obtained when migrating reflections obtained with seismic interferometry rather than with conversions, traditionally used in lithospheric imaging methods. Moreover, reflection data allow the straightforward application of algorithms previously developed in exploration seismology. In particular, the availability of reflection data allows us to extract from it a velocity model using standard multichannel data-processing methods. However, the success of our approach relies mainly on a favourable distribution of earthquakes. In this paper, we investigate how the quality of the reflection response obtained with interferometry is influenced by the distribution of earthquakes and the complexity of the transmitted wavefields. Our analysis shows that a reasonable reflection response could be extracted if (1) the array is approximately aligned with an active zone of earthquakes, (2) different phase responses are used to gather adequate angular illumination of the array and (3) the illumination directions are properly accounted for during processing. We illustrate our analysis using a synthetic data set with similar illumination and source-side reverberation characteristics as field data recorded during the 2000-2001 Laramie broad-band experiment. Finally, we apply our method to the Laramie data, retrieving reflection data. We extract a 2-D velocity model from the reflections and use this model to migrate the data. On the final reflectivity image, we observe a discontinuity in the reflections. We interpret this discontinuity as the Cheyenne Belt, a suture zone between Archean and Proterozoic terranes.

Recent developments in the densification of the seismic network covering the Groningen gas field allow a more detailed study of the connection between induced seismicity and reactivated faults around the gas reservoir at 3 km depth. With the reduction of the average station distance from 20 km to 4–5 km, a probabilistic full-waveform moment tensor inversion procedure could be applied, resulting in both improved hypocenter location accuracy and full moment tensor solutions for events of M≥2.0 recorded in the period 2016–2019. Hypocenter locations as output from the moment tensor inversion are compared to locations from the application of other methods and are found similar within 250 m distance. Moment tensor results show that the double-couple (DC) solutions are in accordance with the known structure, namely normal faulting along 50°–70° dipping faults. Comparison with reprocessed 3D seismic sections, extended to a depth of 6–7 km, demonstrate that (a) most events occur along faults with a small throw and (b) reactivated faults in the reservoir often continue downward in the Carboniferous underburden. From non-DC contributions, the isotropic (ISO) component is dominant and shows consistent negative values, which is expected in a compacting medium. There is some indication that events connected to faults with a large throw (>70  m) exhibit the largest ISO component (40%–50%).

A B S T R A C TProgress in the imaging of the mantle and core is partially limited by the sparse distribution of natural sources; the earthquake hypocenters are mainly along the active lithospheric plate boundaries. This problem can be approached with seismic interferometry. In recent years, there has been considerable progress in the development of seismic interferometric techniques. The term seismic interferometry refers to the principle of generating new seismic responses by cross-correlating seismic observations at different receiver locations. The application of interferometric techniques on a global scale could create sources at locations where no earthquakes occur. In this way, yet unknown responses would become available for the application of travel-time tomography and surface-wave dispersion studies. The retrieval of a dense-enough sampling of source gathers would largely benefit the application of reflection imaging.We derive new elastodynamic representation integrals for global-scale seismic interferometry. The relations are different from other seismic interferometry relations for transient sources, in the sense that they are suited for a rotating closed system like the Earth. We use a correlation of an observed response with a response to which free-surface multiple elimination has been applied to account for the closed system. Despite the fact that the rotation of the Earth breaks source-receiver reciprocity, the seismic interferometry relations are shown to be valid. The Coriolis force is included without the need to evaluate an extra term.We synthesize global-scale earthquake responses and use them to illustrate the acoustic versions of the new interferometric relations. When the sampling of real source locations is dense enough, then both the responses with and without freesurface multiples are retrieved. When we do not take into account the responses from the sources in the direct neighborhood of the seismic interferometry-constructed source location, the response with free-surface multiples can still be retrieved. Even when only responses from sources at a certain range of epicentral distances are available, some events in the Green's function between two receiver locations can still be retrieved. The retrieved responses are not perfect, but the artefacts can largely be ascribed to numerical errors. The reconstruction of internal events -the response as if there was a source and a receiver on (major) contrasts within the model -could possibly be of use for imaging.With modelling it is possible to discover in which region of the correlation panel stationary phases occur that contribute to the retrieval of events. This knowledge opens up a new way of filtering out undesired events and of discovering whether specific events could be retrieved with a given source-receiver configuration. *