Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-pressure elastic properties of dolomite melt supporting carbonate-induced melting in deep upper mantle

Deeply subducted carbonates likely cause low-degree melting of the upper mantle and thus play an important role in the deep carbon cycle. However, direct seismic detection of carbonate-induced partial melts in the Earth’s interior is hindered by our poor knowledge on the elastic properties of carbonate melts. Here we report the first experimentally determined sound velocity and density data on dolomite melt up to 5.9 GPa and 2046 K by in-situ ultrasonic and sink-float techniques, respectively, as well as first-principles molecular dynamics simulations of dolomite melt up to 16 GPa and 3000 K. Using our new elasticity data, the calculated V P /V S ratio of the deep upper mantle (~180–330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated V P /V S ratio of the upper mantle from global seismic observations, supporting the pervasive presence of a low-degree carbonate-rich partial melt (~0.05%) that is consistent with the volatile-induced or redox-regulated initial melting in the upper mantle as argued by petrologic studies. This carbonate-rich partial melt region implies a global average carbon (C) concentration of 80–140 ppm. by weight in the deep upper mantle source region, consistent with the mantle carbon content determined from geochemical studies.

58 GEOSCIENCES↗
NASA NTRS1984

Solar Seismology from Space. a Conference at Snowmass, Colorado

The quality of the ground based observing environment suffers from several degrading factors: diurnal interruptions and thermal variations, atmospheric seeing and transparency fluctuations and adverse weather interruptions are among the chief difficulties. The limited fraction of the solar surface observable from only one vantage point is also a potential limitation to the quality of the data available without going to space. Primary conference goals were to discuss in depth the scientific return from current observations and analyses of solar oscillations, to discuss the instrumental and site requirements for realizing the full potential of the seismic analysis method, and to help bring new workers into the field by collecting and summarizing the key background theory. At the conclusion of the conference there was a clear consensus that ground based observation would not be able to provide data of the quality required to permit a substantial analysis of the solar convection zone dynamics or to permit a full deduction of the solar interior structure.

Ulrich, R. K.↗
DOE OSTI2021

High-resolution Bayesian spatial autocorrelation (SPAC) quasi-3-D Vs model of Utah FORGE site with a dense geophone array

SUMMARY We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2–5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ∼22° to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity.

Geochemistry & Geophysics↗
DOE OSTI2021

InversionNet3D: Efficient and Scalable Learning for 3-D Full-Waveform Inversion

Seismic full-waveform inversion (FWI) techniques aim to find a high-resolution subsurface geophysical model provided with waveform data. Some recent effort in data-driven FWI has shown some encouraging results in obtaining 2-D velocity maps. However, due to high computational complexity and large memory consumption, the reconstruction of 3-D high-resolution velocity maps via deep networks is still a great challenge. Here, in this article, we present InversionNet3D (InvNet3D), an efficient and scalable encoder–decoder network for 3-D FWI. The proposed method employs group convolution in the encoder to establish an effective hierarchy for learning information from multiple sources while cutting down unnecessary parameters and operations at the same time. The introduction of invertible layers further reduces the memory consumption of intermediate features during training and, thus, enables the development of deeper networks with more layers and higher capacity as required by different application scenarios. Experiments on the 3-D Kimberlina dataset demonstrate that InvNet3D achieves state-of-the-art reconstruction performance with lower computational cost and lower memory footprint compared to the baseline.

58 GEOSCIENCES↗
DOE OSTI2024

3D Multiresolution Velocity Model Fusion with Probability Graphical Models

ABSTRACT The variability in spatial resolution of seismic velocity models obtained via tomographic methodologies is attributed to many factors, including inversion strategies, ray-path coverage, and data integrity. Integration of such models, with distinct resolutions, is crucial during the refinement of community models, thereby enhancing the precision of ground-motion simulations. Toward this goal, we introduce the probability graphical model (PGM), combining velocity models with heterogeneous resolutions and nonuniform data point distributions. The PGM integrates data relations across varying resolution subdomains, enhancing detail within low-resolution (LR) domains by utilizing information and prior knowledge from high-resolution (HR) subdomains through a maximum posterior problem. Assessment of efficacy, utilizing both 2D and 3D velocity models—consisting of synthetic checkerboard models and a fault-zone model from Ridgecrest, California—demonstrates noteworthy improvements in accuracy, compared to state-of-the-art fusion techniques. Specifically, we find reductions of 30% and 44% in computed travel-time residuals for 2D and 3D models, respectively, as compared to conventional smoothing techniques. Unlike conventional methods, the PGM’s adaptive weight selection facilitates preserving and learning details from complex, nonuniform HR models and applies the enhancements to the LR background domain.

Geochemistry & Geophysics↗
DOE OSTI2022

BayesMT: A Probabilistic Bayesian Framework for the Seismic Moment Tensor

Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.

58 GEOSCIENCES↗
NASA NTRS1977

Geology of the Imbrium Basin Apennine Mountains and relation to the Apollo 15 landing site

Detailed mapping of the Imbrium Basin Apennine Mountains has emphasized several important factors in the geologic evolution of this region. Pre-Imbrian topography was composed of two types of landforms: (1) the nearby Serenitatis, South Imbrium, Aestuum, and Vaporum basins which produced widely spaced regional lows and a sequence of basin deposits; (2) pre-Imbrian craters that reworked existing basin deposits and produced local topographic lows. The sequence of events during the Imbrium impact can be summarized as follows: intensive structural uplift along the Apennine Front occurred prior to the emplacement of ballistic ejecta. Ejecta became more discontinuous toward the backslope and produced radial texture by sedimentary erosional and depositional processes. Post-ejecta-emplacement downfaulting exposed pre-Imbrian deposits along the front and caused severe seismic disturbances that reworked many of the previous units and accentuated pre-basin structures. Tectonic mapping and orbital geophysical data support the idea that little structural modification has occurred in the Apennines since their formation 3.9 b.y. ago.

Spudis, P.↗
NASA NTRS1993

GPS measurements of crustal deformation associated with the 22 April 1991, Valle de la Estrella , Costa Rica earthquake

Results are presented for analysis of Global Positioning System (GPS) measurements made at sites in Costa Rica in February and July 1991. Significant horizontal and vertical displacements relative to February positions were observed. Differences were found in GPS derived vertical and horizontal displacements as compared to other types of geodetic measurements of uplift in the coastal regions. A slip dislocation model which fits the GPS measured displacement was computed. Differences between these data sets and their associated models were investigated in light of unmodeled slip heterogeneity on the fault and post-seismic displacements.

Lundgren, Paul R.↗
NASA NTRS2016

MarCO: Interplanetary Mission Development on a CubeSat Scale

Shortly after JPL’s Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) mission launches, separates, and commences its cruise phase, two CubeSats will deploy from the launch vehicle’s upper stage and begin independent flight to Mars (Fig. 1). During InSight’s entry, descent, and landing (EDL) sequence, these twin Mars Cube One (MarCO) spacecraft will fly 3,500 km above the Martian surface, recording and relaying InSight UHF radio data to the Deep Space Network (DSN) on Earth1. MarCO is a twin CubeSat mission developed by the NASA Jet Propulsion Laboratory (JPL) to accompany the InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) Mars mission lander. MarCO's primary mission objective is to launch with InSight and independently fly to Mars to serve as a communications relay during InSight's entry, descent, and landing (EDL) phase. MarCO represents a new type of deep space mission: CubeSats at Mars. Building on the development of JPL's first interplanetary CubeSat project, the Interplanetary Nano-Spacecraft Pathfinder in Relevant Environment (INSPIRE), MarCO further refined the approach to hardware, software, and ground architecture development to solve the challenges of quickly building low-budget spacecraft to fly to Mars. The greatest constraint, beyond others typical of CubeSat missions, was time. The duration between MarCO's conception to completion of spacecraft assembly was less than two years - an unprecedented schedule for any planetary mission to date. Through necessity, MarCO has built on previous experience, procedures, systems, and development methodologies, defining a new niche for supporting larger primary missions. The MarCO spacecraft are poised to write a new chapter in deep space exploration. Originally slated to launch and reach Mars in 2016, the InSight mission schedule subsequently slipped to 2018. During the original landing of InSight, Earth would not be in view, and no orbiter around Mars would have been in position to both receive UHF EDL data and simultaneously relay it back to Earth. It was from this obstacle that MarCO was conceived. Regardless of any changes to InSight’s 2018 EDL configuration geometry, MarCO is still expected to fly and serve in the same capacity as originally designed: the first CubeSat mission to Mars. CubeSats have historically been firmly in the domain of universities and small companies. As first conceived, they served as a platform upon which to teach all aspects of the space mission lifecycle. JPL took on this mission type with Interplanetary Nano-Spacecraft Pathfinder in Relevant Environment2 (INSPIRE), moving the concept into a new domain: deep space. Building from the INSPIRE platform and lessons learned, MarCO addressed new challenges in the domain of planetary missions: independent interplanetary flight and navigation, integration with a large-scale mission, long-distance and long-delay communication, short development time, and a small development team. Of these, the greatest constraint was schedule: only 18 months passed from conception of mission concept until delivery of fully assembled and tested flight hardware. Careful selection of mission team, along with extensive use of off-the-shelf equipment, and streamlining automated processes, was essential. This achievement represents the next step in the evolution of CubeSats beyond low-Earth orbit.

