Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic methods”

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 361 records · Page 20

Mathematical Morphological Filtering with a Self-Adaptive Reconstruction Technique and Application to Local Seismic Data

Recorded seismic data are generally contaminated by noise from different sources, which masks the signals of interest. In the seismology community, frequency filtering (FF) is the standard method for noise suppression. However, when the signal of interest and noise share the same frequency band, the latter cannot be filtered out without infringing on the former. We implemented a noise suppression approach based on the mathematical morphology theorem. The method involves compound operations of dilation and erosion using structuring elements of varying lengths and decomposes an input noisy waveform into several time functions with differing characteristics. Further, the filtered waveform is constructed from the time functions using a self-adaptive reconstruction technique. Application to a data set of >4700 local waveforms suggests that the implemented mathematical morphological filtering (MMF) approach is efficient for data with low signal-to-noise ratio (SNR) and significantly outperforms FF in that SNR range. For most of the dataset, FF, machine learning (ML) denoising, and continuous wavelet transform (CWT) thresholding result in higher SNR values compared with the MMF method. However, for ~42% of the waveforms, MMF outperforms FF, and the SNR gain achieved with MMF is as large as ~23 dB. Compared to ML denoising and CWT thresholding, this proportion drops to only ~10%–14%. Our results suggests that in an operational setting, MMF cannot replace the other noise suppression methods; however, signal detection can be improved if MMF is used to supplement them in some scenarios. MMF could help detect signals in problematic low-SNR data, which are currently being missed particularly when using FF alone.

58 GEOSCIENCES↗

A Multi-Instance learning Framework for Seismic Detectors

In this report, we construct and test a framework for fusing the predictions of a ensemble of seismic wave detectors. The framework is drawn from multi-instance learning and is meant to improve the predictive skill of the ensemble beyond that of the individual detectors. We show how the framework allows the use of multiple features derived from the seismogram to detect seismic wave arrivals, as well as how it allows only the most informative features to be retained in the ensemble. The computational cost of the "ensembling" method is linear in the size of the ensemble, allowing a scalable method for monitoring multiple features/transformations of a seismogram. The framework is tested on teleseismic and regional p-wave arrivals at the IMS (International Monitoring System) station in Warramunga, NT, Australia and the PNSU station in University of Utah's monitoring network.

58 GEOSCIENCES↗

Verification of Atmospheric Signals Associated with Major Seismicity by Space and Terrestrial Observations

