Open software

From EuroVolc Wiki
Jump to navigation Jump to search

Access to seismic full wavefield modelling tools (WP23)


VolcSeisSimu: Forward modelling of seismic waves in 2D & 3D heterogeneous media using the Elastic Lattice Method

The infrastructure comprises a software tool (ELM) for the numerical modelling of seismic wave propagation in highly heterogeneous media (including media with fractures). The tool kit also comprises starter heterogeneous volcano models. Whilst full wavefield simulators are becoming more common, a unique aspect of this tool is its ability to handle media with discrete fractures of arbitrary scale, in addition to intrinsic geological heterogeneity.

The Elastic Lattice Method (ELM) is a discrete particle numerical method for simulating seismic wave propagation through complex media in the presence of topography. The model structure consists of particles arranged on a cubic lattice which interact through a central force term and bond-bending force. Particle disturbances are propagated through space by numerically solving their equation of motion. A description of the code, examples, and bench marking against finite difference codes can be found in references [1-6]. The ELM code has both 2D and 3D implementations.

Potential application in volcano observatories

The aim of volcano monitoring is to understand physical processes within a volcano and ultimately help in forecasting potential hazards. One of the primary interests in volcano-seismology is determining volcanic seismic sources, but the generation and propagation of seismic waves throughout volcanic regions is a complex and nonlinear phenomenon controlled by the interaction of many processes. Because volcanoes are highly mechanically heterogeneous they significantly distort seismic wave propagation and hence the source ‘fingerprint’ in seismograms is often hidden by ‘path effects’ created when the wave propagates from to the source to the seismic receiver. Using the VolcSeisSimu tool researchers or monitoring scientists can determine these propagation path effects for arbitrarily heterogeneous models and ‘deconvolve’ them from observed seismograms to better understand and constrain volcano source models. Source inversion is a highly specialised and time consuming research area, but is critical to a better understanding of volcano seismograms in terms of source processes. This scheme offers an alternative 3D method for modelling wave propagation in the presence of strong topography and subsurface heterogeneity, and the flexibility of the method allows for a wide range of possible source mechanisms with no restriction on source geometric shape or on the distribution and degree of elastic heterogeneity.


The manual with information about how to compile and run the code is available here:

The seismic simulation codes are available to download here:


[1] O'Brien, Gareth S., Bean, Christopher J. (2004), A 3D discrete numerical elastic lattice method for seismic wave propagation in heterogeneous media with topography, Geophys. Res. Lett., Vol. 31, No. 14, L14608 doi:10.1029/2004GL020069.

[2] O'Brien, G. S., Bean, C. J. (2004), A discrete numerical method for modeling volcanic earthquake source mechanisms, J. Geophys. Res., Vol. 109, No. B9, B09301 doi:10.1029/2004JB003023.

[3] O'Brien, G. S. (2008), Discrete visco-elastic lattice methods for seismic wave propagation, Geophys. Res. Lett., 35, L02302, doi:10.1029/2007GL032214.

[4] O'Brien G.S., Bean C.J. and Tapamo H. (2009), Dispersion analysis and computational efficiency of Elastic Lattice methods for seismic wave propagation, Computers & Geosciences, Vol. 35, 1768-1775. doi:10.1016/j.cageo.2008.12.004.

[5] O'Brien, G. S. and C.J. Bean (2011), An irregular lattice method for elastic wave propagation, Geophysical Journal International, doi:10.1111/j.1365-246X.2011.05229.x.

[6] O'Brien, G.S., T. Nissen-Meyer and C.J. Bean (2012), A Lattice Boltzmann method for elastic wave propogation in a Poisson solid, Bulletin of the Seismological Soceity of America, VOLUME 103, No. 3, Pages 1224-1234, doi:10.1785/0120110191.

RETREAT - a REal-time TREmor Analysis Tool (WP9)


RETREAT (Smith & Bean, 2020) is a REal-time TREmor Analysis Tool written in python, making use of the ObsPy framework. It performs beamforming using frequency-wavenumber (f-k) analysis, or Least-squares beamforming, on real-time (or optionally archive) seismic array data to calculate the back azimuth and slowness values in a given time window, to aid in the location and analysis of volcanic tremor signals. While primarily intended as a tool for exploiting seismic array data to locate and track volcanic tremor, RETREAT also has the capability to analyze infrasonic array data to track acoustic sources. Release version 2.X now has the ability to process and analyse data from multiple arrays.


