On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160)  Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

/pdf/the-observation-of-gravitational-waves-from-a-binary-black-58srpupgg9.pdf

The Observation of Gravitational Waves from a Binary Black Hole Merger

Today, surface-wave analysis is widely adopted for building near-surface S-wave velocity models. The surface-wave method is under continuous and rapid evolution, also thanks to the lively scientific debate among different disciplines, and interest in the technique has increased significantly during the last decade. A comprehensive review of the literature in the main scientific journals provides historical perspective, methodological issues, applications, and most-promising recent approaches. Higher modes in the inversion and retrieval of lateral variations are dealt with in great detail, and the current scientific debate on these topics is reported. A best-practices guideline is also outlined.

/pdf/surface-wave-analysis-for-building-near-surface-velocity-39ze2co0if.pdf

Surface-wave analysis for building near-surface velocity models — Established approaches and new perspectives

Seismic interferometry involves the crosscorrelation of responses at different receivers to obtain the Green’s function between these receivers. For the simple situation of an impulsive plane wave propagating along the x-axis, the crosscorrelation of the responses at two receivers along the x-axis gives the Green’s function of the direct wave between these receivers. When the source function of the plane wave is a transientas in exploration seismology or a noise signalas in passive seismology, then the crosscorrelation gives the Green’s function, convolved with the autocorrelation of the source function. Direct-wave interferometry also holds for 2D and 3D situations, assuming the receivers are surrounded by a uniform distribution of sources. In this case, the main contributions to the retrieved direct wave between the receivers come from sources in Fresnel zones around stationary points. The main application of direct-wave interferometry is the retrieval of seismic surface-wave responses from ambient noise and the subsequent tomographic determination of the surfacewave velocity distribution of the subsurface. Seismic interferometry is not restricted to retrieving direct waves between receivers. In a classic paper, Claerbout shows that the autocorrelation of the transmission response of a layered medium gives the plane-wave reflection response of that medium. This is essentially 1D reflected-wave interferometry. Similarly, the crosscorrelation of the transmission responses, observed at two receivers, of an arbitrary inhomogeneous medium gives the 3D reflection response of that medium. One of the main applications of reflected-wave interferometry is retrieving the seismic reflection response from ambient noise and imaging of the reflectors in the subsurface. A common aspect of direct- and reflected-wave interferometry is that virtual sources are created at positions where there are only receivers without requiring knowledge of the subsurface medium parameters or of the positions of the actual sources.

/pdf/tutorial-on-seismic-interferometry-part-1-basic-principles-22w81goc8a.pdf

Tutorial on seismic interferometry: Part 1 — Basic principles and applications

Ambient noise tomography has proven to be effective in resolving shallow earth structure. We applied ambient noise tomography on a dense seismic array in Long Beach, California. The array was composed of more than 5200 stations with an average spacing close to 100 m. Three weeks of passive ambient noise were crosscorrelated between each station pair, which resulted in more than 13.5 million crosscorrelations within the area. Clear fundamental-mode Rayleigh waves were observed between 0.5 and 4 Hz, which were most sensitive to structure above 1-km depth. For each station pair, we applied frequency-time analysis to determine the phase traveltime dispersion, and, for each frequency, we applied eikonal tomography to determine the Rayleigh wave phase velocity map. The eikonal tomography accounted for ray bending by tracking the wavefront and allowed uncertainties to be estimated through statistical analysis. The compilation of phase velocity maps was then used to invert for 3D shear velocity structure. The inverted model showed clear correlation with the known geologic features such as the shallow south–north velocity dichotomy and a deeper fast anomaly associated with the Newport-Inglewood fault zone. Our results can potentially be used to complement traditional active source studies.

/pdf/high-resolution-3d-shallow-crustal-structure-in-long-beach-5abdjm0any.pdf

High-resolution 3D shallow crustal structure in Long Beach, California: Application of ambient noise tomography on a dense seismic array

