Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave arrival time”

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 19 records

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗

A Catalog of Temporally Localized Systematic Deviations in Global Body Wave Travel-Time Measurements

Accurate measurements of the arrival times of seismic waves are crucial for seismological analyses such as robust locations of earthquakes, characterization of seismic sources, and high-fidelity imaging of the Earth’s interior. However, these travel-time measurements can sometimes be contaminated by timing errors at the stations which record this data. In this study, we apply a classical approach, based on identifying time-dependence in measured body wave arrival times, to identify these timing errors in a dataset on the order of 107 individual measurements. We find timing deviations at a subset of the stations in our dataset and document the temporal location, extent, and severity of these errors, finding errors at 83 stations, and impacting ~100,000 measurements. This catalog of deviations may enable future investigators to obtain a more accurate dataset through the implementation of quality control measures to eliminate the contaminated data we have identified.

58 GEOSCIENCES↗

Practical Event Location Estimation Algorithm for Power Transmission System Based on Triangulation and Oscillation Intensity

Event location in power systems is quite essential information for system operators to enhance control-room situational awareness capability. Therefore, it is of great importance to develop an event location estimation algorithm for transmission systems with high accuracy. With the development of wide-area measurement system (WAMS) such as FNET/GridEye, and the synchrophasor measurement devices (SMDs) such as frequency disturbance recorders (FDRs), the synchronous measurement data including frequency, voltage amplitude and phase angle can be collected and used for event location estimation. First, the phase angle and rate of change of frequency (RoCoF) trajectories are respectively used for determining two sets of wave arrival time associated with each FDR. Then, a convolutional neural network (CNN) is utilized to determine the wave arrival order to select the more suitable set of wave arrival times for a given case and to perform corresponding modifications. Next, the oscillation intensity associated with each FDR is determined based on phase angle trajectories in the center of inertia (COI) coordinate system. Finally, the multiple criteria for event location estimation are represented. In conclusion, case studies and comparisons between the proposed and previous algorithms using actual and confirmed cases in U.S. power systems are performed to demonstrate the effectiveness and improvement of the proposed algorithm in practical applications.

frequency disturbance recorder (FDR)↗

Seismic structure of the lunar mantle - An overview

The direct P and S wave arrival times from natural lunar seismic events are the most complete and reliable data set for determining the structure of the lunar mantle. A total of 40 events provide sufficiently well-observed arrivals to permit the extraction of structural information. Using this arrival time data set, the average velocities in a two-layered mantle with an assumed crustal structure (from Toksoz et al., 1974) have been obtained. Reflected phases arriving after direct S are then used to calculate the depth of the boundary between the two mantle layers, and to demonstrate that it is probably a complex transition zone approximately 80 km thick between 400 and 480 km depth. To determine velocity gradients in the upper mantle it is required that the model explain the pronounced decrease in shear wave amplitudes and accompanying delay in shear wave arrival times beyond a distance of about 90 deg. The final model is well-constrained.

Goins, N. R.↗

Data-driven Event Identification in the U.S. Power Systems Based on 2D-OLPP and RUSBoosting Trees

Accurate event identification is an essential part of situation awareness ability for power system operators. Therefore, this work proposes an integrated event identification algorithm for power systems. First, to obtain and filter suitable inputs for event identification, an event detection trigger based on the rate of change of frequency (RoCoF) is presented. Then, the wave arrival time difference-based triangulation method considering the anisotropy of wave propagation speed is utilized to estimate the location of the detected event. Next, the two-dimensional orthogonal locality preserving projection (2D-OLPP)-based method, which is suitable for multiple types of measured data, is employed to achieve higher effectiveness in extracting the event features compared with traditional one-dimensional projection and principle component analysis (PCA). Finally, the random undersampling boosted (RUSBoosted) trees-based classifier, which can mitigate the data sample imbalance issue, is utilized to identify the type of the detected event. Furthermore, the proposed approach is demonstrated using the actual measurement data of U.S. power systems from FNET/GridEye. Comparison results show that the proposed event identification algorithm can achieve better performance than existing approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Lunar Seismology: the Internal Structure of the Moon

