Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismologic observations”

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 127 records · Page 7

A Review of Infrasound and Seismic Observations of Sample Return Capsules since the End of the Apollo Era in Anticipation of the OSIRIS-REx Arrival

Advancements in space exploration and sample return technology present a unique opportunity to leverage sample return capsules (SRCs) towards studying atmospheric entry of meteoroids and asteroids. Specifically engineered for the secure transport of valuable extraterrestrial samples from interplanetary space to Earth, SRCs offer unexpected benefits that reach beyond their intended purpose. As SRCs enter the Earth’s atmosphere at hypervelocity, they are analogous to naturally occurring meteoroids and thus, for all intents and purposes, can be considered artificial meteors. Furthermore, SRCs are capable of generating shockwaves upon reaching the lower transitional flow regime, and thus can be detected by strategically positioned geophysical instrumentation. NASA’s OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC is one of only a handful of artificial objects to re-enter the Earth’s atmosphere from interplanetary space since the end of the Apollo era and it will provide an unprecedented observational opportunity. This review summarizes past infrasound and seismic observational studies of SRC re-entries since the end of the Apollo era and presents their utility towards the better characterization of meteoroid flight through the atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects

Researchers at New Mexico Tech, Los Alamos National Laboratory (LANL), Sandia National Laboratory (SNL), and other collaborators are developing a methodology to improve characterization of stress in the subsurface by means of a model-based inversion of reservoir engineering data, time lapse seismic measurements, and microseismicity. Historically, various geophysical techniques have been used in efforts to understand the state of stress in the subsurface through direct and indirect imaging of stress sensitive features (faults and fractures) and observations of transient stress related observations (time variant elastic moduli and microseismicity). Standard direct and indirect techniques for seismic fault and fracture imaging suffer from detectability limits due to reliance on high data quality and multiplicity. Variations in effective elastic properties from inversion of high quality seismic (Vertical Seismic Profile (VSP)) data may be used to infer stress near the wellbore through integration with independent experimental characterization of the stress-velocity relationship. Given a sufficiently robust observation network, microseismic emissions may be inverted to characterize focal mechanisms which, together with supporting assumptions and constraints, inform estimates of in-situ stress. While each of these techniques contributes in part to the characterization of stress within limited spatial and temporal domains, no one method provides an unambiguous stress measurement or a predictive capability over a site scale spatial extent. Challenges associated with solution non-uniqueness, measurement ambiguity, and irregular sampling may be greatly minimized through combination of one or more of independent measurements within a common framework in which realistic geological, hydrodynamic, geomechanical, and seismological constituent process models may act as constraints.

02 PETROLEUM↗

Solar seismology. II - The stochastic excitation of the solar p-modes by turbulent convection

We test the hypothesis that the solar p-modes are stabilized by damping due to turbulent viscosity in the convective zone. Starting from the assumption that the modes are stable, we calculate expectation values for the modal energies. We find that the interaction between a p-mode and the turbulent convection is such that the modal energy tends toward equipartition with the kinetic energy of turbulent eddies whose lifetimes are comparable to the modal period. From the calculated values of the modal energies, we compute rms surface velocity amplitudes. Our predicted rms surface velocities range from 0.01 cm/sec for the fundamental radial mode to 0.6 cm/sec for the radial mode whose period is approximately 5 minutes. The predicted surface velocities for the low order p-modes are much smaller than the velocities inferred from recent observations.

Goldreich, P.↗

Probabilistic Programming for Transportable Source Characterization and Uncertainty Quantification of the North Korean Nuclear Tests 2006–2017

Here, we introduce a transportable technique to determine the yield and depth of burial (DOB) from seismic source spectra of underground nuclear explosions. We demonstrate this technique on the six declared North Korean nuclear tests. This approach derives source spectra in absolute units from regional phase (Pg) amplitudes by correcting the observations for geometric spreading, attenuation, and site amplification. We couple the source spectra and explosion source models with a probabilistic programming framework that integrates deep learning techniques and Bayesian modeling. This approach permits the exchange of information across various data categories to quantify both the data and model uncertainty. This technique stands out as an innovative use of broad‐area propagation models, making it transportable across various geologic settings. This method proves to be effective in scenarios with diverse and/or limited observational data, even when the source depth is unknown. We present new independent estimates of absolute yield and DOB that are consistent with the prior assessments, underscoring the potential of this method in enhancing transportable nuclear explosion monitoring capabilities.