Werne, Thomas↗
NASA NTRS1977

Faulting parameters derived from computer simulation of earthquakes

Seismic source parameters, average displacement, rupture length, and strain energy release are investigated by computer simulation using a coupled massive block model of the sliding along an active fault. Average displacements and energy release vary considerably with the degree of heterogeneity in the friction and elastic parameters used in the model. Strain energy release is determined primarily by the product of dynamic friction, rupture length, and average displacement. Interrelationships among the faulting parameters are consistent with theoretical arguments and experimental data. The variation in the frequency of occurrence of simulation events with strain energy release is different from the variation in the frequency of naturally occurring events with seismic energy.

Cohen, S. C.↗
NASA NTRS1973

Lunar seismic profiling experiment

The Apollo 17 lunar seismic profiling experiment was conducted to record the vibrations of the lunar surface as induced by explosive charges, the thrust of the lunar module ascent engine, and the crash of the lunar module ascent stage. Analysis of the data obtained made it possible to determine the internal characteristics of the lunar crust to a depth of several kilometers. The test equipment used in the experiment is described. Maps showing the location of the geophones and the deployed explosive packages are provided. Samples of the seismic signals recorded by the lunar seismic profiling experiment geophones are included.

Robert L. Kovach↗
DOE OSTI2021

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗
DOE OSTI2022

