Modeling of six-component seismogram envelopes ...

7 downloads 0 Views 5MB Size Report
Circles show locations of the August 2011 Vogtland swarm earthquakes greater magnitude 3. (b) Stacked envelopes from 11 earthquakes for the three ...
Modeling of six-component seismogram envelopes using Radiative Transfer Theory Peter Gaebler1,2, Christoph Sens-Schönfelder1 and Michael Korn2 1

GFZ Potsdam, Section 2.4, Seismology, Potsdam

2

University Leipzig, Institute of Geophysics and Geology, Leipzig

Contact: [email protected]

1 Abstract

4 MC-RTT vs FD simulations of six component envelopes

Measurements of rotational motions of the seismic wave-field have recently shown a high amount of rotational energy also in the coda of seismic events. Especially the rotational motion in the P-wave coda bears interesting information as it can only be excited by scattering of the wave-field at 3D heterogeneities. It indicates the conversion from P to S energy in the Earth’s subsurface and consequently the scattering of high-frequency seismic waves. A suitable method to describe the propagation and scattering of high-frequency waves is the Radiative Transfer Theory (RTT). RTT describes the spatial and temporal distribution of seismic energy emitted from a seismic source. It considers scattering and mode conversions between P, SV and SH polarizations. It also includes the angular dependent scattering pattern derived from the Born approximation. The RTT does not contain any phase information of the seismic waves. Therefore the energy of superimposing waves is considered additive and no interferences between the seismic signals are treated. This makes the RTT a very useful tool for the description of scattered waves. It has recently been shown that the RTT is a powerful method to compute seismogram envelopes for the three translational components of the seismic wave-field. In this study we extend the capabilities of the method to model six-component seismogram envelopes in an elastic medium with randomly distributed velocity and density perturbations. The approach implements Monte Carlo (MC) solutions to the radiative transfer equation. The three additional rotational components can provide independent information about the Earth’s structure and the seismic source. They can for example be used to further constrain scattering properties and thus help to discriminate between intrinsic and scattering attenuation. The results of the MC simulations are verified by comparison with 3D full wave-field finite difference simulations. Six-component seismogram envelopes from the two different approaches are compared. We obtain reasonable agreement not only for the translational components but also for the rotational energy. In conclusion, the RTT is a useful approach to model the six-component seismogram envelopes of highfrequency wave-fields from the initial onset of the direct P-wave to the later part of the S-wave coda in random elastic media.

The results of the Monte Carlo Simulations are verified by comparison with 3D finite difference simulations using a modified version of the finite difference modeling software FDMPI [Bohlen,2002]. The simulations are similiar to the ones decribed by Przybilla and Korn, 2008. We use a model size of 76×76×76km with a compressional isotropic source placed at the the center of the model. The dominant source frequency is 6Hz, mean P and S velocities are 6.0km/s and 3.46km/s. Random velocity fluctuations with an exponential autocorrelation function (eACF) are added to the model. The eACF is described by the fractional fluctuation  and by the correlation distance a. We record translational and rotational components of the wave-field at receivers placed on spheres at distances of 10, 20 and 30km from the source (see figure 3a). Figure 3b and 3c show single traces of translational and rotational motions rotated in the LQT-system recorded at 30km distance in a exponential random medium with  = 0.1 and a = 1. L

L Q

Q

10000

0

10000

20000

30000

30000

30000 20000

20000

ec

t io

n

0 10000

[m

]

10000 20000

20000 30000

30000

x-

d ir

ec

n t io

[m

]

T

d ir

10000 0

T

10000

y-

5

0

• bandpass filtering around fc • rotation from ZNE into LQT-system • stacking of single trace hilbert envelopes

N-S Velocity [m/s]

500

0.00000 0.00005 1000

1500

2000

2500

3000

3500

4000

0.0002 0.0001 0.0000 0.0001 0.0002 500

500 2.0 1e 8 1.5 1.0 0.5 0.0 0.5 1.0 1.5500

0.00010700 0.00010

720

740

760

780

800

0.00005

Figure 4 shows a comparison between total (L+Q+T), radial (L) and transversal component (Q+T) intensities for the modeled translations and rotations from FD and MC-RTT simulations. Simulations were carried out in a random medium with an eACF and a velocity fluctuation  = 0.1. The autocorrelation length is 1km.

0.00000

Conclusions:

0.00005 1000

1500

2000

2500

3000

3500

4000

0.0002 0.0001 0.0000 0.0001 0.0002

