#obspy
These are public posts tagged with #obspy. You can interact with them if you have an account anywhere in the fediverse.
Magnitude 6.1 #earthquake at almost 100 km depth about 100 km East of Antofagasta, #Chile about an hour ago (2025-03-06Z16:21). Nice depth phases. Data from station C1.AC01 (Pan de Azucar, Chile, courtesy Red Sismologica Nacional and delivered via EarthScope using #obspy software) I will never take for granted data like this being available. I am not sure this will always be the case.
https://earthquake.usgs.gov/.../even.../us6000px89/executive
Magnitude 6.0 #earthquake about 25 minutes ago (2024-10-16 07:46:32) - signal on station IU.ANTO (Ankara, Tyrkie) - plot made using free #obspy software.
https://earthquake.usgs.gov/earthquakes/eventpage/us6000nyzk/executive
Magnitude 6.1 #earthquake in the #VolcanoIslands #Japan ~30 minutes ago (2024-02-12Z11:19) Signal on station IU.GUMO (Guam) Figure made using Google colab and #obspy software. Data obtained via #EarthScope
https://earthquake.usgs.gov/earthquakes/eventpage/us7000lym6/executive
The free #obspy python software allows you to make a "mock helicorder" plot of data from seismograph stations from almost anywhere around the world. This is today's datastream from station IU.DAV in the #Philippines. The first big signal is a magnitude 7.6 #earthquake. #aftershocks
![The free #obspy python software allows you to make a "mock helicorder" plot of data from seismograph stations from almost anywhere around the world. This is today's datastream from station IU.DAV in the #Philippines. The first big signal is a magnitude 7.6 #earthquake. #aftershocks
%matplotlib inline
import obspy
from obspy.clients.fdsn import Client
client = Client("IRIS")
from obspy import UTCDateTime
t = UTCDateTime("2023-12-02 12:00:00") #
[ sbef, saft ] = [ 0, 86400 ]
st = client.get_waveforms( "IU", "DAV", "10", "BHZ", t - sbef, t + saft)
st.filter("bandpass", freqmin=1.0, freqmax=4.0, corners=2, zerophase=False)
st.plot(type="dayplot", interval=15, right_vertical_labels=False,
vertical_scaling_range=1.0e5, one_tick_per_line=False,
color=['k', 'c', 'b', 'm'], show_y_UTC_label=True )](https://cdn.masto.host/oslotown/media_attachments/files/111/513/043/896/394/019/original/4390459108efa4c5.png "The free #obspy python software allows you to make a \"mock helicorder\" plot of data from seismograph stations from almost anywhere around the world. This is today's datastream from station IU.DAV in the #Philippines. The first big signal is a magnitude 7.6 #earthquake. #aftershocks
%matplotlib inline
import obspy
from obspy.clients.fdsn import Client
client = Client(\"IRIS\")
from obspy import UTCDateTime
t = UTCDateTime(\"2023-12-02 12:00:00\") #
[ sbef, saft ] = [ 0, 86400 ]
st = client.get_waveforms( \"IU\", \"DAV\", \"10\", \"BHZ\", t - sbef, t + saft)
st.filter(\"bandpass\", freqmin=1.0, freqmax=4.0, corners=2, zerophase=False)
st.plot(type=\"dayplot\", interval=15, right_vertical_labels=False,
vertical_scaling_range=1.0e5, one_tick_per_line=False,
color=['k', 'c', 'b', 'm'], show_y_UTC_label=True )")
Magnitude 6.4 #earthquake in the #Andreanof Islands in the #aleutianislands #Alaska almost exactly one hour ago (2023-10-16Z11:35) Signal on station IU.ADK in the Aleutians obtained from #EarthScope using free #Obspy software.
Quake should be harmless ...
https://earthquake.usgs.gov/earthquakes/eventpage/us6000lfvn/executive
Magnitude 7.2 #earthquake south of the #Alaska Peninsula 45 minutes ago (2023-07-16Z06:48) Signal on station II.KDAK (Kodiak Island, Alaska) made available via #Earthscope and read and plotted using free #obspy software.
https://earthquake.usgs.gov/earthquakes/eventpage/us7000kg30/executive
Magnitude 4.4 #earthquake off the West coast of #Greenland a couple of hours ago (2023-03-21Z06:56:56)
Signal on station IU.SFJD (Søndre Strømfjord, Greenland, made available via #earthscope and viewed using free #obspy software)
High frequency signal so tectonic. Classical Pn, Pg, Sn, Lg sequence for regional signals (about 300 km away)
https://earthquake.usgs.gov/earthquakes/eventpage/us7000jllc/executive
Earthquakes: Cross-correlation detector in #ObsPy is used to match a template from one seismic event to a longer trace from a seismometer, returning events which correlate above a user-defined threshold. Applying this to
#RaspiShake RD884, in Constantine, Cornwall, UK at a threshold of 0.80, I have identified 33 separate events in a swarm of multiplets (or repeaters) from the last week, identical earthquakes produced repeatedly from the same source. 12 of the events from this swarm have also been reported by #BGSSeismology. The largest event has a magnitude of 1.7, so they are all very small, although the M1.7 event was felt at the surface. The #ObsPy tutorial I followed is here: https://docs.obspy.org/master/tutorial/code_snippets/xcorr_detector.html and
#BGSSeismology report the swarm here: https://earthquakes.bgs.ac.uk/earthquakes/recent_uk_events.html.
Magnitude 5.2 #earthquake ca 30 minutes ago just west of #Bali #Indonesia (2022-12-13Z10:38) https://earthquake.usgs.gov/earthquakes/eventpage/us6000j8z3/executive
Signal on station GE.JAGI (Jagag, Java, Indonesia) courtesy of @geofon and processed using free #obspy software from @obspy
A M5.8 #Earthquake 14 km NW of Ciranjang-hilir, Indonesia at 2022-12-08 00:50:58UTC recorded at 109.3° away in SW England and Brittany on #Raspishake and #BGSseismology networks, in the core shadow zone. See: https://earthquake.usgs.gov/earthquakes/eventpage/us6000j7r0/executive. Uses #obspy, #matplotlib & #folium libraries.
Noticias sobre Python científico de la semana, episodio 56 
En resumen: Nuevas versiones de ObsPy y Polars, Python para la industria aseguradora, 40 métodos de detección de anomalías, vuelve la JupyterCon, y oferta Black Friday https://astrojuanlu.substack.com/p/episodio-56 Apoya el noticiero suscribiéndote por correo 
#python #scipy #numpy #jupyter #obspy #pyod #polars #pydata #jupytercon
Magnitude 5.4 #earthquake in the #SolomonIslands It is a bit offshore so hopefully not dangerous - too small to generate any significant tsunami.
Signal at the very close station IU.HNR (about 60 km away) obtained from #IRIS using #obspy software (@obspy)
There are several smaller aftershocks if you look carefully at a longer data segment. I will leave that as an exercise for you!
https://earthquake.usgs.gov/earthquakes/eventpage/us7000is0z/executive
Magnitude 4.5 #earthquake in the remote far North East of #Canada about 90 minutes ago (2022-11-19Z10:24) Signal at station CN.POIN (Pond Inlet, Canada) Data courtesy of #NaturalResourcesCanada and obtained from #IRIS using #obspy software ... only the lines of python you see.
High regional signal with regional P and S phases (about 100 km distance).
https://earthquake.usgs.gov/earthquakes/eventpage/us7000iqz2/origin/phase
I should have acknowledged that the plot is a time frequency representation created using the plot_tfr function in #ObsPy.
@obspy@fosstodon.orgWe're finalizing a new #obspy feature release version 1.4.0, which will get tagged during the next couple of days.
If anybody wants to test it, release candidate is on test.pypi:
ObsPy - a Python framework for seismological observatories.
test.pypi.orgTwo tiny #Earthquakes reported in #Cornwall, UK, on 2022-10-25 in Penpol on the Fal estuary and 2022-10-26 near Portreath on the north coast. Responses were seen on the local #Raspishake network. See: https://earthquakes.bgs.ac.uk/earthquakes/recent_events/20221025184145.html#page=summary & https://earthquakes.bgs.ac.uk/earthquakes/recent_events/20221026192120.html#page=summary. Uses #ObsPy & #MatPlotLib.
M5.2 #Earthquake 23 km ENE of Dipayal, Nepal at 2022-11-12 14:27:05UTC recorded on the Nepali Schools #RaspberryShake network in Nepal. See: https://earthquake.usgs.gov/earthquakes/eventpage/us7000ipbc/executive. Uses #ObsPy, #MatPlotLib & #GoogleEarthPro.
Magnitude 5.3 #earthquake in #Nepal almost exactly one hour ago (2022-11-12Z14:27). It comes four days after a 5.6 which caused significant damage and several fatalities. It is a reminder that earthquake hazard is not dominated by the very largest earthquakes - they are few and far between. The most lives are taken by "medium sized" quakes. Data obtained from #IRIS using #obspy
https://earthquake.usgs.gov/earthquakes/eventpage/us7000ipbc/executive
Responses to today's M7 earthquake from the Fiji region at 2022-11-12 07:09:14UTC recorded in SW England and Brittany, France on the
#raspishake and #BGSseismology networks. See: https://earthquake.usgs.gov/earthquakes/eventpage/us7000ip8h/executive. Uses #obspy, #matplotlib & folium libraries.
