USGS Science⌕ Search

Geology topics

Matthew M. Haney

Publications and source records attributed to Matthew M. Haney.

At least 91 records · Page 5Linked to original sources

Long period seismicity and very long period infrasound driven by shallow magmatic degassing at Mount Pagan, Mariana Islands

Long period (LP) seismicity and very long period infrasound (iVLP) were recorded during continuous degassing from Mount Pagan, Mariana Islands, in July 2013 to January 2014. The frequency content of the LP and iVLP events and delay times between the two arrivals were remarkably stable and indicate nearly co-located sources. Using phase-weighted stacking over similar events to dampen noise, we find that the LP source centroid is located 60 m below and 180 m west of the summit vent. The moment tensor reveals a volumetric source modeled as resonance of a subhorizontal sill intersecting a dike. We model the seismoacoustic wavefields with a coupled earth-air 3-D finite difference code. The ratios of pressure to velocity measured at the infrasound arrays are an order of magnitude larger than the synthetic ratios, so the iVLP is not the result of LP energy transmitting into the atmosphere at its epicenter. Based on crater shape and dimensions determined by structure from motion, we model the iVLP as acoustic resonance of an exponential horn. The source of the continuous plume from gas analysis is shallow magmatic degassing, which repeatedly pressurized the dike-sill portion of the conduit over the 7 months of observation. Periodic gas release caused the geologically controlled sill to partially collapse and resonate, while venting of gas at the surface triggered resonance in the crater. LP degassing only accounts for ~12% of total degassing, indicating that most degassing is relatively aseismic and that multiple active pathways exist beneath the vent.

Mount Pagan↗

Source mechanism of small long-period events at Mount St. Helens in July 2005 using template matching, phase-weighted stacking, and full-waveform inversion

Long-period (LP, 0.5-5 Hz) seismicity, observed at volcanoes worldwide, is a recognized signature of unrest and eruption. Cyclic LP “drumbeating” was the characteristic seismicity accompanying the sustained dome-building phase of the 2004–2008 eruption of Mount St. Helens (MSH), WA. However, together with the LP drumbeating was a near-continuous, randomly occurring series of tiny LP seismic events (LP “subevents”), which may hold important additional information on the mechanism of seismogenesis at restless volcanoes. We employ template matching, phase-weighted stacking, and full-waveform inversion to image the source mechanism of one multiplet of these LP subevents at MSH in July 2005. The signal-to-noise ratios of the individual events are too low to produce reliable waveform-inversion results, but the events are repetitive and can be stacked. We apply network-based template matching to 8 days of continuous velocity waveform data from 29 June to 7 July 2005 using a master event to detect 822 network triggers. We stack waveforms for 359 high-quality triggers at each station and component, using a combination of linear and phase-weighted stacking to produce clean stacks for use in waveform inversion. The derived source mechanism pointsto the volumetric oscillation (~10 m 3 ) of a subhorizontal crack located at shallow depth (~30 m) in an area to the south of Crater Glacier in the southern portion of the breached MSH crater. A possible excitation mechanism is the sudden condensation of metastable steam from a shallow pressurized hydrothermal system as it encounters cool meteoric water in the outer parts of the edifice, perhaps supplied from snow melt.

Washington↗

Monitoring changes in seismic velocity related to an ongoing rapid inflation event at Okmok volcano, Alaska

Okmok is one of the most active volcanoes in the Aleutian Arc. In an effort to improve our ability to detect precursory activity leading to eruption at Okmok, we monitor a recent, and possibly ongoing, GPS-inferred rapid inflation event at the volcano using ambient noise interferometry (ANI). Applying this method, we identify changes in seismic velocity outside of Okmok’s caldera, which are related to the hydrologic cycle. Within the caldera, we observe decreases in seismic velocity that are associated with the GPS-inferred rapid inflation event. We also determine temporal changes in waveform decorrelation and show a continual increase in decorrelation rate over the time associated with the rapid inflation event. Themagnitude of relative velocity decreases and decorrelation rate increases are comparable to previous studies at Piton de la Fournaise that associate such changes with increased production of volatiles and/ormagmatic intrusion within the magma reservoir and associated opening of fractures and/or fissures. Notably, the largest decrease in relative velocity occurs along the intrastation path passing nearest to the center of the caldera. This observation, along with equal amplitude relative velocity decreases revealed via analysis of intracaldera autocorrelations, suggests that the inflation sourcemay be located approximately within the center of the caldera and represent recharge of shallow magma storage in this location. Importantly, there is a relative absence of seismicity associated with this and previous rapid inflation events at Okmok. Thus, these ANI results are the first seismic evidence of such rapid inflation at the volcano.

