Automated Earthquake Response Systems

1
Automated Earthquake Response Systems

• Tom Heaton, Caltech
• Hiroo Kanamori, Caltech
• Egill Hauksson, Caltech
• Georgia Cua, ETH, Switzerland
• Masumi Yamada, Kyoto Univ
• Maren Böse, Caltech

2
Earthquake Alerting a different kind of
prediction

• What if earthquakes were really slow, like the
  weather?
• We could recognize that an earthquake is
  beginning and then broadcast information on its
  development on the news.
• an earthquake on the San Andreas started
  yesterday. Seismologists warn that it may
  continue to strengthen into a great earthquake
  and they predict that severe shaking will hit
  later today.

3
Shakeout Simulation (Aagaard and Graves)

• ..\Documents\2008\ShakeOut_LosAngeles_lowres.mov

4
(No Transcript)
5
If the earthquake is fast, can we be faster?

• Everything must be automated
• Data analysis that a seismologist uses must be
  automated
• Communications must be automated
• Actions must be automated
• Common sense decision making must be automated

6
How would the system work?

• Seismographic Network computers provide estimates
  of the location, size, and reliability of events
  using data available at any instant estimates
  are updated each second
• Each user is continuously notified of updated
  information . Users computer estimates the
  distance of the event, and then calculates an
  arrival time, size, and uncertainty
• An action is taken when the expected benefit of
  the action exceeds its cost
• In the presence of uncertainty, false alarms must
  be expected and managed

7
What we need is a special seismologist

• Someone who has good knowledge of seismology
• Someone who has good judgment
• Someone who works very, very fast
• Someone who doesnt sleep
• We need a Virtual Seismologist

8
Ground motion envelope our definition
Full acceleration time history
Efficient data transmission 3 components each
of Acceleration, Velocity, Displacement, of 9
samples per second
envelope definition max.absolute value over
1-second window
9
Estimating M from ratios of P-wave motions

• P-wave frequency content scales
• with M (Allen and Kanamori, 2003, Nakamura,
  1988)
• Find the linear combination of log(acc) and
  log(disp) that minimizes the variance within
  magnitude-based groups while maximizing
  separation between groups (eigenvalue problem)

• Estimating M from Zad

10

• Voronoi cells are nearest neighbor regions
• If the first arrival is at SRN, the event must
  be within SRNs Voronoi cell
• Green circles are seismicity in week prior to
  mainshock

11
What about Large Earthquakes with Long Ruptures?

• Large events are infrequent, but they have
  potentially grave consequences
• Large events potentially provide the largest
  warnings to heavily shaken regions
• Point source characterizations are adequate for
  Mlt7, but long ruptures (e.g., 1906, 1857) require
  finite fault

12
Strategy to Handle Long Ruptures

• Determine the rupture dimension by using
  high-frequencies to recognize which stations are
  near source
• Determine the approximate slip (and therefore
  instantaneous magnitude) by using low-frequencies
  and evolving knowledge of rupture dimension
• We are using Chi-Chi earthquake data to develop
  and test algorithms

13
10 seconds after origin
20 seconds after origin
Near-field Far-field
Near-field Far-field
14
30 seconds after origin
40 seconds after origin
Near-field Far-field
Near-field Far-field
15

• Once rupture dimension is known
• Obtain approximate slip from long-periods
• Real-time GPS would be very helpful
• Evolving moment magnitude useful for estimating
  probable rupture length
• Magnitude critical for tsunami warning

16
Real-time prediction of ultimate rupture
B?se and Heaton, in prep.
Remaining Rupture Length
slip
Is the rupture on the San Andreas fault?
17
Distributed and Open Seismic Network

• Just in the gedanken phase
• Tens of thousands of inexpensive seismometers
  running on client computers.
• Sensors in buildings, homes, buisinesses
• Data managed by a central site and available to
  everyone.
• It will change the world!

18
Conclusions

• Earthquake warning is an extension of current
  systems that provide rapid information in
  earthquakes (e.g. ShakeMap)
• Significant warning times can be achieved for
  moderate to light shaking
• Heavily damaged areas in moderate events (M 6)
  will not receive warning
• Heavily damaged areas in great earthquakes (M
  7.5) will receive
• Strategies to determine rupture dimension and
  slip look very promising
• User decision making should be based on
  cost/benefit analysis need to develop a
  community that develops optimal responses
• Managing expectations is critical users must
  understand what EEW wont do.

19
Probabilistic Rupture Prediction ? Probabilistic
Ground Shaking
B?se and Heaton, in prep.
20

• We are experimenting with different Linear
  Discriminant analyses to distinguish near-field
  from far-field records

21
Percent of area receiving warning time T or
greater (log N6.89-Mw)
Pseudovelocity cm/sec
Warning time T sec
Heaton, 1985
22
3 sec after initial P detection at SRN

• Prior information
• Voronoi cells
• Gutenberg-Richter

Single station estimate
M, R estimates using 3 sec observations at SRN
No prior information
8 km M4.4

• Prior information
• Voronoi cells
• No Gutenberg-Richter

9 km M4.8
Note star marks actual M, RSRN
23
Virtual Seismologist (VS) method for seismic
early warning

• Bayesian approach to seismic early warning
  designed for regions with distributed seismic
  hazard/risk
• Modeled on back of the envelope methods of
  human seismologists for examining waveform data
• Shape of envelopes, relative frequency content
• Robust analysis
• Capacity to assimilate different types of
  information
• Previously observed seismicity
• State of health of seismic network
• Known fault locations
• Gutenberg-Richter recurrence relationship

24
Data set for learning the envelope
characteristics Most data are from TriNet, but
many larger records are from COSMOS

• 70 events, 2 lt M lt 7.3, R lt 200 km
• Non-linear model estimation (inversion) to
  characterize waveform envelopes for these events
• 30,000 time histories

25
Average Rock and Soil envelopes as functions of
M, R rms horizontal
acceleration
26
Distinguishing between P- and S-waves
27
(No Transcript)
28
(No Transcript)
29
(No Transcript)
30
(No Transcript)
31
(No Transcript)
32
(No Transcript)
33

• Sum of 9 point source envelopes
• Vertical acceleration

34
horizontal acceleration ampl rel. to ave. rock
site
Vertical P-wave acceleration ampl rel. to ave.
rock site
horizontal velocity ampl rel. to ave. rock site
vertical P-wave velocity ampl rel. to ave. rock
site
35
Strategy for acceleration envelopes

• High-frequency energy is proportional to rupture
  are (Brune scaling)
• Sum envelopes from 10-km patches