Observations from the last twenty years suggest the existence of electromagnetic (EM) phenomena during or preceding some earthquakes [Haykawa et a!, 2004; Pulinets at al, 1999,2004, 2006, Ouzounov et all 2007 and Liu et all 20041. Both our previous studies [Pulinets at al, 2005, 2006, Ouzounov et al, 2006, 20071 and the latest review by the Earthquake Remote Precursor Sensing panel [ERPS; 2003- 20051; have shown that there were precursory atmospheric TIR signals observed on the ground and in space associated with several recent earthquakes. [Tramutoli at al, 2005, 2006, Cervone et al, 2006, Ouzounov et all 2004,2006JT.o study these signals, we applied both multi parameter statistical analysis and data mining methods that require systematic measurements from an Integrated Sensor Web of observations of several physical and environmental parameters. These include long wave earth infra-red radiation, ionospheric electrical and magnetic parameters, temperature and humidity of the boundary layer, seismicity and may be associated with major earthquakes. Our goal is to verify the earthquake atmospheric correlation in two cases: (i) backward analysis - 2000-2008 hindcast monitoring of multi atmospheric parameters over the Kamchatka region, Russia ; and (ii) forward real-time alert analysis over different seismo-tectonic regions for California, Turkey, Taiwan and Japan. Our latest results, from several post-earthquake independent analyses of more then 100 major earthquakes, show that joint satellite and some ground measurements, using an integrated web, could provide a capability for observing pre-earthquake atmospheric signals by combining the information from multiple sensors into a common framework. Using our methodology, we evaluated and compared the observed signals preceding the latest M7.9 Sichuan earthquake (0511212008), M8.0 earthquake in Peru (0811512007), M7.6 Kashmir earthquake (1010812005) and M9.0 Sumatra earthquake (1212812004). We found evidence of the systematic appearance of both atmospheric and ionospheric anomalies preceding most of the major events during the period of our analysis 2001 -2008.

Taylor, Patrick↗

SeismoGen: Seismic Waveform Synthesis Using GAN With Application to Seismic Data Augmentation

Abstract Detecting earthquake arrivals within seismic time series can be a challenging task. Visual, human detection has long been considered the gold standard but requires intensive manual labor that scales poorly to large data sets. In recent years, automatic detection methods based on machine learning have been developed to improve the accuracy and efficiency. However, the accuracy of those methods relies on access to a sufficient amount of high‐quality labeled training data, often tens of thousands of records or more. We aim to resolve this dilemma by answering two questions: (1) provided with a limited amount of reliable labeled data, can we use them to generate additional, realistic synthetic waveform data? and (2) can we use those synthetic data to further enrich the training set through data augmentation, thereby enhancing detection algorithms? To address these questions, we use a generative adversarial network (GAN), a type of machine learning model which has shown supreme capability in generating high‐quality synthetic samples in multiple domains. Once trained, our GAN model is capable of producing realistic seismic waveforms of multiple labels (noise and event classes). Applied to real Earth seismic data sets in Oklahoma, we show that data augmentation from our GAN‐generated synthetic waveforms can be used to improve earthquake detection algorithms in instances when only small amounts of labeled training data are available.

Wang, Tiantong↗

Back-Azimuth Estimation of Air-to-Ground Coupled Infrasound from Transverse Coherence Minimization

We present the transverse coherence minimization method (TCM)—an approach to estimate the back–azimuth of infrasound signals that are recorded on an infrasound microphone and a colocated three–component seismometer. Accurate back–azimuth information is important for a variety of monitoring efforts, but it is currently only available for infrasound arrays and for seismoacoustic sensor pairs separated by 10 s of meters. Our TCM method allows for the analysis of colocated sensor pairs, sensors located within a few meters of each other, which may extend the capabilities of existing seismoacoustic networks and supplement operating infrasound arrays. This approach minimizes the coherence of the transverse component of seismic displacement with the infrasound wave to estimate the infrasound back–azimuth. After developing an analytical model, we investigate seismoacoustic signals from the August 2012 Humming Roadrunner experiment and the 26 May 2021 eruption of Great Sitkin Volcano, Alaska, U.S.A., at the ranges of 6.5–185 km from the source. We discuss back–azimuth estimates and potential sources of deviation (1°–15°), such as local terrain effects or deviation from common analytical models. This practical method complements existing seismoacoustic tools and may be suitable for routine application to signals of interest.

58 GEOSCIENCES↗

NASA Tech Briefs, November 2012

The topics include: Visual System for Browsing, Analysis, and Retrieval of Data (ViSBARD); Time-Domain Terahertz Computed Axial Tomography NDE System; Adaptive Sampling of Time Series During Remote Exploration; A Tracking Sun Photometer Without Moving Parts; Surface Temperature Data Analysis; Modular, Autonomous Command and Data Handling Software with Built-In Simulation and Test; In-Situ Wire Damage Detection System; Amplifier Module for 260-GHz Band Using Quartz Waveguide Transitions; Wideband Agile Digital Microwave Radiometer; Buckyball Nucleation of HiPco Tubes; FACT, Mega-ROSA, SOLAROSA; An Integrated, Layered-Spinel Composite Cathode for Energy Storage Applications; Engineered Multifunctional Surfaces for Fluid Handling; Polyolefin-Based Aerogels; Adjusting Permittivity by Blending Varying Ratios of SWNTs; Gravity-Assist Mechanical Simulator for Outreach; Concept for Hydrogen-Impregnated Nanofiber/Photovoltaic Cargo Stowage System; DROP: Durable Reconnaissance and Observation Platform; Developing Physiologic Models for Emergency Medical Procedures Under Microgravity; Spectroscopic Chemical Analysis Methods and Apparatus; Low Average Sidelobe Slot Array Antennas for Radiometer Applications; Motion-Corrected 3D Sonic Anemometer for Tethersondes and Other Moving Platforms; Water Treatment Systems for Long Spaceflights; Microchip Non-Aqueous Capillary Electrophoresis (MicronNACE) Method to Analyze Long-Chain Primary Amines; Low-Cost Phased Array Antenna for Sounding Rockets, Missiles, and Expendable Launch Vehicles; Mars Science Laboratory Engineering Cameras; Seismic Imager Space Telescope; Estimating Sea Surface Salinity and Wind Using Combined Passive and Active L-Band Microwave Observations; A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing; Scalable SCPPM Decoder; QuakeSim 2.0; HURON (HUman and Robotic Optimization Network) Multi-Agent Temporal Activity Planner/Scheduler; MPST Software: MoonKommand

Source record↗

Small-Array Location Capabilities Using the Seismometer to Investigate Ice and Ocean Structure (SIIOS): Implications for an Ocean World Lander

Ocean worlds have thick icy shells covering subsurface oceans. Due to the potential habitability of the subsurface ocean, Europa has become a target for a potential lander mission. Seismology is the preeminent method for constraining the thickness of an icy shell. The Seismometer to Investigate Ice and Ocean Structure (SIIOS) uses flight-candidate instrumentation to develop approaches for seismic studies of icy bodies. The SIIOS team deployed small aperture seismic arrays on Gulkana Glacier in 2017 and in Northwest Greenland in 2018.

Marusiak, A. G.↗

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗

Seismic monitoring and high-frequency noise using arrays in the Arabian Peninsula

Seismic arrays provide useful tools for regional seismic monitoring. Two small-aperture, regionally-oriented arrays, QWAR and HQAR, were deployed in Saudi Arabia and Oman in 2012 and 2016. We utilize a time-domain beampacking method, similar to frequency-wavenumber analysis, to examine the performance of the arrays in terms of slowness and azimuthal bias and event detection capabilities. Additionally, we investigate persistent ambient noise sources recorded by the arrays. We find that the arrays provide slowness vectors with biases comparable with similar-sized arrays of the International Monitoring System in other locations around the world. At QWAR, regional events of magnitude 3.0 and above are detected a majority of the time, and as magnitudes increase to 4.0 and above, the detection rate is greater than 82%. Finally, strong noise generation, primarily with slownesses characteristic of Lg waves, is found in the directions of the northern Arabian/Persian Gulf and southern Red Sea and may be a factor in event detection capabilities.

58 GEOSCIENCES↗

Proximal Observations of Epicentral Infrasound Generated by Shallow Low-Magnitude Earthquakes in the Permian Basin, West Texas

Infrasound generated by earthquakes and explosions is generally detected at receivers at epicentral distances of hundreds to thousands of kilometers. However, proximal (<50 km) observations are especially important for low-magnitude earthquakes and low-yield explosions that may not generate signals capable of being detected at great ranges. Here, in this study, we present on the signals detected on an infrasound array 3 km away from two M L 2.9 earthquakes in the Permian Basin of west Texas. Local infrasound (LIS), generated at receivers during the passage of seismic waves, was detected following each earthquake. Epicentral infrasound (EIS), created at or near the epicenter and propagating away as a sound wave, was also detected. Array processing methods show that the EIS signals arrive from the same direction as the earthquake epicenters and at acoustic speeds. To our knowledge, these are the first observations of laterally propagating EIS at proximal ranges following an earthquake of any magnitude.

58 GEOSCIENCES↗

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗

CO 2 rock physics modeling for reliable monitoring of geologic carbon storage

Monitoring, verification, and accounting (MVA) are crucial to ensure safe and long-term geologic carbon storage. Seismic monitoring is a key MVA technique that utilizes seismic data to infer elastic properties of CO 2 -saturated rocks. Reliable accounting of CO 2 in subsurface storage reservoirs and potential leakage zones requires an accurate rock physics model. However, the widely used CO 2 rock physics model based on the conventional Biot-Gassmann equation can substantially underestimate the influence of CO 2 saturation on seismic waves, leading to inaccurate accounting. We develop an accurate CO 2 rock physics model by accounting for both effects of the stress dependence of seismic velocities in porous rocks and CO 2 weakening on the rock framework. We validate our CO 2 rock physics model using the Kimberlina-1.2 model (a previously proposed geologic carbon storage site in California) and create time-lapse elastic property models with our new rock physics method. We compare the results with those obtained using the conventional Biot-Gassmann equation. Our innovative approach produces larger changes in elastic properties than the Biot-Gassmann results. Using our CO 2 rock physics model can replicate shear-wave speed reductions observed in the laboratory. Our rock physics model enhances the accuracy of time-lapse elastic-wave modeling and enables reliable CO 2 accounting using seismic monitoring.

58 GEOSCIENCES↗

Computational Tools and Workflows for Quantitative Risk Assessment and Decision Support for Geologic Carbon Storage Sites: Progress and Insights from the U.S. DOE’s National Risk Assessment Partnership

The 2005 Intergovernmental Panel on Climate Change (IPCC) Special Report on CCS raised the profile of CO2 capture and storage (CCS) as an important technology for reducing greenhouse gas (GHG) emissions. CCS is now recognized as a key component of most climate change mitigation scenarios. Since publication of that report the international research, development, and deployment (RD&D) community has advanced key technical aspects, clarified regulatory requirements, explored value chain and infrastructure solutions, and developed incentive paradigms to enable and promote large-scale deployment of CCS. These efforts have included research to better characterize geologic storage resources, to improve injection performance and storage efficiency, to assess and manage subsurface environmental risks, and to advance monitoring technologies to assure system conformance. These efforts have helped to build confidence in the viability of geologic carbon storage (GCS), but stakeholder concerns about long-term risks and liability associated with GCS remain a hurdle to broad acceptance and large-scale deployment of CCS. Since 2010, the U.S. DOE’s National Risk Assessment Partnership (NRAP) – a research collaboration between five contributing national laboratories – has worked to establish and demonstrate methods and tools to quantify and manage the subsurface environmental risks associated with GCS, amidst uncertainty. This work supports the Office of Fossil Energy and Carbon Management Carbon Transport and Storage Program’s goal of advancing safe and secure commercial-scale GCS deployment. To address the technical challenge of simulating the physical response of the GCS site to large-scale CO2 injection, NRAP has adopted an approach that relies on coupling computationally efficient reduced-order and/or data-driven proxy models of important system components (i.e., storage reservoir, sealing caprock, leakage pathways, intermediate formations, overlying groundwater aquifers, and the atmosphere) in integrated assessment framework. That integrated model of the physical system is complemented with fit-for purpose functionality to support site characterization and risk-related decisions. The recently released NRAP Phase II toolset includes the Open-Source Integrated Assessment Model (NRAP-Open-IAM) for evaluation of trends in leakage risk and potential impact, tools to support monitoring design optimization (Designs for Risk Evaluation and Management – DREAM v3.0 and Passive Seismic Monitoring Tool - PSMT), and tools for state of stress evaluation (State-of-Stress Analysis Tool - SOSAT) and forecasting induced seismicity risk. The NRAP team has also released a pair of reports describing conceptual workflows to incorporate physics-based, quantitative risk assessment into many of the design, planning, operation, and closure decisions for GCS projects. An online catalogue highlights published studies where these tools and methods are demonstrated. In this presentation, the utility of these products to assess risks and address key stakeholder questions will be highlighted through examples, and related insights about the safety and security of geologic carbon storage in qualified storage sites will be discussed. The prospect of rapid, large-scale deployment of GCS technology to aggressively reduce anthropogenic CO2 emissions requires careful consideration of interference between multiple commercial-scale storage projects within a basin. Going forward, NRAP is expanding and adapting site-scale risk quantification tools and methods to enable assessment of risks and inform management decisions for basin-scale deployment. Increasingly, this work will leverage next-generation approaches for surrogate modelling, fast prediction, and advanced visualization enabled by machine learning and artificial intelligence to promote virtual learning, scenario evaluation, and augment risk-based decision making.

quantitative risk assessment, geologic carbon stor↗

On the Feasibility of Geophysical Methods for CO 2 Monitoring in the North Dakota CarbonSAFE Project

Conference presentation at American Institute of Chemical Engineers (AIChE) Annual Meeting, Phoenix, AZ, November 13–18, 2022. An analysis of rock physics crossplots and diagnostic models suggests that the baseline seismic data may be sensitive to changes in lithology and porosity of the Broom Creek Formation and that these changes may be detectable by prestack seismic inversion.

20 FOSSIL-FUELED POWER PLANTS↗

Wellbore integrity assessment with casing-based advanced sensing

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage 1 and waste disposal. After installation, well casing and cement are subject to mechanical stress due to near-well pressure changes and fluid induced corrosion. This is exacerbated for geothermal wells where produced fluid is at high temperature and corrosive. The current state-of-the-art technologies for wellbore integrity assessments are an array of cased hole logging tools. Wireline deployed acoustic, electromagnetic and mechanical tools are all available to inspect steel casing corrosion and casing-cement bond and these tools can provide high-resolution assessment of borehole conditions. They are intrusive, however, in terms of borehole preparation and interruption to the normal operation of the wells, and not suitable for high temperature or highly deviated well deployments. In addition, these measurements are performed infrequently due to high cost, and are therefore incapable of providing frequent data to better predict borehole degradation trajectory, which can help provide early warning of potential borehole failures. For this project we are developing a suite of novel, non-invasive, casing based tools for wellbore integrity assessment, combining fast/low cost screening with higher-precision investigation. Our approach is based on monitoring the response of the casing when energized at the wellhead, thereby interrogating the casing without well intervention. Lab, field and numerical approaches are used in our study. During the early stage of the research, we focus on numerical simulations, which have shown the sensitivity of the low frequency electromagnetic (EM) signals to changes in borehole depths and have successfully tested the concept at a field site with different length well casings. Initial seismic modeling efforts have also demonstrated our capability to simulate seismic tube wave and seismic field alterations due to borehole breakage and associated fluid leakage. Further numerical, laboratory and field experiments are underway for additional technology sensitivity analysis, particularly the transient EM/Seismic reflectometry methods, data acquisition optimization, and numerical simulation improvements.

Wilt, Michael↗

Deep learning multiphysics network for imaging CO 2 saturation and estimating uncertainty in geological carbon storage

Multiphysics inversion exploits different types of geophysical data that often complement each other and aims to improve overall imaging resolution and reduce uncertainties in geophysical interpretation. Despite the advantages, traditional multiphysics inversion is challenging because it requires a large amount of computational time and intensive human interactions for preprocessing data and finding trade-off parameters. These issues make it nearly impossible for traditional multiphysics inversion to be applied as a real-time monitoring tool for geological carbon storage. In this paper, we present a deep learning (DL) multiphysics network for imaging CO 2 saturation in real time. The multiphysics network consists of three encoders for analysing seismic, electromagnetic and gravity data and shares one decoder for combining imaging capabilities of the different geophysical data for better predicting CO 2 saturation. The network is trained on pairs of CO 2 label models and multiphysics data so that it can directly image CO 2 saturation. Here we use the bootstrap aggregating method to enhance the imaging accuracy and estimate uncertainties associated with CO 2 saturation images. Using realistic CO 2 label models and multiphysics data derived from the Kimberlina CO 2 storage model, we evaluate the performance of the deep learning multiphysics network and compare its imaging results to those from the deep learning single-physics networks. Our modelling experiments show that the deep learning multiphysics network for seismic, electromagnetic, and gravity data not only improves the imaging accuracy but also reduces uncertainties associated with CO 2 saturation images. Our results also suggest that the deep learning multiphysics network for the non-seismic data (i.e., electromagnetic and gravity) can be used as an effective low-cost monitoring tool in between regular seismic monitoring.

58 GEOSCIENCES↗

FY 2026 Midyear Report: Seismic Monitoring of Underground Vibration Sources Using Distributed Acoustic Sensing and Seismometers

Safeguards-relevant temporal changes in underground facilities can be observed using geophysical monitoring techniques. Seismic waves, in particular, provide valuable insights into subsurface activities and can serve as an important tool for detecting anomalous events that may indicate containment breaches at geological repositories. This midyear report summarizes ongoing efforts to automatically and rapidly detect and locate anomalous vibration signals that could be indicative of potential containment breaches. Previous work during FY25 focused on compiling continuous seismic datasets from two underground sites and developing a database of continuous waveforms and ground-truth event data derived from multiple sensing modalities. Building on this foundation, we are adapting anomaly detection and geolocation algorithms to explore methods for monitoring underground activities using two relatively low-maintenance sensing technologies: a dense surface geophone array deployed at the Pleasant Gap mine in Pennsylvania, and a three-dimensional fiber-optic cable array for distributed acoustic sensing (DAS) installed in the subsurface at the Sanford Underground Research Facility (SURF) in South Dakota. This report summarizes work conducted during the first two quarters of FY26, during which we refined a dynamic power spectral density (PSD)-based detector, applied it independently to each geophone station, and then combined the per‑station detections with density-based spatial clustering of applications with noise (DBSCAN) to cluster events and produce spatial maps over a nine‑day interval. In addition, we outline plans for a field trial at the Waste Isolation Pilot Plant (WIPP) in New Mexico to compare traditional seismic monitoring approaches with DAS techniques and to evaluate the benefits of combined data analysis. Activities during the past two quarters have included the preparation and submission of a Field Test Plan to WIPP for approval, as well as submission to headquarters for review and feedback.

58 GEOSCIENCES↗

IC reports (2021)

We requested HPC support to continue research on the seismic waves generated by impacts. We used the following codes: (1) the Hybrid Optimization Software Suite (HOSS), developed at LANL. HOSS is based on a combined Finite and Discrete Element Method (FDEM). New material models are developed for the sedimentary rocks. HOSS has been recently benchmarked to iSale and FLAG codes. SPECFEM3D is an open-source code developed since the last 90s. It won the Gordon Bell award for best performance in 2003, was finalist again in 2008 for a run at 0.16 petaflops on 149,784 cores on the ‘Jaguar’ Cray system at Oak Ridge National Laboratory. It also won the BULL Joseph Fourier supercomputing award in 2010.; SW4 is a 4th-order finite difference code developed at LLNL which is currently actively developed to handle complex 3D models and to be ported on future exascale platforms. We assessed our need to a total of 3.9M CPU-hrs for year 1 and and 3.1 M for year 2.

79 ASTRONOMY AND ASTROPHYSICS↗