06.08.2013 Views

Lawrence et al., 2012: DRAFT: JGR: Ambient Noise Numerical ...

Lawrence et al., 2012: DRAFT: JGR: Ambient Noise Numerical ...

Lawrence et al., 2012: DRAFT: JGR: Ambient Noise Numerical ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Lawrence</strong> <strong>et</strong> <strong>al</strong>., <strong>2012</strong>: <strong>DRAFT</strong>: <strong>JGR</strong>: <strong>Ambient</strong> <strong>Noise</strong> Numeric<strong>al</strong> Ev<strong>al</strong>uation<br />

We explore the param<strong>et</strong>er domain within which equation 10 is v<strong>al</strong>id to characterize the<br />

appropriate ambient seismic field. The three dimension sc<strong>al</strong>es, the source-length sc<strong>al</strong>e Ls,<br />

the receiver-length sc<strong>al</strong>e Lx and the source-receiver separation Δsx reflect the radiation<br />

aperture with the ratios Ls/Δsx and Ls/Lx (i.e, Figure 1). Another factor to consider is the<br />

number of sources, Ns, relative to the magnitude distribution, Ms. The larger the range of<br />

source amplitudes used to generate the synth<strong>et</strong>ic spectra the fewer the sources dominate<br />

any given seismogram. Consequently the fraction provides a measure of<br />

the number of sources contributing to the spectra. Fin<strong>al</strong>ly, Due to geom<strong>et</strong>ric<strong>al</strong> spreading,<br />

the range of source-to-receiver distances <strong>al</strong>so plays an important role in the number of<br />

dominant sources: [max(rsx)-min(rsx)].<br />

For a given scenario, we generate synth<strong>et</strong>ic ambient noise seismograms cut into<br />

1800s time windows, as found to be optim<strong>al</strong> for fastest convergence of the NCF toward<br />

stable NCF [Seats <strong>et</strong> <strong>al</strong>., 2011]. We then compute the coherencies (in frequency domain)<br />

or NCFs (in time domain) and stack Re[γxy(ω)] over <strong>al</strong>l station pairs to estimate the<br />

coherency at various distances. For each frequency, we estimate phase velocity and<br />

attenuation coefficients and we fin<strong>al</strong>ly ev<strong>al</strong>uate the match b<strong>et</strong>ween the synth<strong>et</strong>ic, or<br />

“observed” coherency, and the predicted from equation 2. Table 1 summarizes the<br />

param<strong>et</strong>ers to simulate the various scenarios.<br />

The synth<strong>et</strong>ic noise seismograms are a good visu<strong>al</strong> control of similarity b<strong>et</strong>ween the<br />

generated noise and the true ambient seismic field. Figure 2 shows examples of the<br />

generated ambient seismic field with different noise source distribution and the resulting<br />

variability in the sign<strong>al</strong>s. Due to the stochastic nature of our source generator, we can<br />

describe the statistics related to seismograms, which change as a function of source-<br />

receiver geom<strong>et</strong>ry and number of sources per time window.<br />

Figure 2 presents example synth<strong>et</strong>ic noise seismograms for the four source<br />

distributions illustrated in Figure 1, namely circular, linear, uniform, and uniform shell.<br />

We imposed 128 sources per time window (Ns=128) because it provides a reasonable<br />

b<strong>al</strong>ance b<strong>et</strong>ween synth<strong>et</strong>ics seismograms with identifiable phases above the noise and<br />

those without. With 128 sources per time window, Δsx=500km, Ls=1000km, and<br />

Lx=100km, only 14 % of the time windows contain events that have a short-term<br />

11<br />

Marine Denolle 12/31/12 7:10 PM<br />

Comment: Do we play with this<br />

param<strong>et</strong>er?

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!