Alaska↗

A comparison of methods to estimate seismic phase delays--Numerical examples for coda wave interferometry

Time-shift estimation between arrivals in two seismic traces before and after a velocity perturbation is a crucial step in many seismic methods. The accuracy of the estimated velocity perturbation location and amplitude depend on this time shift. Windowed cross correlation and trace stretching are two techniques commonly used to estimate local time shifts in seismic signals. In the work presented here, we implement Dynamic Time Warping (DTW) to estimate the warping function – a vector of local time shifts that globally minimizes the misfit between two seismic traces. We illustrate the differences of all three methods compared to one another using acoustic numerical experiments. We show that DTW is comparable to or better than the other two methods when the velocity perturbation is homogeneous and the signal-to-noise ratio is high. When the signal-to-noise ratio is low, we find that DTW and windowed cross correlation are more accurate than the stretching method. Finally, we show that the DTW algorithm has better time resolution when identifying small differences in the seismic traces for a model with an isolated velocity perturbation. These results impact current methods that utilize not only time shifts between (multiply) scattered waves, but also amplitude and decoherence measurements. DTW is a new tool that may find new applications in seismology and other geophysical methods (e.g., as a waveform inversion misfit function).

Geophysical Journal International↗

Point spread functions for earthquake source imaging: An interpretation based on seismic interferometry

Recently, various methods have been proposed and applied for earthquake source imaging, and theoretical relationships among the methods have been studied. In this study, we make a follow-up theoretical study to better understand the meanings of earthquake source imaging. For imaging problems, the point spread function (PSF) is used to describe the degree of blurring and degradation in an obtained image of a target object as a response of an imaging system. In this study, we formulate PSFs for earthquake source imaging. By calculating the PSFs, we find that waveform source inversion methods remove the effect of the PSF and are free from artifacts. However, the other source imaging methods are affected by the PSF and suffer from the effect of blurring and degradation due to the restricted distribution of receivers. Consequently, careful treatment of the effect is necessary when using the source imaging methods other than waveform inversions. Moreover, the PSF for source imaging is found to have a link with seismic interferometry with the help of the source-receiver reciprocity of Green’s functions. In particular, the PSF can be related to Green’s function for cases in which receivers are distributed so as to completely surround the sources. Furthermore, the PSF acts as a low-pass filter. Given these considerations, the PSF is quite useful for understanding the physical meaning of earthquake source imaging.

Geophysical Journal International↗

Non-perturbational surface-wave inversion: A Dix-type relation for surface waves

We extend the approach underlying the well-known Dix equation in reflection seismology to surface waves. Within the context of surface wave inversion, the Dix-type relation we derive for surface waves allows accurate depth profiles of shear-wave velocity to be constructed directly from phase velocity data, in contrast to perturbational methods. The depth profiles can subsequently be used as an initial model for nonlinear inversion. We provide examples of the Dix-type relation for under-parameterized and over-parameterized cases. In the under-parameterized case, we use the theory to estimate crustal thickness, crustal shear-wave velocity, and mantle shear-wave velocity across the Western U.S. from phase velocity maps measured at 8-, 20-, and 40-s periods. By adopting a thin-layer formalism and an over-parameterized model, we show how a regularized inversion based on the Dix-type relation yields smooth depth profiles of shear-wave velocity. In the process, we quantitatively demonstrate the depth sensitivity of surface-wave phase velocity as a function of frequency and the accuracy of the Dix-type relation. We apply the over-parameterized approach to a near-surface data set within the frequency band from 5 to 40 Hz and find overall agreement between the inverted model and the result of full nonlinear inversion.

Geophysics↗

Tracking changes in volcanic systems with seismic Interferometry

