SEARCH · Search NASA
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.
We study spectra output from spatially proximal DPRK underground nuclear explosions from data challenged by several constraints.
Distributed acoustic sensing (DAS) has a demonstrated potential for wide-scale and continuous in situ monitoring of near-surface environmental and anthropogenic processes. DAS is attractive for development as a multi-geophysical observatory due to the prevalence of existing fiber infrastructure in regions with environmental, cultural, or strategic significance. To evaluate the efficacy of this technology for monitoring of polar environmental processes, we collected DAS data from a 37-km long section of seafloor telecommunications fiber located on the continental shelf of the Beaufort Sea, Alaska. This experiment spanned eight, one-week, seasonally-distributed periods across two years. This was the first ever deployment of seafloor DAS beneath sea ice, and the first deployment in any marine environment to span multiple seasons. We recorded a variety of environmental and anthropogenic signals with demonstrable utility for the study of sea ice dynamics and tracking of ocean vessels and ice-traversing vehicles.
Explore the source record for details and available documents.
The Terrestrial Planet Finder (TPF) was proposed as a mission concept to the 2000 Decadal Survey, and received a very high ranking amongst the major initiatives that were then reviewed. As proposed, it was a formation flying array of four 3-m class mid-infrared telescopes, linked together as an interferometer. Its science goal was to survey 150 nearby stars for the presence of Earth-like planets, to detect signs of life or habitability, and to enable revolutionary advances in high angular resolution astrophysics. The Decadal Survey Committee recommended that $200M be invested to advance TPF technology development in the Decade of 2000-2010. This paper presents the results of NASA's investment.
ABSTRACT Recent determinations of the mean free path of ionizing photons (mfp) in the intergalactic medium (IGM) at z = 6 are lower than many theoretical predictions. In order to gain insight, we investigate the evolution of the mfp in our new massive fully coupled radiation-hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDa III). CoDa III’s scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the IGM during reionization. The simulation was performed with ramses-cudaton on Summit, and used 13 1072 processors coupled to 24 576 GPUs, making it the largest reionization simulation, and largest ever ramses simulation. A superior agreement with global constraints on reionization is obtained in CoDa III over Cosmic Dawn II (CoDa II), especially for the evolution of the neutral hydrogen fraction and the cosmic photoionization rate, thanks to an improved calibration, later end of reionization (z = 5.6), and higher spatial resolution. Analysing the mfp, we find that CoDa III reproduces the most recent observations very well, from z = 6 to z = 4.6. We show that the distribution of the mfp in CoDa III is bimodal, with short (neutral) and long (ionized) mfp modes, due to the patchiness of reionization and the coexistence of neutral versus ionized regions during reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from 0.1 Mpc h−1 at z = 7 to ∼10 Mpc h−1 at z = 5.2. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at z > 13 (z < 5).
To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.
Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply scattered or reflected waves. However, a challenge arises when multiply scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural or induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (TNN also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The TNN with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. Furthermore, these results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.
When an acoustic wave strikes a topographic feature, some of its energy is scattered. Sensors on the ground cannot capture these scattered signals when they propagate at high angles. We report observations of upwardly-scattered acoustic waves prior to refraction back to the ground, intercepting them with a set of balloon-borne infrasound microbarometers in the lower stratosphere over northern Sweden. We show that these scattered waves generate a coda whose presence can be related to topography beneath balloons and low-altitude acoustic ducts. The inclination of the coda signals changes systematically with time, as expected from waves arriving from scatterers successively closer to receivers. The codas are present when a temperature inversion channels infrasound from a set of ground chemical explosions along the ground, but are absent following the inversion's dissipation. Since scattering partitions energy away from the main arrival, these observations imply a mechanism of amplitude loss that had previously been inaccessible to measurement. As such, these data and results allow for a better comprehension of interactions between atmospheric infrasound propagation and the solid earth
When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.
During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.
Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply-scattered or reflected waves. However, a challenge arises when multiply-scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The twin neural network with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. These results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.
Recent research indicates that 8 to 12% of the global oil and gas production methane emissions could be attributed to ultra-emitters, which result in high concentration ‘hotspots’ near point sources. Identifying these emissions in near real time provides useful information to the policy makers and private industry, who are working to reduce their impact. To meet this need, scientists are increasingly analyzing satellite data from the TROPOspheric Monitoring Instrument (TROPOMI) instrument aboard ESA’s Sentinel 5-Precursor mission. While direct analysis of TROPOMI level 2 swath data has been successful in identifying some large emission events, identifying hotpots is challenging because of the imaging noise due to a variety of artifacts and limits in daily coverage. Here we explore possible methodologies to detect methane hotspots using a new, gap-filled, and temporally continuous methane product from NASA’s Goddard Earth Observing System (GEOS) Constituent Data Assimilation System (CoDAS), which assimilates column averaged methane mole fractions from the TROPOMI with capabilities to assimilate other remote sensing measurements. The CoDAS has been expanded from a heritage of stratospheric composition and carbon dioxide assimilation allowing for the support of regional modeling, validation with non-coincident operations, and merging variety of datasets. The current work mainly explores Observing System Simulation Experiments with methane GEOS simulations without assimilation to prepare the groundwork for further experiments with assimilated TROPOMI. First, known hotspots based on the known inventory are identified to demonstrate the capability of the system to point out emission hotspots. In the next step, a variety of machine learning techniques such as Self-Organizing Maps and Deep Learning are explored to automate detection of the plumes. Finally, a few approaches to quantify emissions from the identified hotspots are presented and are evaluated against the inventory. The effort is directed toward a future evaluation of the CoDAS based methane monitoring system’s ability to successfully detect and quantify hotspots.
Ozone and water vapor in the lower stratosphere are important trace gases for atmospheric chemistry and radiative budget. The Stratospheric Aerosol and Gas Experiment (SAGE) missions have been crucial in monitoring the stratospheric ozone loss and the subsequent recovery as well as the trends in water vapor linked to surface temperature trends. The SAGE III instrument aboard the International Space Station (ISS) continues the SAGE mission record, with high vertical resolution profiles of ozone and water vapor available since mid 2017. The NASA GEOS Earth system model has the new capability to assimilate multi-constituents from ground and space-based instruments using the GEOS Constituent Data Assimilation System (CoDAS). The recently released MERRA-2 Stratospheric Composition Reanalysis with Aura MLS (M2-SCREAM) assimilates version 4.2 MLS ozone, water vapor and other chemically-reactive species with the NASA GEOS model coupled to a stratospheric-only chemistry mechanism and transport constrained to the MERRA-2 reanalysis. While the number of solar occultation observations a day from SAGE III/ISS is about 1% of the total number of profiles observed globally by MLS, the chemical timescales of ozone and water vapor in the lower stratosphere are long enough that the SAGE III/ISS data may provide a useful constraint on the assimilated product. Using the same GEOS CoDAS configuration as M2-SCREAM, we will present a series of experiments to investigate if ozone and water vapor trends are consistent with the assimilation of SAGE observations with and without Aura retrievals, and to determine if the assimilation of SAGE observations produces a steady product for trend analysis, especially as the end of the Aura mission nears. In our experiments, assimilating only SAGE III/ISS water vapor profiles results in water vapor fields more consistent with experiments that assimilate MLS v5; however, in the polar regions SAGE III/ISS observations are not available and the values are unconstrained. We are encouraged by the positive benefit assimilating the less frequent SAGE III/ISS observations has on stratospheric composition. Sensitivity experiments such as these will allow us to assess the added value of SAGE data for continued monitoring of the stratospheric composition for climate and ozone recovery assessments.
Abstract During routine operations monitoring Oklahoma earthquakes, we found that certain earthquakes occurred closely both in space and time and had overlapping phases at the recording stations. Through further scrutiny and analysis, we determined that rather than being distinctly different earthquakes, some of the earthquakes exhibited multiphase arrivals and longer than expected coda due to unique ray paths that encounter impedance contrasts such as at the sedimentary rock-basement. Of course, some of these events truly were distinct events, which we term overlapping earthquakes, for which perceived coda duration overlaps and obscures the phase arrivals of the second event due to the source proximity in both time and space. We detail our classification scheme to separate the local earthquakes in Oklahoma as single, overlapping earthquakes, or those associated with multiphase arrivals. We forward model seismic wave propagation in a 2D crustal model and develop a methodology that utilizes waveform correlation to distinguish phases from overlapping earthquakes to those from crustal reverberations. Duration analysis shows a more elongated duration, qualitatively similar to the duration produced by overlapping earthquakes, at the sites where multiphase arrivals are observed.
Discriminating low-yield underground nuclear explosions from small earthquakes is a key task in monitoring nuclear test ban treaties. P/S amplitude ratios have been an effective discriminant for moderate-sized events recorded at regional distances, but it is unclear if they are as effective in discriminating small seismic events recorded at local distances (<150 km). The difference between local magnitude (M L ) and coda duration magnitude (M c ) has been proposed as a new discriminant that may complement P/S amplitude ratios at local distances. Here, in this work, we calculate high-frequency (up to ∼4 Hz) synthetic seismograms at epicentral distances of 0–30 km in realistic models of the Salt Lake basin (Utah, United States) to better understand how variations in source type and depth affect M L −M c values. The Earth models incorporate simplified 1D and deterministic 3D structures, small-wavelength stochastic velocity perturbations, and surface topography. Coda waves are enhanced for the more complicated models compared to the base 1D model, but still underpredict observed durations by about a factor of two, which results in overprediction of amplitude to duration ratios (i.e., M L −M c values) for a near-surface explosion and a 7 km deep earthquake. For both source types, the predicted M L and M c values decrease as source depth increases, and M L −M c shows only minor variation with depth; however, M L −M c is on average ∼0.5 units smaller for explosions than earthquakes. This finding may imply that M L −M c has sensitivity to source type, in addition to being a depth discriminant, but more modeling is needed given the limitations of the current study. Future modeling should incorporate higher-frequency (≳5 Hz) simulations over a larger distance range (0–150 km), where M L and M c are commonly measured, while honoring low shear velocities (<300 m/s) near the surface and sampling a wider range of earthquake and explosion source mechanisms.
The Origins, Spectral Interpretation, Resource Identification, and Security‐Regolith Explorer sample return capsule (SRC) re‐entered the Earth’s atmosphere at hypersonic speeds from interplanetary space on 24 September 2023. The current work reports on 18 ground‐based acoustic sensors deployed at Wendover Airport, the same location that the Genesis and Stardust SRC re‐entries were recorded. Four different sensors (Chaparral Physics, Gem, Wilson Engineering Research and Development [WERD], and RedVox) were deployed in close proximity to compare their performance. All the sensors captured an N‐wave signal associated with the SRC re‐entry shock wave followed by a broadband coda. The Chaparral Physics array served as the high‐fidelity reference measurement. The N‐wave signal had a peak‐to‐peak amplitude of 4.07 Pa with a fundamental frequency of 4.98 Hz from 167.5° measured clockwise from north, nearly perpendicular to the SRC trajectory. In addition, high coherence in the coda was shown to be associated reflections from the surrounding mountains. In general, the more economical sensors (Gem, WERD, and RedVox) produced results that were consistent with these observations and sensor specifications. Beamforming with these single sensors arranged as an array showed agreement with the high‐fidelity array to within a couple of degrees. Furthermore, the current high‐fidelity results were compared with the measurements during the Genesis and Stardust SRC re‐entries. All three entries produced a broadband fundamental peak at a frequency that was inversely related to the SRC diameter as well as evidence of reflections from the surrounding topography.
Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.
Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.