Z Rotation rate [rad/s]

E-W Velocity [m/s]

0.00005

0.00010700 0.00010

720

740

760

780

800

0.00005 0.00000 0.00005 1000

1500

2000

1000

1500

2000

2500

3000

3500

4000

2500

3000

3500

4000

time[s]

0.00010700 2.0 1e 9 1.5 1.0 0.5 0.0 0.5 1.0 1.5 2.0700

720

740

720

740

760

780

800

760

780

800

time[s]

Figure 1: Top three traces show ground velocities, bottom trace shows rotation rate around the vertical axis. P-wave arrival is marked in orange.

3 RTT - Radiative Transfer Theory RTT describes the the energy flux from a radiating surface into a unit solid angle around direction n, this quantity is called specific intensity I(n, r). It depends on position r, as well as on the direction of the energy flux n. Phenomenologically the transport problem can be regarded as a conservation of energy flux in a direction n and at position r. Figure 2 illustrates the conservation of energy flux in the radiative transfer equations.

10

1

10

0

10

3

10

2

10

1

10

0

0

10

3

10

2

4

2

6

8

10

12

14

16

4

2

6

8

10

12

14

10

3

10

2

10

1

10

0

18

ra dia l rota tions (inte ns ity )

Z Velocity [m/s]

0.00010

2

16

10

3

10

2

10

1

10

0

18

tra ns v e rs a l rota tions (inte ns ity )

• determination of seimic phase velocities from collocated translational and rotational measurements

10

FD 10km dis ta nce MC 10km dis ta nce FD 20km dis ta nce MC 20km dis ta nce FD 30km dis ta nce MC 30km dis ta nce

0

ra dia l tra ns la tions (inte ns ity )

• recording particles in LQT-system • stacking of single trace hilbert envelopes • simulate FD calculations by → correcting for source time function → correcting for wandering effect

Tra ns la tiona l m otions

3

tota l rota tions (inte ns ity )

FD-simulations processing

10

tota l tra ns la tions (inte ns ity )

To compare the six-component envelopes from the FD and MC simulations following data processing steps were done:

• use of rotation sensor as wave polarization filter

0.0002 0.0001 0.0000 0.0001 0.0002

5

0

20

15

10

tim e [s ]

20

transversal axis recorded at 30km distance. Vertical orange lines indicate P and S arrival times for waves propagating vP , vS .

MC-RTT-simulations processing

2011/03/11 Tohoku Earthquake, Magnitude 9.0, Bandpass 0.1-5Hz

15

10

tim e [s ]

Figure 3: a) Receivers placed at 10, 20, 30km distance b) Radial and transversal translational motions recorded at 30km distance c) Rotations around the radial and the

• incorporation of rotational motions into the seismic inverse problem

Figure 1 shows translational and vertical rotational ground motion induced by the March 2011 Tohoku earthquake recorded at the Fundamentalstation Wettzell, Germany. Rotational motions around a vertical axis can only be excited by horizontally polarized shear waves. The high amount of energy visible in the vertical rotation rate at the time of the P-wave arrival therefore clearly indicates a conversion from P to S energy. This can be achieved by scattering of the high-frequency seismic wave-field at local or regional small scale heterogeneities.

LQT ROTATIONS

20000

tra ns v e rs a l tra ns la tions (inte ns ity )

Additionally to the three components of translational motion the seismic wave-field is described by six components of strain and three components of rotational motion [Aki and Richards, 2002]. Due to the steady increase in sensor sensitivity it is nowadays possible to measure rotational motions in the seismic wave field induced by earthquakes and even by ambient seismic noise with very high accuracy. There are several applications using the records of rotational motions, for example:

c

LQT TRANSLATIONS

30000

z-d ire c tio n [m ]

2 Rotational Motions in Seismology

b

a

10

3

10

2

Rota tiona l m otions

0

2

4

6

8

10

12

14

16

18

0

2

4

6

8

10

12

14

16

18

Comparison between MC-RTT and FD 10 10 envelopes show a very good agreement from the initial P-wave onset until late parts of the S-wave coda. Modelling of 10 10 translational motions agrees with previous 4 14 4 14 0 2 6 8 10 12 16 18 0 2 6 8 10 12 16 18 tim e [s ] tim e [s ] results by Przybilla and Korn, 2008. Figure 4: Comparison between MC-RTT and FD simulations The results of this study clearly demonstrate that also the three rotational components of motion can be modeled using RTT. We therefore conclude, that the RTT approach can be used to analyze six-component-envelopes of the high-frequency seismic wave-field. 1

