Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Feasibility of Geophysical Monitoring of Carbon-Sequestrated Deep Saline Aquifers

As carbon dioxide (CO{sub 2}) is sequestered from the bottom of a brine reservoir and allowed to migrate upward, the effects of the relative permeability hysteresis due to capillary trapping and buoyancy driven migration tend to make the reservoir patchy saturated with different fluid phases over time. Seismically, such a patchy saturated reservoir induces an effective anisotropic behavior whose properties are primarily dictated by the nature of the saturation of different fluid phases in the pores and the elastic properties of the rock matrix. By combining reservoir flow simulation and modeling with seismic modeling, it is possible to derive these effective anisotropic properties, which, in turn, could be related to the saturation of CO{sub 2} within the reservoir volume any time during the post-injection scenario. Therefore, if time-lapse seismic data are available and could be inverted for the effective anisotropic properties of the reservoir, they, in combination with reservoir simulation could potentially predict the CO{sub 2} saturation directly from the time-lapse seismic data. It is therefore concluded that the time-lapse seismic data could be used to monitor the carbon sequestrated saline reservoirs. But for its successful implementation, seismic modeling and inversion methods must be integrated with the reservoir simulations. In addition, because CO{sub 2} sequestration induces an effective anisotropy in the sequestered reservoir and anisotropy is best detected using multicomponent seismic data compared to single component (P-wave) data, acquisition, processing, and analysis is multicomponent seismic data is recommended for these time-lapse studies. Finally, a successful implementation of using time-lapse seismic data for monitoring the carbon sequestrated saline reservoirs will require development of a robust methodology for inverting multicomponent seismic data for subsurface anisotropic properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Identification of Natural Oscillation Modes for Purposes of Seismic Assessment and Monitoring of HPP Dams

The paper puts forward a method for processing data from detailed seismic assessments of HPP dams (dynamic tests). A detailed assessment (hundreds of observation points in dam galleries) is performed with consideration of operating dam equipment and the microseismic noise. It is shown that dynamic oscillation characteristics (natural oscillation frequencies and modes in the main dam axes, the velocities of propagation of elastic waves with given polarization, and so on.) can be determined with sufficient accuracy by using complex transfer functions and pulse characteristics. Monitoring data is processed using data from a detailed assessment, taking account of identified natural oscillation modes and determined ranges of natural frequencies. The spectra of characteristic frequencies thus obtained are used to choose substitution models and estimate the elastic characteristics of the “dam – rock bed” construction system, viz., the modulus of elasticity (the Young modulus), the Poisson ratio, the dam section stiffness with respect to shear, tension and compression and the elastic characteristics of the rock foundation.

42 ENGINEERING↗
DOE OSTI2015

sRecovery Act: Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO 2 Storage

