Using ObsPy to plot earthquake seismogram from HKO data

1. install PyWEED from anaconda command prompt. This will also install ObsPy

2. find a report  from HKO Earthquake

 utilising ObsPy , you can run the FDSN web services

3, run the following simple python script (wait for 10 seconds or maybe more)

example: 2022/01/21 16:00 ( UTC) 

my python script is uploaded to gitlab, only broadband sensors are selected.

btw, data is only available for date after 2009

note: 

network: HK

station: HKPS , intro from Dr Wong Wing Tak

instrument list of HKPS

three letters meaning；

1st letter : (response)

V = Very Long Period

H = High Broad Band

B = Broad Band

L = Long Period

2nd letter: (instrument)
C = ?

E = ?

H = High gain seismometer

K = ?

M = mass position seismometer

N =accelerometer

3rd letter:(orientation)
Z = Vertical
N = North-South
E = East-West

the time on HKO is in HKT(UTC+8), so you need to convert to UTC

Streckeisen STS-2/Quanterra 330HR Linear Phase: Broadband Sensor, response from 120 seconds to about 10 Hz

Episensor 200 Hz: Surface accelerometer

Steven Gibbons twitter

Krischer

Online URL builder example:

https://service.iris.edu/irisws/timeseries/1/query?net=HK&sta=HKPS&cha=BHZ&start=2022-09-09T05:08:00&end=2022-09-09T05:10:00&format=plot&loc=--

Amaseis (windows) sw for display the "worm" lines

OpendTect