Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic source”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Seismic velocity modeling and earthquake relocations for the southern Nevada National Security Site (NNSS), with a focus on Rock Valley

The Rock Valley Direct Comparison (RV/DC) is the third phase of the Source Physics Experiment (SPE; Snelson et al., 2014), a project aimed at improving seismic discrimination between explosive sources and natural earthquakes. Phases I and II of SPE focused on studying the generation of shear waves in endmember geologies (granite and alluvium, respectively), while Phase III aims to detonate a chemical explosion co-located with the source region of the 1993 Rock Valley earthquake sequence, located in the Rock Valley Fault Zone (RVFZ), in the southern portion of the Nevada Test Site (now the Nevada National Security Site, or NNSS).

58 GEOSCIENCES↗
DOE OSTI2024

Seismic Tomography 2024

Seismic tomography is the most abundant source of information about the internal structure of the Earth at scales ranging from a few meters to thousands of kilometers. It constrains the properties of active volcanoes, earthquake fault zones, deep reservoirs and storage sites, glaciers and ice sheets, or the entire globe. It contributes to outstanding societal problems related to natural hazards, resource exploration, underground storage, and many more. The recent advances in seismic tomography are being translated to nondestructive testing, medical ultrasound, and helioseismology. Nearly 50 yr after its first successful applications, this article offers a snapshot of modern seismic tomography. Focused on major challenges and particularly promising research directions, it is intended to guide both Earth science professionals and early-career scientists. The individual contributions by the coauthors provide diverse perspectives on topics that may at first seem disconnected but are closely tied together by a few coherent threads: multiparameter inversion for properties related to dynamic processes, data quality, and geographic coverage, uncertainty quantification that is useful for geologic interpretation, new formulations of tomographic inverse problems that address concrete geologic questions more directly, and the presentation and quantitative comparison of tomographic models. Finally, it remains to be seen which of these problems will be considered solved, solved to some extent, or practically unsolvable over the next decade.

58 GEOSCIENCES↗
DOE OSTI2025

Effect of Time Window and Spectral Measurement Options on Empirical Green’s Function Analysis Using DAS Array and Seismic Stations

The recorded seismic waveform is a convolution of event source term, path term, and station term. Removing high-frequency attenuation due to path effect is a challenging problem. Empirical Green’s function (EGF) method uses nearly collocated small earthquakes to correct the path and station terms for larger events recorded at the same station. However, this method is subject to variability due to many factors. Here, we focus on three events that were well recorded by the seismic network and a rapid response distributed acoustic sensing (DAS) array. Using a suite of high-quality EGF events, we assess the influence of time window, spectral measurement options, and types of data on the spectral ratio and relative source time function (RSTF) results. Increased number of tapers (from 2 to 16) tends to increase the measured corner frequency and reduce the source complexity. Extended long time window (e.g., 30 s) tends to produce larger variability of corner frequency. The multitaper algorithm that simultaneously optimizes both target and EGF spectra produces the most stable corner-frequency measurements. The stacked spectral ratio and RSTF from the DAS array are more stable than two nearby seismic stations, and are comparable to stacked results from the seismic network, suggesting that DAS array has strong potential in source characterization.

58 GEOSCIENCES↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗
DOE OSTI2024

SMART Deliverable 6.1.2a: Application of the ORION tool to the IBDP Carbon Storage Site

Forecasting and managing potential induced seismic activity is one of the challenges facing commercialscale geologic carbon sequestration (GCS), as well as other geologic energy extraction and byproduct disposal technologies. Historically, the process to develop robust, science-based forecasts of induced seismicity has required an integrated effort from experts in seismology, geomechanics, and reservoir engineering to manage data, develop and evaluate models of subsurface processes, and to calibrate and interpret the results from a range of models to understand site behavior relative to prescribed standards and in the context of uncertainty in geologic characterization data, forecasting models, and operational scenario uncertainty. The Operational Forecasting of Induced Seismicity (ORION) toolkit is an open-source, observation-based forecasting toolkit that is being co-developed by two U.S. DOE-funded initiatives: the National Risk Assessment Partnership (NRAP) and the Scienceinformed Machine Learning for Accelerating Real Time Decisions in Subsurface Applications (SMART) Initiative. ORION is designed to provide functionality to support decision making about seismic hazard analysis and risk management for GCS stakeholders ranging from the public to site operators to expert seismologists. The tool, which is written as open-source code in the Python programming language, is composed of a desktop graphical user interface (GUI) and an underlying forecasting engine. The forecasting engine uses available reservoir properties, well and fluid injection scenario details, and observed seismic catalog data as inputs to produce a set of temporal and spatio-temporal seismic forecasts.

