advances in earthquake location and tomography

61
Advances in Earthquake Location and Tomography William Menke Lamont-Doherty Earth Observatory Columbia University

Upload: ciqala

Post on 07-Jan-2016

21 views

Category:

Documents


0 download

DESCRIPTION

Advances in Earthquake Location and Tomography. William Menke Lamont-Doherty Earth Observatory Columbia University. Outline. Part 1: Advantage of using differential arrival times to locate earthquakes Part 2: Simultaneous earthquake location and tomography - PowerPoint PPT Presentation

TRANSCRIPT

Page 1: Advances in Earthquake Location and Tomography

Advances in Earthquake Locationand

Tomography

William Menke

Lamont-Doherty Earth Observatory

Columbia University

Page 2: Advances in Earthquake Location and Tomography

Part 1: Advantage of using differential arrival times to locate earthquakes

Part 2: Simultaneous earthquake location and tomography

Part 3: In depth analysis of the special case of unknown origin time

Outline

Page 3: Advances in Earthquake Location and Tomography

Part 1

Advantage of using differential arrival times to locate earthquakes

Page 4: Advances in Earthquake Location and Tomography
Page 5: Advances in Earthquake Location and Tomography

that was the recent Gulf of Mexico earthquake,

by the way …

Page 6: Advances in Earthquake Location and Tomography

Locating an earthquakerequires knowing the

seismic velocity structure

accurately

Page 7: Advances in Earthquake Location and Tomography
Page 8: Advances in Earthquake Location and Tomography
Page 9: Advances in Earthquake Location and Tomography
Page 10: Advances in Earthquake Location and Tomography

What’s the best way to represent 3 dimensional structure

Best for what?

compatibility with data sources

ease of visualization and editing

facilitating calculation

Page 11: Advances in Earthquake Location and Tomography

Overall organization into interfaces

Small-scale organization into tetrahedra

Linear interpolation within tetrahedra implying rays that are circular arcs

Page 12: Advances in Earthquake Location and Tomography

seismometer earthquake

Page 13: Advances in Earthquake Location and Tomography

Location Errors: = 0.5 degree = 55 km = 30 miles

Note: this preliminary calculation used data from a limited number of stations

Page 14: Advances in Earthquake Location and Tomography

Two parallel approaches

work to improve earth model

design earthquake location techniques that are as insensitive to model as possible

Page 15: Advances in Earthquake Location and Tomography

Waves from earthquake first arrived in Palisades NY at 15:00:32 on Sept

10, 2006

Page 16: Advances in Earthquake Location and Tomography
Page 17: Advances in Earthquake Location and Tomography

Arrival Time ≠Travel Time

Q: a car arrived in town after traveling for an half an hour at sixty miles an hour. Where did it start?

A. Thirty miles away

Q: a car arrived in town at half past one, traveling at sixty miles an hour. Where did it start?

A. Are you crazy?

Page 18: Advances in Earthquake Location and Tomography

Suppose you contour arrival timeon surface of earth

Earthquake’s (x,y) is center

of bullseye

but what about its depth?

Page 19: Advances in Earthquake Location and Tomography

Earthquake’s depth related to

curvature of arrival time at

origin

Deep

Shallow

Page 20: Advances in Earthquake Location and Tomography
Page 21: Advances in Earthquake Location and Tomography
Page 22: Advances in Earthquake Location and Tomography

Courtesty of Felix Walhhauser, LDEO

Earthquakes in Long Valley Caldera, California located with absolute traveltimes

Page 23: Advances in Earthquake Location and Tomography

Courtesty of Felix Walhhauser, LDEO

Earthquakes in Long Valley Caldera, California located with differential traveltimes

Page 24: Advances in Earthquake Location and Tomography

differential arrival time = difference in arrival times

Page 25: Advances in Earthquake Location and Tomography

T = arrival time

TT = travel time

To = Origin Time (start time of earthquake)

mean origin time cancels out

Page 26: Advances in Earthquake Location and Tomography

Station i

Page 27: Advances in Earthquake Location and Tomography
Page 28: Advances in Earthquake Location and Tomography

Very accurate DT’s !

Page 29: Advances in Earthquake Location and Tomography

A technical question for Applied Math types …

Are differential arrival times as calculated by cross-correlation less correlated than implied by the formula

They seem to be.

If so, the this is another advantage of using the method

Page 30: Advances in Earthquake Location and Tomography
Page 31: Advances in Earthquake Location and Tomography

How does differential arrival time vary spatially?

Depends strongly on this angle

Page 32: Advances in Earthquake Location and Tomography

In a 3 dimensional homogeneous box …

maximum

meanminimum

If you can identify the line AB, then you can locate earthquakes

Page 33: Advances in Earthquake Location and Tomography

as long as you have more than two earthquakes