The direct P and S wave arrival times are the primary data set that can be measured on the seismograms of natural lunar seismic events. Polarization filtering techniques allow the enhancement of secondary body wave arrivals and record curves to identify the secondary phases and deduce structural information. Finally, shear wave amplitude vs. distance curves yield information on the location and magnitude of seismic velocity gradients in the interior. The results of these analyses show that the moon appears to have a two-layer crust at all four seismic stations: a 20 km upper crust that seems to be constant at all sites and a lower crust that is 40 km thick at stations 12 and 14 (mare), 55 + or - 10 km at station 16 (highland), and tentatively either 40 km or 70 km at station 15. The lower mantle extends from 480 km to at least 1100 km depth which is the maximum depth of penetration of all but a few seismic waves used as data. No definitive evidence for or against a lunar core exists.

Goins, N. R.↗

STILGAR: Subsurface Models for Graymont Pleasant Gap Mine

The detection, location, and monitoring of underground structures are of great importance to national and global security. Tunnels and voids generate seismic signatures detectable at the surface, but using non-invasive seismic data to image near-surface presents several challenges in real-world applications. In this report, we describe the use of a dense surface seismic deployment to generate subsurface models of the Graymont Pleasant Gap mine - a single-layer mine with a complex structure embedded in a high-velocity P-wave limestone bedrock. Our approach consists of three key methods. We use P-wave arrival times from local blast events to perform a tomography inversion with the tomoTD method, constructing a P-wave velocity model of the subsurface. We model the layer above the mine using Rayleigh wave ellipticity and inversion techniques. We leverage ongoing anthropogenic activities to identify and locate noise sources both on the surface and within the subsurface. With this integrated approach we aim to overcome the challenges and enhance our ability to non-invasively characterize underground structures, contributing to improved seismic monitoring techniques.

58 GEOSCIENCES↗

The lunar interior - A summary report

The complete seismic data set collected by the Apollo network contains about 40 events which provide significant structural information on the lunar interior. The seismograms from these events yield a set of direct wave arrival times that constitutes the most reliable information on the seismic structure of the moon. Secondary data include possible reflected arrivals from crustal and mantle interfaces, an apparent shear wave shadow zone for surface events beginning at about 90 deg distance, and the shear wave amplitude decay with distance. Analysis of these data give well-constrained and stable average velocity values for the upper and lower mantle regions independent of most assumptions. The upper-lower mantle transition can begin no shallower than 400 km depth and may represent a compositional change although the effects of increased temperature cannot be ruled out.

Goins, N. R.↗

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗

Evaluation of a preliminary regional Earth model through comparison of synthetic and observed waveform data

In this report, we document the process related to developing a regional geologic model of a 605 x 1334 km area centered around Utah and encompassing surrounding states. This model is developed to test the effect that composition of a model has on the generation of synthetic data with the intent of using this information to improve upon full waveform moment tensor inversions. We compare observed data from three seismic events and five stations to the synthetic data generated by a preliminary model derived from a geologic framework model (GFM) developed by the USGS. The synthetic data and observed data comparisons indicate that our preliminary model performs well at smaller offset distances in the northern and central sections of the model. However, the southern stations consistently display synthetic data P- and S-wave arrival times that do not match the observed data arrival times, indicating that the velocity structure of the southern part of the model especially is inaccurate.

58 GEOSCIENCES↗

Denoising Seismograms in the Time Domain Using a Deep Learning Model

Deep learning has emerged as a transformative tool for enhancing the extraction of reliable information from seismograms, addressing the increasing demand for precise and efficient seismic data analysis. We introduce an innovative encoder–decoder deep learning model, named WaveDenoiser, designed for noise reduction in the time domain, thereby eliminating the need for spectrogram computations that have been used for existing deep learning tools and significantly improving processing speed. Utilizing the benchmark dataset that is Stanford Earthquake Dataset, we developed three models of varying sizes: base, medium, and large. Notably, the large (referred to as WaveDenoiser) model demonstrated superior performance, achieving a median signal‐to‐noise ratio improvement of 8.8 dB on in‐distribution unseen data (in the same geographic region) and 7.7 dB on out‐distribution unseen data (in a new geographic region), outpacing both the base and medium models. Further evaluation of the WaveDenoiser model revealed a reduction in median arrival‐time errors by 0.02 s for P waves and 0.01 s for S waves when processing waveforms prior to phase picking using PhaseNet on in‐distribution unseen data. When tested on out‐distribution unseen data, the model also effectively reduced the P‐wave median arrival‐time error by 0.02 and 0.01 s in median arrival‐time error for S waves. Importantly, the application of WaveDenoiser resulted in a significant reduction of phase picking outliers by 1.1% to 3.6% for both P and S waves. In addition, we achieved over five times acceleration in processing speed compared with the seisBench implementation of DeepDenoiser. Our findings underscore the potential of WaveDenoiser as a powerful tool for improving seismic data analysis and processing efficiency.

