Skip to content

Releases: SeismicSource/sourcespec

Release v1.8

07 May 09:45
Compare
Choose a tag to compare

This long overdue release brings many improvements, new features and bugfixes
gradually introduced during the last year.

Release highlights:

  • New file formats for events and spectra
  • New configuration options to better specifify velocity and density models
  • Better support for P-wave inversion and teleseismic events
  • Radiation pattern correction from focal mechanism in radiated energy computation
  • Option for travel time-based signal window length
  • More options to control the calculation of source radius and stress drop
  • Improved estimation of radiated energy
  • New source parameter: apparent stress
  • Fix for map tiles not plotted anymore

This release requires at least Python 3.7.

Warning: the SQLite database used by this version is not compatible with
previous versions. You will need to upgrade your old database manually or using
source_spec -u DATABASE_FILE_NAME.

Check out the full Changelog

Release v1.7

31 Mar 14:46
Compare
Choose a tag to compare

This release improves trace processing through the use of modern routines for
instrument correction, optional baseline removal, a new clipping detection
algorithm and a better definition of signal and noise time windows.

A new plot is introduced, "stacked spectra" which allows to compare all the
spectra at once (and easily detect problematic stations 😉).
Also, a new command line tool, plot_sourcepars, allows making aggregate
plots of source parameters for many events (starting from the SQLite database).

New config file parameters have been added, while some have been removed.
Please run source_spec -U CONFIG_FILE_NAME to update your old config file.

As always, many bugfixes and improvements have been made in this release.
Thanks to all the users who took time to write and ask questions (by mail
or using the official SourceSpec [Discussions]).

Big kudos to Kris Vanneste [@krisvanneste] who helped all along the development
with code review and testing and submitted pull requests on noise windows and
clipping detection.

Check out the full Changelog.

Release v1.6

02 Aug 16:49
Compare
Choose a tag to compare

This release introduces several modifications to the config file.
You will need to upgrade your old config files manually or using
source_spec -U CONFIG_FILE_NAME.

This release requires at least Python 3.6.

A lot of effort has been devoted to improve the documentation.
Please check it out on https://sourcespec.readthedocs.io/