58 GEOSCIENCES↗

On Geomagnetism and Paleomagnetism

A statistical description of Earth's broad scale, core-source magnetic field has been developed and tested. The description features an expected, or mean, spatial magnetic power spectrum that is neither "flat" nor "while" at any depth, but is akin to spectra advanced by Stevenson and McLeod. This multipole spectrum describes the magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Natural variations of core multipole powers about their mean values are to be expected over geologic time and are described via trial probability distribution functions that neither require nor prohibit magnetic isotropy. The description is thus applicable to core-source dipole and low degree non-dipole fields despite axial dipole anisotropy. The description is combined with main field models of modem satellite and surface geomagnetic measurements to make testable predictions of: (1) the radius of Earth's core, (2) mean paleomagnetic field intensity, and (3) the mean rates and durations of both dipole power excursions and durable axial dipole reversals. The predicted core radius is 0.7% above the 3480 km seismologic value. The predicted root mean square paleointensity (35.6 mu T) and mean Virtual Axial Dipole Moment (about 6.2 lx 1022 Am(exp 2)) are within the range of various mean paleointensity estimates. The predicted mean rate of dipole power excursions, as defined by an absolute dipole moment <20% of the 1980 value, is 9.04/Myr and 14% less than obtained by analysis of a 4 Myr paleointensity record. The predicted mean rate of durable axial dipole reversals (2.26/Myr) is 2.3% more than established by the polarity time-scale for the past 84 Myr. The predicted mean duration of axial dipole reversals (5533 yr) is indistinguishable from an observational value. The accuracy of these predictions demonstrates the power and utility of the description, which is thought to merit further development and testing. It is suggested that strong stable stratification of Earth's uppermost outer core leads to a geologically long interval of no dipole reversals and a very nearly axisymmetric field outside the core. Statistical descriptions of other planetary magnetic fields are outlined.

Voorhies, Coerte V.↗

Estimation of the orientation of stress in the Earth’s crust without earthquake or borehole data

Abstract Mechanical stress acting in the Earth’s crust is a fundamental property that is important for a wide range of scientific and engineering applications. The orientation of maximum horizontal compressive stress can be estimated by inverting earthquake source mechanisms and measured directly from borehole-based measurements, but large regions of the continents have few or no observations. Here we present an approach to determine the orientation of maximum horizontal compressive stress by measuring stress-induced anisotropy of nonlinear susceptibility, which is the derivative of elastic modulus with respect to strain. Laboratory and Earth experiments show that nonlinear susceptibility is azimuthally dependent in an anisotropic stress field and is maximum in the orientation of maximum horizontal compressive stress. We observe this behavior in the Earth—in Oklahoma and New Mexico, U.S.A, where maximum nonlinear susceptibility coincides with the orientation of maximum horizontal compressive stress measured using traditional methods. Our measurements use empirical Green’s functions and solid-earth tides and can be applied at different temporal and spatial scales.

58 GEOSCIENCES↗

The nature of the earth's core

The properties of the earth's core are overviewed with emphasis on seismologically determined regions and pressures and seismologically measured density, elastic wave velocities, and gravitational acceleration. Attention is given to solid-state convection of the inner core, and it is noted that though seismological results do not conclusively prove that the inner core is convective, the occurrence and magnitude of seismic anisotropy are explained by the effects of solid-state convection. Igneous petrology and geochemistry of the inner core, a layer at the base of the mantle and contact metasomatism at the core-mantle boundary, and evolution of the core-mantle system are discussed. It is pointed out that high-pressure melting experiments indicate that the temperature of the core is ranging from 4500 to 6500 K, and a major implication of such high temperature is that the tectonics and convection of the mantle, as well as the resulting geological processes observed at the surface, are powered by heat from the core. As a result of the high temperatures, along with the compositional contrast between silicates and iron alloy, the core-mantle boundary is considered to be most chemically active region of the earth.

Jeanloz, Raymond↗

Seismotectonic evidence for subduction beneath the Eastern Greater Caucasus