Page 34: Advances in Earthquake Location and Tomography

In a vertically-stratified earth, rays are bent back up to the surface, so both Points A and B are on the surface.

The pattern of differnetial traveltime is more complicated …

ray

wavefront

Page 35: Advances in Earthquake Location and Tomography

The same idea works …

p q

Page 36: Advances in Earthquake Location and Tomography

Patterns of differential arrival time

C

CC

C

C

B

C

B BA

A A

B

A

Can you guess the orientation of the two sources in these six cases?

Page 37: Advances in Earthquake Location and Tomography

This pattern an be seen in actual data, in this case from a pair of earthquakes on the San Andreas Fault

Boxes: differential arrival times observed at particular stations

Shading: theoretical calculation for best-fitting locations of the earthquake pair

C

A

B

Page 38: Advances in Earthquake Location and Tomography

Anotherexample …

Page 39: Advances in Earthquake Location and Tomography

What is the practical advantageof using differential arrival times

to locate earthquakes

My approach is toexamine the statistics of location errorsusing numerical simulations

Compare the result of usingabsolute arrival time data

Anddifferential arrival time data

Whenthe data are noise

Orthe earth structure is poorly known

Page 40: Advances in Earthquake Location and Tomography

Geometry of the numerical experiment …

Page 41: Advances in Earthquake Location and Tomography

Effect of noisy data(10 milliseconds of measurement error)

absolute data

absolute data

differential data

differential data

Page 42: Advances in Earthquake Location and Tomography

Effect of near surface heterogeneities(1 km/s of velocity variation with a scale length of 5 km)

absolute data

differential data differential

dataabsolute data

Page 43: Advances in Earthquake Location and Tomography

Both absolute locations and relative locations of earthquakes are improved by using differential arrival time data

when arrival times are nosily measured andwhen near-surface earth structure is poorly

modeled

Relative location errors can be just a few meters even when errors are “realistically large”

Page 44: Advances in Earthquake Location and Tomography

Part 2

Simultaneous earthquake location and tomography

Page 45: Advances in Earthquake Location and Tomography
Page 46: Advances in Earthquake Location and Tomography
Page 47: Advances in Earthquake Location and Tomography
Page 48: Advances in Earthquake Location and Tomography

simultaneous earthquake location and tomography?

Many earthquakes with unknown X, Y, Z, To

Unknown velocity structure

Solve for everything

Using either

absolute arrival timesor

differential arrival times

Page 49: Advances in Earthquake Location and Tomography

A numerical test

11 stations

50 earthquakeson fault zone

Heterogeneitynear fault zone only

Page 50: Advances in Earthquake Location and Tomography

True earthquake locationsAnd fault zone heterogenity( 1 km/s)

Reconstructed earthquake locationsAnd fault zone heterogenity, using noise free differential data

Note the amplitude of the “signal” is only 1 ms, so noise might be a problem.

Page 51: Advances in Earthquake Location and Tomography

Reality Check: How big is the Signal?

How much better are the data fit?

When the earth structure is allowed to vary

compared with holding a simple, layered

earth structure fixed?

Answer: 0.7 milliseconds, for a dataset that has traveltimes of a few seconds

Need very precise measurements!

Page 52: Advances in Earthquake Location and Tomography

Part 3

Is Joint Tomography/Earthquake Location

Really Possible ?

Study a simplified version of the problem

In depth analysis of the special case of unknown origin timebut known location

Page 53: Advances in Earthquake Location and Tomography
Page 54: Advances in Earthquake Location and Tomography
Page 55: Advances in Earthquake Location and Tomography

Station 1 2 3 4

Event 1

Event 2

Event 3

Page 56: Advances in Earthquake Location and Tomography

If you can …

Then that structure is indistinguishable from a perturbation in origin time!

Page 57: Advances in Earthquake Location and Tomography

Case of sources near bottom of the model

This velocity perturbation causes constant travel time perturbation for a station on the surface anywhere in the grey box for the event at but zero traveltime perturbation for all the sources at !

Page 58: Advances in Earthquake Location and Tomography

Case of sources near top of model

This velocity perturbation causes constant travel time perturbation for a station on the surface anywhere in the grey box for the event at but zero traveltime perturbation for all the sources at !

Page 59: Advances in Earthquake Location and Tomography

But you can always find such structures!

And they often look ‘geologically interesting’

Yet their presence of absence in an area cannot be proved or disproved by the tomography.

Page 60: Advances in Earthquake Location and Tomography

Summary

Part 1: Earthquake location with differential data is the way to go!

Part 2: Simultaneous tomography / earthquake location possible with differential data, but requires high-precision data.

Part 3: Coupled Tomography/Location is extremely nonunique and extremely likely to fool you.

Page 61: Advances in Earthquake Location and Tomography