Seismology Software

Seismic Attenuation
At Boston University (Geoff Abers) we developed attenuation spectral fitting routines using the Matlab-Antelope interface . These matlab scripts obtain origin and station information from an Antelope database, import waveforms and fit multi-taper spectra for seismic moment, corner frequency, and attenuation operator (t*). Currently, these scripts require BOTH P AND S picks for a particular earthquake/station combination. This is in the process of being fixed along with the ability to auto-window the frequency band of the spectrum based on a preset SNR.

The short story: This suite of matlab scripts, one tcl script, and a couple fortran routines determines attenuation operators (t*), moments, and corner frequencies from an Antelope database of waveforms, origins, and arrivals. Using the developed schema (atten_schema), the attenuation results are saved to an Antelope-style database. Spectra are also saved as text files (for now) and are accessed through the specdisc table.
Progress is being made to convert entirely to matlab scripts, and then port all routines to be accessed through a Tcl/Tk interface.

Modifications to Antelope dbevents: aeic_dbevents.tar (6.8Mb)
Updated: June 18, 2004

Notable changes:
-uncommented the parts to show dbpick

-maxage command line argument for number of days of events to display (from niko)

-added ability to play sound files when a new event comes in (when database is updated). this only happens when the waveforms are being displayed. I generated some simple sound files (.au) for numbers, seismic region numbers (SRN), and geographic region numbers (GRN). GRN only for SRN=1 (alaska). others can easily be generated from here , but extra logic needs to be added. Sound files are kept in directory 'soundpath' as designated in parameter file. Also, an "alert" sound is played before the earthquake location text, as designated by 'soundfile' in pf. Variables 'alerton' for audible alert and 'repeaton' for repeated audible alert (although this isn't operational yet-idea is to have popup window that needs to be closed before audible sound is stopped), both in pf. Added "Alert" pulldown menu in menubar to control these parameters on the fly.

-QDDS control. We needed a way to perform some quality control on events sent to QDDS(ANSS recenteqs webpage). Added another line when you right-click on event (on map) to run a script that deletes an event from the ANSS website, only deleting from the web site, not from the database. For our operations, the database opened with dbevents must be the same database that sends events to QDDS because we generate QDDS event ids from evid, which can differ between our four acquisition systems. "qddson" in pf allows for user to perform this function. Also, "Control" pulldown in menubar controls this variable. Bottom three lines in "Control" menu don't do anything yet. Not sure how useful all this is outside of AEIC.

-Latest feature: click on event and display detailed event map (probably the same one as generated by dbrecenteqs) with xv. Rt-click on event, "Show Detail Map" in pulldown menu. Pf "webmapsdb" points to database with webmaps table. Pf "xvgeometry" places image on screen in desired spot. Caveat: evids must be the same in webmapsdb as in dbname for events. Also, i run into trouble when i point dbevents at an origin table managed by dbrecenteqs.

Things i've thought about adding, but haven't gotten to it:
-when "Show Stations", be able to click on station and get station information (sta, lat, lon, etc.)
---> Done (6/18/2004)
-click on event in list and get same information and menu as when you right-click on event in map.
-add elog
HOME