SUMMARY We generated high-resolution 3-D seismic velocity models as well as a relocated earthquake catalogue across the eastern Greater Caucasus and Kura basins. This work was done using data from the recently upgraded Republic Seismological Survey Center's (RSSC) seismic network. We generated our tomographic images of crustal velocity structure in Azerbaijan using double-difference inversions (i.e. tomoDD and hypoDD). Earthquake catalogues from the RSSC between 2011 and 2016 were used; these catalogues include absolute arrival times of 103 288 P- and 120 952 S-wave traveltime picks for 7574 events recorded at 35 stations in Azerbaijan. Beginning with a layered, 1-D velocity model that was estimated using VELEST, we inverted simultaneously for relative location, Vp and Vs on a 3-D grid with dimensions 670 × 445 × 45 km, with a uniform grid spacing of 55 × 55 × 5 km for all of Azerbaijan. We observe that the relocated hypocentres cluster into two depth ranges, at the surface and at depth, that appear to correspond to major fault zones and the top of a subducting plate. Additionally, we note intermediate depth seismicity (∼50–60 km) beneath the Kura Basin, and a northward deepening of earthquake depths. Seismic velocities vary significantly throughout the study region; we observe very slow velocities throughout the Kura Basin between 5 and 15 km, and elevated velocities at 20–35 km. The wholesale velocity structure and seismic structure of Kura Basin strongly mirrors that of the Caspian Sea, which suggests that the geodynamics of the Caspian continue westwards into Azerbaijan. The key results of this study suggest that the northward subduction observed in the Caspian Sea continues beneath the Eastern Greater Caucasus, as well as provides evidence for active faulting along the southern margin of the mountain range.

Gunnels, Michael↗

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗

Earth rotation and core topography

The NASA Geodynamics program has as one of its missions highly accurate monitoring of polar motion, including changes in length of day (LOD). These observations place fundamental constraints on processes occurring in the atmosphere, in the mantle, and in the core of the planet. Short-timescale (t less than or approx 1 yr) variations in LOD are mainly the result of interaction between the atmosphere and the solid earth, while variations in LOD on decade timescales result from the exchange of angular momentum between the mantle and the fluid core. One mechanism for this exchange of angular momentum is through topographic coupling between pressure variations associated with flow in the core interacting with topography at the core-mantel boundary (CMB). Work done under another NASA grant addressing the origin of long-wavelength geoid anomalies as well as evidence from seismology, resulted in several models of CMB topography. The purpose of work supported by NAG5-819 was to study further the problem of CMB topography, using geodesy, fluid mechanics, geomagnetics, and seismology. This is a final report.

Hager, Bradford H.↗

Observations and 3D MHD Modeling of Small-scale Solar Prominence Oscillations

Small scaleoscillations were recently observed in prominences with high resolution Hinode/SOT and IRIS SJI in unprecedented detailand identified as fast magnetosonic waves. The high-resolution observations provide evidence of flows of cool material, fluid instabilities such as Kelvin-Helmholtz and Rayleigh-Taylor, and the spectroscopic information provides evidence of Doppler shifts, heating, and cooling of the prominence material. Recently nonlinear fast magnetosonic waves were observed withperiods in the range of 5–11 minutes and wavelengths in the plane of the sky (POS) of ~2000kmand flows, and the flows in narrow threads withPOS speed ∼16–46 km/s. Thenonlinear fast magnetosonic waves and flows were modeled using 2.5D MHD (Ofman & Kucera 2020). Here, we extend the studyof the nonlinear fast magnetosonic waves using more realistic 3D MHD model. The model includes the 3D structure of the density and the magnetic fieldas well as mode coupling anddemonstrate the propagating density compressions associated with the nonlinear fast magnetosonicwaves, as well as the guided modesin the pillar. The 3DMHD modelingimprove understanding ofthe POSeffects in limb prominence oscillations observations.The results of the modeling are useful for coronal seismology of the prominence structure.

Leon Ofman↗

Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks

Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.

47 OTHER INSTRUMENTATION↗

Forecasting induced seismicity in Oklahoma using machine learning methods