The detection and evaluation of time-dependent changes at volcanoes form the foundation upon which successful volcano monitoring is built. Temporal changes at volcanoes occur over all time scales and may be obvious (e.g., earthquake swarms) or subtle (e.g., a slow, steady increase in the level of tremor). Some of the most challenging types of time-dependent change to detect are subtle variations in material properties beneath active volcanoes. Although difficult to measure, such changes carry important information about stresses and fluids present within hydrothermal and magmatic systems. These changes are imprinted on seismic waves that propagate through volcanoes. In recent years, there has been a quantum leap in the ability to detect subtle structural changes systematically at volcanoes with seismic waves. The new methodology is based on the idea that useful seismic signals can be generated “at will” from seismic noise. This means signals can be measured any time, in contrast to the often irregular and unpredictable times of earthquakes. With seismic noise in the frequency band 0.1–1 Hz arising from the interaction of the ocean with the solid Earth known as microseisms, researchers have demonstrated that cross-correlations of passive seismic recordings between pairs of seismometers yield coherent signals (Campillo and Paul 2003; Shapiro and Campillo 2004). Based on this principle, coherent signals have been reconstructed from noise recordings in such diverse fields as helioseismology (Rickett and Claerbout 2000), ultrasound (Weaver and Lobkis 2001), ocean acoustic waves (Roux and Kuperman 2004), regional (Shapiro et al. 2005; Sabra et al. 2005; Bensen et al. 2007) and exploration (Draganov et al. 2007) seismology, atmospheric infrasound (Haney 2009), and studies of the cryosphere (Marsan et al. 2012). Initial applications of ambient seismic noise were to regional surface wave tomography (Shapiro et al. 2005). Brenguier et al. (2007) were the first to use ambient noise tomography (ANT) to map the 3D structure of a volcanic interior (at Piton de la Fournaise). Subsequent studies have imaged volcanoes with ANT at Okmok (Masterlark et al. 2010), Toba (Stankiewicz et al. 2010), Katmai (Thurber et al. 2012), Asama (Nagaoka et al. 2012), Uturuncu (Jay et al. 2012), and Kilauea (Ballmer et al. 2013b). In addition, Ma et al. (2013) have imaged a scatterer in the volcanic region of southern Peru by applying array techniques to ambient noise correlations. Prior to and in tandem with the development of ANT, researchers discovered that repeating earthquakes, which often occur at volcanoes, could be used to monitor subtle time-dependent changes with a technique known as the doublet method or coda wave interferometry (CWI) (Poupinet et al. 1984; Roberts et al. 1992; Ratdomopurbo and Poupinet 1995; Snieder et al. 2002; Pandolfi et al. 2006; Wegler et al. 2006; Martini et al. 2009; Haney et al. 2009; De Angelis 2009; Nagaoka et al. 2010; Battaglia et al. 2012; Erdem and Waite 2005; Hotovec-Ellis et al. 2014). Chaput et al. (2012) have also used scattered waves from Strombolian eruption coda at Erebus volcano to image the reflectivity of the volcanic interior with body wave interferometry. However, CWI in its original form was limited in that repeating earthquakes, or doublets, were not always guaranteed to occur. With the widespread use of noise correlations in seismology following the groundbreaking work by Campillo and Paul (2003) and Shapiro et al. (2005), it became evident that the nature of the ambient seismic field, due to its oceanic origin, enabled the continuous monitoring of subtle, time-dependent changes at both fault zones (Wegler and Sens-Schönfelder 2007; Brenguier et al. 2008b; Wegler et al. 2009; Sawazaki et al. 2009; Tatagi et al. 2012) and volcanoes (Sens-Schönfelder and Wegler 2006; Brenguier et al. 2008a) without the need for repeating earthquakes. Seismic precursors to eruptions based on ambient noise we

Book chapter↗

The 2013 eruption of Pavlof Volcano, Alaska: a spatter eruption at an ice- and snow-clad volcano

The 2013 eruption of Pavlof Volcano, Alaska began on 13 May and ended 49 days later on 1 July. The eruption was characterized by persistent lava fountaining from a vent just north of the summit, intermittent strombolian explosions, and ash, gas, and aerosol plumes that reached as high as 8 km above sea level and on several occasions extended as much as 500 km downwind of the volcano. During the first several days of the eruption, accumulations of spatter near the vent periodically collapsed to form small pyroclastic avalanches that eroded and melted snow and ice to form lahars on the lower north flank of the volcano. Continued lava fountaining led to the production of agglutinate lava flows that extended to the base of the volcano, about 3–4 km beyond the vent. The generation of fountain-fed lava flows was a dominant process during the 2013 eruption; however, episodic collapse of spatter accumulations and formation of hot spatter-rich granular avalanches was a more efficient process for melting snow and ice and initiating lahars. The lahars and ash plumes generated during the eruption did not pose any serious hazards for the area. However, numerous local airline flights were cancelled or rerouted, and trace amounts of ash fall occurred at all of the local communities surrounding the volcano, including Cold Bay, Nelson Lagoon, Sand Point, and King Cove.