The EGS Collab Project – Stimulations at Two Depths

The EGS Collab project, supported by the US Department of Energy, is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to address challenges in implementing enhanced geothermal systems (EGS). Data and observations from the field tests are compared to simulations to understand processes and build confidence in numerical modeling of the processes. We have completed Experiment 1 (of 3), which examined hydraulic fracturing in a well-characterized underground fractured phyllite test bed at a depth of approximately 1.5 km at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. Testbed characterization included fracture mapping, borehole acoustic and optical televiewers, full waveform sonic, conductivity, resistivity, temperature, campaign p- and s-wave investigations and electrical resistance tomography. Borehole geophysical techniques including passive seismic, continuous active source seismic monitoring, electrical resistance tomography, fiber-based distributed strain, distributed temperature, and distributed acoustic monitoring, were used to carefully monitor stimulation events and flow tests. More than a dozen stimulations and nearly one year of flow tests were performed. Quality data and detailed observations were collected and analyzed during stimulation and water flow tests using ambient temperature and chilled water. We achieved adaptive control of the tests using real-time monitoring and rapid dissemination of data and near-real-time simulation. More detailed numerical simulation was performed to answer key experimental design questions, forecast fracture propagation trajectories and extents, and analyze and evaluate results. Data are freely available from the Geothermal Data Repository. Experiment 2 examines the potential for hydraulic shearing in amphibolite at a depth of about 1.25 km at SURF. This site has a different set of stress and fracture conditions than Experiment 1. The Experiment 2 testbed consists of nine subhorizontal boreholes configured in two fans of two boreholes which surround the testbed and contain grouted-in electrical resistance tomography, seismic sensors, active seismic sources and distributed fiber sensors. A “five-spot” set of test wells that extends from a custom mined alcove includes an injection well and four production/monitoring wells. The testbed was characterized geophysically and hydrologically, and three stimulations have been performed using the Step-Rate Injection Method for Fracture In-Situ Properties (SIMFIP) tool to measure strains, and a new strain quantifying tool (downhole robotic strain analysis tool -DORSA) was deployed in a monitoring hole during stimulation. Real-time data were broadcast during stimulations to allow real-time response to arising issues.

EGS Collab, Enhanced Geothermal Systems, EGS, fiel↗
NASA NTRS1988

Relationships of earthquakes (and earthquake-associated mass movements) and polar motion as determined by Kalman filtered, Very-Long-Baseline-Interferometry

A Kalman filter was designed to yield optimal estimates of geophysical parameters from Very Long Baseline Interferometry (VLBI) group delay data. The geophysical parameters are the polar motion components, adjustments to nutation in obliquity and longitude, and a change in the length of day parameter. The VLBI clock (and clock rate) parameters and atmospheric zenith delay parameters are estimated simultaneously. Filter background is explained. The IRIS (International Radio Interferometric Surveying) VLBI data are Kalman filtered. The resulting polar motion estimates are examined. There are polar motion signatures at the times of three large earthquakes occurring in 1984 to 1986: Mexico, 19 September, 1985 (Magnitude M sub s = 8.1); Chile, 3 March, 1985 (M sub s = 7.8); and Taiwan, 14 November, 1986 (M sub s = 7.8). Breaks in polar motion occurring about 20 days after the earthquakes appear to correlate well with the onset of increased regional seismic activity and a return to more normal seismicity (respectively). While the contribution of these three earthquakes to polar motion excitations is small, the cumulative excitation due to earthquakes, or seismic phenomena over a Chandler wobble damping period may be significant. Mechanisms for polar motion excitation due to solid earth phenomena are examined. Excitation functions are computed, but the data spans are too short to draw conclusions based on these data.