Oklahoma earthquakes in the past decade have been mostly associated with wastewater injection. Here we use a machine learning technique—the Random Forest to forecast induced seismicity rate in Oklahoma based on injection-related parameters. We split the data into training (2011.01–2015.05) and test (2015.06–2020.12) periods. The model forecasts seismicity rate during the test period based on input features, including operational parameters (injection rate and pressure), geological information (depth to basement), and modeled pore pressure and poroelastic stress. The results show overall good match with observed seismicity rate (adjusted R 2 of 0.75). The model shows that pore pressure rate and poroelastic stressing rates are the two most important features in forecasting. The absolute values of pore pressure and poroelastic stress, and the injection rate itself, are less important than the stressing rates. These findings further emphasize that temporal changes of stressing rates would lead to significant changes in seismicity rates.

58 GEOSCIENCES↗

Nonlinear Fast Magnetosonic Waves in Solar Prominence Pillars

We investigate the properties of nonlinear fast magnetosonic (NFM) waves in a solar prominence, motivated by recent high-resolution and high-cadence Hinode/SOT observations of small-scale oscillations in a prominence pillar. As an example, we analyze the details of the 2012 February 14 Hinode/SOT observations of quasi-periodic propagating features consistent with NFM waves, imaged in emission in Ca II and in the far blue wing of Hα. We perform wavelet analysis and find oscillations in the 1-3 min period range. Guided by these observations, we model the NFM waves with a three-dimensional magnetohydrodynamics (3D MHD) model, extending previous 2.5D MHD studies. The new model includes the structure of the high-density, low-temperature material of the prominence pillar embedded in the hot corona, in both potential and non-force-free sheared magnetic field configurations. The nonlinear model demonstrates the effects of mode coupling and the propagating density compressions associated with linear and NFM waves. The guided fast magnetosonic waves, together with density compressions and currents, are reproduced in the 3D pillar structure. We demonstrate or the first time the dynamic effects of the Lorentz force due to the magnetic shear in the non-force-free field on the pillar structure and on the propagation of the waves. The insights gained from the 3D MHD modeling are useful for improving coronal seismology of prominence structures that exhibit fast MHD wave activity.

Leon Ofman↗

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.↗

Regolith, rock and fluid distributions at the Upper Colorado River Basin via a multicomponent seismic imaging approach

In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.

58 GEOSCIENCES↗

Constraints on the Lunar Interior Through Comparison of Observed Apollo Seismic Travel Time With Simulated Lunar Interiors

Understanding the current state of the lunar interior is key to evaluating its formation history and evolution. The analyses of seismic waves recorded by 5 years of the Apollo passive seismic network, from 1972-1977, have provided some fundamental constraints on the lunar interior structure. Here, we show that comparisons between (a) calculated seismic velocities of simulated lunar interiors and (b) measured Apollo seismic velocities can resolve new solutions. Specifically, we show that comparison of seismic travel times can be used to constrain the lower bound of possible selenotherms.

Moon↗

The Effect of Sonic Booms on Earthquake Warning Systems

Several aerospace companies are designing quiet supersonic business jets for service over the United States. These aircraft have the potential to increase the occurrence of mild sonic booms across the country. This leads to interest among earthquake warning (EQW) developers and the general seismological community in characterizing the effect of sonic booms on seismic sensors in the field, their potential impact on EQW systems, and means of discriminating their signatures from those of earthquakes. The SonicBREWS project (Sonic Boom Resistant Earthquake Warning Systems) is a collaborative effort between Seismic Warning Systems, Inc. (SWS) and NASA Dryden Flight Research Center. This project aims to evaluate the effects of sonic booms on EQW sensors. The study consists of exposing high-sample-rate (1000 sps) triaxial accelerometers to sonic booms with overpressures ranging from 10 to 600 Pa in the free field and the built environment. The accelerometers record the coupling of the sonic boom to the ground and surrounding structures, while microphones record the acoustic wave above ground near the sensor. Sonic booms are broadband signals with more high-frequency content than earthquakes. Even a 1000 sps accelerometer will produce a significantly aliased record. Thus the observed peak ground velocity is strongly dependent on the sampling rate, and increases as the sampling rate is reduced. At 1000 sps we observe ground velocities that exceed those of P-waves from M L 3 earthquakes at local distances, suggesting that sonic booms are not negligible for EQW applications. We present the results of several experiments conducted under SonicBREWS showing the effects of typical-case low amplitude sonic booms and worst-case high amplitude booms. We show the effects of various sensor placements and sensor array geometries. Finally, we suggest possible avenues for discriminating sonic booms from earthquakes for the purposes of EQW.

Wurman, Gilead↗