SEISAN Earthquake Analysis Software Lars Ottemöller and Jens

SEISAN Earthquake Analysis Software Lars Ottemöller and Jens Havskov (University of Bergen) Peter Voss (GEUS)

Background • • • Development started in the late 80 s Large number of software contributors Steady improvements, now Version 9. 1 Basics unchanged Supported operating systems: Solaris, Windows, Linux and Mac. OS

Who uses SEISAN • Geographically widely used (Norway, UK, France, Spain, . . . , India, Iran, Philippines, . . . , Central and South America, South and East Africa, . . . ) • 15% of data sent to ISC in SEISAN format • National and regional seismic networks • Volcano observatories • Researchers

What can you do with SEISAN • Import/Export data • Store data (waveform and parametric) • Observatory routine: – Analyse earthquake data (signal processing, phase identification, . . . ) – Compute earthquake parameters (location, magnitude, . . . ) – Fault plane solution – Seismicity maps – Produce bulletin and catalogue search • Research platform, several additional data processing tools

Philosophy • • • Open source Open system Non-commercial Multi-platform Integrate essential tools and data storage into one package • Provide documentation • Provide support • Provide training

SEISAN Help • Manual (PDF and HTML) • Training document • Seisan Webpage http: //www. uib. no/rg/geodyn/artikler/201 0/02/software • Mailing list: seisan@geo. uib. no • Authors: jens. havskov@geo. uib. no lars. ottemoller@geo. uib. no pv@geus. dk