Alaska↗

Surface‐wave Green’s tensors in the near field

We demonstrate the connection between theoretical expressions for the correlation of ambient noise Rayleigh and Love waves and the exact surface‐wave Green’s tensors for a point force. The surface‐wave Green’s tensors are well known in the far‐field limit. On the other hand, the imaginary part of the exact Green’s tensors, including near‐field effects, arises in correlation techniques such as the spatial autocorrelation (SPAC) method. Using the imaginary part of the exact Green’s tensors from the SPAC method, we find the associated real part using the Kramers–Kronig relations. The application of the Kramers–Kronig relations is not straightforward, however, because the causality properties of the different tensor components vary. In addition to the Green’s tensors for a point force, we also derive expressions for a general point moment tensor source.

Bulletin of the Seismological Society of America↗

Backprojection of volcanic tremor

Backprojection has become a powerful tool for imaging the rupture process of global earthquakes. We demonstrate the ability of backprojection to illuminate and track volcanic sources as well. We apply the method to the seismic network from Okmok Volcano, Alaska, at the time of an escalation in tremor during the 2008 eruption. Although we are able to focus the wavefield close to the location of the active cone, the network array response lacks sufficient resolution to reveal kilometer-scale changes in tremor location. By deconvolving the response in successive backprojection images, we enhance resolution and find that the tremor source moved toward an intracaldera lake prior to its escalation. The increased tremor therefore resulted from magma-water interaction, in agreement with the overall phreatomagmatic character of the eruption. Imaging of eruption tremor shows that time reversal methods, such as backprojection, can provide new insights into the temporal evolution of volcanic sources.

Alaska↗

Distinguishing high surf from volcanic long-period earthquakes

Repeating long-period (LP) earthquakes are observed at active volcanoes worldwide and are typically attributed to unsteady pressure fluctuations associated with fluid migration through the volcanic plumbing system. Nonvolcanic sources of LP signals include ice movement and glacial outburst floods, and the waveform characteristics and frequency content of these events often make them difficult to distinguish from volcanic LP events. We analyze seismic and infrasound data from an LP swarm recorded at Pagan volcano on 12–14 October 2013 and compare the results to ocean wave data from a nearby buoy. We demonstrate that although the events show strong similarity to volcanic LP signals, the events are not volcanic but due to intense surf generated by a passing typhoon. Seismo-acoustic methods allow for rapid distinction of volcanic LP signals from those generated by large surf and other sources, a critical task for volcano monitoring.

Pagan Volcano↗

Shaking up volcanoes

Most volcanic eruptions that occur shortly after a large distant earthquake do so by random chance. A few compelling cases for earthquake-triggered eruptions exist, particularly within 200 km of the earthquake, but this phenomenon is rare in part because volcanoes must be poised to erupt in order to be triggered by an earthquake (1). Large earthquakes often perturb volcanoes in more subtle ways by triggering small earthquakes and changes in spring discharge and groundwater levels (1, 2). On page 80 of this issue, Brenguier et al. (3) provide fresh insight into the interaction of large earthquakes and volcanoes by documenting a temporary change in seismic velocity beneath volcanoes in Honshu, Japan, after the devastating Tohoku-Oki earthquake in 2011.

Science↗

Ambient seismic noise interferometry in Hawai'i reveals long-range observability of volcanic tremor