58 GEOSCIENCES↗
DOE OSTI2026

Yield Estimates of three Historical Atmospheric Nuclear Explosions at Lop Nor, from Reduced Order Models of Regionally Recorded Rayleigh Waveforms and One Seismic Station: A Brief Communication

We use a reduced order model (ROM) for Rayleigh waveforms sourced by large atmospheric explosions along with data collected from a single seismic station (TLG) to estimate the yields of three atmospheric nuclear tests with order-megaton (MT) yields. These historical tests were conducted by China between 1973 and 1980 (CHIC 15, CHIC 16, and CHIC 26). We use our ROM to estimate yields of 2.34 MT (CHIC 15), 0.51 MT (CHIC 16), and 1.15 (CHIC 26) in the 20.5 s period band that are consistent with publicly accepted values of 3 MT (CHIC 15), 0.50 MT (CHIC 16), and 0.730 MT (CHIC 26).

58 GEOSCIENCES↗
DOE OSTI2021

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗
DOE OSTI2022

Analysis of the Source Physics Experiment Dry Alluvium Geology Ball Drop Sources

Two balls were dropped approximately 150 m down an open borehole as part of an investigation into the aftershocks from a 51 ton TNT-equivalent chemical explosion. The ball drop sources, a concrete-filled basketball and a bowling ball, were intended to mimic falling debris within an underground explosion borehole and cavity. Falling debris or sloughing of the borehole and cavity walls is thought to be one of the sources of explosion aftershocks. The signals from the ball drop sources were recorded by the same seismic network that monitored the aftershocks of the chemical explosion, so the ball drop recordings were used as templates to search through the continuous stream of aftershocks and detect correlated signals as a proportion of all aftershocks. The search resulted in only four similar event recordings out of the thousands of detected aftershocks. A synthetic test showed that the small, high-frequency ball drop templates were limited to producing highly correlated detections for events occurring in a small region around the impact location at 150 m depth. The high-frequency templates were necessary since the small ball drop signals were only above the noise at frequencies greater than 10 Hz. Unfortunately, the impact location is where the borehole was filled with sand in preparation for the subsequent chemical explosion in the test series, making it 150 m above the explosion created cavity centered at 300 m depth. This puts the ball drop templates 150 m away from the cavity, where many more events from falling debris would be expected. Further synthetic tests confirmed events emanating from the cavity 150 m below the impact site would have correlated with the lower frequency template and been detected had the template been usable below 10 Hz. Therefore, we conclude that the four detected events were most probably falling debris within the open borehole that impacted the floor at 150 m depth just as the ball drop sources had, and that this type of explosion aftershock mechanism comprises less than 1% of all explosion aftershocks.

58 GEOSCIENCES↗
DOE OSTI2024

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

58 GEOSCIENCES↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗
DOE OSTI2026

Evolution of storage monitoring – update in response to commercial and regulatory drivers

