Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysical”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Investigating the Relationship Between Bolide Entry Angle and Apparent Direction of Infrasound Signal Arrivals

Infrasound sensing offers critical capabilities for detecting and geolocating bolide events globally. However, the observed back azimuths, directions from which infrasound signals arrive at stations, often differ from the theoretical expectations based on the bolide’s peak brightness location. For objects with shallow entry angles, which traverse longer atmospheric paths, acoustic energy may be emitted from multiple points along the trajectory, leading to substantial variability in back azimuth residuals. This study investigates how the entry angle of energetic bolides affects the back azimuth deviations, independent of extrinsic factors such as atmospheric propagation, station noise, and signal processing methodologies. A theoretical framework, the Bolide Infrasound Back-Azimuth EXplorer Model (BIBEX-M), was developed to compute predicted back azimuths solely from geometric considerations. The model quantifies how these residuals vary as a function of source-to-receiver distance, revealing that bolides entering at shallow angles, e.g., 10°, can produce average residuals of 20°, with deviations reaching up to 46° at distances below 1000 km, and remaining significant even at 5000 km (up to 8°). In contrast, bolides with steeper entry angles, e.g., > 60°, show smaller deviations, typically under 5° at 1000 km and diminishing to less than ~1° beyond 5000 km. These findings attest to the need for careful interpretation when evaluating signal detections and estimating bolide locations. This work is not only pertinent to bolides but also to other high-energy, extended-duration atmospheric phenomena such as space debris and reentry events, where similar geometric considerations can influence infrasound arrival directions.

Acoustics↗
DOE OSTI2025

Effect of a Fine-Scale Layered Structure of the Atmosphere on Infrasound Signals from Fragmenting Meteoroids

We investigate the influence of a fine-scale (FS) layered structure in the atmosphere on the propagation of infrasound signals generated by fragmenting meteoroids. Using a pseudo-differential parabolic equation (PPE) approach, we model broadband acoustic signals from point sources at altitudes of 35–100 km. The presence of FS fluctuations in the stratosphere (37–45 km) and the lower thermosphere (100–120 km) modifies ray trajectories, causing multiple arrivals and prolonged signal durations at ground stations. In particular, meteoroids fragmenting at 80–100 km can produce two distinct thermospheric arrivals beyond 150 km range, while meteoroids descending to 50 km or below yield weak, long-lived arrivals within the acoustic shadow zone via antiguiding propagation and diffraction. Comparison with observed infrasound data confirms that FS-layered inhomogeneities can account for multi-arrival “N-waves,” broadening potential interpretations of meteoroid signals. The results also apply to other atmospheric-entry objects, such as sample return capsules, emphasizing how FS structure impacts shock wave propagation. In conclusion, our findings advance understanding of wavefield evolution in a layered atmosphere and have broad relevance for global infrasound monitoring of diverse phenomena (e.g., re-entry capsules, rocket launches, and large-scale explosions).

Aeroacoustics↗
DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Closure of Earth’s Global Seasonal Cycle of Energy Storage

Abstract The global seasonal cycle of energy in Earth’s climate system is quantified using observations and reanalyses. After removing long-term trends, net energy entering and exiting the climate system at the top of the atmosphere (TOA) should agree with the sum of energy entering and exiting the ocean, atmosphere, land, and ice over the course of an average year. Achieving such a balanced budget with observations has been challenging. Disagreements have been attributed previously to sparse observations in the high-latitude oceans. However, limiting the local vertical integration of new global ocean heat content estimates to the depth to which seasonal heat energy is stored, rather than integrating to 2000 m everywhere as done previously, allows closure of the global seasonal energy budget within statistical uncertainties. The seasonal cycle of energy storage is largest in the ocean, peaking in April because ocean area is largest in the Southern Hemisphere and the ocean’s thermal inertia causes a lag with respect to the austral summer solstice. Seasonal cycles in energy storage in the atmosphere and land are smaller, but peak in July and September, respectively, because there is more land in the Northern Hemisphere, and the land has more thermal inertia than the atmosphere. Global seasonal energy storage by ice is small, so the atmosphere and land partially offset ocean energy storage in the global integral, with their sum matching time-integrated net global TOA energy fluxes over the seasonal cycle within uncertainties, and both peaking in April.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Neutron capture signatures determined from gadolinium isotope compositions of uranium ore concentrates using MC-ICPMS

We present Gd isotope compositions on 25 well-characterized uranium ore concentrates (UOCs). About half the UOCs have isotope depletions in 157 Gd coupled with anticorrelated excesses in 158 Gd due to neutron capture effects. UOCs from older ore bodies and with higher U contents show larger neutron capture effects than UOCs from younger ore bodies with lower U contents. These Gd data are correlated with their previously measured Sm isotope compositions, and we use the Sm data to estimate the neutron fluence of these samples. In conclusion, this work demonstrates how Gd isotope signatures can be employed as a new tool for nuclear forensics.

Earth Sciences↗
DOE OSTI2020

Improved sampling technique to collect natural gas from hydrate-bearing pressure cores