Full documentation on how to download, install and run the software is available here:

See also the accompanying publication Smith and Bean, 2020 for further description.


Eibl, Eva P. S., Bean, C.J., Vogfjörd, K.S., Ying, Y., Lokmer, I., Möllhoff, M., O’Brien, G.S., & Pálsson, F. (2017a). Tremor-rich shallow dyke formation followed by silent magma flow at Bárðarbunga in Iceland. Nature Geoscience volume 10, pages 299–304, doi:10.1038/ngeo2906

Eibl, E. P. S., Bean, C. J., Jónsdóttir, I., Höskuldsson, A., Thordarson, T., Coppola, D., Witt, T., and Walter, T. R. (2017b), Multiple coincident eruptive seismic tremor sources during the 2014–2015 eruption at Holuhraun, Iceland, J. Geophys. Res. Solid Earth, 122, 2972– 2987, doi:10.1002/2016JB013892.

Smith, P. J. and Bean, C. J., (2020), RETREAT: A REal-Time TREmor Analysis Tool for Seismic Arrays, With Applications for Volcano Monitoring. Front. Earth Sci. 8:586955. doi: doi:10.3389/feart.2020.586955

VolUnD - VOLcano UNrest Detection Toolkit (WP9)


VolUnD is a VOLcano UNrest Detection Toolkit written in Python 3 using the PyTorch library. It performs anomaly detection of volcanic historical data in an unsupervised way, leveraging the hypothesis that normal activity dominates the event distribution and can therefore be treated as the main source of information for training a model able to identify deviations as anomalies. These anomalies could be attributed to changes in the volcano state.


Further information and full documentation on how to download, install and run the software are available here:


Raffaele Mineo, Paolo Spadaro, Flavio Cannavò, Simone Palazzo, Andrea Cannata & Concetto Spampinato. (2021). EUROVOLC-ML/VolUnD-Toolkit. Zenodo. doi:10.5281/zenodo.5537253

WebOBS (WP10)


WebObs [1] is an integrated web-based system for data monitoring and networks management. Seismological and volcanological observatories have common needs and often common practical problems for multi disciplinary data monitoring applications. In fact, access to integrated data in real-time and estimation of uncertainties are keys for an efficient interpretation, but instruments variety, heterogeneity of data sampling and acquisition systems lead to difficulties that may hinder crisis management. In the Guadeloupe observatory, we have developed in the last 15 years an operational system that attempts to answer the questions in the context of a pluri-instrumental observatory. Based on a single computer server, open source scripts (with few free binaries) and a Web interface, the system proposes:

  • an extended database for networks management, stations and sensors (maps, station file with log history, technical characteristics, meta-data, photos and associated documents);
  • a web-form interfaces for manual data input/editing and export (like geochemical analysis, some of the deformation measurements, …); routine data processing with dedicated automatic scripts for each technique, production of validated data outputs, static graphs on preset moving time intervals, possible e-mail alarms, sensors and station status based on data validity;
  • in the special case of seismology, a multichannel continuous stripchart associated with EarthWorm/SeisComP acquisition chain, event classification database, automatic shakemap reports, regional catalog with associated hypocenter maps.
  • a statistical GNSS data modelling application [2] that provides simple Mogi and complex pCDM sources (whose application was highlighted in Eurovolc deliverable x.x)

WebObs is presently fully functional and used in a dozen observatories worldwide, but the documentation is mostly incomplete. We hope to shortly finish the main user’s manual. If you are in a hurry, please contact the project coordinator and we will be happy to help you to install it.


Further information and full documentation on how to download, install and run the software are available here:

Source code and binary packages repository are available here:


[1] Beauducel F., D. Lafon, X. Béguin, J.-M. Saurel, A. Bosson, D. Mallarino, P. Boissier, C. Brunet, A. Lemarchand, C. Anténor-Habazac, A. Nercessian, A. A. Fahmi (2020), WebObs: The volcano observatories missing link between research and real-time monitoring, Frontiers in Earth Sciences, doi:10.3389/feart.2020.00048

[2] Beauducel F., A. Peltier, A. Villié, W. Suryanto (2020), Mechanical imaging of a volcano plumbing system from unsupervised GNSS modeling, Geophysical Research Letters, doi:10.1029/2020GL089419