TCS AH provides access to software dealing with specified scientific problems, designed in the area of anthropogenic hazards. All of these services consist of software packages and visualization services, both applied to selected data subsets that a user has pre-defined and uploaded to the personal workspace. The users have therefore the possibility to proceed to either episode-oriented research (application of various methodologies focusing on a particular episode), or method-oriented research (perform a certain methodology at several datasets in order to compare results with each other and test the efficiency of selected approaches).

Services to be implemented are grouped within six blocks:

  • Basic services for data integration and handling;
  • Services for physical models of stress/strain changes over time and space as driven by geo-resource production;
  • Services for analysing geophysical signals;
  • Services to extract the relation between technological operations and observed induced seismic/deformation;
  • Services to quantitative probabilistic assessments of anthropogenic seismic hazard  - statistical properties of anthropogenic seismic series and their dependence on time-varying anthropogenesis; ground motion prediction equations; stationary and time-dependent probabilistic seismic hazard estimates, related to time-changeable technological factors inducing the seismic process;
  • Simulator for Multi-hazard/multi-risk assessment in ExploRation/exploitation of GEoResources (MERGER) - numerical estimate of the occurrence probability of chains of events or processes impacting the environment.

TCS AH will also serve the public sector expert knowledge and background information. In order to fulfill this aim the services for outreach, dissemination & communication will be implemented.

The following applications are already implemented:

Anderson-Darling test for exponentiality of inter-event time

This Application performs the Anderson-Darling test to study the null hypothesis that the inter-event times of a given series of events are drawn from an exponential distribution, alternatively that they are drawn from a Weibull distribution.


The autocorrelation (serial correlation) is a tool to find repeating patterns in a signal. Autocorrelation of a random process is the correlation between values of the process at different times, as a function of the two times or of the time lag.

Coefficient of Randomness

A statistical test to express changes in seismicity patterns by means of an objective criterion. The Matsumura Test provides the parameter of Matsumura, value of the parameter in 5% and 95% percentile and the patterns classification as regular, completely random, and clustered.

Cross Correlation

The Cross Correlation is a measure of similarity of two data series as a function of the displacement of one relative to the other.

Effective Stress Drop Estimate

Estimate of the scaling between total seismic moment and areal extent of earthquake cluster, which depends on the average stress drop released during the seismic sequence.

Estimation of source parameters in time-varying production parameters geometry

To estimate time-dependent values of seismic hazard parameters: the activity rate, the Gutenberg-Richter b-value, the return period and the exceedance probability for a prescribed area. The parameter values for the time moment t are estimated optionally either from events that occurred in dt time units preceding t, where dt is kept constant, or from n last events before t, where n is kept constant. Four magnitude distribution estimation methods are supported: maximum likelihood using the unbounded Gutenberg-Richter relation based model, maximum likelihood using the upper-bounded Gutenberg-Richter relation based model, unbounded non-parametric kernel estimation, upper-bounded non-parametric kernel estimation. The upper limit of magnitude distribution is evaluated using the Kijko-Sellevol generic formula.

FOCI application

This application makes use the P-wave first arrivals from the waveforms recorded by the seismic stations and resolves seismic moment tensors using various decomposition schemes. Three moment tensors are calculated: unconstrained (full) moment tensor solution, deviatoric moment tensor solution (no volumetric change in the source), and double-couple moment tensor solution (no volumetric change in the source, no linear dipole). In addition to the moment tensor components various source parameters are calculated as well (fault plane solutions, P/T/B axes orientation, scalar seismic moment etc).


This application implements TRMLOC localization/relocalization software which performs efficiently the inversion of seismic (acoustic) first arrival time onsets for hypocenter location using the probabilistic inverse theory approach. The application estimates the maximum likelihood hypocenter location, enables other hypocenter location estimators to be calculated and it allows an advanced analysis of location (inversion) uncertainties, providing in addition, meta-characteristics of the data and solutions.

Priestley-Subba Rao (PSR) Test