This study focuses on evaluating the feasibility and suitability of using the Jurassic/Triassic (J/TR) sediments of the South Georgia Rift basin (SGR) for CO 2 storage in southern South Carolina and southern Georgia The SGR basin in South Carolina (SC), prior to this project, was one of the least understood rift basin along the east coast of the U.S. In the SC part of the basin there was only one well (Norris Lightsey #1) the penetrated into J/TR. Because of the scarcity of data, a scaled approach used to evaluate the feasibility of storing CO 2 in the SGR basin. In the SGR basin, 240 km (~149 mi) of 2-D seismic and 2.6 km2 3-D (1 mi2) seismic data was collected, process, and interpreted in SC. In southern Georgia 81.3 km (~50.5 mi) consisting of two 2-D seismic lines were acquired, process, and interpreted. Seismic analysis revealed that the SGR basin in SC has had a very complex structural history resulting the J/TR section being highly faulted. The seismic data is southern Georgia suggest SGR basin has not gone through a complex structural history as the study area in SC. The project drilled one characterization borehole (Rizer # 1) in SC. The Rizer #1 was drilled but due to geologic problems, the project team was only able to drill to 1890 meters (6200 feet) instead of the proposed final depth 2744 meters (9002 feet). The drilling goals outlined in the original scope of work were not met. The project was only able to obtain 18 meters (59 feet) of conventional core and 106 rotary sidewall cores. All the conventional core and sidewall cores were in sandstone. We were unable to core any potential igneous caprock. Petrographic analysis of the conventional core and sidewall cores determined that the average porosity of the sedimentary material was 3.4% and the average permeability was 0.065 millidarcy. Compaction and diagenetic studies of the samples determined there would not be any porosity or permeability at depth in SC. In Georgia there appears to be porosity in the J/TR section based on neutron log porosity values. The only zones in Rizer #1 that appear to be porous were fractured diabase units where saline formation water was flowing into the borehole. Two geocellular models were created for the SC and GA study area. Flow simulation modeling was performed on the SC data set. The injection simulation used the newly acquired basin data as well as the Petrel 3-D geologic model that included geologic structure. Due to the new basin findings as a result of the newly acquired data, during phase two of the modeling the diabase unit was used as reservoir and the sandstone units were used as caprock. Conclusion are: 1) the SGR basin is composed of numerous sub-basins, 2) this study only looked at portions of two sub-basins, 3) in SC, 30 million tonnes of CO 2 can be injected into the diabase units if the fracture network is continuous through the units, 4) due to the severity of the faulting there is no way of assuring the injected CO 2 will not migrate upward into the overlying Coastal Plain aquifers, 5) in Georgia there appears to porous zones in the J/TR sandstones, 6) as in SC there is faulting in the sub-basin and the seismic suggest the faulting extends upward into the Coastal Plain making that area not suitable for CO 2 sequestration, 7) the complex faulting observed at both study areas appear to be associated with transfer fault zones (Heffner 2013), if sub-basins in the Georgia portion of the SGR basin can be located that are far away from the transfer fault zones there is a strong possibility of sequestering CO 2 in these areas, and 9) the SGR basin covers area in three states and this project only studied two small areas so there is enormous potential for CO 2 sequestration in other portions the basin and further research needs to be done to find these areas.

58 GEOSCIENCES↗
DOE OSTI2019

The initial appraisal of buried DAS system in CO2CRC Otway Project: the comparison of buried standard fibre-optic and helically wound cables using 2D imaging

This study aims to assess the ability of shallow distributed acoustic sensing (DAS) to serve as a cost-effective seismic sensor array for permanent monitoring applications. To this end, as part of the CO2CRC seismic monitoring program, a fibre-optic DAS array was deployed alongside a permanently buried geophone array at the Otway Project site (Victoria, Australia). The DAS array consisted of a standard commercially available tactical fibre-optic cable, which was deployed in 0.8m deep trenches. A custom-designed helically wound (HW) cable was also deployed in one of the DAS trenches for comparison of the cable designs. Simultaneous acquisition of the seismic data was carried out using 3000 vibroseis source points and geophones, DAS standard and HW cables. For initial assessment of the seismic images acquired with DAS and to compare different cable designs, preliminary 2D seismic reflection processing is conducted on both DAS cables and geophone data along a single 2D line. The geophone data processing guided processing of the DAS data. Several shallow structures (100-450ms) and some important reflectors at 450-600ms are observed on the final DAS images. Comparison of the two different DAS cable types demonstrated that seismic imaging would benefit DAS technology. However, the benefit of utilising HW cable is modest compared with the standard cable. The workflows and results of this study pave the way for processing of the 3D seismic data set acquired with the DAS array, as well as further detailed analysis of the DAS cables and the system itself.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

West Flank Coso FORGE: Seismic Reflection

PDFs of seismic reflection profiles 101,110, 111 local to the West Flank FORGE site. 45 line kilometers of seismic reflection data are processed data collected in 2001 through the use of vibroseis trucks. The initial analysis and interpretation of these data was performed by Unruh et al. (2001). Optim processed these data by inverting the P-wave first arrivals to create a 2-D velocity structure. Kirchhoff images were then created for each line using velocity tomograms (Unruh et al., 2001).

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

Magnetotelluric Data Collected in 2016 over the San Emidio Geothermal Field in Nevada

This data set includes the magnetotelluric (MT) data collected from October 21 to November 9, 2016 over the San Emidio geothermal field in Nevada by Quantec Geoscience USA Inc. on behalf of US Geothermal Inc. as part of a project entitled "A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography". This data set includes descriptions of the instrumentation, data acquisition and processing procedures, as well as the final processed data and digital archive formats. A total of 81 MT locations were surveyed (52 profile sites, and 29 MT sites). Data were processed and inspected for quality assurance on site, and reviewed daily by the geophysicist in charge of the project.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Estimation of anisotropy parameters in organic-rich shale: Rock physics forward modeling approach

Anisotropy analysis becomes an important step in processing and interpretation of seismic data. One of the most important things in anisotropy analysis is anisotropy parameter estimation which can be estimated using well data, core data or seismic data. In seismic data, anisotropy parameter calculation is generally based on velocity moveout analysis. However, the accuracy depends on data quality, available offset, and velocity moveout picking. Anisotropy estimation using seismic data is needed to obtain wide coverage of particular layer anisotropy. In anisotropic reservoir, analysis of anisotropy parameters also helps us to better understand the reservoir characteristics. Anisotropy parameters, especially ε, are related to rock property and lithology determination. Current research aims to estimate anisotropy parameter from seismic data and integrate well data with case study in potential shale gas reservoir. Due to complexity in organic-rich shale reservoir, extensive study from different disciplines is needed to understand the reservoir. Shale itself has intrinsic anisotropy caused by lamination of their formed minerals. In order to link rock physic with seismic response, it is necessary to build forward modeling in organic-rich shale. This paper focuses on studying relationship between reservoir properties such as clay content, porosity and total organic content with anisotropy. Organic content which defines prospectivity of shale gas can be considered as solid background or solid inclusion or both. From the forward modeling result, it is shown that organic matter presence increases anisotropy in shale. The relationships between total organic content and other seismic properties such as acoustic impedance and Vp/Vs are also presented.

58 GEOSCIENCES↗
DOE OSTI2014

2D Seismic Reflection Data across Central Illinois

In a continuing collaboration with the Midwest Geologic Sequestration Consortium (MGSC) on the Evaluation of the Carbon Sequestration Potential of the Cambro-Ordovician Strata of the Illinois and Michigan Basins project, Schlumberger Carbon Services and WesternGeco acquired two-dimensional (2D) seismic data in the Illinois Basin. This work included the design, acquisition and processing of approximately 125 miles of (2D) seismic reflection surveys running west to east in the central Illinois Basin. Schlumberger Carbon Services and WesternGeco oversaw the management of the field operations (including a pre-shoot planning, mobilization, acquisition and de-mobilization of the field personnel and equipment), procurement of the necessary permits to conduct the survey, post-shoot closure, processing of the raw data, and provided expert consultation as needed in the interpretation of the delivered product. Three 2D seismic lines were acquired across central Illinois during November and December 2010 and January 2011. Traversing the Illinois Basin, this 2D seismic survey was designed to image the stratigraphy of the Cambro-Ordovician sections and also to discern the basement topography. Prior to this survey, there were no regionally extensive 2D seismic data spanning this section of the Illinois Basin. Between the NW side of Morgan County and northwestern border of Douglas County, these seismic lines ran through very rural portions of the state. Starting in Morgan County, Line 101 was the longest at 93 miles in length and ended NE of Decatur, Illinois. Line 501 ran W-E from the Illinois Basin – Decatur Project (IBDP) site to northwestern Douglas County and was 25 miles in length. Line 601 was the shortest and ran N-S past the IBDP site and connected lines 101 and 501. All three lines are correlated to well logs at the IBDP site. Originally processed in 2011, the 2D seismic profiles exhibited a degradation of signal quality below ~400 millisecond (ms) which made interpretation of the Mt. Simon and Knox sections difficult. The data quality also gradually decreased moving westward across the state. To meet evolving project objectives, in 2012 the seismic data was re-processed using different techniques to enhance the signal quality thereby rendering a more coherent seismic profile for interpreters. It is believed that the seismic degradation could be caused by shallow natural gas deposits and Quaternary sediments (which include abandoned river and stream channels, former ponds, and swamps with peat deposits) that may have complicated or changed the seismic wavelet. Where previously limited by seismic coverage, the seismic profiles have provided valuable subsurface information across central Illinois. Some of the interpretations based on this survey included, but are not limited to: - Stratigraphy generally gently dips to the east from Morgan to Douglas County. - The Knox Supergroup roughly maintains its thickness. There is little evidence for faulting in the Knox. However, at least one resolvable fault penetrates the entire Knox section. - The Eau Claire Formation, the primary seal for the Mt. Simon Sandstone, appears to be continuous across the entire seismic profile. - The Mt. Simon Sandstone thins towards the western edge of the basin. As a result, the highly porous lowermost Mt. Simon section is absent in the western part of the state. - Overall basement dip is from west to east. - Basement topography shows evidence of basement highs with on-lapping patterns by Mt. Simon sediments. - There is evidence of faults within the lower Mt. Simon Sandstone and basement rock that are contemporaneous with Mt. Simon Sandstone deposition. These faults are not active and do not penetrate the Eau Claire Shale. It is believed that these faults are associated with a possible failed rifting event 750 to 560 million years ago during the breakup of the supercontinent Rodinia.

Smith, Valerie↗
DOE OSTI2015

IDC Reengineering Phase 2 Project Scope

The Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) operates the International Data Centre (IDC) to support treaty verification activities by Member States. The IDC collects, stores, processes, and distributes data from seismic, hydroacoustic, infrasound, and radionuclide stations. A significant component of the IDC is the seismic, hydroacoustic, and infrasound (SHI) processing system. The IDC has recognized the need to reengineer the SHI processing system and has defined a three-phase reengineering effort. IDC Reengineering Phase 1 (RP1) started in 2011 with the goal to significantly modernize SHI automatic and interactive data processing software, focusing on enhancements to individual components of the system. Reengineering Phase 2 (RP2), started in 2014, addresses the task of specifying a unified architecture for all SHI software, across processing stages, to pave the way for further software development and sustainment in the future. A key objective of RP2 is an architecture sufficient to provide a basis for a cost estimate for the development or enhancement of the software components and subsystems. This report provides a brief description of the scope of the IDC Reengineering Phase 2 project. It describes the goals and objectives of reengineering, the system definition, and the technical scope of the system.This report expands upon the “Scope for IDC Reengineering Phase 2”, July 2014, ECS-DIS-WGB-43-PTS-MATERIAL-12-ADD-1, produced by the IDC.