[1] Understanding the geotectonic evolution of the southeastern Tibetan plateau requires knowledge about the structure of the lithosphere. Using data from 77 broadband stations in SW China, we invert Rayleigh wave phase velocity dispersion curves from ambient noise interferometry (T =1 0–40 s) and teleseismic surface waves (T =2 0–150 s) for 3‐D heterogeneity and azimuthal anisotropy in the lithosphere to ∼150 km depth. Our surface wave array tomography reveals (1) deep crustal zones of anomalously low shear wave speed and (2) substantial variations with depth of the pattern of azimuthal anisotropy. Upper crustal azimuthal anisotropy reveals a curvilinear pattern around the eastern Himalayan syntaxis, with fast directions generally parallel to the main strike slip faults. The mantle pattern of azimuthal anisotropy is different from that in the crust and varies from north to south. The tomographically inferred 3‐D variation in azimuthal anisotropy helps constrain the source region of shear wave splitting. South of ∼26°N (off the high plateau) most of the observed splitting can be accounted for by upper mantle anisotropy, but for stations on the plateau proper (with thick crust) crustal anisotropy cannot be ignored. On long wavelengths, the pattern of azimuthal anisotropy in the crust differs from that in the mantle. This is easiest explained if deformation varies with depth. The deep crustal zones of low shear wave speed (and, presumably, mechanical strength) may represent loci of ductile deformation. But their lateral variation suggests that in SE Tibet (localized) crustal channel flow and motion along the major strike slip faults are both important.

/pdf/heterogeneity-and-anisotropy-of-the-lithosphere-of-se-tibet-3mpfrp5ys7.pdf

Heterogeneity and anisotropy of the lithosphere of SE Tibet from surface wave array tomography

One application of seismic interferometry is to retrieve the impulse response (Green's function) from crosscorrelation of ambient seismic noise. Various researchers show results for retrieving the surface-wave part of the Green's function. However, reflection retrieval has proven more challenging. We crosscorrelate ambient seismic noise, recorded along eight parallel lines in the Sirte basin east of Ajdabeya, Libya, to obtain shot gathers that contain reflections. We take advantage of geophone groups to suppress part of the undesired surface-wave noise and apply frequency-wavenumber filtering before crosscorrelation to suppress surface waves further. After comparing the retrieved results with data from an active seismic exploration survey along the same lines, we use the retrieved reflection data to obtain a migrated reflection image of the subsurface.

/pdf/reflection-images-from-ambient-seismic-noise-17u8e35gi4.pdf

Reflection images from ambient seismic noise

The successful surface waves retrieval in solid-Earth seismology using long-time correlations and subsequent tomographic images of the crust have sparked interest in extraction of subsurface information from noise in the exploration seismology. Subsurface information in exploration seismology is usually derived from body-wave reflections > 1 Hz, which is challenging for utilization of ambient noise. We use 11 h of noise recorded in the Sirte basin, Libya. First, we study the characteristics of the noise. We show that the bulk of the noise is composed of surface waves at frequencies below 6 Hz. Some noise panels contain nearly vertically traveling events. We further characterize these events using a beamforming algorithm. From the beamforming, we conclude that these events represent body-wave arrivals with a fairly rich azimuthal distribution. Having body-wave arrivals in the noise is a prerequisite for body-wave reflections retrieval. We crosscorrelate and sum the recorded ambient-noise panels to retrieve common-source gathers, following two approaches - using all the noise and using only noise panels containing body-wave arrivals likely to contribute to the reflections retrieval. Comparing the retrieved gathers with active seismic data, we show that the two-way traveltimes at short offsets of several retrieved events coincide with those of reflections in the active data and thus correspond to apexes of reflections. We then compare retrieved stacked sections of the subsurface from both approaches with the active-data stacked section and show that the reflectors are consistent along a line. The results from the second approach exhibit the reflectors better. © 2013. American Geophysical Union. All Rights Reserved.

/pdf/seismic-exploration-scale-velocities-and-structure-from-1hdlwp6ahz.pdf

Seismic exploration‐scale velocities and structure from ambient seismic noise (>1 Hz)

https://repository.tudelft.nl/islandora/object/uuid%3Ae23dede0-cc57-464e-ad42-77b49a1bb6e1/datastream/OBJ/download

Extraction of P-wave reflections from microseisms

In 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.

/pdf/high-resolution-lithospheric-imaging-with-seismic-11i899bu0q.pdf

Xander Campman

Papers

Tutorial on seismic interferometry: Part 1 — Basic principles and applications

Reflection images from ambient seismic noise

Seismic exploration‐scale velocities and structure from ambient seismic noise (>1 Hz)

Extraction of P-wave reflections from microseisms

High-resolution lithospheric imaging with seismic interferometry