The use of seismic noise interferometry to retrieve Green's functions and the analysis of volcanic tremor are both useful in studying volcano dynamics. Whereas seismic noise interferometry allows long-range extraction of interpretable signals from a relatively weak noise wavefield, the characterization of volcanic tremor often requires a dense seismic array close to the source. We here show that standard processing of seismic noise interferometry yields volcanic tremor signals observable over large distances exceeding 50 km. Our study comprises 2.5 yr of data from the U.S. Geological Survey Hawaiian Volcano Observatory short period seismic network. Examining more than 700 station pairs, we find anomalous and temporally coherent signals that obscure the Green's functions. The time windows and frequency bands of these anomalous signals correspond well with the characteristics of previously studied volcanic tremor sources at Pu'u 'Ō'ō and Halema'uma'u craters. We use the derived noise cross-correlation functions to perform a grid-search for source location, confirming that these signals are surface waves originating from the known tremor sources. A grid-search with only distant stations verifies that useful tremor signals can indeed be recovered far from the source. Our results suggest that the specific data processing in seismic noise interferometry—typically used for Green's function retrieval—can aid in the study of both the wavefield and source location of volcanic tremor over large distances. In view of using the derived Green's functions to image heterogeneity and study temporal velocity changes at volcanic regions, however, our results illustrate how care should be taken when contamination by tremor may be present.

Hawai'i↗

Source characterization for an explosion during the 2009 eruption of Redoubt Volcano from very-long-period seismic waves

The 2009 eruption of Redoubt produced several very-long-period (VLP) signals associated with explosions. We invert for the source location and mechanism of an explosion at Redoubt volcano using waveform methods applied to broadband recordings. Such characterization of the source carries information on the geometry of the conduit and the physics of the explosion process. Inversions are carried out assuming the volcanic source can be modeled as a point source, with mechanisms described by a) a set of 3 orthogonal forces, b) a moment tensor consisting of force couples, and c) both forces and moment tensor components. We find that the source of the VLP seismic waves during the explosion is well-described by either a combined moment/force source located northeast of the crater and at an elevation of 1.6 km ASL or a moment source at an elevation of 800 m to the southwest of the crater. The moment tensors for the solutions with moment and force and moment-only share similar characteristics. The source time functions for both moment tensors begin with inflation (pressurization) and execute two cycles of deflation-reinflation (depressurization–repressurization). Although the moment/force source provides a better fit to the data, we find that owing to the limited coverage of the broadband stations at Redoubt the moment-only source is the more robust and reliable solution. Based on the moment-only solution, we estimate a volume change of 19,000 m 3 and a pressure change of 7 MPa in a dominant sill and an out-of-phase volume change of 5000 m 3 and pressure change of 1.8 MPa in a subdominant dike at the source location. These results shed new light on the magmatic plumbing system beneath Redoubt and complement previous studies on Vulcanian explosions at other volcanoes.

Alaska↗

Seismic observations of Redoubt Volcano, Alaska - 1989-2010 and a conceptual model of the Redoubt magmatic system

Seismic activity at Redoubt Volcano, Alaska, has been closely monitored since 1989 by a network of five to ten seismometers within 22 km of the volcano's summit. Major eruptions occurred in 1989-1990 and 2009 and were characterized by large volcanic explosions, episodes of lava dome growth and failure, pyroclastic flows, and lahars. Seismic features of the 1989-1990 eruption were 1) weak precursory tremor and a short, 23-hour-long, intense swarm of repetitive shallow long-period (LP) events centered 1.4 km below the crater floor, 2) shallow volcano-tectonic (VT) and hybrid earthquakes that separated early episodes of dome growth, 3) 13 additional swarms of LP events at shallow depths precursory to many of the 25 explosions that occurred over the more than 128 day duration of eruptive activity, and 4) a persistent cluster of VT earthquakes at 6 to 9 km depth. In contrast the 2009 eruption was preceded by a pronounced increase in deep-LP (DLP) events at lower crustal depths (25 to 38 km) that began in mid-December 2008, two months of discontinuous shallow volcanic tremor that started on January 23, 2009, a strong phreatic explosion on March 15, and a 58-hour-long swarm of repetitive shallow LP events. The 2009 eruption consisted of at least 23 major explosions between March 23 and April 5, again accompanied by shallow VT earthquakes, several episodes of shallow repetitive LP events and dome growth continuing until mid July. Increased VT earthquakes at 4 to 9 km depth began slowly in early April, possibly defining a mid-crustal magma source zone. Magmatic processes associated with the 2009 eruption seismically activated the same portions of the Redoubt magmatic system as the 1989-1990 eruption, although the time scales and intensity vary considerably among the two eruptions. The occurrence of precursory DLP events suggests that the 2009 eruption may have involved the rise of magma from lower crustal depths. Based on the evolution of seismicity during the 1989-1990 and 2009 eruptions the Redoubt magmatic system is envisioned to consist of a shallow system of cracks extending 1 to 2 km below the crater floor, a magma storage or source region at roughly 3 to 9 km depth, and a diffuse magma source region at 25 to 38 km depth. Close tracking of seismic activity allowed the Alaska Volcano Observatory to successfully issue warnings prior to many of the hazardous explosive events that occurred in 2009.