58 GEOSCIENCES↗
DOE OSTI2019

Frequency-Bessel Transform Method for Effective Imaging of Higher-Mode Rayleigh Dispersion Curves From Ambient Seismic Noise Data

It has been widely recognized that the cross-correlation function of ambient seismic noise data recorded at two stations approximates to the part of Greens function between two stations. Therefore, the cross-correlation function should include higher modes, aside from the fundamental mode. However, the problem of measuring or extracting overtones from ambient seismic noise data remains. In this paper, we propose the frequency-Bessel transform method (F-J method) for extracting the dispersion curves of higher modes from ambient seismic noise data. We then assess the validity, accuracy, and applicability of the F-J method by conducting extensive numerical simulations and processing the observed ambient seismic noise data of the USArray. As demonstrated in this study, the F-J method is a convenient, practical, and accurate method for extracting the dispersion curves of multimodes from ambient seismic noise data and therefore has significant potentiality in the field of ambient seismic noise tomography.

58 GEOSCIENCES↗
DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗
DOE OSTI2013

Rotation-Enabled 7-Degree of Freedom Seismometer for Geothermal Resource Development. Phase 1 Final Report

Under this Department of Energy (DOE) grant, A-Tech Corporation d.b.a. Applied Technology Associates (ATA), seeks to develop a seven-degree-of-freedom (7-DOF) seismic measurement tool for high-temperature geothermal applications. The Rotational-Enabled 7-DOF Seismometer includes a conventional tri-axial accelerometer, a conventional pressure sensor or hydrophone, and a tri-axial rotational sensor. The rotational sensing capability is novel, based upon ATA's innovative research in rotational sensing technologies. The geothermal industry requires tools for high-precision seismic monitoring of crack formation associated with Enhanced Geothermal System (EGS) stimulation activity. Currently, microseismic monitoring is conducted by deploying many seismic tools at different depth levels along a 'string' within drilled observation wells. Costs per string can be hundreds of thousands of dollars. Processing data from the spatial arrays of linear seismometers allows back-projection of seismic wave states. In contrast, a Rotational-Enabled 7-DOF Seismometer would simultaneously measure p-wave velocity, s-wave velocity, and incident seismic wave direction all from a single point measurement. In addition, the Rotational-Enabled 7-DOF Seismometer will, by its nature, separate p- and s-waves into different data streams, simplifying signal processing and facilitating analysis of seismic source signatures and geological characterization. By adding measurements of three additional degrees-of-freedom at each level and leveraging the information from this new seismic observable, it is likely that an equally accurate picture of subsurface seismic activity could be garnered with fewer levels per hole. The key cost savings would come from better siting of the well due to increased information content and a decrease in the number of confirmation wells drilled, also due to the increase in information per well. Improved seismic tools may also increase knowledge, understanding, and confidence, thus removing some current blocks to feasibility and significantly increasing access to potential geothermal sites. During the Phase 1 effort summarized in this final report, the ATA Team modeled and built two TRL 3 proof-of-concept test units for two competing rotational sensor technologies. The two competing technologies were based on ATA's angular rate and angular displacement measurement technologies; Angular rate: ATA's Magnetohydrodynamic Angular Rate Sensor (Seismic MHD); and Angular displacement: ATA's Low Frequency Improved Torsional Seismometer (LFITS). In order to down-select between these two technologies and formulate a go / no go decision, the ATA Team analyzed and traded scientific performance requirements and market constraints against sensor characteristics and components, acquiring field data where possible to validate the approach and publishing results from these studies of rotational technology capability. Based on the results of Phase 1, the ATA Team finds that the Seismic MHD (SMHD) technology is the best choice for enabling rotational seismometry and significant technical potential exists for micro-seismic monitoring using a downhole 7-DOF device based on the SMHD. Recent technical papers and field data confirm the potential of rotational sensing for seismic mapping, increasing confidence that cost-reduction benefits are achievable for EGS. However, the market for geothermal rotational sensing is small and undeveloped. As a result, this report recommends modifying the Phase 2 plan to focus on prototype development aimed at partnering with early adopters within the geothermal industry and the scientific research community. The highest public benefit will come from development and deployment of a science-grade SMHD rotational seismometer engineered for geothermal downhole conditions and an integrated test tool for downhole measurements at active geothermal test sites.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Utah FORGE - Development of a Reservoir Seismic Velocity Model and Seismic Resolution Study