High quality gas compositional data are an important factor in interpreting the genetic source of natural gas hosted in hydrate-bearing sediments and other subsurface systems. In order to accurately characterize the composition of gas samples degassed from hydrate-bearing pressure cores, one must use a reproducible sampling technique that minimizes artifacts of the sampling process. Herein, we review sediment core degassing techniques and compare data obtained from a commonly used degassing approach, which we term the standard quantitative degassing (SQD) technique, to our newly developed modified quantitative degassing (MQD) method designed to minimize atmospheric contamination and gas-water interactions. The SQD method allows sample gas to interact with water in a bubbling chamber, which we hypothesize could alter the gas composition following mixing with water or dissolved gases in the bubbling chamber. Whereas, the MQD method allows for the collection of sample gas prior to the bubbling chamber. To compare the SQD and MQD methods, we performed a side-by-side comparison of noble (He, Ne, Ar, Kr, and Xe), major (H 2 , N 2 , O 2 , and CO 2 ), and hydrocarbon (CH 4 , C 2 H 6 , C 3 H 8 , i-C 4 H 10 , C 4 H 10 , i-C 5 H 12 , C 5 H 12 ) gas concentrations and select isotopic compositions obtained using both sample collection techniques. Gas samples were collected from hydrate-bearing pressure cores recovered and maintained under hydrate stable conditions from the northern Gulf of Mexico during the UT-GOM 2 -1 Expedition. The MQD method displayed significantly lower concentrations of atmospheric gases, higher proportions of hydrocarbon gases, lower ratios of C 1 /C 2 + , and heavier stable carbon and hydrogen isotopes of methane than the SQD method. These results demonstrate that the MQD method reduced air contamination and minimized alteration of the hydrocarbon gases. Finally, we conclude this method may be important for future work that seeks to determine the composition of natural gas from pressure cores using quantitative degassing experiments, especially those seeking to measure major (e.g., N 2 ) and noble gases.

58 GEOSCIENCES↗
DOE OSTI2021

Systematic investigations on iron cycling in phosphorus/siderophore systems: Synergism or antagonism?

Synergisms between microbial exudates on Fe (hydr)oxide dissolution as an effective Fe acquisition pathway have been recently addressed and vigorously debated. However, Fe liberation mechanisms and where siderophores and phosphorus (P) coexist received little attentions. Current study systematically investigated ferrihydrite dissolution in the presence of desferrioxamine B (DFOB) (a kind of fungally-derived siderophores) and inorganic/organic phosphorus (orthophosphate, Pi; myo-inositol hexaphosphate, IHP), as a function of solute pH, reaction time and reagent content. Reacted solids were characterized by N 2 -BET adsorption, zeta (ζ) potential analysis, sequential extraction analysis (SEDEX), field emission scanning electron microscopy (FESEM), high-resolution transmission electron microscopy (HRTEM), X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), micro Raman spectroscopy and attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR). Our results indicate that upon reaction with P-only or (DFOB + P) systems, interfacial complexation partially switched from monolayer bidentate-binuclear surface complexes to ternary complexes, or laterally transformed into amorphous Fe–P precipitates. The Fe-Pi complex precipitated more readily under acidic conditions, and Fe-IHP complex preferentially nucleated in neutral-alkaline environments. Phosphorus slightly promoted Fe release from minerals and fixation to the leached layer or interfacial liquid zone initially, but subsequently prevented further attacks from protons and DFOB. The co-effects of P and DFOB likely correspond to two successive scenarios: 1) DFOB is preferentially attracted toward ferrihydrite surfaces by negative electrical fields induced by adsorbed phosphorus and can act synergistically with labile P–Fe complexes, resulting in intensive temporal dissolution of Fh; 2) subsequent Fe shuttling to DFOB can be prohibited by stabilized, passive P/Fe–P layers. Finally, our results emphasize the antagonism between P compounds and siderophores (i.e., DFOB here) on ferrihydrite dissolution to improve understanding of the biologically-mediated Fe cycling in natural systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2021

Aspects of rare earth element enrichment in Allegheny Plateau coals, Pennsylvania, USA

The Allegheny Plateau in western Pennsylvania, western Maryland, and northern West Virginia has relatively flat-lying strata with southwest-northeast-trending anticlines cored by splay faults originating from the thrust faults and southeast-northwest-trending lineaments. The lineaments define both the structural relief of the anticlines and the coal rank. Coal rank in the Pennsylvanian Allegheny Group coals is medium volatile bituminous in the region of less-prominent structural relief and low volatile bituminous in the region with more pronounced anticlines. Zinc and Cl exhibit high values in coals in the vicinity of the lineaments. Although coal rank was not a parameter in the Principal components analysis, the highest rank coals clustered together in a quadrant defined, in part, by the Zn concentration. Zinc was selected as a parameter due to its possible link with hydrothermal mineralization and Zr and TiO 2 /Al 2 O 3 were selected based on their presumed association with detrital mineral influences. The highest rare earths + Y values in benches of the Lower Kittanning coal were found at sites close to the lineaments. Negative or negligible Ce anomalies (δCe) and strongly positive Gd anomalies (δGd) suggest marine or hydrothermal influences and negative Eu anomalies (δEu) along with high TiO 2 /Al 2 O 3 values suggest a terrestrial input. In conclusion, the plots of δCe vs. δGd and δCe vs. δEu show that the low volatile bituminous coals are segregated from the other coals, suggesting that the factors that influenced coal rank also influenced the distribution of the rare earth elements.

58 GEOSCIENCES↗