Carbon Capture and Storage (CCS) is in transition from first-of-a kind projects and research-orientated pilots to commercially-motivated applications. Monitoring results from many newly developed and planned large scale commercial projects are limited; however, it is worthwhile to assess their evolution and consider new strategies as part of an effort to assess and document best practices. Commercial monitoring is targeted to activities that comply with regulatory drivers and de-risk investments. Commercial monitoring also supports accounting that storage has occurred and is tied to project financing. It deals with long time frames and large volumes injected into multiple wells and multiple projects in favorable areas. We see developing trends toward reproducible workflows that systematically reduce risks and clarify expectations for oversight and long-term surveillance. Monitoring techniques showing increasing trends include injection zone pressure as a history-matching and compliance tool. To reduce cost and environmental impact of time-lapse seismic data collection, deploying new approaches and tools, such as use of fibre and installed sources are increasingly applied. Concern over the risk of induced seismicity by regulatory bodies and the general public has increased, which has also resulted in increased monitoring. Some techniques used in the early research phases have been sidelined or used only in restricted applications. For example, geochemical analyses in the injection zone as well as the environment are now being deployed less than it was in research-oriented programs, except in the US where it is required by the permitting process. Expectations of frequent area-wide near surface monitoring have also decreased.

25 ENERGY STORAGE↗
DOE OSTI2021

The First Detection of an Earthquake From a Balloon Using Its Acoustic Signature

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground-based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon-borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo-acoustic phenomenology demonstrates that seismic activity is detectable from a high-altitude platform on Earth, and that Rayleigh wave-induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus.

58 GEOSCIENCES↗
DOE OSTI2020

Play Fairway Analysis: Structurally Controlled Geothermal Systems in the Eastern Great Basin Extensional Regime, Utah

A research team with membership from the University of Utah/Energy & Geoscience Institute, the University of Utah/Dept. of Geology & Geophysics, and the Utah Geological Survey, undertook a play fairway analysis (PFA) for geothermal resources in the Eastern Great Basin (EGB) extensional tectonic regime of western Utah. This is a high-priority region for geothermal exploration because active Basin and Range (B&R) extension with volcanism having a N-S strike is superimposed upon pre-existing E-W belts of plutonic rocks and large-scale structural lineaments. Cumulative heat flow along the N-S strike of the state totals approximately 5 GWt above background stable interior. Three electricity producing power plants currently exist with substantial potential for increase. Succinctly, our PFA approach aims to resolve potential sources of heat and permeability in the region, which are the two principal criteria for establishing a geothermal resource. An initial Phase 1 was carried out using only existing geoscientific data in the area. Criteria selected for focusing heat potential include direct heat flow measurements in boreholes, magnetotelluric (MT) low resistivity anomalies, fluid/gas geochemistry, and proximity to recent volcanic eruptions. Permeability is established through geological structures (fault density, critically stressed areas, seismicity, gravity), and MT low resistivity anomalies. In Phase II of this PFA project, additional geological, geophysical and geochemical data were acquired and analysis carried out primarily over promising composite common risk (CCR) areas initially identified in Phase I in order to focus prospectivity and prepare for drilling recommendations. These prospects are near the Twin Peaks rhyolite field, high heat flow areas north of the producing Cove Fort system, and geophysical structure beneath the Crater Knoll area off the northeast flank of the Mineral Mountains. Additional data included MT site fill-in, structural mapping and analysis using high-resolution imagery, gravity and on-ground mapping, Nodal 3C passive seismic collection, and passive 3He surveying. Heat source and permeability potential are again expressed in terms of their individual common risk segment (CRS) maps, with a color scheme using green for most favorable (low risk) and red for least favorable (high risk). Diverse data types are united through the technique of probability kriging, which establishes prospectivity thresholds for each data type and then computes probability of exceeding that threshold over the PFA area. Modified CRS and CCRS maps are compared to those of Phase I to highlight tighter prospectivity focus. In doing so, the prospectivity threshold for heat was increased significantly to narrow the targeting. In the final Phase III of this project where a recommended deep thermal gradient hole was sited, additional geophysical, geochemical and geological field collection and analysis was carried out to refine drill hole targeting. This includes prospect-scale MT, gravity, structure, passive seismic deployment (Cove Fort area), and detailed 3He isotope profiling. It was the recommendation of the DOE Technical Monitoring Team (TMT) that one or more holes be sited in the north-ern Cove Fort area where legacy TG gradient holes showed high cumulative heat flow. These were to be of moderate depth, 2000-3000 feet, to reach the geothermal fluid table expected to start in excess of 1000 feet depth. The drilling organization stipulated by the DOE/GTO was that of the USGS Research Drilling Program (RDP) centered in Las Vegas, NV. A detailed well plan, appended to this report, was developed principally by Dr. Ben Barker consulting to University of Utah, Dr. Steve Pye on the DOE TMT, Mr. Steven Crawford of the USGS-RDP, and the project PI Phil Wannamaker. However, temperature and possible H 2 S at the systems lead to cancellation of the drilling last-minute as this appeared outside the experience base of the USGS-RDP. We hope to have the opportunity to revisit the test drilling and expand the Play Fairway Analysis of this region at some point in the future.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Wastewater Disposal Has Not Significantly Altered the Regional Stress State in Southern Kansas