Alaska↗

Catalog of earthquake hypocenters at Alaskan volcanoes: January 1 through December 31, 2012

Between January 1 and December 31, 2012, the Alaska Volcano Observatory located 4,787 earthquakes, of which 4,211 occurred within 20 kilometers of the 33 volcanoes monitored by a seismograph network. There was significant seismic activity at Iliamna, Kanaga, and Little Sitkin volcanoes in 2012. Instrumentation highlights for this year include the implementation of the Advanced National Seismic System Quake Monitoring System hardware and software in February 2012 and the continuation of the American Recovery and Reinvestment Act work in the summer of 2012. The operational highlight was the removal of Mount Wrangell from the list of monitored volcanoes. This catalog includes hypocenters, magnitudes, and statistics of the earthquakes located in 2012 with the station parameters, velocity models, and other files used to locate these earthquakes.

Alaska↗

Casual instrument corrections for short-period and broadband seismometers

Of all the filters applied to recordings of seismic waves, which include source, path, and site effects, the one we know most precisely is the instrument filter. Therefore, it behooves seismologists to accurately remove the effect of the instrument from raw seismograms. Applying instrument corrections allows analysis of the seismogram in terms of physical units (e.g., displacement or particle velocity of the Earth’s surface) instead of the output of the instrument (e.g., digital counts). The instrument correction can be considered the most fundamental processing step in seismology since it relates the raw data to an observable quantity of interest to seismologists. Complicating matters is the fact that, in practice, the term “instrument correction” refers to more than simply the seismometer. The instrument correction compensates for the complete recording system including the seismometer, telemetry, digitizer, and any anti‐alias filters. Knowledge of all these components is necessary to perform an accurate instrument correction. The subject of instrument corrections has been covered extensively in the literature (Seidl, 1980; Scherbaum, 1996). However, the prospect of applying instrument corrections still evokes angst among many seismologists—the authors of this paper included. There may be several reasons for this. For instance, the seminal paper by Seidl (1980) exists in a journal that is not currently available in electronic format and cannot be accessed online. Also, a standard method for applying instrument corrections involves the programs TRANSFER and EVALRESP in the Seismic Analysis Code (SAC) package (Goldstein et al. , 2003). The exact mathematical methods implemented in these codes are not thoroughly described in the documentation accompanying SAC.

Seismological Research Letters↗

Extension of the spatial autocorrelation (SPAC) method to mixed-component correlations of surface waves

Using ambient seismic noise for imaging subsurface structure dates back to the development of the spatial autocorrelation (SPAC) method in the 1950s. We present a theoretical analysis of the SPAC method for multicomponent recordings of surface waves to determine the complete 3 × 3 matrix of correlations between all pairs of three-component motions, called the correlation matrix. In the case of isotropic incidence, when either Rayleigh or Love waves arrive from all directions with equal power, the only non-zero off-diagonal terms in the matrix are the vertical–radial (ZR) and radial–vertical (RZ) correlations in the presence of Rayleigh waves. Such combinations were not considered in the development of the SPAC method. The method originally addressed the vertical–vertical (ZZ), RR and TT correlations, hence the name spatial autocorrelation. The theoretical expressions we derive for the ZR and RZ correlations offer additional ways to measure Rayleigh wave dispersion within the SPAC framework. Expanding on the results for isotropic incidence, we derive the complete correlation matrix in the case of generally anisotropic incidence. We show that the ZR and RZ correlations have advantageous properties in the presence of an out-of-plane directional wavefield compared to ZZ and RR correlations. We apply the results for mixed-component correlations to a data set from Akutan Volcano, Alaska and find consistent estimates of Rayleigh wave phase velocity from ZR compared to ZZ correlations. This work together with the recently discovered connections between the SPAC method and time-domain correlations of ambient noise provide further insights into the retrieval of surface wave Green’s functions from seismic noise.

Geophysical Journal International↗