A test of stationarity to check, if statistical properties do not change with time. Stationarity in a strict sense is the strongest form of stationarity. It means that the joint statistical distribution of any collection of the time series variates never depends on time, the mean, variance and any moment of any variate is the same whichever variate. The second order stationarity is often consider for normal use. In this case a constant mean, a constant variance, and an autocovariance that does not depend on time. To check if a time series is stationary in second order mean a Priestley-Subba Rao (PSR) Test is used.

Spectral Analysis

P and S waves spectral levels and corner frequencies using Snoke alghorithm (Snoke, 1987).

Stationary Hazard: Exceedance Probability

This application uses the magnitude distribution and activity rate calculated from a seismic catalog to estimate the Exceedance Probability (EP). The EP of magnitude M in the time period of length T is the probability of an earthquake of magnitude M or greater to occur in T. The EP values are calculated optionally either for a fixed time period T and a series of magnitude values or for a fixed magnitude M and a series of period length values.

Stationary Hazard: Maximum Credible Magnitude

This application uses the magnitude distribution and activity rate from calculated a seismic catalog to estimate the Maximum Credible Magnitude (MCM). The MCM for the time period of length T is the magnitude value whose mean return period is T. The MCM values are calculated for a series of periods of length values.

Stationary Hazard: Mean Return Period

This application uses the magnitude distribution and activity rate calculated from a seismic catalog to estimate the Mean Return Period (MRP). The MRP of magnitude M is the average elapsed time between the consecutive earthquakes of magnitude M. The MRPs are calculated for a series of magnitude values.

Fracture Network Models - Mechanical Stresses

Interactively plot stress patterns associated with the hydraulic fracturing of a shale gas well. The data is based on the physical and geological conditions at the Preese Hall Episode site.

Front Advance Histograms

Application for drawing histograms and cumulative histograms of events count and seismic energy for selected bins, the width of which corresponds to a constant value of the mining front advance. Calculations are performed from the beginning of the excavation along the axis of the excavation.

Seismic Activity with Front Advance

Application allowing for visualization of the activity changes in time relatively to the front location

Source Localization
This application implements TRMLOC localization/relocalization software which performs efficiently the inversion of seismic (acoustic) first arrival time onsets for hypocenter location using the probabilistic inverse theory approach.

Moment Tensor Inversion
This application implements a command line application, focimt for estimating the moment tensor, which is a standard for description of earthquake kinematic source processes in the whole range of magnitudes.

Earthquake Spectra
This application determines seismic source parameters based on spectral analysis of the P and S waves. Two independent spectral parameters, i.e. low frequency spectral level and corner frequency, are estimated with  Snoke’s fitting algorithm. Assuming Brune's source model and spectral parameters.

Completeness Magnitude Estimation
The completeness magnitude of a given catalog, Mc, defined as the lowest magnitude at which all seismic events in a selected space-time volume are recorded by a network, is calculated by this service.

Source Size Distribution Functions | Stationary Seismic Hazard
This application implements several functions for estimating the probability density function and cumulative distribution function of magnitude.

Time-Dependent Seismic Hazard (in mining front surroundings)
This application comprises functions which are dedicated to episodes associated with mining induced seismicity. Its purpose is to evaluate seismic hazard parameters in connection with the evolution of mining operations and therefore to detect a causative relationship between seismic events and front advance.

Time-Dependent Seismic Hazard (in a selected area)
This application is a generalization of time-dependent seismic hazard, which can be performed for all episodes. Its purpose is to evaluate the temporal evolution of the seismic hazard parameters and therefore to detect a causative relationship between seismic events and production/technology activity.

Stress Inversion
This application implements the software package MSATSI which allows performing stress inversion using earthquake focal mechanisms. The purpose of this service is to determine the stress axes orientation and relative stress magnitude by inverting earthquake focal mechanisms.

Visualization include several generic tools implemented in the platform. The tools can be used to illustrate, e.g. the evolution of seismic activity together with the production/technological processes related to each episode. Various plotting options are provided, allowing the users to identify potential patterns and links between seismicity-technology. Depending on the inducing technology of the episode and the available data, the following additional visualization tools are offered: (i) visualization of mining front advance; (ii) seismic activity with water level, water volume, injection rate or wellhead pressure; (iii) injection rate and wellhead pressure histograms; (iv) focal mechanisms demonstration; (v) integrated episode data and 3D visualization.