Automatic Moment Tensor solution for SeisComP3

Transcription

1 Second European Conference on Earthquake Engineering and Seismology, Istanbul Aug , 2014 Automatic Moment Tensor solution for SeisComP3 Nikolaos Triantafyllis1, Efthimios Sokos2 and Aristidis Ilias1 1 University of Patras, Computer Engineering and Informatics Department, Rio 26504, Patras, Greece, triantafyl@ceid.upatras.gr, aristeid@ceid.upatras.gr 2 University of Patras, Geology Department, Rio 26504, Patras, Greece, esokos@upatras.gr

2 Introduction Moment Tensors are important for studies like shakemap generation, tsunami warnings, ground motion evaluation and more Importance of automatic, quick and reliable Moment Tensor solution

3 Scisola Open-source python based software User friendly GUI Supports: Automatic Moment Tensor calculation of events provided by SeisComP3 in real-time Easy solution overview Quick solution revision Extensive configuration

4 SeisComP3 Likely the most widely distributed software package for seismological purposes Evolved within approximately the last 10 years Its use ranges from pure acquisition or real-time data exchange over Internet to a fully featured real-time earthquake monitoring Provides real-time waveform data through Seedlink protocol

5 ISOLA ISOLA moment tensor algorithm: Point source iterative deconvolution method Full wavefield is taken into account Discrete wavenumber method for local or regional distances Moment tensor is found by least-square minimization of misfit between observed and synthetic waveforms Position and time of centroid is optimized through grid search Used in a manual way in GI-NOA, UPSL and other institutes worldwide More info in Sokos & Zahradnik, 2008 and Sokos & Zahradnik, 2013

6 Architecture

7 SeisComP3 Interconnection

8 Automatic Procedure

9 Automatic Procedure

10 Origin Triggering Automatic Mode Watch SeisComP3 through scevtls Retrieve origin's info through scxmldump Manual Mode Execution through python script (mostly for testing purposes)

11 Automatic Procedure

12 Station Selection based on Distance Retrieves Stations/Streams Info from scisola database Filters Stations/Streams by certain type (e.g. HHN, HNE) Removes Blacklisted Stations/Streams defined by the user Calculates distance and azimuth of stations according to epicenter Selects Stations/Streams according to distance rule, defined by the user (e.g. 3.5 mw distance 100 km)

13 Automatic Procedure

14 Bad Station Data Filtering Filters unavailable stations/streams according to seedlink Retrieves records from seedlink in mseed format Filters stations/streams with gaps Filters stations/streams with clipping

15 Automatic Procedure

16 Station Data Correction Rotates stations/streams automatically Corrects stations/streams data Removes instrumental effect Aligns according to origin time Cuts to predefined duration Resamples

17 Automatic Procedure

18 Station Selection based on Azimuth Station distribution according to epicenter location 8 sectors, each of 45 Minimum number of sectors defined by the user Maximum stations per sector defined by the user

19 Automatic Procedure

20 Green's Functions Computation 1D crustal model defined by the user Computation of elementary seismograms (delta time function) Time window length of inversion procedure (tl) is defined by the user according to tl rule (e.g. 3.5 mw 4.5 tl = sec) Centroid horizontal position fixed at epicenter location while centroid depth is grid searched The number of trial sources above and below automatic depth estimation is defined by the user

21 Automatic Procedure

22 Inversion Computation Source inversion predefined to deviatoric type Inversion frequency band is defined by the user according to frequency rule (e.g. 3.5 mw 4.5 frequencies = [0.04, 0.05, 0.08, 0.09] Hz) Time grid search for centroid time is defined by the user

23 Automatic Procedure

24 Result Plotting Generates text files with results and the final focal mechanism Generates a map containing the epicenter location and contributing stations' location Generates a misfit of observed and synthetic waveforms plot Generates correlation plots containing focal mechanisms on different depth and time and the best solutions for the different depths Generates contributing streams plot Saves results to scisola database

25 Revise Procedure User can manually remove streams that have been selected by the automatic procedure from the inversion User can change inversion frequencies

26 Watcher

27 Automatic Vs Manual

28 Automatic Vs Manual Comparison of 46 manual MT solutions by GI-NOA, with 46 automatic by scisola Comparison using the Kagan angle metric Locations red dots

29 Automatic Vs Manual

30 Automatic Vs Manual 34/46 = similar solutions 74% success

31 Automatic Vs Manual

32 Automatic Vs Manual 39/46 = similar mw 85% success origins' difference not more than 0.3 of units

33 Automatic Vs Manual

34 Automatic Vs Manual 35/46 = similar c.depth 76% success

35 Automatic Vs Manual Results: Automatic Vs Manual 74% similar solutions Automatic MTs recognize the size and depth of the seismic source with adequate accuracy a few minutes after the event. Important for quick estimation of ground motions or tsunami hazard Average time of the automatic procedure is 5 minutes Quick revision, in just a few minutes, provides highly accurate solutions

36 Future Improvements Advanced methods and artificial intelligence techniques for stations/streams selection Implementation of multiple crustal models based on earthquake s location Advanced signal processing methods to be applied to seismic waveforms (to avoid various problems e.g. disturbances, noise and data transmission problems) Search of centroid in 3D grid surrounding hypocenter Pre-calculated Green s functions (for faster performance) Estimation of solution s quality

37 GUI

38 Log-in screen scisola DB (MySQL) SeisComP3 DB (MySQL/PostgreSQL)

39 Log-in screen databases scisola DB (MySQL) SeisComP3 DB (MySQL/PostgreSQL)

40 Main Screen

41 Main Screen main tabs

42 Main Screen main buttons

43 Main Screen latest 20 origins

44 Main Screen

45 Main Buttons start/stop watcher settings screen search screen update latest 20 origins about screen

46 Main Buttons start/stop watcher settings screen search screen update latest 20 origins about screen

47 Settings Screen

48 Settings Screen settings tabs

49 Settings Screen edit scisola stations

50 Settings Screen

51 Settings Screen

52 Settings Screen

53 Settings Screen

54 Settings Screen click to edit

55 Main Screen

56 Main Screen double click to review

57 Review Screen

58 Review Screen review tabs

59 Review Screen revision

60 Review Screen double click to remove stream revision

61 Review Screen delete origin

62 Review Screen

63 Review Screen scroll to zoom in/out

64 Review Screen

65 Review Screen

66 Review Screen

67 Review Screen

68 Review Screen

69 Review Screen

70 Main Buttons start/stop watcher settings screen search screen update latest 20 origins about screen

71 Search Screen shows origins in user defined datetime (GMT) range

72 Main Screen searched origins

73 Main Screen clicking on minimize

74 For More... source code and manual can be found here: