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 145 records · Page 8

Integrating Atmospheric Specifications into Seismoacoustic Event Localization

This report investigates the integration of infrasound and seismic data to improve event localization accuracy, specifically focusing on a surface explosion at the Utah Training and Testing Range (UTTR). Utilizing the Seismoacoustic Bayesian Event Locator (SABEL) framework, we incorporated atmospheric specifications derived from Ground to Space (G2S) profiles to enhance celerity-range priors. Our analysis revealed that while the combination of infrasound and seismic observations significantly reduced localization uncertainty, challenges remained, particularly with returns at distances less than 200 km from the source and the influence of specific observations on location estimates. The results indicate that broader celerity distributions, such as those from Blom et al. (2020), facilitate better alignment with ground truth locations compared to narrower models. Overall, this work demonstrates the promise of seismoacoustic approaches in refining event localization and highlights the need for further exploration of celerity-range models to ensure reliable outcomes.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

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↗

Constraints on lunar structure

A brief review is given of the constraints placed on lunar structure and composition by seismic data and density models. Bounds on the crustal velocity structure in Mare Cognitum are derived using travel-time data from artificial impacts, and a velocity model is determined on the basis of synthetic seismograms. It is shown that the velocities of P- and S-waves in the mantle can be fixed by a least-squares analysis of arrival times from meteor impacts and moonquakes, and that lunar density can be determined from the seismic structure, mean density, and moment of inertia. Olivine-pyroxene mixtures and certain olivine-rich compositions are found to be consistent with the seismic-velocity and density limits. Maximum radii are calculated for pure Fe and pure FeS cores, and it is concluded that the possibility of an ancient lunar magnetic dynamo may have to be reevaluated in the light of these figures.

Dainty, A. M.↗

Inverting travel times with a triplication

A method based on the use of constrained spline fits is used to overcome the difficulties arising when body-wave data in the form of T-delta are reduced to the tau-p form in the presence of cusps. In comparison with unconstrained spline fits, the method proposed here tends to produce much smoother models which lie approximately in the middle of the bounds produced by the extremal method. The method is noniterative and, therefore, computationally efficient. The method is applied to the lunar seismic data, where at least one triplication is presumed to occur in the P-wave travel-time curve. It is shown, however, that because of an insufficient number of data points for events close to the antipode of the center of the lunar network, the present analysis is not accurate enough to resolve the problem of a possible lunar core.

Jarosch, H. S.↗

Nowcasting Earthquakes With Stochastic Simulations: Information Entropy of Earthquake Catalogs

Earthquake nowcasting has been proposed as a means of tracking the change in large earthquake potential in a seismically active area. The method was developed using observable seismic data, in which probabilities of future large earthquakes can be computed using Receiver Operating Characteristic methods. Furthermore, analysis of the Shannon information content of the earthquake catalogs has been used to show that there is information contained in the catalogs, and that it can vary in time. So an important question remains, where does the information originate? In this paper, we examine this question using stochastic simulations of earthquake catalogs. Our catalog simulations are computed using an Earthquake Rescaled Aftershock Seismicity (“ERAS”) stochastic model. This model is similar in many ways to other stochastic seismicity simulations, but has the advantage that the model has only 2 free parameters to be set, one for the aftershock (Omori-Utsu) time decay, and one for the aftershock spatial migration away from the epicenter. Generating a simulation catalog and fitting the two parameters to the observed catalog such as California takes only a few minutes of wall clock time. While clustering can arise from random, Poisson statistics, we show that significant information in the simulation catalogs arises from the “non-Poisson” power-law aftershock clustering, implying that the practice of de-clustering observed catalogs may remove information that would otherwise be useful in forecasting and nowcasting. We also show that the nowcasting method provides similar results with the ERAS model as it does with observed seismicity.

58 GEOSCIENCES↗

Internal friction quality-factor Q under confining pressure