Other software • SAC (http: //www. iris. edu/manuals/sac/index. htm) • Geotool (CTBTO) • Seismic Handler (http: //www. szgrf. bgr. de/sh-doc/index. html) • Giant/Pitsa (http: //www. geo. unipotsdam. de/forschung/software/giant. html) • … • Link on Orefus and IRIS

Data Types • Waveform data • Calibration data • Parametric data

Waveform data • Created by digitising analogue output from seismometer and time stamping • Consists of amplitude and time • Translates into ground motion using instrument response information • Written out by data acquisition system • Waveformats understood directly by SEISAN: SEISAN, SEED, mini. SEED, SAC, GSE

Waveform data files • Single or multiple channels • Data identifier: Station code, Channel code, Network code, Location code • Timing information • Data values • Used for: Recording, analysis, archival, data exchange, …

Waveform data into SEISAN • Conversion of data created by data logger • Conversion of data provided by another agency (IRIS, ORFEUS, . . . ) • Data written out by near real-time processing system (Earthworm, Seislog, Seis. Com. P, . . . ) • Automated data collection: SEISNET

Waveform Conversion Examples • Rdseed: • Wavetool: • … SEED -> mini. SEED -> SAC … (mini)SEED <-> SEISAN SAC <-> SEISAN GSE <-> SEISAN

Parametric data • Derived from waveform data • Observations: – Arrival times and residuals – Amplitudes – Polarities • Computed: – Hypocenter location – Magnitudes – Fault plane solution

SEISAN Directory Structure 2005 2006 BER 2005 INF WOR Workspace BER DAT Documentation TEST CAL Parameter Files WAV REA Response Files C: seismo SEISAN-TOP 01 02 03. . . Parametric files Waveform files

Main SEISAN programs • EEV: Database tool • MULPLT: Seismogram analysis • HYPOCENTER: Location of local, regional and teleseismic earthquakes • EPIMAP: Seismicity maps • FOCMEC: Fault plane solution • SELECT: Extract parametric data

EEV

Main SEISAN programs • EEV: Database tool • MULPLT: Seismogram analysis • HYPOCENTER: Location of local, regional and teleseismic earthquakes • EPIMAP: Seismicity maps • FOCMEC: Fault plane solution • SELECT: Extract parametric data

MULPLT

MULPLT

MULPLT

Main SEISAN programs • EEV: Database tool • MULPLT: Seismogram analysis • HYPOCENTER: Location of local, regional and teleseismic earthquakes • EPIMAP: Seismicity maps • FOCMEC: Fault plane solution • SELECT: Extract parametric data

HYPOCENTER MAGNITUDES All main magnitude scales can be used: - Local magnitude Ml - Coda magnitude Mc - Body wave magnitide mb, m. B - Surface wave magnitude Ms, MS - Moment magnitude Mw

HYPOCENTER

Main SEISAN programs • EEV: Database tool • MULPLT: Seismogram analysis • HYPOCENTER: Location of local, regional and teleseismic earthquakes • EPIMAP: Seismicity maps • FOCMEC: Fault plane solution • SELECT: Extract parametric data

EPIMAP

Main SEISAN programs • EEV: Database tool • MULPLT: Seismogram analysis • HYPOCENTER: Location of local, regional and teleseismic earthquakes • EPIMAP: Seismicity maps • FOCMEC: Fault plane solution • SELECT: Extract parametric data

FOCMEC

Main SEISAN programs • EEV: Database tool • MULPLT: Seismogram analysis • HYPOCENTER: Location of local, regional and teleseismic earthquakes • EPIMAP: Seismicity maps • FOCMEC: Fault plane solution • SELECT: Extract parametric data

SELECT

Format Examples GSE Waveform DATA_TYPE WAVEFORM GSE 2. 0 WID 2 2007/06/04 17: 34: 54. 000 ESK HHZ 0. 10 E+01 1. 000 -1. 0 0. 0 STA 2 55. 31650 -3. 20517 DAT 2 -986 -988 -990 -987 -992 -984 -982 -992 -984. . . INT 40000 100. 000000 0. 261 0. 000

Format Examples GSE Calibration DATA_TYPE RESPONSE GSE 2. 0 CAL 2 ESK HHZ CMG-3 0. 13 E+00 1. 100. 00000 2005/06/01 16: 00 PAZ 2 1 V 1. 12701317 E+03 5 3 Laplace transform -0. 148031845+00 -0. 11309730 E+04 0. 0000 E+00 -0. 10053090 E+04 0. 0000 E+00 -0. 50265500 E+03 0. 00000000 E+00 0. 0000 E+00 DIG 2 2 6. 29326620 E+05 100. 00000 Reftek

Format Examples GSE Phases EVENT 00000001 Date Time rms OT_Error 2007/04/28 07: 18: 11. 3 0. 70 +- 2. 32 Sta TFO 1 MENF APA SKP Dist 9. 18 Ev. Az 338. 0 129. 0 298. 0 m m mc m Latitude Longitude Smajor Sminor Az 51. 0730 0. 0 Phase P S P P 1. 1670 0. 0 0 Date 2007/04/28 2007/04/28 Depth Err Ndef Nsta Gap mdist Mag 1 Err 0. 0 5. 6 156 124 65 0. 04 11. 85 ML 4. 3 41 +-0. 2 +- Time 07: 18: 12. 6 07: 18: 13. 8 07: 18: 27. 2 07: 18: 38. 7 07: 18: 34. 2 07: 18: 37. 2 TRes 0. 1 0. 3 0. 2 0. 3 -0. 2 0. 3 Azim Az. Res Slow N Mag 2 Err SRes Def T T T N Mag 3 Err N Author Quality ID 00000001 m i se SNR 9. 2 3. 8 3. 5 Amp Per Mag 1 Mag 2 Arr ID 00000001 00000002 00000003 00000004 00000005 00000006

Format Examples SEISAN Waveform 21107 155 CWF HH EDI HH GAL 1 HH HPK HH KPL HH LRW HH PGB 1 HH CWF Z Z Z Z 6 4 17 29 59. 000 0. 99 0. 99 HH Z 107 155 1201. 01 1201. 01 6 1202. 000 CWF HH EDI HH GAL 1 HH HPK HH KPL HH LRW HH PGB 1 HH N N N N 4 17 29 59. 990 0. 99 0. 00 1201. 01 1202. 00 CWF HH EDI HH GAL 1 HH HPK HH KPL HH LRW HH PGB 1 HH E E E E 0. 99 0. 00 100. 00 120102 1201. 01 1202. 00 4 … -4231 -4200 -4217 -4231 -4217 -4178 -4195 -4220 -4198 -4208 -4239 -4207 -4178 -4198 -4223 -4201 -4221 -4238 -4198 -4183 -4197 -4227 -4202 -4223 -4230 -4195 -4172 -4186 -4229 -4217 -4216 -4200 -4178 -4196 -4218 -4202 -4233 -4231 -4192 -4179 -4192 -4206 -4205 -4236 -4235 -4183 -4185 -4196

Format Examples SEISAN Calibration ESK S Z 97 044 2 13 17 30 0. 000 55. 3167 -3. 2050 263 Mk 3; 500 V/m/s; fm. A/M(5)40 d. B; Flt. 002&40 Hz@-1&6 p; VME A/D 2048 c/V RMY 13/02/97 1. 000 0. 700 499. 251 40. 0002048. 000. 459 E+09 0. 002 -1. 000 40. 000 6. 000 0. 000. 100. 200. 300. 400. 500. 600. 700. 800. 900 1. 00. 140 E-02. 112 E-01. 377 E-01. 888 E-01. 171. 286. 435. 609. 801 1. 00 -1. 701 -1. 864 -2. 026 -2. 192 -2. 366 -2. 545 -2. 727 -2. 905 -3. 077 3. 047 1. 10 1. 20 1. 60 2. 00 2. 40 3. 00 4. 00 5. 00 6. 00 8. 00 1. 20 1. 39 2. 10 2. 73 3. 32 4. 18 5. 60 7. 00 8. 40 11. 2 2. 902 2. 771 2. 380 2. 129 1. 954 1. 765 1. 542 1. 371 1. 226 0. 970 10. 0 12. 0 14. 0 17. 0 20. 0 23. 0 26. 0 30. 0 40. 0 50. 0 14. 0 16. 8 19. 6 23. 8 28. 0 32. 2 36. 3 41. 4 39. 6 17. 8 0. 738 0. 516 0. 297 -0. 028 -0. 358 -0. 697 -1. 052 -1. 565 -3. 107 1. 917

Format Examples SEISAN Phases 2007 8 1 0413 43. 84 L 56. 463 -6. 531 12. 2 BGS 8 0. 4 1. 4 LBGS 1 GAP=270 1. 08 6. 3 29. 0 12. 4 0. 1431 E+03 -0. 1777 E+03 -0. 1814 E+02 E XNEAR 150. 0 XFAR 250. 0 SDEP 10. 0 3 2007 -08 -01 -0414 -21 S. MOQ___014 6 2007 -08 -01 -0412 -50 S. KYL___011 6 ACTION: REE 07 -08 -02 11: 53 OP: gdf STATUS: ID: 20070801041441 I STAT SP IPHASW D HRMM SECON CODA AMPLIT PERI AZIMU VELO SNR AR TRES W DIS CAZ 7 KAR 1 SZ EP 413 54. 93 95 -0. 0910 66. 6 40 KPL HN ES 414 14. 81 67 -0. 0710 111 28 KPL HN AML 414 15. 16 13. 6 0. 19 111 28 KPL HE AML 414 15. 63 19. 1 0. 18 111 28 KAC SZ EP 414 5. 68 67 0. 1310 138 32 EAB SZ EP 414 5. 55 67 -0. 2410 139 102 PGB 1 HE ES 414 24. 01 67 0. 3710 147 119 PGB 1 HE AML 414 25. 71 3. 3 0. 24 147 119 PGB 1 HN AML 414 26. 70 3. 8 0. 30 147 119 ELO SZ EP 414 9. 49 54 -1. 03 8 174 89 EDI HE ES 414 38. 60 54 -0. 43 3 216 105 EBH SZ EP 414 12. 70 54 0. 39 6 188 96 KPL HZ EP 414 2. 16 67 0. 4010 111 28

SEISAN introductory training course • • • Demonstration Installation Seisan basics Phase picking Earthquake location Magnitudes Fault plane solution Earthquake spectra Participant’s data

Earthworm and SEISAN • Earthworm is a modular automated processing system • Functionality: – – – real-time data import/export archive data signal processing: decimate, filter, . . . helicorder and spectrogram plots detect seismic events locate and determine magnitude • Data output can be in SEISAN format

Earthworm Data Acquisition Data Import Shared Memory Association Verification Picker WAVERING Wave Server SEISAN Data Archival format PICKRING Data Processing Location HYPORING Event Trigger Waveform and parametric detection files