P-waves↗

Discrimination of Seismic Events (2006–2020) in North Korea Using P/Lg Amplitude Ratios from Regional Stations and a Bivariate Discriminant Function

Two events of magnitude (m b ) 3.6–3.8 occurred in southern North Korea (NK) on 27 June 2019 and 11 May 2020. Although these events were located ~330–400 km from the known nuclear test site, the fact that they occurred within the territory of NK, a country with a recent history of underground nuclear tests, made them events of interest for the monitoring community. In this work, we used P/Lg ratios from regional stations to categorize seismic events that occurred in NK from 2006 to May 2020, including these two recent events, the six declared NK nuclear tests, and the cavity collapse and triggered earthquakes that followed the 3 September 2017 nuclear explosion. We were able to separate the cavity collapse from the population of nuclear explosions. However, based on P/Lg ratios, the distinction between the earthquakes and the cavity collapse is ambiguous. The performed discriminant analyses suggest that combining Pg/Lg and Pn/Lg ratios results in improved discriminant power compared with any of the ratio types alone. We used the two ratio types jointly in a quadratic discriminant function and successfully classified the six declared nuclear tests and the triggered earthquakes that followed the September 2017 explosion. Our analyses also confirm that the recent southern events of June 2019 and May 2020 are both tectonic earthquakes that occurred naturally.

58 GEOSCIENCES↗

Initial Results from the Variable Intensity Sonic Boom Propagation Database

An extensive sonic boom propagation database with low- to normal-intensity booms (overpressures of 0.08 lbf/sq ft to 2.20 lbf/sq ft) was collected for propagation code validation, and initial results and flight research techniques are presented. Several arrays of microphones were used, including a 10 m tall tower to measure shock wave directionality and the effect of height above ground on acoustic level. A sailplane was employed to measure sonic booms above and within the atmospheric turbulent boundary layer, and the sailplane was positioned to intercept the shock waves between the supersonic airplane and the ground sensors. Sailplane and ground-level sonic boom recordings were used to generate atmospheric turbulence filter functions showing excellent agreement with ground measurements. The sonic boom prediction software PCBoom4 was employed as a preflight planning tool using preflight weather data. The measured data of shock wave directionality, arrival time, and overpressure gave excellent agreement with the PCBoom4-calculated results using the measured aircraft and atmospheric data as inputs. C-weighted acoustic levels generally decreased with increasing height above the ground. A-weighted and perceived levels usually were at a minimum for a height where the elevated microphone pressure rise time history was the straightest, which is a result of incident and ground-reflected shock waves interacting.

Haering, Edward A., Jr.↗

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗

Delivery and application of precise timing for a traveling wave powerline fault locator system

The Bonneville Power Administration (BPA) has successfully operated an in-house developed powerline fault locator system since 1986. The BPA fault locator system consists of remotes installed at cardinal power transmission line system nodes and a central master which polls the remotes for traveling wave time-of-arrival data. A power line fault produces a fast rise-time traveling wave which emanates from the fault point and propagates throughout the power grid. The remotes time-tag the traveling wave leading edge as it passes through the power system cardinal substation nodes. A synchronizing pulse transmitted via the BPA analog microwave system on a wideband channel sychronizes the time-tagging counters in the remote units to a different accuracy of better than one microsecond. The remote units correct the raw time tags for synchronizing pulse propagation delay and return these corrected values to the fault locator master. The master then calculates the power system disturbance source using the collected time tags. The system design objective is a fault location accuracy of 300 meters. BPA's fault locator system operation, error producing phenomena, and method of distributing precise timing are described.

Street, Michael A.↗

Rock Valley Accelerated Weight Drop Preliminary P-wave Tomographic Model

An active source experiment using an accelerated weight drop was conducted in Rock Valley, Nevada National Security Site, during the spring of 2021 in order to characterize the shallow seismic structure of the region. P-wave first arrival travel times picked from this experiment were used to construct a preliminary 3-D compressional wave speed model over an area that is roughly 4 km wide east-west and 8 km north-south to a depth of about 500-600 m below the surface, but with primary data concentration along the transects of the experimental lines. The preliminary model shows good correlation with basic geology and surface features, but geological interpretation is not the focus of this report. We describe the methods used in the tomographic inversion of the data and show results from this preliminary P-wave model.

