/SeismoRMS

A simple Jupyter Notebook example for getting the RMS of a seismic signal (from PSDs)

Primary LanguageHTMLEuropean Union Public License 1.1EUPL-1.1

Ground Motion Displacement RMS vs Time

an example simple tutorial for getting seismic data, computing the power spectral densities, extracting the RMS and plotting

DOI

Required:

  • python
  • obspy (and its dependencies)
  • pandas
  • jupyter
  • notebook
  • tqdm

this should be easy to set up in a conda env: conda create -c conda-forge -n covid python=3.7 obspy pandas jupyter notebook tqdm

Author: Thomas Lecocq @seismotom, Fred Massin @fmassin

Run it interactively on mybinder.org

Original Example:

This was the output of the original code shared on Twitter end of March 2020. The following data shows the effect of the Social Distancing measures from the Belgian Government (2020-03-16 at midnight, and 2020-03-18 at midday):

Example image from this code:

Current Example:

The code has evolved and includes several new useful plots for interpreting the time series.

Standard ObsPy PPSD plots:

Standard PPSD:

Standard PPSD_spectrogram:

New plots:

RMS time series: RMS_timeseries:

Changes per day of week & time of day, before and after lockdown: RMS_daily_changes:

Changes per day of week & time of day, before and after lockdown, visualized as a 24-hour clock: RMS_hourly:

Colormapped clock plot, each record of the disc is 1 day: RMS_clockmap:

Colormapped grid plot, each column of the grid is 1 day: RMS_gridmap: