2017 - Gabrielle Simard-Moore - Calculating the Epicenter of an Earthquake with Python

Transcript

1. Calculating the Epicenter of an Earthquake with Python Gabrielle Simard-Moore
2. From Seismograph to Epicenter • 3 Seismic Stations • P- and S-wave arrival time difference • Distance to Earthquake • Draw 3 circles • Intersection = Epicenter
3. Seismic Data
4. Getting the Distance and Magnitude
6. Tech Stack • Python • Numpy • Datetime • Math • Flask • Ajax/Jquery • Google Maps JS API
7. Assumptions • Wave velocity in CA crustal rocks P-wave velocity ~ 6.34 km/s (varies by +/- 2+) S-wave velocity ~ 3.67 km/s (varies by +/- 1+) • Homogeneity of materials (very unrealistic for a heterogeneous lithosphere)
8. 3 stations
9. Many steps
10. Degrees to Radians
11. Geodetic (latitude, longitude) to ECEF X, Y, Z coordinates to Numpy array ECEF = “Earth centered, Earth fixed”
12. Vector Normalization
13. Final Trilateration point: Epicenter
14. ECEF XYZ to Lat, Lon
15. Front end
16. Thanks! github.com/gsimore/earthquakes gsimardmoore@gmail.com