This is data from and a final report on the development of a 3D velocity model for the larger FORGE area and on the seismic resolution in the stimulated fracture volume at the bottom of well 16A-32. The velocity model was developed using RMS velocities of the seismic reflection survey and seismic velocity logs from borehole measurements as an input model. To improve the accuracy of the model in the shallow subsurface, travel times phase arrivals of the direct propagating P-waves were determined from the seismic reflection data, using PhaseNet, a deep-neural-network-based seismic arrival time picking method. The travel times were subsequently inverted using the input velocity model. The seismic resolution study used borehole and surface seismic sensors as well as the seismicity observed during the April 2022 stimulation experiment to estimate the seismic resolution in the activated fracture reservoir. The data contain a 3D P- and S-wave velocity model for the larger FORGE area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗
DOE OSTI2020

Seismology in the Cloud: A New Streaming Workflow

Data-intensive research in seismology is experiencing a recent boom, driven in part by large volumes of available data and advances in the growing field of data science. However, there are significant barriers to processing large data volumes, such as long retrieval times from data repositories, complex data management, and limited computational resources. New tools and platforms have reduced the barriers to entry for scientific cluster computing, including the maturation of the commercial cloud as an accessible instrument for research. Here, we build a customized research cluster in the cloud to test a new workflow for large-scale seismic analysis, in which data are processed as a stream (retrieved on-the-fly and acted upon without storing), with data from the Incorporated Research Institutions for Seismology Data Management Center. We use this workflow to deploy a spectral peak detection algorithm over 5.6 TB of compressed continuous seismic data from 2074 stations of the USArray Transportable Array EarthScope network. Using a 50-node cluster in the cloud, we completed the noise survey in 80 hr, with an average data throughput of 1.7 GB per minute. By varying cluster sizes, we find the scaling of our analysis to be sublinear, due to a combination of algorithmic limitations and data center response times. The cloud-based streaming workflow represents an order-of-magnitude increase in acquisition and processing speed compared to a traditional download-store-process workflow, and offers the additional benefits of employing a flexible, accessible, and widely used computing architecture. It is limited, however, due to its reliance on Internet transfer speeds and data center service capacity, and may not work well for repeated analyses or those for which even higher data throughputs are needed. These research applications will require a new class of cloud-native approaches in which both data and analysis are in the cloud.

58 GEOSCIENCES↗