Preisig, Joseph Richard Mark↗
DOE OSTI2024

Multi-phenomenology Yield Characterization

This report serves as the first delivery of a four-year applied science effort to transform and advance the error bounds for the yield estimate of an explosion. Each year’s delivery will be in this form, culminating in the submission of this work for peer review to a scientific journal. Importantly, the yearly progress reports can then also be viewed as expanding drafts working towards a formal journal article submission. For the first tranche of funding, we collaborated with Air Force Technical Applications Center (AFTAC) scientists to identify unclassified real-world data that demonstrate and validate our advanced error propagation methods. Collaboration includes visits to AFTAC and telecons. For this development, we illustrate the fusion of seismic, acoustic, optical, and surface effect signatures from an explosion. The mathematics and code being adapted to this specific application (Williams et al., 2021) involves physics models of multiple sensor signatures. We have also identified related physics models and have integrated them into code. Current methods of underground explosion yield estimation for the Threshold Test Ban Treaty (TTBT) have served the US treaty monitoring mission well for decades. A research objective of the Defense Nuclear Nonproliferation Research and Development (DNN R&D) office of the National Nuclear Security Administration (NNSA) has always been to provide new technical capabilities for monitoring lower thresholds. The general model and error propagation code to be developed in this project is based on significant advances in error modeling and propagation needed to analyze data at lower detection thresholds. The second tranche of funding for this project began on May 1, 2022, and planned work for the second tranche includes: i) completing the integration of physical model code into the general error model framework; this code accommodates a wide range of linear/nonlinear source models, fixed/ random effects, and frequentist/Bayesian analyses (the purpose of which is not to dictate to users how to analyze data, but instead to allow users the maximum flexibility in their work); ii) illustrative application of code to identified data, and; iii) initial planning with AFTAC researchers on delivery of code to the Common Development Environment at AFTAC, and continued writing of the planned final journal article submission (year two of this progress report), with particular emphasis on descriptions of data identified for this effort.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Utah FORGE: Composite 3D Seismic Velocity Model

This is a composite 3D seismic velocity that was constructed from compiled information from several local studies regarding seismic velocities and structural information. This seismic velocity model is provided in NonLinLoc format (slow_len), which is readily usable in NonLinLoc software. Other model formats and versions of the model can be produced using the Python script provided with this data set. Details on how the model was created and prior velocity and structural information was used is provided in the accompanying documentation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

On Magnitudes and Inferred Yields of the 39 Underground Nuclear Test Explosions at the Novaya Zemlya Test Site

ABSTRACT We review the previously compiled seismic magnitudes for underground nuclear explosions at the Union of the Soviet Socialist Republics Novaya Zemlya Test Site, and analyze these data in the context of published details about individual tests and official Russian reports of annual cumulative explosive energy release at the test site. We derive relationships between teleseismic body wave magnitude mb and explosive yield Y by predicting the yield of each test using mb and minimizing the difference between predicted and reported annual cumulative yield. A linear relationship between mb and logY with a slope of 0.82 provides a good fit to the annual total yield, with a root mean square (rms) misfit of 0.095 logY units. A spectral magnitude relationship based on the explosion source model of Mueller and Murphy (1971) with a single-free-parameter provides a slightly better fit. For the 10 tests for which individual yields can be directly deduced from the official reports, the rms misfit of reported yields is less than 0.1 logY units or about 25%. More complex relationships that differentiate between the northern and southern test areas, or take account of variations in the scaled depth of burial, do not provide substantially better fits to the reported yields.

Geochemistry & Geophysics↗
DOE OSTI2025

Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS): Task 4.0 Deliverable – Geologic Analysis Report

Project OASIS (Optimizing Alabama’s CO 2 Storage in Shelby County) is a geologic and reservoir characterization study designed to evaluate deep saline formations for potential long-term carbon dioxide (CO 2 ) storage in central Alabama near the National Carbon Capture Center (NCCC) and Alabama Power’s Plant Gaston. The project centers on the Cambro-Ordovician Knox Group and underlying strata such as the Conasauga and Rome Formations, which were investigated through the drilling of two stratigraphic test wells to obtain electronic well logs, core, and sidewall core plugs. These data provide direct measurements of porosity, permeability, and lithologic variability critical for reservoir characterization. Complementing the well program, a limited 2D seismic survey was acquired to help select the site for Westover #2, and a more regional Seismic Exchange (SEI) seismic survey was licensed and interpreted to define structural and stratigraphic frameworks in a new Static Earth Model (SEM), map reservoir continuity, and identify potential sealing intervals. Integrated with geologic and reservoir modeling, these datasets form the basis for evaluating storage capacity, injectivity, and containment.

20 FOSSIL-FUELED POWER PLANTS↗