Input/output

  • QuakeML output (when using QuakeML input)
  • Command line option --run-id to provide a string identifying the
    current run (see pull request [#6])
  • Write SourceSpec version to parfile
  • Write SourceSpec version and run complete time to SQLite file
  • Write author and agency info (if specified) to output files, figures and
    HTML report
  • HTML page for misfit plots (when using grid search or importance sampling)
  • Station table in HTML report is now sortable (and its header remains fixed)!
  • Reduce PNG figures file size, while improving their resolution 😃
  • Removed option to read event information and traces from a pickle file
    (rarely used)

Processing

  • Support for P-wave spectral inversion (see pull request [#9])
  • It is now possible to provide different vp and vs velocities, close to the
    source and close to the stations (see the new config options above and
    issue [#5])
  • Possibility to choose a geometrical spreading model between
    (see issue [#8]):
    • rⁿ (default: n=1 – body waves)
    • Boatwright et al. (2002): "r" below a cutoff distance, frequency-dependent
      above the cutoff distance
  • Use travel time to compute quality factor from t* (and viceversa)
    (see issue [#5])
  • Compute travel time from pick and origin time, when possible
    (see issue [#10])
  • Warn if noise window ends after P or S arrival

Post-Inversion

  • Subtract the integral of noise spectrum from the integral of signal spectrum
    when computing radiated energy, under the hypothesis that energy is additive
    and noise is stationary

Config file

  • Config parameter paz has been removed and merged into station_metadata
  • Config parameters vp and vs have been renamed to vp_source
    and vs_source (see issue [#5])
  • New, optional, config parameter vp_stations and vs_stations
  • Config parameter pre_p_time and pre_s_time have been renamed to
    noise_pre_time and signal_pre_time, respectively (see pull request [#9])
  • Config parameter rps_from_focal_mechanism renamed to
    rp_from_focal_mechanism (see pull request [#9])
  • New config parameter: geom_spread_model (see issue [#8])
  • Config parameters PLOT_SHOW, PLOT_SAVE and PLOT_SAVE_FORMAT are now
    lowercase (plot_show, plot_save and plot_save_format)
  • New, optional, general config parameters for specyfing author and agency
    information. This information is written to output files and figures,
    if specified
  • New config parameter, event_url, to link the event page from the HTML
    report
  • Removed DEBUG config parameter
  • Parameters from GENERAL PARAMETERS section reorganized into a new
    section called TRACE AND METADATA PARAMETERS
  • Some parameters from INVERSION PARAMETERS moved into a new section
    called SPECTRAL MODEL PARAMETERS

Bugfixes

  • Fix for not working weighting options: frequency and no_weight
  • Fix for negative weights occasionally generated by interpolation
  • Fix bug when event coordinates are written into sqlite as binary blobs

Release v1.5

22 May 20:16
Compare
Choose a tag to compare

This is a pretty big release coming after several months of work on a large
dataset of more than 5000 events in Mayotte.

Please read through the changelog to discover all the improvements and new
features.

You will need to update your old config files via
source_spec -U CONFIG_FILE_NAME

Note that v1.5 is no more compatible with Python 2!

  • Input/output:
    • Write output files into a subdirectory of OUTDIR, whose name is
      the event id
    • Support for HYPOINVERSE-2000 output files
    • Removed autodetection of hypo71 file paths (specific to CRL case)
    • SQLite output: added radiated energy, weighted averages,
      errors on parameters, number of observations, hypocentral location
      and origin time
    • Removed -C argument to apply station correction to spectra. Now spectra
      are automatically corrected if residuals_filepath is specified in the
      configuration file
    • Save an additional event parameter to output files: average quality factor
    • Save additional station parameters to output files:
      source radius, Brune stress drop, source radius, quality factor
    • Mark outliers in .out file and in html report
    • Colored console output for log messages! (Not supported on Windows)
  • Processing:
    • New parameter for setting the width of the spectral smoothing window
      in terms of frequency decades: spectral_smooth_width_decades
    • Compute spectral weights after spectral correction (when a station
      residuals file is specified via residuals_filepath)
    • Removed configuration parameter trace_format
    • New configuration parameter sensitivity to provide a constant sensor
      sensitivity (flat response curve), which ovverides any response curve
      provided in metadata
    • New parameter for manually specifying trace units: trace_units
      (defaults to auto)
    • New approach for trace clipping detection (requires just one configuration
      parameter, named clip_max_percent)
      • Check for trace clipping only in the processing window
      • Use histogram of samples to detect clipping
    • Fix for wrong component used for 'SV' spectra (#3)
  • Inversion:
    • New config option: Mw_0_variability. Allowed variability around Mw_0
      during the main inversion. Previously hardcoded to 0.1
    • New inversion methods for grid sampling:
      • grid search (very slow!)
      • importance sampling of the misfit grid using k-d tree
        (faster, but less accurate)
    • Fix for Basin-hopping algorithm not running
  • Post-Inversion:
    • New set of post-inversion parameters to reject certain inversion results,
      per-station: pi_fc_min_max, pi_t_star_min_max, pi_bsd_min_max,
      pi_misfit_max
    • Reject inversion results when inverted fc is within 10% of fc_min or
      fc_max
    • Fix: use logarithmic error width for weighted logarithmic averages
      • previous way of computing weighted logarithmic averages was not correct!
    • Option to reject outliers using the IQR (interquartile range) method:
      parameter nIQR
    • Support for non symmetric error on station spectral parameters
    • Compute additional, per-station parameters: source radius, Brune stress
      drop and quality factor
    • Compute errors for all station parameters
    • Compute weighted averages for all event parameters (except radiated energy)
    • Compute spectral residuals using weighted average spectral parameters
  • Plotting:
    • Source parameter box plots to evaluate parameter dispersion across stations
      and visually detect outliers
    • Misfit plot (2D and 1D) when using grid sampling
    • cartopy removed as installation dependency, since it is not easily
      installable via pip
    • Use GSHHS database to draw coastlines.
      • New config option: plot_coastline_resolution
    • Correctly show circles on maps with diagonal smaller than 1 km
    • Fix plotting map colorbar on Matplotlib 3.5
    • Make average and errorbar lines more visible on map colorbar
    • Fix for error on plotting fc map, when only one station is available
    • Fix trace plot scaling for traces with larger signal outside the plot
      window
    • Do not plot 'H' spectrum if there is only one instrument component
      (since it will coincide with the only component)
    • Plot uncorrected spectrum when station correction is used

Release v1.4

15 Oct 11:44
Compare
Choose a tag to compare
  • New config option rps_from_focal_mechanism to compute station-specific
    S-wave radiation pattern from focal mechanism, if a focal mechanism is
    available in the QuakeML file
  • Plot the focal mechanism on maps, if it is available
  • Change default inversion algorithm to TNC (truncated
    Newton algorithm)
  • Config option dataless renamed to station_metadata
  • Config option traceids renamed to traceid_mapping_file
  • Config options ignore_stations and use_stations renamed to
    ignore_traceids and use_traceids, respectively
  • Support for 2D NonLinLoc grids (via nllgrid >= 1.4.1)
  • Possibility of using a generic DEFAULT NonLinLoc time grid
  • Added cartopy as an installation dependency
  • Fixed: nllgrid was always requested at runtime
  • Fixed: gracefully handle the case when there is no internet connection and
    map tiles cannot be downloaded
  • Fixed (Windows): suppress colored terminal output, which is not supported
  • Fixed (Windows): it is now possible to relaunch the same run, without
    having to delete the output directory first
  • Fixed (Windows): use same timezone names than Linux and macOS

Release v1.3

07 Oct 08:57
Compare
Choose a tag to compare
  • HTML reports
  • Option to provide station-specific spectral windowing

Release v1.2

07 Oct 08:59
Compare
Choose a tag to compare
  • Use python-versioneer to manage version numbers

Release v1.1

07 Oct 08:59
Compare
Choose a tag to compare
  • Bug fixes:
    • Accept band code C for broadband seismometers sampled at >=250 Hz
    • Require cartopy>=0.18 for compatibility with matplotlib>=3.3

Release v1.0

07 Oct 09:00
Compare
Choose a tag to compare
  • Simplification of time window parameters:
    • an unique window length, win_length, is used for time-domain S/N ratio
      calculation and spectral analysis
    • the old parameters, noise_win_length and s_win_length, are no more
      supported and generate an error
  • Reorganized config file and improved inline documentation within the file
    • removed unused option: fc_0
    • removed pre_filt option
    • post-processing check on fc bounds has been removed
  • Option to specify non standard instrument codes (e.g., L for
    acceleration)
  • Option to specify filter frequencies for local magnitude computation
  • Plotting improvements:
    • Show S/N ratio on trace plots
    • Show spectral S/N ratio on spectrum plots
    • Option to show/hide ignored traces/spectra (low S/N)
  • Bug fixes:
    • Pay attention to location code when building average spectra
    • Plotting: avoid overlapping traces with different location code
    • Plotting: avoid overlapping spectra with different location code

Release v0.9

07 Oct 09:01
Compare
Choose a tag to compare
  • Support for QuakeML input and StationXML
  • Support for Python 3.5
  • Only compatible with ObsPy >= 1.1.0
  • Project reorganization:
    • Project renamed to SourceSpec
    • ssp_residuals renamed to source_residuals
    • New installable package (e.g., via pip)
  • Spectra are smoothed in log-freq (no more Konno-Ohmachi)
  • Inversion is performed in a log-freq space
  • Option to invert for t_star_0 on the plateau level
  • Traces are filtered before computing S/N ratio
  • Trace clipping detection
  • Traces are always plotted, even if no inversion is performed
  • Use by default a global model for theoretical travel time calculation
  • Possibility of using NonLinLoc travel time grids (requires nllgrid)
  • New options for P and S arrival time tolerance
  • New option for maximum epicentral distance for trace processing
  • Possibility of using a NonLinLoc model grid for obtaining vs at the source
    and at the receiver (requires nllgrid)
  • Use log10 of weighting function, since the inversion is done in magnitude
    units
  • Use json format for traceid correction file
  • Save config file to output dir (only for source_spec)
  • Save run completion time into output file
  • Logarithmic average and logarithmic (asymmetric) error bars for Mo, fc and
    source radius
  • Computation of radiated energy
  • Station-specific filters
  • New parameters: gap_max, overlap_max
  • Add legend to spectral plots
  • Add event information, code version and run completion time to plots
  • Multifigure plotting for traces and spectra (for large number of stations)
  • New option to plot a station map, color-coded by station magnitude
    (requires Cartopy)
  • Refactoring of local magnitude computation code:
    • Wood-Anderson amplitude is computed from the whole trace converted to W-A
      and not converting only the min and max peaks (which is the default in
      ObsPy)
    • Trace windowing is computed on HF envelopes
    • New option for custom Ml coefficients
    • Remove unnecessary Ml options
  • Code cleaning and optimization:
    • Switch from optparse (deprecated) to argparse
    • Code style fixes and refactoring
  • New option to update a config file from previous versions
  • BUGFIX: Fix a major bug in reading hypo pick file
  • BUGFIX: Fix for pick file not read in source_model