42 ENGINEERING↗

Field Evaluations of Sampling, Interviewing, and Flight Tracking of NASA's Low Boom Flight Demonstrator Aircraft

The first year’s effort identified sampling and interviewing as the principal risks to assessment of prompt reactions to overflights producing low-amplitude sonic booms. It also 1) established the utility of geo-information system-based route planning for LBFD flight missions, 2) developed and demonstrated a prototype of a geographically-distributed, Internet-enabled instrumentation system capable of wide-area tracking of LBFD aircraft in near-real time. The latter system permits synchronizing the conduct of interviews in multiple overflown communities with arrival times of shock waves at interviewing sites; and of measuring, archiving, and processing their acoustic signatures. Means were also recommended for constructing representative, telephone-based samples of eligible respondents living in households within carpet boom corridors adjacent to LBFD flight tracks, and for conducting interviews with cross-sectional (independent) samples of such respondents about their prompt reactions to exposure to low-amplitude sonic booms. A detailed study design was prepared and accepted by NASA for a set of single-contact attempt telephone interviews with a nationally representative sample of households. The study design focused on testing automated and live agent interview completion rates obtainable without callbacks. A minimal (two monitoring station) version of the aircraft tracking system was built and installed near a civil airport in a successful demonstration of the system’s ability to detect and track aircraft movements. The field exercise also demonstrated the ability of the system to capture the acoustic emissions of departing aircraft, and to serve aircraft position and sound level information to remote, geographically-distributed analysts in near-real time. Upon approval of OMB and IRB of the detailed study plan, a stratified, nationally representative sample of landline and wireless telephone-subscribing households was constructed. A total of 12,734 telephone interview contact attempts of the sort required by a straightforward cross-sectional study design were then made. These contact attempts demonstrated the impracticality of conducting a time-critical, cross-sectional study of prompt community response to low-amplitude sonic booms by means of “independent” (single contact attempt per respondent for each LBFD flight mission) telephone samples of respondents. The observed interview completion rates for these single telephone contact attempts were so low (~ 1% to 3% for automated and live agent interviews, respectively) that: 1) the representativeness of collected opinions would be susceptible to intuitive challenge as inadequate, even absent conclusive evidence of non-representativeness. Refuting challenges to representativeness would have to demonstrate that the composition of the actual sample did not differ from that of the target population, a task that is tantamount to proving a negative; 2) the information required to refute allegations of non-representativeness would require a questionnaire considerably lengthier than that required simply to determine the prevalence of boom-induced startle and annoyance. Such a questionnaire would have to inquire about potentially sensitive and intrusive matters, including respondents’ age, gender, education, employment, home ownership, income, ethnicity, family size, and other demographic factors; and 3) unreasonable numbers of attempts would be required to re-contact households with unsuccessful initial contact attempts, given the limited time available for doing so. For example, if about 500 completed interviews were desired in a supersonically overflown community, approximately 50,000 automated interview attempts would have to be made within ten to fifteen minutes of each LBFD overflight. Such large numbers of contact attempts could well exceed the numbers of households available for interview in areas of similar boom exposure levels in some communities near LBFD flight tracks. Such large numbers of interviews could be cost-effectively undertaken only by means of automated (i.e., outgoing interactive voice response) interviewing, a data collection method ill-suited for complex and sensitive questionnaire items. The infeasibility of independent sampling for evaluating prompt responses to LBFD overflights in a cross-sectional study is due in large part to simple non-response: that is, potential respondents – particularly those contacted on wireless telephones – refusing to answer calls with unfamiliar caller IDs. It is also due in part, however, to 1) the lack of time to attempt to contact the same respondent more than once within a few minutes after the arrival of a shock wave at the respondent’s location; and 2) the need to place calls during weekday/daytime hours, when response rates are notably lower than during evenings and weekends. Despite the poor interview completion rates achieved under the above constraints, cross sectional assessments of delayed reactions to LBFD overflights could still be feasible, if multiple attempts could be made to contact respondents during evening and weekend time periods, over extended time periods. Detailed plans for a longitudinal (panel) sample were developed as an alternative to a cross sectional sample design.

Fidell, Sanford↗