It has been found in previous studies that small amounts of adsorbed volatiles can have a profound effect on the internal friction quality-factor Q of rocks and other porous media. Pandit and Tozer (1970) have suggested that the laboratory-measured Q of volatile-free rocks should be similar to the in situ seismic Q values of near-surface lunar rocks which according to Latham et al. (1970) are in the range of 3000-5000. Observations of dramatic increases in Q with outgassing up to values approaching 2000 in the seismic frequency range confirm this supposition. Measurements under confining pressures with the sample encapsulated under hard vacuum are reported to aid in the interpretation of seismic data obtained below the lunar surface. It has been possible to achieve in the experiments Q values just under 2000 at about 1 kbar for a terrestrial analog of lunar basalt. It was found that a well-outgassed sample maintains a high Q whereas one exposed to moisture maintains a low Q as the confining pressure is raised to 2.5 kbar. This result suggests that volatiles can indeed affect Q when cracks are partially closed and the high lunar seismic Q values reported are concomitant with very dry rock down to depths of at least 50 km.

Tittmann, B. R.↗

Lunar Seismometer and Burial System

Beginning in 1969, Apollo successfully deployed a long-lived network of seismometers on the Moon. Seismic studies provide definitive knowledge of internal planetary structure, and analysis of the Apollo seismic data has contributed to the magma ocean hypothesis for initial terrestrial planetary differentiation [Wieczoreket al., 2006]. While the general model is widely accepted, details such as mantle composition, stratification and possible overturn, lateral structure, and thermal inhomogeneities remain unresolved. The Moon experiences moonquakes at varying depths [Nakamura, 1983]. Shallow quakes are relatively large but rare, similar to terrestrial intra-plate earthquakes. Deeper quakes are comparatively smaller but more frequent, occurring periodically according to the tidal cycle. On the Moon, the lack of an atmosphere enables seismic experiments to potentially constrain meteorite impact flux, which informs cratering rates assumed throughout the solar system. The large diurnal temperature variation between day and night also induces thermal moonquakes, which may contribute to regolith production [Duennebier& Sutton, 1974; Weber et al., 2017]. Still, many questions remain regarding the frequency and distribution of natural moonquakes. This translates into an incomplete understanding of the Moon’s hemispherical dichotomies in crustal thickness, mare volcanism, seismicity, and the distribution of heat-producing elements. The Planetary Decadal Survey (National Research Council, 2013) identifies a New Frontiers Lunar Geophysical Network (LGN) mission to answer such questions.

Bailey, S. H.↗

Some examples of deep structure of the Archean from geophysics

The development of Archean crust remains as one of the significant problems in earth science, and a major unknown concerning Archean terrains is the nature of the deep crust. The character of crust beneath granulite terrains is especially fascinating because granulites are generally interpreted to represent a deep crustal section. Magnetic data from this area can be best modeled with a magnetized wedge of older Archean rocks (granulitic gneisses) underlying the younger Archean greenstone terrain. The dip of the boundary based on magnetic modeling is the same as the dip of the postulated thrust-fault reflection. Thus several lines of evidence indicate that the younger Archean greenstone belt terrain is thrust above the ancient Minnesota Valley gneiss terrain, presumably as the greenstone belt was accreted to the gneiss terrain, so that the dipping reflection represents a suture zone. Seismic data from underneath the granulite-facies Minnesota gneiss terrain shows abundant reflections between 3 and 6 s, or about 9 to 20 km. These are arcuate or dipping multicyclic events indicative of layering.

Smithson, S. B.↗

GRAIL Refinements to Lunar Seismic Structure

Gravity field measurements are perhaps the most numerous of the indirect observations relevant to the Moon's internal structure. Multiple recent missions have mapped the global lunar gravity field, each one improving upon the resolution of the last. The details of the Moon's deepest structure, including the parameters that define the lunar core, however, were still largely unaddressed by pre-GRAIL gravity measurements, which were not high enough resolution to resolve the tidal coefficients at a sufficient degree of accuracy. Current constraints on core size and state arise from other indirect measurements, including lunar laser ranging, magnetic induction studies, and analyses of elemental abundances in depth-derived mare basalts. These inferences vary widely, but when considered together with structure models derived from the seismic data gathered during the Apollo missions, a schematic of the lunar interior containing a partially molten deepest mantle layer overlying molten outer and solid inner core layers was obtained. Seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL mission's high-resolution measurements of the lunar gravity field are being used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL's constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) therefore complement seismic investigations. This work focuses on expanding our knowledge of the Moon's internal structure using joint gravity and seismic analyses, which will improve constraints on the deep lunar mantle and core.

Weber, Renee C.↗

Probabilistic Seismic Hazard Assessment for Azerbaijan

Probabilistic Seismic Hazard Assessments (PSHA) underpin the determination of seismic loads in most contemporary seismic provisions of building codes around the world. Modern building codes are migrating towards using the entire uniform hazard spectrum at a range of vibration periods (typically up to 4 s or 10 s) rather than a single peak value such as peak ground acceleration (PGA), thus requiring a larger range of PSHA outputs. The hazard maps in the current version of the building code of Azerbaijan (2011) are in terms of intensity and peak ground acceleration. Recognizing the need for an up-to-date seismic hazard assessment in the country, Seismic Cooperation Program (SCP) under the Lawrence Livermore National Laboratory (LLNL) undertook, in coordination with the Republican Seismic Survey Center of the Azerbaijan National Academy of Sciences and the Azerbaijan Scientific Research Institute of Construction and Architecture, a PSHA study that reflects new seismic data recorded locally and recent research conducted nationally and regionally since the last update to the building code. The PSHA framework for this project was designed to help develop a new earthquake catalogue, to incorporate a novel characterization of ground motions that specifically reflects the attenuation characteristics in the eastern Caucasus, and to generate hazard information in a form that is useful for an update of the current building code or the development of a new building code for Azerbaijan. The project also aimed to provide training and support for the local seismologists and engineers related to the seismic hazard models, probabilistic seismic hazard results and their use towards changes in the building code.

58 GEOSCIENCES↗

Advanced Shipboard Communications Demonstrations with ACTS

For ships at sea. satellites provide the only option for high data rate (HDR), long haul communications. Furthermore the demand for HDR satellite communications (SATCOM) for military and commercial ships. and other offshore platforms is increasing. Presently the bulk of this maritime HDR SATCOM connectivity is provided via C-band and X-band. However, the shipboard antenna sizes required to achieve a data rate of, say T 1 (1.544 Mbps) with present C-/X-band SATCOM systems range from seven to ten feet in diameter. This limits the classes of ships to which HDR services can be provided to those which are large enough to accommodate the massive antennas. With its high powered K/Ka-band spot beams, the National Aeronautics and Space Administration's (NASA) Advanced Communications Technology Satellite (ACTS) was able to provide T I and higher rate services to ships at sea using much smaller shipboard antennas. This paper discusses three shipboard HDR SATCOM demonstrations that were conducted with ACTS between 1996 and 1998. The first demonstration involved a 2 Mbps link provided to the seismic survey ship MN Geco Diamond equipped with a 16-inch wide, 4.5-inch tall, mechanically steered slotted waveguide array antenna developed by the Jet Propulsion Laboratory. In this February 1996 demonstration ACTS allowed supercomputers ashore to process Geco Diamond's voluminous oceanographic seismic data in near real time. This capability allowed the ship to adjust its search parameters on a daily basis based on feedback from the processed data, thereby greatly increasing survey efficiency. The second demonstration was conducted on the US Navy cruiser USS Princeton (CG 59) with the same antenna used on Geco Diamond. Princeton conducted a six-month (January-July 1997) Western Hemisphere solo deployment during which time T1 connectivity via ACTS provided the ship with a range of valuable tools for operational, administrative and quality-of-life tasks. In one instance, video teleconferencing (VTC) via ACTS allowed the ship to provide life-saving emergency medical aid, assisted by specialists ashore. to a fellow mariner - the Master of a Greek cargo ship. The third demonstration set what is believed to be the all-time SATCOM data rate record to a ship at sea, 45 Mbps in October 1998. This Lake Michigan (Chicago area) demonstration employed one of ACTS' fixed beams and involved the smallest of the three vessels, the 45-foot Bayliner M/V Entropy equipped with a modified commercial-off-the-shelf one-meter antenna. A variety of multi-media services were provided to Entropy through a stressing range of sea states. These three demonstrations provided a preview of the capabilities that could be provided to future mariners on a more routine basis when K/Ka-band SATCOM systems are widely deployed.

Axford, Roy A.↗

Determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution, appendix 2

This thesis reviews the technique established to clear channels in the Power Spectral Estimate by applying linear combinations of well known window functions to the autocorrelation function. The need for windowing the auto correlation function is due to the fact that the true auto correlation is not generally used to obtain the Power Spectral Estimate. When applied, the windows serve to reduce the effect that modifies the auto correlation by truncating the data and possibly the autocorrelation has on the Power Spectral Estimate. It has been shown in previous work that a single channel has been cleared, allowing for the detection of a small peak in the presence of a large peak in the Power Spectral Estimate. The utility of this method is dependent on the robustness of it on different input situations. We extend the analysis in this paper, to include clearing up to three channels. We examine the relative positions of the spikes to each other and also the effect of taking different percentages of lags of the auto correlation in the Power Spectral Estimate. This method could have application wherever the Power Spectrum is used. An example of this is beam forming for source location, where a small target can be located next to a large target. Other possibilities extend into seismic data processing. As the method becomes more automated other applications may present themselves.

Ioup, George E.↗

Global correlation of topographic heights and gravity anomalies

The short wavelength features were obtained by subtracting a calculated 24th-degree-and-order field from observed data written in 1 deg x 1 deg squares. The correlation between the two residual fields was examined by a program of linear regression. When run on a worldwide scale over oceans and continents separately, the program did not exhibit any correlation; this can be explained by the fact that the worldwide autocorrelation function for residual gravity anomalies falls off much faster as a function of distance than does that for residual topographic heights. The situation was different when the program was used in restricted areas, of the order of 5 deg x 5 deg square. For 30% of the world,fair-to-good correlations were observed, mostly over continents. The slopes of the regression lines are proportional to apparent densities, which offer a large spectrum of values that are being interpreted in terms of features in the upper mantle consistent with available heat-flow, gravity, and seismic data.

Roufosse, M. C.↗

Storage Field Development Plan: One Earth Energy

This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, and operational plans, and costs associated with the proposed One Earth Sequestration Carbon Capture and Storage (OES-CCS) site in McLean County, Illinois, near Gibson City. The proposed storage complex, known as the Mt. Simon Storage Complex, comprises the Cambrian Mt. Simon Sandstone reservoir and the primary seal, the Cambrian Eau Claire Formation. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Ordovician St. Peter Sandstone. Geologic characterization of the Mt. Simon Storage Complex at the OES-CCS site was performed by the Illinois Storage Corridor CarbonSAFE Phase III project, which also prepared and submitted three UIC Class VI applications to construct three injection wells; the permit applications were submitted and are in the federal EPA review process. A characterization well, OEE #1, was drilled to collect site-specific data. These data were analyzed and used to develop the UIC Class VI applications. The OEE #1 well will be converted to an in-zone monitoring (IZM) well for the injection phase. The proposed buildout for the OES-CCS site includes (1) three injection wells (OES #1, OES #2, and OES #3), (2) two IZM wells, (3) two above confining zone (ACZ) monitoring wells, one of which will be used to monitor the lowermost USDW, (4) capture and compression facilities, and (5) transportation facilities, i.e., pipelines. A pre-operational testing program was proposed in the Class VI permit application and will be employed at the site pending approval. Additional pre-injection (baseline), syn-injection, and post-injection monitoring and site care procedures will be followed by OES to ensure that injection activities are protective of human health and the environment. Injection is scheduled to begin in 2025, distributed across the three injection wells in accordance with the permit operating conditions. One Earth Sequestration intends to inject up to 90 million tonnes of CO 2 over a period of approximately 20 years. Injection will begin at approximately 0.5 million tonnes of CO 2 annually and ramp up to a maximum of 4.5 million tonnes annually. Daily injection rates are expected to range from 1,400 to 1,500 tonnes per day initially and reach a maximum of approximately 4,225 tonnes per day, depending on site geology and injectivity at each injection well location, and CO 2 availability. The costs associated with the OES-CCS project include pre-operational costs (e. g. additional seismic data acquisition and well drilling), capture and transportation facility and equipment costs, predicted field operating expenditures (OpEx), and decommissioning and post-injection site care (PISC) costs. The risks associated with project activities, such as site construction, injection operations, and verification of secure storage were evaluated, and mitigation strategies proposed to alleviate those risks.

