SEARCH · Engineering Papers
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.
Hybrid Optimization Software Suite (HOSS) is based on finite-discrete element modeling to capture the complex physics of explosions. Calibration of near-source modeling for underground explosions in geophysical models is presented, showing the effects of the amplitude of radial and vertical seismic components.
The thickness of the crust, the thickness of the basal (intermediate or lower) crustal layer, and the average velocity at the top of the mantle have been mapped using all available deep-penetrating seismic-refraction profiles in the conterminous United States and surrounding border areas. These profiles are indexed to their literature data sources. The more significant long wavelength anomalies on the three maps are briefly discussed and analyzed. An attempt to use Bouguer gravity to validate mantle structure was inconclusive.
A relatively high-amplitude magnetic anomaly directly detected with the Apollo 15 subsatellite magnetometer and centered near the crater Gerasimovich on the southeastern lunar far side is found to correlate with the location of a conspicuous Reiner Gamma-type swirl marking visible on a Zond 8 photograph. Examinations of available direct and indirect orbital magnetics measurements demonstrate that most strong anomalies occur in areas where morphologically similar markings are concentrated. Even though photogeologic studies indicate an impact-related origin for the swirls, both the swirls and their associated strong anomalies tend to exist preferentially in or near areas that have been seismically modified. Modeling of improved vector magnetic anomaly maps is used to infer 28 independent bulk directions of magnetization for relatively strong and isolated lunar magnetic anomaly sources.
Here, in the study of subsurface seismic imaging, solving the acoustic wave equation is a pivotal component in existing models. The advancement of deep learning (DL) enables solving partial differential equations (PDEs), including wave equations, by applying neural networks to identify the mapping between the inputs and the solution. This approach can be faster than traditional numerical methods when numerous instances are to be solved. Previous works that concentrate on solving the wave equation by neural networks consider either a single velocity model or multiple simple velocity models, which is restricted in practice. Instead, inspired by the idea of operator learning, this work leverages the Fourier neural operator (FNO) to effectively learn the frequency domain seismic wavefields under the context of variable velocity models. We also propose a new framework paralleled FNO (PFNO) for efficiently training the FNO-based solver given multiple source locations and frequencies. Numerical experiments demonstrate the high accuracy of both FNO and PFNO with complicated velocity models in the OpenFWI datasets. Furthermore, the cross-dataset generalization test verifies that PFNO adapts to out-of-distribution velocity models. Finally, PFNO admits higher computational efficiency on large-scale testing datasets than the traditional finite-difference method. The aforementioned advantages endow the FNO-based solver with the potential to build powerful models for research on seismic waves.
SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.
The Construction Resource Utilization Explorer (CRUX) is a technology maturation project for the U.S. National Aeronautics and Space Administration to provide enabling technology for lunar and planetary surface operations (LPSO). The CRUX will have 10 instruments, a data handling function (Mapper - with features of data subscription, fusion, interpretation, and publication through geographical information system [GIs] displays), and a decision support system DSS) to provide information needed to plan and conduct LPSO. Six CRUX instruments are associated with an instrumented drill to directly measure regolith properties (thermal, electrical, mechanical, and textural) and to determine the presence of water and other hydrogen sources to a depth of about 2 m (Prospector). CRUX surface and geophysical instruments (Surveyor) are designed to determine the presence of hydrogen, delineate near subsurface properties, stratigraphy, and buried objects over a broad area through the use of neutron and seismic probes, and ground penetrating radar. Techniques to receive data from existing space qualified stereo pair cameras to determine surface topography will also be part of the CRUX. The Mapper will ingest information from CRUX instruments and other lunar and planetary data sources, and provide data handling and display features for DSS output. CRUX operation will be semi-autonomous and near real-time to allow its use for either planning or operations purposes.
Geothermal systems in the Basin and Range Province of the western United States probably differ in many respects from geothermal systems already discovered in other parts of the world because of the unique tectonic setting. To investigate this, a study of the geothermal occurrences at Fly Ranch, approximately 100 miles north of Reno, Nevada, has been undertaken. Ample evidence for a geothermal system exists in this area, including the surface expression of heat flow in the form of hot springs, an extensive area of low electrical resistivity, and a high level of seismicity along faults bounding the thermal area. However, geophysical and geological studies have not yet provided evidence for a local heat source at depth. Additional detailed geophysical and geological studies, as well as drilling, must be completed before the geothermal system can be described fully.
High efficiency dynamic Radioisotope Power Systems (RPS) could be mission enabling for low power space applications such as small probes, landers rovers, and communication repeaters. These applications would contain science instruments and be distributed across planetary surfaces or near objects of interest where solar flux is insufficient for using solar cells. Small RPS could be used to provide power for sensing radiation, temperature, pressure, seismic activity, and other measurements of interest to planetary scientists. Small RPS would use fractional versions of the General Purpose Heat Source (GPHS) or Light Weight Radioisotope Heater Units (LWRHU), to heat power conversion technologies. Dynamic power systems are capable of three to four times higher conversion efficiency compared to static power conversion technologies, and would provide an equal amount of power using less fuel or more power using an equal amount of fuel. Providing spacecraft with more power could decrease duty cycling of basic functions and, therefore, increase the quality and abundance of science data. NASA GRC is developing a low power dynamic RPS that would convert heat from multiple LWRHU to one watt of usable direct current electric power for spacecraft instrumentation and communication. The power system could be used to charge batteries or capacitors for higher power burst usage. The initial design, called Small Stirling Technology Exploration Power (smallSTEP), is around 3 kg, 11 cm diameter X 32 cm long, and converts 8 watts of heat to one watt of electricity using a Stirling convertor. This low power conversion system represents a new class of RPS with power levels two orders of magnitude lower than prototypes currently being developed for space applications under NASA contracts. Development of the 1-watt RPS includes maturation of convertor and controller designs, performance evaluation of an evacuated metal foil insulation, and development of system interfaces. Initial demonstration of the subsystems has been completed in a laboratory environment and a higher fidelity system is being pursued for demonstration in relevant environments for use on small spacecraft needed to carry out future space science missions.
High efficiency dynamic Radioisotope Power Systems (RPS) could be mission enabling for low power space applications such as small probes, landers rovers, and communication repeaters. These applications would contain science instruments and be distributed across planetary surfaces or near objects of interest where solar flux is insufficient for using solar cells. Small RPS could be used to provide power for sensing radiation, temperature, pressure, seismic activity, and other measurements of interest to planetary scientists. Small RPS would use fractional versions of the General Purpose Heat Source (GPHS) or Light Weight Radioisotope Heater Units (LWRHU), to heat power conversion technologies. Dynamic power systems are capable of three to four times higher conversion efficiency compared to static power conversion technologies, and would provide an equal amount of power using less fuel or more power using an equal amount of fuel. Providing spacecraft with more power could decrease duty cycling of basic functions and, therefore, increase the quality and abundance of science data. NASA Glenn Research Center (GRC) is developing a low power dynamic RPS that would convert heat from multiple LWRHU to one watt of usable direct current electric power for spacecraft instrumentation and communication. The power system could be used to charge batteries or capacitors for higher power burst usage. The initial design, called Small Stirling Technology Exploration Power (smallSTEP), is around 3 kg, 11 cm diameter X 32 cm long, and converts 8 watts of heat to one watt of electricity using a Stirling convertor. This low power conversion system represents a new class of RPS with power levels two orders of magnitude lower than prototypes currently being developed for space applications under NASA contracts. Development of the 1-watt RPS includes maturation of convertor and controller designs, performance evaluation of an evacuated metal foil insulation, and development of system interfaces. Initial demonstration of the subsystems has been completed in a laboratory environment and a higher fidelity system is being pursued for demonstration in relevant environments for use on small spacecraft needed to carry out future space science missions.
Background: Oleaginous yeast are prodigious producers of oleochemicals, offering alternative and secure sources for applications in foodstuff, skincare, biofuels, and bioplastics. Nitrogen starvation is the primary strategy used to induce oil accumulation in oleaginous yeast as part of a global stress response. While research has demonstrated that post-translational modifications (PTMs), including phosphorylation and protein cysteine thiol oxidation (redox PTMs), are involved in signaling pathways that regulate stress responses in metazoa and algae, their role in oleaginous yeast remain understudied and unexplored. Results: Towards linking the yeast oleaginous phenotype to protein function, we integrated lipidomics, redox proteomics, and phosphoproteomics to investigate Rhodotorula toruloides under nitrogen-rich and starved conditions over time. Our lipidomics results unearthed interactions involving sphingolipids and cardiolipins with ER stress and mitophagy. Our redox and phosphoproteomics data highlighted the roles of the AMPK, TOR, and calcium signaling pathways in regulation of lipogenesis, autophagy, and oxidative stress response. As a first, we also demonstrated that lipogenic enzymes including fatty acid synthase are modified as a consequence of shifts in cellular redox states due to nutrient availability. Conclusions: We conclude that lipid accumulation is largely a consequence of carbon rerouting and autophagy governed by changes to PTMs, and not increases in the abundance of enzymes involved in central carbon metabolism and fatty acid biosynthesis. Our systems-level approach sets the stage for acquiring multidimensional data sets for protein structural modeling and predicting the functional relevance of PTMs using Artificial Intelligence/Machine Learning (AI/ML). Coupled to those bioinformatics approaches, the putative PTM switches that we delineate will enable advanced metabolic engineering strategies to decouple lipid accumulation from nitrogen limitation.
The python routines are outlined in detail to perform the methods and results in the manuscript under review in the journal Geophysical Research Letters titled “Seismic features predict ground motions during repeating caldera collapse sequence” with LA-UR-23-33345. All routines are written in open source python and were applied to publicly available data sets. The codes formats the data into the appropriate structure required to train a boosted tree regression model. Other codes produce figure results.
Stressing and rupture of a locked zone adjacent to a creeping fault segment was studied with special reference to strength heterogeneity depthwise and along-strike. The resulting precursory temporal and spatial variations of surface strain rate profiles were compared to geodetic measurements on the San Andreas fault in central California. Crustal deformation in great California earthquake cycles was also studied with special reference to the temporal decay of strain rate observed since the 1957 and 1906 great earthquakes, and comtemporary surface strain rate and velocity profiles at several locations along the San Andreas. The effect of viscoelastic response in the deep aseismic shear zone on the surface deformation behavior was examined. Work was begun on a fundamental reformulation of the crustal deformation problem focusing on the crustal deformation process affected by deep aseismic slip as the slip zone progresses toward an instability and as deep seismic slip continues postseismically, the 3-D nature of the problem due to geometry and material heterogeneity, and the time-dependent source coming from the lithosphere/astenospheric coupling process.
The design, deployment, and operation of the Apollo 16 passive seismic experiment (PSE) are discussed. Since activation, all elements of the PSE have operated as planned, with the exception of the sensor thermal control system. Significant progress in the measurement of meteoroid flux in near-earth space has been made, along with dilineation of active moonquake source regions. The data obtained indicate that moonquakes are concentrated at great depth (800 to 1000 km) and that the apparent disparity between meteoroid flux estimtes based on lunar crater counts and those from earth-based observations can be resolved by seismic measurements in favor of the lower flux indicated by the crater count method. The results obtained from the PSE are summarized and their significance is discussed in detail.
This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.
Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.
Seismic data recorded at the Tonto Forest Seismological Observatory in Arizona and the Uinta Basin Seismological Observatory in Utah were used to compare the frequency of occurrence, severity, and spectral content of ground motions resulting from earthquakes, and other natural and man-made sources with the motions generated by sonic booms. A search of data recorded at the two observatories yielded a classification of over 180,000 earthquake phase arrivals on the basis of frequency of occurrence versus maximum ground velocity. The majority of the large ground velocities were produced by seismic surface waves from moderate to large earthquakes in the western United States, and particularly along the Pacific Coast of the United States and northern Mexico. A visual analysis of raw film seismogram data over a 3-year period indicates that local and regional seismic events, including quarry blasts, are frequent in occurrence, but do not produce ground motions at the observatories comparable to either the large western United States earthquakes or to sonic booms. Seismic data from the Nevada Test Site nuclear blasts were used to derive magnitude-distance-sonic boom overpressure relations.
We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.
Four areas of investigation, each dealing with the measurement of a particular geophysical property, are discussed. These properties are the gravity field, seismicity, magnetism, and heat flow. All are strongly affected by conditions, past or present, in the planetary interior; their measurement is the primary source of information about planetary interiors.