0

1

0

5 Real data application Figure 5b shows stacked hilbert envelopes for 11 earthquakes at 10 stations (see figure 5a). The earthquake station distance ranges from 100 to 160km. Additionally to the three translational envelopes (black) we measure a vertical rotation envelope (red) using data from the ringlaser instrument located at the Fundamentalstation Wettzell. These four component records will be inverted for following medium parameters using MC-RTT simulations: • fractional fluctuation  • correlation length a • quality factor for intrinsic S-wave attenuation QS

a

10˚

11˚

12˚

13˚

b

14˚

Z-com pone nt

R-com pone nt

170

T-com pone nt

170

170

GR_GRC3

GR_GRC3 160

GR_GRC3

160

160

GR_GRC1

GR_GRC1 150

GR_GRC1

150

150

Figure 2: Energy conservation in the radiative transfer equations. Illustrated is the energy flux from a radiating surface (RS) into a unit solid angle (USA) around a direction

GR_GRA2 GR_GRA4

The Radiative Transfer Equations are solved using Monte Carlo techniques. In this study a so called Particle Counting Method is used. The specific intensity I(n, r) is modeled by a number density of particles N (n, r) moving into direction n and loctated a position r. Particles are emitted from a source following a fixed radiation pattern and are described by following parameters:

GR_GRB4

GR_GRB2

140

GR_GRB5

130

BW_RLAS GR_WET GR_GRB2

120

49˚

GR_GRC1 GR_GRC3

110

During the course of the simulation run particles can undergo scattering processes at medium heterogeneities following the Born approximation. This processes include mode conversion and a change of propagation direction. To account for attenuation particles loose energy (I P , IθS , IφS ) as they are propagating through the medium. When no scattering events occur particles move through the medium according to ray theory. This includes the interaction with interfaces (reflection, transmission, mode conversions) or the propagation in depth-dependent velocity structures. Snapshots of the energy field are stored in an four-dimensional-array consisting of three spatial and one temporal dimension.

130

GR_WET GR_GRB2

120

GR_GRB4

GR_GRA4

GR_GRA4

GR_GRA2

100

11˚

130

GR_WET GR_GRB2

120

12˚

13˚

14˚

GR_GRB4 GR_GRA4

GR_GRA2

100

GR_GRA3

10˚

GR_GRB5

110

GR_GRB4

GR_GRA3

90

0

20

40

60

tim e [s ]

80

GR_GRA2

100

GR_GRA3

90

• velocity vector: vx, vy , vz • Stokes Vector: (I P , IθS , IφS , U, V )

GR_GRB5

110

• mode: P, S • coordinates in space and time: x, y, z, t

140

GR_WET

GR_GRB5 GR_GRC4

49˚

GR_GRC4

s ta tion dis ta nc e [km ]

sources (D). Figure redrawn from Sens-Schönfelder et al., 2009.

140

s ta tion dis ta nc e [km ]

GR_GRA3

GR_GRC4

GR_GRC4

50˚ s ta tion dis ta nc e [km ]

n. The total flux is reduced by scattering in other directions (A) and by dissipation (B). A flux increase results from scattering from other directions into direction n and from

50˚

100

90

0

20

40

60

tim e [s ]

80

100

0

20

40

60

tim e [s ]

80

100

Figure 5: (a) Triangles indicate location of Gräfenberg array stations plus the location of the Fundamentalstation Wettzell. Circles show locations of the August 2011 Vogtland swarm earthquakes greater magnitude 3. (b) Stacked envelopes from 11 earthquakes for the three translational components of the shown stations (black). Rotation rate around the vertical axis measured with the ringlaser in Wettzell is shown in red. The blue and orange line indicate the expected arrival of the P respectively the S-wave, assuming vP = 5.6km/s and vS = 3.4km/s.

6 References • Aki and Richards. 2002. Quantitative Seismology. University Science Books. • Bohlen. 2002. Parallel 3-D viscoelastic finite difference seismic modelling. Computers and Geosciences. • Przybilla and Korn. 2008. Monte Carlo simulation of radiative energy transfer in continuous elastic random media-three-component envelopes and numerical validation. Geophysical Journal International. • Sens-Schönfelder, Margerin and Campillo. 2009. Laterally heterogeneous scattering explains Lg blockage in the Pyrenees. Journal of Geophysical Research.