Monte Carlo simulation of seismogram envelopes in scattering media
|
|
- Victor Barton
- 5 years ago
- Views:
Transcription
1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 105, NO. B3, PAGES , MARCH 10, 2000 Monte Carlo simulation of seismogram envelopes in scattering media Kazuo Yoshimoto Department of Geophysics, Graduate School of Science, Tohoku University, Sendai, Japan Abstract. The analysis of the seismogram coda envelopes of local and regional earthquakes is one of the most effective strategies for investigating the heterogeneous lithospheric structure characterized by the seismic scattering and attenuation. In order to synthesize the coda envelope we introduce a numerical scheme called the direct simulation Monte Carlo method, which has been used in the field of the kinetic theory of gases. Because of the simplicity of the algorithm the method has several advantages over previous methods in terms of the flexibility of the numerical calculation to incorporate various factors required to construct realistic seismogram envelopes. On the basis of coda envelope simulations, including multiple scattering, we show that an increase of seismic velocity with depth severely affects the shape of the coda envelope. The effects of ray bending due to the velocity increase at the Moho and the reflection at the free surface are clearly found in the synthesized envelope for a shallow earthquake. Our simulation demonstrates that the amplitude of the envelope is magnified by stagnation of seismic energy at shallow depths due to the positive velocity gradient with depth. Because of this effect, for an a priori assumption of a homogeneous velocity model the measurement of the scattering coefficient by conventional methods may be overestimated. 1. Introduction of the seismic energy just beneath the free surface, as suggested by Hayakawa [1998] for the single isotropic scattering Coda waves of local and regional earthquakes are consid- model. In order to quantify this effect on seismogram enveered to be incoherent successive arrivals of the scattered waves lopes of local earthquakes it is necessary to trace the ray paths from random heterogeneoustructures [Aki end Chouet, 1975]. Aki [1980] indicated that coda waves have the same attenuation characteristics as S waves in the lithosphere. Array observations of high-frequency seismic waves reported that the coda waves consist of scattered body waves, except for the surface waves generated by local near-surface weathered layers [e.g., l/emon et el., 1998]. On the basis of these characteristics the analysis of the seismogram coda envelope has been one of the most effective tools to investigate the heterogeneity of the of seismic waves in three dimensions by using a realistic velocity structure model. The Monte Carlo method plays an important role in investigating the effect of velocity structure on coda envelopes. Adopting an analytical integral representation for the distance between scattering points, Hoshiba [1997] developed an efficient algorithm for a simple structure model and studied the characteristics of coda envelopes in layered media. In order to calculate seismogram envelopes for a realistic Earth model in lithosphere. Many methods for synthesizing coda envelopes which seismic velocity varies complexly we introduce a direct based on the energy transport approach have been proposed for this purpose [e.g., Gusev end Abubekirov, 1987, 1996; Hoshibe, 1993; Zeng et el., 1991; Nishimure et el., 1997; Yoshimoto et el., 1997a, b]. By incorporating the effect of multiple seismic scattering, the analysis of coda envelopes at large lapse simulation Monte Carlo (DSMC) method, which has been used in the field of the kinetic theory of gases [e.g., Bird, 1976, 1994; Bellomo, 1995]. The DSMC method, which utilizes a finite difference scheme for ray tracing, has an advantage over previous methods in terms of flexibility and can be applied to times made it possible to separately measure the scattering three-dimensional (3-D) (i.e., laterally varying) velocity struccoefficient and the intrinsic attenuation of the lithosphere [e.g., Meyede et el., 1992; Hoshibe, 1993; Jin et el., 1994]. However, tures. In this paper, we synthesize seismogram envelopes for a shallow local earthquake using a realistic velocity structure. there are few studies that adopt a realistic model of seismic Instead of the layered structures which were used by Hoshiba velocity structure [Hoshibe, 1997; Mergerin et el., 1998]. The use of an inappropriate model, such as a homogeneous velocity structure, for the envelope analysis may cause significant error [1997] and Margetin et al. [1998], we adopt a lithospheric model with a positive velocity gradient with depth to evaluate the effect of seismic ray bending on coda envelopes. in estimated medium parameters of the heterogeneous lithosphere (e.g., scattering coefficient). 2. Method In general, seismic velocity in the lithosphere increases with depth [e.g., Christensen end Moony, 1995]. In such a structure, We apply the DSMC method to synthesize seismogram envelopes in 3-D scattering media. The treatment of seismic not only ray paths of direct waves but also those of scattered waves is completely acoustic in this paper (i.e., only S waves waves bend upward. This causes concentration and stagnation are considered acoustically). For simplicity, let us suppose a Copyright 2000 by the American Geophysical Union. Paper number 1999JB /00/1999JB D medium in which the scattering coefficient # is spatially uniform and the seismic velocity changes only vertically. The scattering coefficient is also known as the turbidity and the 6153
2 6154 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE (a) (b) Source z So u rc Particle i x Receiver Scat,t/ering Free Surface Torus Volume x Particle time interval At the particle at depth z propagates over a distance of v(z)at, and the probability of the occurrence of scattering along this path is represented by 1 - exp[-v(z)at/ l] --- v(z)at/l, where v(z) is the seismic velocity at the depth and l -- #-1 is the mean free path of the scattering medium. This probability can be rewritten as g v(z)at by using the scattering coefficient. Under the condition of v(z)at << #- 1 the process of propagation is approximately characterized by successive independent mathematical operations that represent free propagation and scattering, respectively [Bird, 1967, 1994; Nanbu, 1980]. This is called the principle of uncoupling. In the case that the particle is scattered in this time interval its motion is calculated as the sequence that the particle first moves over the distance v(z)at without scattering and then changes its propagation direction by scattering. The occurrence of scattering along the distance of v(z)at is determined by whether the following inequality is satisfied: gv(z)at > U3, (2) where U3 is the uniform random deviate between 0.0 and 1.0. Since the probability distribution of the scattering direction is assumed to be isotropic, the propagation direction of the scattered particle is also evaluated by using (1); a newly generated set of parameters is used for the next step of calculation. Introducing the x'-z plane which includes ray path as shown in Figure la, the vector representing the propagation of the particle in time interval At is (Ax', Az) = [v(z)at sin O, v(z)at cos O], (3) where O represents the angle of particle propagation measured from the downward vertical and satisfies the following equation [Cerveny and Rayindra, 1971]' dv(z) A O = d-- -- At sin O. (4) Figure 1. Schematic illustration of the propagation of energy particle. The location of the source is taken as the origin of the Cartesian coordinate. (a) A spherical volume at a receiver to count the number of energy particles. (b) A torus volume just beneath the free surface to count the number of energy particles at a receiver at epicentral distance A. scattering cross section per unit volume [Dainty, 1981]. We assume that the source radiation and the scattering of seismic waves are isotropic, which implies a spherical symmetry of these phenomena. Random perturbation of seismic velocity in space is not considered here. We shoot an energy particle from a point source (Figure la). Because of the assumption of isotropy, takeoff angle 0 and azimuthal angle 0 of the particle are given randomly as arccos (1-2U1), 2'rrU2, (1) NO. I l SmaRT I ß I SHOO,N I '1 I,.,,, ua NO I YES (gvat> 1' J Sq.1 I where U and U 2 are the uniform random deviates between 0.0 and 1.0. Random variables U and U2 are generated independently. Equation (1) ensures uniform probability distribution of the shooting direction on the focal sphere. The energy particle, which represents a seismic wavelet of unit energy, propagates in the 3-D medium following seismic ray theory. The location of the particle at a certain lapse time is calculated on the basis of a finite difference scheme. For a small T YES Figure 2. Flow chart of the numerical calculation of seismogram envelope.
3 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE 6155 (a) lo 2o 3o 4o 5o Model 1" Uødel11' Model 2 I ' Model 1 I I I ' ' ' dia for every time step, unlike the methods of Hoshiba [1997] and Margerin et al. [1998]. This step-by-step calculation requires additional computational work. However, because of the simplicity of the algorithm the DSMC method can be easily applied to a structure with a complex spatial variation of velocity and is executable on a personal computer with small on-board memory. A computational investigation of the accuracy of numerical calculation by this method is presented in the appendix. 3. Numerical Simulations 3.1. Seismogram Envelopes in a Realistic Velocity Structure We synthesize seismogram envelopes in a scattering medium with a vertical velocity gradient by using a model of S wave velocity structure beneath the Tohoku district, Japan [Hase- $ wave Velocity (km/s) (b) Epicentral Distance (km) 10'7 _o 40.o loo 1.o u 1/ ' -... / -. ß '- G 5o Figure 3. (a) Veloci structure models of S waves used in the synthesis. (b) Schematic illustration of the free propagation of energy particle. Model 1 is used for the calculation of ray path. =50km The left-hand term represents the increment of 0 in time interval At. The initial value of 0 is set to be 0 at the location of the source. The location of the energy particle at a certain lapse time t is evaluated by iterating the calculation described above. In this computation the highest order of scattering possible is the total number of iterations because a scattering is allowed only once in each time interval At. Under the condition of #v(z)at << 1 a small time step gives a good approximation for the tracing of the particle in a scattering media. A small volume element A V is considered at a receiver. In the case where the particle is included in this volume element, we count it as an arrival of a seismic energy packet. On the basis of the Monte Carlo scheme we shoot many energy particles of number N from the source. The energy density at the receiver at a time t is evaluated as n(t)/(nav), where n(t) is the number of particles included in the volume element at the receiver. We refer to the time history of the seismic energy density at a receiver as the seismogram envelope throughout this paper. The size and shape of the volume element affect the result; however, this effect is not severe for the coda portion where the spatial and temporal variation of energy density is small. Time-space distribution of energy density is also evaluated from the same calculation for a set of spatially distributed receivers. Figure 2 shows the flow chart of the seismogram envelope synthesis. The DSMC traces particle motion in scattering me- A=lOOkm A= 150km A=200km \ 35km/s Figure 4. An example of the synthetic seismogram envelopes for model 1 in Figure 3a. Envelopes of the surface receivers at epicentral distances from 10 to 200 km are exhibited. Sampling interval of the data is 1 s. The lapse time is measured from the source origin time. The dashed line indicates apparent velocity of 4.35 km s -, which is the S wave velocity of the uppermost mantle. The source depth is 15 km. The scattering coefficient is 0.01 km -. The parameters used in the simulation are as follows: At = 0.2 s, N = 107 and AV = 1.0 x x 103 km 3 (varying with epicentral distance as A V c A).
4 , 6156 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE. '10-5 I,,,,,,,,I,,, 10-6, 10 '7 LU ß 10.8, I 10 ' 2 = 50 kn "'"" >' I 10-6 L 10.7 ß 10-8 E z lo ,J II... : '5 r 10 ' LU ß 10.8 E ' '6 10 ' A=I I i i i i i i1[ i i i i i i i i i i ij i i i A =, Figure 5. Seismogram envelopes for different velocity structures. The solid, dashed and dotted lines indicate the envelopes for the velocity models of model 1, 2, and 3 in Figure 3a, respectively. Envelopes at epicentral distances of 50, 100, 150, and 200 km are displayed on Figures 5a, 5b, 5c, and 5d, respectively. Sampling interval of the data is 1 s. The source depth is 15 km. The scattering coefficient is 0.01 km -. The parameters used in the simulation are as follows: At = 0.2 s, N = 107, and AV = 5.0 x 102, 1.0 x 103, 1.5 x 103, and 2.0 x 103 km 3 (at epicentral distances 50, 100, 150, and 200 km, for lapse times smaller than 200 s). For lapse times larger than 200 s, a receiver volume 4 times as large as the original is used to minimize the ripple of envelope. gawa et al., 1978]. The model, which is used for routine hypo- on the envelopes is a result of the finite number of energy central determination Tohoku University, has an S wave particles. It takes 4 hours of CPU time on a 270-MHz Sun velocity of 3.15 km s - at the free surface and 4.35 km s - in Ultra5 to calculate the 20 envelopes on Figure 4, where 20% of the uppermost mantle, as shown in Figure 3a. The velocity has the time is consumed by generating random numbers. The maximum gradient at the crust-mantle transition at a depth of apparent velocity of the first arrival increases with epicentral -30 km. For simplicity, a reflection from the transition struc- distance, as is expected for the direct wave. The shape of the ture which includes a discontinuity of the first derivative of envelope is impulsive at small epicentral distances broadening velocity is not modeled in the synthesis. We assume isotropic as the distance increases. At epicentral distances smaller than scattering and a constant scattering coefficient of 0.01 km km the first part of the envelope corresponds to the direct irrespective of depth. Intrinsic attenuation is not considered arrival of S waves, which propagate through the crust. In conhere. We put a point source of isotropic radiation at a depth of trast, the shape of the first arrival at larger distances affected 15 km so that the effect of velocity structure on seismogram by the velocity structure of the crust-mantle transition, and the envelopes is emphasized. Examples of free (nonscattering) effect of surface reflection dominates. The dominant impulsive propagation of energy particles are schematically shown in signal seen at epicentral distances larger than 150 km reflects Figure 3b. For the free surface we assume total reflection the arrival of the ss phase. according to ray theory. By multiple reflections at the free At large lapse time the decay of envelope amplitude is surface an energy particle that is shot in a nearly horizontal smooth irrespective of epicentral distance. This indicates that direction from the source propagates through the crust, that is, the characteristic structure of the crust-mantle transition does a crust with a positive velocity gradient behaves as a waveguide. not cause any short-time fluctuations of envelope amplitude In the numerical calculations of particle tracing, we adopt a except around the direct arrivals. We note that the small decay time step of At = 0.2 s, which is roughly one hundredth of the rate of amplitude in the coda is due to the exclusion of the mean free time l/v(z) and is small enough to ensure the intrinsic attenuation in the synthesis. The intrinsic attenuation stability of the calculation (see Figure A2). By usin geomet- can easily be incorporated into the DSMC method in a rical symmetry with respect to the epicenter we used a torus straightforward manner. volume just beneath the free surface to count the number of energy particles at a receiver (Figure lb). The size of the torus, 3.2. Effect of a Velocity Gradient on which increases with the epicentral distance, chosen to be Seismogram Envelopes large enough to stabilize ripple of the envelopes but small Figure 5 is a comparative plot of envelopes for the three enough to retain temporal resolution. velocity models: models 1, 2, and 3 in Figure 3a. Models 2 and Figure 4 exhibitsynthesized envelopes for surface receivers 3 are homogeneous velocity models and take 3.15 and 4.35 km at epicentral distances from 10 to 200 km. The amplitude s -, which correspond to the velocity at the free surface and sampled every 1 s is shown on a logarithmic scale. The ripple the uppermost mantle of model 1, respectively. Envelopes at
5 . YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE 6157 the epicentral distances of 50, 100, 150, and 200 km are calculated with the same parameters as those of Figure 4. A nonimpulsive shape of the seismogram envelope around the direct arrival is most prominent for model 1 with the velocity gradient. This is due to the free surface reflection and the ray bending caused by the velocity structure of the crustmantle transition mentioned in section 3.1. For each velocity model, envelope amplitudes at different epicentral distances are almost identical at large lapse times (see, for example, amplitudes at a lapse time of 500 s). As seen in Figures 5a-5d, the envelope shows a systematic difference in amplitude between the homogeneous velocity models (models 2 and 3). The amplitude at a fixed lapse time is consistently smaller for the high-velocity model. This is explained by the fast spatial spreading of the seismic energy in the high-velocity model. The normalization of the lapse time by using the seismic velocities of these models substantially diminishes this difference (see equation (A1) in the appendix). From just after the onset of envelope, as seen in Figure 5a, the exponent of the decay rate of amplitude is almost equal to t-.s, which is predicted by the long-time limit of the diffusion solution for the infinite medium. 10 '4 10 's 10 ' s ' '9 _ A :50km (a) 100km Model 1 L. J5Okm... I... I... I... I... I... I... I... I I... I... I... I... I... I... I... h,,,,,,,,h,,,,,,,, (b) Model 1' [... I'"' '1... I... I... I... I... I... I '4 I... I... I... I... I... I... I... I... I... I... I As depicted by the seismogram envelope of model 1 in Figure 5, the positive velocity gradient with depth increases the Model 1" amplitude of the envelope. This is the same tendency reported 10-6 by Hoshiba [1997] and Margerin et al. [1998] for layered structures. The amplitude calculated for the velocity gradient model 10'7 is roughly twice that for the homogeneous velocity models at 10-8 large lapse times. The temporal decay rate of amplitude is not identical with those of the homogeneous velocity models; how- 10 ' I ever, it is roughly characterized by t-.5 as the lapse time increases. Our result clearly indicates that upward ray bending due to the positive velocity gradient as well as the multiple reflection in the layer structure causes amplification of the seismogram envelope. We investigate how a small variation of a velocity structure Figure 6. Seismogram envelopes for different velocity structures. Envelopes at epicentral distances of 50, 100, 150, and affects the shape of seismogram envelopes. For this purpose, 200 km are shown. Sampling interval of the data is 1 s. Figures we use two models: one is a model with a low-velocity zone in 6a, 6b, and 6c are calculated using the velocity models of the lower crust (model 1') and the other has a high positive models 1, 1', and 1" in Figure 3a, respectively. The source velocity gradient near the free surface (model 1"). These mod- depth is 15 km. Set of parameters, except for At = 0.05 s, used els are constructed from model 1. Figure 6 shows seismogram in the calculation is identical with that used in Figure 4. envelopes for these three velocity models. The parameters, except for At = 0.05 s, used in the synthesis are identical with those of Figure 4. Except for the onset of envelopes, the shape of envelopes is almost identical among these models. The difference of amplitudes at a lapse time of 200 s is <10%. This implies that a small difference of velocity structure is easily hidden from the shape of coda envelopes. temporal rate of the energy leak decreases with lapse time; however, it is not identical among the velocity models. The smallest rate is found for the model with positive velocity gradient (model 1). In this model, >10% of the energy radiated from the source is seen at depths shallower than 15 km even at a lapse time of 200 s, twice the value for the homoge Spatial Energy Distribution neous velocity models. We call this phenomenon "stagnation." Figure 7, showing the depth distribution of the energy versus the lapse time, gives an intuitive insight into the dependence of the seismogram envelopes on velocity structure. The ratio of the energy included in the three layers bounded by depths 0-15, and 30-o km is plotted on the vertical axis (i.e., horizontal distribution of the energy is not reflected in Figure 7). For example, at a lapse time of 100 s, Figure 7a depicts that 21%, 13%, and 66% of energy is included in the top, middle, and bottom layers, respectively. As the lapse time increases, the seismic energy gradually leaks out of the crust for all the velocity models (Figures 7a, 7b, and 7c). The leak starts at a lapse time of -4 s when the wave front of the direct waves reaches the Moho. In general, the The dependence of energy stagnation on source depth is found in Figures 7a, 7d, and 7e. Depth of source is 15, 40, and 80 km for Figures 7a, 7d, and 7e, respectively. Since the source is located in the mantle for Figures 7d and 7e, the energy in the top and middle layer increases just after the arrival of the wave front at the lower boundary of each layer. As lapse time increases, energy propagation in the opposite direction due to the reflection of the free surface decreases the energy in these layers. The distribution of energy at short lapse times is quite different for each case; however, it is worthy to compare them at lapse times larger than several tens of seconds. The amount of energy stagnating in the crust, the top layer, and middle layer is clearly large for the shallow source: e.g., at a lapse time
6 6158 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE '?.!:!::::...:;:::.i::..:.:?....: i5.. i!?i:..: :::.i..:!????.! o ( ~ O Figure 7. Depth distribution of seismic energy versus lapse time. The ratio of seismic energy included in three different depth ranges (0-15, 15-30, 30-oo km) is shown. Figures 7a, 7b, and 7c are calculated using the velocity models of models 1, 2, and 3 in Figure 3a, respectively. Figures 7d and 7e depict the depth distribution of energy in model 1 for the sources with depths of 40 and 80 km, respectively. Set of parameters used in the calculation is identical with those used in Figure 5. of 80 s the energy of the event with a depth of 15 km is twice as large as that of 80 km. This indicates that the positive velocity gradient with depth effectively behaves as a waveguide for a shallow source Analysis of Seismogram Envelopes Our synthesis demonstrates that the effect of a positive velocity gradient with depth on the envelope amplitude is not negligible; a significant enlargement of envelope amplitude due to the stagnation of energy at shallow depths occurs for a scattering medium. It is not easy to discriminate between this phenomenon and that due to strong scattering. This implies that the analysis of observed seismogram envelopes of local and regional earthquakes might result in an overestimate of the scattering coefficient when a homogeneous velocity model is assumed. To investigate this, we carry out a numerical investigation of estimation of the scattering parameter by inversion analysis of seismogram envelopes using the theory of this paper. The first 20 s of the seismogram envelopes at epicentral distances of 10, 30, 50, 70, and 90 km in Figure 4 are used as the data of the inversion. We adopt a linearized iterative technique which is based on a singular value decomposition method. Assuming a homogeneous velocity model with an S wave velocity of 3.5 km s -, we seek a model parameter # which minimizes the squared residual between the seismogram envelopes of the data and this model. The residual of each sampling point is normalized by the amplitude of the envelope. A finite differ- ence approximation is used to calculate the differentials of seismogram envelopes by perturbing the # value 5%. Irrespective of an initial value of #, the squared residual becomes sufficiently small after several iterations. The estimated value is not sensitive to the initial value and is almost 60% larger than # km-, which is used to calculate the data envelopes. We apply the multiple lapse time window analysis [Hoshiba, 1993] to the envelopes at epicentral distances between 10 and 100 km in Figure 4. This method, which employs time windows for the analysis of seismogram envelopes, is commonly used in investigating the scattering coefficient # on the basis of the analysis of the time-space distribution of coda wave energy. In this numerical test, following Hoshiba [1993], we use three time windows of 15 s, which are applied from the onset of S waves without overlap. Let us representhe integral of seismic energy in each time window with the normalization of the envelope amplitude at lapse time of 60 s as E (r), E2(r), and E3(r), where r is the hypocentral distance. Adopting the homogenous velocity model used previously, we carried out the multiple lapse time window analysis to seek the # value which minimizes the squared log residual between data and calculation. In the grid search analysis, # value is varied from 0.0 to 2.0 x 10-2 km - by a step of 5.0 x 10-3 km -. The best fit result is shown in Figure 8. The residual takes a minimum value, which is one half of that for the true value, at# = 1.35 x 10-2 km -. Figure 8 shows good coincidence between data and calculation. However, our result indicates that the multiple lapse time
7 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE r r 2El(r) the field measurement of scattering coefficient under the assumption of homogeneous velocity structure might be overestimated as a result. Our result indicates that an appropriate structure model is indispensable for the synthesis of realistic seismogram envelopes r To evaluate the DSMC method, we investigate the accuracy of the numerical calculation by comparing it with an analytic r r2e3(r) solution, and varying the calculation parameters. Since there is at present no analytic solution, or Green function, of the multiple isotropic scattering model for a 3-D medium in which g best-' 1.35x 10-2 km'l velocity change spatially, we use the analytic solution obtained '0 4'0 6'0 ' 8'0 ' 1 )0 by Zeng et al. [1991] for a homogeneous velocity structure in an infinite volume. We synthesize seismogram envelopes for the 3-D infinite scattering medium characterized by the parameters in section 2 and compare them with those of analytic Hypocentral Distance (km) solutions. Source radiation is assumed to be temporally impul- Figure 8. Plots of normalized energy corrected for geo- sive and azimuthally isotropic. For the syntheses in this section metrical spreading by 4 rr 2. Circles, triangles, and diamonds represent energy integrals for the time windows of 0-15, 15-30, and s after the onset of S waves. The solid line shows the best fit result, which is calculated using # = 1.35 x 10-2 ' a) 0.5 >, 0.4 p =0.8 window analysis may overestimate # value under an inappropriate (homogeneous) structure model. The use of a precise velocity model is essential for envelope analysis. We believe that the analysis of seismogram envelopes is important for investigation of the global and regional structure of the lithosphere in connection with the characteristics of short-wavelengtheterogeneities. The DSMC method, which is very flexible in incorporating the effects of free-surface reflection, intrinsic attenuation, nonisotropic source radiation, and nonisotropic scattering on the seismogram envelope, would be a very effective tool for this study. 4. Conclusions A numerical algorithm called the DSMC method is introduced to synthesize seismogram envelopes in a scattering medium with nonuniform seismic velocity. We checked the capability of this method by comparing its numerical outputs and analytic solutions for a multiple isotropic scattering model (see the appendix). The result of the inspection indicates high precision and numerical stability of the DSMC method. Because of the simplicity of the algorithm it is easy to incorporate flee-surface reflection, intrinsic attenuation, nonisotropic source radiation, and nonisotropic scattering as well as the depth dependency of seismic wave velocity. The calculation can be executed on a personal computer with small on-board memory. Using the DSMC method, we investigated the characteristics of the seismogram envelope in a structure with a positive velocity gradient with depth. The velocity gradient alters not only envelope amplitude but also its temporal shape. The en- Appendix: Computational Accuracy and Stability o z ß 0.07 ß 0.06 >, 0.05 o) ' Normalized Time! i i i i i i i p u.i 0.02 O.Ol o Z ' O. 002 E O Normalized Time Normalized Time Figure A1. Normalized energy density at normalized distances of p = 0.8, 1.6, and 3.2. The result from the DSMC method (the solid line) and that from the analytic calculation velope shape of a shallow earthquake is affected by the flee- (the dashed line) are compared. The former is the average of surface reflection and by the ray bending due to the velocity 20 independent numerical simulations. Standard deviation is change at the crust-mantle transition. The stagnation of energy represented by the dotted lines. The solid and dotted lines lie at shallow depths is clearly found in the synthesis. Since this on top of the dashed line except for p = 3.2. The parameters used phenomenon is not easily discriminated from strong scattering, in the simulation are A, = 0.01, N = 106, and AV =
8 6160 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE o.o?, o.o?j I >,0.05' 0.06' I >,0.05J 0.06'J [ I A ' = ' O.04-J 0.03,. c e 0.03l 0.02l W o.02j [ O.OlJ O.01J 11 o o.oo I O0'00 'J Normalized Time Normalized Time c 0.07' ( 0.06 >,0.05 b z = e c 0.03 iil 0.02 E O.01 o 0.00 o 3 6 f 8 Normalized Time ' >,0.05- o.o4l ce o.o3l E o.oo Normalized Time Figure A2. Numerical investigation of the dependence of the result of the simulation on AT. Normalized energy density is calculated for p = 1.6. The result from the DSMC method (the solid line) and that from the analytic calculation (the dashed line) are compared. The former is the average of 20 independent numerical simulations. Standardeviation is represented by the dotted line. The solid and dotted lines lie on top of the dashed line for AT and 0.1. The set of parameters used in the calculation is identical with those used in Figure A1. we use the normalized time T and the normalized distance p following Sato [1995]: T = gvt, p = gr. (A1) The parameter r is the distance between the source and the receiver. the coarse sampling interval and is not critical. The shapes of the envelopes are identical until Av exceeds 0.1. As the value of Av increases beyond 0.1, the DSMC method exhibits slightly larger outputs compared with the reference. However, the DSMC method does not show any numerical instability. This result indicates that the upper limit of the time step interval in the calculation for accurate results is about one-tenth of the mean free time of the scattering medium. Figure A1 shows seismogram envelopes at different normalized distances: p - 0.8, 1.6, and 3.2. The average and plus or minus one standard deviation for 20 simulations are repre- Acknowledgments. The author is grateful to Haruo Sato for his sented by the solid line and the dotted line, respectively. The valuable discussions and suggestive comments. The author thanks number of particles used in each simulation is 106. A spherical Masakazu Ohtake for his thoughtful reviews and suggestions for volume with the radius of 0.1 is used to estimate the spatial greatly improving the manuscript. Thanks also to Shigeo Kinoshita for his kind encouragement for this study. Critical and helpful reviews of density of energy particles at the receiver. The amplitude of this manuscript by Anton Dainty, Stephen Gao, and Mitsuyuki seismogram envelopes calculated by the DSMC method is Hoshiba are gratefully acknowledged. For this study the author used identical with that calculated by the method of Zeng et al. the computer system of the Earthquake Research Institute, University [1991](the dashed line), except for ripples of the envelope at of Tokyo. p = 3.2. The ripples scatter around the analytic solution; however, they do not show any apparent bias within plus or minus one standard deviation. Shooting a larger number of particles References substantially diminishes the ripples. In addition, the geometric Aki, K., Scattering and attenuation of shear waves in the lithosphere, symmetry with respecto the point source may also be used for J. Geophys. Res., 85, , this purpose: the use of a concentric spherical shell with radius Aki, K., and B. Chouet, Origin of coda waves: Source, attenuation and scattering effects, J. Geophys. Res., 80, , p = 3.2 as the receiver volume element effectively reduces the Bellomo, N. (Ed.) Lecture Notes on the Mathematical Theory of the ripples. Boltzmann Equation, World Sci., River Edge, N.J., The DSMC method employs the decoupling principle, ex- Bird, G. A., The velocity distribution function within a shock wave, J. pressed by the inequality AT << 1, to calculate the motion of Fluid Mech., 30, , Bird, G. A., Molecular Gas Dynamics, Oxford Univ. Press, New York, energy particles. Next, we investigate the dependence of syn thesized envelopes on AT. Figure A2 displays the seismogram Bird, G. A., Molecular Gas Dynamics and the Direct Simulation of Gas envelopes at p = 1.6 for the numerical conditions AT = 0.01, Flows, Oxford Univ. Press, New York, , 0.2, and 0.4. After synthesizing 20 envelopes for each Cerveny, V., and R. Ravindra, Theory of Seismic Head Waves, Univ. of value, their average (the solid line) and standard deviation (the Toronto Press, Toronto, Ont., Christensen, N. I., and W. D. Mooney, Seismic velocity structure and dotted line) are calculated. One million particles are used in composition of the continental crust: A global view, J. Geophys. Res., the calculation. The analytic solution (the dashed line) is plot- 100, , ted for comparison. The difference around the onset is due to Dainty, A.M., A scattering model to explain seismic Q observations in
9 YOSHIMOTO: MONTE CARLO SIMULATION OF SEISMOGRAM ENVELOPE 6161 the lithosphere between 1 and 30 Hz, Geophys. Res. Lett., 8, equation, I, Monocomponent gases, J. Phys. Soc. Jpn., 49, , 1128, Gusev, A. A., and I. R. Abubakirov, Monte-Carlo simulation of record Nishimura, T., M. Fehler, W. S. Baldridge, P. Roberts, and L. Steck, envelope of a near earthquake, Phys. Earth Planet. Inter., 49, 30-36, Heterogeneoustructure around the Jemez volcanic field, New Mex ico, USA, as inferred from the envelope inversion of active- Gusev, A. A., and I. R. Abubakirov, Simulated envelopes of non- experiment seismic data, Geophys. J. Int., 131, , isotropically scattered body waves as compared to observed ones: Sato, H., Formulation of the multiple non-isotropic scattering process Another manifestation of fractal heterogeneity, Geophys. J. Int., 127, in 3-D space on the basis of energy transportheory, Geophys. J. Int., 49-60, , , Hasegawa, A., N. Umino, and A. Takagi, Double-planed structure of Vernon, F. L., G. L. Pavlis, T. J. Owens, D. E. McNamara, and P. N. the deep seismic zone in the northeastern Japan arc, Tectonophysics, Anderson, Near-surface scattering effects observed with a high- 47, 43-58, frequency phased array at Pinyon Flats, California, Bull. Seismol. Hayakawa, T., Derivation of the envelope Green function in scattering Soc. Am., 88, , media with the depth-dependent velocity structure, and its applica- Yoshimoto, K., H. Sato, and M. Ohtake, Three-component seismotion to the analysis of source rupture process (in Japanese), master's gram envelope synthesis in randomly inhomogeneous semi-infinite thesis, 94 pp., Tohoku Univ., Sendai, Japan, media based on the single scattering approximation, Phys. Earth Hoshiba, M., Separation of scattering attenuation and intrinsic absorp- Planet. Inter., 104, 37-61, 1997a. tion in Japan using the multiple lapse time window analysis of full Yoshimoto, K., H. Sato, and M. Ohtake, Short-wavelength crustal seismogram envelope, J. Geophys. Res., 98, 15,809-15,824, heterogeneities in the Nikko area, central Japan, revealed from the Hoshiba, M., Seismic coda wave envelope in depth-dependent S wave three-component seismogram envelope analysis, Phys. Earth Planet. velocity structure, Phys. Earth Planet. Inter., 104, 15-22, Inter., 104, 63-73, 1997b. Jin, A., K. Mayeda, D. Adams, and K. Aki, Separation of intrinsic and Zeng, Y., F. Su, and K. Aki, Scattering wave energy propagation in a scattering attenuation in southern California using TERRAscope random isotropic scattering medium, 1, Theory, J. Geophys. Res., 96, data, J. Geophys. Res., 99, 17,835-17,848, , Margerin, L., M. Campillo, and B. V. Tiggelen, Radiative transfer and diffusion of waves in a layered medium: New insight into coda Q, K. Yoshimoto, Department of Geophysics, Graduate School of Sci- Geophys. J. Int., 134, , ence, Tohoku University, Aoba-ku, Sendai , Japan. Mayeda, K., S. Koyanagi, M. Hoshiba, K. Aki, and Y. Zeng, A com- tohoku.ac.jp) parative study of scattering, intrinsic, and coda Q- for Hawaii, Long Valley, and central California between 1.5 and 15.0 Hz, J. Geophys. Res., 97, , (Received January 25, 1999; revised October 19, 1999; Nanbu, K., Direct simulation scheme derived from the Boltzmann accepted December 6, 1999.)
Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationIntrinsic and Scattering Seismic Attenuation in W. Greece
Pure appl. geophys. 153 (1998) 703 712 0033 4553/98/040703 10 $ 1.50+0.20/0 Intrinsic and Scattering Seismic Attenuation in W. Greece G-AKIS TSELENTIS 1 Abstract Intrinsic (Q 1 i ) and scattering (Q 1
More informationSEISMIC WAVE PROPAGATION AND SCATTERING IN THE HETEROGENEOUS EARTH
Haruo Sato Tohoku University, Japan Michael C. Fehler Massachusetts Institute of Technology, U.S.A. Takuto Maeda The University of Tokyo, Japan SEISMIC WAVE PROPAGATION AND SCATTERING IN THE HETEROGENEOUS
More informationScattering and intrinsic attenuation structure in Central Anatolia, Turkey using BRTR (PS-43) array
Scattering and intrinsic attenuation structure in Central Anatolia, Turkey using BRTR (PS-43) array CTBT: Science & Technology 2011 Korhan Umut SEMIN Nurcan Meral OZEL B.U. Kandilli Observatory & Earthquake
More informationUnified Explanation of Envelope Broadening and Maximum-Amplitude. Decay of High-Frequency Seismograms based on the Envelope
submitted to J. Geophys. Res. Unified Explanation of Envelope Broadening and Maximum-Amplitude Decay of High-Frequency Seismograms based on the Envelope Simulation using the Markov Approximation: Forearc
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2010) 182, 988 1000 doi: 10.1111/j.1365-246X.2010.04657.x Envelope broadening characteristics of crustal earthquakes in northeastern Honshu, Japan Jayant
More informationSUMMARY. fractured reservoirs in the context of exploration seismology.
Modeling scattering of cross-well seismic waves using Radiative Transfer Theory Josimar A. Da Silva Jr, Oleg V. Poliannikov and Michael Fehler, Earth Resources Laboratory / M.I.T SUMMARY We model P and
More informationMODELING OF HIGH-FREQUENCY WAVE RADIATION PROCESS ON THE FAULT PLANE FROM THE ENVELOPE FITTING OF ACCELERATION RECORDS
MODELING OF HIGH-FREQUENCY WAVE RADIATION PROCESS ON THE FAULT PLANE FROM THE ENVELOPE FITTING OF ACCELERATION RECORDS Yasumaro KAKEHI 1 SUMMARY High-frequency (higher than 1 Hz) wave radiation processes
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
DEVELOPMENT OF REGIONAL PHASE TOMOGRAPHIC ATTENUATION MODELS FOR EURASIA Thorne Lay 1, Xiao-Bi Xie 1, and Xiaoning (David) Yang 2 University of California, Santa Cruz 1 and Los Alamos National Laboratory
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
ANALYSIS AND SIMULATION OF THREE-DIMENSIONAL SCATTERING DUE TO HETEROGENEOUS CRUSTAL STRUCTURE AND SURFACE TOPOGRAPHY ON REGIONAL PHASES; MAGNITUDE AND DISCRIMINATION Arben Pitarka 1, Don V. Helmberger
More informationMain Menu. Summary. Introduction
Kyosuke Okamoto *, JSPS Research Fellow, Kyoto University; Ru-shan Wu, University of California, Santa Cruz; Hitoshi Mikada, Tada-nori Goto, Junichi Takekawa, Kyoto University Summary Coda-Q is a stochastic
More information2009 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
DEVELOPMENT OF REGIONAL PHASE TOMOGRAPHIC ATTENUATION MODELS FOR EURASIA Thorne Lay 1, Xiao-Bi Xie 1, and Xiaoning (David) Yang 2 University of California, Santa Cruz 1 and Los Alamos National Laboratory
More informationHaruhisa N. (Fig. + ) *+ Graduate School of Environmental Studies, Nagoya University, Furo-cho, Chikusa-ku, Nagoya.0. 20*+ Japan.
/- (,**2) 0,+/,,+ Source Mechanism and Seismic Velocity Structure of Source Region of Deep Low-frequency Earthquakes beneath Volcanoes: Case Studies of Mt Iwate and Mt Fuji Haruhisa N AKAMICHI + +3 (Fig
More informationSimulating the Envelope of Scalar Waves in 2D Random Media Having Power-Law Spectra of Velocity Fluctuation
Bulletin of the Seismological Society of America, Vol. 93, No., pp. 4 5, February 3 Simulating the Envelope of Scalar Waves in D Random Media Having Power-Law Spectra of Velocity Fluctuation by Tatsuhiko
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationSeismic Scattering in the Deep Earth
Seismic Scattering in the Deep Earth Peter Shearer IGPP/SIO/U.C. San Diego September 2, 2009 Earthquake Research Institute Mantle mixing calculations Heterogeneity is likely at all scales Davies (2002)
More informationOrigin of Coda Waves: Earthquake Source Resonance
Origin of Coda Waves: Earthquake Source Resonance Yinbin Liu Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia, Vancouver, BC, Canada Email: yliu@eoas.ubc.ca Abstract
More informationSeismic Wave Propagation and Scattering in the Heterogeneous Earth : Second Edition
Seismic Wave Propagation and Scattering in the Heterogeneous Earth : Second Edition Bearbeitet von Haruo Sato, Michael C. Fehler, Takuto Maeda 1. Auflage 2012. Buch. xvi, 496 S. Hardcover ISBN 978 3 642
More informationTomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima
1 2 3 Auxiliary materials for Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area 4 5 6 7 8 9 Ping Tong 1,2, Dapeng Zhao 1 and Dinghui Yang 2 [1] {Department of Geophysics,
More informationThe global short-period wavefield modelled with a Monte Carlo seismic phonon method
Geophys. J. Int. (24) 158, 113 1117 doi: 1.1111/j.1365-246X.24.2378.x The global short-period wavefield modelled with a Monte Carlo seismic phonon method Peter M. Shearer 1 and Paul S. Earle 2 1 Institute
More informationImaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges
2017 SCEC Proposal Report #17133 Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges Principal Investigator Zhongwen Zhan Seismological Laboratory, California
More informationEstimation of Coda Wave Attenuation Quality Factor from Digital Seismogram Using Statistical Approach
Science and Technology 2012, 2(1): 1-7 DOI: 10.5923/j.scit.20120201.01 Estimation of Coda Wave Attenuation Quality Factor from Digital Seismogram Using Statistical Approach Jwngsar Brahma School of Petroleum
More informationTeleseismic receiver function using stacking and smoothing of multi seismic-records at a single station
Earthq Sci (2012)25: 75 81 75 doi:10.1007/s11589-012-0833-7 Teleseismic receiver function using stacking and smoothing of multi seismic-records at a single station Yi Yang and Fuhu Xie Earthquake Administration
More informationANALYTICAL STUDY ON RELIABILITY OF SEISMIC SITE-SPECIFIC CHARACTERISTICS ESTIMATED FROM MICROTREMOR MEASUREMENTS
ANALYTICAL STUDY ON RELIABILITY OF SEISMIC SITE-SPECIFIC CHARACTERISTICS ESTIMATED FROM MICROTREMOR MEASUREMENTS Boming ZHAO 1, Masanori HORIKE 2 And Yoshihiro TAKEUCHI 3 SUMMARY We have examined the site
More informationSimulation of earthquake rupture process and strong ground motion
Simulation of earthquake rupture process and strong ground motion Takashi Miyatake (1) and Tomohiro Inoue (2) (1) Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo, Tokyo, 113-0032, Japan
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationTHREE-DIMENSIONAL FINITE DIFFERENCE SIMULATION OF LONG-PERIOD GROUND MOTION IN THE KANTO PLAIN, JAPAN
THREE-DIMENSIONAL FINITE DIFFERENCE SIMULATION OF LONG-PERIOD GROUND MOTION IN THE KANTO PLAIN, JAPAN Nobuyuki YAMADA 1 And Hiroaki YAMANAKA 2 SUMMARY This study tried to simulate the long-period earthquake
More informationBROADBAND STRONG MOTION SIMULATION OF THE 2004 NIIGATA- KEN CHUETSU EARTHQUAKE: SOURCE AND SITE EFFECTS
Third International Symposium on the Effects of Surface Geology on Seismic Motion Grenoble, France, 30 August - 1 September 2006 Paper Number: 105 BROADBAND STRONG MOTION SIMULATION OF THE 2004 NIIGATA-
More informationEvidence for P P asymmetrical scattering at near podal distances
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl052179, 2012 Evidence for P P asymmetrical scattering at near podal distances Wenbo Wu, 1,2 Sidao Ni, 2 and Xiangfang Zeng 1 Received 3 May 2012;
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
MODELING TRAVEL-TIME CORRELATIONS BASED ON SENSITIVITY KERNELS AND CORRELATED VELOCITY ANOMALIES William L. Rodi 1 and Stephen C. Myers 2 Massachusetts Institute of Technology 1 and Lawrence Livermore
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More informationby Xiao-Bi Xie and Thorne Lay
Bulletin of the Seismological Society of America, Vol. 17, No. 1, pp. 22 33, February 217, doi: 1.1785/121623 E Effects of Laterally Varying Mantle Lid Velocity Gradient and Crustal Thickness on Pn Geometric
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/1131692/dc1 Supporting Online Material for Localized Temporal Change of the Earth s Inner Core Boundary This PDF file includes: Materials and Methods Figs. S1 to S3
More informationTracing rays through the Earth
Tracing rays through the Earth Ray parameter p: source receiv er i 1 V 1 sin i 1 = sin i 2 = = sin i n = const. = p V 1 V 2 V n p is constant for a given ray i 2 i 3 i 4 V 2 V 3 V 4 i critical If V increases
More informationReceived May 11, 1989; Accepted July 10, * To whom correspondence should be addresseḍ
J. Phys. Earth, 38, 163-177, 1990 Teleseismic P- Wave Travel Time and Amplitude Anomalies Observed in Hokkaido Region, Japan Ichiro Nakanishi1,* and Yoshinobu Motoya2 1Department of Geophysics and 2Research
More informationEARTH STRUCTURE & DYNAMICS EARTHQUAKE SEISMOLOGY PRACTICALS. G.R. Foulger
1 EARTH STRUCTURE & DYNAMICS EARTHQUAKE SEISMOLOGY PRACTICALS G.R. Foulger 1. A large earthquake is recorded well at a three-component seismic station in Hawaii (coordinates 205 E, 20 N). The epicentral
More informationSeismic Scattering in the Deep Earth
Seismic Scattering in the Deep Earth Peter M. Shearer Institute of Geophysics and Planetary Physics Scripps Institution of Oceanography University of California, San Diego La Jolla, CA 92093-0225 October
More informationBasic Ray Tracing. Rick Aster and Sue Bilek. October 3, 2003
Basic Ray Tracing Rick Aster and Sue Bilek October 3, 3 A key observation that we can make about a seismic signal is its arrival time. From systematic observations of arrival times, we can deduce useful
More informationComparison of crustal and upper mantle heterogeneity in different time periods: Indonesian subduction zone to northern Australia
Earthq Sci (2014) 27(1):47 55 DOI 10.1007/s11589-013-0059-3 RESEARCH PAPER Comparison of crustal and upper mantle heterogeneity in different time periods: Indonesian subduction zone to northern Australia
More informationDOUBLE-PLANED STRUCTURE OF INTERMEDIATE- DEPTH SEISMIC ZONE AND THERMAL STRESS IN THE DESCENDING PLATE. (Received December 20, 1983)
J. Phys. Earth, 31, 329-347, 1983 DOUBLE-PLANED STRUCTURE OF INTERMEDIATE- DEPTH SEISMIC ZONE AND THERMAL STRESS IN THE DESCENDING PLATE Hiroyuki HAMAGUCHI, Kazuhiko GOTO, and Ziro SUZUKI Geophysical Institute,
More informationLETTER Earth Planets Space, 56, , 2004
LETTER Earth Planets Space, 56, 921 925, 2004 Inversion of the high-frequency source radiation of M6.8 Avachinsky Gulf, Kamchatka, earthquake using empirical and theoretical envelope Green functions Anatoly
More informationSome aspects of seismic tomography
Some aspects of seismic tomography Peter Shearer IGPP/SIO/U.C. San Diego September 7, 2009 Earthquake Research Institute Part 1: Global Tomography P velocity perturbations 200 km 1000 km 2700 km MIT 2006
More informationGlobal surface-wave tomography
Global surface-wave tomography Lapo Boschi (lapo@erdw.ethz.ch) October 7, 2009 Love and Rayleigh waves, radial anisotropy Whenever an elastic medium is bounded by a free surface, coherent waves arise that
More informationContents of this file
Geophysical Research Letters Supporting Information for Intraplate volcanism controlled by back-arc and continental structures in NE Asia inferred from trans-dimensional ambient noise tomography Seongryong
More information4((F'~) 2) = ~ = (2)
Bulletin of the Seismological Society of America, Vol. 74, No. 5, pp. 1615-1621, October 1984 AVERAGE BODY-WAVE RADIATION COEFFICIENTS BY DAVID M. BOORE AND JOHN BOATWRIGHT ABSTRACT Averages of P- and
More informationStress associated coda attenuation from ultrasonic waveform measurements
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L09307, doi:10.1029/2007gl029582, 2007 Stress associated coda attenuation from ultrasonic waveform measurements Meng-Qiu Guo 1 and Li-Yun Fu 1 Received 8 February
More informationComposite memory variables for viscoelastic synthetic seismograms
Geophys. J. Int. (1995) 121,634-639 RESEARCH NOTE Composite memory variables for viscoelastic synthetic seismograms Tong Xu George A. McMechan Center for Lithospheric Studies, The University of Texas at
More informationDistortion of the apparent S-wave radiation pattern in the high-frequency wavefield: Tottori-Ken Seibu, Japan, earthquake of 2000
Geophys. J. Int. (2009) 178, 950 961 doi: 10.1111/j.1365-246X.2009.04210.x Distortion of the apparent S-wave radiation pattern in the high-frequency wavefield: Tottori-Ken Seibu, Japan, earthquake of 2000
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
TRANSITION ZONE WAVE PROPAGATION: CHARACTERIZING TRAVEL-TIME AND AMPLITUDE INFORMATION Peter M. Shearer and Jesse F. Lawrence University of California San Diego, Institute of Geophysics and Planetary Physics
More informationAn intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L24619, doi:10.1029/2004gl021228, 2004 An intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake Changjiang
More informationWaveform search for the innermost inner core
Waveform search for the innermost inner core Vernon F. Cormier 1 and Anastasia Stroujkova 1,2 University of Connecticut Storrs, CT 06269-3046 Abstract Waveforms of the PKIKP seismic phase in the distance
More informationVelocity Changes of Seismic Waves and Monitoring Stress in the Crust
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 12,**- pp.,+3,,0 + * +, +, Velocity Changes of Seismic Waves and Monitoring Stress in the Crust Muneyoshi Furumoto + *, Yoshihiro Hiramatsu + and Takashi Satoh,
More informationSEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS
3 Underwater propagation 3. Ray acoustics 3.. Relevant mathematics We first consider a plane wave as depicted in figure. As shown in the figure wave fronts are planes. The arrow perpendicular to the wave
More informationRadiation pattern in homogeneous and transversely isotropic attenuating media
Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010
More informationDynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake
Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 016 Mw 7.0 Kumamoto Earthquake Heng-Yi Su 1 *, Aitaro Kato 1 Department of Earth Sciences, National Central University, Taoyuan
More informationTomographic imaging of P wave velocity structure beneath the region around Beijing
403 Doi: 10.1007/s11589-009-0403-9 Tomographic imaging of P wave velocity structure beneath the region around Beijing Zhifeng Ding Xiaofeng Zhou Yan Wu Guiyin Li and Hong Zhang Institute of Geophysics,
More informationHigh Resolution Imaging of Fault Zone Properties
Annual Report on 1998-99 Studies, Southern California Earthquake Center High Resolution Imaging of Fault Zone Properties Yehuda Ben-Zion Department of Earth Sciences, University of Southern California
More informationSmall-scale heterogeneities in the oceanic lithosphere inferred from guided waves
GEOPHYSICAL RESEARCH LETTERS, VOL. 40, 1708 1712, doi:10.1002/grl.50330, 2013 Small-scale heterogeneities in the oceanic lithosphere inferred from guided waves Azusa Shito, 1 Daisuke Suetsugu, 1 Takashi
More informationSeismic ray path variations in a 3D global velocity model
Physics of the Earth and Planetary Interiors 141 (2004) 153 166 Seismic ray path variations in a 3D global velocity model Dapeng Zhao, Jianshe Lei Geodynamics Research Center, Ehime University, Matsuyama
More informationTHE EFFECT OF DIRECTIVITY ON THE STRESS PARAMETER DETERMINED FROM GROUND MOTION OBSERVATIONS
Bulletin of the Seismological Society of America, Vol. 79, No. 6, pp. 1984-1988, December 1989 THE EFFECT OF DIRECTIVITY ON THE STRESS PARAMETER DETERMINED FROM GROUND MOTION OBSERVATIONS BY DAVID M. BOORE
More informationboundaries with additional record sections, as emphasized in Fig. S2. The observations at the
Data used to Sample African Anomaly. The great circle paths of the source-receiver combinations used in this study are shown in Fig. S1. The event information is given in Table S1. Abrupt Changes across
More informationACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT
ACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT Edward H FIELD 1 And SCEC PHASE III WORKING GROUP 2 SUMMARY Probabilistic seismic hazard analysis
More informationSUPPLEMENTARY INFORMATION
Inability of additional parameters to resolve the Rayleigh-Love discrepancy Radial anisotropy is introduced to resolve the Rayleigh-Love misfit discrepancy that exists across large regions of the western
More informationFrequency-dependent attenuation of S-waves in the Kanto region, Japan
Earth Planets pace, 61, 67 75, 2009 Frequency-dependent attenuation of -waves in the Kanto region, Japan Kazuo Yoshimoto and Mariko Okada International Graduate chool of Arts and ciences, Yokohama City
More informationThree Dimensional Simulations of Tsunami Generation and Propagation
Chapter 1 Earth Science Three Dimensional Simulations of Tsunami Generation and Propagation Project Representative Takashi Furumura Authors Tatsuhiko Saito Takashi Furumura Earthquake Research Institute,
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA
MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA Ali Fatehi and Keith D. Koper Saint Louis University Sponsored by the Air Force Research Laboratory ABSTRACT Contract No. FA8718-06-C-0003 We have used data
More informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
More informationNegative repeating doublets in an aftershock sequence
LETTER Earth Planets Space, 65, 923 927, 2013 Negative repeating doublets in an aftershock sequence X. J. Ma and Z. L. Wu Institute of Geophysics, China Earthquake Administration, 100081 Beijing, China
More informationGEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi: /2004gl020366, 2004
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi:10.1029/2004gl020366, 2004 Characteristic seismic activity in the subducting plate boundary zone off Kamaishi, northeastern Japan, revealed by precise
More informationShort Note Constituent Energy of Regional Seismic Coda
Bulletin of the Seismological Society of America, Vol. 98, No. 1, pp. 454 462, February 28, doi: 1.1785/127121 Short Note Constituent Energy of Regional Seismic Coda by Tae-Kyung Hong and William Menke
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2014) Geophysical Journal International Advance Access published February 21, 2014 doi: 10.1093/gji/ggu031 Strong seismic wave scattering in the low-velocity
More informationGlobal propagation of body waves revealed by cross-correlation analysis of seismic hum
GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 1691 1696, doi:1.12/grl.5269, 213 Global propagation of body waves revealed by cross-correlation analysis of seismic hum K. Nishida 1 Received 1 January 213; revised
More informationSeismogram Interpretation. Seismogram Interpretation
Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for
More informationUPPER MANTLE ATTENUATION STRUCTURE BENEATH THE EASTERN HOKKAIDO, JAPAN AND ITS EFFECTS ON STRONG GROUND MOTIONS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 914 UPPER MANTLE ATTENUATION STRUCTURE BENEATH THE EASTERN HOKKAIDO, JAPAN AND ITS EFFECTS ON STRONG
More informationProbing Mid-Mantle Heterogeneity Using PKP Coda Waves
Probing Mid-Mantle Heterogeneity Using PKP Coda Waves Michael A.H. Hedlin and Peter M. Shearer Cecil H. and Ida M. Green Institute of Geophysics and Planetary Physics Scripps Institution of Oceanography,
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara DEFINITION OF NORMALIZED ENERGY DENSITY AND ITS APPLICATION
More informationTime dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations?
Physics of the Earth and Planetary Interiors 122 (2000) 221 228 Time dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations? Xiaodong Song Department of Geology,
More informationSPATIAL DISTRIBUTION OF STRONG GROUND MOTION CONSIDERING ASPERITY AND DIRECTIVITY OF FAULT
SPATIAL DISTRIBUTION OF STRONG GROUND MOTION CONSIDERING ASPERITY AND DIRECTIVITY OF FAULT Shunroku YAMAMOTO SUMMARY Waveform simulations of the 995 Hyogo-ken Nanbu earthquake were carried out to study
More informationEstimation of deep fault geometry of the Nagamachi-Rifu fault from seismic array observations
Earth Planets Space,,, Estimation of deep fault geometry of the Nagamachi-Rifu fault from seismic array observations Ayako Nakamura, Youichi Asano, and Akira Hasegawa Research Center for Prediction of
More informationarxiv:physics/ v2 [physics.geo-ph] 18 Aug 2003
Is Earthquake Triggering Driven by Small Earthquakes? arxiv:physics/0210056v2 [physics.geo-ph] 18 Aug 2003 Agnès Helmstetter Laboratoire de Géophysique Interne et Tectonophysique, Observatoire de Grenoble,
More informationUnphysical negative values of the anelastic SH plane wave energybased transmission coefficient
Shahin Moradi and Edward S. Krebes Anelastic energy-based transmission coefficient Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient ABSTRACT Computing reflection
More informationvolcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1
volcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1 Università della Calabria, Cosenza, Italy 2 Istituto Nazionale di Geofisica e Vulcanologia Osservatorio Vesuviano,
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationA GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR
A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR Annemarie CHRISTOPHERSEN 1 And Euan G C SMITH 2 SUMMARY This paper considers the distribution of aftershocks in space, abundance, magnitude and time. Investigations
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/326/5949/112/dc1 Supporting Online Material for Global Surface Wave Tomography Using Seismic Hum Kiwamu Nishida,* Jean-Paul Montagner, Hitoshi Kawakatsu *To whom correspondence
More informationMultifractal Analysis of Seismicity of Kutch Region (Gujarat)
P-382 Multifractal Analysis of Seismicity of Kutch Region (Gujarat) Priyanka Midha *, IIT Roorkee Summary The geographical distribution of past earthquakes is not uniform over the globe. Also, in one region
More informationSeismic Velocity Structure in the Crust and Upper Mantle beneath Northern Japan
J. Phys. Earth, 42, 269-301, 1994 Seismic Velocity Structure in the Crust and Upper Mantle beneath Northern Japan Hiroki Miyamachi,l,* Minoru Kasahara,2 Sadaomi Suzuki,2,** Kazuo Tanaka,3 and Akira Hasegawa
More informationLateral variation of the D 00 discontinuity beneath the Cocos Plate
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15612, doi:10.1029/2004gl020300, 2004 Lateral variation of the D 00 discontinuity beneath the Cocos Plate T. Lay Earth Sciences Department, University of California,
More informationCorrelogram Analyses of Seismograms by Means of. By Keiiti AKI Geophysical Institute, Faculty of Science, Tokyo University.
JOURNAL OF PHYSICS OF THE EARTH, VOL. 4, No. 2, 1956 71 Correlogram Analyses of Seismograms by Means of a Simple Automatic Computer. By Keiiti AKI Geophysical Institute, Faculty of Science, Tokyo University.
More informationRECIPE FOR PREDICTING STRONG GROUND MOTIONS FROM FUTURE LARGE INTRASLAB EARTHQUAKES
RECIPE FOR PREDICTING STRONG GROUND MOTIONS FROM FUTURE LARGE INTRASLAB EARTHQUAKES T. Sasatani 1, S. Noguchi, T. Maeda 3, and N. Morikawa 4 1 Professor, Graduate School of Engineering, Hokkaido University,
More informationSynthetic sensitivity analysis of high frequency radiation of 2011 Tohoku-Oki (M W 9.0) earthquake
Earthq Sci (214) 27(4):355 364 DOI 1.17/s11589-14-88-6 RESEARCH PAPER Synthetic sensitivity analysis of high frequency radiation of 211 Tohoku-Oki (M W 9.) earthquake Haoran Meng Yongshun John Chen Received:
More information3D waveform simlation in Kobe of the 1995 Hyogoken-Nanbu earthquake by FDM using with discontinuous grids
3D waveform simlation in Kobe of the 1995 Hyogoken-Nanbu earthquake by FDM using with discontinuous grids S. Aoi National Research Institute for Earth Science and Disaster Prevention H. Sekiguchi, T. Iwata
More informationJOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109, B06306, doi: /2003jb002761, 2004
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003jb002761, 2004 Imaging the fault zones of the 2000 western Tottori earthquake by a new inversion method to estimate three-dimensional distribution
More informationRadiative Transfer of Seismic Waves
Radiative Transfer of Seismic Waves L. Margerin CEREGE, CNRS, Aix en Provence, France Atelier sur les ondes élastiques, Col de Porte, 13 janvier 2009 En collaboration avec: N. Le Bihan, M. Campillo, E.
More informationNonlinear site response from the 2003 and 2005 Miyagi-Oki earthquakes
LETTER Earth Planets Space, 58, 1593 1597, 2006 Nonlinear site response from the 2003 and 2005 Miyagi-Oki earthquakes Kenichi Tsuda and Jamison Steidl Department of Earth Science and Institute for Crustal
More informationON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL
1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI
More informationSeismic scatterers within subducting slab revealed from ambient noise autocorrelation
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl053321, 2012 Seismic scatterers within subducting slab revealed from ambient noise autocorrelation Yoshihiro Ito 1 and Katsuhiko Shiomi 2 Received
More informationSeismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone
IJMS 2017 vol. 4 (2): 49-54 International Journal of Multidisciplinary Studies (IJMS) Volume 4, Issue 2, 2017 DOI: http://doi.org/10.4038/ijms.v4i2.22 Seismic Activity near the Sunda and Andaman Trenches
More informationInvestigation of long period amplifications in the Greater Bangkok basin by microtremor observations
Proceedings of the Tenth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Pacific 6-8 November 2015, Sydney, Australia Investigation of long period amplifications in the Greater
More informationImprovement in the Fault Boundary Conditions for a Staggered Grid Finite-difference Method
Pure appl. geophys. 63 (6) 977 99 33 553/6/9977 DOI.7/s-6-8- Ó Birkhäuser Verlag, Basel, 6 Pure and Applied Geophysics Improvement in the Fault Boundary Conditions for a Staggered Grid Finite-difference
More information