Meteor in Russia – Seismic Signal ?

Very early this morning, a meteor lit up the skies of Russia, somewhere close to Ekaterinburg.


Same as for the Korean Boom Boom, we wanted to have a look at the seismic data of a close-by station to check if something was visible. The closest IRIS seismic station to this place is II.ARU (see location), at a distance of roughly 135 km.


One BHZ component is available on an STS-1 seismometer.To know the timestamp of the “event”, i took a guess that the overlay on the video was correct, and once corrected for the timezone, gave 03:26 UTC. Then, the easy part, getting data from IRIS :

from obspy.iris import Client
comp ="BHZ"
client = Client()
t1 = UTCDateTime("2013-02-15 03:00:00")
st = client.getWaveform("II", "ARU", "00", comp, t1 , t1 + 3600)

one gets:


or, high-pass filtered above 1 Hz:


There are two clear signals at 23 min and 38 min, but I have *no* idea if those are related to the meteor, its sonic boom(s), a potential impact with the ground, any-thing! If you have ideas, please share in the comments section below !

Finally, two close-ups of the higher amplitude stuff:

figure_3 figure_4


A M5.7 earthquake occurred in Tonga at 03:02:21.1, using obspy’s implementation of TauP, I get:

from obspy.taup.taup import getTravelTimes
from obspy.core import UTCDateTime
import pyproj
slat =56.43
slon =58.56

t0 = UTCDateTime("2013-02-15 03:02:21.1")
elat = -19.77
elon = 174.07
edep = 77

gc = pyproj.Geod(ellps='WGS84')
az1,az2,d2 = gc.inv(elon,elat,slon,slat)
km2deg = 40075.016686 / 360. #This is WGS84 norm
d = d2/(km2deg*1e3)
tt = getTravelTimes(d, edep, model='ak135')
if tt:
    for tti in tt:
        print tti['phase_name'],t0+tti['time']

arrival times:

     Pdiff 2013-02-15T03:17:27.424402Z
    pPdiff 2013-02-15T03:17:48.467493Z
    sPdiff 2013-02-15T03:17:56.730005Z
     PKPdf 2013-02-15T03:21:03.121484Z
     PKiKP 2013-02-15T03:21:03.211206Z
    pPKPdf 2013-02-15T03:21:24.910425Z
    pPKiKP 2013-02-15T03:21:24.993921Z
    sPKPdf 2013-02-15T03:21:33.009302Z
    sPKiKP 2013-02-15T03:21:33.094141Z
        PP 2013-02-15T03:22:29.301294Z
     SKPdf 2013-02-15T03:24:31.683740Z
     SKiKP 2013-02-15T03:24:32.054101Z
     PKSdf 2013-02-15T03:24:39.782373Z
     SKSac 2013-02-15T03:27:53.448022Z
     SKSdf 2013-02-15T03:28:08.284570Z
    pSKSac 2013-02-15T03:28:23.035425Z
    sSKSac 2013-02-15T03:28:31.218286Z
    pSKSdf 2013-02-15T03:28:38.175683Z
    sSKSdf 2013-02-15T03:28:46.273706Z
    SKKSac 2013-02-15T03:29:22.002222Z
     Sdiff 2013-02-15T03:30:14.318506Z
    pSdiff 2013-02-15T03:30:41.697168Z
    sSdiff 2013-02-15T03:30:50.556787Z
    PKKPbc 2013-02-15T03:31:10.561548Z
    PKKPab 2013-02-15T03:31:12.351709Z
    PKKPdf 2013-02-15T03:31:42.770532Z
        SP 2013-02-15T03:32:11.134790Z
        PS 2013-02-15T03:32:20.137964Z
    SKKPbc 2013-02-15T03:34:51.262353Z
    PKKSbc 2013-02-15T03:34:59.426904Z
    SKKPdf 2013-02-15T03:35:11.194726Z
    PKKSdf 2013-02-15T03:35:19.292139Z
    SKKSac 2013-02-15T03:38:36.127588Z
    SKKSdf 2013-02-15T03:38:47.644678Z
        SS 2013-02-15T03:38:50.680566Z
    P'P'df 2013-02-15T03:39:54.882959Z
    S'S'ac 2013-02-15T03:53:42.594873Z
    S'S'df 2013-02-15T03:54:13.323633Z

and, plotting those:


doesn’t explain all the peaks visible… Ideas, guys ?


2 thoughts on “Meteor in Russia – Seismic Signal ?”

  1. Cool stuff! Nice how one can get so easily real measurement data! Did you get some feedback on your 2 events in the meantime?

Leave a Reply

Your email address will not be published. Required fields are marked *