Missing info in event.pkl for events before 2011

[okling]

Hi Kasra,

I encountered another issue when I was updating my dataset. I'm using the following example command:

obspyDMT --datapath 2010-2019_EQ --min_date 2010-06-16 --max_date 2010-06-18 --net="CH" --station_rect="1/22/40/52" --loc "*" --cha "HHZ" --event_rect="170/180/-40/-30" --min_mag 6 --preset 0 --offset 3600 --sampling_rate 20 --instrument_correction --corr_unit "DIS" --pre_filt "0.01, 0.02, 10, 20" --event_catalog "NEIC_USGS" --data_source all --syngine --syngine_bg_model iasp91_2s --req_parallel 

However, when I read the event pickle file, some of the important information such as focal mechanism was missing:

>>> import obspy
>>> import pickle
>>> import os
>>> ev_load = open("/2010-2019_EQ/20100617_130646.a/info/event.pkl","rb")
>>> ev_pkl = pickle.load(ev_load)
>>> ev_load.close()
>>> print (ev_pkl)
OrderedDict([('number', 1), ('latitude', -33.168), ('longitude', 179.719), ('depth', 170.4), ('datetime', UTCDateTime(2010, 6, 17, 13, 6, 46, 590000)), ('magnitude', 6.0), ('magnitude_type', 'mwb'), ('author', 'None'), ('event_id', '20100617_130646.a'), ('origin_id', ResourceIdentifier(id="quakeml:earthquake.usgs.gov/archive/product/origin/usp000he6x/us/1415324502367/product.xml")), ('focal_mechanism', False), ('source_duration', False), ('flynn_region', 'NAN'), ('t1', UTCDateTime(2010, 6, 17, 13, 6, 46, 590000)), ('t2', UTCDateTime(2010, 6, 17, 14, 6, 46, 590000))])

Hence, the synthetics retrieved through Syngine using obspyDMT don't seem to be correct. Is this related to how obspyDMT read the NEIC product?

Thanks again for your help.

Cheers,
Angel

[kasra-hosseini]

Hi Angel, I ran your command in issue #55 , and it seems that the focal mechanism is retrieved correctly:

import obspy
import pickle
import os
ev_pkl = pickle.load(open("./Alaska_M6.5/EVENTS-INFO/event_list_pickle", "rb"))

Output:

[OrderedDict([('number', 1),
              ('latitude', 51.4234),
              ('longitude', -178.0262),
              ('depth', 33.93),
              ('datetime', 2018-08-15T21:56:56.310000Z),
              ('magnitude', 6.5),
              ('magnitude_type', 'mww'),
              ('author', 'None'),
              ('event_id', '20180815_215656.a'),
              ('origin_id',
               quakeml:earthquake.usgs.gov/archive/product/origin/us1000ga0z/us/1540657323040/product.xml),
              ('focal_mechanism',
               [6.1296e+18,
                -6.0987e+18,
                -3.09e+16,
                4.0674e+18,
                3.3288e+18,
                -2.0849e+18]),
              ('source_duration', False),
              ('flynn_region', 'NAN'),
              ('t1', 2018-08-15T21:56:56.310000Z),
              ('t2', 2018-08-15T23:56:56.310000Z)])]

[kasra-hosseini]

I also get the same results using:

ev_pkl = pickle.load(open("./Alaska_M6.5/20180815_215656.a/info/event.pkl", "rb"))

OUTPUT:

OrderedDict([('number', 1),
             ('latitude', 51.4234),
             ('longitude', -178.0262),
             ('depth', 33.93),
             ('datetime', 2018-08-15T21:56:56.310000Z),
             ('magnitude', 6.5),
             ('magnitude_type', 'mww'),
             ('author', 'None'),
             ('event_id', '20180815_215656.a'),
             ('origin_id',
              quakeml:earthquake.usgs.gov/archive/product/origin/us1000ga0z/us/1540657323040/product.xml),
             ('focal_mechanism',
              [6.1296e+18,
               -6.0987e+18,
               -3.09e+16,
               4.0674e+18,
               3.3288e+18,
               -2.0849e+18]),
             ('source_duration', False),
             ('flynn_region', 'NAN'),
             ('t1', 2018-08-15T21:56:56.310000Z),
             ('t2', 2018-08-15T23:56:56.310000Z)])

[okling]

Hi Kasra, thank you for checking. But the event you show here is in 2018. Could you try my example command or an event before 2011? The info is still missing after updating obspyDMT.