1 vertically integrated seismic analysis stuart russell computer science division, uc berkeley nimar...
TRANSCRIPT
![Page 1: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/1.jpg)
1
Vertically Integrated Seismic Analysis
Stuart RussellComputer Science Division, UC Berkeley
Nimar Arora, Erik Sudderth, Nick Hay
![Page 2: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/2.jpg)
2
Outline
Seismic event monitoring as probabilistic inference Vertically integrated probability models …
Connect events to sensor data and everything in between Associate events and detections optimally Automatically take nondetections into account May improve low-amplitude detection and noise rejection
Inference using MCMC (poster) Empirical estimation of model components Preliminary experimental results
![Page 3: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/3.jpg)
3
Bayesian model-based learning
Generative approach P(world) describes prior over what is (source), also over
model parameters, structure P(signal | world) describes sensor model (channel) Given new signal, compute
P(world | signal) ~ P(signal | world) P(world)
Learning Adapt model parameters or structure to improve fit Operates continuously as data are acquired and analyzed
Substantial recent advances in modeling capabilities, general-purpose inference algorithms
![Page 4: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/4.jpg)
4
Generative model for IDC arrival data
Events occur in time and space with magnitude Natural spatial distribution a mixture of Fisher-Binghams Man-made spatial distribution uniform Time distribution Poisson with given spatial intensity Magnitude distribution Gutenberg-Richter Aftershock distribution (not yet implemented)
Travel time according to IASPEI91 model+corrections Detection depends on magnitude, distance, station* Detected azimuth, slowness w/ empirical residuals False detections with station-dependent distribution
![Page 5: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/5.jpg)
5
Seismic event
Travel times
Seismic event
Travel times
Station 1picks
Station 2picks
Generative structure
Detected atStation 1?
Detected atStation 2?
Station 1noise
Station 2 noise
![Page 6: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/6.jpg)
6
Inference
MCMC (Markov chain Monte Carlo) (see poster S31B-1713 for details)
Efficient sampling of hypothetical worlds (events, travel times, detections, noise, etc.)
Converges to true posterior given evidence Key point: computing posterior probabilities takes the algorithm off the table; to get better answers, either Improve the model, or Add more sensors
![Page 7: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/7.jpg)
7
Vertical integration: Detection
Basic idea: analyzing each signal separately throws away information. Multiple weak signals are mutually reinforcing via a higher-
level hypothesis Multiple missing signals indicate that other “detections” may
be coincidental noise Simple example: K sensors record either
Independent noise drawn from N[0,1] Common signal drawn from N[0,1-] + independent N[0,] noise
Separate detectors fail completely!
Joint detection succeeds w.p. 1 as 0 or K Travel time accuracy affects detection capability!
![Page 8: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/8.jpg)
8
STA/LTA Threshold
![Page 9: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/9.jpg)
9
![Page 10: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/10.jpg)
10
![Page 11: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/11.jpg)
11
![Page 12: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/12.jpg)
12
Outline
Seismic event monitoring as probabilistic inference Vertically integrated probability models …
Connect events to sensor data and everything in between Associate events and detections optimally Automatically take nondetections into account May improve low-amplitude detection and noise rejection
Inference using MCMC (poster) Empirical estimation of model components Preliminary experimental results
![Page 13: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/13.jpg)
13
Seismic event
Travel times
Seismic event
Travel times
Station 1picks
Station 2picks
Generative structure
Detected atStation 1?
Detected atStation 2?
Station 1noise
Station 2 noise
![Page 14: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/14.jpg)
14
![Page 15: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/15.jpg)
15
![Page 16: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/16.jpg)
16
Seismic event
Travel times
Seismic event
Travel times
Station 1picks
Station 2picks
Generative structure
Detected atStation 1?
Detected atStation 2?
Station 1noise
Station 2 noise
![Page 17: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/17.jpg)
17
![Page 18: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/18.jpg)
18
![Page 19: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/19.jpg)
19
Seismic event
Travel times
Seismic event
Travel times
Station 1picks
Station 2picks
Generative structure
Detected atStation 1?
Detected atStation 2?
Station 1noise
Station 2 noise
![Page 20: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/20.jpg)
20
Overall Pick Error
![Page 21: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/21.jpg)
21
WRA Pick Error
![Page 22: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/22.jpg)
22
Overall IASPEI Error
![Page 23: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/23.jpg)
23
WRA - IASPEI Error
![Page 24: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/24.jpg)
24
Seismic event
Travel times
Seismic event
Travel times
Station 1picks
Station 2picks
Generative structure
Detected atStation 1?
Detected atStation 2?
Station 1noise
Station 2 noise
![Page 25: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/25.jpg)
25
Overall Azimuth Error
![Page 26: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/26.jpg)
26
WRA - Azimuth Error
![Page 27: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/27.jpg)
27
Seismic event
Travel times
Seismic event
Travel times
Station 1picks
Station 2picks
Generative structure
Detected atStation 1?
Detected atStation 2?
Station 1noise
Station 2 noise
![Page 28: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/28.jpg)
28
![Page 29: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/29.jpg)
29
Analyzing Performance
Min-cost max-cardinality matching where edges exist between prediction and ground truth events within 50 seconds and 5 degrees.
Precision – percentage of predictions that match. Recall – percentage of ground truths that match. F1 – harmonic mean of precision and recall. Error – average distance between matching events.
(Cost of matching / size of matching)
![Page 30: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/30.jpg)
30
Evaluation vs LEB (human experts)
F1 Precision/
Recall
Error/S.D.
(km)
Average Log-likelihood
SEL3 (IDC Automated)
55.6 46.2 / 69.7 98 / 119 _
VISA (Best Start) 80.4 70.9 / 92.9 100 / 117 -1784
VISA (SEL3 Start) 55.2 44.3 / 73.4 104 / 124 -1791
VISA (Back projection Start)
50.6 49.1 / 52.0 126 / 139 -1818
![Page 31: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/31.jpg)
31
INFERENCE EXAMPLE
model2_71_223.avi
![Page 32: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/32.jpg)
32
Summary
Vertically integrated probability models Connect events, transmission, detection,
association Information flows in all directions, reinforcing or
rejecting local hypotheses to form a global solution Better travel time model => better signal detection Nondetections automatically play a role Local sensor models calibrated continuously with no
need for ground truth May give more reliable detection and localization of
lower-magnitude events
![Page 33: 1 Vertically Integrated Seismic Analysis Stuart Russell Computer Science Division, UC Berkeley Nimar Arora, Erik Sudderth, Nick Hay](https://reader038.vdocument.in/reader038/viewer/2022110322/56649d215503460f949f7076/html5/thumbnails/33.jpg)
33
Ongoing Work
More sophisticated MCMC design Add more phases and phase relabeling Extend model all the way down to waveforms Evaluation using data from high-density
networks (Japan Meteorological Agency, some regions within ISC data)