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.
Item purpose – This poster will be presented at the SSA meeting Internal peer review status – This item has completed internal peer review. HQ review status – This item is being sent to HQ for review
A spectral-based correction adapted from ISO 11204 [1] is investigated here to mitigate the effects of ambient noise contamination of sonic booms recorded by noise monitors during field tests. The algorithm from sections 5.4.2 and 7 of ISO 11204 is customized using six nonstandard correction grades in addition to implementing the two standardized grades described in ISO 11204. The six nonstandard grades allow for more aggressive correction of the levels of the sonic boom spectrum when they are proximate to the ambient spectrum. Of the eight correction grades, the most aggressive custom correction grade, termed “Custom F” here, performed best under the conditions that were studied. Consequently, “Custom F” is recommended for use when processing in-field recordings of X-59 sonic booms. To evaluate the eight correction grades, mock X-59 acoustic recordings were generated using predicted ground-level X-59 sonic booms. These ground-level waveforms were created by Doebler [2] by propagating nearfield CFD solutions of the X-59 C612A configuration to the ground using PCBoom [3] (please see Acknowledgements and Refs. [4, 5, 2, 6]). To generate the mock acoustic recordings from those ground-level waveforms, they were further modified using turbulence filters developed during the NASA SonicBAT effort [7], post boom noise audio clips from the NASA SonicBAT tests [7], and ambient noise from the NASA QSF18 test [8, 9]. These mock acoustic recordings enabled evaluation of the ambient noise mitigation methods since the proper loudness levels of the X-59 waveforms in absence of ambient noise are known. Specifically, these known levels provide a benchmark against which the corrected loudness levels are compared, where the corrected loudness levels are computed when ambient noise is present within the waveforms. Importantly, similar analyses using in-field recordings are not possible since the proper loudness levels of the sonic boom waveforms in absence of ambient noise are unknowable when analyzing in-field recordings. Consequently, if additional analyses of ambient noise mitigation methods are needed, then use of mock recordings like those used here is recommended.
Characterizing and monitoring geologic formations around a borehole are crucial for energy and environmental applications. However, conventional wireline sonic logging usually cannot be used in high-temperature environments nor is the tool feasible for long-term monitoring. We introduce and evaluate the feasibility of a source-free distributed-acoustic-sensing (DAS) logging method based on borehole DAS ambient noise. Our new logging method provides a next-generation borehole imaging tool. The tool is source free because it uses ever-present ambient noises as sources and does not need a borehole sonic source that cannot be easily re-inserted into a borehole after well completion for time-lapse monitoring. The receivers of our source-free DAS logging tool are fiber optic cables cemented behind casing, enabling logging in harsh, high-temperature environments, and eliminating the receiver repeatability issue of conventional wireline sonic logging for time-lapse monitoring. We analyze a borehole DAS ambient noise dataset to obtain root-mean-squares (RMS) amplitudes and use these amplitudes to infer subsurface elastic properties. We find that the ambient noise RMS amplitudes correlate well with anomalies in conventional logging data. The source-free DAS logging tool can advance our ability to characterize and monitor subsurface geologic formations in an efficient and cost-effective manner, particularly in high-temperature environments such as geothermal reservoirs. Further validation of the source-free DAS logging method using other borehole DAS ambient noise data would enable the new logging tool for wider applications.
Research in planetary seismology is fundamentally constrained by a lack of data. Seismo-logical science products of future missions can typically only be informed by theoretical signal/noise characteristics of the environment or likely Earth-analogues. Although objectives can be re-assessed after some initial data-collection upon lander arrival, transfer of high-resolution data back to Earth is costly on lander power usage. Over the last several years, development of GPU computing techniques and open-source high-level APIs have led to rapid advances in deep learning within the fields of computer vision, natural language processing, and collaborative filtering. These techniques are actively being adapted in seismology for a variety of tasks, including: earthquake detection, seismic phase discrimination, and ground-motion prediction. Until the recent detection of mars quakes during the Mars InSight mission, the only other measurements of seismicity recorded outside of Earth was on the Moon during the Apollo missions between 1969 to 1977. These unique data sets have been periodically revisited using new seismological methods, including ambient noise interferometry and Hidden Markov Models. Our objective is to develop a deep learning seismic detector and use it to catalog moonquakes from the Apollo 17 Lunar Seismic Profiling Experiment (LSPE) and compare the results with those obtained by other methods. Additionally, we will assess the accuracy tradeoff between using a training set of lunar data and one composed of Earth seismicity. In this document, we present preliminary results using a prototype classifier trained on a small set of earthquakes that was able to obtain detections for LSPE moonquakes with a greater accuracy than a recent study using Hidden Markov Models.
As part of the Cyclotron Road program, Zanskar Geothermal & Minerals, Inc. investigated the application of micro-earthquake and ambient noise seismology methods to imaging and characterizing the structural characteristics and hydrothermal flux of subsurface faults. Significant advances in what could be resolved were enabled by two major developments in seismology: 1) the availability of large-n arrays of low-cost seismometers, and 2) the availability of increased computational power and semi-automated data reduction algorithms. In tandem, these advances may improve the signal-to-noise ratio and spatial precision of the data collected and enable higher-resolution characterization of subsurface fracture systems and their spatio-temporal evolution. These tools supported efforts to reduce dry-hole risk and to improve wellfield productivity for geothermal resource development. In particular, two applications of these advances were evaluated: 1) fracture-seismic imaging, which was used to detect ambient emissions from fluid-filled fractures, and 2) reservoir tomography, which used information about travel paths, source locations, and source parameters of micro-earthquakes to identify areas of enhanced permeability. Integration of these methods provided guidance for siting wells and served as prior constraints for reservoir models, informing forecasts of power potential and production and injection strategies aimed at minimizing temperature decline and improving overall resource productivity.
Data-intensive research in seismology is experiencing a recent boom, driven in part by large volumes of available data and advances in the growing field of data science. However, there are significant barriers to processing large data volumes, such as long retrieval times from data repositories, complex data management, and limited computational resources. New tools and platforms have reduced the barriers to entry for scientific cluster computing, including the maturation of the commercial cloud as an accessible instrument for research. Here, we build a customized research cluster in the cloud to test a new workflow for large-scale seismic analysis, in which data are processed as a stream (retrieved on-the-fly and acted upon without storing), with data from the Incorporated Research Institutions for Seismology Data Management Center. We use this workflow to deploy a spectral peak detection algorithm over 5.6 TB of compressed continuous seismic data from 2074 stations of the USArray Transportable Array EarthScope network. Using a 50-node cluster in the cloud, we completed the noise survey in 80 hr, with an average data throughput of 1.7 GB per minute. By varying cluster sizes, we find the scaling of our analysis to be sublinear, due to a combination of algorithmic limitations and data center response times. The cloud-based streaming workflow represents an order-of-magnitude increase in acquisition and processing speed compared to a traditional download-store-process workflow, and offers the additional benefits of employing a flexible, accessible, and widely used computing architecture. It is limited, however, due to its reliance on Internet transfer speeds and data center service capacity, and may not work well for repeated analyses or those for which even higher data throughputs are needed. These research applications will require a new class of cloud-native approaches in which both data and analysis are in the cloud.
Time-distance helioseismology is a method of ambient noise imaging using the solar oscillations. The basic realization that led to time-distance helioseismology was that the temporal cross correlation of the signals at two 'surface' (or photospheric) locations should show a feature at the time lag corresponding to the subsurface travel time between the locations. The temporal cross correlation, as a function of the location separation, is the Fourier transform of the spatio-temporal power spectrum of the solar oscillations, a commonly used function in helioseismology. It is therefore likely the characteristic ridge structure of the correlation function had been seen before without appreciation of its significance. Travel times are measured from the cross correlations. The times are sensitive to a number of important subsurface solar phenomena. These include sound speed variations, flows, and magnetic fields. There has been much interesting progress in the 17 years since the first paper on this subject (Duvall et al., Nature, 1993, 362, 430-432). This progress will be reviewed in this paper.
Near-surface seismic structure, particularly the shear wave velocity (V s ), can strongly affect local site response, and should be accurately estimated for ground motion prediction during seismic hazard assessment. The Imperial Valley (California), occupying the southern end of the Salton Trough, is a seismically active basin with thick surficial lacustrine sedimentary deposits. In this study, we utilize ambient noise records and local earthquake events for high-resolution near-surface characterization and site effect estimation with an unlit fiber-optic telecommunication infrastructure (dark fiber) in Imperial Valley by using the distributed acoustic sensing (DAS) technique. We apply ambient noise interferometry to retrieve coherent surface waves from DAS records, and evaluate performances of three different surface wave methods on DAS ambient noise dispersion imaging. We develop a quality control workflow to improve the dispersion measurement of noisy portions of the DAS data set by using a data selection strategy. Using the joint inversion of both the fundamental mode and higher overtones of Rayleigh waves, a high resolution two-dimensional (2D) V s structure down to 70 m depth is obtained. We successfully achieve an improved V s 30 (the time-averaged shear-wave velocity in the top 30 m) model with higher spatial-resolution and reliability compared to the existing community model for the area. We also explore the potential for utilizing DAS earthquake events for site amplification estimation. The preliminary results reveal a clear anti-correlation between the approximated site response and the V s 30 profile. In conclusion, our results indicate the potential utility of DAS deployed on dark fiber for near-surface characterization in appropriate contexts.
Explore the source record for details and available documents.
This report details seismic ambient noise analysis to improve seismic velocity estimates of the southeastern Nevada National Security Site (NNSS). We compare two different methods for estimating surface wave dispersion curves from ambient noise cross-correlations: frequency time analysis (FTAN) and Aki’s cross-spectral method (XSpec). We find that XSpec performs better for our local-regional dataset and frequency content. Using phase velocity estimates from XSpec, we build a preliminary phase velocity dispersion dataset, which contains data for 1,054 station pairs and 8,905 discrete phase velocity measurements for periods between 0.5 and 13 s. This phase velocity dataset will be incorporated into an updated local-regional P and S wave velocity model of the southeastern NNSS in the future and is expected to improve upon shallow velocity estimates.
Research studies of residents' responses to noise consist of interviews with samples of individuals who are drawn from a number of different compact study areas. The statistical techniques developed provide a basis for those sample design decisions. These techniques are suitable for a wide range of sample survey applications. A sample may consist of a random sample of residents selected from a sample of compact study areas, or in a more complex design, of a sample of residents selected from a sample of larger areas (e.g., cities). The techniques may be applied to estimates of the effects on annoyance of noise level, numbers of noise events, the time-of-day of the events, ambient noise levels, or other factors. Methods are provided for determining, in advance, how accurately these effects can be estimated for different sample sizes and study designs. Using a simple cost function, they also provide for optimum allocation of the sample across the stages of the design for estimating these effects. These techniques are developed via a regression model in which the regression coefficients are assumed to be random, with components of variance associated with the various stages of a multi-stage sample design.
A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.
A new physics-based method called Fundamental Rotorcraft Acoustic Modeling from Experiments (FRAME) is used to demonstrate the change in rotor harmonic noise of a helicopter operating at different ambient conditions. FRAME is based upon a non-dimensional representation of the governing acoustic and performance equations of a single rotor helicopter. Measured external noise is used together with parameter identification techniques to develop a model of helicopter external noise that is a hybrid between theory and experiment. The FRAME method is used to evaluate the main rotor harmonic noise of a Bell 206B3 helicopter operating at different altitudes. The variation with altitude of Blade-Vortex Interaction (BVI) noise, known to be a strong function of the helicopter s advance ratio, is dependent upon which definition of airspeed is flown by the pilot. If normal flight procedures are followed and indicated airspeed (IAS) is held constant, the true airspeed (TAS) of the helicopter increases with altitude. This causes an increase in advance ratio and a decrease in the speed of sound which results in large changes to BVI noise levels. Results also show that thickness noise on this helicopter becomes more intense at high altitudes where advancing tip Mach number increases because the speed of sound is decreasing and advance ratio increasing for the same indicated airspeed. These results suggest that existing measurement-based empirically derived helicopter rotor noise source models may give incorrect noise estimates when they are used at conditions where data were not measured and may need to be corrected for mission land-use planning purposes.
Understanding creep is essential for assessing seismic hazards along the creeping Hayward Fault, which threatens the San Francisco area. Current 1-D monitoring methods, like creepmeters and theodolites, provide limited measurements with low spatial resolution. InSAR offers higher 2-D spatial sampling but lower temporal resolution. We introduce a novel method using coda wave interferometry (CWI) to monitor creep in 2-D with high resolution. By calculating changes in ambient seismic noise travel times, we relate these measurements to creep offsets, bridging the gap between existing methods. Applied to seismic data from 2021 around the Hayward Fault, our CWI-derived offsets aligned with creepmeter measurements, revealing continuous and episodic creep. A 2-D regionalization identified a spatial boundary between these regimes, corroborating previous findings. This approach could enhance seismic monitoring and improve understanding of fault mechanics and earthquake risks.
A new methodology is developed for the construction of helicopter source noise models for use in mission planning tools from experimental measurements of helicopter external noise radiation. The models are constructed by employing a parameter identification method to an assumed analytical model of the rotor harmonic noise sources. This new method allows for the identification of individual rotor harmonic noise sources and allows them to be characterized in terms of their individual non-dimensional governing parameters. The method is applied to both wind tunnel measurements and ground noise measurements of two-bladed rotors. The method is shown to match the parametric trends of main rotor harmonic noise, allowing accurate estimates of the dominant rotorcraft noise sources to be made for operating conditions based on a small number of measurements taken at different operating conditions. The ability of this method to estimate changes in noise radiation due to changes in ambient conditions is also demonstrated.
Using a complex psychomotor task performed for 50 minutes in the presence of low frequency noise, high frequency noise, or ambient noise, annoyance ratings were obtained for noises of various frequencies by the method of magnitude estimation. The results suggest that high frequency noise affects female performance to a greater extent than male performance. Contrasted to these performance effects, the sexes did not differ in their annoyance ratings. A monotonically increasing relationship between annoyance and noise frequency was found (except for a decrease in annoyance at 8,000 Hz). It is concluded that both performance and annoyance responses may need to be assessed in certain situations to adequately describe human reaction to noise.
Propagation of sonic booms through turbulence reduces mean sonic boom perception metric levels and also causes considerable variability. NASA’s PCBoom suite of sonic boom acoustic propagation modules includes an approximate method for accounting for the effects of turbulence on traditional N-wave sonic booms. The current implementation is ineffective for shaped sonic booms or low-booms, and it also has limited values for turbulence and ambient input parameters. NASA’s future X-59 low-boom community noise surveys require an accurate estimate of the effects of turbulence in regions across the USA, so the module must be improved. This work presents the methods of selecting which ambient and turbulence parameters should be included in an improved PCBoom turbulence module. Turbulence and ambient data were collected from two atmospheric model databases, the Climate Forecast System Version 2 and European Centre for Medium-Range Weather Forecast Reanalysis Version 5 (ERA5), hourly from 7 AM to 7 PM local time for 10 years at 19 locations across the USA. A fully-factorial propagation analysis using these parameters would be exceedingly computationally expensive. Instead, a central composite design was chosen resulting in 45 combinations of ambient and turbulence parameters. These 45 cases effectively sample the space balancing computational burden.
SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.