09 BIOMASS FUELS↗

Thermoelastic Properties of Iron-Rich Ringwoodite and the Deep Mantle Aerotherm of Mars

The Martian mantle is considered to have a higher Fe/Mg ratio than the Earth's mantle. Ringwoodite, γ-(Mg,Fe) 2 SiO 4 , is likely the dominant polymorph of olivine in the core-mantle boundary (CMB) region of Mars. We synthesized anhydrous iron-rich ringwoodite with molar Mg/(Mg + Fe) = 0.44 and determined its thermal equation of state up to 35 GPa and 750 K by synchrotron X-ray diffraction. Using a third order Birch-Murnaghan equation of state, we obtain K T0 = 182 (3) GPa, K' = 4.6 (2), and α 0 = 3.18 (6) × 10 -5 K -1 . Using these results and an updated mineralogical model with an iron-rich composition of Mg/(Mg + Fe) = 0.75 for the Martian mantle, we estimate ~1900 K for the temperature of the D1000 seismic discontinuity inside Mars. The resulting adiabat predicts a warm aerotherm, which could explain the presence of partial melt at the CMB of Mars recently detected with seismic data from the 2019 InSight mission.

58 GEOSCIENCES↗

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗

Array-Based Seismic Measurements of OSIRIS-REx’s Re-Entry

The return home of the OSIRIS-REx spacecraft in September 2023 marked only the fifth time that an artificial object entered the Earth’s atmosphere at interplanetary velocities. Although rare, such events serve as valuable analogs for natural meteoroid re-entries; enabling study of hypersonic dynamics, shock wave generation, and acoustic-to-seismic coupling. Here, in this study, we report on the signatures recorded by a dense (100 m scale) 11-station array located almost directly underneath the capsule’s point of peak atmospheric heating in northern Nevada. Seismic data are presented, which allow inferences to be made about the shape of the shock wave’s footprint on the surface, the capsule’s trajectory, and its flight parameters.

58 GEOSCIENCES↗

Shallow Crustal Structure and Site Response at Fort Greely, Alaska

Fort Greely, Alaska is located in proximity to numerous faults capable of generating large (MW > 6.5) and damaging earthquakes. However, little knowledge exists about the noise properties, shallow crustal structure and site response of the area. We address this by examining ambient noise conditions, Rayleigh wave group wavespeed dispersion, and site response measurements using new and existing seismic data from December 2020 – April 2021 for four sites along a north-south line—stations GREN, GRES, PS09, K24K operated by the Alaska Earthquake Center. We find that ambient noise levels are relatively high on base compared to other places in Central Alaska, especially above 1 Hz, while noise levels in the microseismic band are typical of the region. Short-period (0.5 – 2Hz) Rayleigh wave group wavespeeds appear to decrease with increasing sediment thickness and increase when going southward from GREN to K24K. The predominant site frequency related to sediment thickness is lowest beneath PS09 and highest beneath K24K. We use the group velocities and site response values to propose a rough north-south layered model for the region.

58 GEOSCIENCES↗