Wastewater disposal is primarily responsible for the increased seismicity rate since ~2013 in southern Kansas. Previous work that used shear-wave splitting (SWS) in southern Kansas interpreted an ~90° temporal rotation in the fast polarization direction and attributed it to increased pore pressures resulting from fluid injection. However, this interpreted rotation coincided with a change in the stations used to make the SWS measurements. We investigate the temporal variability of fast azimuths in southern Kansas by making SWS measurements on earthquake families with similar source–receiver paths recorded on a stable local seismic network. We select high-quality SWS measurements by investigating the stability of results across 65 different frequency bands between 0.5 and 15 Hz. We find that the fast polarization direction in southern Kansas is relatively constant with an average east-northeast (~N79°E) orientation between 2014 and 2017. Our fast polarization measurements are primarily a reflection of the maximum principal horizontal stress direction (SHmax). We observe a slight spatial change in SHmax to the northeast (~N55°E) near the Nemaha ridge in Oklahoma. However, we do not observe any significant temporal rotation of SHmax or variation in delay time (i.e., crack density) in southern Kansas, contrary to the earlier study. The previously interpreted ~90° rotation may either be a reflection of a very local stress change or a misinterpretation of SWS results potentially due to the use of inconsistent source–receiver paths. Our SWS measurements cover the period of peak wastewater disposal and seismicity rates and suggest an absence of significant temporal rotations in the local anisotropy and stress orientations associated with wastewater disposal.

Geochemistry & Geophysics↗
DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
DOE OSTI2022

Report on the LLNL Global Seismic Waveform Tomography Modeling Project

Earth models have important applications for seismic event monitoring including the location and characterization of potential underground nuclear tests. Seismic waves generated by earthquakes and man-made events are altered by variations in Earth’s subsurface properties; and these distortions effect the arrival times of energy packets (phases) and obscure the characteristics of the original energy source (event mechanism). Three-dimensional models of Earth’s seismic properties can help explain and predict the distortions to the seismic wavefield and reveal the properties of a source that generated the waves, including the event location using the model-predicted timing of the phases and source mechanism using the modelpredicted waveform characteristics.

58 GEOSCIENCES↗
DOE OSTI2024

AmpSuite

Seismic amplitudes offer vital information about explosion source characteristics, including discrimination and yield estimation. To take advantage of this, we developed an interactive Python package to measure, control data quality, generate broad area propagation models and perform discrimination and estimate yield. Propagation models are essential in support of transportable yield and broad area discrimination. The key benefit of this package will be its ability to continuously integrate data and new techniques. The AmpSuite framework will provide standardized, repeatable, and accurate model generation and characterization routines. The capability is crucial for monitoring agencies tasked with rapid and high-quality seismic event characterization. The AmpSuite software includes a series of independent modules to perform: • Direct Phase Amplitude Measurement and Storage • Coda Envelope Measurement and Storage • Data Quality Control • New Propagation Model Developments • Seismic Discrimination and Analysis • Yield Estimation and supporting utility software. The AmpSuite software provides comprehensive solutions for monitoring agencies seeking to optimize model generation and event analysis within a contemporary Python framework. Stakeholders (AFTAC) have begun to move towards the Python language for scientific analysis as a new workforce emerges.

Alfaro, Richard↗