diff --git a/paper/paper.md b/paper/paper.md index 23ee52d..9437e36 100644 --- a/paper/paper.md +++ b/paper/paper.md @@ -39,22 +39,42 @@ The Statewide California Earthquake Center (SCEC) Ground Motion Simulation Valid Each code included in the GMSV Toolkit was redesigned to work in a standalone package so that its full capabilities are accessible for the validation of any simulated seismograms, without the need to install other software. For each code packaged in the GMSV Toolkit, we include both a command-line interface and a Python API. Each tool can therefore be run directly from the shell, allowing codes to be scripted together in a user-specified way, or with the Python API for use with other Python programs. We preserve full compatibility with the original codes from the BBP, and reproduce metric computation, statistics, and plot results obtained in previous BBP studies. -# Mathematics +# Time Series Comparisons -Single dollars ($) are required for inline mathematics e.g. $f(x) = e^{\pi/x}$ +The SCEC GMSV Toolkit provides a collection of seismogram processing, comparison, and plotting tools, such as: filtering, rotation, baseline correction, resampling, two, and three component comparison plots, etc. -Double dollars make self-standing equations: +Figure 1: Comparing two velocity timeseries (top left); side-by-side comparison of acceleration, velocity, and displacement timeseries (top-right); side-by-side acceleration, FAS, and PSA comparison (bottom). -$$\Theta(x) = \left\{\begin{array}{l} -0\textrm{ if } x < 0\cr -1\textrm{ else} -\end{array}\right.$$ +# Metric Computation -You can also use plain \LaTeX for equations -\begin{equation}\label{eq:fourier} -\hat f(\omega) = \int_{-\infty}^{\infty} f(x) e^{i\omega x} dx -\end{equation} -and refer to \autoref{eq:fourier} from text. +The SCEC GMSV Toolkit contains most validation metrics available in the Broadband Platform, including: +*Arias Intensity and Arias Intensity Duration +*RotD50 PSA +*Fourier Amplitude Spectra (FAS) +*RZZ 2015 metrics +*GMM computations for RotD50 PSA (NGA-West 1 and NGA-West 2 + +# PSA Validation + +Figure 2: In the PSA validation workflow, acceleration timeseries go through the RotDXX module where RotD50 is computed. This results is compared station by station with the plot_rotdxx module and aggregated across all stations and compared against a second data set with the PSA GoF tool. The PSA GoF Plot module generates different PSA GoF plots so users can see how two datasets match. + +Figure 3: PSA values for the two horizontal components and RotD50 (right). The green line shows calculated data and the blue line shows PSA values for the recorded seismogram. The two vertical lines (purple and red) come from the station list and indicate limits for which the recorded data is valid for this particular station (top). Calculated fourier amplitude spectra (FAS) for the two horizontal components, along with Smoothed Effective Amplitude Spectra (SEAS) (right). + +# Statistical Computation and Plot Generation + +Multiple tools available for combined calculated per-station metrics and generating multiple Goodness-of-Fit comparison plots: +*Pseudo-Spectral Acceleration (PSA) GoF +*Fourier Amplitude Spectra (FAS) GoF +*Anderson 2004 GoF +*GMM GoF comparisons + +Figure 4: Map PSA GoF plot (top left), showing per-station comparison of simulated against recorded data, useful for checking bias related to station location and directivity effects. Vs30 PSA GoF plot (bottom left), showing the same dataset plotted against each station Vs30 value. PSA GoF plot (right) showing how simulated data compares against recorded data at different periods. The solid red line shows the mean, the narrow band is the 90% confidence interval of the mean, and the wide band shows the standard deviation centered around the mean. Note that for periods between 1s and 5s the red line is close to zero, indicating the two sets are quite close. For periods < 1.0s (higher frequencies) the red line dips below zero, signaling the simulated data is overpredicting the recorded data. + +# GMSV Toolkit Summary + +*Makes well-verified and useful software available to a broad seismological and engineering community, leveraging over a decade of scientific, engineering, and software development by dozens of SCEC contributors +*Provides a standalone tool with a command-line interface and Python APIs +*Available on GitHub along with extensive documentation and examples # Citations